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

doi:10.1016/j.jcp.2018.05.030

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

Full Record
Other Related Research

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.  doi: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.  doi: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.  doi: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

DOI: 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

Save / Share:

Export Metadata

Save to My Library

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
DOI: 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
DOI: 10.1088/0004-637X/737/1/13
Observation of a Kelvin-Helmholtz Instability in a High-Energy-Density Plasma on the Omega Laser

Journal Article Harding, E C ; Hansen, J F ; Hurricane, O A ; ... - Physical Review Letters, vol. 103, N/A, July 24, 2009, pp. 045005

A laser initiated experiment is described in which an unstable plasma shear layer is produced by driving a blast wave along a plastic surface with sinusoidal perturbations. In response to the vorticity deposited and the shear flow established by the blast wave, the interface rolls up into large vortices characteristic of the Kelvin-Helmholtz (KH) instability. The experiment used x ray radiography to capture the first well-resolved images of KH vortices in a high-energy-density plasma, and possibly the first images of transonic shocks generated by large-scale structures in a shear layer. The physical processes governing the evolution of a stratified fluidmore »« less
DOI: 10.1103/PhysRevLett.103.045005

Full Text Available
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
DOI: 10.1088/0004-637X/777/1/48
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
DOI: 10.1006/jcph.1998.5989

Similar Records