A Finite Volume approximation of the Navier-Stokes equations with nonlinear filtering stabilization

We consider a Leray model with a nonlinear differential low-pass filter for the simulation of incompressible fluid flow at moderately large Reynolds number (in the range of a few thousands) with under-refined meshes. For the implementation of the model, we adopt the three-step algorithm Evolve-Filter-Relax (EFR). The Leray model has been extensively applied within a Finite Element (FE) framework. Here, we propose to combine the EFR algorithm with a computationally efficient Finite Volume (FV) method. Our approach is validated against numerical data available in the literature for the 2D flow past a cylinder and against experimental measurements for the 3D fluid flow in an idealized medical device, as recommended by the U.S. Food and Drug Administration. We will show that for similar levels of mesh refinement FV and FE methods provide significantly different results. Through our numerical experiments, we are able to provide practical directions to tune the parameters involved in the model. Furthermore, we are able to investigate the impact of mesh features (element type, non-orthogonality, local refinement, and element aspect ratio) and the discretization method for the convective term on the agreement between numerical solutions and experimental data.

[1]  Steven Deutsch,et al.  Assessment of CFD Performance in Simulations of an Idealized Medical Device: Results of FDA’s First Computational Interlaboratory Study , 2012 .

[2]  O. Reynolds.  III. An experimental investigation of the circumstances which determine whether the motion of water shall be direct or sinuous, and of the law of resistance in parallel channels , 1883, Proceedings of the Royal Society of London.

[3]  A. Kolmogorov.  The local structure of turbulence in incompressible viscous fluid for very large Reynolds numbers , 1991, Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences.

[4]  C. Rhie,et al.  Numerical Study of the Turbulent Flow Past an Airfoil with Trailing Edge Separation , 1983 .

[5]  V. M. Tikhomirov,et al.  Dissipation of Energy in Isotropic Turbulence , 1991 .

[6]  Julia S. Mullen,et al.  Filter-based stabilization of spectral element methods , 2001 .

[7]  Franck Nicoud,et al.  Large-Eddy Simulation of Turbulence in Cardiovascular Flows , 2018 .

[8]  Hrvoje Jasak,et al.  A tensorial approach to computational continuum mechanics using object-oriented techniques , 1998 .

[9]  Julia S. Mullen,et al.  Filtering techniques for complex geometry fluid flows , 1999 .

[10]  Christophe Duwig,et al.  On the implementation of low-dissipative Runge–Kutta projection methods for time dependent flows using OpenFOAM® , 2014 .

[12]  Hrvoje Jasak,et al.  Error analysis and estimation for the finite volume method with applications to fluid flows , 1996 .

[13]  Leo G. Rebholz,et al.  Improved accuracy in regularization models of incompressible flow via adaptive nonlinear filtering , 2012 .

[14]  Christophe Prud'Homme,et al.  High order finite element simulations for fluid dynamics validated by experimental data from the fda benchmark nozzle model , 2017, 1701.02179.

[15]  Leo G. Rebholz,et al.  Numerical study of a regularization model for incompressible flow with deconvolution-based adaptive nonlinear filtering , 2013 .

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

[17]  P. Moin,et al.  Eddies, streams, and convergence zones in turbulent flows , 1988 .

[18]  Todd Allen Simons,et al.  Large eddy simulation of turbulent flows using finite volume methods with structured, unstructured, and zonal embedded grids , 1998 .

[19]  F. Nicoud,et al.  Using singular values to build a subgrid-scale model for large eddy simulations , 2011 .

[20]  J. Smagorinsky,et al.  GENERAL CIRCULATION EXPERIMENTS WITH THE PRIMITIVE EQUATIONS , 1963 .

[21]  C. Fureby,et al.  Numerical Simulation of an Oscillating Cylinder Using Large Eddy Simulation and Implicit Large Eddy Simulation , 2012 .

[22]  A Veneziani,et al.  Validation of an open source framework for the simulation of blood flow in rigid and deformable vessels , 2013, International journal for numerical methods in biomedical engineering.

[23]  G. Rozza,et al.  POD-Galerkin reduced order methods for CFD using Finite Volume Discretisation: vortex shedding around a circular cylinder , 2017, 1701.03424.

[24]  Kameswararao Anupindi,et al.  Large Eddy Simulation of FDA’s Idealized Medical Device , 2013, Cardiovascular engineering and technology.

[25]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

[26]  A. Quaini,et al.  On the Sensitivity to Model Parameters in a Filter Stabilization Technique for Advection Dominated Advection-Diffusion-Reaction Problems , 2018, Lecture Notes in Computational Science and Engineering.

[27]  Martin Kronbichler,et al.  Modern discontinuous Galerkin methods for the simulation of transitional and turbulent flows in biomedical engineering: A comprehensive LES study of the FDA benchmark nozzle model. , 2019, International journal for numerical methods in biomedical engineering.

[28]  Po-Lin Hsu,et al.  On the representation of effective stress for computing hemolysis , 2019, Biomechanics and modeling in mechanobiology.

[29]  A. Gosman,et al.  Solution of the implicitly discretised reacting flow equations by operator-splitting , 1986 .

[30]  Zhu Wang,et al.  Approximate Deconvolution Reduced Order Modeling , 2015, 1510.02726.

[31]  高等学校計算数学学報編輯委員会編.  高等学校計算数学学報 = Numerical mathematics , 1979 .

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

[33]  Leo G. Rebholz,et al.  Modular Nonlinear Filter Stabilization of Methods for Higher Reynolds Numbers Flow , 2012 .

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

[35]  D. Spalding,et al.  A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows , 1972 .

[36]  Traian Iliescu,et al.  A Bounded Artificial Viscosity Large Eddy Simulation Model , 2008, SIAM J. Numer. Anal..

[37]  Steven Deutsch,et al.  Multilaboratory particle image velocimetry analysis of the FDA benchmark nozzle model to support validation of computational fluid dynamics simulations. , 2011, Journal of biomechanical engineering.

[38]  M. Darwish,et al.  The Finite Volume Method in Computational Fluid Dynamics: An Advanced Introduction with OpenFOAM® and Matlab , 2015 .

[39]  A. Dunca,et al.  On the Stolz-Adams Deconvolution Model for the Large-Eddy Simulation of Turbulent Flows , 2006, SIAM J. Math. Anal..

[40]  Jean-Luc Guermond,et al.  From Suitable Weak Solutions to Entropy Viscosity , 2011, J. Sci. Comput..

[41]  G. Rozza,et al.  Finite volume POD-Galerkin stabilised reduced order methods for the parametrised incompressible Navier–Stokes equations , 2017, Computers & Fluids.

[42]  Jean Leray,et al.  Sur le mouvement d'un liquide visqueux emplissant l'espace , 1934 .

[43]  A. Goulas,et al.  On the order of accuracy of the divergence theorem (Green-Gauss) method for calculating the gradient in finite volume methods , 2017 .

[44]  Gábor Janiga,et al.  Large eddy simulation of the FDA benchmark nozzle for a Reynolds number of 6500 , 2014, Comput. Biol. Medicine.

[45]  P. Lax,et al.  Systems of conservation laws , 1960 .

[46]  R. F. Warming,et al.  Upwind Second-Order Difference Schemes and Applications in Aerodynamic Flows , 1976 .

[47]  Annalisa Quaini,et al.  Deconvolution‐based nonlinear filtering for incompressible flows at moderately large Reynolds numbers , 2016 .

[48]  K. Carlson,et al.  Turbulent Flows , 2020, Finite Analytic Method in Flows and Heat Transfer.

[49]  M. Olshanskii,et al.  A connection between filter stabilization and eddy viscosity models , 2013, 1302.4487.

[50]  J. P. V. Doormaal,et al.  ENHANCEMENTS OF THE SIMPLE METHOD FOR PREDICTING INCOMPRESSIBLE FLUID FLOWS , 1984 .

[51]  Franck Nicoud,et al.  About the numerical robustness of biomedical benchmark cases: Interlaboratory FDA's idealized medical device , 2017, International journal for numerical methods in biomedical engineering.

[52]  M. Darwish,et al.  The Finite Volume Method , 2016 .

[53]  A. W. Vreman.  An eddy-viscosity subgrid-scale model for turbulent shear flow: Algebraic theory and applications , 2004 .

[54]  Aslak W. Bergersen,et al.  The FDA nozzle benchmark: “In theory there is no difference between theory and practice, but in practice there is” , 2018, International journal for numerical methods in biomedical engineering.

[55]  Traian Iliescu,et al.  An evolve‐then‐filter regularized reduced order model for convection‐dominated flows , 2015, 1506.07555.

[56]  D. K. Walters,et al.  Laminar, Turbulent, and Transitional Simulations in Benchmark Cases with Cardiovascular Device Features , 2013 .