Hostname: page-component-848d4c4894-8kt4b Total loading time: 0 Render date: 2024-06-16T21:28:16.380Z Has data issue: false hasContentIssue false

Solution Reconstruction on Unstructured Tetrahedral Meshes Using P1-Conservative Interpolation

Published online by Cambridge University Press:  08 July 2016

Biao Peng*
Affiliation:
Department of Aerodynamics, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
Chunhua Zhou*
Affiliation:
Department of Aerodynamics, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
Junqiang Ai*
Affiliation:
The First Aircraft Institute, Aviation Industry Corporation of China, Xi'an 710089, China
*
*Corresponding author. Email:pb_nuaa@126.com (B. Peng), chzhou@nuaa.edu.cn (C. H. Zhou), avic603ajq@163.com (J. Q. Ai)
*Corresponding author. Email:pb_nuaa@126.com (B. Peng), chzhou@nuaa.edu.cn (C. H. Zhou), avic603ajq@163.com (J. Q. Ai)
*Corresponding author. Email:pb_nuaa@126.com (B. Peng), chzhou@nuaa.edu.cn (C. H. Zhou), avic603ajq@163.com (J. Q. Ai)
Get access

Abstract

This paper extends an algorithm of P1-conservative interpolation on triangular meshes to tetrahedral meshes and thus constructs an approach of solution reconstruction for three-dimensional problems. The conservation property is achieved by local mesh intersection and the mass of a tetrahedron of the current mesh is calculated by the integral on its intersection with the background mesh. For each current tetrahedron, the overlapped background tetrahedrons are detected efficiently. A mesh intersection algorithm is proposed to construct the intersection of a current tetrahedron with the overlapped background tetrahedron and mesh the intersection region by tetrahedrons. A localization algorithm is employed to search the host units in background mesh for each vertex of the current mesh. In order to enforce the maximum principle and avoid the loss of monotonicity, correction of nodal interpolated solution on tetrahedral meshes is given. The performance of the present solution reconstruction method is verified by numerical experiments on several analytic functions and the solution of the flow around a sphere.

Type
Research Article
Copyright
Copyright © Global-Science Press 2016 

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Jiao, X. and Heath, M., Common-refinement-based data transfer between non-matching meshes in multiphysics simulations, Int. J. Numer. Methods Eng., 61(14) (2004), pp. 24022427.Google Scholar
[2]Ringler, T. and Randall, D., A potential entropy and energy conserving numerical scheme for solution of the shallow-water equations on a geodesic grid, Monthly Weather Review, 130(5) (2002), pp. 13971410.Google Scholar
[3]Thuburn, J., Some conservation issues for the dynamical cores of NWP and climate models, J. Comput. Phys., 227(7) (2007), pp. 37153730.Google Scholar
[4]Cullen, M., Modeling atmospheric flows, Acta Numerica, 16 (2007), pp. 67154.Google Scholar
[5]Alauzet, F. and Mehrenberger, M., P1-conservative solution interpolation on unstructured triangular meshes, Int. J. Numer. Methods Eng., 84 (2010), pp. 15521588.Google Scholar
[6]Alauzet, F., Frey, P., George, P. and Mohammadi, B., 3D transient fixed point mesh adaptation for time-dependent problems: application to CFD simulations, J. Comput. Phys., 222 (2007), pp. 592623.Google Scholar
[7]Dukowicz, J. and Kodis, J., Accurate conservative remapping (rezoning) for arbitrary Lagrangian-Eulerian computations, SIAM J. Sci. Statist. Comput., 8(3) (1987), pp. 305321.Google Scholar
[8]Garimella, R., Kucharik, M. and Shashkov, M., An efficient linearity and bound preserving conservative interpolation (remapping) on polyhedral meshes, Comput. Fluids, 36(2) (2007), pp. 224237.Google Scholar
[9]Margolin, L. and Shashkov, M., Second-order sign-preserving conservative interpolation (remapping) on general grids, J. Comput. Phys., 184(1) (2003), pp. 266298.Google Scholar
[10]Jones, P. W., First- and second-order conservative remapping schemes for grids in spherical coordinates, Monthly Weather Review, 127(9) (1999), pp. 22042210.Google Scholar
[11]Cheng, J. and Shu, C. W., A high order accurate conservative remapping method on staggered meshes, Appl. Numer. Math., 58(7) (2008), pp. 10421060.Google Scholar
[12]Grandy, J., Conservative remapping and regions overlays by intersecting polyhedral, J. Comput. Phys., 148(2) (1999), pp. 433466.Google Scholar
[13]Dukowicz, J. K., Conservative rezoning (remapping) for general quadrilateral meshes, J. Comput. Phys., 54 (1984), pp. 411424.Google Scholar
[14]Geuzaine, C., Meys, B., Henrotte, F., Dular, P. and Legros, W., A Galerkin projection method for mixed finite elements, IEEE Transactions on Magnetics, 35(3) (1999), pp. 14381441.Google Scholar
[15]Farrell, P., Piggott, M., Pain, C. and Gorman, G., Conservative interpolation between unstructured meshes via supermesh construction, Comput. Methods Appl. Mech. Eng., 198 (2009), pp. 26322642.Google Scholar
[16]Farrell, P. E. and Maddison, J. R., Conservative interpolation between volume meshes by local Galerkin projection, Comput. Methods Appl. Mech. Eng., 200 (2011), pp. 89100.Google Scholar
[17]Zhou, C. H. and Ai, J. Q., Mesh adaptation for simulation of unsteady flow with moving immersed boundaries, Int. J. Numer. Methods Fluids, 72 (2013), pp. 453477.Google Scholar
[18]Frey, P. and George, P., Mesh Generation, Application to Finite Elements (2nd edn), ISTE Ltd and Wiley: New York, 2008.Google Scholar
[19]Löhner, R., Applied CFD Techniques, An Introduction Based on Finite Element Methods, Wiley: New York, 2001.Google Scholar
[20]Preparata, F. P. and Shamos, M. I., Computational Geometry: An Introduction, Second ed., Springer, 1985.Google Scholar
[21]Seidel, R., Convex Hull Computations, in: Goodman, J. E. and óRourke, J. (Eds.), Handbook of Discrete and Computational Geometry, Chapman & Hall/CRC, 2004, pp. 495512.Google Scholar