Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

Estimating uncertainties in statistics computed from direct numerical simulation

Journal Article · · Physics of Fluids
DOI:https://doi.org/10.1063/1.4866813· OSTI ID:2565029

Rigorous assessment of uncertainty is crucial to the utility of direct numerical simulation (DNS) results. Uncertainties in the computed statistics arise from two sources: finite statistical sampling and the discretization of the Navier–Stokes equations. Due to the presence of non-trivial sampling error, standard techniques for estimating discretization error (such as Richardson extrapolation) fail or are unreliable. This work provides a systematic and unified approach for estimating these errors. First, a sampling error estimator that accounts for correlation in the input data is developed. Then, this sampling error estimate is used as part of a Bayesian extension of Richardson extrapolation in order to characterize the discretization error. These methods are tested using the Lorenz equations and are shown to perform well. These techniques are then used to investigate the sampling and discretization errors in the DNS of a wall-bounded turbulent flow at Reτ ≈ 180. Both small (Lx/δ × Lz/δ = 4π × 2π) and large (Lx/δ × Lz/δ = 12π × 4π) domain sizes are investigated. For each case, a sequence of meshes was generated by first designing a “nominal” mesh using standard heuristics for wall-bounded simulations. These nominal meshes were then coarsened to generate a sequence of grid resolutions appropriate for the Bayesian Richardson extrapolation method. In addition, the small box case is computationally inexpensive enough to allow simulation on a finer mesh, enabling the results of the extrapolation to be validated in a weak sense. For both cases, it is found that while the sampling uncertainty is large enough to make the order of accuracy difficult to determine, the estimated discretization errors are quite small. This indicates that the commonly used heuristics provide adequate resolution for this class of problems. However, it is also found that, for some quantities, the discretization error is not small relative to sampling error, indicating that the conventional wisdom that sampling error dominates discretization error for this class of simulations needs to be reevaluated.

Sponsoring Organization:
USDOE
Grant/Contract Number:
FC52-08NA28615
OSTI ID:
2565029
Journal Information:
Physics of Fluids, Journal Name: Physics of Fluids Journal Issue: 3 Vol. 26; ISSN 1070-6631
Publisher:
American Institute of PhysicsCopyright Statement
Country of Publication:
United States
Language:
English

References (28)

A Critical Evaluation of the Resolution Properties of B-Spline and Compact Finite Difference Methods journal December 2001
�ber die partiellen Differenzengleichungen der mathematischen Physik journal December 1928
Spectral methods for the Navier-Stokes equations with one infinite and two periodic directions journal October 1991
Higher order B-spline collocation at the Greville abscissae journal January 2005
Review of code and solution verification procedures for computational simulation journal May 2005
Bayesian uncertainty analysis with applications to turbulence modeling journal September 2011
Scaling of the energy spectra of turbulent channels journal January 1999
Turbulence statistics in fully developed channel flow at low Reynolds number journal April 1987
Spectra of the very large anisotropic scales in turbulent channels journal January 2003
Scaling of the velocity fluctuations in turbulent channels up to Reτ=2003 journal January 2006
Dissipation and enstrophy in isotropic turbulence: Resolution effects and scaling in direct numerical simulations journal April 2008
Reynolds number effects on the Reynolds-stress budgets in turbulent channels journal October 2008
Three Curious Properties of the Sample Variance and Autocovariance for Stationary Processes with Unknown Mean journal November 1993
B-spline Methods in Fluid Dynamics journal March 2003
emcee : The MCMC Hammer
  • Foreman-Mackey, Daniel; Hogg, David W.; Lang, Dustin
  • Publications of the Astronomical Society of the Pacific, Vol. 125, Issue 925 https://doi.org/10.1086/670067
journal March 2013
Bayesian uncertainty quantification applied to RANS turbulence models journal December 2011
What are we learning from simulating wall turbulence?
  • Jiménez, Javier; Moser, Robert D.
  • Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, Vol. 365, Issue 1852 https://doi.org/10.1098/rsta.2006.1943
journal January 2007
Automatic spectral analysis with time series models journal April 2002
Finite sample criteria for autoregressive order selection journal January 2000
Comments on the performance of maximum entropy algorithms journal November 1978
Commentary on the denominator recursion for Burg's block algorithm journal January 1986
Pure order recursive least-squares ladder algorithms journal August 1986
Practical Aspects of the Spectral Analysis of Irregularly Sampled Data With Time-Series Models journal May 2009
On Direct Numerical Simulation of Turbulent Flows journal March 2011
Petascale direct numerical simulation of turbulent channel flow on up to 786K cores
  • Lee, Myoungkyu; Malaya, Nicholas; Moser, Robert D.
  • Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis on - SC '13 https://doi.org/10.1145/2503210.2503298
conference January 2013
DIRECT NUMERICAL SIMULATION: A Tool in Turbulence Research journal January 1998
On the Calculation of Filter Coefficients for Maximum Entropy Spectral Analysis journal February 1974
Ensemble samplers with affine invariance journal January 2010

Similar Records

Comparison of Full-Field and Integrated CFD Convergence Based on Richardson Extrapolation
Conference · Thu Aug 01 00:00:00 EDT 2024 · OSTI ID:2438863

A Method for Treating Discretization Error in Nondeterministic Analysis
Conference · Tue Jan 26 23:00:00 EST 1999 · OSTI ID:3300

Quantifying the uncertainty introduced by discretization and time-averaging in two-fluid model predictions
Journal Article · Wed Jul 12 00:00:00 EDT 2017 · AIChE Journal · OSTI ID:1440343

Related Subjects