Unified semi-analytical wall boundary conditions applied to 2-D incompressible SPH

This work aims at improving the 2-D incompressible SPH model (ISPH) by adapting it to the unified semi-analytical wall boundary conditions proposed by Ferrand et al. [10]. The ISPH algorithm considered is as proposed by Lind et al. [25], based on the projection method with a divergence-free velocity field and using a stabilising procedure based on particle shifting. However, we consider an extension of this model to Reynolds-Averaged Navier-Stokes equations based on the k-@e turbulent closure model, as done in [10]. The discrete SPH operators are modified by the new description of the wall boundary conditions. In particular, a boundary term appears in the Laplacian operator, which makes it possible to accurately impose a von Neumann pressure wall boundary condition that corresponds to impermeability. The shifting and free-surface detection algorithms have also been adapted to the new boundary conditions. Moreover, a new way to compute the wall renormalisation factor in the frame of the unified semi-analytical boundary conditions is proposed in order to decrease the computational time. We present several verifications to the present approach, including a lid-driven cavity, a water column collapsing on a wedge and a periodic schematic fish-pass. Our results are compared to Finite Volumes methods, using Volume of Fluids in the case of free-surface flows. We briefly investigate the convergence of the method and prove its ability to model complex free-surface and turbulent flows. The results are generally improved when compared to a weakly compressible SPH model with the same boundary conditions, especially in terms of pressure prediction.

[1]  Roland W. Lewis,et al.  A variational formulation based contact algorithm for rigid boundaries in two-dimensional SPH applications , 2004 .

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

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

[4]  Damien Violeau,et al.  On the maximum time step in weakly compressible SPH , 2014, J. Comput. Phys..

[5]  Daniel J. Price.  Smoothed particle hydrodynamics and magnetohydrodynamics , 2010, J. Comput. Phys..

[6]  Michael Dumbser,et al.  A new 3D parallel SPH scheme for free surface flows , 2009 .

[7]  Nikolaus A. Adams,et al.  Assessment of Implicit Large-Eddy Simulation with a conservative immersed interface method for turbulent cylinder flow , 2010 .

[8]  Benedict D. Rogers,et al.  Study of differential operators in the context of the semi-analytical wall boundary conditions , 2012 .

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

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

[11]  Dominique Laurence,et al.  A LINEARISED TURBULENT PRODUCTION IN THE k-ε MODEL FOR ENGINEERING APPLICATIONS , 2002 .

[12]  Christophe Kassiotis,et al.  Unified semi-analytical wall boundary conditions in SPH: analytical extension to 3-D , 2014, Numerical Algorithms.

[13]  Nikolaus A. Adams,et al.  An incompressible multi-phase SPH method , 2007, J. Comput. Phys..

[14]  Bertrand Alessandrini,et al.  Specific pre/post treatments for 3-D SPH applications through massive HPC simulations , 2009 .

[15]  S. Shao,et al.  Incompressible SPH simulation of wave breaking and overtopping with turbulence modelling , 2006 .

[16]  Dominique Laurence,et al.  Incompressible separated flows simulations with the smoothed particle hydrodynamics gridless method , 2005 .

[17]  Rui Xu,et al.  Comparisons of weakly compressible and truly incompressible algorithms for the SPH mesh free particle method , 2008, J. Comput. Phys..

[18]  Javier Bonet,et al.  Dynamic refinement and boundary contact forces in SPH with applications in fluid flow problems , 2007 .

[19]  Nikolaus A. Adams,et al.  SPH simulations of flow around a periodic array of cylinders confined in a channel , 2011 .

[20]  B. Launder,et al.  Mathematical Models of turbulence , 1972 .

[21]  Yih-Ferng Peng,et al.  Transition in a 2-D lid-driven cavity flow , 2003 .

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

[23]  Damien Violeau,et al.  ONE AND TWO-EQUATIONS TURBULENT CLOSURES FOR SMOOTHED PARTICLE HYDRODYNAMICS , 2004 .

[24]  J. Monaghan.  Smoothed particle hydrodynamics , 2005 .

[25]  Joe J. Monaghan,et al.  SPH particle boundary forces for arbitrary boundaries , 2009, Comput. Phys. Commun..

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

[27]  Robert C. Armstrong,et al.  Viscoelastic flow of polymer solutions around a periodic, linear array of cylinders: comparisons of predictions for microstructure and flow fields , 1998 .

[28]  Afzal Suleman,et al.  SPH with the multiple boundary tangent method , 2009 .

[29]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[30]  Bertrand Alessandrini,et al.  Normal flux method at the boundary for SPH , 2009 .

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

[32]  Salvatore Marrone,et al.  Fast free-surface detection and level-set function definition in SPH solvers , 2010, J. Comput. Phys..

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

[34]  Yuichi Matsuo,et al.  Direct Numerical Simulation of a Fully Developed Turbulent Channel Flow With Respect to the Reynolds Number Dependence , 2001 .

[35]  Rui Xu,et al.  Accuracy and stability in incompressible SPH (ISPH) based on the projection method and a new approach , 2009, J. Comput. Phys..

[36]  Alexandre J. Chorin,et al.  On the Convergence of Discrete Approximations to the Navier-Stokes Equations , 1969 .

[37]  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..

[38]  Nikolaus A. Adams,et al.  A conservative immersed interface method for Large-Eddy Simulation of incompressible flows , 2010, J. Comput. Phys..

[39]  F. Archambeau,et al.  Code Saturne: A Finite Volume Code for the computation of turbulent incompressible flows - Industrial Applications , 2004 .

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

[41]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..