Wall‐modeled large‐eddy simulation in a finite element framework. (28th October 2019)

• Record Type: Journal Article
• Authors: Owen, Herbert; Chrysokentis, Georgios; Avila, Matias; Mira, Daniel; Houzeaux, Guillaume; Borrell, Ricard; Cajas, Juan Carlos; Lehmkuhl, Oriol
• Is Part Of: International journal for numerical methods in fluids. Volume 92:Number 1(2020)
• Page Start: 20
• Page End: 37
• Publication Date: 2019-10-28
• Subjects: exchange location -- finite elements -- large‐eddy simulation -- law of the wall -- wall modeling; Fluid dynamics -- Mathematics -- Periodicals; 532
• DOI: 10.1002/fld.4770
• Languages: English
• ISSNs: 0271-2091

Abstract

Summary: This work studies the implementation of wall modeling for large‐eddy simulation in a finite element context. It provides a detailed description of how the approach used by the finite volume and finite differences communities is adapted to the finite element context. The new implementation is as simple and easy to implement as the classical finite element one, but it provides vastly superior results. In the typical approach used in finite elements, the mesh does not extend all the way to the wall, and the wall stress is evaluated at the first grid point, based on the velocity at the same point. Instead, we adopt the approach commonly used in finite differences, where the mesh covers the whole domain and the wall stress is obtained at the wall grid point, with the velocity evaluated at the first grid point off the wall. The method is tested in a turbulent channel flow at R e τ =2003, a neutral atmospheric boundary layer flow, and a flow over a wall‐mounted hump, with significant improvement in the results compared to the standard finite element approach. Additionally, we examine the effect of evaluating the input velocity further away from the wall, as well as applying temporal filtering on the wall‐model input. Abstract : In the typical wall‐modeling approach used in finite elements, the mesh does not extend all the way to the wall, and the wall stress is evaluated at the first grid point, based on the velocity at the same point. Instead, we adopt the approach commonly used in finite differences, where the mesh covers the whole domain and the wall stress is obtained at the wall grid point, with the velocity evaluated at the first grid point off the wall.