Lagrangian discontinuous Galerkin hydrodynamic methods in axisymmetric coordinates

Liu, Xiaodong; Morgan, Nathaniel R.; Burton, Donald E.

doi:10.1016/j.jcp.2018.06.073

Title: Lagrangian discontinuous Galerkin hydrodynamic methods in axisymmetric coordinates

Full Record
Other Related Research

Abstract

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 »« less

Authors:

Liu, Xiaodong ^[1];

^[1];

^[1]

Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

Publication Date:: Tue Jul 03 00:00:00 EDT 2018

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

OSTI Identifier:: 1459834

Alternate Identifier(s):: OSTI ID: 1532565

Report Number(s):: LA-UR-18-21449
Journal ID: ISSN 0021-9991; TRN: US1901810

Grant/Contract Number:: AC52-06NA25396

Resource Type:: Accepted Manuscript

Journal Name:: Journal of Computational Physics

Additional Journal Information:: Journal Volume: 373; Journal ID: ISSN 0021-9991

Publisher:: Elsevier

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

Citation Formats


                    Liu, Xiaodong, Morgan, Nathaniel R., and Burton, Donald E. Lagrangian discontinuous Galerkin hydrodynamic methods in axisymmetric coordinates.  United States: N. p., 2018. 
Web.  doi:10.1016/j.jcp.2018.06.073.

Copy to clipboard


                    Liu, Xiaodong, Morgan, Nathaniel R., & Burton, Donald E. Lagrangian discontinuous Galerkin hydrodynamic methods in axisymmetric coordinates.  United States.  https://doi.org/10.1016/j.jcp.2018.06.073

Copy to clipboard


                    Liu, Xiaodong, Morgan, Nathaniel R., and Burton, Donald E. Tue .  
"Lagrangian discontinuous Galerkin hydrodynamic methods in axisymmetric coordinates".  United States.  https://doi.org/10.1016/j.jcp.2018.06.073.  https://www.osti.gov/servlets/purl/1459834.

Copy to clipboard


                    
@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         = {Tue Jul 03 00:00:00 EDT 2018},

  month        = {Tue Jul 03 00:00:00 EDT 2018}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (Publisher)

Accepted Manuscript (DOE)

Publisher's Version of Record

https://doi.org/10.1016/j.jcp.2018.06.073

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 14 works

Citation information provided by
Web of Science

Save / Share:

Export Metadata

Save to My Library

Similar Records in DOE PAGES and OSTI.GOV collections:

A high-order Lagrangian discontinuous Galerkin hydrodynamic method for quadratic cells using a subcell mesh stabilization scheme

Journal Article Liu, Xiaodong ; Morgan, Nathaniel R. ; Burton, Donald E. - Journal of Computational Physics

We present here a Lagrangian discontinuous Galerkin (DG) hydrodynamic method that is up to third-order accurate using subcell mesh stabilization (SMS) for compressible flows on quadratic meshes in two-dimensional (2D) Cartesian coordinates. Similar to the second-order accurate Lagrangian DG method with linear meshes, the physical evolution equations for the specific volume, velocity, and specific total energy are discretized using a modal DG method with Taylor series polynomials. The Riemann velocity at the vertices of a curvilinear cell, and the corresponding surface forces, are calculated by solving a multidirectional approximate Riemann problem. Curvilinear cells (e.g., quadratic quadrilateral meshes in this work)more »« less
Cited by 24
https://doi.org/10.1016/j.jcp.2019.02.008

Full Text Available
A Godunov-like point-centered essentially Lagrangian hydrodynamic approach

Journal Article Morgan, Nathaniel R. ; Waltz, Jacob I. ; Burton, Donald E. ; ... - Journal of Computational Physics

We present an essentially Lagrangian hydrodynamic scheme suitable for modeling complex compressible flows on tetrahedron meshes. The scheme reduces to a purely Lagrangian approach when the flow is linear or if the mesh size is equal to zero; as a result, we use the term essentially Lagrangian for the proposed approach. The motivation for developing a hydrodynamic method for tetrahedron meshes is because tetrahedron meshes have some advantages over other mesh topologies. Notable advantages include reduced complexity in generating conformal meshes, reduced complexity in mesh reconnection, and preserving tetrahedron cells with automatic mesh refinement. A challenge, however, is tetrahedron meshesmore »« less
Cited by 18
https://doi.org/10.1016/j.jcp.2014.10.048
A point-centered arbitrary Lagrangian Eulerian hydrodynamic approach for tetrahedral meshes

Journal Article Morgan, Nathaniel R. ; Waltz, Jacob I. ; Burton, Donald E. ; ... - Journal of Computational Physics

We present a three dimensional (3D) arbitrary Lagrangian Eulerian (ALE) hydrodynamic scheme suitable for modeling complex compressible flows on tetrahedral meshes. The new approach stores the conserved variables (mass, momentum, and total energy) at the nodes of the mesh and solves the conservation equations on a control volume surrounding the point. This type of an approach is termed a point-centered hydrodynamic (PCH) method. The conservation equations are discretized using an edge-based finite element (FE) approach with linear basis functions. All fluxes in the new approach are calculated at the center of each tetrahedron. A multidirectional Riemann-like problem is solved atmore »« less
Cited by 16
https://doi.org/10.1016/j.jcp.2015.02.024

Full Text Available
A Lagrangian discontinuous Galerkin hydrodynamic method

Journal Article Liu, Xiaodong ; Morgan, Nathaniel Ray ; Burton, Donald E. - Computers and Fluids

Here, we present a new Lagrangian discontinuous Galerkin (DG) hydrodynamic method for solving the two-dimensional gas dynamic equations on unstructured hybrid meshes. The physical conservation laws for the momentum and total energy are discretized using a DG method based on linear Taylor expansions. Three different approaches are investigated for calculating the density variation over the element. The first approach evolves a Taylor expansion of the specific volume field. The second approach follows certain finite element methods and uses the strong mass conservation to calculate the density field at a location inside the element or on the element surface. The thirdmore »« less
Cited by 32
https://doi.org/10.1016/j.compfluid.2017.12.007

Full Text Available
RAM: a Relativistic Adaptive Mesh Refinement Hydrodynamics Code

Technical Report Zhang, Wei-Qun ; /KIPAC, Menlo Park ; MacFadyen, Andrew I ; ...

The authors have developed a new computer code, RAM, to solve the conservative equations of special relativistic hydrodynamics (SRHD) using adaptive mesh refinement (AMR) on parallel computers. They have implemented a characteristic-wise, finite difference, weighted essentially non-oscillatory (WENO) scheme using the full characteristic decomposition of the SRHD equations to achieve fifth-order accuracy in space. For time integration they use the method of lines with a third-order total variation diminishing (TVD) Runge-Kutta scheme. They have also implemented fourth and fifth order Runge-Kutta time integration schemes for comparison. The implementation of AMR and parallelization is based on the FLASH code. RAM ismore »« less
https://doi.org/10.2172/890471

Full Text Available

Similar Records