skip to main content

DOE PAGESDOE PAGES

This content will become publicly available on July 3, 2019

Title: Lagrangian discontinuous Galerkin hydrodynamic methods in axisymmetric coordinates

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 ofmore » 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.« less
Authors:
 [1] ; ORCiD logo [1] ; ORCiD logo [1]
  1. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Publication Date:
Report Number(s):
LA-UR-18-21449
Journal ID: ISSN 0021-9991
Grant/Contract Number:
AC52-06NA25396
Type:
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 373; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Research Org:
Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Sponsoring Org:
USDOE National Nuclear Security Administration (NNSA); LANL Laboratory Directed Research and Development (LDRD) Program
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; Lagrangian hydrodynamics; discontinuous Galerkin; cell-centered hydrodynamics; axisymmetric coordinates; true volume approach; area-weighted approach
OSTI Identifier:
1459834

Liu, Xiaodong, Morgan, Nathaniel R., and Burton, Donald E.. Lagrangian discontinuous Galerkin hydrodynamic methods in axisymmetric coordinates. United States: N. p., Web. doi:10.1016/j.jcp.2018.06.073.
Liu, Xiaodong, Morgan, Nathaniel R., & Burton, Donald E.. Lagrangian discontinuous Galerkin hydrodynamic methods in axisymmetric coordinates. United States. doi:10.1016/j.jcp.2018.06.073.
Liu, Xiaodong, Morgan, Nathaniel R., and Burton, Donald E.. 2018. "Lagrangian discontinuous Galerkin hydrodynamic methods in axisymmetric coordinates". United States. doi:10.1016/j.jcp.2018.06.073.
@article{osti_1459834,
title = {Lagrangian discontinuous Galerkin hydrodynamic methods in axisymmetric coordinates},
author = {Liu, Xiaodong and Morgan, Nathaniel R. and Burton, Donald E.},
abstractNote = {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.},
doi = {10.1016/j.jcp.2018.06.073},
journal = {Journal of Computational Physics},
number = ,
volume = 373,
place = {United States},
year = {2018},
month = {7}
}