Freely floating-body simulation by a 2D fully nonlinear numerical wave tank

Nonlinear interactions between large waves and freely floating bodies are investigated by a 2D fully nonlinear numerical wave tank (NWT). The fully nonlinear 2D NWT is developed based on the potential theory, MEL/material-node time-marching approach, and boundary element method (BEM). A robust and stable 4th-order Runge–Kutta fully updated time-integration scheme is used with regriding (every time step) and smoothing (every five steps). A special φ_n-η type numerical beach on the free surface is developed to minimize wave reflection from end-wall and wave maker. The acceleration-potential formulation and direct mode-decomposition method are used for calculating the time derivative of velocity potential. The indirect mode-decomposition method is also independently developed for cross-checking. The present fully nonlinear simulations for a 2D freely floating barge are compared with the corresponding linear results, Nojiri and Murayama’s (Trans. West-Jpn. Soc. Nav. Archit. 51 (1975)) experimental results, and Tanizawa and Minami’s (Abstract for the 6th Symposium on Nonlinear and Free-surface Flow, 1998) fully nonlinear simulation results. It is shown that the fully nonlinear results converge to the corresponding linear results as incident wave heights decrease. A noticeable discrepancy between linear and fully nonlinear simulations is observed near the resonance area, where the second and third harmonic sway forces are even bigger than the first harmonic component causing highly nonlinear features in sway time series. The surprisingly large second harmonic heave forces in short waves are also successfully reproduced. The fully updated time-marching scheme is found to be much more robust than the frozen-coefficient method in fully nonlinear simulations with floating bodies. To compare the role of free-surface and body-surface nonlinearities, the body-nonlinear-only case with linearized free-surface condition was separately developed and simulated.     

Numerical wave tank computations of nonlinear motions of two-dimensional arbitrarily shaped free floating bodies

Time domain simulations of nonlinear motions of two-dimensional floating bodies in waves are presented. The so called `body exact' approach is adopted in a numerical wave tank. A new scheme for pressure evaluation on the wetted hull is developed and systematically used with good results in terms of accuracy and stability. Strongly flared geometries are successfully handled even at very large amplitude motions. The validation of the code is carried out according to the 20th ITTC Seakeeping Committee recommendations through internal checks of consistency and through comparisons with available experimental data.     

水下平板与波浪相互作用完全非线性数值模拟(英文)

利用完全非线性数值波浪水槽技术研究水下平板与波浪的相互作用。假定水下平板厚度极薄、刚性,位于有限水深并且非常接近自由水面。应用四阶龙格库塔方法追踪每一时刻的波面形状,采用阻尼层来吸收反射波以保证算法的稳定性,同时引入平滑和重组的方法抑制自由表面控制点的较高梯度。通过对波浪与浮动圆柱相互作用的数值模拟证实了数值波浪水槽方法的有效性,计算结果与线性理论吻合良好。在波浪数值水槽方法中引入造波板模拟波浪产生并与水下平板发生相互作用,应用傅立叶解析方法对波面变形、波浪力作了分析。结果表明在板非常接近自由水面的情况下会表现出现很强的非线性,揭示了线性理论的局限性。     

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.     

Fully nonlinear numerical wave tank (NWT) simulations and wave run-up prediction around 3-D structures

A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed.     

Bragg reflection in a fully nonlinear numerical wave tank based on boundary integral equation method

We investigated the phenomenon of Bragg reflection of submerged structures in a 2D fully nonlinear numerical wave tank (NWT) based on the boundary integral equation method (BIEM). This model was validated by comparing not only the free surface elevations with that of the analytic solution of Stokes’ second-order wave theory, but also the reflection coefficients of submerged bars with that from other sources. The results of the present model show that the free surface nonlinear effect on the reflection coefficient of the primary resonance reduces significantly for all of the submerged bars considered. Finally, a case study is presented to demonstrate the reflecting capacity and overall performance of various submerged bars. Results indicate that sinusoidal bar has the maximum reflection capacity at the primary resonance, but the trapezoidal submerged bar is suggested as the better option for the practical convenience of coastal underwater construction.     

完全非线性波浪与二维固定结构物作用的IGN-BEM耦合分析模型

为高效准确地对完全非线性波浪与二维固定结构物的相互作用进行模拟分析,建立了二维完全非线性时域耦合模型。耦合模型将计算域划分为靠近结构物的内域和远离结构物的外域,每个区域均采用满足完全非线性自由水面边界条件的波浪模型进行求解。在内域使用Laplace方程描述流体运动并采用高阶边界元法（BEM）对其进行求解;而在没有结构物的外域,波浪运动的控制方程为Irrotational Green-Naghdi（IGN）方程并采用有限元法（FEM）对其进行求解。内域和外域通过一段重叠区域进行耦合,从而实现模型间变量的传递。首先利用耦合模型分别对规则波的传播、直墙前立波的生成以及相关物理模型试验进行模拟,数值结果与精确解和试验结果的良好吻合验证了耦合模型耦合方式的合理性以及处理非线性问题的准确性;然后使用耦合模型模拟分析了波浪与固定结构物间的相互作用,并将结果与线性解析解以及完全非线性BEM模型的结果进行了对比分析,进一步证明了耦合模型的正确性与高效性。     

Direct time domain analysis of floating structures with linear and nonlinear mooring stiffness in a 3D numerical wave tank

In this paper, motion response of a moored floating structure interacting with a large amplitude and steep incident wave field is studied using a coupled time domain solution scheme. Solution of the hydrodynamic boundary value problem is achieved using a three-dimensional numerical wave tank (3D NWT) approach based upon a form of Mixed-Eulerian–Lagrangian (MEL) scheme. In the developed method, nonlinearity arising due to incident wave as well as nonlinear hydrostatics is completely captured while the hydrodynamic interactions of radiation and diffraction are determined at every time step based on certain simplifying approximations. Mooring lines are modelled as linear as well as nonlinear springs. The horizontal tension for each individual mooring line is obtained from the nonlinear load-excursion plot of the lines computed using catenary theory, from which the linear and nonlinear line stiffness are determined. Motions of three realistic floating structures with different mooring systems are analyzed considering various combinations of linear and approximate nonlinear hydrodynamic load computations and linear/nonlinear mooring line stiffness. Results are discussed to bring out the influence and need for consideration of nonlinearities in the hydrodynamics and hydrostatics as well as the nonlinear modelling of the line stiffness.     

Hydroelastic analysis of fully nonlinear water waves with floating elastic plate via multiple knot B-splines

Hydroelastic analysis of fully nonlinear water waves with the floating elastic plate is a hard mission. Especially, the behavior of the wave would be more complex when water wave encounter the floating elastic plate. In this paper, the meshless numerical method is devoted to solve such a problem. Fundamental solution method is applied to approximate the velocity potential in the fluid domain. When the water wave encounters the plate, the wave function would not be enough smooth in the edge of plate compared to the other points. Hence, to analyze numerically the behavior of wave, the solution space should include the basis functions that are not enough smooth in the edge of plate. Moreover, to decrease computational cost significantly, the basis functions had better to have local compact support. The multiple knot B-spline basis functions are suitable that contain both properties. The number of repeated knots, the degree of B-spline and the spatial points are challengeable that are discussed in this paper. The results are in good agreement with those obtained from other numerical works.