High-order simulations of isothermal flows using the local anisotropic basis function method (LABFM)

Mesh-free methods have significant potential for simulations of flows in complex geometries, with the difficulties of domain discretisation greatly reduced. However, many mesh-free methods are limited to low order accuracy. In order to compete with conventional mesh-based methods, high order accuracy is essential. The Local Anisotropic Basis Function Method (LABFM) is a meshfree method introduced in King et al., J. Comput. Phys. 415:109549 (2020), which enables the construction of highly accurate difference operators on disordered node discretisations. Here, we introduce a number of developments to LABFM, in the areas of basis function construction, stencil optimisation, stabilisation, variable resolution, and high order boundary conditions. With these developments, direct numerical simulations of the Navier Stokes equations are possible at extremely high order (up to 10 order in characteristic node spacing internally). We numerically solve the isothermal compressible Navier Stokes equations for a range of geometries: periodic and channel flows, flows past a cylinder, and porous media. Excellent agreement is seen with analytical solutions, published numerical results (using a spectral element method), and experiments. The potential of the method for direct numerical simulations in complex geometries is demonstrated with simulations of subsonic and transonic flows through an inhomogeneous porous media at pore Reynolds numbers up to Rep = 968.

[1]  Abouzied M. A. Nasar,et al.  High order difference schemes using the local anisotropic basis function method , 2019, J. Comput. Phys..

[2]  Mostafa Abbaszadeh,et al.  Simulation flows with multiple phases and components via the radial basis functions-finite difference (RBF-FD) procedure: Shan-Chen model , 2020 .

[3]  Thierry Poinsot,et al.  Comparison of outflow boundary conditions for subsonic aeroacoustic simulations , 2012 .

[4]  P. Alam.  ‘A’ , 2021, Composites Engineering: An A–Z Guide.

[5]  Stefano Sibilla,et al.  An algorithm to improve consistency in Smoothed Particle Hydrodynamics , 2015 .

[6]  Sourabh V. Apte,et al.  Modeling Turbulent Flows in Porous Media , 2020, Annual Review of Fluid Mechanics.

[7]  Harsha K. Chelliah,et al.  DNS Investigation of Cavity Stabilized Premixed Turbulent Ethylene-Air Flame , 2018 .

[8]  Qianlong Liu,et al.  Nonreflecting boundary conditions based on nonlinear multidimensional characteristics , 2007 .

[9]  D. Rempfer.  On Boundary Conditions for Incompressible Navier-Stokes Problems , 2006 .

[10]  James Demmel,et al.  LAPACK Users' Guide, Third Edition , 1999, Software, Environments and Tools.

[11]  Mehdi Dehghan,et al.  An upwind local radial basis functions-differential quadrature (RBF-DQ) method with proper orthogonal decomposition (POD) approach for solving compressible Euler equation , 2017, Engineering Analysis with Boundary Elements.

[12]  J. Morris,et al.  Modeling Low Reynolds Number Incompressible Flows Using SPH , 1997 .

[13]  Alessandro Reali,et al.  Novel finite particle formulations based on projection methodologies , 2011 .

[14]  Bengt Fornberg,et al.  Stabilization of RBF-generated finite difference methods for convective PDEs , 2011, J. Comput. Phys..

[15]  Bengt Fornberg,et al.  On the role of polynomials in RBF-FD approximations: II. Numerical solution of elliptic PDEs , 2017, J. Comput. Phys..

[16]  Axel Brandenburg,et al.  Computational aspects of astrophysical MHD and turbulence , 2001, Advances in Nonlinear Dynamos.

[17]  José M. Domínguez,et al.  A versatile algorithm for the treatment of open boundary conditions in Smoothed particle hydrodynamics GPU models , 2018, Computer Methods in Applied Mechanics and Engineering.

[18]  Kunihiko Taira,et al.  Two-dimensional compressible viscous flow around a circular cylinder , 2015, Journal of Fluid Mechanics.

[19]  N. D. Waters,et al.  Unsteady flow of an elastico-viscous liquid , 1970 .

[20]  B. Weigand,et al.  Direct numerical simulation of convective heat transfer in porous media , 2019, International Journal of Heat and Mass Transfer.

[21]  Wing Kam Liu,et al.  Meshfree and particle methods and their applications , 2002 .

[22]  V. Stovin,et al.  Estimating drag coefficient for arrays of rigid cylinders representing emergent vegetation , 2018, Journal of Hydraulic Research.

[23]  Peter Hansbo,et al.  CutFEM: Discretizing geometry and partial differential equations , 2015 .

[24]  Furen Ming,et al.  A novel non-reflecting boundary condition for fluid dynamics solved by smoothed particle hydrodynamics , 2018, Journal of Fluid Mechanics.

[25]  W. Dehnen,et al.  Improving convergence in smoothed particle hydrodynamics simulations without pairing instability , 2012, 1204.2471.

[26]  C. Ross Ethier,et al.  A high-order discontinuous Galerkin method for the unsteady incompressible Navier-Stokes equations , 2007, J. Comput. Phys..

[27]  H. Herwig,et al.  Numerical investigation of the possibility of macroscopic turbulence in porous media: a direct numerical simulation study , 2015, Journal of Fluid Mechanics.

[28]  Xiaozhe Hu,et al.  A compatible high-order meshless method for the Stokes equations with applications to suspension flows , 2016, J. Comput. Phys..

[29]  Wing Kam Liu,et al.  Reproducing kernel particle methods for structural dynamics , 1995 .

[30]  Xiaozhe Hu,et al.  Compact moving least squares: An optimization framework for generating high-order compact meshless discretizations , 2016, J. Comput. Phys..

[31]  K. Suga,et al.  Direct numerical simulation of turbulence over anisotropic porous media , 2017, Journal of Fluid Mechanics.

[32]  T. Poinsot.  Boundary conditions for direct simulations of compressible viscous flows , 1992 .

[33]  Christopher A. Kennedy,et al.  Improved boundary conditions for viscous, reacting, compressible flows , 2003 .

[34]  Wing Kam Liu,et al.  Reproducing kernel particle methods , 1995 .

[35]  Romesh C. Batra,et al.  Modified smoothed particle hydrodynamics method and its application to transient problems , 2004 .

[36]  Grady B. Wright,et al.  A Robust Hyperviscosity Formulation for Stable RBF-FD Discretizations of Advection-Diffusion-Reaction Equations on Manifolds , 2020, SIAM J. Sci. Comput..

[37]  Onset of Vortex Shedding in a Periodic Array of Circular Cylinders , 2006 .

[38]  W.,et al.  Time Dependent Boundary Conditions for Hyperbolic Systems , 2003 .

[39]  Sahil Garg,et al.  Meshfree Methods: A Comprehensive Review of Applications , 2017 .

[40]  G. Karniadakis,et al.  Spectral/hp Element Methods for Computational Fluid Dynamics , 2005 .

[41]  Boyce E. Griffith,et al.  Immersed Methods for Fluid-Structure Interaction. , 2020, Annual review of fluid mechanics.

[42]  Benedict D. Rogers,et al.  High-order velocity and pressure wall boundary conditions in Eulerian incompressible SPH , 2020, J. Comput. Phys..

[43]  A. Kuznetsov,et al.  Turbulence modeling for flows in wall bounded porous media: An analysis based on direct numerical simulations , 2017 .

[44]  P. Hopkins.  A new class of accurate, mesh-free hydrodynamic simulation methods , 2014, 1409.7395.

[45]  Gianluca Iaccarino,et al.  IMMERSED BOUNDARY METHODS , 2005 .

[46]  M. Carpenter,et al.  Several new numerical methods for compressible shear-layer simulations , 1994 .

[47]  Emmanuel Motheau,et al.  Navier–Stokes Characteristic Boundary Conditions Using Ghost Cells , 2017 .

[48]  Jing Tang Xing,et al.  Shape adaptive RBF-FD implicit scheme for incompressible viscous Navier–Strokes equations , 2014 .

[49]  Sudarshan Tiwari,et al.  On meshfree GFDM solvers for the incompressible Navier–Stokes equations , 2017, 1701.03427.

[50]  Hong G. Im,et al.  Characteristic boundary conditions for simulations of compressible reacting flows with multi-dimensional, viscous and reaction effects , 2007 .

[51]  George Em Karniadakis,et al.  Spectral element simulations of laminar and turbulent flows in complex geometries , 1989 .

[52]  Bengt Fornberg,et al.  Solving PDEs with radial basis functions * , 2015, Acta Numerica.

[53]  Xiaozhe Hu,et al.  A spatially adaptive high-order meshless method for fluid–structure interactions , 2019, Computer Methods in Applied Mechanics and Engineering.

[54]  C. Shu,et al.  An upwind local RBF-DQ method for simulation of inviscid compressible flows , 2005 .

[55]  Bengt Fornberg,et al.  Fast generation of 2-D node distributions for mesh-free PDE discretizations , 2015, Comput. Math. Appl..

[56]  P. Sagaut,et al.  Simulation of a viscous compressible flow past a circular cylinder with high-order discontinuous Galerkin methods , 2002 .

[57]  J. Bonet,et al.  Variational and momentum preservation aspects of Smooth Particle Hydrodynamic formulations , 1999 .

[58]  Qian Wang,et al.  AUGEM: Automatically generate high performance Dense Linear Algebra kernels on x86 CPUs , 2013, 2013 SC - International Conference for High Performance Computing, Networking, Storage and Analysis (SC).

[59]  T. Fries.  Higher-order conformal decomposition FEM (CDFEM) , 2018 .

[60]  S. Apte,et al.  Relative performance of body fitted and fictitious domain simulations of flow through fixed packed beds of spheres , 2013 .

[61]  Steven A. Orszag,et al.  On the Elimination of Aliasing in Finite-Difference Schemes by Filtering High-Wavenumber Components , 1971 .

[62]  Varun Shankar,et al.  Hyperviscosity-Based Stabilization for Radial Basis Function-Finite Difference (RBF-FD) Discretizations of Advection-Diffusion Equations , 2018, J. Comput. Phys..

[63]  Nikolaus A. Adams,et al.  A generalized wall boundary condition for smoothed particle hydrodynamics , 2012, J. Comput. Phys..

[64]  B. Fornberg,et al.  Radial basis function interpolation: numerical and analytical developments , 2003 .

[65]  Robert D. Guy,et al.  Immersed Boundary Smooth Extension (IBSE): A high-order method for solving incompressible flows in arbitrary smooth domains , 2016, J. Comput. Phys..

[66]  D. Lohse,et al.  From Rayleigh–Bénard convection to porous-media convection: how porosity affects heat transfer and flow structure , 2020, Journal of Fluid Mechanics.

[67]  Luis Gavete,et al.  Solving second order non-linear elliptic partial differential equations using generalized finite difference method , 2017, J. Comput. Appl. Math..

[68]  Yukie Tanino,et al.  Laboratory Investigation of Mean Drag in a Random Array of Rigid, Emergent Cylinders , 2008 .

[69]  Richard Shock,et al.  Numerical study of flow past an impulsively started cylinder by the lattice-Boltzmann method , 2004, Journal of Fluid Mechanics.

[70]  Yi Wang,et al.  Characteristic boundary conditions for direct simulations of turbulent counterflow flames , 2005 .

[71]  Madeleine Coutanceau,et al.  The early stage of development of the wake behind an impulsively started cylinder for 40 < Re < 104 , 1980, Journal of Fluid Mechanics.

[72]  P. A. Smith,et al.  Impulsively started flow around a circular cylinder by the vortex method , 1988, Journal of Fluid Mechanics.

[73]  R. Lewis,et al.  Low-storage, Explicit Runge-Kutta Schemes for the Compressible Navier-Stokes Equations , 2000 .

[74]  Pavel B. Bochev,et al.  A High-Order Staggered Meshless Method for Elliptic Problems , 2017, SIAM J. Sci. Comput..

[75]  L. Gavete,et al.  Solving parabolic and hyperbolic equations by the generalized finite difference method , 2007 .

[76]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[77]  Tim Colonius,et al.  The immersed boundary method: A projection approach , 2007, J. Comput. Phys..

[78]  Steven J. Lind,et al.  High-order Eulerian incompressible smoothed particle hydrodynamics with transition to Lagrangian free-surface motion , 2016, J. Comput. Phys..

[79]  Giancarlo Sangalli,et al.  Particle Methods for a 1D Elastic Model Problem: Error Analysis and Development of a Second-Order Accurate Formulation , 2010 .