An implicit boundary approach for viscous compressible high Reynolds flows using a hybrid remeshed particle hydrodynamics method

Abstract We present an implicit boundary particle method with background mesh adaptation. We use a Brinkman penalisation to represent the boundary of the domain and a remeshed particle method to simulate viscous flow with high Reynolds numbers. A penalty term is added to the Navier-Stokes equations to impose the boundary conditions. The boundary conditions are enforced to a specific precision with no need to modify the numerical method or change the grid, achieving an implicit approach for flow around complex boundaries/geometries simulation. The main idea of the Brinkman penalisation method is to model the solid obstacle as a porous medium. The governing equations for the compressible fluid and penalised Navier-Stokes for the porous medium are solved simultaneously, without need for interface conditions. The accuracy of the method is tested for a number of benchmark problems starting with simple cases such as a periodic laminar flow in a channel (Poiseuille flow), to more complex problems such as the lid-driven cavity with high Reynolds numbers, the reflection and propagation of a compressible shock wave, and finally swirling flow in a model two-stroke diesel engine. The remeshed particle-mesh method coupled with Brinkman penalisation provides a good quality simulation and the results are in agreement with analytical or reference solutions.

[1]  M. Lai,et al.  An Immersed Boundary Method with Formal Second-Order Accuracy and Reduced Numerical Viscosity , 2000 .

[2]  G. Cottet,et al.  A level-set formulation of immersed boundary methods for fluid–structure interaction problems , 2004 .

[3]  Anthony Peter Whitworth,et al.  A new prescription for viscosity in smoothed particle hydrodynamics. , 1996 .

[4]  Yiannis Ventikos,et al.  Dual pulsating or steady slot jet cooling of a constant heat flux surface , 2003 .

[5]  C. Peskin Flow patterns around heart valves: A numerical method , 1972 .

[6]  Philippe Angot,et al.  Analysis of singular perturbations on the Brinkman problem for fictitious domain models of viscous flows , 1999 .

[7]  L. Sirovich,et al.  Modeling a no-slip flow boundary with an external force field , 1993 .

[8]  J. Monaghan Simulating Free Surface Flows with SPH , 1994 .

[9]  Stéphane P. A. Bordas,et al.  Corotational Cut Finite Element Method for real-time surgical simulation: application to needle insertion simulation , 2017, Computer Methods in Applied Mechanics and Engineering.

[10]  C. L. Tien,et al.  Boundary and inertia effects on flow and heat transfer in porous media , 1981 .

[11]  Philippe Angot,et al.  A penalization method to take into account obstacles in incompressible viscous flows , 1999, Numerische Mathematik.

[12]  L. Lucy A numerical approach to the testing of the fission hypothesis. , 1977 .

[13]  J. Monaghan Particle methods for hydrodynamics , 1985 .

[14]  James L. Beck,et al.  Convection in a box of porous material saturated with fluid , 1972 .

[15]  M. Kaviany Principles of heat transfer in porous media , 1991 .

[16]  H. Brinkman A calculation of the viscous force exerted by a flowing fluid on a dense swarm of particles , 1949 .

[17]  Guido Morgenthal,et al.  An immersed interface method for the Vortex-In-Cell algorithm , 2007 .

[18]  P. Cleary,et al.  Conduction Modelling Using Smoothed Particle Hydrodynamics , 1999 .

[19]  Jens Honore Walther,et al.  PIV study of the effect of piston position on the in-cylinder swirling flow during the scavenging process in large two-stroke marine diesel engines , 2013 .

[20]  Yuri Bazilevs,et al.  An immersogeometric variational framework for fluid-structure interaction: application to bioprosthetic heart valves. , 2015, Computer methods in applied mechanics and engineering.

[21]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[22]  Petros Koumoutsakos,et al.  An immersed boundary-lattice-Boltzmann method for the simulation of the flow past an impulsively started cylinder , 2008, J. Comput. Phys..

[23]  J. Morris Simulating surface tension with smoothed particle hydrodynamics , 2000 .

[24]  Oleg V. Vasilyev,et al.  An Adaptive Wavelet Collocation Method for Fluid-Structure Interaction at High Reynolds Numbers , 2005, SIAM J. Sci. Comput..

[25]  J. Ferziger,et al.  A ghost-cell immersed boundary method for flow in complex geometry , 2002 .

[26]  M. Berger,et al.  An Adaptive Version of the Immersed Boundary Method , 1999 .

[27]  G. Cottet,et al.  A LEVEL SET METHOD FOR FLUID-STRUCTURE INTERACTIONS WITH IMMERSED SURFACES , 2006 .

[28]  P. Moin,et al.  Numerical Simulation of Turbulent Flows , 1984 .

[29]  J. Trulsen,et al.  Regularized Smoothed Particle Hydrodynamics: A New Approach to Simulating Magnetohydrodynamic Shocks , 2001 .

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

[31]  Oleg V. Vasilyev,et al.  A Brinkman penalization method for compressible flows in complex geometries , 2007, J. Comput. Phys..

[32]  Martin Robinson,et al.  Direct numerical simulation of decaying two‐dimensional turbulence in a no‐slip square box using smoothed particle hydrodynamics , 2011, 1101.2240.

[33]  T. Lund,et al.  A Non-Body Conformal Grid Method for Simulation of Compressible Flows with Complex Immersed Boundaries , 2004 .

[34]  Damien Violeau,et al.  Numerical modelling of complex turbulent free‐surface flows with the SPH method: an overview , 2007 .

[35]  J. Monaghan,et al.  Smoothed particle hydrodynamics: Theory and application to non-spherical stars , 1977 .

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

[37]  R. Wooding Steady state free thermal convection of liquid in a saturated permeable medium , 1957, Journal of Fluid Mechanics.

[38]  J. Monaghan SPH compressible turbulence , 2002, astro-ph/0204118.

[39]  Markus H. Gross,et al.  Particle-based fluid simulation for interactive applications , 2003, SCA '03.

[40]  D. D. Zeeuw,et al.  An adaptively refined Cartesian mesh solver for the Euler equations , 1993 .

[41]  Petros Koumoutsakos,et al.  Remeshed smoothed particle hydrodynamics for the simulation of viscous and heat conducting flows , 2002 .

[42]  J. W. Purvis,et al.  Prediction of critical Mach number for store configurations , 1979 .

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

[44]  J. Trulsen,et al.  Regularized smoothed particle hydrodynamics with improved multi-resolution handling , 2005 .

[45]  John A. Evans,et al.  An Isogeometric design-through-analysis methodology based on adaptive hierarchical refinement of NURBS, immersed boundary methods, and T-spline CAD surfaces , 2012 .

[46]  Stéphane Bordas,et al.  Three‐dimensional remeshed smoothed particle hydrodynamics for the simulation of isotropic turbulence , 2018 .

[47]  S. Cummins,et al.  An SPH Projection Method , 1999 .

[48]  S. Shao,et al.  INCOMPRESSIBLE SPH METHOD FOR SIMULATING NEWTONIAN AND NON-NEWTONIAN FLOWS WITH A FREE SURFACE , 2003 .

[49]  Eugene Fiume,et al.  Depicting fire and other gaseous phenomena using diffusion processes , 1995, SIGGRAPH.

[50]  S. Miyama,et al.  Numerical Simulation of Viscous Flow by Smoothed Particle Hydrodynamics , 1994 .

[51]  S. Biringen,et al.  Numerical Simulation of a Cylinder in Uniform Flow , 1996 .

[52]  Petros Koumoutsakos,et al.  An immersed boundary method for smoothed particle hydrodynamics of self-propelled swimmers , 2008, J. Comput. Phys..

[53]  C. Peskin The immersed boundary method , 2002, Acta Numerica.

[54]  Yuri Bazilevs,et al.  A new formulation for air-blast fluid–structure interaction using an immersed approach: part II—coupling of IGA and meshfree discretizations , 2017 .

[55]  O. M. Phillips,et al.  Flow and Reactions in Permeable Rocks , 1991 .

[56]  P. Koumoutsakos,et al.  A Lagrangian particle level set method. , 2005 .

[57]  Elie Rivoalen,et al.  JOSEPHINE: A parallel SPH code for free-surface flows , 2012, Comput. Phys. Commun..

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

[59]  K. Meyer,et al.  Large eddy simulations of the influence of piston position on the swirling flow in a model two-stroke diesel engine , 2014 .

[60]  M. D. Salas,et al.  Euler calculations for multielement airfoils using Cartesian grids , 1986 .

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

[62]  Joe F. Thompson,et al.  Boundary-fitted coordinate systems for numerical solution of partial differential equations—A review , 1982 .

[63]  Benedict D. Rogers,et al.  Numerical Modeling of Water Waves with the SPH Method , 2006 .