Comparison of different iterative schemes for ISPH based on Rankine source solution

Smoothed Particle Hydrodynamics (SPH) method has a good adaptability for the simulation of free surface flow problems. There are two forms of SPH. One is weak compressible SPH and the other one is incompressible SPH (ISPH). Compared with the former one, ISPH method performs better in many cases. ISPH based on Rankine source solution can perform better than traditional ISPH, as it can use larger stepping length by avoiding the second order derivative in pressure Poisson equation. However, ISPH_R method needs to solve the sparse linear matrix for pressure Poisson equation, which is one of the most expensive parts during one time stepping calculation. Iterative methods are normally used for solving Poisson equation with large particle numbers. However, there are many iterative methods available and the question for using which one is still open. In this paper, three iterative methods, CGS, Bi-CGstab and GMRES are compared, which are suitable and typical for large unsymmetrical sparse matrix solutions. According to the numerical tests on different cases, still water test, dam breaking, violent tank sloshing, solitary wave slamming, the GMRES method is more efficient than CGS and Bi-CGstab for ISPH method.

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

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

[3]  S. Shao,et al.  Corrected Incompressible SPH method for accurate water-surface tracking in breaking waves , 2008 .

[4]  D. R. Fokkema,et al.  BICGSTAB( L ) FOR LINEAR EQUATIONS INVOLVING UNSYMMETRIC MATRICES WITH COMPLEX , 1993 .

[5]  A. Colagrossi,et al.  Numerical simulation of interfacial flows by smoothed particle hydrodynamics , 2003 .

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

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

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

[9]  Gui-Rong Liu,et al.  Restoring particle consistency in smoothed particle hydrodynamics , 2006 .

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

[11]  Masashi Kashiwagi,et al.  Numerical simulation of violent sloshing by a CIP-based method , 2006 .

[12]  Martin B. van Gijzen,et al.  IDR(s): A Family of Simple and Fast Algorithms for Solving Large Nonsymmetric Systems of Linear Equations , 2008, SIAM J. Sci. Comput..

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

[14]  Afzal Suleman,et al.  A robust weakly compressible SPH method and its comparison with an incompressible SPH , 2012 .

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

[16]  Seiji Fujino,et al.  GPBiCG(m, l): a hybrid of BiCGSTAB and GPBiCG methods with efficiency and robustness , 2002 .

[17]  H. Schwaiger.  An implicit corrected SPH formulation for thermal diffusion with linear free surface boundary conditions , 2008 .

[18]  Judith A. Vogel,et al.  Flexible BiCG and flexible Bi-CGSTAB for nonsymmetric linear systems , 2007, Appl. Math. Comput..

[19]  Wen-yang Duan,et al.  Incompressible SPH Based on Rankine Source Solution for Water Wave Impact Simulation , 2015 .

[20]  Juntao Zhou,et al.  MLPG_R Method for Numerical Simulation of 2D Breaking Waves , 2009 .

[21]  Murray Rudman,et al.  Comparative study on the accuracy and stability of SPH schemes in simulating energetic free-surface flows , 2012 .

[22]  Songdong Shao,et al.  Incompressible SPH simulation of water entry of a free‐falling object , 2009 .

[23]  R. Mittal,et al.  An efficient method for constructing an ILU preconditioner for solving large sparse nonsymmetric linear systems by the GMRES method , 2003 .

[24]  Pep Español,et al.  Incompressible smoothed particle hydrodynamics , 2007, J. Comput. Phys..

[25]  G. Golub,et al.  Gmres: a Generalized Minimum Residual Algorithm for Solving , 2022 .

[26]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[27]  C. K. Thornhill,et al.  Part IV. An experimental study of the collapse of liquid columns on a rigid horizontal plane , 1952, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[28]  Hitoshi Gotoh,et al.  Enhancement of stability and accuracy of the moving particle semi-implicit method , 2011, J. Comput. Phys..

[29]  Xin Liu,et al.  An improved incompressible SPH model for simulation of wave–structure interaction , 2013 .

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

[31]  Andreas G. Boudouvis,et al.  Bifurcation detection with the (un)preconditioned GMRES(m) , 2004 .

[32]  P. Sonneveld.  CGS, A Fast Lanczos-Type Solver for Nonsymmetric Linear systems , 1989 .

[33]  D. Graham,et al.  Simulation of wave overtopping by an incompressible SPH model , 2006 .

[34]  Bertrand Alessandrini,et al.  An improved SPH method: Towards higher order convergence , 2007, J. Comput. Phys..

[35]  Q. W. Ma,et al.  Incompressible SPH method based on Rankine source solution for violent water wave simulation , 2014, J. Comput. Phys..

[36]  Hitoshi Gotoh,et al.  GPU-acceleration for Moving Particle Semi-Implicit method , 2011 .

[37]  S. Koshizuka,et al.  Moving-Particle Semi-Implicit Method for Fragmentation of Incompressible Fluid , 1996 .

[38]  Yoshiaki Oka,et al.  Numerical analysis of fragmentation mechanisms in vapor explosions , 1999 .

[39]  K. Morita,et al.  An improved MPS method for numerical simulations of convective heat transfer problems , 2006 .

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

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