A 3D fully coupled analysis of nonlinear sloshing and ship motion

This study investigates the coupling effects of six degrees of freedom in ship motion with fluid oscillation inside a three-dimensional rectangular container using a novel time domain simulation scheme. During the time marching, the tank-sloshing algorithm is coupled with the vessel-motion algorithm so that the influence of tank sloshing on vessel motions and vice versa can be assessed. Several factors influencing the dynamic behavior of tank–liquid system due to moving ship are also investigated. These factors include container parameters, environmental settings such as the significant wave height, current velocity as well as the direction of wind, wave and flow current acting on the ship. The nonlinear sloshing is studied using a finite element model whereas nonlinear ship motion is simulated using a hybrid marine control system. Computed roll response is compared with the existing results, showing fair agreement. Although the two hull forms and the sea states are not identical, the numerical result shows the same trend of the roll motion when the anti-rolling tanks are considered. Thus, the numerical approach presented in this paper is expected to be very useful and realistic in evaluating the coupling effects of nonlinear sloshing and 6-DOF ship motion.     

Resonant heave motion of semisubmersible vessels

This paper is concerned with nonlinear resonant heave motion of a semisubmersible vessel at the survival draft. Due to the small potential of the hulls at deep draft the resonant motion is governed almost entirely by nonlinear drag forces on the hull and bracing members.Long period ocean swell following a storm can excite semisubmersiblesat resonance in heave and delay the resumption of drilling activities.The paper presents a simple closed form solution for determining the heave motion at the natural period which agrees well with other theoretical and experimental results.     

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.     

大陆张裂三维有限元的数值模拟

采用纯数值模拟的方法,研究了大陆岩石圈在拉张和底部隆起作用下的张裂现象,编制了三维粘弹性有限元程序,结果表明:(1)单纯拉张或底部隆起时不会造成岩石圈块体层间的相互滑动,不会产生不对称共轭扩张边缘,块体的中间部分总是首先被拉薄而分离,应力应变转移集中到两侧,张裂的不对称性取决于周边的约束或加载条件,张裂共轭边缘的不对称性与周边块体的作用及其它地质改造作用密切相关。(2)拉张与底部隆起同时作用时,将会造成层间的相互滑动,且应力较大。(3)介质的塑性或流变性对张裂及应力有很大的影响,如对位移产生一定影响,同时应力值降低,使层间产生相对运动,张裂以一种渐变的形式出现。     

Maneuvering performance of high-speed ships with effect of roll motion

Equations of yaw, sway, roll and rudder motions are formulated to represent realistic maneuvering behavior of high-speed ships such as destroyers. Important coupling terms between yaw, sway, roll and rudder were included on the basis of recent captive model test results of a high-speed ship. A series of computer runs was made by using equations of yaw, sway, roll and rudder motions. Results indicate substantial coupling effects between yaw, roll, and rudder, which introduce changes in maneuvering characteristics and reduce course stability in high-speed operation. These effects together with relatively small GM (which is typical for certain high-speed ships) produce large rolling motions in a seaway as frequently observed in actual operations. Results of digital simulations and captive model tests clearly indicate the major contributing factors to such excessive rolling motions at sea.     

Interactions between nonlinear water waves and non-wall-sided 3D structures

Interactions between water waves and non-wall-sided cylinders are analyzed based on velocity potential theory with fully nonlinear boundary conditions on the free surface and the body surface. The finite element method (FEM) is adopted together with a 3D mesh generated through an extension of a 2D Delaunay grid on a horizontal plane along the depth. The linear matrix equation for the velocity potential is constructed by imposing the governing equation and boundary conditions through the Galerkin method and is solved through an iterative method. By imposing the gradient of the potential equal to the velocity, the Galerkin method is used again to obtain the velocity field in the fluid domain. Simulations are made for bottom mounted and truncated cylinders with flare in a numerical tank. Periodic waves and wave groups are generated by a piston type wave maker mounted on one end of the tank. Results are obtained for forces, wave profiles and wave runups. Further simulations are made for a cylinder with flare subjected to forced motion in otherwise still open water. Results are provided for surge and heave motion in different amplitudes, and for a body moving in a circular path in the horizontal plane. Comparisons are made in several cases with the results obtained from the second order solution in the time domain.     

A comparative linear and nonlinear ship motion study using 3-D time domain methods

This paper reports on 3D time domain seakeeping computations at different level of modeling of nonlinear effects. Four successively improved levels of computations are considered: (i) a linear computation, (ii) Froude-Krylov nonlinear computation, (iii) body nonlinear computation where the perturbation potential is computed based on the instantaneous hull under the mean free-surface, and finally (iv) the body-exact nonlinear computation, where the perturbation potential is determined based on the wetted hull under the incident wave profile after suitable mapping of the hull into a computational domain. Computations are carried out for a Wigley hull (having less ‘geometric’ nonlinearity due to vertical sides), and a S175 hull at different forward speeds. The results are obtained for both regular and irregular waves.     

On the numerical prediction of seakeeping and of structural loads of high-speed vessels

The applicability of an efficient 3D numerical method for the evaluatfion of the seakeeping performance and structural loads of advanced high-speed vessels is investigated. A semi-empirical approach to account for the viscous flow effects on ship's behavior at speed is presented and discussed. Numerical results are compared with model tests for two high-speed monohull vessels with satisfactory agreement.     

群岛峡道浓盐水排放扩散的三维数值模拟

采用三维潮流温盐数学模型对六横岛附近海域的水动力环境及盐度场分布进行数值模拟,并以实测资料进行验证,所建立模型可以较好地反映该海域潮流动力特征及盐度场的分布情况。利用验证后的三维数学模型,对海水淡化工程浓盐水排放后盐度场分布进行计算,将排放后的盐度分布与工程前进行对比分析,并推算盐升面积及垂向盐度增量。结果表明,海水淡化工程排放的浓盐水对六横岛海域盐度分布的影响主要集中在排水口附近的底层,影响区域呈带状分布,最大盐度增量为1.2左右,且排水口附近海域出现盐度垂向分层。     

Shoreline motion in nonlinear shallow water coastal models

Different shoreline boundary conditions for numerical models of the Non-Linear Shallow Water Equations based on Godunov-type schemes are compared. The study focuses on the Peregrine and Williams [Peregrine, D.H., Williams, S.M., 2001. Swash overtopping a truncated plane beach. Journal of Fluid Mechanics 440, 391–399.] problem of a single bore collapsing on a slope. This is considered the best test to assess performances of the shoreline boundary treatments in terms of all the parameters of interest in swash zone modelling. Emphasis is given to the shoreline trajectory and flow velocity modelling. A mismatch of the velocity at the early stage of the motion is highlighted. Most of the tested techniques perform similarly in terms of maximum run-up, the backwash phase is critical in all cases. Starting from the Brocchini et al. [Brocchini, M., Bernetti, R., Mancinelli, A., Albertini, G., 2001. An efficient solver for nearshore flows based on the WAF method. Coastal Engineering 43(2), 105–129.] shoreline boundary treatment, a simple technique that improves the accuracy of velocity predictions is also developed. A sensitivity analysis of the domain resolution and the threshold value of the water depth that defines a wet cell is also presented.     

An investigation on the vertical motion sickness characteristics of a high-speed catamaran ferry

Low frequent motions of vessel may cause motion sickness in rough seas. These undesirable effects induce fatigue of crews during the navigation. The motion sickness is always an important criterion for the high-speed craft design. Modern ferry designs have been marketed with a great emphasis on the seakeeping performance. This research has been carried out by investigating the results on the vertical motion sickness incidence (MSI) study for a 40 m wave-piecing catamaran at seas. The primary purpose of this research is to investigate the vertical motion sickness characteristics of a high-speed catamaran ferry. Two mathematical models, three-dimensional translating–pulsating source distribution technique and three-dimensional pulsating source distribution technique, are used for predicting the vertical acceleration responses of the wave-piecing catamaran in oblique waves. The comparison between numerical predictions and experimental data shows a good agreement except that around the pitch resonance region in FP vertical acceleration motions. Based on the experimental observation, the discrepancies may be caused by the nonlinear effects of centre bow during large pitch motions in waves. The comfort assessments are based on the ISO-2631/1997 standard with the hydrodynamic analysis for determining the acceleration levels in different locations on the vessel. The effects of seating location, wave heading and duration of motion exposure on seasickness are discussed.     

Fully nonlinear 3D interaction of bubble dynamics and a submerged or floating structure

This paper is concerned with the interaction of bubbles, a submerged or floating structure, and free surface waves. A three-dimensional fully nonlinear model has been developed based on the coupling of the boundary integral method (BIM) for bubble dynamics and free surface waves and the finite element method for structure deformation. The present method is well validated by comparing the numerical results with the experimental data. Three structure characteristics, including fixed, rigidly moving and flexible, are investigated separately to determine their influence on bubble dynamics. For a free-floating structure, the free surface causes not only a larger reduction in peak pressure for a rigid structure compared with a fixed body but also the modification of the bubble period and structural response. The interaction between a bubble and a flexible structure, in the absence of a free surface, is simulated. Both the rigid motion and the deformation at the local structure appear in the simulation. The effect of the structural thickness on the reduction in peak pressure is also considered.     

三维完全非线性波浪水槽的数值模拟

用有限元求解拉普拉斯方程,建立了三维完全非线性数值波浪水槽.跟踪流体自由表面的方法为满足完全非线性自由表面条件的半拉格朗日法,对离散单元采用20节点的六面体二次等参数单元.并把数值计算结果与水面初始升高产生箱体内流体运动解析解和二阶斯托克斯波理论解进行了对比,结果表明该模型是稳定的、守恒的,能精确模拟非线性波浪的产生和传播.     

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

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

非线性freak波分析及模拟理论探讨

概述了国内外学者对freak波的认识，并从各种观测资料分析总结了freak波的特征。引入完全非线性控制方程Schroedinger方程的1＋1模式，对freak波进行模拟，为在实验室造波或计算研究提供一种方法。     

完全非线性深水波的数值模拟

基于势流理论,并结合深水波质点运动从水面向下呈e指数衰减的特性,建立了完全非线性数值变深水槽模型,通过实时模拟活塞式造波机运动来产生波浪.采用时域高阶边界元法进行模拟,利用混合欧拉-拉格朗日方法和四阶Runge-Kutta方法追踪流体瞬时水面,应用镜像格林函数消除了水槽两个侧面的积分,在水槽末端布置人工阻尼层来消除反射...     

浅水中浮体在波浪作用下运动的三维时域计算模型

港口中系泊船在波浪作用下运动问题的本质是浅水波浪与浮体的相互作用。与深水情况不同,浅水问题应当考虑水底、水域边界的影响及浅水波浪自身的特性,单一模型很难实现该模拟过程。为此,建立了Boussinesq方程计算入射波和Laplace方程计算散射波的全时域组合计算模型。有限元法求解的Boussinesq方程能使入射波充分考虑到水底、水域边界的影响和浅水波浪的特性;散射波被线性化,采用边界元法求解,并以浮体运动时的物面条件为入射波和散射波求解的匹配条件。该方法为完全的时域方法,计算网格不随时间变动,计算过程较为方便。通过与实验及其他数值方法的结果进行比较,验证了本模型对非线性波面、浮体的运动都有比较理想的计算结果,显示了本模型对非线性问题具有较好的计算能力。     

非均匀风影响下风力机三维气动粘性流场的数值模拟

基于面向对象的开源软件OpenFOAM,选择美国国家新能源实验室(NREL)Phase VI风力机为对象,对以往研究较少的非均匀来流风速作用下风力机三维气动粘性流场进行数值模拟。采用较为接近于真实情况的指数型风剖面,计算了轮毂处风速分别为5、10、15和25 m/s四种工况下的叶片表面压力分布、叶片的推力、尾涡等气动力数据,并与均匀来流风速下的风力机气动力学性能进行详细的对比,探讨非均匀风剖面对风力机流场结构和流动特性影响的物理现象和规律。     

波浪作用下缆船拖带系统非线性运动数值模拟

基于船舶操纵性运动方程和拖缆的三维动力学运动方程,提出了被拖点位置匹配的方法,建立了拖船—拖缆—被拖船系统整体非线性拖带动力学模型。为了考察被拖船航向稳定性与横向稳性的关系以及波浪载荷作用的影响,被拖船采用水平面四自由度运动方程,并引入了波浪的作用力和力矩。拖船采用PD控制方法较真实地模拟了拖船航向改变的运动过程。对一个拖船—拖缆—被拖船系统(5 000 t的拖船和3 000 t的被拖船)在时域内进行了规则波浪作用下拖带运动的模拟,计算结果表明被拖带船舶在波浪中运动呈现运动稳定、不稳定和临界状态3种可能的特性。根据模拟计算结果,认为波浪中拖带航向稳定是被拖带船舶保持稳性的必要条件。