skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Temporal slow-growth formulation for direct numerical simulation of compressible wall-bounded flows

; ; ;
Publication Date:
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
Grant/Contract Number:
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Physical Review Fluids
Additional Journal Information:
Journal Volume: 2; Journal Issue: 8; Related Information: CHORUS Timestamp: 2017-08-08 10:09:31; Journal ID: ISSN 2469-990X
American Physical Society
Country of Publication:
United States

Citation Formats

Topalian, Victor, Oliver, Todd A., Ulerich, Rhys, and Moser, Robert D.. Temporal slow-growth formulation for direct numerical simulation of compressible wall-bounded flows. United States: N. p., 2017. Web. doi:10.1103/PhysRevFluids.2.084602.
Topalian, Victor, Oliver, Todd A., Ulerich, Rhys, & Moser, Robert D.. Temporal slow-growth formulation for direct numerical simulation of compressible wall-bounded flows. United States. doi:10.1103/PhysRevFluids.2.084602.
Topalian, Victor, Oliver, Todd A., Ulerich, Rhys, and Moser, Robert D.. 2017. "Temporal slow-growth formulation for direct numerical simulation of compressible wall-bounded flows". United States. doi:10.1103/PhysRevFluids.2.084602.
title = {Temporal slow-growth formulation for direct numerical simulation of compressible wall-bounded flows},
author = {Topalian, Victor and Oliver, Todd A. and Ulerich, Rhys and Moser, Robert D.},
abstractNote = {},
doi = {10.1103/PhysRevFluids.2.084602},
journal = {Physical Review Fluids},
number = 8,
volume = 2,
place = {United States},
year = 2017,
month = 8

Journal Article:
Free Publicly Available Full Text
This content will become publicly available on August 7, 2018
Publisher's Accepted Manuscript

Save / Share:
  • A new class of compact-reconstruction weighted essentially non-oscillatory (CRWENO) schemes were introduced (Ghosh and Baeder in SIAM J Sci Comput 34(3): A1678–A1706, 2012) with high spectral resolution and essentially non-oscillatory behavior across discontinuities. The CRWENO schemes use solution-dependent weights to combine lower-order compact interpolation schemes and yield a high-order compact scheme for smooth solutions and a non-oscillatory compact scheme near discontinuities. The new schemes result in lower absolute errors, and improved resolution of discontinuities and smaller length scales, compared to the weighted essentially non-oscillatory (WENO) scheme of the same order of convergence. Several improvements to the smoothness-dependent weights, proposed inmore » the literature in the context of the WENO schemes, address the drawbacks of the original formulation. This paper explores these improvements in the context of the CRWENO schemes and compares the different formulations of the non-linear weights for flow problems with small length scales as well as discontinuities. Simplified one- and two-dimensional inviscid flow problems are solved to demonstrate the numerical properties of the CRWENO schemes and its different formulations. Canonical turbulent flow problems—the decay of isotropic turbulence and the shock-turbulence interaction—are solved to assess the performance of the schemes for the direct numerical simulation of compressible, turbulent flows« less
  • Computational analyses of fluid flow through packed pebble bed domains using the Reynolds-averaged NavierStokes framework have had limited success in the past. Because of a lack of high-fidelity experimental or computational data, optimization of Reynolds-averaged closure models for these geometries has not been extensively developed. In the present study, direct numerical simulation was employed to develop a high-fidelity database that can be used for optimizing Reynolds-averaged closure models for pebble bed flows. A face-centered cubic domain with periodic boundaries was used. Flow was simulated at a Reynolds number of 9308 and cross-verified by using available quasi-DNS data. During the simulations,more » low-frequency instability modes were observed that affected the stationary solution. Furthermore, these instabilities were investigated by using the method of proper orthogonal decomposition, and a correlation was found between the time-dependent asymmetry of the averaged velocity profile data and the behavior of the highest energy eigenmodes.« less
  • A B-spline based numerical method on a zonal embedded grid has been developed. The method is aimed at reducing the computational requirements for large eddy simulations (LES) and direct numerical simulations (DNS) of wall-bounded turbulent flows. The objective is to reduce the number of grid points required to resolve the near-wall eddies without placing a large number of grid points in the outer layers. DNS and LES calculations of a turbulent channel flow were performed on a grid with a zone near the wall that was refined in all three directions. The results from the zonal grid calculations show goodmore » agreement with previously published numerical and experimental results obtained for the same flow conditions. The zonal grid calculations required only a fraction of the CPU time required for the single zone grid calculation with the same near-wall grid density. In addition, the memory requirements for the zonal grid calculations are significantly reduced. 26 refs., 19 figs., 1 tab.« less
  • Capabilities and accuracy issues in Lagrangian tracking of heavy particles in velocity fields obtained from large-eddy simulations (LES) of wall-bounded turbulent flows are reviewed. In particular, it is shown that, if no subgrid scale (SGS) model is added to the particle motion equations, particle preferential concentration and near-wall accumulation are significantly underestimated. Results obtained with SGS modeling for the particle motion equations based on approximate deconvolution are briefly recalled. Then, the error purely due to filtering in particle tracking in LES flow fields is singled out and analyzed. The statistical properties of filtering errors are characterized in turbulent channel flowmore » both from an Eulerian and a Lagrangian viewpoint. Implications for stochastic SGS modeling in particle motion equations are briefly outlined.« less
  • The pressure Hessian tensor plays a key role in shaping the behavior of the velocity gradient tensor, and in turn, that of many incumbent non-linear processes in a turbulent flow field. In compressible flows, the role of pressure Hessian is even more important because it represents the level of fluid-thermodynamic coupling existing in the flow field. In this work, we first perform a direct numerical simulation-based study to clearly identify, isolate, and understand various important inviscid mechanisms that govern the evolution of the pressure Hessian tensor in compressible turbulence. The ensuing understanding is then employed to introduce major improvements tomore » the existing Lagrangian model of the pressure Hessian tensor (the enhanced Homogenized Euler equation or EHEE) in terms of (i) non-symmetric, non-isentropic effects and (ii) improved representation of the anisotropic portion of the pressure Hessian tensor. Finally, we evaluate the new model extensively by comparing the new model results against known turbulence behavior over a range of Reynolds and Mach numbers. Indeed, the new model shows much improved performance as compared to the EHEE model.« less