×

zbMATH — the first resource for mathematics

High-order curvilinear finite elements for axisymmetric Lagrangian hydrodynamics. (English) Zbl 1290.76061
Summary: In this paper we present an extension of our general high-order curvilinear finite element approach for solving the Euler equations in a Lagrangian frame to the case of axisymmetric problems. The numerical approximation of these equations is important in a number of applications of compressible shock hydrodynamics and the reduction of 3D problems with axial symmetry to 2D computations provides a significant computational advantage. Unlike traditional staggered-grid hydrodynamics (SGH) methods, which use the so-called “area-weighting” scheme, we formulate our semi-discrete axisymmetric conservation laws directly in 3D and reduce them to a 2D variational form in a meridian cut of the original domain. This approach is a natural extension of the high-order curvilinear finite element framework we have developed for 2D and 3D problems in Cartesian geometry, leading to a rescaled momentum conservation equation which includes new radial terms in the pressure gradient and artificial viscosity forces. We show that this approach exactly conserves energy and we demonstrate via computational examples that it also excels at preserving symmetry in problems with symmetric initial conditions. The results also illustrate that our computational method does not produce spurious symmetry breaking near the axis of rotation, as is the case with many area-weighted approaches.

MSC:
76M10 Finite element methods applied to problems in fluid mechanics
65M60 Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs
Software:
GLVis; BLAST; MFEM
PDF BibTeX XML Cite
Full Text: DOI
References:
[1] Dobrev VA, Kolev TzV, Rieben RN. High order curvilinear finite element methods for Lagrangian hydrodynamics. SIAM J Sci Comput, submitted for publication.
[2] Ockendon, H.; Ockendon, J., Waves and compressible flow, Texts in applied mathematics, vol. 47, (2004), Springer-Verlag · Zbl 1041.76001
[3] Benson, D. J., Computational methods in Lagrangian and Eulerian hydrocodes, Comput Methods Appl Mech Eng, 99, 235-394, (1992) · Zbl 0763.73052
[4] Zel’dovich, Ya. B.; Raizer, Yu. P., Physics of shock waves and high-temperature hydrodynamic phenomena, (1967), Academic Press New York
[5] Castor, J., Radiation hydrodynamics, (2004), Cambridge University Press Cambridge, UK
[6] Dobrev, V. A.; Ellis, T. E.; Kolev, Tz. V.; Rieben, R. N., Curvilinear finite elements for Lagrangian hydrodynamics, Int J Numer Methods Fluids, 65, 11-12, 1295-1310, (2010) · Zbl 1255.76075
[7] Wilkins, M. L., Methods in computational physics, vol. 3, (1964), Academic Press, [chapter Calculation of elastic-plastic flow]
[8] Whalen, P., Algebraic limitations on two dimensional hydrodynamics simulations, J Comput Phys, 124, 1, 46-54, (1996) · Zbl 0849.76079
[9] Caramana, E. J.; Burton, D. E.; Shashkov, M. J.; Whalen, P. P., The construction of compatible hydrodynamics algorithms utilizing conservation of total energy, J Comput Phys, 146, 227-262, (1998) · Zbl 0931.76080
[10] Margolin, L.; Shashkov, M., Using a curvilinear grid to construct symmetry-preserving discretizations for Lagrangian gas dynamics, J Comput Phys, 149, 389-417, (1999) · Zbl 0936.76057
[11] Barlow, A., A compatible finite element multi-material ALE hydrodynamics algorithm, Int J Numer Methods Fluids, 56, 953-964, (2007) · Zbl 1169.76030
[12] Scovazzi, G.; Christon, M. A.; Hughes, T. J.R.; Shadid, J. N., Stabilized shock hydrodynamics: I. A Lagrangian method, Comput Methods Appl Mech Eng, 196, 4-6, 923-966, (2007) · Zbl 1120.76334
[13] Scovazzi, G.; Love, E.; Shashkov, M., Multi-scale Lagrangian shock hydrodynamics on Q1/P0 finite elements: theoretical framework and two-dimensional computations, Comput Methods Appl Mech Eng, 197, 1056-1079, (2008) · Zbl 1169.76396
[14] Lipnikov, K.; Shashkov, M., A framework for developing a mimetic tensor artificial viscosity for Lagrangian hydrocodes on arbitrary polygonal meshes, J Comput Phys, 229, 20, 7911-7941, (2010) · Zbl 1425.76186
[15] Matejovic, P.; Adamik, V., A one-point integration quadrilateral with hourglass control in axisymmetric geometry, Comput Methods Appl Mech Eng, 70, 3, 301-320, (1988) · Zbl 0658.73051
[16] Matejovic, P., Quadrilateral with high coarse-mesh accuracy for solid mechanics in axisymmetric geometry, Comput Methods Appl Mech Eng, 88, 2, 241-258, (1991) · Zbl 0745.73016
[17] Maire, P-H., A high-order cell-centered Lagrangian scheme for two-dimensional compressible fluid flows on unstructured meshes, J Comput Phys, 228, 7, 2391-2425, (2009) · Zbl 1156.76434
[18] Maire, P-H., A high-order cell-centered Lagrangian scheme for compressible fluid flows in two-dimensional cylindrical geometry, J Comput Phys, 228, 6882-6915, (2009) · Zbl 1261.76021
[19] Sedov, L. I., Similarity and dimensional methods in mechanics, (1993), CRC Press
[20] Burton DE. Multidimensional discretizations of conservation laws for unstructured polyhedral grids. Technical report UCRL-JC-118306. Lawrence Livermore National Laboratory; 1994.
[21] Caramana, E. J.; Shashkov, M. J., Elimination of artificial grid distortion and hourglass-type motions by means of Lagrangian subzonal masses and pressures, J Comput Phys, 142, 2, 521-561, (1998) · Zbl 0932.76068
[22] BLAST: high-order finite element Lagrangian hydrocode. <http://www.llnl.gov/CASC/blast>.
[23] MFEM: modular parallel finite element methods library, <http://mfem.googlecode.com>.
[24] GLVis: OpenGL visualization tool. <http://glvis.googlecode.com>.
[25] Dukowicz, J. K.; Meltz, B. J., Vorticity errors in multidimensional Lagrangian codes, J Comput Phys, 99, 1, 115-134, (1992) · Zbl 0743.76058
[26] Galera, S.; Maire, P-H.; Breil, J., A two-dimensional unstructured cell-centered multi-material ALE scheme using VOF interface reconstruction, J Comput Phys, 229, 16, 5755-5787, (2010) · Zbl 1346.76105
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.