Streamline integration as a method for two-dimensional elliptic grid generation
- Institute for Ion Physics and Applied Physics, Universität Innsbruck, A-6020 Innsbruck (Austria)
- Numerical Analysis group, Universität Innsbruck, A-6020 Innsbruck (Austria)
We propose a new numerical algorithm to construct a structured numerical elliptic grid of a doubly connected domain. Our method is applicable to domains with boundaries defined by two contour lines of a two-dimensional function. Furthermore, we can adapt any analytically given boundary aligned structured grid, which specifically includes polar and Cartesian grids. The resulting coordinate lines are orthogonal to the boundary. Grid points as well as the elements of the Jacobian matrix can be computed efficiently and up to machine precision. In the simplest case we construct conformal grids, yet with the help of weight functions and monitor metrics we can control the distribution of cells across the domain. Our algorithm is parallelizable and easy to implement with elementary numerical methods. We assess the quality of grids by considering both the distribution of cell sizes and the accuracy of the solution to elliptic problems. Among the tested grids these key properties are best fulfilled by the grid constructed with the monitor metric approach. - Graphical abstract: - Highlights: • Construct structured, elliptic numerical grids with elementary numerical methods. • Align coordinate lines with or make them orthogonal to the domain boundary. • Compute grid points and metric elements up to machine precision. • Control cell distribution by adaption functions or monitor metrics.
- OSTI ID:
- 22622303
- Journal Information:
- Journal of Computational Physics, Journal Name: Journal of Computational Physics Vol. 340; ISSN JCTPAH; ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
The adaptive mesh method based on HOC difference scheme for convection diffusion equations with boundary layers
Adaptive grid generation by elliptic equations with grid control at all of the boundaries