Technical Brief

The Role of Forcing and Eddy Viscosity Variation on the Log-Layer Mismatch Observed in Wall-Modeled Large-Eddy Simulations

[+] Author and Article Information
Rey DeLeon

Department of Mechanical Engineering,
University of Idaho Moscow,
Moscow, ID 83844
e-mail: anthony.rey.deleon@gmail.com

Inanc Senocak

Department of Mechanical Engineering and
Materials Science,
University of Pittsburgh,
Pittsburgh, PA 15261
e-mail: senocak@pitt.edu

1Present address: ANSYS, Inc. Park City, UT 84098.

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received March 5, 2018; final manuscript received September 24, 2018; published online November 16, 2018. Assoc. Editor: Moran Wang.

J. Fluids Eng 141(5), 054501 (Nov 16, 2018) (6 pages) Paper No: FE-18-1151; doi: 10.1115/1.4041562 History: Received March 05, 2018; Revised September 24, 2018

We investigate the role of eddy viscosity variation and the effect of zonal enforcement of the mass flow rate on the log-layer mismatch problem observed in turbulent channel flows. An analysis of the mean momentum balance shows that it lacks a degree-of-freedom (DOF) when eddy viscosity is large, and the mean velocity conforms to an incorrect profile. Zonal enforcement of the target flow rate introduces an additional degree-of-freedom to the mean momentum balance, similar to an external stochastic forcing term, leading to a significant reduction in the log-layer mismatch. We simulate turbulent channel flows at friction Reynolds numbers of 2000 and 5200 on coarse meshes that do not resolve the viscous sublayer. The second-order turbulence statistics agree well with the direct numerical simulation benchmark data when results are normalized by the velocity scale extracted from the filtered velocity field. Zonal enforcement of the flow rate also led to significant improvements in skin friction coefficients.

Copyright © 2019 by ASME
Your Session has timed out. Please sign back in to continue.


Nikitin, N. , Nicoud, F. , Wasistho, B. , Squires, K. , and Spalart, P. , 2000, “ An Approach to Wall Modeling in Large-Eddy Simulations,” Phys. Fluids, 12(7), pp. 1629–1632. [CrossRef]
Germano, M. , 2004, “ Properties of the Hybrid RANS/LES Filter,” Theor. Comput. Fluid Dyn., 17(4), pp. 225–231. [CrossRef]
Hamba, F. , 2006, “ A Hybrid RANS/LES Simulation of High-Reynolds-Number Channel Flow Using Additional Filtering at the Interface,” Theor. Comput. Fluid Dyn., 20(2), pp. 89–101. [CrossRef]
Keating, A. , and Piomelli, U. , 2006, “ A Dynamic Stochastic Forcing Method as a Wall-Layer Model for Large-Eddy Simulation,” J. Turbul., 7(12), pp. 1–24.
Davidson, L. , and Billson, M. , 2006, “ Hybrid LES-RANS Using Synthesized Turbulent Fluctuations for Forcing in the Interface Region,” Int. J. Heat Fluid Flow, 27(6), pp. 1028–1042. [CrossRef]
Larsson, J. , Lien, F. , and Yee, E. , 2007, “ The Artificial Buffer Layer and the Effects of Forcing in Hybrid LES-RANS,” Int. J. Heat Fluid Flow, 28(6), pp. 1443–1459. [CrossRef]
Davidson, L. , 2016, “ Zonal PANS: Evaluation of Different Treatments of the Rans-Les Interface,” J. Turbul., 17(3), pp. 274–307. [CrossRef]
Yang, X. I. A. , Park, G. I. , and Moin, P. , 2017, “ Log-Layer Mismatch and Modeling of the Fluctuating Wall Stress in Wall-Modeled Large-Eddy Simulations,” Phys. Rev. Fluids, 2, p. 104601. [CrossRef]
Baggett, J. S. , 1998, “ On the Feasibility of Merging LES with RANS for the Near-wall Region of Attached Turbulent Flows,” Annual Research Briefs, Center for Turbulence Research, Stansford, CA, pp. 267–277.
Piomelli, U. , Balaras, E. , Pasinato, H. , Squires, K. D. , and Spalart, P. R. , 2003, “ The Inner–Outer Layer Interface in Large-Eddy Simulations With Wall-Layer Models,” Int. J. Heat Fluid Flow, 24(4), pp. 538–550. [CrossRef]
Meneveau, C. , Lund, T. , and Cabot, W. , 1996, “ A Lagrangian Dynamic Subgrid-Scale Model of Turbulence,” J. Fluid Mech., 319(1), pp. 353–85. [CrossRef]
Schumann, U. , 1975, “ Subgrid Scale Model for Finite Difference Simulations of Turbulent Flows in Plane Channels and Annuli,” J. Comput. Phys., 18(4), pp. 376–404. [CrossRef]
Grötzbach, G. , 1987, “ Direct Numerical and Large Eddy Simulation of Turbulent Channel Flows,” Encyclopedia Fluid Mechanics, Vol. 6, N. P. Cheremisinoff, ed., Gulf Publishing, West Orange, NJ.
Lee, J. , Cho, M. , and Choi, H. , 2013, “ Large Eddy Simulations of Turbulent Channel and Boundary Layer Flows at High Reynolds Number With Mean Wall Shear Stress Boundary Condition,” Phys. Fluids, 25(11), p. 110808. [CrossRef]
Cabot, W. , and Moin, P. , 2000, “ Approximate Wall Boundary Conditions in the Large-Eddy Simulation of High Reynolds Number Flow,” Flow Turbul. Combust., 63(1–4), pp. 269–291. [CrossRef]
Jacobsen, D. , and Senocak, I. , 2011, “ A Full-Depth Amalgamated Parallel 3D Geometric Multigrid Solver for GPU Clusters,” AIAA Paper No. 2011-946.
Jacobsen, D. A. , and Senocak, I. , 2013, “ Multi-Level Parallelism for Incompressible Flow Computations on GPU Clusters,” Parallel Comput., 39(1), pp. 1–20. [CrossRef]
Benocci, C. , and Pinelli, A. , 1990, “ The Role of the Forcing Term in the Large Eddy Simulation of Equilibrium Channel Flow,” Engineering Turbulence Modeling and Experiments, W. Rodi and E. Ganic , eds., Elsevier, New York, pp. 287–296.
Lee, M. , and Moser, R. , 2015, “ Direct Numerical Simulation of Turbulent Channel Flow Up to Reτ ≈ 5200,” J. Fluid Mech., 774(7), pp. 395–415. [CrossRef]
Pope, S. B. , 2000, Turbulent Flows, Cambridge University Press, Cambridge, UK.
Hoyas, S. , and Jimenez, J. , 2006, “ Scaling of the Velocity Fluctuations in Turbulent Channels Up to Reτ = 2003,” Phys. Fluids, 18(1), p. 011702. [CrossRef]
Dean, R. , 1978, “ Reynolds Number Dependence of Skin Friction and Other Bulk Flow Variables in Two-Dimensional Rectangular Duct Flow,” ASME J. Fluids Eng., 100(2), pp. 215–223. [CrossRef]
Tuerke, F. , and Jiménez, J. , 2013, “ Simulations of Turbulent Channels With Prescribed Velocity Profiles,” J. Fluid Mech., 723, pp. 587–603. [CrossRef]


Grahic Jump Location
Fig. 1

Split forcing demarcations. (a) Sketch of three forcing regions in a channel flow. (b) The normalized mean eddy viscosity profile obtained from a RANS-LES simulation of Reτ = 2000 turbulent channel flow. The dashed line represents the height of the RANS-LES interface, hRL, and the first split forcing demarcation. The dotted line represents the second split forcing demarcation, hsp.

Grahic Jump Location
Fig. 2

Nondimensional mean velocity profiles. (a) Log-layer mismatch reproduced using the single forcing approach. Reτ = 2000. (b) Reτ = 5200 (shifted up by 3.0) and Reτ = 2000 simulated with the tri-split approach. All results are normalized by uτ,w obtained from the wall shear stress.

Grahic Jump Location
Fig. 3

Attainment of a constant velocity scale in the LES region. Profiles of local friction velocity normalized by the friction velocity at the wall are shown. Both Reτ = 2000 and Reτ = 5200 are simulated with the tri-split forcing approach. Dashed vertical lines, hRL; dotted vertical lines, hsp.

Grahic Jump Location
Fig. 4

Profiles of RMS velocity fluctuations and the effect of velocity scale choice on the results. (a) Reτ = 2000 and (b) Reτ = 5200. Streamwise and spanwise components shifted up by 2.0 and 1.0, respectively. Dashed vertical lines, hRL; dotted vertical lines, hsp.

Grahic Jump Location
Fig. 5

Visualization of instantaneous streamwise velocity for Reτ = 2000 at y/δ ≈ 0.13. (a) Single forcing and (b) tri-split forcing. All velocities and lengths are normalized by the bulk velocity and channel half height, respectively.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In