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.     

A corrected symmetric SPH method to simulate viscoelastic free surface flows based on the PTT model

In this work, a corrected symmetric and periodic density reinitialized SPH (CSPDR‐SPH) method is proposed and extended to simulate the viscoelastic free surface flows based on the Phan–Thien–Tanner model. The improvements mainly lie in deriving a corrected symmetric kernel gradient, and combining it with a periodic density reinitialization procedure. In addition, a simple artificial viscosity and a simple artificial stress form are adopted. Thus, the CSPDR‐SPH method has higher accuracy and better stability than the SPH method, and conserves both linear and angular momentums. The consistency and convergence of the CSPDR‐SPH method are justified by approximating a function in one and two dimensions. The merits of CSPDR‐SPH method are demonstrated by several benchmarks. The simple flow in a two‐dimensional channel is investigated to show the capability of the CSPDR‐SPH method to simulate the viscoelastic free surface flow. Then the CSPDR‐SPH method is extended to simulate the impacting drop problem. Numerical results show that the CSPDR‐SPH method can precisely capture the viscoelastic free surface. The Reynolds number, Weissenberg number and elongation parameter have remarkable effect on the flows. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

Computation of turbulent free‐surface flows around modern ships

This paper presents the calculated results for three classes of typical modern ships in modelling of ship‐generated waves. Simulations of turbulent free‐surface flows around ships are performed in a numerical water tank, based on the FINFLO‐RANS SHIP solver developed at Helsinki University of Technology. The Reynolds‐averaged Navier–Stokes (RANS) equations with the artificial compressibility and the non‐linear free‐surface boundary conditions are discretized by means of a cell‐centred finite‐volume scheme. The convergence performance is improved with the multigrid method. A free surface is tracked using a moving mesh technology, in which the non‐linear free‐surface boundary conditions are given on the actual location of the free surface. Test cases recommended are a container ship, a US Navy combatant and a tanker. The calculated results are compared with the experimental data available in the literature in terms of the wave profiles, wave pattern, and turbulent flow fields for two turbulence models, Chien's low Reynolds number k–εmodel and Baldwin–Lomax's model. Furthermore, the convergence performance, the grid refinement study and the effect of turbulence models on the waves have been investigated. Additionally, comparison of two types of the dynamic free‐surface boundary conditions is made. Copyright © 2003 John Wiley& Sons, Ltd.     

Introduction of turbulent model in a mixed finite volume/finite element method

The aim of this work is to present a new numerical method to compute turbulent flows in complex configurations. With this in view, a k-? model with wall functions has been introduced in a mixed finite volume/finite element method. The numerical method has been developed to deal with compressible flows but is also able to compute nearly incompressible flows. The physical model and the numerical method are first described, then validation results for an incompressible flow over a backward-facing step and for a supersonic flow over a compression ramp are presented. Comparisons are performed with experimental data and with other numerical results. These simulations show the ability of the present method to predict turbulent flows, and this method will be applied to simulate complex industrial flows (flow inside the combustion chamber of gas turbine engines). The main goal of this paper is not to test turbulence models, but to show that this numerical method is a solid base to introduce more sophisticated turbulence model.     

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.     

Nonlinear turbulence models for predicting strong curvature effects

Prediction of the characteristics of turbulent flows with strong streamline curvature, such as flows in turbomachines, curved channel flows, flows around airfoils and buildings, is of great importance in engineering applications and poses a very practical challenge for turbulence modeling. In this paper, we analyze qualitatively the curvature effects on the structure of turbulence and conduct numerical simulations of a turbulent Uduct flow with a number of turbulence models in order to assess their overall performance. The models evaluated in this work are some typical linear eddy viscosity turbulence models, nonlinear eddy viscosity turbulence models （NLEVM） （quadratic and cubic）, a quadratic explicit algebraic stress model （EASM） and a Reynolds stress model （RSM） developed based on the second-moment closure. Our numerical results show that a cubic NLEVM that performs considerably well in other benchmark turbulent flows, such as the Craft, Launder and Suga model and the Huang and Ma model, is able to capture the major features of the highly curved turbulent U-duct flow, including the damping of turbulence near the convex wall, the enhancement of turbulence near the concave wall, and the subsequent turbulent flow separation. The predictions of the cubic models are quite close to that of the RSM, in relatively good agreement with the experimental data, which suggests that these models may be employed to simulate the turbulent curved flows in engineering applications.     

Incompressible SPH simulation of wave breaking and overtopping with turbulence modelling

In this paper a truly incompressible version of the smoothed particle hydrodynamics (SPH) method is presented to investigate the surface wave overtopping. SPH is a pure Lagrangian approach which can handle large deformations of the free surface with high accuracy. The governing equations are solved based on the SPH particle interaction models and the incompressible algorithm of pressure projection is implemented by enforcing the constant particle density. The two‐equation k–ε model is an effective way of dealing with the turbulence and vortices during wave breaking and overtopping and it is coupled with the incompressible SPH numerical scheme. The SPH model is employed to reproduce the experiment and computations of wave overtopping of a sloping sea wall. The computations are validated against the experimental and numerical data found in the literatures and good agreement is observed. Besides, the convergence behaviour of the numerical scheme and the effects of particle spacing refinement and turbulence modelling on the simulation results are also investigated in further detail. The sensitivity of the computed wave breaking and overtopping on these issues is discussed and clarified. Copyright © 2005 John Wiley & Sons, Ltd.     

Turbulent liquid metal flow in rectangular shaped contraction nozzles for target applications

Owing to the high beam power densities envisaged in advanced nuclear targets, liquid metal-operated free surface targets are conceived as one feasible option. There, the free surface is formed by an adequately shaped upstream located nozzle. Target boundary conditions necessitate a detailed knowledge on the turbulent flow in contraction nozzles in order to identify turbulence models accurately predicting experimental findings within the velocity range of interest for nuclear target and hence can then act as design optimisation tools. In this context, a combined experimental and numerical study is conducted on the basis of the turbulent flow in the contraction nozzle of the Super-FRS target. Two aspects determining the turbulent flow in the nozzle have been investigated. The first is a potential relaminarisation of the boundary layer caused by the acceleration within the contraction and the second is a development of the secondary flows due to the pressure gradient in the rectangular shaped ducts cross-section. Regarding the three different turbulence models investigated here only the V2F model exhibited the capability to predict the relaminarisation of the turbulent boundary layer both qualitatively and quantitatively. All turbulence models are able to predict the development of secondary flows induced by pressure gradients in transverse direction with an acceptable accuracy.     

Numerical simulation of turbulent free surface flow with two‐equation k–ε eddy‐viscosity models

This paper presents a finite difference technique for solving incompressible turbulent free surface fluid flow problems. The closure of the time‐averaged Navier–Stokes equations is achieved by using the two‐equation eddy‐viscosity model: the high‐Reynolds k–ε (standard) model, with a time scale proposed by Durbin; and a low‐Reynolds number form of the standard k–ε model, similar to that proposed by Yang and Shih. In order to achieve an accurate discretization of the non‐linear terms, a second/third‐order upwinding technique is adopted. The computational method is validated by applying it to the flat plate boundary layer problem and to impinging jet flows. The method is then applied to a turbulent planar jet flow beneath and parallel to a free surface. Computations show that the high‐Reynolds k–ε model yields favourable predictions both of the zero‐pressure‐gradient turbulent boundary layer on a flat plate and jet impingement flows. However, the results using the low‐Reynolds number form of the k–ε model are somewhat unsatisfactory. Copyright © 2004 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.     

Parallel large eddy simulation of turbulent flow around MIRA model using linear equal‐order finite element method

A parallel large eddy simulation code that adopts domain decomposition method has been developed for large‐scale computation of turbulent flows around an arbitrarily shaped body. For the temporal integration of the unsteady incompressible Navier–Stokes equation, fractional 4‐step splitting algorithm is adopted, and for the modelling of small eddies in turbulent flows, the Smagorinsky model is used. For the parallelization of the code, METIS and Message Passing Interface Libraries are used, respectively, to partition the computational domain and to communicate data between processors. To validate the parallel architecture and to estimate its performance, a three‐dimensional laminar driven cavity flow inside a cubical enclosure has been solved. To validate the turbulence calculation, the turbulent channel flows at Re_τ = 180 and 1050 are simulated and compared with previous results. Then, a backward facing step flow is solved and compared with a DNS result for overall code validation. Finally, the turbulent flow around MIRA model at Re = 2.6 × 10⁶ is simulated by using approximately 6.7 million nodes. Scalability curve obtained from this simulation shows that scalable results are obtained. The calculated drag coefficient agrees better with the experimental result than those previously obtained by using two‐equation turbulence models. Copyright © 2007 John Wiley & Sons, Ltd.     

Effects of free‐stream turbulence on large‐scale coherent structures of separated boundary layer transition

It has been well established that large‐scale structures, usually called coherent structures, exist in many transitional and turbulent flows. The topology and range of scales of those large‐scale structures vary from flow to flow such as counter‐rotating vortices in wake flows, streaks and hairpin vortices in turbulent boundary layer. There has been relatively little study of large‐scale structures in separated and reattached transitional flows. Large‐eddy simulation (LES) is employed in the current study to investigate a separated boundary layer transition under 2% free‐stream turbulence on a flat plate with a blunt leading edge. The Reynolds number based on the inlet free stream velocity and the plate thickness is 6500. A dynamic subgrid‐scale model is employed to compute the subgrid‐scale stresses more accurately in the current transitional flow case. Flow visualization has shown that the Kelvin–Helmholtz rolls, which have been so clearly visible under no free‐stream turbulence (NFST) are not as apparent in the present study. The Lambda‐shaped vortical structures which can be clearly seen in the NFST case can hardly be identified in the free‐stream turbulence (FST) case. Generally speaking, the effects of free‐stream turbulence have led to an early breakdown of the boundary layer, and hence increased the randomization in the vortical structures, degraded the spanwise coherence of those large‐scale structures. Copyright © 2005 John Wiley & Sons, Ltd.     

无壁面参数低雷诺数非线性涡黏性模式研究

建立了一个低雷诺数的非线性涡黏性湍流模式,该模式的一个显著特性是它不包含壁面参数（如y^ ,n等）,因而特别适用于复杂几何流场的计算,本模式在几种包括回流、分离、激波等典型流动中进行了验证,结果令人满意。     

Computation of free surface waves around an arbitrary body by a Navier-Stokes solver using the psuedocompressibility technique

A Navier-Stokes equation solver is developed for computing free surface wave and viscous flow around an arbitrary body, in which a free surface model is introduced into the pseudocompressibility solution. The governing equations are classified in a vectorial form, with primitive variables, and a block diagonal system is generated by the discretization of an implicit factorization method. A moving grid system fitted to both the free surface and body surface is generated by an effective cubic spline fitting technique. Two zero-equation turbulence models, namely the Cebeci-Smith model and the Baldwin-Lomax model, are used for turbulent calculations. Numerical simulations are carried out for the free surface viscous flows generated by a submerged hydrofoil and a ship model. Computed results are in reasonable agreement with measurements.     

A gas‐kinetic scheme for the simulation of turbulent flows on unstructured meshes

This paper presents a numerical method for simulating turbulent flows via coupling the Boltzmann BGK equation with Spalart–Allmaras one equation turbulence model. Both the Boltzmann BGK equation and the turbulence model equation are carried out using the finite volume method on unstructured meshes, which is different from previous works on structured grid. The application of the gas‐kinetic scheme is extended to the simulation of turbulent flows with arbitrary geometries. The adaptive mesh refinement technique is also adopted to reduce the computational cost and improve the efficiency of meshes. To organize the unstructured mesh data structure efficiently, a non‐manifold hybrid mesh data structure is extended for polygonal cells. Numerical experiments are performed on incompressible flow over a smooth flat plate and compressible turbulent flows around a NACA 0012 airfoil using unstructured hybrid meshes. These numerical results are found to be in good agreement with experimental data and/or other numerical solutions, demonstrating the applicability of the proposed method to simulate both subsonic and transonic turbulent flows. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical analysis of turbulent structure in compound meandering open channel by algebraic Reynolds stress model

The turbulent flow in a compound meandering channel with a rectangular cross section is one of the most complicated turbulent flows, because the flow behaviour is influenced by several kinds of forces, including centrifugal forces, pressure‐driven forces and shear stresses generated by momentum transfer between the main channel and the flood plain. Numerical analysis has been performed for the fully developed turbulent flow in a compound meandering open‐channel flow using an algebraic Reynolds stress model. The boundary‐fitted coordinate system is introduced as a method for coordinate transformation in order to set the boundary conditions along the complicated shape of the meandering open channel. The turbulence model consists of transport equations for turbulent energy and dissipation, in conjunction with an algebraic stress model based on the Reynolds stress transport equations. With reference to the pressure–strain term, we have made use of a modified pressure–strain term. The boundary condition of the fluctuating vertical velocity is set to zero not only for the free surface, but also for computational grid points next to the free surface, because experimental results have shown that the fluctuating vertical velocity approaches zero near the free surface. In order to examine the validity of the present numerical method and the turbulent model, the calculated results are compared with experimental data measured by laser Doppler anemometer. In addition, the compound meandering open channel is clarified somewhat based on the calculated results. As a result of the analysis, the present algebraic Reynolds stress model is shown to be able to reasonably predict the turbulent flow in a compound meandering open channel. Copyright © 2005 John Wiley & Sons, Ltd.     

An implementation of the Spalart–Allmaras DES model in an implicit unstructured hybrid finite volume/element solver for incompressible turbulent flow

In this paper, we describe an implicit hybrid finite volume (FV)/element (FE) incompressible Navier–Stokes solver for turbulent flows based on the Spalart–Allmaras detached eddy simulation (SA‐DES). The hybrid FV/FE solver is based on the segregated pressure correction or projection method. The intermediate velocity field is first obtained by solving the original momentum equations with the matrix‐free implicit cell‐centered FV method. The pressure Poisson equation is solved by the node‐based Galerkin FE method for an auxiliary variable. The auxiliary variable is closely related to the real pressure and is used to update the velocity field and the pressure field. We store the velocity components at cell centers and the auxiliary variable at vertices, making the current solver a staggered‐mesh scheme. The SA‐DES turbulence equation is solved after the velocity and the pressure fields have been updated at the end of each time step. The same matrix‐free FV method as the one used for momentum equations is used to solve the turbulence equation. The turbulence equation provides the eddy viscosity, which is added to the molecular viscosity when solving the momentum equation. In our implementation, we focus on the accuracy, efficiency and robustness of the SA‐DES model in a hybrid flow solver. This paper will address important implementation issues for high‐Reynolds number flows where highly stretched elements are typically used. In addition, some aspects of implementing the SA‐DES model will be described to ensure the robustness of the turbulence model. Several numerical examples including a turbulent flow past a flat plate and a high‐Reynolds number flow around a high angle‐of‐attack NACA0015 airfoil will be presented to demonstrate the accuracy and efficiency of our current implementation. Copyright © 2008 John Wiley & Sons, Ltd.     

A correction for balancing discontinuous bed slopes in two‐dimensional smoothed particle hydrodynamics shallow water modeling

In this paper, a smoothed particle hydrodynamics (SPH) numerical model for the shallow water equations (SWEs) with bed slope source term balancing is presented. The solution of the SWEs using SPH is attractive being a conservative, mesh‐free, automatically adaptive method without special treatment for wet‐dry interfaces. Recently, the capability of the SPH–SWEs numerical scheme with shock capturing and general boundary conditions has been used for predicting practical flooding problems. The balance between the bed slope source term and fluxes in shallow water models is desirable for reliable simulations of flooding over bathymetries where discontinuities are present and has received some attention in the framework of Finite Volume Eulerian models. The imbalance because of the source term resulting from the calculation of the the water depth is eradicated by means of a corrected mass, which is able to remove the error introduced by a bottom discontinuity. Two different discretizations of the momentum equation are presented herein: the first one is based on the variational formulation of the SWEs in order to obtain a fully conservative formulation, whereas the second one is obtained using a non‐conservative form of the free‐surface elevation gradient. In both formulations, a variable smoothing length is considered. Results are presented demonstrating the corrections preserve still water in the vicinity of either 1D or 2D bed discontinuities and provide close agreement with 1D analytical solutions for rapidly varying flows over step changes in the bed. The method is finally applied to 2D dam break flow over a square obstacle where the balanced formulation improves the agreement with experimental measurements of the free surface. Copyright © 2012 John Wiley & Sons, Ltd.