Abstract:
The paper presents a numerical study of the convergence order of the projection-characteristic СРР (Cubic Polynomial Projection) method for solving a three-dimensional stationary transport equation on unstructured tetrahedral meshes. The method is based on a characteristic approach to solving the transport equation, has a minimal stencil within a single tetrahedron and a high (third) order of approximation. Unlike classical grid-characteristic methods, in this method, the final numerical approach is constructed not on the basis of interpolation operators of some order of approximation, but on the basis of orthogonal projection operators on the functional space used to approximate the solution. The base scheme is a one-dimensional scheme referred to the Hermitian cubic Interpolation СIP (Cubic Interpolation Polynomial) scheme. The use of interpolation operators is often implemented to sufficiently smooth functions. However, even if the exact solution has sufficient smoothness, some types of illumination of tetrahedra lead to the appearance of non-smooth grid solutions. The transition to orthogonal projectors solves two problems: firstly, the problem of the appearance of angular directions that are coplanar with the faces of the cells, and secondly, the problem of the appearance of non-smooth numerical solutions in the faces of the mesh cell. The convergence result is compared with the theoretical estimates obtained for the first time by one of the authors of this work. The third order of convergence of the method is shown, provided that the solution is sufficiently smooth and the absorption coefficient in the cells is near to constant.
Keywords:
transport equation, unstructured meshes, characteristic method.
Citation:
E. N. Aristova, G. O. Astafurov, “A third-order projection-characteristic method for solving the transport equation on unstructed grids”, Mat. Model., 35:11 (2023), 79–93; Math. Models Comput. Simul., 16:2 (2024), 208–216
\Bibitem{AriAst23}
\by E.~N.~Aristova, G.~O.~Astafurov
\paper A third-order projection-characteristic method for solving the transport equation on unstructed grids
\jour Mat. Model.
\yr 2023
\vol 35
\issue 11
\pages 79--93
\mathnet{http://mi.mathnet.ru/mm4506}
\crossref{https://doi.org/10.20948/mm-2023-11-06}
\transl
\jour Math. Models Comput. Simul.
\yr 2024
\vol 16
\issue 2
\pages 208--216
\crossref{https://doi.org/10.1134/S2070048224020066}
Linking options:
https://www.mathnet.ru/eng/mm4506
https://www.mathnet.ru/eng/mm/v35/i11/p79
This publication is cited in the following 2 articles:
E. N. Aristova, G. O. Astafurov, “Vysokotochnaya skhema dlya uravneniya perenosa v zadache neitronnoi zaschity”, Preprinty IPM im. M. V. Keldysha, 2024, 013, 21 pp.
E. N. Aristova, N. I. Karavaeva, I. R. Ivashkin, “Monotonizatsiya modifitsirovannoi skhemy s ermitovoi interpolyatsiei dlya chislennogo resheniya neodnorodnogo uravneniya perenosa s pogloscheniem”, Preprinty IPM im. M. V. Keldysha, 2024, 065, 40 pp.