Lagrangian discontinuous Galerkin hydrodynamic methods in axisymmetric coordinates
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
In this paper, we present new Lagrangian discontinuous Galerkin (DG) hydrodynamic methods for compressible flows on unstructured meshes in axisymmetric coordinates. The physical evolution equations for the specific volume, velocity, and specific total energy are discretized using a modal DG method with linear Taylor series polynomials. Two different approaches are used to discretize the evolution equations – the first one is the true volume approach and the second one is the area-weighted approach. For the true volume approach, the DG equations are derived using the true 3D volume that is consistent with the geometry conservation law (GCL). The Riemann velocity at the nodes on the surface of the element, and the corresponding surface forces, are calculated by solving a multidirectional approximate Riemann problem using surfaces areas for axisymmetric coordinates. This true volume approach conserves mass, momentum, and total energy and satisfies the GCL. However, it can not preserve spherical symmetry on an equal-angle polar grid with 1D radial flows. For the area-weighted approach, the DG equations are based on the 2D Cartesian geometry that is rotated about the axis of symmetry using a single, element average radius. With this approach, the Riemann velocity at the nodes on the surface of the element, and the corresponding surface forces, are calculated by solving a multidirectional approximate Riemann problem in 2D Cartesian geometry. This area-weighted approach, in the limit of an infinitesimal mesh size, conserves physical momentum, and physical total energy. The area-weighted approach preserves spherical symmetry on an equal-angle polar grid for 1D radial flows, but it does not satisfy the GCL. Finally, a suite of test problems are calculated to demonstrate stable mesh motion, the expected second order accuracy of these methods, and that the new area-weighted DG method preserves spherical symmetry on 1D radial flow problems with equal-angle polar meshes.
- Research Organization:
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
- Sponsoring Organization:
- USDOE National Nuclear Security Administration (NNSA); LANL Laboratory Directed Research and Development (LDRD) Program
- Grant/Contract Number:
- AC52-06NA25396
- OSTI ID:
- 1459834
- Alternate ID(s):
- OSTI ID: 1532565
- Report Number(s):
- LA-UR-18-21449; TRN: US1901810
- Journal Information:
- Journal of Computational Physics, Vol. 373; ISSN 0021-9991
- Publisher:
- ElsevierCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
Similar Records
A Godunov-like point-centered essentially Lagrangian hydrodynamic approach
A point-centered arbitrary Lagrangian Eulerian hydrodynamic approach for tetrahedral meshes