SPH computation of plunging waves using a 2‐D sub‐particle scale (SPS) turbulence model

The paper presents a 2‐D large eddy simulation (LES) modelling approach to investigate the properties of the plunging waves. The numerical model is based on the smoothed particle hydrodynamics (SPH) method. SPH is a mesh‐free Lagrangian particle approach which is capable of tracking the free surfaces of large deformation in an easy and accurate way. The Smagorinsky model is used as the turbulence model due to its simplicity and effectiveness. The proposed 2‐D SPH–LES model is applied to a cnoidal wave breaking and plunging over a mild slope. The computations are in good agreement with the documented data. Especially the computed turbulence quantities under the breaking waves agree better with the experiments as compared with the numerical results obtained by using the k–ε model. The sensitivity analyses of the SPH–LES computations indicate that both the turbulence model and the spatial resolution play an important role in the model predictions and the contributions from the sub‐particle scale (SPS) turbulence decrease with the particle size refinement. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical study of breaking waves by a two‐phase flow model

A two‐phase flow model, which solves the flow in the air and water simultaneously, is presented for modelling breaking waves in deep and shallow water, including wave pre‐breaking, overturning and post‐breaking processes. The model is based on the Reynolds‐averaged Navier–Stokes equations with the k ?ε turbulence model. The governing equations are solved by the finite volume method in a Cartesian staggered grid and the partial cell treatment is implemented to deal with complex geometries. The SIMPLE algorithm is utilised for the pressure‐velocity coupling and the air‐water interface is modelled by the interface capturing method via a high resolution volume of fluid scheme. The numerical model is validated by simulating overturning waves on a sloping beach and over a reef, and deep‐water breaking waves in a periodic domain, in which good agreement between numerical results and available experimental measurements for the water surface profiles during wave overturning is obtained. The overturning jet, air entrainment and splash‐up during wave breaking have been captured by the two‐phase flow model, which demonstrates the capability of the model to simulate free surface flow and wave breaking problems.Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical modelling of complex turbulent free‐surface flows with the SPH method: an overview

The gridless smoothed particle hydrodynamics (SPH) method is now commonly used in computational fluid dynamics (CFD) and appears to be promising in predicting complex free‐surface flows. However, increasing flow complexity requires appropriate approaches for taking account of turbulent effects, whereas some authors are still working without any turbulence closure in SPH. A review of recently developed turbulence models adapted to the SPH method is presented herein, from the simplistic point of view of a one‐equation model involving mixing length to more sophisticated (and thus realistic) models like explicit algebraic Reynolds stress models (EARSM) or large eddy simulation (LES). Each proposed model is tested and validated on the basis of schematic cases for which laboratory data, theoretical or numerical solutions are available in the general field of turbulent free‐surface incompressible flows (e.g. open‐channel flow and schematic dam break). They give satisfactory results, even though some progress should be made in the future in terms of free‐surface influence and wall conditions. Recommendations are given to SPH users to apply this method to the modelling of complex free‐surface turbulent flows. Copyright © 2006 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 landslide impulsive waves by incompressible smoothed particle hydrodynamics

An incompressible‐smoothed particle hydrodynamics (I‐SPH) formulation is presented to simulate impulsive waves generated by landslides. The governing equations, Navier–Stokes equations, are solved in a Lagrangian form using a two‐step fractional method. Landslides in this paper are simulated by a submerged mass sliding along an inclined plane. During sliding, both rigid and deformable landslides mass are considered. The present numerical method is examined for a rigid wedge sliding into water along an inclined plane. In addition solitary wave generated by a heavy box falling inside water, known as Scott Russell wave generator, which is an example for simulating falling rock avalanche into artificial and natural reservoirs, is simulated and compared with experimental results. The numerical model is also validated for gravel mass sliding along an inclined plane. The sliding mass approximately behaves like a non‐Newtonian fluid. A rheological model, implemented as a combination of the Bingham and the general Cross models, is utilized for simulation of the landslide behaviour. In order to match the experimental data with the computed wave profiles generated by deformable landslides, parameters of the rheological model are adjusted and the numerical model results effectively match the experimental results. The results prove the efficiency and applicability of the I‐SPH method for simulation of these kinds of complex free surface problems. Copyright © 2007 John Wiley & Sons, Ltd.     

An improved KGF‐SPH with a novel discrete scheme of Laplacian operator for viscous incompressible fluid flows

The kernel gradient free (KGF) smoothed particle hydrodynamics (SPH) method is a modified finite particle method (FPM) which has higher order accuracy than the conventional SPH method. In KGF‐SPH, no kernel gradient is required in the whole computation, and this leads to good flexibility in the selection of smoothing functions and it is also associated with a symmetric corrective matrix. When modeling viscous incompressible flows with SPH, FPM or KGF‐SPH, it is usual to approximate the Laplacian term with nested approximation on velocity, and this may introduce numerical errors from the nested approximation, and also cause difficulties in dealing with boundary conditions. In this paper, an improved KGF‐SPH method is presented for modeling viscous, incompressible fluid flows with a novel discrete scheme of Laplacian operator. The improved KGF‐SPH method avoids nested approximation of first order derivatives, and keeps the good feature of ‘kernel gradient free’. The two‐dimensional incompressible fluid flow of shear cavity, both in Euler frame and Lagrangian frame, are simulated by SPH, FPM, the original KGF‐SPH and improved KGF‐SPH. The numerical results show that the improved KGF‐SPH with the novel discrete scheme of Laplacian operator are more accurate than SPH, and more stable than FPM and the original KGF‐SPH. Copyright © 2015 John Wiley & Sons, Ltd.     

Incompressible SPH simulation of water entry of a free‐falling object

In this paper an incompressible smoothed particle hydrodynamics (Incom‐SPH) model is used to simulate the interactions between the free surface flow and a moving object. Incom‐SPH method is a two‐step semi‐implicit hydrodynamic formulation of the SPH algorithm and is capable of accurately treating the free surface deformations and impact forces during the solid–fluid interactions. For a free‐falling object, its motion is tracked by an additional Lagrangian algorithm based on Newton's law to couple with the Incom‐SPH program. The developed model is employed to investigate the water entry of a free‐falling wedge. The accuracy of the computations is validated by the good agreement in predicting the relevant hydrokinematic and hydrodynamic parameters. Finally, a numerical test is performed to study the influence of spatial resolution on the water entry features. The Incom‐SPH modeling coupled with the solid–fluid interaction algorithm could provide a promising computational tool to predict the slamming problems in coastal and offshore engineering. Copyright © 2008 John Wiley & Sons, Ltd.     

海堤越浪的数值模拟

基于RANS方程和两方程湍流模型,采用有限体积法,将人射波波场作为人工的分布源项加人动量方程,提出了适用于VOF方法的源造波一消波技术。通过对行波及驻波的计算,分别考察了数值波浪水槽前端及末端消波段的有效性。在本文建立的数值波浪水槽内对规则波在海堤上爬高和越浪过程进行了数值模拟,并将计算结果与现有实验结果进行了比较。验证计算结果表明,数值模拟结果较好地复演了海堤越浪过程。为了研究模型尺度对越浪量的影响,文中设计了一组满足重力相似但具有不同几何比尺的数值实验模型。系列数值实验结果表明,若按重力相似换算越浪量,计算结果与实验预报值间的偏差随模型比尺的增大和堤前波浪破碎强度的增强而增大,建议在进行越浪物理模型实验时需进一步考虑模型比尺对原型预报值的影响。     

A 2D+t SPH model to study the breaking wave pattern generated by fast ships

A 2D+t approach is applied to study the wave pattern generated by high speed slender ships with a sharp stem. This allows approximating the ship motion through a set of equations which are mathematically equivalent to those governing the unsteady 2D free-surface flow generated by a deformable body in the vertical plane transverse to the ship. To describe the body deformation, a proper modeling of the solid boundaries is needed. To this purpose the fixed ghost particle technique recently developed by the authors for 2D SPH schemes is adopted. Then, an in-depth analysis of the bow breaking wave evolution is provided focusing on the different features of the plunging jet as function of the ship forward velocity. The comparison with experimental data proves the proposed SPH scheme to be robust and accurate.     

Improvement of moving particle semi‐implicit method for simulation of progressive water waves

Precise simulation of the propagation of surface water waves, especially when involving breaking wave, takes a significant place in computational fluid dynamics. Because of the strong nonlinear properties, the treatment of large surface deformation of free surface flow has always been a challenging work in the development of numerical models. In this paper, the moving particle semi‐implicit (MPS) method, an entirely Lagrangian method, is modified to simulate wave motion in a 2‐D numerical wave flume preferably. In terms of consecutive pressure distribution, a new and simple free surface detection criterion is proposed to enhance the free surface recognition in the MPS method. In addition, a revised gradient model is deduced to diminish the effect of nonuniform particle distribution and then to reduce the numerical wave attenuation occurring in the original MPS model. The applicability and stability of the improved MPS method are firstly demonstrated by the calculation of hydrostatic problem. It is revealed that these modifications are effective to suppress the pressure oscillation, weaken the local particle clustering, and boost the stability of numerical algorithm. It is then applied to investigate the propagation of progressive waves on a flat bed and the wave breaking on a mild slope. Comparisons with the analytical solutions and experimental results indicate that the improved MPS model can give better results about the profiles and heights of surface waves in contrast with the previous MPS models. Copyright © 2017 John Wiley & Sons, Ltd.     

Finite volume computation of incompressible turbulent flows in general co-ordinates on staggered grids

A brief review of the computation of incompressible turbulent flow in complex geometries is given. A 2D finite volume method for the calculation of turbulent flow in general curvilinear co-ordinates is described. This method is based on a staggered grid arrangement and the contravariant flux componets are chosen as primitive variables. Turbulence is modelled either by the standard k–ε model or by a k–ε model based on RNG theory. Convection is approximated with central differences for the mean flow quantities and a TVD-type MUSCL scheme for the turbulence equations. The sensitivity of the method to the grid properties is investigated. An application of this method to a complex turbulent flow is presented. The results of computations are compared with experimental data and other numerical solutions and are found to be satisfactory.     

Numerical analysis of breaking waves using the moving particle semi-implicit method

The numerical method used in this study is the moving particle semi-implicit (MPS) method, which is based on particles and their interactions. The particle number density is implicitly required to be constant to satisfy incompressibility. A semi-implicit algorithm is used for two-dimensional incompressible non-viscous flow analysis. The particles whose particle number densities are below a set point are considered as on the free surface. Grids are not necessary in any calculation steps. It is estimated that most of computation time is used in generation of the list of neighboring particles in a large problem. An algorithm to enhance the computation speed is proposed. The MPS method is applied to numerical simulation of breaking waves on slopes. Two types of breaking waves, plunging and spilling breakers, are observed in the calculation results. The breaker types are classified by using the minimum angular momentum at the wave front. The surf similarity parameter which separates the types agrees well with references. Breaking waves are also calculated with a passively moving float which is modelled by particles. Artificial friction due to the disturbed motion of particles causes errors in the flow velocity distribution which is shown in comparison with the theoretical solution of a cnoidal wave. © 1998 John Wiley & Sons, Ltd.     

An SPH model for free surface flows with moving rigid objects

This paper presents a computational model for free surface flows interacting with moving rigid bodies. The model is based on the SPH method, which is a popular meshfree, Lagrangian particle method and can naturally treat large flow deformation and moving features without any interface/surface capture or tracking algorithm. Fluid particles are used to model the free surface flows which are governed by Navier–Stokes equations, and solid particles are used to model the dynamic movement (translation and rotation) of moving rigid objects. The interaction of the neighboring fluid and solid particles renders the fluid–solid interaction and the non‐slip solid boundary conditions. The SPH method is improved with corrections on the SPH kernel and kernel gradients, enhancement of solid boundary condition, and implementation of Reynolds‐averaged Navier–Stokes turbulence model. Three numerical examples including the water exit of a cylinder, the sinking of a submerged cylinder and the complicated motion of an elliptical cylinder near free surface are provided. The obtained numerical results show good agreement with results from other sources and clearly demonstrate the effectiveness of the presented meshfree particle model in modeling free surface flows with moving objects. Copyright © 2013 John Wiley & Sons, Ltd.     

Conservation of circulation in SPH for 2D free‐surface flows

In this paper, we study how accurately the Smoothed Particle Hydrodynamics (SPH) scheme accounts for the conservation and the generation of vorticity and circulation, in a low viscosity, weakly compressible, barotropic fluid in the context of free‐surface flows. We consider a number of simple examples to clarify the processes involved and the accuracy of the simulations. The first example is a differentially rotating fluid where the integration path for the circulation becomes progressively more complicated, whereas the structure of the velocity field remains simple. The second example is the collision of two rectangular regions of fluid. We show that SPH accurately predicts the time variation of the circulation as well as the total vorticity for selected domains advected by the fluid. Finally, a breaking wave is considered. For such a problem we show how the dynamics of the vorticity generated by the breaking process is captured by the SPH model. Copyright © 2012 John Wiley & Sons, Ltd.     

A level set technique applied to unsteady free surface flows

An unsteady Navier–Stokes solver for incompressible fluid is coupled with a level set approach to describe free surface motions. The two‐phase flow of air and water is approximated by the flow of a single fluid whose properties, such as density and viscosity, change across the interface. The free surface location is captured as the zero level of a distance function convected by the flow field. To validate the numerical procedure, two classical two‐dimensional free surface problems in hydrodynamics, namely the oscillating flow in a tank and the waves generated by the flow over a bottom bump, are studied in non‐breaking conditions, and the results are compared with those obtained with other numerical approaches. To check the capability of the method in dealing with complex free surface configurations, the breaking regime produced by the flow over a high bump is analyzed. The analysis covers the successive stages of the breaking phenomenon: the steep wave evolution, the falling jet, the splash‐up and the air entrainment. In all phases, numerical results qualitatively agree with the experimental observations. Finally, to investigate a flow in which viscous effects are relevant, the numerical scheme is applied to study the wavy flow past a submerged hydrofoil. Copyright © 2001 John Wiley & Sons, Ltd.     

An efficient numerical tank for non-linear water waves,based on the multi-subdomain approach with BEM

An efficient 2D non-linear numerical wave tank called LONGTANK has been developed based on a multi-subdomain (MSD) approach combined with the conventional boundary element method (BEM). The multi-subdomain approach aims at optimized matrix diagonalization, thus minimizing the computing time and reserved storage. The CPU per time step in LONGTANK simulation is found to increase only linearly with the number of surface nodes, which makes LONGTANK highly efficient especially when simulating long-time wave evolutions in space. Appropriate treatment of special points on the boundary ensures high resolution in LONGTANK simulation beyond initial deformation and breaking, which allows detailed study of breaking criterion, breaker morphology, breaking dissipation, vorticity generation, etc. Detailed numerical implementation has been given with demonstration of LONGTANK simulations.     

A non‐linear k–ε model with realizability for prediction of flows around bluff bodies

The incompressible flow around bluff bodies (a square cylinder and a cube) is investigated numerically using turbulence models. A non‐linear k–ε model, which can take into account the anisotropy of turbulence with less CPU time and computer memory then RSM or LES, is adopted as a turbulence model. In tuning of the model coefficients of the non‐linear terms are adjusted through the examination of previous experimental studies in simple shear flows. For the tuning of the coefficient in the eddy viscosity (=C_μ), the realizability constraints are derived in three types of basic 2D flow patterns, namely, a simple shear flow, flow around a saddle and a focal point. C_μ is then determined as a function of the strain and rotation parameters to satisfy the realizability. The turbulence model is first applied to a 2D flow around a square cylinder and the model performance for unsteady flows is examined focussing on the period and the amplitude of the flow oscillation induced by Karman vortex shedding. The applicability of the model to 3D flows is examined through the computation of the flow around a surface‐mounted cubic obstacle. The numerical results show that the present model performs satisfactorily to reproduce complex turbulent flows around bluff bodies. Copyright © 2003 John Wiley & Sons, Ltd.     

Calculation of three‐dimensional turbulent flow with a finite volume multigrid method

A numerical method for the efficient calculation of three‐dimensional incompressible turbulent flow in curvilinear co‐ordinates is presented. The mathematical model consists of the Reynolds averaged Navier–Stokes equations and the k–ε turbulence model. The numerical method is based on the SIMPLE pressure‐correction algorithm with finite volume discretization in curvilinear co‐ordinates. To accelerate the convergence of the solution method a full approximation scheme‐full multigrid (FAS‐FMG) method is utilized. The solution of the k–ε transport equations is embedded in the multigrid iteration. The improved convergence characteristic of the multigrid method is demonstrated by means of several calculations of three‐dimensional flow cases. Copyright © 1999 John Wiley & Sons, Ltd.     

Optimized incompressible smoothed particle hydrodynamics methods and validations

The solution of the Poisson's equation used by the incompressible smoothed particle hydrodynamics (ISPH) methods for estimating the pressure field is expensive in CPU time. The CPU time, consumed by the inversion of the operator ∇(1/ρ∇) and the estimation of the right hand side of the Poisson's equation, increases with the number N of particles used in a purely Lagrangian framework. In this work, this default of ISPH methods is overcome by solving the Poisson's equation on a Cartesian grid. This SPH-mesh coupling is equivalent to the particle in cell method. In a first step, in order to analyze its efficiency, the optimized version of two ISPH methods (divergence free and density invariant) is compared with the standard weakly compressible SPH method through two benchmarks of incompressible bidimensional flows characterized by the Reynolds number Re, Lamb-Oseen vortex (10 ≤Re≤ 100) and lid-driven cavity flow (100 ≤Re≤ 1000). In a second step, the numerical results obtained by the three SPH methods are compared to laboratory experimental data of a dam break flow in order to show the performance of the SPH-mesh coupling in a practical and complex flow problem. As in the configuration of the experimental setup, the numerical results are obtained for a Reynolds number Re = 3.8 10⁶.