用MPS方法数值分析孤立波与平板相互作用问题

表面孤立波与海洋结构物的相互作用常伴随有砰击、上浪等现象,是船舶与海洋工程领域需要研究的重要问题之一。该文采用基于无网格半隐式移动粒子MPS法开发的MLParticle-SJTU求解器,对孤立波与平板结构相互作用问题进行数值计算分析研究。模拟了二维孤立波对平板的冲击作用过程,得到的平板受力与实验结果吻合良好,验证了MLParticle-SJTU求解器在处理此类问题上的可行性。通过对孤立波与平板相互作用过程的三维模拟,观察到了有别于二维空间下的波浪在平板上表面演化过程,即呈现U型演化过程,说明对此类问题进行三维数值模拟的必要性。     

THREE-DIMENSIONAL NUMERICAL SIMULATION OF WAVE SLAMMING ON AN OPEN STRUCTURE

The three-dimensional numerical simulation of wave slamming on an open structure in the splash zone is carried out based upon the Volume Of Fluid (VOF) method.A wave basin is established by solving the...     

Dynamic analysis of wave slamming on plate with elastic support

Experiments are conducted to investigate the dynamic response of a plate with elastic support under a regular wave slamming. The statistical analysis results obtained in different model testing cases are presented. The theoretical analysis of the plate vibrations(including the forced and free vibrations) is performed. Four characteristic stages of the plate vibration accelerations between two consecutive wave impacts are identified. The submergence durations of the plate during the wave action and the hydro-elastic effects are discussed. Finally, some useful conclusions are drawn.     

Numerical simulation of hydro-elastic problems with smoothed particle hydrodynamics method

Violent free surface flows with strong fluid-solid interactions can produce a tremendous pressure load on structures, resu-lting in elastic and even plastic deformations. Modeling hydro-elastic problem...     

A 3-D numerical study of solitary wave interaction with vertical cylinders using a parallelised particle-in-cell solver

This paper aims to provide a better understanding of the interaction between solitary waves and vertical circular cylinders. This is achieved via process based numerical modelling using the parallel particle-in-cell based incompressible flow solver PICIN. The numerical model solves the Navier-Stokes equations for free-surface flows and incorporates a Cartesian cut cell method for fluid-structure interaction. Solitary waves are generated using a piston-type wave paddle. The PICIN model is first validated using a test case that involves solitary wave scattering by a single vertical cylinder. Comparisons between the present results and experimental data show good agreement for the free surface elevations around the cylinder and the horizontal wave force on the cylinder. The model is then employed to investigate solitary wave interaction with a group of eleven vertical cylinders. The wave run-up and wave forces on the cylinders are discussed.     

Towards development of enhanced fully-Lagrangian mesh-free computational methods for fluid-structure interaction

Simulation of incompressible fluid flow-elastic structure interactions is targeted by using fully-Lagrangian mesh-free computational methods. A projection-based fluid model(moving particle semi-implicit(MPS)) is coupled with either a Newtonian or a Hamiltonian Lagrangian structure model(MPS or HMPS) in a mathematically-physically consistent manner. The fluid model is founded on the solution of Navier-Stokes and continuity equations. The structure models are configured either in the framework of Newtonian mechanics on the basis of conservation of linear and angular momenta, or Hamiltonian mechanics on the basis of variational principle for incompressible elastodynamics. A set of enhanced schemes are incorporated for projection-based fluid model(Enhanced MPS), thus, the developed coupled solvers for fluid structure interaction(FSI) are referred to as Enhanced MPS-MPS and Enhanced MPS-HMPS. Besides, two smoothed particle hydrodynamics(SPH)-based FSI solvers, being developed by the authors, are considered and their potential applicability and comparable performance are briefly discussed in comparison with MPS-based FSI solvers. The SPH-based FSI solvers are established through coupling of projection-based incompressible SPH(ISPH) fluid model and SPH-based Newtonian/Hamiltonian structure models, leading to Enhanced ISPH-SPH and Enhanced ISPH-HSPH. A comparative study is carried out on the performances of the FSI solvers through a set of benchmark tests, including hydrostatic water column on an elastic plate,high speed impact of an elastic aluminum beam, hydroelastic slamming of a marine panel and dam break with elastic gate.     

A NEW NUMERICAL WAVE FLUME COMBINING THE 0-1 TYPE BEM AND THE VOF METHOD

A new coupling numerical wave model, based on both the Boundary Element Method (BEM) and the Volume Of Fluid (VOF) method, is established by taking advantages of the both methods to solve the wave-structure interaction problems. In this model, the wave transformation in front of structures is calculated by the 0-1 type BEM, and the intense wave motions near the structures are calculated by the VOF method. In this paper, the characteristics of the BEM and the VOF method are discussed first, and then the coupling treatments are described in detail. In the end, the accuracy and the validity of the coupling model are examined by comparing the numerical results with experiment results and other numerical results available for the interactions between regular waves with a monolayer horizontal plate.     

基于CIP法的波物相互作用强非线性数值模拟

极端波浪环境作用下,浮体易发生上浪拍击等强非线性现象,这一问题也是波浪与浮体相互作用研究的难点之一。该文基于高阶CIP(Constrained Interpolation Profile)法,采用THINC(Tangent of Hyperbola for Interface Capturing)数值技术精确捕捉自由液面,建立了二维多相流强非线性数值水槽模型。通过模拟线性波和孤立波的生成、传播过程,对该模型的收敛性、稳定性及造波性能进行了全面的验证与分析。利用该数值模型模拟和分析了极端波浪对浮体的强非线性作用。     

基于Bragg共振原理水面多重竖直板消浪性能实验研究

该文采用实验方法探索研究了表面波经过多重穿透水面的刚性薄板结构物传播时,透射波的波高变化。将多重穿透水面的刚性薄板竖直固定及等间距布置。基于Bragg共振条件,板间距与波长的一半相当。改变板浸入水面下的深度、板的数目和波陡,进行了实验研究。实验时采用浪高仪测量模型前后波形,整理出入射波幅和透射波幅,绘出透射系数tK随无因次波浪参数kh变化的曲线。结果显示,当水波波长接近板间距两倍左右时,出现透射系数较低的现象。     

NUMERICAL SIMULATIONS OF NONLINEAR WAVES OVERTOPPING AN OBSTRUCTION

1.　ＩＮＴＲＯＤＵＣＴＩＯＮＴｈｅｐｒｏｂｌｅｍｏｆｗａｖｅｓｏｖｅｒｔｏｐｐｉｎｇｏｂｓｔｒｕｃｔｉｏｎｓｉｓｏｆｃｏｎｓｉｄｅｒａｂｌｅｉｍｐｏｒｔａｎｃｅｉｎｏｆｆｓｈｏｒｅ,ｈｙｄｒａｕｌｉｃａｎｄｃｏａｓｔａｌｅｎｇｉｎｅｅｒｉｎｇ.Ｔｈｅｓ...     

NUMERICAL SIMULATION OF SOLITARY WAVE RUN-UP AND OVER-TOPPING USING BOUSSINESQ-TYPE MODEL

In this article,the use of a high-order Boussinesq-type model and sets of laboratory experiments in a large scale flume of breaking solitary waves climbing up slopes with two inclinations are presented to study the shoreline behavior of breaking and non-breaking solitary waves on plane slopes.The scale effect on run-up height is briefly discussed.The model simulation capability is well validated against the available laboratory data and present experiments.Then,serial numerical tests are conducted to study the shoreline motion correlated with the effects of beach slope and wave nonlinearity for breaking and non-breaking waves.The empirical formula proposed by Hsiao et al.for predicting the maximum run-up height of a breaking solitary wave on plane slopes with a wide range of slope inclinations is confirmed to be cautious.Furthermore,solitary waves impacting and overtopping an impermeable sloping seawall at various water depths are investigated.Laboratory data of run-up height,shoreline motion,free surface elevation and overtopping discharge are presented.Comparisons of run-up,run-down,shoreline trajectory and wave overtopping discharge are made.A fairly good agreement is seen between numerical results and experimental data.It elucidates that the present depth-integrated model can be used as an efficient tool for predicting a wide spectrum of coastal problems.     

NUMERICAL MODELING OF WAVE EVOLUTION AND RUNUP IN SHALLOW WATER

Based on the Navier-Stokes (N-S) equations for viscous, incompressible fluid and the VOF method, 2-D and 3-D Numerical Wave Tanks (NWT) for nonlinear shallow water waves are built. The dynamic mesh technique is applied, which can save computational resources dramatically for the simulation of solitary wave propagating at a constant depth. Higher order approximation for cnoidal wave is employed to generate high quality waves. Shoaling and breaking of solitary waves over different slopes are simulated and analyzed systematically. Wave runup on structures is also investigated. The results agree very well with experimental data or analytical solutions.     

非淹没刚性植被影响下孤立波在岸滩上的爬高

为探究海岸刚性植被对海啸波的削弱效应,通过物理模型试验和数值模拟研究了孤立波与不同坡度岸滩上非淹没刚性植被的相互作用问题。物理模型试验在波浪槽中进行,测试了不同的入射波波高、植被密度和岸滩坡度对孤立波爬高的影响,并运用物理模型试验数据校核改进后的Boussinesq方程,得到植被的拖曳力系数。结果表明:拖曳力系数随植被密度的增大而增大,随坡度增大而减小;植被后的无量纲透射波高和无量纲岸滩爬高随着无量纲的入射波高的增大而减小,随着植被密度的增大而减小;当岸滩坡度增大时,无量纲透射波高增大而无量纲爬高并无显著差异。最后根据回归分析得出了岸滩爬高与相对入射波高、植被密度和岸滩坡度的幂函数型经验关系式。     

三孤立波的斜坡爬高特征

该文采用推板式造波机和改进的Goring造波方法在波浪水槽中生成了由三个等波高孤立波组成的波列,首先在坡度为1/10的斜坡上进行了多个孤立波连续爬高的实验研究。实验测量了等波高双孤立波和三孤立波在等水深波浪水槽中传播和在斜坡上的爬高过程,给出了沿程多点浪高仪测量的波面时间序列和高速摄像记录的水线在斜坡上爬高与回落过程影像。分析了来波相对波高和无量纲波峰间距对每个孤立波爬高放大系数的影响。在基于RANS方程和VOF方法的数值波浪水槽中,模拟了三孤立波爬高过程,给出了最大爬高与破碎特征。实验和数值模拟结果表明,孤立波在斜坡上回落时会产生强劲的回流水流,该反向流动导致后随孤立波破碎,消耗波动能量,使得后随孤立波的爬高放大系数明显小于没有回落水流影响的首个孤立波的爬高放大系数。当孤立的波峰间较大时,若来波相对波高较小,后两个的孤立波爬高基本相同;若来波相对波高较大时,第三个孤立波的爬高大于第二个孤立波的爬高。     

AN OPTIMIZATION OF EIGENFUNCTION EXPANSION METHOD FOR THE INTERACTION OF WATER WAVES WITH AN ELASTIC PLATE

The hydroelastic interaction of an incident wave with a semi-infinite horizontal elastic plate floating on a homogenous fluid of finite depth is analyzed using the eigenfunction expansion method. The fluid is assumed to be inviscid and incompressible and the wave amplitudes are assumed to be small. A two-dimensional problem is formulated within the framework of linear potential theory. The fluid domain is divided into two regions, namely an open water region and a plate-covered region. In this paper, the orthogonality property of eigenfunctions in the open water region is used to obtain the set of simultaneous equations for the expansion coefficients of the velocity potentials and the edge conditions are included as a part of the equation system. The results indicate that the thickness and the density of plate have almost no influence on the reflection and transmission coefficients. Numerical analysis shows that the method proposed here is effective and has higher convergence than the previous results.     

NUMERICAL STUDY OF MULTIDIRECTIONAL FOCUSING WAVE RUN-UP ON A VERTICAL SURFACE-PIERCING CYLINDER

In this article,a numerical model based on improved Boussinesq equations and the Finite Element Method(FEM)with unstructured triangular elements is proposed and verified by experimental results for the focusing wave group interaction with a vertical surface-piercing cylinder.The multidirectional focusing waves with different directional spreading parameters and their interaction with a vertical surface-piercing cylinder are numerically simulated.The directional spreading parameter influences the wave form,that is,a wider directional spreading would induce a narrower wave form in space.The wave run-up on the cylinder increases with the increase of the wave directional spreading parameter,but when the parameter is greater than 30-40,the wave run-up changes only slightly and the influence of the directional spreading can be neglected.A larger diameter cylinder may have a larger wave run-up at the front face.But at the back side,it is other way round.     

部分可渗透直立圆柱的消波作用研究

本文导出了线性波浪作用于上部可渗透的直立刚性圆柱体上的水动力的解析解。在线性有势流的假定下，用特征函数展开方法求出了内部和外部流场的解。给出了不同波浪要素及结构尺寸参数下的数值结果。由结果可知，提高可渗透性可以有效消减波幅。     

3-D simulation of transient flow patterns in a corridor-shaped air-cushion surge chamber based on computational fluid dynamics

The 3-D characteristics of the water-air flow patterns in a corridor-shaped air-cushion surge chamber during hydraulic transients need to be considered in the shape optimization. To verify the reliability of the water-air two-phase model, namely, the volume of fluid model, the process of charging water into a closed air chamber is successfully simulated. Using the model, the 3-D flow characteristics under the load rejection and acceptance conditions within the air-cushion surge chamber of a specific hydropower station are studied. The free surface waves, the flow patterns, and the pressure changes during the surge wave process are analyzed in detail. The longitudinal flow of water in the long corridor-shaped surge chamber is similar to the open channel flow with respect to the wave propagation, reflection and superposition characteristics. The lumped parameters of the 3-D numerical simulation agree with the results of a 1-D calculation of hydraulic transients in the whole water conveying system, which validates the 3-D method. The 3-D flow structures obtained can be applied to the shape optimization of the chamber.     

NUMERICAL SIMULATION OF PROPAGATION AND BREAKING PROCESSES OF A FOCUSED WAVES GROUP

A two-dimensional numerical wave flume is developed to study the focused waves group propagation and the consequent breaking processes. The numerical model is based on the Reynolds-Averaged Navier-Stokes (RANS) equations, with the standard k-ε turbulence model to simulate the turbulence effects. To track the complicated and broken free-surface, the Volume Of Fluid (VOF) method is employed. The numerical model combines the "Partial Cell Treatment (PCT)" method with the "Locally Relative Stationary (LRS)" concept to treat the moving wave paddle so that various waves can be generated directly in a fixed Cartesian grid system. The theoretical results of the linear and nonlinear waves are used to validate the numerical wave flume firstly, and then a plunging breaking wave created by a focused waves group is simulated. The numerical results are compared to the experimental data and other simulation results, with very good agreements. The turbulence intensity, the flow field and the energy dissipation in the breaking processes are analyzed based on the numerical results. It is shown that the present numerical model is efficient and accurate for studying the waves group generation, the waves packet propagation, and the wave breaking processes.     

Application of CLEAR-VOF method to wave and flow simulations

A two-dimensional numerical model based on the Navier-Stokes equations and computational Lagrangian-Eulerian advection remap-volume of fluid (CLEAR-VOF) method was developed to simulate wave and flow problems. The Navier-Stokes equations were discretized with a three-step finite element method that has a third-order accuracy. In the CLEAR-VOF method, the VOF function F was calculated in the Lagrangian manner and allowed the complicated free surface to be accurately captured. The propagation of regular waves and solitary waves over a flat bottom, and shoaling and breaking of solitary waves on two different slopes were simulated with this model, and the numerical results agreed with experimental data and theoretical solutions. A benchmark test of dam-collapse flow was also simulated with an unstructured mesh, and the capability of the present model for wave and flow simulations with unstructured meshes, was verified. The results show that the model is effective for numerical simulation of wave and flow problems with both structured and unstructured meshes.