Wall modeling via function enrichment within a high‐order DG method for RANS simulations of incompressible flow

Summary We present a novel approach to wall modeling for RANS within the discontinuous Galerkin method. Wall functions are not used to prescribe boundary conditions as usual but they are built into the function space of the numerical method as a local enrichment, in addition to the standard polynomial component. The Galerkin method then automatically finds the optimal solution among all shape functions available. This idea is fully consistent and gives the wall model vast flexibility in separated boundary layers or high adverse pressure gradients. The wall model is implemented in a high-order discontinuous Galerkin solver for incompressible flow complemented by the Spalart–Allmaras closure model. As benchmark examples we present turbulent channel flow starting from Reτ=180 and up to Reτ=100,000 as well as flow past periodic hills at Reynolds numbers based on the hill height of ReH=10,595 and ReH=19,000. This article is protected by copyright. All rights reserved.

[1]  Luca Heltai,et al.  The deal . II Library , Version 8 . 4 . 0 , 2016 .

[2]  Javier Jiménez,et al.  Spectra of the very large anisotropic scales in turbulent channels , 2003 .

[3]  John Kim,et al.  DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOWS UP TO RE=590 , 1999 .

[4]  Martin Kronbichler,et al.  A high-order semi-explicit discontinuous Galerkin solver for 3D incompressible flow with application to DNS and LES of turbulent channel flow , 2016, J. Comput. Phys..

[5]  Ralf Hartmann,et al.  Adjoint Consistency Analysis of Discontinuous Galerkin Discretizations , 2007, SIAM J. Numer. Anal..

[6]  Paul Steinmann,et al.  On the $h$-adaptive PUM and the $hp$-adaptive FEM approaches applied to PDEs in quantum mechanics , 2016, 1612.02305.

[7]  S. Orszag,et al.  High-order splitting methods for the incompressible Navier-Stokes equations , 1991 .

[8]  K. Hanjalic,et al.  Compound Wall Treatment for RANS Computation of Complex Turbulent Flows and Heat Transfer , 2007 .

[9]  David Wells,et al.  The deal.II Library, Version 8.4 , 2016, J. Num. Math..

[10]  D. Arnold.  An Interior Penalty Finite Element Method with Discontinuous Elements , 1982 .

[11]  R. Moser,et al.  Direct numerical simulation of turbulent channel flow up to $\mathit{Re}_{{\it\tau}}\approx 5200$ , 2014, Journal of Fluid Mechanics.

[12]  Koen Hillewaert,et al.  Development of the discontinuous Galerkin method for high-resolution, large scale CFD and acoustics in industrial geometries , 2013 .

[13]  B. Launder,et al.  The numerical computation of turbulent flows , 1990 .

[14]  P. Spalart.  Comments on the feasibility of LES for wings, and on a hybrid RANS/LES approach , 1997 .

[15]  Jaime Peraire,et al.  Navier-Stokes Solution Using Hybridizable Discontinuous Galerkin methods , 2011 .

[16]  R. B. Dean.  Reynolds Number Dependence of Skin Friction and Other Bulk Flow Variables in Two-Dimensional Rectangular Duct Flow , 1978 .

[17]  B. Rivière,et al.  Improved energy estimates for interior penalty, constrained and discontinuous Galerkin methods for elliptic problems. Part I , 1999 .

[18]  R. Simpson.  A model for the backflow mean velocity profile , 1983 .

[19]  Charbel Farhat,et al.  A discontinuous enrichment method for the finite element solution of high Péclet advection-diffusion problems , 2009 .

[20]  R. B. Dean.  A Single Formula for the Complete Velocity Profile in a Turbulent Boundary Layer , 1976 .

[21]  Jochen Fröhlich,et al.  Highly resolved large-eddy simulation of separated flow in a channel with streamwise periodic constrictions , 2005, Journal of Fluid Mechanics.

[22]  Sj Steven Hulshoff,et al.  POD enriched boundary models and their optimal stabilisation , 2015 .

[23]  Michael Manhart,et al.  Flow over periodic hills: an experimental study , 2011 .

[24]  T. Belytschko,et al.  The extended/generalized finite element method: An overview of the method and its applications , 2010 .

[25]  Paul A. Durbin,et al.  Limiters and wall treatments in applied turbulence modeling , 2009 .

[26]  Valerio D'Alessandro,et al.  A Spalart-Allmaras turbulence model implementation in a discontinuous Galerkin solver for incompressible flows , 2013, J. Comput. Phys..

[27]  Suad Jakirlić,et al.  Extending the bounds of ‘steady’ RANS closures: Toward an instability-sensitive Reynolds stress model , 2015 .

[28]  Yan Chao,et al.  A Spalart–allmaras Turbulence Model Implementation for High-order Discontinuous Galerkin Solution of the Reynolds-averaged Navier-stokes Equations , 2016 .

[29]  Tobias Knopp,et al.  A grid and flow adaptive wall-function method for RANS turbulence modelling , 2006, J. Comput. Phys..

[30]  C. Farhat,et al.  The Discontinuous Enrichment Method , 2000 .

[31]  I. Babuska,et al.  The partition of unity finite element method: Basic theory and applications , 1996 .

[32]  D. Turner,et al.  ANALYTIC ENRICHMENT OF FINITE ELEMENT FORMULATIONS FOR CAPTURING BOUNDARY LAYER BEHAVIOR , 2010 .

[33]  P. Spalart.  A One-Equation Turbulence Model for Aerodynamic Flows , 1992 .

[34]  D. Mavriplis,et al.  Robust Computation of Turbulent Flows Using a Discontinuous Galerkin Method , 2012 .

[35]  Katharina Kormann,et al.  A generic interface for parallel cell-based finite element operator application , 2012 .

[36]  G. Iaccarino,et al.  Near-wall behavior of RANS turbulence models and implications for wall functions , 2005 .

[37]  Javier Jiménez,et al.  Scaling of the velocity fluctuations in turbulent channels up to Reτ=2003 , 2006 .

[38]  Ted Belytschko,et al.  Elastic crack growth in finite elements with minimal remeshing , 1999 .

[39]  D. Spalding.  A Single Formula for the “Law of the Wall” , 1961 .

[40]  Wolfgang A. Wall,et al.  A new approach to wall modeling in LES of incompressible flow via function enrichment , 2015, J. Comput. Phys..

[41]  Erik Burman,et al.  Numerical Approximation of Large Contrast Problems with the Unfitted Nitsche Method , 2011 .