Anisotropic adaptive stabilized finite element solver for RANS models

Aerodynamic characteristics of various geometries are predicted using a finite element formulation coupled with several numerical techniques to ensure stability and accuracy of the method. First, an edge based error estimator and anisotropic mesh adaptation are used to detect automatically all flow features under the constraint of a fixed number of elements, thus controlling the computational cost. A Variational MultiScale stabilized finite element method is employed to solve the incompressible Navier-Stokes equations. Finally, the Spalart-Allmaras turbulence model is solved using the Streamline Upwind Petrov-Galerkin (SUPG) method. This paper is meant to show that the combination of anisotropic unsteady mesh adaptation with stabilized finite element methods provides an adequate framework for solving turbulent flows at high Reynolds numbers. The proposed method was validated on several test cases by confrontation with literature of both numerical and experimental results, in terms of accuracy on the prediction of the drag and lift coefficients as well as their evolution in time for unsteady cases. This article is protected by copyright. All rights reserved.

[1]  Endre Süli,et al.  The Residual-Free-Bubble Finite Element Method on Anisotropic Partitions , 2007, SIAM J. Numer. Anal..

[2]  Thierry Coupez,et al.  Stabilized finite element method for incompressible flows with high Reynolds number , 2010, J. Comput. Phys..

[3]  Tayfun E. Tezduyar,et al.  Finite element stabilization parameters computed from element matrices and vectors , 2000 .

[4]  Simona Perotto,et al.  New anisotropic a priori error estimates , 2001, Numerische Mathematik.

[5]  T. Hughes,et al.  Isogeometric analysis : CAD, finite elements, NURBS, exact geometry and mesh refinement , 2005 .

[6]  Gerd Kunert,et al.  Edge residuals dominate a posteriori error estimates for linear finite element methods on anisotropic triangular and tetrahedral meshes , 2000, Numerische Mathematik.

[7]  Wagdi G. Habashi,et al.  Multiscale finite element method applied to the Spalart–Allmaras turbulence model for 3D detached-eddy simulation , 2012 .

[8]  Rolf Rannacher,et al.  Finite Element Methods for the Incompressible Navier-Stokes Equations , 2000 .

[9]  Ekkehard Ramm,et al.  Stabilized finite element formulation for incompressible flow on distorted meshes , 2009 .

[10]  Santiago Badia,et al.  Analysis of a Stabilized Finite Element Approximation of the Transient Convection-Diffusion Equation Using an ALE Framework , 2006, SIAM J. Numer. Anal..

[11]  Frédéric Hecht,et al.  New Progress in Anisotropic Grid Adaptation for Inviscid and Viscous Flows Simulations , 1995 .

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

[13]  Alvaro L. G. A. Coutinho,et al.  Edge‐based finite element implementation of the residual‐based variational multiscale method , 2009 .

[14]  Forrester T. Johnson,et al.  Modi cations and Clari cations for the Implementation of the Spalart-Allmaras Turbulence Model , 2011 .

[15]  Michael J. Hemsch,et al.  Statistical Analysis of Computational Fluid Dynamics Solutions from the Drag Prediction Workshop , 2004 .

[16]  Jürgen Fuhrmann,et al.  Guermond : " Theory and Practice of Finite Elements " , 2017 .

[17]  Frédéric Alauzet,et al.  Time accurate anisotropic goal-oriented mesh adaptation for unsteady flows , 2012, J. Comput. Phys..

[18]  L. Davidson,et al.  Large Eddy Simulation of Flow Past a Square Cylinder: Comparison of Different Subgrid Scale Models , 2000 .

[19]  J. Guermond,et al.  Theory and practice of finite elements , 2004 .

[20]  Pascal Frey,et al.  Anisotropic mesh adaptation for CFD computations , 2005 .

[21]  Antony Jameson,et al.  A perspective on computational algorithms for aerodynamic analysis and design , 2001 .

[22]  Volker John,et al.  Reference values for drag and lift of a two‐dimensional time‐dependent flow around a cylinder , 2004 .

[23]  Elie Hachem,et al.  On optimal simplicial 3D meshes for minimizing the Hessian-based errors , 2016 .

[24]  Thierry Coupez,et al.  Anisotropic meshing with time-stepping control for unsteady convection-dominated problems , 2015 .

[25]  Simona Perotto,et al.  Stabilized Finite Elements on Anisotropic Meshes: A Priori Error Estimates for the Advection-Diffusion and the Stokes Problems , 2003, SIAM J. Numer. Anal..

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

[27]  Kenji Shimada,et al.  APPLICATION OF A MODIFIED k–ε MODEL TO THE PREDICTION OF AERODYNAMIC CHARACTERISTICS OF RECTANGULAR CROSS-SECTION CYLINDERS , 2002 .

[28]  Nicola Parolini,et al.  Fluid-structure interaction problems in free surface flows: application to boat dynamics , 2008 .

[29]  E. D. Obasaju,et al.  An experimental study of pressure fluctuations on fixed and oscillating square-section cylinders , 1982, Journal of Fluid Mechanics.

[30]  Thomas J. R. Hughes,et al.  What are C and h ?: inequalities for the analysis and design of finite element methods , 1992 .

[31]  J. Remacle,et al.  Anisotropic adaptive simulation of transient flows using discontinuous Galerkin methods , 2005 .

[32]  R. Rannacher,et al.  Benchmark Computations of Laminar Flow Around a Cylinder , 1996 .

[33]  Yan Bao,et al.  Numerical prediction of aerodynamic characteristics of prismatic cylinder by finite element method with Spalart-Allmaras turbulence model , 2011 .

[34]  T. Coupez,et al.  Metric construction by length distribution tensor and edge based error for anisotropic adaptive meshing , 2011, J. Comput. Phys..

[35]  Simona Perotto,et al.  Anisotropic error estimates for elliptic problems , 2003, Numerische Mathematik.

[36]  Frédéric Alauzet,et al.  Adaptation de maillage anisotrope en trois dimensions : applications aux simulations instationnaires en mécanique des fluides , 2003 .

[37]  T. Hughes.  Multiscale phenomena: Green's functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of stabilized methods , 1995 .

[38]  Christopher L. Rumsey,et al.  Effective Inflow Conditions for Turbulence Models in Aerodynamic Calculations , 2007 .

[39]  Alfio Quarteroni,et al.  Mathematical models and numerical simulations for the America"s Cup , 2005 .

[40]  Sanjay Mittal,et al.  On the performance of high aspect ratio elements for incompressible flows , 2000 .