Dynamic loads on submerged bodies in a viscous numerical wave tank at small KC numbers

In this paper, an improved version of the MAC method (SIMAC), recently developed at the University of Trieste, is employed for the study of the wave generation and propagation into a numerical wave tank and for the evaluation of dynamic loads over submerged fixed bodies.In the first part of the work, a numerical wave tank was developed. A pneumatic wave-maker at the left-hand side of the tank was implemented by the use of a pressure perturbation at the free surface. The pressure varies in time with a sinusoidal law. Grid sensitivity tests, checks on mass conservation and the Fourier analysis of the waves which propagate in the tank showed the effectiveness of SIMAC in treating such problems. The wave-maker was then calibrated.In the second part of the work, the dynamic loads over submerged square and rectangular cylinders were evaluated. The time records of the horizontal and vertical forces which act over the body were then treated using the Morison equation in order to derive the inertial and damping coefficients. The analysis was carried out for KC numbers ranging between 0.447 and 3.58. Numerical results satisfactorily tallied with experimental data. The analysis of the velocity field near the body evidenced the influence of vortex generation and vortex shedding on the coefficients of inertial forces.     

Wave propagation in a fully nonlinear numerical wave tank: A desingularized method

The problem of wave propagation in a fully nonlinear numerical wave tank is studied using desingularized boundary integral equation method coupled with mixed Eulerian–Lagrangian formulation. The present method is employed to solve the potential flow boundary value problem at each time step. The fourth-order predictor–corrector Adams–Bashforth–Moulton scheme is used for the time-stepping integration of the free surface boundary conditions. A damping layer near the end-wall of wave tank is added to absorb the outgoing waves with as little wave reflection back into the wave tank as possible. The saw-tooth instability is overcome via a five-point Chebyshev smoothing scheme. The model is applied to several wave propagations including solitary, irregular and random incident waves.     

Linear irregular wave generation in a numerical wave tank

In the design of any floating or fixed marine structure, it is vital to test models in order to understand the fluid/structure interaction involved. A relatively inexpensive method, compared to physical model testing, of achieving this is to numerically model the structure and the wave conditions in a numerical wave tank. In this paper, a methodology for accurately replicating measured ocean waves in a numerical model at full scale is detailed. A Fourier analysis of the measured record allows the wave to be defined as a summation of linear waves and, therefore, Airy's linear wave theory may be used to input the wave elevation and associated water particle velocities. Furthermore, a structure is introduced into the model to display the ability of the model to accurately predict wave–structure interaction. A case study of three individual measured waves, which are recorded at the Atlantic marine energy test site, off the west coast of Ireland, is also presented. The accuracy of the model to replicate the measured waves and perform wave–structure interaction is found to be very high. Additionally, the absolute water particle velocity profile below the wave from the numerical model is compared to a filtered analytical approximation of the measured wave at a number of time-steps and is in very good agreement.     

一个二维数值波浪水槽的改进

对Drimer,Agnon及Segre研发的数值波浪水槽(简称DAS),首先通过修正DAS中某些积分计算的错误,建立了第一个修正模型(MDAS1),使得模型更加准确稳定。其次,将DAS模型自由表面线性元近似替换为更为合理准确的三阶元近似,建立了第二个修正模型(MDAS2),使得模型计算的波面更合理可靠。平底水槽中波群的生成及传播算例表明:与DAS结果相比,MDAS1对波面没有显著影响,而采用三阶元的MDAS2对波面有显著影响。与Hansen和Svendsen的实测资料相比较证实了两个修正模型的有效性。     

基于物理造波和主动消波的黏性数值波浪水槽开发

基于推板造波理论和摇板造波理论,在Open FOAM平台上采用重叠网格技术建立黏性数值波浪水槽,并使用一种结合SIMPLE算法和PISO算法的PIMPLE算法对数值模型进行求解。利用开发的数值模型通过数值收敛性测试和网格独立性测试分别重点研究了时间步长、库朗数和网格尺寸对数值精度和计算效率的影响。并对比研究了此数值模型分别嵌入层流模型和湍流模型的计算精度和计算效率。实现的规则波和二阶有限振幅波与理论结果和试验结果吻合,验证了此黏性数值波浪水槽的造波和主动消波功能。基于二维数值波浪水槽,进一步研究了三维数值造波,数值计算结果与理论结果吻合良好。研究结果不仅验证了重叠网格在二维和三维两相流体域中求解运动物体与流场交互的可靠性和正确性,而且为使用此黏性数值波浪水槽解决更复杂的海洋工程问题提供了依据。     

Determining the onset and strength of unforced wave breaking in a numerical wave tank

A numerical wave tank is used to investigate the onset and strength of unforced wave breaking, and the waves have three types of initial spectra： constant amplitude spectrum, constant steepness spectrum and Pierson-Moscowitz spectrum. Numerical tests are performed to validate the model results. Then, the onset of wave breaking is discussed with geometric, kinematic, and dynamic breaking criteria. The strength of wave breaking, which is always characterized by the fractional energy loss and breaking strength coefficient, is studied for different spectra. The results show how the energy growth rate is better than the initial wave steepness on estimating the fractional energy losses as well as breaking strength coefficient.     

不规则波在孔隙介质礁坪上传播过程的波高和增水变化数值研究

采用$\sigma $坐标系统下以体积平均的雷诺时间平均方程作为控制方程的三维非静压模型,对随机波浪在带有孔隙介质的岛礁地形上的传播过程进行了模拟,重点分析了礁坪上方波高和增水的变化。通过与多个组次工况的物理模型实验数据进行对比,结果显示,本文模型能很好地模拟波浪在孔隙介质上传播演化的过程,与实验结果吻合程度很高。分析结果表明,相比于光滑底床,孔隙介质的存在造成破碎点附近波高平均下降12%,礁坪上方波高平均下降28%。对于平均水位,孔隙底床条件下的最大减水幅值减小了43%,同时礁坪上方增水幅值上升6%。另外,孔隙率在0.47～0.87范围内变化时,对礁坪上方平均水位的变化基本无影响。     

数值水槽内的阻尼消波和波浪变形计算

采用时域内对波面运动位置追踪的边界元方法,建立了一种非线性波浪变形计算的三维数值模式,并借助Spongelayer阻尼消波和Sommerfeld放射条件相匹配的处理方式,开发了一条三维非线性波的数值造波水槽,进而对水槽内的Stokes波进行了波浪变形计算。     

Linear and nonlinear irregular waves and forces in a numerical wave tank

A time-domain higher-order boundary element scheme was utilized to simulate the linear and nonlinear irregular waves and diffractions due to a structure. Upon the second-order irregular waves with four Airy wave components being fed through the inflow boundary, the fully nonlinear boundary problem was solved in a time-marching scheme. The open boundary was modeled by combining an absorbing beach and the stretching technique. The proposed numerical scheme was verified by simulating the linear regular and irregular waves. The scheme was further applied to compute the linear and nonlinear irregular wave diffraction forces acting on a vertical truncated circular cylinder. The nonlinear results were also verified by checking the accuracy of the nonlinear simulation.     

Experimental and numerical investigation of non-predictive phase-control strategies for a point-absorbing wave energy converter

Phase control may substantially increase the power absorption in point-absorber wave energy converters. This study deals with validation of dynamic models and latching control algorithms for an oscillating water column (OWC) inside a fixed vertical tube of small circular cross-section by small-scale testing. The paper describes experimental and numerical results for the system's dynamics, using simple and practical latching control techniques that do not require the prediction of waves or wave forces, and which will be relevant to any type of point-absorbing devices.In the experimental set-up, the upper end of the tube was equipped with an outlet duct and a shut-off valve, which could be controlled to give a latching of the inner free surface movement. The pressure drop through the open valve is used as a simplified measure of the energy extraction. The control was realized by using the real-time measurement signals for the inner and outer surface displacement.A mathematical model of the system was established and applied in numerical simulation. In the case the OWC's diameter is much smaller than the wavelength and the wave amplitude much smaller than the draft, the free surface movement inside the tube can be described as an oscillating weightless piston. For this hydrodynamic problem an analytical solution is known. In addition, the mathematical model includes the effects of viscous flow losses, the air compressibility inside the chamber and the pressure drop across the valve. Experimental results were used to calibrate some of the model parameters, and the total model was formulated as a coupled system of six non-linear, first-order differential equations. Time-domain integration was used to simulate the system in order to test the control strategies and compare with experimental results.     

基于MPS无网格方法的孤立波传播的数值模拟

在MPS无网格方法中,引进预定候选粒子集概念用以生成邻接粒子集矩阵,使该部分的机时耗费缩短为引进前的1/11;采用Bi-CGSTAB方法求解压力泊松方程,显著地提高了求解速度.模拟了孤立波在数值波浪水槽中的传播及其与直墙作用时的爬升、回落过程,结果表明模拟波面结果与解析值及实测结果基本相符,针对不同波高的孤立波计算得到的墙前最大爬升值与实测结果也是一致的.     

Numerical simulation of solitary wave propagation based on MPS method

The concept of candidate particle set is introduced in the MPS gridless numerical method to generate neighboring particle set matrix, which can reduce the CPU time to 1/11 of that before introduction. The Bi-CGSTAB (bi-conjugate gradient stabilized) algorithm is applied to solving the Poisson pressure equation, by which the solving speed is significantly accelerated. The process of solitary waves propagating over a numerical flume and interacting with a vertical wall is simulated. The simulated results of water surface elevation are in good agreement with the analytical solution as well as the measured data. The predicted maximum values of the run-up of solitary waves with various relative incident wave heights agree well with the measured results.     

Effects of mesoscale eddies on the internal solitary wave propagation

The mesoscale eddy and internal wave both are phenomena commonly observed in oceans. It is aimed to investigate how the presence of a mesoscale eddy in the ocean affects wave form deformation of the internal solitary wave propagation. An ocean eddy is produced by a quasi-geostrophic model in f-plane, and the one-dimensional nonlinear variable-coefficient extended Korteweg-de Vries (eKdV) equation is used to simulate an internal solitary wave passing through the mesoscale eddy field. The results suggest that the mode structures of the linear internal wave are modified due to the presence of the mesoscale eddy field. A cyclonic eddy and an anticyclonic eddy have different influences on the background environment of the internal solitary wave propagation. The existence of a mesoscale eddy field has almost no prominent impact on the propagation of a smallamplitude internal solitary wave only based on the first mode vertical structure, but the mesoscale eddy background field exerts a considerable influence on the solitary wave propagation if considering high-mode vertical structures. Furthermore, whether an internal solitary wave first passes through anticyclonic eddy or cyclonic eddy, the deformation of wave profiles is different. Many observations of solitary internal waves in the real oceans suggest the formation of the waves. Apart from topography effect, it is shown that the mesoscale eddy background field is also a considerable factor which influences the internal solitary wave propagation and deformation.     

孤立波爬坡的二维数值模拟

波浪在滩地上以及遇海岸工程后传播发生变形、爬坡等现象,对其进行数值模拟具有广阔的工程应用背景.应用基于Boltzmann方程的KFVS(kinetic flux vector splitting)格式求解二维浅水方程,同时采用干底Riemann解模拟动边界问题.模型模拟了孤立波在滩地上传播变形、爬坡的过程,以及孤立波在滩地上遇圆柱后绕射、变形和爬坡的过程.计算结果与实验结果非常吻合,表明模型具有较大的推广应用价值.     

基于高阶边界元的三维数值波浪港池

初步建立了一个基于高阶边界元的三维数值波浪港池,港池具有造波和消波功能。采用高阶边界元16节点四边形单元和基于二阶显式泰勒展开的混合欧拉-拉格朗日时间步进求解带自由表面的完全非线性势流方程。模型中对于影响数值精度的问题作了细致的处理。数值计算结果表明本港池可以用来模拟非线性波浪的传播,具有很高的数值精度和稳定性。     

Comparative study on the accuracy of solitary wave generations in an ISPH-based numerical wave flume

In the present study, six solitary wave generations by different mathematical approximations are investigated using a piston type wave maker at dimensionless amplitudes ranging from 0.1 to 0.6 and two water depths. Incompressible smoothed particle hydrodynamics is used to simulate solitary wave propagation along the fixed depth channel. The present numerical results are compared with analytical results and experimental data in terms of free surface displacements, fluid particle velocity, phase speed, paddle motion, etc. The present mesh-free numerical results of wave profile variations over time proved that “Rayleigh” has the lowest relative wave height variation. However, its solitary wave has notable phase lead, while “Third order” and “Ninth order” have the least wave lags. Furthermore, the record of present numerical free surface elevation at different distances and the loss of amplitude of the main pulse showed that regarding both of them, “Ninth order” has supremacy over five others. Considering the numerical velocity components of generated solitary wave, “Third order” and “Ninth order” trace analytical results more accurately than other four ones, whereas “Rayleigh” is the most accurate one in predicting the maximum runup. Finally, the paddle motion, its velocity, and displacement, as well as phase speed and outskirts decay coefficient are also compared and discussed intensely.     

基于二阶时间步进法数值波浪水槽研究

基于势流理论和时域高阶边界元方法,建立了三维完全非线性数值波浪水槽模型.利用源造波法产生入射波浪,应用五阶斯托克斯波理论给定波浪速度;采用混合欧拉-拉格朗日方法追踪流体瞬时水面,将二阶泰勒级数展开法应用于更新下一时间步的波面和速度势;通过加速势的方法准确计算自由水面速度的法向导数和物面速度势的时间导数.对完全非线性波浪进行了模拟,得到了稳定的波形.当波浪非线性较小时,与四阶Runge- Kutta法(RK4)计算结果和五阶斯托克斯波理论解均吻合良好;随着波浪非线性的增大,计算结果误差逐渐增大.通过数值试验分析,在满足精度要求的基础上,本方法计算时间略大于四阶Runge- Kutta法的四分之一,大大减少了计算量.     

基于HUST-Ship的船舶耐波性数值试验水池参数化建模方法研究

基于自研的HUST-Ship黏性流CFD求解器,提出了一种船舶耐波性数值试验水池参数化建模方法和网格生成技术。在自研的CFD软件平台中,根据波浪参数(规则波的波长、波高和不规则波的特征周期、有效波高)自动生成船舶耐波性数值试验水池。以结构化重叠网格技术为基础,分别对船舶与水池进行网格划分并进行网格的组装、重叠,实现船舶在耐波性数值试验水池中的数值“试验”,对船舶耐波性进行数值预报。根据波浪种类分为规则波数值试验水池与不规则波数值试验水池,可对各类波浪条件实现精确可靠的数值造波,模拟船舶在各种波浪环境中运动响应、耐波性事件与非线性运动特性。采用参数化建模方法生成的耐波性数值试验水池能够提高建模和网格生成效率,通过对国际标模进行数值预报并与船模试验结果对比,验证了参数化建模方法生成的耐波性数值试验水池对船舶在波浪中运动性能的计算精度。     

珊瑚岛礁孤立波爬坡的平面二维数值模拟研究

为了探究岛屿周围珊瑚礁在抵御海啸灾害中的作用,采用激波捕捉类Boussinesq模型FUNWAVE-TVD,对孤立波在理想化三维岛礁地形上的传播及爬坡开展了现场尺度的平面二维数值模拟,分析了入射波高、礁坪水深、礁坪宽度、礁前斜坡坡度、礁后斜坡坡度、珊瑚礁糙率对岛屿四周孤立波爬高分布的影响。结果表明,珊瑚礁的存在总体上可有效降低岛屿四周孤立波的最大爬坡高度;入射波高、礁坪水深、礁坪宽度、珊瑚礁糙率是影响珊瑚岛礁四周孤立波爬坡分布的主要因素,岛礁四周最大爬坡高度会随入射波高和礁坪水深的增大、礁坪宽度和珊瑚礁糙率的减小而不断增大;当礁坪水深增大到一定程度时,珊瑚礁主要会对岛屿背浪面的爬高失去影响,而当礁坪宽度和珊瑚礁糙率减小至一定程度时,会出现岛礁四周最大爬高高于无珊瑚礁时爬高的现象;礁后斜坡的变缓会使岛礁周围的最大爬高有所减小,而礁前斜坡坡度对珊瑚岛礁周围的最大爬高几乎没有影响。     

孤立波在岸礁地形上传播与爬坡数值模拟研究

为准确模拟孤立波在岸礁地形上的传播和爬坡,采用基于完全非线性Boussinesq方程开发的Funwave-TVD模型,探究模型的可行性,并利用验证后的模型进一步研究岸礁各地形因素对孤立波爬高的影响。研究结果表明:模型能准确模拟孤立波在岸礁陡变地形上的传播及变形,摩擦系数对礁前陡坡及礁坪上的波浪传播模拟影响不大,但对爬坡预测的敏感性较强;模型空间步长可适当增大,提高计算效率;随着礁坪宽度的增大以及礁后斜坡的变缓,孤立波爬坡高度下降明显,而礁前陡坡坡度变化对孤立波爬坡高度影响不大。