A mixed/hybrid Dirichlet formulation for wall-bounded flow problems including turbulent flow

Abstract In this study, a new method for the weak imposition of Dirichlet boundary conditions in convection-dominated convection–diffusion and wall-bounded turbulent flow problems is developed. It is based on an embedded Dirichlet formulation which utilises an additional discontinuous stress field for imposing boundary conditions. Values at outflow boundaries may not be sufficiently controllable on underresolved meshes, e.g., when using Nitsche’s method for weak enforcement of Dirichlet boundary conditions. Hence, a consistent modification of an earlier-proposed constant-parameter embedded Dirichlet formulation that allows to regain control if required is introduced. A separate treatment of boundary conditions in wall-normal and tangential direction is enabled. This split is physically motivated, and it turns out to be crucial for a successful application of the method. On the one hand, it allows for including a formulation based on the law of the wall in tangential direction. On the other hand, a consistent modification derived from theoretical considerations for convection–diffusion problems is applied in wall-normal direction, preventing unphysical penetration of walls. Numerical results obtained for turbulent channel flow at modest friction Reynolds number Re τ = 395 employing meshes with and without adequate resolution of the boundary layer are presented. Superior performance of the proposed method compared to strongly-imposed boundary conditions is demonstrated by an error reduction of up to 50% for mean velocity results in the core of the channel on coarse discretisations.

[1]  S. Mittal,et al.  Incompressible flow computations with stabilized bilinear and linear equal-order-interpolation velocity-pressure elements , 1992 .

[2]  T. Hughes,et al.  Isogeometric variational multiscale modeling of wall-bounded turbulent flows with weakly enforced boundary conditions on unstretched meshes , 2010 .

[3]  P. Moin,et al.  Approximate Wall Boundary Conditions in the Large-Eddy Simulation of High Reynolds Number Flow , 2000 .

[4]  U. Schumann Subgrid Scale Model for Finite Difference Simulations of Turbulent Flows in Plane Channels and Annuli , 1975 .

[5]  Peter Hansbo,et al.  A velocity pressure streamline diffusion finite element method for Navier-Stokes equations , 1990 .

[6]  Erik Burman,et al.  Stabilized finite element methods for the generalized Oseen problem , 2007 .

[7]  G. Hulbert,et al.  A generalized-α method for integrating the filtered Navier–Stokes equations with a stabilized finite element method , 2000 .

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

[9]  U. Piomelli,et al.  Wall-layer models for large-eddy simulations , 2008 .

[10]  Wolfgang A. Wall,et al.  Time-dependent subgrid scales in residual-based large eddy simulation of turbulent channel flow , 2010 .

[11]  M. Kronbichler,et al.  An algebraic variational multiscale-multigrid method for large eddy simulation of turbulent flow , 2010 .

[12]  Pavel B. Bochev,et al.  On the Finite Element Solution of the Pure Neumann Problem , 2005, SIAM Rev..

[13]  T. Hughes,et al.  Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows , 2007 .

[14]  Victor M. Calo,et al.  Weak Dirichlet Boundary Conditions for Wall-Bounded Turbulent Flows , 2007 .

[15]  Wolfgang A. Wall,et al.  An embedded Dirichlet formulation for 3D continua , 2010 .

[16]  L. Franca,et al.  On an Improved Unusual Stabilized Finite Element Method for theAdvective-Reactive-Diffusive Equation , 1999 .

[17]  Wolfgang A. Wall,et al.  Residual‐based variational multiscale methods for laminar, transitional and turbulent variable‐density flow at low Mach number , 2011 .

[18]  J. Deardorff A numerical study of three-dimensional turbulent channel flow at large Reynolds numbers , 1970, Journal of Fluid Mechanics.

[19]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[20]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[21]  Victor M. Calo,et al.  Residual-based multiscale turbulence modeling: Finite volume simulations of bypass transition , 2005 .

[22]  Ramon Codina,et al.  A symmetric method for weakly imposing Dirichlet boundary conditions in embedded finite element meshes , 2012 .

[23]  P. Sagaut Large Eddy Simulation for Incompressible Flows , 2001 .

[24]  Wolfgang A. Wall,et al.  An algebraic variational multiscale-multigrid method for large-eddy simulation: generalized-α time integration, Fourier analysis and application to turbulent flow past a square-section cylinder , 2011 .

[25]  J. Fröhlich,et al.  Hybrid LES/RANS methods for the simulation of turbulent flows , 2008 .

[26]  Peter Hansbo,et al.  A velocity-pressure streamline diffusion finite element method for the incompressible Navier-Stokes equation , 1990 .

[27]  U. Piomelli Wall-layer models for large-eddy simulations , 2008 .

[28]  R. Codina Stabilization of incompressibility and convection through orthogonal sub-scales in finite element methods , 2000 .

[29]  Wolfgang A. Wall,et al.  An algebraic variational multiscale-multigrid method for large-eddy simulation of turbulent variable-density flow at low Mach number , 2010, J. Comput. Phys..

[30]  J. Nitsche Über ein Variationsprinzip zur Lösung von Dirichlet-Problemen bei Verwendung von Teilräumen, die keinen Randbedingungen unterworfen sind , 1971 .

[31]  Thomas J. R. Hughes,et al.  Finite element modeling of blood flow in arteries , 1998 .

[32]  R. Codina,et al.  Time dependent subscales in the stabilized finite element approximation of incompressible flow problems , 2007 .

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

[34]  Thomas J. R. Hughes,et al.  Weak imposition of Dirichlet boundary conditions in fluid mechanics , 2007 .