High-order velocity and pressure wall boundary conditions in Eulerian incompressible SPH

Abstract High-order velocity and pressure boundary conditions are presented in Eulerian incompressible smoothed particle hydrodynamics (ISPH). While the high-order convergence of Eulerian ISPH has been demonstrated by the authors for periodic internal flows using Gaussian kernels this was limited by first to second-order accuracy for cases with solid boundaries. Since the SPH interpolation method is numerically robust there is potential for obtaining high-order accuracy in topologically complex domains with robust high-order accurate boundary conditions. In this paper high-order finite-difference extrapolation methods at solid boundaries are developed in Eulerian ISPH to allow for enforcement of the Dirichlet boundary condition for velocity and the Neumann boundary condition for pressure with high-order accuracy. Convergence up to fourth-order is demonstrated for 2-D Taylor-Couette flow and 3-D simulations of Taylor-Couette cellular flow structures are used to demonstrate accuracy and robustness. The order of accuracy may be extended to even higher-order using the analysis presented. Compact fourth-order Wendland-type kernels have also been derived to reduce the particle support region thereby lowering computational effort without loss of high-order convergence. The proposed formulation is therefore entirely high order.

[1]  José M. Domínguez,et al.  Local uniform stencil (LUST) boundary condition for arbitrary 3-D boundaries in parallel smoothed particle hydrodynamics (SPH) models , 2019, Computers & Fluids.

[2]  Mauro De Marchis,et al.  A coupled Finite Volume–Smoothed Particle Hydrodynamics method for incompressible flows , 2016 .

[3]  Jose L. Cercos-Pita,et al.  A Boundary Integral SPH Formulation --- Consistency and Applications to ISPH and WCSPH --- , 2012 .

[4]  M. Lastiwka,et al.  Truncation error in mesh‐free particle methods , 2006 .

[5]  Bertrand Alessandrini,et al.  Violent Fluid-Structure Interaction simulations using a coupled SPH/FEM method , 2010 .

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

[7]  Dimos Poulikakos,et al.  High order interpolation and differentiation using B-splines , 2004 .

[8]  I. Babuska.  The finite element method with Lagrangian multipliers , 1973 .

[9]  S. J. Lind,et al.  Incompressible smoothed particle hydrodynamics for free-surface flows: A generalised diffusion-based algorithm for stability and validations for impulsive flows and propagating waves , 2012, J. Comput. Phys..

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

[11]  James J. Feng,et al.  Pressure boundary conditions for computing incompressible flows with SPH , 2011, J. Comput. Phys..

[12]  A. Chorin.  The Numerical Solution of the Navier-Stokes Equations for an Incompressible Fluid , 2015 .

[13]  M. Gómez-Gesteira,et al.  Boundary conditions generated by dynamic particles in SPH methods , 2007 .

[14]  Alistair Revell,et al.  Flexible slender body fluid interaction: Vector-based discrete element method with Eulerian smoothed particle hydrodynamics , 2019, Computers & Fluids.

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

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

[17]  F. Brezzi.  On the existence, uniqueness and approximation of saddle-point problems arising from lagrangian multipliers , 1974 .

[18]  Jie Shen,et al.  An overview of projection methods for incompressible flows , 2006 .

[19]  Alistair Revell,et al.  Eulerian weakly compressible smoothed particle hydrodynamics (SPH) with the immersed boundary method for thin slender bodies , 2019, Journal of Fluids and Structures.

[20]  A. Colagrossi,et al.  δ-SPH model for simulating violent impact flows , 2011 .

[21]  T. Mullin.  Mutations of steady cellular flows in the Taylor experiment , 1982, Journal of Fluid Mechanics.

[22]  L. Brookshaw,et al.  A Method of Calculating Radiative Heat Diffusion in Particle Simulations , 1985, Publications of the Astronomical Society of Australia.

[23]  Jie Shen,et al.  On the error estimates for the rotational pressure-correction projection methods , 2003, Math. Comput..

[24]  Lloyd N. Trefethen,et al.  The Exponentially Convergent Trapezoidal Rule , 2014, SIAM Rev..

[25]  Jean-Luc Guermond,et al.  On the approximation of the unsteady Navier–Stokes equations by finite element projection methods , 1998, Numerische Mathematik.

[26]  Benedict D. Rogers,et al.  An Eulerian-Lagrangian incompressible SPH formulation (ELI-SPH) , 2018 .

[27]  Frans N. van de Vosse,et al.  An approximate projec-tion scheme for incompressible ow using spectral elements , 1996 .

[28]  R. Fatehi,et al.  Error estimation in smoothed particle hydrodynamics and a new scheme for second derivatives , 2011, Comput. Math. Appl..

[29]  D. Violeau,et al.  Smoothed particle hydrodynamics (SPH) for free-surface flows: past, present and future , 2016 .

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

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

[32]  Hitoshi Gotoh,et al.  Current achievements and future perspectives for projection-based particle methods with applications in ocean engineering , 2016 .

[33]  Dominique Laurence,et al.  Unified semi‐analytical wall boundary conditions for inviscid, laminar or turbulent flows in the meshless SPH method , 2013 .

[34]  Xiaohu Guo,et al.  New massively parallel scheme for Incompressible Smoothed Particle Hydrodynamics (ISPH) for highly nonlinear and distorted flow , 2018, Comput. Phys. Commun..

[35]  James J. Feng,et al.  A particle-based model for the transport of erythrocytes in capillaries , 2009 .

[36]  J. Monaghan,et al.  Extrapolating B splines for interpolation , 1985 .

[37]  Holger Wendland,et al.  Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree , 1995, Adv. Comput. Math..

[38]  Hitoshi Gotoh,et al.  Comparative study on accuracy and conservation properties of two particle regularization schemes and proposal of an optimized particle shifting scheme in ISPH context , 2017, J. Comput. Phys..

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