改进的移动粒子半隐式法模拟楔形体入水砰击

移动粒子半隐式法(Moving-Particle Semi-Implicit Method,MPS)是一种新的基于拉格朗日(Lagrange)理念的无网格方法,适用于模拟自由液面的大变形和水流的喷射现象。用基于大涡模拟的改进MPS法首先模拟了矩形体的入水砰击,砰击压力的计算结果证明了这种方法的正确性,然后模拟了楔形体的匀速入水砰击,并与实验结果进行了对比,验证了大涡模拟改进MPS法在砰击问题中的适用性。     

Numerical Simulation of Sloshing Using the MPS-FSI Method with Large Eddy Simulation

A numerical model has been developed to study sloshing of turbulent flow in a tank with elastic baffles. The Moving-Particle Semi-implicit method(MPS) is a kind of meshless Lagrangian calculation method. The large eddy simulation(LES) approach is employed to model the turbulence by using the Smagorinsky Sub-Particle Scale(SPS)closure model. This paper uses MPS-FSI method with LES to simulate the interaction between free surface flow and a thin elastic baffle in sloshing. Then, the numerical model is validated, and the numerical solution has good agreement with experimental data for sloshing in a tank with elastic baffles. Furthermore, under external excitations,the MPS is applied to viscous laminar flow and turbulent flow, with both the deformation of elastic baffles and the wave height of the free surface are compared with each other. Besides, the impact pressure with/without baffles and wave height of free surface are investigated and discussed in detail. Finally, preliminary simulations are carried out in the damage problem of elastic baffles, taking the advantage of the MPS-FSI method in computations of the fluid–structure interaction with large deformation.     

波浪与抛石潜堤相互作用的MPS法数值模型研究

为了研究波浪与抛石潜堤相互作用过程中大自由表面变形和堤内渗流等强非线性紊流运动问题,利用改进的MPS法,建立了模拟波浪与抛石潜堤相互作用的MPS法数值计算模型。模型将抛石潜堤假定为均质多孔介质,采用Drew的二相流运动方程描述多孔介质内外的流体运动;通过在动量方程中增加非线性阻力项,并引入亚粒子尺度紊流模型,模拟波浪与可渗结构物相互作用过程中的紊流运动。选取“U”型管中多孔介质内渗流过程和孤立波与可渗潜堤相互作用两个典型的渗流问题,通过将数值计算结果与理论解和实测值的对比分析,对所提出的MPS法紊流渗流模型的模拟精度进行验证。结果表明:基于改进的MPS法构建的垂向二维紊流渗流模型可以很好地再现“U”型管中多孔介质内渗流以及波浪作用下可渗潜堤内外的复杂流场,显著缓解流-固界面处的压力震荡与粒子分布不均匀问题,实现了较高的模拟精度。     

MPS法高精度自由表面识别方法研究

为了提高粒子法即移动粒子半隐式法(Moving Particle Semi-implicit method, MPS法)中自由表面粒子的识别精度,降低由粒子误判引起的非物理压力振荡,对MPS法的自由表面识别方法进行了改进。在原始自由表面判别标准的基础上增加辅助判别条件,提出2种新的自由表面判别法即压力判别法和填充率判别法。利用对静水问题和溃坝流问题的模拟计算,对比分析选用不同自由表面判别法得到的数值计算结果,揭示粒子识别精度对压力计算的重要影响。研究结果表明:新提出的压力判别法和填充率判别法可以有效地提高自由表面粒子的识别精度,减轻压力计算中的非物理压力振荡现象,从而提高压力计算的稳定性以及整体数值计算的模拟精度。     

三维自由面流动模拟中GPU并行计算技术

MPS(Moving Particle Semi-implicit)法能够有效地处理溃坝、晃荡等自由面大变形流动问题。在三维MPS方法中,粒子数量的急剧增加会导致其计算效率的降低并限制其在大规模流动问题中的应用。基于自主开发的MPS求解器MLParticleSJTU,本文对求解过程中耗时最多的邻居粒子搜寻和泊松方程求解两个模块采用了GPU并行加速,详细探讨了CPU+GPU策略。以三维晃荡和三维溃坝这两种典型的自由面大变形流动为例,比较了CPU+GPU相对于MLParticle-SJTU串行求解时的加速情况,结果表明CPU+GPU在邻居粒子和泊松方程这两个模块中的加速比最高能达到十倍左右。此外,采用CPU+GPU并行能够较准确地模拟溃坝、晃荡等自由面大变形问题。     

Numerical simulation of wave interaction with porous structures using an improved smoothed particle hydrodynamic method

A smoothed particle hydrodynamic (SPH) model is developed to simulate wave interaction with porous structures. The mean flow outside the porous structures is obtained by solving Reynolds Averaged Navier–Stokes (RANS) equations and the turbulence field is calculated by a large eddy simulation (LES) model. The porous flow is described by the spatially averaged Navier–Stokes type equations with the resistance effect of the porous media being represented by an empirical frictional source term. The interface boundaries between the porous flow and the outside flow are modeled by means of specifying a transition zone along the interface. The model is validated against other available numerical results and experimental data for wave damping over porous seabed with different levels of permeability. The validated model is then employed to investigate wave breaking over a submerged porous breakwater and good agreements between the SPH model results and the experimental data are obtained in terms of free surface displacement. In addition the predicted velocity, vorticity and pressure fields near the porous breakwater and in the breaking wave zone are also analyzed.     

MPS方法数值模拟液舱晃荡问题

基于无网格粒子法MPS方法(moving particle semi-implicit method)研究了液舱晃荡问题。针对二维矩形液舱晃荡问题进行了数值验证,结果表明MPS方法能够很好地计算晃荡产生的拍击压力。同时将MPS方法应用到带隔板的液舱晃荡问题计算中,分析了二维和三维带隔板液舱晃荡问题。计算结果表明:隔板的存在很大程度地限制了流体的水平运动,隔板附近出现了自由面的翻卷、破碎和融合现象,MPS方法能够很好地模拟这些流动现象。计算得到的波高与实验测得的波高吻合较好,表明MPS方法模拟带隔板的晃荡问题具有一定的可靠性。     

多消浪室开孔沉箱消浪特性的数值模拟研究

本文基于雷诺平均的Navier-Stokes方程和k-ε模型求解湍流流动,采用流体体积法(Volume of Fluid,VOF)追踪自由表面运动,建立无反射波浪数值水槽,对多消浪室开孔沉箱的消浪特性进行数值模拟研究。将单消浪室和多消浪室开孔沉箱反射系数和结构前波面分布的数值分析结果与物理模型试验结果进行对比验证,两者符合良好。利用数值算例,研究多消浪室开孔沉箱的反射特性以及开孔结构附近的速度场和湍流强度分布。分析结果表明:波浪与开孔沉箱相互作用时,涡旋和湍动主要分布在开孔墙和消浪室内部自由表面附近;与单消浪室开孔沉箱相比,多消浪室开孔沉箱可以更有效的耗散波浪能量,降低结构的反射系数。本文分析结果可为开孔沉箱结构的工程设计提供参考依据。     

Numerical simulation of solitary and randomwave propagation through vegetation based on VOFmethod

A vertical two-dimensional numerical model has been applied to solving the Reynolds Averaged Navier- Stokes （RANS} equations in the simulation of current and wave propagation through vegetated and non- vegetated waters. The k-e model is used for turbulence closure of RANS equations. The effect of vegeta- tion is simulated by adding the drag force of vegetation in the flow momentum equations and turbulence model. To solve the modified N-S equations, the finite difference method is used with the staggered grid system to solver equations. The Youngs＇ fractional volume of fluid （VOF） is applied tracking the free sur- face with second-order accuracy. The model has been tested by simulating dam break wave, pure current with vegetation, solitary wave runup on vegetated and non-vegetated channel, regular and random waves over a vegetated field. The model reasonably well reproduces these experimental observations, the model- ing approach presented herein should be useful in simulating nearshore processes in coastal domains with vegetation effects.     

Investigation of the Effects of Baffles on the Shallow Water Sloshing in A Rectangular Tank Using A 2D Turbulent ISPH Method

Liquid sloshing is a common phenomenon in the liquid tanks transportation. Liquid waves lead to fluctuating forces on the tank wall. Uncontrolled fluctuations lead to large forces and momentums. Baffles can control these fluctuations. A numerical method, which has been widely used to model this phenomenon, is Smoothed Particle Hydrodynamics(SPH). The Lagrangian nature of this method makes it suitable for simulating free surface flows. In the present study, an accurate Incompressible Smoothed Particle Hydrodynamics(ISPH) method is developed and improved using the kernel gradient correction tensors, particle shifting algorithms, k–ε turbulence model, and free surface particle detectors. Comparisons with the experimental data approve the ability of the present algorithm for simulating shallow water sloshing. The main aim of this study is to investigate the effects of the vertical baffle on the damping of liquid sloshing. Results show that baffles number has a major role in sloshing fluctuation damping.     

N-S方程的数值解法及其在水波动力学中应用的综述

对于描述流体运动的基本方程N-S方程,现已开发出了多种不同的数值求解方法。其中以MAC法和VOF法为代表的网格数值方法较为成熟,并被逐渐用于实际工程的水动力研究中：以SPH和MPS为代表的粒子方法则刚刚起步．具有很大的发展空间。针对N-S方程中的能量耗散问题,采用雷诺时均化的方程形式仍是目前解决水动力问题的主要途径,但需要引入相应湍流模型,以封闭方程。本文首先对N-S方程求解方法的发展过程进行简要回顾,对一些重要的计算方法进行评述,简要介绍其在水波动力学中的应用;然后对时均化的N-S方程的主要封闭模式进行总结介绍;最后对数值求解N-S方程的发展趋势进行展望。     

Propagation of solitary waves over a submerged permeable breakwater

We study the interactions between a non-breaking solitary wave and a submerged permeable breakwater experimentally and numerically. The particle image velocimetry (PIV) technique is employed to measure instantaneous free surface displacements and velocity fields in the vicinity of a porous dike. The porous medium, consisting of uniform glass spheres, is mounted on the seafloor. Due to the limited size of each field of view (FOV) for high spatial resolution purposes, four FOVs are set in order to form a continuous flow field around the structure. Quantitative mean properties are obtained by ensemble averaging 30 repeated instantaneous measurements. The Reynolds decomposition method is then adopted to separate the velocity fluctuations for each trial to estimate the turbulent kinetic energy. In addition, a highly accurate two-dimensional model with the volume of fluid interface tracking technique is used to simulate an idealized volume-averaged porous medium. The model is based on the Volume-Averaged Reynolds Averaged Navier–Stokes equations coupled with the non-linear k–ε turbulence closure solver. Comparisons are performed between measurements and numerical results for the time histories of the free surface elevation recorded by wave gauges and the spatial distributions of free surface displacement with the corresponding velocity and turbulent kinetic energy around the permeable object imaged by the PIV system. Fairly good agreements are obtained. It is found that the measured and modeled turbulent intensities on the weather side are much larger than those on the lee side of the object, and that the magnitude of the turbulent intensity increases with increasing wave height of a solitary wave at a constant water depth. The verified numerical model is then used to estimate the energy reflection, transmission and dissipation using the energy integral method by varying the aspect ratio and the grain size of the permeable obstacle.     

Large-eddy simulation of the flow past a circular cylinder at sub- to super-critical Reynolds numbers

Large-eddy simulation of turbulent flow past a circular cylinder at sub- to super-critical Reynolds numbers is performed using a high-fidelity orthogonal curvilinear grid solver. Verification studies investigate the effects of grid resolution, aspect ratio and convection scheme. Monotonic convergence is achieved in grid convergence studies. Validation studies use all available experimental benchmark data. Although the grids are relatively large and fine enough for sufficiently resolved turbulence near the cylinder, the grid uncertainties are large indicating the need for even finer grids. Large aspect ratio is required for sub-critical Reynolds number cases, whereas small aspect ratio is sufficient for critical and super-critical Reynolds number cases. All the experimental trends were predicted with reasonable accuracy, in consideration the large facility bias, age of most of the data, and differences between experimental and computational setup in particular free stream turbulence and roughness. The largest errors were for under prediction of turbulence separation.     

Modified moving particle method for modeling wave interaction with multi layered porous structures

Modified Moving Particle method in Porous media (MMPP) is introduced in this study for simulating a flow interaction with porous structures. By making use of the sub-particle scale (SPS) turbulence model, a unified set of equations are introduced for the entire computational domain and a proper boundary treatment is suggested at the interfaces between fluid and the porous media. Similar to the Incompressible Smoothed Particle Hydrodynamic (ISPH) method, a robust two-step semi-implicit scheme is utilized to satisfy the incompressibility criterion. By means of the introduced model, different flow regimes through multi-layered porous structures with arbitrary shapes can be simulated and there is no need to implement calibration factors.The developed MMPP model is then validated via simulating the experiments of Liu et al. (1999) i.e. linear and turbulent flows through porous dams and the experiments of Sakakiyama and Liu (2001) i.e. wave overtopping on a caisson breakwater protected by multi layered porous materials. Good agreements between numerical and laboratory data present the ability of the introduced model in simulating various flow regimes through multi-layered porous structures. It is concluded that the turbulent flow is an important issue particularly at the interface between the free fluid and porous media and consequently, the accuracy of the previous Lagrangian models that were based on neglecting the turbulence effect can be improved significantly by means of the present model. In addition, to satisfy the continuity criteria in the SPH models, it is necessary to modify density of particles in accordance with their porosity.     

粗颗粒两相流计算方法研究

针对具有颗粒粒级跨度大且粒径大特点的固液两相流系统,如深海采矿系统等,开发新型粗颗粒—均质浆体计算模型,该模型采用欧拉—拉格朗日计算法,大粒径颗粒视为固相粗颗粒,其余视为细颗粒与海水组成液相均质浆体。固—液相间滑移速度v_(slip)作为粗细颗粒划分标准,通过联立颗粒雷诺数Re_p,液相对颗粒绕流阻力系数C_d,颗粒粒径d_p等参数,得到滑移速度v_(slip)与颗粒粒径d_p的相关函数,获得粗细颗粒的划分粒径,进而得到两相流控制方程。利用CFD软件,使用该方法对粗颗粒固—液两相流垂直管提升系统进行仿真计算,并将所获得的仿真结果与已发表成果进行对比,发现数据较为吻合,从一定程度上验证了该计算模型的准确性。     

Large-scale turbulence under a solitary wave

The structure of large-scale turbulence under a broken solitary wave on a 1 in 50 plane slope was studied. Three-component velocity measurements were taken at different heights above a smooth bed in the middle surf zone using an acoustic Doppler velocimeter. The measured data showed that turbulent velocity components were well correlated in the middle part of the water column. The velocity correlations could be produced by an oblique vortex similar to the obliquely descending eddy observed previously by other investigators. The vertical distributions of the relative values of the components of the Reynolds stress tensor showed that the structure of turbulence evolved continuously between the free surface and the bottom. The evolution was related to transition from two-dimensional to three-dimensional flow structures and the effect of the solid bottom on flow structures. Time histories of measured turbulent kinetic energy and turbulence stresses showed episodic turbulent events near the free surface but more sporadic turbulence in the lower layer. Large or intense turbulent events were found to have short duration and time lag relative to the wave crest point. These events also maintained good correlations between the turbulence velocity components close to the bottom.Instantaneous turbulent velocity fields were measured near the bottom at the same cross-shore location by using a stereoscopic particle image velocimetry system. These measurements showed that the near-bed flow field was characterized by large-scale, coherent flow structures that were the sources of most of the turbulent kinetic energy and turbulence stresses. The types of organized flow structures observed included vortices and downbursts of turbulence descending directly from above, lateral spreading of turbulent fluid along the bed, and formation of vortices in shear layers between fluid streams. A common feature of the organized flow structures near the bed was the large turbulence velocities in the longitudinal and transverse directions, which reflected the influence of a solid bottom on the breaking-wave-generated turbulence arriving at the bed.     

New modified gradient models for MPS method applied to free-surface flow simulations

In this paper, two modified pressure gradient models based on Taylor series expansion are proposed to enhance the higher order source term MPS (MPS-HS) method. The modified models consist of gradient correction matrices applied to the existing (base) pressure gradient models. To validate the modified pressure gradient models first hydrostatic pressure test is simulated and compared to both the base and modified MPS methods. Using the modified models are shown to reduce unphysical pressure oscillations observed in the base models. Second, an evolution of an elliptical drop in a 2D flow field is examined and shown to verify the models. Third, the proposed models illustrated appropriate stability and consistency properties against analytical solutions when an altered gravitational acceleration was superimposed to the hydrostatic pressure test. In addition, an improved performance is observed when Higher order Laplacian (HL) and Error-Compensating Source (ECS) of the Poisson Pressure Equation (PPE) schemes are coupled with the modified pressure gradient models compared to coupling them with the base gradient models. Finally, the modified MPS methods enhanced performances are validated in a free-surface flow simulation for a dam break problem with impact pressure, and a violent sloshing flow in a rectangular tank when compared to the base MPS methods against an existing experimental data.     

Modeling of Motion of Particle Clouds Formed by Dumping Dredged Material

The motion of particle clouds formed by dumping dredged material into quiescent waters is experimentally and numerically studied. In the numerical model, the particle phase is modeled by the dispersion model, and turbulence is calculated by the large eddy simulation. The governing equations, including the filtered Navier-Stokes equations and mass transport equation, are solved based on the operator-splitting algorithm and an implicit cubic spline interpolation scheme. The eddy viscosity is evaluated by the modified Smagorinsky model including the buoyancy term. Comparisons of main flow characteristics, including shape, size, average density excess, moving speed and the amount of particles deposited on the bed, between experimental and computational results show that the numerical model well predicts the motion of the cloud from the falling to spreading stage. The effects of silt-fence on the motion of the particle cloud are also investigated.     

A Numerical Study on the Structure and Dynamics of Turbulence in Oscillatory Bottom Layer Flow over A Flat or Rippled Bed

1 .IntroductionUnder the influence of surface waves ,sandripples often appear on beaches . Whenthe amplitudeof water oscillationis sufficientlylarge ,vortices are formed onthe lee of every sand ripple crest . A-mongthese vortices ,the most important are t…     

基于MPS方法模拟低雷诺数方柱绕流问题

柱体绕流问题是流体力学领域一个非常经典的问题。当流体流经柱体时,由于黏性的存在,会发生许多复杂的流动现象,如流动分离、涡旋周期性生成与脱落等,经常被作为标准验证算例。同时,柱体绕流广泛存在于实际工程中,并在一定工况下可能对工程产生巨大危害,因此对柱体绕流进行深入研究具有重要意义。研究中,拟将一种无网格类方法——半隐式移动粒子方法(moving particle semi-implicit method,简称MPS)引入到柱体绕流问题的数值研究中,并对不同雷诺数下二维方柱绕流问题进行数值模拟。首先,使用基于MPS方法自主开发的MLParticle-SJTU求解器,结合入口边界条件和出口边界条件,模拟了雷诺数Re分别为40、200和1 000时均匀来流条件下的方柱绕流。随后,将模拟的绕流结果与文献中试验和数值计算结果进行了对比,结果吻合较好,并且在雷诺数为200和1 000时,可以清晰地捕捉到方柱尾流中的卡门涡街现象,验证了MPS方法在柱体绕流问题模拟上的有效性和适用性。