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 »
- Authors:
-
- Univ. of California, Davis, CA (United States)
- 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}
}
Web of Science
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
On the Finite-Time Splash and Splat Singularities for the 3-D Free-Surface Euler Equations
journal, December 2013
- Coutand, Daniel; Shkoller, Steve
- Communications in Mathematical Physics, Vol. 325, Issue 1
The numerical simulation of two-dimensional fluid flow with strong shocks
journal, April 1984
- Woodward, Paul; Colella, Phillip
- Journal of Computational Physics, Vol. 54, Issue 1
The method of fractional steps for conservation laws
journal, September 1980
- Crandall, Michael; Majda, Andrew
- Numerische Mathematik, Vol. 34, Issue 3
Weighted Essentially Non-oscillatory Schemes on Triangular Meshes
journal, March 1999
- Hu, Changqing; Shu, Chi-Wang
- Journal of Computational Physics, Vol. 150, Issue 1
Efficient Implementation of Weighted ENO Schemes
journal, June 1996
- Jiang, Guang-Shan; Shu, Chi-Wang
- Journal of Computational Physics, Vol. 126, Issue 1
Nonoscillatory Central Schemes for Multidimensional Hyperbolic Conservation Laws
journal, November 1998
- Jiang, Guang-Shan; Tadmor, Eitan
- SIAM Journal on Scientific Computing, Vol. 19, Issue 6
Theory of the Rayleigh-Taylor instability
journal, August 1991
- Kull, H. J.
- Physics Reports, Vol. 206, Issue 5
Comparison of Several Difference Schemes on 1D and 2D Test Problems for the Euler Equations
journal, January 2003
- Liska, Richard; Wendroff, Burton
- SIAM Journal on Scientific Computing, Vol. 25, Issue 3
Compressibility effects on the Rayleigh–Taylor instability growth between immiscible fluids
journal, January 2004
- Livescu, D.
- Physics of Fluids, Vol. 16, Issue 1
Errors for calculations of strong shocks using an artificial viscosity and an artificial heat flux
journal, September 1987
- Noh, W. F.
- Journal of Computational Physics, Vol. 72, Issue 1
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
Investigation of the Character of the Equilibrium of an Incompressible Heavy Fluid of Variable Density
journal, November 1882
- Rayleigh,
- Proceedings of the London Mathematical Society, Vol. s1-14, Issue 1
A space–time smooth artificial viscosity method for nonlinear conservation laws
journal, February 2013
- Reisner, J.; Serencsa, J.; Shkoller, S.
- Journal of Computational Physics, Vol. 235
Revisiting Wall Heating
journal, August 2000
- Rider, William J.
- Journal of Computational Physics, Vol. 162, Issue 2
An overview of Rayleigh-Taylor instability
journal, July 1984
- Sharp, D. H.
- Physica D: Nonlinear Phenomena, Vol. 12, Issue 1-3
Resolution of high order WENO schemes for complicated flow structures
journal, April 2003
- Shi, Jing; Zhang, Yong-Tao; Shu, Chi-Wang
- Journal of Computational Physics, Vol. 186, Issue 2
High-order Finite Difference and Finite Volume WENO Schemes and Discontinuous Galerkin Methods for CFD
journal, March 2003
- Shu, Chi-Wang
- International Journal of Computational Fluid Dynamics, Vol. 17, Issue 2
Formation of singularities in three-dimensional compressible fluids
journal, December 1985
- Sideris, Thomas C.
- Communications in Mathematical Physics, Vol. 101, Issue 4
A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws
journal, April 1978
- Sod, Gary A.
- Journal of Computational Physics, Vol. 27, Issue 1
Time dependent boundary conditions for hyperbolic systems
journal, January 1987
- Thompson, Kevin W.
- Journal of Computational Physics, Vol. 68, Issue 1
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
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