The effect of artificial bulk viscosity in simulations of forced compressible turbulence

Campos, A.; Morgan, B.

doi:10.1016/j.jcp.2018.05.030

Title: The effect of artificial bulk viscosity in simulations of forced compressible turbulence

Abstract

The use of an artificial bulk viscosity for shock stabilization is a common approach employed in turbulence simulations with high-order numerics. The effect of the artificial bulk viscosity is analyzed in the context of large eddy simulations by using as a test case simulations of linearly-forced compressible homogeneous turbulence (Petersen and Livescu, 2010 [12]). This case is unique in that it allows for the specification of a priori target values for total dissipation and ratio of solenoidal to dilatational dissipation. A comparison between these target values and the true predicted levels of dissipation is thus used to investigate the performance of the artificial bulk viscosity. Results show that the artificial bulk viscosity is effective at achieving stable solutions, but also leads to large values of artificial dissipation that outweigh the physical dissipation caused by fluid viscosity. An alternate approach, which employs the artificial thermal conductivity only, shows that the dissipation of dilatational modes is entirely due to the fluid viscosity. However, this method leads to unwanted Gibbs oscillations around the shocklets. The use of shock sensors that further localize the artificial bulk viscosity did not reduce the amount of artificial dissipation introduced by the artificial bulk viscosity. Finally, an improvedmore »« less

Authors:

^[1]; Morgan, B. ^[1]

Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)

Publication Date:: 2018-05-17

Research Org.:: Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)

Sponsoring Org.:: USDOE National Nuclear Security Administration (NNSA)

OSTI Identifier:: 1455410

Alternate Identifier(s):: OSTI ID: 1568906

Report Number(s):: LLNL-JRNL-739678
Journal ID: ISSN 0021-9991; 893182; TRN: US1901234

Grant/Contract Number:: AC52-07NA27344

Resource Type:: Accepted Manuscript

Journal Name:: Journal of Computational Physics

Additional Journal Information:: Journal Volume: 371; Journal Issue: C; Journal ID: ISSN 0021-9991

Publisher:: Elsevier

Country of Publication:: United States

Language:: English

Subject:: 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; 97 MATHEMATICS AND COMPUTING; arti cial viscosity; compressible turbulence; large-eddy simulation

Citation Formats


                    Campos, A., and Morgan, B. The effect of artificial bulk viscosity in simulations of forced compressible turbulence.  United States: N. p., 2018. 
Web.  https://doi.org/10.1016/j.jcp.2018.05.030.

Copy to clipboard


                    Campos, A., & Morgan, B. The effect of artificial bulk viscosity in simulations of forced compressible turbulence.  United States.  https://doi.org/10.1016/j.jcp.2018.05.030

Copy to clipboard


                    Campos, A., and Morgan, B. Thu .  
"The effect of artificial bulk viscosity in simulations of forced compressible turbulence".  United States.  https://doi.org/10.1016/j.jcp.2018.05.030.  https://www.osti.gov/servlets/purl/1455410.

Copy to clipboard


                    
@article{osti_1455410,

  title        = {The effect of artificial bulk viscosity in simulations of forced compressible turbulence},

  author       = {Campos, A. and Morgan, B.},

  abstractNote = {The use of an artificial bulk viscosity for shock stabilization is a common approach employed in turbulence simulations with high-order numerics. The effect of the artificial bulk viscosity is analyzed in the context of large eddy simulations by using as a test case simulations of linearly-forced compressible homogeneous turbulence (Petersen and Livescu, 2010 [12]). This case is unique in that it allows for the specification of a priori target values for total dissipation and ratio of solenoidal to dilatational dissipation. A comparison between these target values and the true predicted levels of dissipation is thus used to investigate the performance of the artificial bulk viscosity. Results show that the artificial bulk viscosity is effective at achieving stable solutions, but also leads to large values of artificial dissipation that outweigh the physical dissipation caused by fluid viscosity. An alternate approach, which employs the artificial thermal conductivity only, shows that the dissipation of dilatational modes is entirely due to the fluid viscosity. However, this method leads to unwanted Gibbs oscillations around the shocklets. The use of shock sensors that further localize the artificial bulk viscosity did not reduce the amount of artificial dissipation introduced by the artificial bulk viscosity. Finally, an improved forcing function that explicitly accounts for the role of the artificial bulk viscosity in the budget of turbulent kinetic energy was explored.},

  doi          = {10.1016/j.jcp.2018.05.030},

  journal      = {Journal of Computational Physics},

  number       = C,

  volume       = 371,

  place        = {United States},

  year         = {2018},

  month        = {5}

}

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.05.030

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 2 works

Citation information provided by
Web of Science

Figures / Tables:

Figure 1: Time history of turbulent dissipations for simulations that use artificial bulk viscosity only. Blue, $n$ = 64³; green, $n$ = 128³; red, $n$ = 256³; yellow, $n$ = 512₃. (For interpretation of the references to color in this figure, the reader is referred to the web version ofmore »

All figures and tables (9 total)

Save / Share:

Export Metadata

Save to My Library

Figures / Tables found in this record:

Figures/Tables have been extracted from DOE-funded journal article accepted manuscripts.

Similar Records in DOE PAGES and OSTI.GOV collections:

Assessment of high-resolution methods for numerical simulations of compressible turbulence with shock waves

Journal Article Johnsen, Eric ; Larsson, Johan ; Bhagatwala, Ankit V ; ... - Journal of Computational Physics

Flows in which shock waves and turbulence are present and interact dynamically occur in a wide range of applications, including inertial confinement fusion, supernovae explosion, and scramjet propulsion. Accurate simulations of such problems are challenging because of the contradictory requirements of numerical methods used to simulate turbulence, which must minimize any numerical dissipation that would otherwise overwhelm the small scales, and shock-capturing schemes, which introduce numerical dissipation to stabilize the solution. The objective of the present work is to evaluate the performance of several numerical methods capable of simultaneously handling turbulence and shock waves. A comprehensive range of high-resolution methodsmore »« less
https://doi.org/10.1016/j.jcp.2009.10.028
COMPARING NUMERICAL METHODS FOR ISOTHERMAL MAGNETIZED SUPERSONIC TURBULENCE

Journal Article Kritsuk, Alexei G ; Collins, David ; Norman, Michael L ; ... - Astrophysical Journal

Many astrophysical applications involve magnetized turbulent flows with shock waves. Ab initio star formation simulations require a robust representation of supersonic turbulence in molecular clouds on a wide range of scales imposing stringent demands on the quality of numerical algorithms. We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine popular astrophysical MHD methods actively used to model star formation. The set of nine codes includes: ENZO, FLASH, KT-MHD, LL-MHD, PLUTO, PPML, RAMSES, STAGGER, and ZEUS. These applications employ a variety of numerical approaches, including both split and unsplit,more »« less
https://doi.org/10.1088/0004-637X/737/1/13
Formulations of artificial viscosity for multi-dimensional shock wave computations

Journal Article Caramana, E J ; Shashkov, M J ; Whalen, P P - Journal of Computational Physics

In this paper the authors present a new formulation of the artificial viscosity concept. Physical arguments for the origins of this term are given and a set of criteria that any proper functional form of the artificial viscosity should satisfy is enumerated. The first important property is that by definition a viscosity must always be dissipative, transferring kinetic energy into internal energy, and must never act as a false pressure. The artificial viscous force should be Galilean invariant and vary continuously as a function of the criterion used to determine compression and expansion, and remain zero for the latter case.more »« less
https://doi.org/10.1006/jcph.1998.5989
TURBULENCE IN THE INTERGALACTIC MEDIUM: SOLENOIDAL AND DILATATIONAL MOTIONS AND THE IMPACT OF NUMERICAL VISCOSITY

Journal Article Zhu, Weishan ; Gu, Qiusheng ; Feng, Long-long ; ... - Astrophysical Journal

We use a suite of cosmological hydrodynamical simulations, run by two fixed grid codes, to investigate the properties of solenoidal and dilatational motions of the intergalactic medium (IGM) and the impact of numerical viscosity on turbulence in an ΛCDM universe. The codes differ only in the spatial difference discretization. We find that (1) The vortical motion grows rapidly since z = 2 and reaches ∼10 km s{sup –1}-90 km s{sup –1} at z = 0. Meanwhile, the small-scale compressive ratio r{sub CS} drops from 0.84 to 0.47, indicating comparable vortical and compressive motions at z = 0. (2) Power spectramore »« less
https://doi.org/10.1088/0004-637X/777/1/48
Entropy-based artificial viscosity stabilization for non-equilibrium Grey Radiation-Hydrodynamics

Journal Article Delchini, Marc O., E-mail: delchinm@email.tamu.edu ; Ragusa, Jean C., E-mail: jean.ragusa@tamu.edu ; Morel, Jim - Journal of Computational Physics

The entropy viscosity method is extended to the non-equilibrium Grey Radiation-Hydrodynamic equations. The method employs a viscous regularization to stabilize the numerical solution. The artificial viscosity coefficient is modulated by the entropy production and peaks at shock locations. The added dissipative terms are consistent with the entropy minimum principle. A new functional form of the entropy residual, suitable for the Radiation-Hydrodynamic equations, is derived. We demonstrate that the viscous regularization preserves the equilibrium diffusion limit. The equations are discretized with a standard Continuous Galerkin Finite Element Method and a fully implicit temporal integrator within the MOOSE multiphysics framework. The methodmore »« less
https://doi.org/10.1016/J.JCP.2015.04.039

Similar Records