Finite difference methods for continuum mechanics problems on triangular and tetrahedral grids

Yu. N. Karamzin ; I. V. Popov ; S. V. Polyakov

Matematičeskoe modelirovanie, Tome 15 (2003) no. 11, pp. 3-12

Résumé

The solution of 2-D and 3-D continuum mechanics problems with help of finite difference schemes on unstructured grids are discussed in the report. The first problem is the discretization of computational domain. To solve it we use some new algorithms for decomposition of 2-D and 3-D multiply connected non-convex domains in “rough” triangles and tetrahedrons with consequent application of special procedures for its refinement. To this end the discrete local analyzis of grid points for given contours is fulfilled. Then the “rough” triangular or tetrahedral grid is constructed and after its refinement we get the grid with given properties satisfying Delaunay criteria. For grid refinement we propose two original algorithms. The first one is connected with initial generation of adaptive grid to domain boundary and to the solution of the problem. The second one gives the opportunity to construct a grid using the given boundary and domain mesh points. The first algorithm has $O(N)$ and the second $O(N\sqrt N)$ operation, where $N$ is the final number of triangles (tetrahedrons). The second problem discussed in the report is the construction of finite difference schemes on such unstructured grids. We propose original finite difference schemes with higher order of accuracy for parabolic equations and for transport equation. In these schemes the values of all functions are prescribed only to the mesh points. That gives the possibility to solve the problems with arbitrary types of boundary conditions and to use the adaptive grids without interpolation procedure.