DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A space-time smooth artificial viscosity method with wavelet noise indicator and shock collision scheme, Part 2: The 2-D case

Abstract

This is the second part to our companion paper. Herein, we generalize to two space dimensions the C-method developed in for adding localized, space-time smooth artificial viscosity to nonlinear systems of conservation laws that propagate shock waves, rarefaction waves, and contact discontinuities. For gas dynamics, the C-method couples the Euler equations to scalar reaction-diffusion equations, which we call C-equations, whose solutions serve as space-time smooth artificial viscosity indicators for shocks and contacts. We develop a high-order numerical algorithm for gas dynamics in 2-D which can accurately simulate the Rayleigh-Taylor (RT) instability with Kelvin-Helmholtz (KH) roll-up of the contact discontinuity, as well as shock collision and bounce-back. Solutions to our C-equations not only indicate the location of the shocks and contacts, but also track the geometry of the evolving fronts. This allows us to implement both directionally isotropic and anisotropic artificial viscosity schemes, the latter adding diffusion only in directions tangential to the evolving front. We additionally produce a novel shock collision indicator function, which naturally activates during shock collision, and then smoothly deactivates. Moreover, we implement a high-frequency 2-D wavelet-based noise detector together with an efficient and localized noise removal algorithm. To test the methodology, we use a highly simplifiedmore » WENO-based discretization scheme. We provide numerical results for some classical 2-D test problems, including the RT problem, the Noh problem, a circular explosion problem from the Liska & Wendroff [13] review paper, the Sedov blast wave problem, the double Mach 10 reflection test, and a shock-wall collision problem. In particular, we show that our artificial viscosity method can eliminate the wall-heating phenomenon for the Noh problem, and thereby produce an accurate, non-oscillatory solution, even though our simplified WENO-type scheme fails to run for this problem.« less

Authors:
 [1];  [2]; ORCiD logo [1]
  1. Univ. of California, Davis, CA (United States)
  2. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Publication Date:
Research Org.:
Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA), Office of Defense Nuclear Nonproliferation
OSTI Identifier:
1501807
Alternate Identifier(s):
OSTI ID: 1635993
Report Number(s):
LA-UR-18-24730
Journal ID: ISSN 0021-9991
Grant/Contract Number:  
89233218CNA000001
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 387; Journal Issue: C; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; Shocks

Citation Formats

Ramani, Raaghav, Reisner, Jon, and Shkoller, Steve. A space-time smooth artificial viscosity method with wavelet noise indicator and shock collision scheme, Part 2: The 2-D case. United States: N. p., 2019. Web. doi:10.1016/j.jcp.2019.02.048.
Ramani, Raaghav, Reisner, Jon, & Shkoller, Steve. A space-time smooth artificial viscosity method with wavelet noise indicator and shock collision scheme, Part 2: The 2-D case. United States. https://doi.org/10.1016/j.jcp.2019.02.048
Ramani, Raaghav, Reisner, Jon, and Shkoller, Steve. Thu . "A space-time smooth artificial viscosity method with wavelet noise indicator and shock collision scheme, Part 2: The 2-D case". United States. https://doi.org/10.1016/j.jcp.2019.02.048. https://www.osti.gov/servlets/purl/1501807.
@article{osti_1501807,
title = {A space-time smooth artificial viscosity method with wavelet noise indicator and shock collision scheme, Part 2: The 2-D case},
author = {Ramani, Raaghav and Reisner, Jon and Shkoller, Steve},
abstractNote = {This is the second part to our companion paper. Herein, we generalize to two space dimensions the C-method developed in for adding localized, space-time smooth artificial viscosity to nonlinear systems of conservation laws that propagate shock waves, rarefaction waves, and contact discontinuities. For gas dynamics, the C-method couples the Euler equations to scalar reaction-diffusion equations, which we call C-equations, whose solutions serve as space-time smooth artificial viscosity indicators for shocks and contacts. We develop a high-order numerical algorithm for gas dynamics in 2-D which can accurately simulate the Rayleigh-Taylor (RT) instability with Kelvin-Helmholtz (KH) roll-up of the contact discontinuity, as well as shock collision and bounce-back. Solutions to our C-equations not only indicate the location of the shocks and contacts, but also track the geometry of the evolving fronts. This allows us to implement both directionally isotropic and anisotropic artificial viscosity schemes, the latter adding diffusion only in directions tangential to the evolving front. We additionally produce a novel shock collision indicator function, which naturally activates during shock collision, and then smoothly deactivates. Moreover, we implement a high-frequency 2-D wavelet-based noise detector together with an efficient and localized noise removal algorithm. To test the methodology, we use a highly simplified WENO-based discretization scheme. We provide numerical results for some classical 2-D test problems, including the RT problem, the Noh problem, a circular explosion problem from the Liska & Wendroff [13] review paper, the Sedov blast wave problem, the double Mach 10 reflection test, and a shock-wall collision problem. In particular, we show that our artificial viscosity method can eliminate the wall-heating phenomenon for the Noh problem, and thereby produce an accurate, non-oscillatory solution, even though our simplified WENO-type scheme fails to run for this problem.},
doi = {10.1016/j.jcp.2019.02.048},
journal = {Journal of Computational Physics},
number = C,
volume = 387,
place = {United States},
year = {Thu Mar 07 00:00:00 EST 2019},
month = {Thu Mar 07 00:00:00 EST 2019}
}

Journal Article:

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

Save / Share:

Works referenced in this record:

On Validating an Astrophysical Simulation Code
journal, November 2002

  • Calder, A. C.; Fryxell, B.; Plewa, T.
  • The Astrophysical Journal Supplement Series, Vol. 143, Issue 1
  • DOI: 10.1086/342267

On the Finite-Time Splash and Splat Singularities for the 3-D Free-Surface Euler Equations
journal, December 2013


The numerical simulation of two-dimensional fluid flow with strong shocks
journal, April 1984


The method of fractional steps for conservation laws
journal, September 1980

  • Crandall, Michael; Majda, Andrew
  • Numerische Mathematik, Vol. 34, Issue 3
  • DOI: 10.1007/BF01396704

Weighted Essentially Non-oscillatory Schemes on Triangular Meshes
journal, March 1999

  • Hu, Changqing; Shu, Chi-Wang
  • Journal of Computational Physics, Vol. 150, Issue 1
  • DOI: 10.1006/jcph.1998.6165

Efficient Implementation of Weighted ENO Schemes
journal, June 1996

  • Jiang, Guang-Shan; Shu, Chi-Wang
  • Journal of Computational Physics, Vol. 126, Issue 1
  • DOI: 10.1006/jcph.1996.0130

Nonoscillatory Central Schemes for Multidimensional Hyperbolic Conservation Laws
journal, November 1998


Theory of the Rayleigh-Taylor instability
journal, August 1991


Comparison of Several Difference Schemes on 1D and 2D Test Problems for the Euler Equations
journal, January 2003


Compressibility effects on the Rayleigh–Taylor instability growth between immiscible fluids
journal, January 2004


Errors for calculations of strong shocks using an artificial viscosity and an artificial heat flux
journal, September 1987


A contribution to the great Riemann solver debate
journal, March 1994

  • Quirk, James J.
  • International Journal for Numerical Methods in Fluids, Vol. 18, Issue 6
  • DOI: 10.1002/fld.1650180603

Investigation of the Character of the Equilibrium of an Incompressible Heavy Fluid of Variable Density
journal, November 1882


A space–time smooth artificial viscosity method for nonlinear conservation laws
journal, February 2013


Revisiting Wall Heating
journal, August 2000


An overview of Rayleigh-Taylor instability
journal, July 1984


Resolution of high order WENO schemes for complicated flow structures
journal, April 2003


High-order Finite Difference and Finite Volume WENO Schemes and Discontinuous Galerkin Methods for CFD
journal, March 2003


Formation of singularities in three-dimensional compressible fluids
journal, December 1985

  • Sideris, Thomas C.
  • Communications in Mathematical Physics, Vol. 101, Issue 4
  • DOI: 10.1007/BF01210741

A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws
journal, April 1978


Time dependent boundary conditions for hyperbolic systems
journal, January 1987


Experimental study of Rayleigh–Taylor instability: Low Atwood number liquid systems with single-mode initial perturbations
journal, May 2001

  • Waddell, J. T.; Niederhaus, C. E.; Jacobs, J. W.
  • Physics of Fluids, Vol. 13, Issue 5
  • DOI: 10.1063/1.1359762

A Method for the Numerical Calculation of Hydrodynamic Shocks
journal, March 1950

  • VonNeumann, J.; Richtmyer, R. D.
  • Journal of Applied Physics, Vol. 21, Issue 3
  • DOI: 10.1063/1.1699639