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

Title: Scalar excursions in large-eddy simulations

Abstract

Here, the range of values of scalar fields in turbulent flows is bounded by their boundary values, for passive scalars, and by a combination of boundary values, reaction rates, phase changes, etc., for active scalars. The current investigation focuses on the local conservation of passive scalar concentration fields and the ability of the large-eddy simulation (LES) method to observe the boundedness of passive scalar concentrations. In practice, as a result of numerical artifacts, this fundamental constraint is often violated with scalars exhibiting unphysical excursions. The present study characterizes passive-scalar excursions in LES of a shear flow and examines methods for diagnosis and assesment of the problem. The analysis of scalar-excursion statistics provides support of the main hypothesis of the current study that unphysical scalar excursions in LES result from dispersive errors of the convection-term discretization where the subgrid-scale model (SGS) provides insufficient dissipation to produce a sufficiently smooth scalar field. In the LES runs three parameters are varied: the discretization of the convection terms, the SGS model, and grid resolution. Unphysical scalar excursions decrease as the order of accuracy of non-dissipative schemes is increased, but the improvement rate decreases with increasing order of accuracy. Two SGS models are examined, themore » stretched-vortex and a constant-coefficient Smagorinsky. Scalar excursions strongly depend on the SGS model. The excursions are significantly reduced when the characteristic SGS scale is set to double the grid spacing in runs with the stretched-vortex model. The maximum excursion and volume fraction of excursions outside boundary values show opposite trends with respect to resolution. The maximum unphysical excursions increase as resolution increases, whereas the volume fraction decreases. The reason for the increase in the maximum excursion is statistical and traceable to the number of grid points (sample size) which increases with resolution. In contrast, the volume fraction of unphysical excursions decreases with resolution because the SGS models explored perform better at higher grid resolution.« less

Authors:
 [1];  [1]
  1. California Inst. of Technology (CalTech), Pasadena, CA (United States)
Publication Date:
Research Org.:
California Institute of Technology (CalTech), Pasadena, CA (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA), Office of Defense Programs (DP); US Air Force Office of Scientific Research (AFOSR)
OSTI Identifier:
1343122
Alternate Identifier(s):
OSTI ID: 1359309
Grant/Contract Number:  
NA0002382; FA9550-12-1-0461
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 327; 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; large eddy simulation; numerical methods; numerical model error; scalar mixing; turbulent mixing

Citation Formats

Matheou, Georgios, and Dimotakis, Paul E. Scalar excursions in large-eddy simulations. United States: N. p., 2016. Web. doi:10.1016/j.jcp.2016.08.035.
Matheou, Georgios, & Dimotakis, Paul E. Scalar excursions in large-eddy simulations. United States. https://doi.org/10.1016/j.jcp.2016.08.035
Matheou, Georgios, and Dimotakis, Paul E. Wed . "Scalar excursions in large-eddy simulations". United States. https://doi.org/10.1016/j.jcp.2016.08.035. https://www.osti.gov/servlets/purl/1343122.
@article{osti_1343122,
title = {Scalar excursions in large-eddy simulations},
author = {Matheou, Georgios and Dimotakis, Paul E.},
abstractNote = {Here, the range of values of scalar fields in turbulent flows is bounded by their boundary values, for passive scalars, and by a combination of boundary values, reaction rates, phase changes, etc., for active scalars. The current investigation focuses on the local conservation of passive scalar concentration fields and the ability of the large-eddy simulation (LES) method to observe the boundedness of passive scalar concentrations. In practice, as a result of numerical artifacts, this fundamental constraint is often violated with scalars exhibiting unphysical excursions. The present study characterizes passive-scalar excursions in LES of a shear flow and examines methods for diagnosis and assesment of the problem. The analysis of scalar-excursion statistics provides support of the main hypothesis of the current study that unphysical scalar excursions in LES result from dispersive errors of the convection-term discretization where the subgrid-scale model (SGS) provides insufficient dissipation to produce a sufficiently smooth scalar field. In the LES runs three parameters are varied: the discretization of the convection terms, the SGS model, and grid resolution. Unphysical scalar excursions decrease as the order of accuracy of non-dissipative schemes is increased, but the improvement rate decreases with increasing order of accuracy. Two SGS models are examined, the stretched-vortex and a constant-coefficient Smagorinsky. Scalar excursions strongly depend on the SGS model. The excursions are significantly reduced when the characteristic SGS scale is set to double the grid spacing in runs with the stretched-vortex model. The maximum excursion and volume fraction of excursions outside boundary values show opposite trends with respect to resolution. The maximum unphysical excursions increase as resolution increases, whereas the volume fraction decreases. The reason for the increase in the maximum excursion is statistical and traceable to the number of grid points (sample size) which increases with resolution. In contrast, the volume fraction of unphysical excursions decreases with resolution because the SGS models explored perform better at higher grid resolution.},
doi = {10.1016/j.jcp.2016.08.035},
journal = {Journal of Computational Physics},
number = C,
volume = 327,
place = {United States},
year = {Wed Aug 31 00:00:00 EDT 2016},
month = {Wed Aug 31 00:00:00 EDT 2016}
}

Journal Article:

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

Save / Share:

Works referenced in this record:

An Intercomparison of Large-Eddy Simulations of the Stable Boundary Layer
journal, February 2006

  • Beare, Robert J.; Macvean, Malcolm K.; Holtslag, Albert A. M.
  • Boundary-Layer Meteorology, Vol. 118, Issue 2
  • DOI: 10.1007/s10546-004-2820-6

The Effects of Numerical Dissipation in Large Eddy Simulations
journal, October 2000


A further study of numerical errors in large-eddy simulations
journal, January 2003


Direct numerical simulation of stationary homogeneous stratified sheared turbulence
journal, March 2012


Large-Eddy Simulation of Stratified Turbulence. Part I: A Vortex-Based Subgrid-Scale Model
journal, May 2014

  • Chung, Daniel; Matheou, Georgios
  • Journal of the Atmospheric Sciences, Vol. 71, Issue 5
  • DOI: 10.1175/JAS-D-13-0126.1

The mixing transition in RayleighTaylor instability
journal, January 1999


The mixing transition in turbulent flows
journal, April 2000


Turbulent Mixing
journal, January 2005


Some consequences of the boundedness of scalar fluctuations
journal, November 1990

  • Dimotakis, Paul E.; Miller, Paul L.
  • Physics of Fluids A: Fluid Dynamics, Vol. 2, Issue 11
  • DOI: 10.1063/1.857666

Revisiting Numerical Errors in Direct and Large Eddy Simulations of Turbulence: Physical and Spectral Spaces Analysis
journal, December 2001

  • Fedioun, Ivan; Lardjane, Nicolas; Gökalp, Iskender
  • Journal of Computational Physics, Vol. 174, Issue 2
  • DOI: 10.1006/jcph.2001.6939

An Analysis of Numerical Errors in Large-Eddy Simulations of Turbulence
journal, April 1996


Numerical Calculation of Time-Dependent Viscous Incompressible Flow of Fluid with Free Surface
journal, January 1965

  • Harlow, Francis H.; Welch, J. Eddie
  • Physics of Fluids, Vol. 8, Issue 12
  • DOI: 10.1063/1.1761178

Hybrid tuned center-difference-WENO method for large eddy simulations in the presence of strong shocks
journal, March 2004


Comparison of conservative and rotational forms in large eddy simulation of turbulent channel flow
journal, August 1987


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

On the Effect of Numerical Errors in Large Eddy Simulations of Turbulent Flows
journal, March 1997


A stable and accurate convective modelling procedure based on quadratic upstream interpolation
journal, June 1979


New Trends in Large-Eddy Simulations of Turbulence
journal, January 1996


Weighted Essentially Non-oscillatory Schemes
journal, November 1994

  • Liu, Xu-Dong; Osher, Stanley; Chan, Tony
  • Journal of Computational Physics, Vol. 115, Issue 1
  • DOI: 10.1006/jcph.1994.1187

Strained spiral vortex model for turbulent fine structure
journal, January 1982


Large-eddy simulation of mixing in a recirculating shear flow
journal, March 2010

  • Matheou, Georgios; Bonanos, Aristides M.; Pantano, Carlos
  • Journal of Fluid Mechanics, Vol. 646
  • DOI: 10.1017/S0022112009992965

A recycling method for the large-eddy simulation of plumes in the atmospheric boundary layer
journal, May 2015


Large-Eddy Simulation of Stratified Turbulence. Part II: Application of the Stretched-Vortex Model to the Atmospheric Boundary Layer
journal, December 2014

  • Matheou, Georgios; Chung, Daniel
  • Journal of the Atmospheric Sciences, Vol. 71, Issue 12
  • DOI: 10.1175/JAS-D-13-0306.1

On the Fidelity of Large-Eddy Simulation of Shallow Precipitating Cumulus Convection
journal, September 2011

  • Matheou, Georgios; Chung, Daniel; Nuijens, Louise
  • Monthly Weather Review, Vol. 139, Issue 9
  • DOI: 10.1175/2011MWR3599.1

Large-eddy simulations of turbulent mixing layers using the stretched-vortex model
journal, February 2011


A comparison of second- and sixth-order methods for large-eddy simulations
journal, May 2002


Spectral large-eddy simulation of isotropic and stably stratified turbulence
journal, June 1992


A vortex-based subgrid stress model for large-eddy simulation
journal, August 1997

  • Misra, Ashish; Pullin, D. I.
  • Physics of Fluids, Vol. 9, Issue 8
  • DOI: 10.1063/1.869361

Suitability of Upwind-Biased Finite Difference Schemes for Large-Eddy Simulation of Turbulent Flows
journal, August 1997

  • Mittal, Rajat; Moin, Parviz
  • AIAA Journal, Vol. 35, Issue 8
  • DOI: 10.2514/2.253

Spectral Analysis of Large-Eddy Simulations of the Convective Boundary Layer
journal, December 1988


An introduction to SPH
journal, January 1988


Fully Conservative Higher Order Finite Difference Schemes for Incompressible Flow
journal, June 1998

  • Morinishi, Y.; Lund, T. S.; Vasilyev, O. V.
  • Journal of Computational Physics, Vol. 143, Issue 1
  • DOI: 10.1006/jcph.1998.5962

A high-resolution Petrov–Galerkin method for the 1D convection–diffusion–reaction problem
journal, January 2010

  • Nadukandi, Prashanth; Oñate, Eugenio; Garcia, Julio
  • Computer Methods in Applied Mechanics and Engineering, Vol. 199, Issue 9-12
  • DOI: 10.1016/j.cma.2009.10.009

On the Elimination of Aliasing in Finite-Difference Schemes by Filtering High-Wavenumber Components
journal, September 1971


LES approach for high Reynolds number wall-bounded flows with application to turbulent channel flow
journal, November 2008

  • Pantano, C.; Pullin, D. I.; Dimotakis, P. E.
  • Journal of Computational Physics, Vol. 227, Issue 21
  • DOI: 10.1016/j.jcp.2008.04.015

Discretization errors in large eddy simulation: on the suitability of centered and upwind-biased compact difference schemes
journal, August 2004

  • Park, Noma; Yoo, Jung Yul; Choi, Haecheon
  • Journal of Computational Physics, Vol. 198, Issue 2
  • DOI: 10.1016/j.jcp.2004.01.017

Ten questions concerning the large-eddy simulation of turbulent flows
journal, January 2004


A vortex-based model for the subgrid flux of a passive scalar
journal, September 2000


Reynolds stresses and one‐dimensional spectra for a vortex model of homogeneous anisotropic turbulence
journal, May 1994

  • Pullin, D. I.; Saffman, P. G.
  • Physics of Fluids, Vol. 6, Issue 5
  • DOI: 10.1063/1.868240

Solution of the Advective Equation by Upstream Interpolation with a Cubic Spline
journal, January 1976


High Accuracy Compact Schemes and Gibbs' Phenomenon
journal, December 2004

  • Sengupta, Tapan K.; Ganerwal, G.; Dipankar, Anurag
  • Journal of Scientific Computing, Vol. 21, Issue 3
  • DOI: 10.1007/s10915-004-1317-2

Quasi-Monotone Advection Schemes Based on Explicit Locally Adaptive Dissipation
journal, June 1998


MPDATA: An edge-based unstructured-grid formulation
journal, July 2005

  • Smolarkiewicz, Piotr K.; Szmelter, Joanna
  • Journal of Computational Physics, Vol. 206, Issue 2
  • DOI: 10.1016/j.jcp.2004.12.021

Iterated upwind schemes for gas dynamics
journal, January 2009

  • Smolarkiewicz, Piotr K.; Szmelter, Joanna
  • Journal of Computational Physics, Vol. 228, Issue 1
  • DOI: 10.1016/j.jcp.2008.08.008

Spectral methods for the Navier-Stokes equations with one infinite and two periodic directions
journal, October 1991

  • Spalart, Philippe R.; Moser, Robert D.; Rogers, Michael M.
  • Journal of Computational Physics, Vol. 96, Issue 2
  • DOI: 10.1016/0021-9991(91)90238-G

Semi-Lagrangian Integration Schemes for Atmospheric Models—A Review
journal, September 1991


An improved bounded semi-Lagrangian scheme for the turbulent transport of passive scalars
journal, September 2014


A physical-space version of the stretched-vortex subgrid-stress model for large-eddy simulation
journal, July 2000

  • Voelkl, Tobias; Pullin, D. I.; Chan, Daniel C.
  • Physics of Fluids, Vol. 12, Issue 7
  • DOI: 10.1063/1.870429

Comparision of Numerical Schemes in Large-Eddy Simulation of the Temporal Mixing Layer
journal, February 1996


Two-Dimensional Semi-Lagrangian Transport with Shape-Preserving Interpolation
journal, January 1989


A particle-gridless hybrid method for incompressible flows
journal, June 1999


On maximum-principle-satisfying high order schemes for scalar conservation laws
journal, May 2010


Works referencing / citing this record:

Turbulence Structure in a Stratocumulus Cloud
journal, October 2018


LES of Subsonic Reacting Mixing Layers
journal, December 2019

  • Kartha, Anand; Subbareddy, Pramod K.; Candler, Graham V.
  • Flow, Turbulence and Combustion, Vol. 104, Issue 4
  • DOI: 10.1007/s10494-019-00066-4

Turbulent shear-layer mixing: initial conditions, and direct-numerical and large-eddy simulations
journal, August 2019

  • Sharan, Nek; Matheou, Georgios; Dimotakis, Paul E.
  • Journal of Fluid Mechanics, Vol. 877
  • DOI: 10.1017/jfm.2019.591

Galilean invariance of shallow cumulus convection large-eddy simulations
journal, February 2021