A remedy for numerical oscillations in weakly compressible smoothed particle hydrodynamics

Weakly Compressible Smoothed Particle Hydrodynamics (WCSPH) can lead to non‐physical oscillations in the pressure and density fields when simulating incompressible flow problems. This in turn may result in tensile instability and sometimes divergence. In this paper, it is shown that this difficulty originates from the specific form of spatial discretization used for the pressure term when solving the mass conservation equation. After describing the pressure–velocity decoupling problem associated with the so‐called colocated grid methods, a modified approach is presented that overcomes this problem using a different discretization scheme for the second derivative of pressure. The modified scheme is employed for solving a number of benchmark problems including both single‐phase and two‐phase test cases. Copyright © 2010 John Wiley & Sons, Ltd.     

Towards a transparent boundary condition for compressible Navier–Stokes equations

A new artificial boundary condition for two‐dimensional subsonic flows governed by the compressible Navier–Stokes equations is derived. It is based on the hyperbolic part of the equations, according to the way of propagation of the characteristic waves. A reference flow, as well as a convection velocity, is used to properly discretize the terms corresponding to the entering waves. Numerical tests on various classical model problems, whose solutions are known, and comparisons with other boundary conditions (BCs), show the efficiency of the BC. Direct numerical simulations of more complex flows over a dihedral plate are simulated, without creation of acoustic waves going back in the flow. Copyright © 2001 John Wiley & Sons, Ltd.     

A complete boundary integral formulation for compressible Navier–Stokes equations

A complete boundary integral formulation for compressible Navier–Stokes equations with time discretization by operator splitting is developed using the fundamental solutions of the Helmholtz operator equation with different order. The numerical results for wall pressure and wall skin friction of two‐dimensional compressible laminar viscous flow around airfoils are in good agreement with field numerical methods. Copyright © 2004 John Wiley & Sons, Ltd.     

A consistent incompressible SPH method for internal flows with fixed and moving boundaries

An improved incompressible smoothed particle hydrodynamics (ISPH) method is presented, which employs first‐order consistent discretization schemes both for the first‐order and second‐order spatial derivatives. A recently introduced wall boundary condition is implemented in the context of ISPH method, which does not rely on using dummy particles and, as a result, can be applied more efficiently and with less computational complexity. To assess the accuracy and computational efficiency of this improved ISPH method, a number of two‐dimensional incompressible laminar internal flow benchmark problems are solved and the results are compared with available analytical solutions and numerical data. It is shown that using smaller smoothing lengths, the proposed method can provide desirable accuracies with relatively less computational cost for two‐dimensional problems. Copyright © 2015 John Wiley & Sons, Ltd.     

A high‐order scheme for the incompressible Navier–Stokes equations with open boundary condition

We present a numerical scheme to solve the incompressible Navier–Stokes equations with open boundary condition. After replacing the incompressibility constraint by the pressure Poisson equation, the key is how to give an appropriate boundary condition for the pressure Poisson equation. We propose a new boundary condition for the pressure on the open boundary. Some numerical experiments are presented to verify the accuracy and stability of scheme. Copyright © 2013 John Wiley & Sons, Ltd.     

光滑粒子法中的一种新的核函数

分析了传统的核函数产生压缩失稳现象的原因,提出了消除这种压缩失稳现象的一种新的核函数。采用改进的光滑粒子法,对几种常用的核函数进行了一维应变波的对比计算。结果表明:所提出新的核函数在应力波计算中不但保证了计算精度,还能有效地消除压缩失稳。     

Accurate particle splitting for smoothed particle hydrodynamics in shallow water with shock capturing

The solution for the shallow water equations using smoothed particle hydrodynamics is attractive, being a mesh‐free, automatically adaptive method without special treatment for wet–dry interfaces. However, the relatively new method is limited by the variable kernel size or smoothing length being inversely proportional to water depth causing poor resolution at small depths. Boundary conditions at solid walls have also not been well resolved. To solve the resolution problem in small depths, a particle splitting procedure was developed (conveniently into seven particles), which conserves mass and momentum by varying the smoothing length, velocity and acceleration of each refined particle. This improves predictions in the shallowest depths where the error associated with splitting is reduced by one order of magnitude in comparison to other published works. To provide good shock capturing behaviour, particle interactions are treated as a Riemann problem with Monotone Upstream‐centred Scheme for Conservation Laws (MUSCL) reconstruction providing stability. For solid boundaries, the recent modified virtual boundary particle method was developed further to enable the zeroth moment to be accurately conserved where the smoothing length of particles is changing rapidly during particle splitting. The resulting method is applied to the one‐dimensional and the two‐dimensional axisymmetric wet‐bed dam break problems showing close agreement with analytical solutions, demonstrating the need for particle splitting. To demonstrate wetting and drying in a more complex case, the scheme is applied to oscillating water in a two‐dimensional parabolic basin and produces good agreement with the analytical solution. The method is finally applied to the European Concerted Action on DAm break Modelling dam‐break test case representative of realistic conditions and good predictions are made of experimental measurements with a 40% reduction in the computational time when particle splitting is employed. The overall method has thus become quite sophisticated but its generality and versatility will be attractive for various shallow water problems. Copyright © 2011 John Wiley & Sons, Ltd.     

A critical investigation of smoothed particle hydrodynamics applied to problems with free‐surfaces

In this paper, an in‐depth study of SPH method, in its original weakly compressible version, is achieved on dedicated 2D and 3D free‐surface flow test cases. These rather critical prototype problems shall constitute suitable test cases to get through when building a free‐surface SPH model. The present work aims at investigating various numerical aspects of this method, often little mentioned in literature. In particular, a great care is paid to the dynamic part of the solution, which is critical to the local hydrodynamic load prediction. The role of numerical errors in the development of acoustic frequencies in the pressure signals is discussed, as well as the influence of the choice of the sound velocity. On the shown test problems, it is also evidenced that some numerical tools are crucial to ensure the robustness and accuracy of the standard SPH method. The convergence of our model is heuristically proved on these nonlinear prototype tests, showing at the same time the very satisfactory level of accuracy reached. Through these tests, some other numerical specificities of the SPH method are discussed, such as the self‐redistribution of the particles occurring during the Lagrangian evolution. A higher order model is also proposed, and its advantages and drawbacks are discussed. Copyright © 2013 John Wiley & Sons, Ltd.     

Three‐dimensional wave impact on a rigid structure using smoothed particle hydrodynamics

Navier–Stokes computations of a wave–structure interaction are performed with the aim of assessing the potential of smoothed particle hydrodynamics to accurately estimate impact loading time history. A three‐dimensional dam‐break flow with a rectangular column located downstream is considered. The net force and impulse exerted on the column is monitored throughout the simulation with the results correlating well with existing experimental data. Initial and boundary conditions and numerical parameters are varied and their effect on the column load investigated. The column load is found to be most sensitive to the choice of boundary treatment. Copyright © 2011 John Wiley & Sons, Ltd.     

A comparative study of truly incompressible and weakly compressible SPH methods for free surface incompressible flows

In this paper, the performance of the incompressible SPH (ISPH) method and an improved weakly compressible SPH (IWCSPH) method for free surface incompressible flows are compared and analyzed. In both methods, the Navier–Stokes equations are solved, and no artificial viscosity is used. The ISPH algorithm in this paper is based on the classical SPH projection method with common treatments on solid boundaries and free surfaces. The IWCSPH model includes some advanced corrective algorithms in density approximation and solid boundary treatment (SBT). In density approximation, the moving least squares (MLS) approach is applied to re‐initialize density every several steps to obtain smoother and more stable pressure fields. An improved coupled dynamic SBT algorithm is implemented to obtain stable pressure values near solid wall areas and, thus, to minimize possible numerical oscillations brought in by the solid boundaries. Three representative numerical examples, including a benchmark test for hydrostatic pressure, a dam breaking problem and a liquid sloshing problem, are comparatively analyzed with ISPH and IWCSPH. It is demonstrated that the present IWCSPH is more attractive than ISPH in modeling free surface incompressible flows as it is more accurate and more stable with comparable or even less computational efforts. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical simulation of soil–water interaction using smoothed particle hydrodynamics (SPH) method

An application of smoothed particle hydrodynamics (SPH) to simulation of soil–water interaction is presented. In this calculation, water is modeled as a viscous fluid with week compressibility and soil is modeled as an elastic–perfectly plastic material. The Mohr–Coulomb failure criterion is applied to describe the stress states of soil in the plastic flow regime. Dry soil is modeled by one-phase flow while saturated soil is modeled by separate water and soil phases. Interaction between soil and water is taken into account by means of pore water pressure and seepage force. Simulation tests of soil excavation by a water jet are calculated as a challenging example to verify the broad applicability of the SPH method. The excavations are carried out in two different soil models, one is dry soil and the other is fully saturated soil. Numerical results obtained in this paper have shown that the gross discontinuities of soil failure can be simulated without any difficulties. This supports the feasibility and attractiveness of this a new approach in geomechanics applications. Advantages of the method are robustness, conceptual simplicity and relative ease of incorporating new physics.     

基于拟流体模型的SPH新方法及其在弹丸超高速碰撞薄板中的应用

引入颗粒动力学理论(拟流体模型)建立了适用于超高速碰撞的SPH新方法。将超高速碰撞中处于损伤状态的碎片等效为拟流体,在描述其运动过程中引入了碎片间相互作用和气体相对碎片的作用。采用该方法对球形弹丸超高速碰撞薄板形成碎片云的过程进行了数值模拟,得到了弹坑直径、外泡碎片云和内核碎片云的形状、分布,并与使用传统SPH方法、自适应光滑粒子流体动力学(ASPH)方法的模拟结果进行对比,结果显示：新方法在内核碎片云形状和分布上计算结果更加准确。同时对Whipple屏超高速碰撞问题进行了研究,分析了不同撞击速度下防护屏弹坑尺寸及舱壁损伤特性等特性,计算结果与实验吻合较好且符合Whipple防护结构的典型撞击极限曲线。

     

Time‐dependent flow across a step: the slip with friction boundary condition

The paper studies numerically the slip with friction boundary condition in the time‐dependent incompressible Navier–Stokes equations. Numerical tests on two‐ and three‐dimensional channel flows across a step using this boundary condition on the bottom wall are performed. The influence of the friction parameter on the flow field is studied and the results are explained according to the physics of the flow. Due to the stretching and tilting of vortices, the three‐dimensional results differ in many respects from the two‐dimensional ones. Copyright © 2005 John Wiley & Sons, Ltd.     

Unified semi‐analytical wall boundary conditions for inviscid,laminar or turbulent flows in the meshless SPH method

Wall boundary conditions in smoothed particle hydrodynamics (SPH) is a key issue to perform accurate simulations. We propose here a new approach based on a renormalising factor for writing all boundary terms. This factor depends on the local shape of a wall and on the position of a particle relative to the wall, which is described by segments (in two‐dimensions), instead of the cumbersome fictitious or ghost particles used in most existing SPH models. By solving a dynamic equation for the renormalising factor, we significantly improve traditional wall treatment in SPH, for pressure forces, wall friction and turbulent conditions. The new model is demonstrated for cases including hydrostatic conditions for still water in a tank of complex geometry and a dam break over triangular bed profile with sharp angle where significant improved behaviour is obtained in comparison with the conventional boundary techniques. The latter case is also compared with a finite volume and volume‐of‐fluid scheme. The performance of the model for a two‐dimensional laminar flow in a channel is demonstrated where the profiles of velocity are in agreement with the theoretical ones, demonstrating that the derived wall shear stress balances the pressure gradient. Finally, the performance of the model is demonstrated for flow in a schematic fish pass where both the velocity field and turbulent viscosity fields are satisfactorily reproduced compared with mesh‐based codes. Copyright © 2012 John Wiley & Sons, Ltd.     

A method for computing compressible,highly stratified flows in astrophysics based on operator splitting

A new numerical approach based on consistent operator splitting is presented for computing compressible, highly stratified flows in astrophysics. The algorithm is particularly designed to search for steady or almost steady solutions for the time-dependent Navier–Stokes equations, describing viscous flow under the influence of a strong gravitational field. The algorithm proposed is multidimensional and works in Cartesian, cylindrical or spherical co-ordinates. It uses a second-order finite volume scheme with third-order upwinding and a second-order time discretization. An adaptive time step control and monotonic multilevel grid distribution has been incorporated to speed up convergence. This method has been incorporated into a hydrodynamical code by which, for the first time, for two-dimensional models the dynamics of the boundary layer in the accretion disk around a compact star could be computed over the whole viscous time scale. © 1998 John Wiley & Sons, Ltd.     

On the stability and dissipation of wall boundary conditions for compressible flows

Characteristic formulations for boundary conditions have demonstrated their effectiveness to handle inlets and outlets, especially to avoid acoustic wave reflections. At walls, however, most authors use simple Dirichlet or Neumann boundary conditions, where the normal velocity (or pressure gradient) is set to zero. This paper demonstrates that there are significant differences between characteristic and Dirichlet methods at a wall and that simulations are more stable when using walls modelled with a characteristic wave decomposition. The derivation of characteristic methods yields an additional boundary term in the continuity equation, which explains their increased stability. This term also allows to handle the two acoustic waves going towards and away from the wall in a consistent manner. Those observations are confirmed by stability matrix analysis and one‐ and two‐dimensional simulations of acoustic modes in cavities. Copyright © 2009 John Wiley & Sons, Ltd.     

A multi‐energy‐level lattice Boltzmann model for the compressible Navier–Stokes equations

In this paper, we propose a new lattice Boltzmann model for the compressible Navier–Stokes equations. The new model is based on a three‐energy‐level and three‐speed lattice Boltzmann equation by using a method of higher moments of the equilibrium distribution functions. As the 25‐bit model, we obtained the equilibrium distribution functions and the compressible Navier–Stokes equations with the second accuracy of the truncation errors. The numerical examples show that the model can be used to simulate the shock waves, contact discontinuities and supersonic flows around circular cylinder. The numerical results are compared with those obtained by traditional method. Copyright © 2007 John Wiley & Sons, Ltd.     

Efficient computation of the flow around single fluid particles using an artificial boundary condition

Two‐phase flows around fluid particles are often considered to be in infinite domains, to avoid influence of the domain walls. Numerical simulations, however, must be modeled with a bounded domain, thus introducing artificial boundaries. Modeling of fluid flow in a domain with such artificial boundaries requires a careful choice of suitable boundary conditions. Slip boundary conditions for example can have a large impact on the computational results if the domain is chosen to be too small, because they model impermeable walls. This paper introduces an artificial boundary condition for simulations of the flow around single rising or settling fluid particles based on the approximated decay behavior of the velocity and the pressure field in the surrounding liquid. This is applied to the simulation of rising gas bubbles in systems with a Reynolds number of up to 50, and the outcome is compared with experimental results and simulations with slip boundary condition. It is found that domain size can be reduced by a factor of about two compared with slip boundary conditions without loss of accuracy. Copyright © 2014 John Wiley & Sons, Ltd.     

A promising boundary element formulation for three‐dimensional viscous flow

In this paper, a new set of boundary‐domain integral equations is derived from the continuity and momentum equations for three‐dimensional viscous flows. The primary variables involved in these integral equations are velocity, traction, and pressure. The final system of equations entering the iteration procedure only involves velocities and tractions as unknowns. In the use of the continuity equation, a complex‐variable technique is used to compute the divergence of velocity for internal points, while the traction‐recovery method is adopted for boundary points. Although the derived equations are valid for steady, unsteady, compressible, and incompressible problems, the numerical implementation is only focused on steady incompressible flows. Two commonly cited numerical examples and one practical pipe flow problem are presented to validate the derived equations. Copyright © 2004 John Wiley & Sons, Ltd.     

No‐slip consistent immersed boundary particle tracking to simulate impaction filtration in porous media

In this paper, we present a new method for simulating the motion of a disperse particle phase in a carrier gas through porous media. We assume a sufficiently dilute particle‐laden flow and compute, independently of the disperse phase, the steady laminar fluid velocity using the immersed boundary method. Given the velocity of the carrier gas, the equations of motion for the particles experiencing the Stokes drag force are solved to determine their trajectories. The ‘no‐slip consistent’ particle tracking algorithm avoids possible numerical filtration of very small particles due to the nonzero velocity field at the solid–fluid interface introduced by the immersed boundary method. This physically consistent tracking allows a reliable estimation of the filtration efficiency of porous filters due to inertial impaction. We illustrate and test our new approach for model porous media consisting of a structured array of aligned rectangular fibers, arranged in line and staggered. In the staggered geometry, the effect of the residual velocity at the solid–fluid interface is significant for particles with low inertia. Without adopting the developed no‐slip consistent numerical method, an artificial numerical filtration is observed, which becomes dominant for small enough particles. For both the in line and the staggered geometries, the filtration rate depends quite strongly and non monotonically on the particle inertia. This is expressed most clearly in the staggered arrangement in which a very strong increase in the filtration efficiency is observed at a well‐defined critical droplet size, corresponding to a qualitative change in the dominant particle paths in the porous medium. Copyright © 2013 John Wiley & Sons, Ltd.