skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A positivity-preserving semi-implicit discontinuous Galerkin scheme for solving extended magnetohydrodynamics equations

Authors:
; ORCiD logo;
Publication Date:
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1228388
Grant/Contract Number:
NA0001836; NA0001855
Resource Type:
Journal Article: Published Article
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 278; Journal Issue: C; Related Information: CHORUS Timestamp: 2015-03-19 03:52:33; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Country of Publication:
United States
Language:
English

Citation Formats

Zhao, Xuan, Yang, Yang, and Seyler, Charles E. A positivity-preserving semi-implicit discontinuous Galerkin scheme for solving extended magnetohydrodynamics equations. United States: N. p., 2014. Web. doi:10.1016/j.jcp.2014.08.044.
Zhao, Xuan, Yang, Yang, & Seyler, Charles E. A positivity-preserving semi-implicit discontinuous Galerkin scheme for solving extended magnetohydrodynamics equations. United States. doi:10.1016/j.jcp.2014.08.044.
Zhao, Xuan, Yang, Yang, and Seyler, Charles E. Mon . "A positivity-preserving semi-implicit discontinuous Galerkin scheme for solving extended magnetohydrodynamics equations". United States. doi:10.1016/j.jcp.2014.08.044.
@article{osti_1228388,
title = {A positivity-preserving semi-implicit discontinuous Galerkin scheme for solving extended magnetohydrodynamics equations},
author = {Zhao, Xuan and Yang, Yang and Seyler, Charles E.},
abstractNote = {},
doi = {10.1016/j.jcp.2014.08.044},
journal = {Journal of Computational Physics},
number = C,
volume = 278,
place = {United States},
year = {Mon Dec 01 00:00:00 EST 2014},
month = {Mon Dec 01 00:00:00 EST 2014}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1016/j.jcp.2014.08.044

Citation Metrics:
Cited by: 8works
Citation information provided by
Web of Science

Save / Share:
  • Shallow water equations with a non-flat bottom topography have been widely used to model flows in rivers and coastal areas. An important difficulty arising in these simulations is the appearance of dry areas where no water is present, as standard numerical methods may fail in the presence of these areas. These equations also have still water steady state solutions in which the flux gradients are nonzero but exactly balanced by the source term. In this paper we propose a high order discontinuous Galerkin method which can maintain the still water steady state exactly, and at the same time can preservemore » the non-negativity of the water height without loss of mass conservation. A simple positivity-preserving limiter, valid under suitable CFL condition, will be introduced in one dimension and then extended to two dimensions with rectangular meshes. Numerical tests are performed to verify the positivity-preserving property, well-balanced property, high order accuracy, and good resolution for smooth and discontinuous solutions.« less
  • Here, a space and time third-order discontinuous Galerkin method based on a Hermite weighted essentially non-oscillatory reconstruction is presented for the unsteady compressible Euler and Navier–Stokes equations. At each time step, a lower-upper symmetric Gauss–Seidel preconditioned generalized minimal residual solver is used to solve the systems of linear equations arising from an explicit first stage, single diagonal coefficient, diagonally implicit Runge–Kutta time integration scheme. The performance of the developed method is assessed through a variety of unsteady flow problems. Numerical results indicate that this method is able to deliver the designed third-order accuracy of convergence in both space and time,more » while requiring remarkably less storage than the standard third-order discontinous Galerkin methods, and less computing time than the lower-order discontinous Galerkin methods to achieve the same level of temporal accuracy for computing unsteady flow problems.« less
  • A comparative study of two classes of third-order implicit time integration schemes is presented for a third-order hierarchical WENO reconstructed discontinuous Galerkin (rDG) method to solve the 3D unsteady compressible Navier-Stokes equations: — 1) the explicit first stage, single diagonally implicit Runge-Kutta (ESDIRK3) scheme, and 2) the Rosenbrock-Wanner (ROW) schemes based on the differential algebraic equations (DAEs) of Index-2. Compared with the ESDIRK3 scheme, a remarkable feature of the ROW schemes is that, they only require one approximate Jacobian matrix calculation every time step, thus considerably reducing the overall computational cost. A variety of test cases, ranging from inviscid flowsmore » to DNS of turbulent flows, are presented to assess the performance of these schemes. Here, numerical experiments demonstrate that the third-order ROW scheme for the DAEs of index-2 can not only achieve the designed formal order of temporal convergence accuracy in a benchmark test, but also require significantly less computing time than its ESDIRK3 counterpart to converge to the same level of discretization errors in all of the flow simulations in this study, indicating that the ROW methods provide an attractive alternative for the higher-order time-accurate integration of the unsteady compressible Navier-Stokes equations.« less
  • The performance of a nonstationary iterative scheme for solving the large system of linear equations arising in neutron diffusion calculations is reported. The method is compared with successive overrelaxation and conjugate gradient methods. Practical aspects of its implementation on computers are considered.