考虑定常兴波影响的辐射问题去奇Rankine源方法数值计算

随着中高速船舶在海洋资源开发中的广泛应用,航速对船舶在波浪中摇荡运动的影响受到重视。本文基于势流理论,采用非均匀有理B样条(NURBS)表达物面,应用去奇Rankine源方法求解流场速度势;推导了兴波问题的二阶近似非线性自由面边界条件,以及考虑二阶近似非线性兴波影响的辐射速度势的线性自由面和物面边界条件。本文编程首先数值计算了Wigley船型的近似非线性兴波问题,将计算的线性和非线性兴波阻力系数与文献结果做了比较;然后数值计算了计及线性定常兴波影响的下潜圆球有航速辐射问题,将计算的附加质量和阻尼系数与文献结果比较。本文计算结果与文献结果吻合程度较好,验证了本文数值计算方法的有效性。文中还给出了兴波波形和辐射波形。     

船体在波浪中的非线性水动压力

本文提出了一个预报船体在波浪中大幅运动时非线性水动压力场的二维时域理论。船体扰动势用时域自由面格林函数和在入射波下的瞬时湿表面上的分布源求解；与非线性水动压力场相匹配的船体运动用差分法求得。为提高计算效率和避免数值过程发散，采用了改进的数值模型和方案。通过线性理论计算与模型试验结果的比较，指出了线性切片理论在预报水动压力场时的不足，水动压力与波高的非线性关系及正负水动压力沿船体表面的分布在Wigley船的计算比较中得到了说明。初步计算表明，该理论的实用化发展前景是令有鼓舞的。相应的计算机程序可在PC机上运行。     

港湾内非线性波浪对船剖面作用的时域计算模型(英文)

文中用边界元数值求解复速度势模拟了波浪对船剖面的时域作用力,采用满足运动学和动力学条件的半拉朗格朗日法跟踪流体的自由表面.推导了适用本数值模型的辐射边界条件,此边界条件和海绵层边界条件组成的出口边界条件,可以使波浪很好地被吸收或传播出去.进行了相关物理模型实验,数值计算结果和实验结果的对比表明此模型能精确模拟波浪对船的作用力.     

数值波浪水池中船舶水动力系数测试与分析技术

以CFD理论为基础,建立了数值波浪水池,给出了一种基于三维数值波浪水池的船舶水动力系数测试与分析技术。对Wigley-III船模在数值波浪水池中受迫振荡进行了数值模拟,计算分析了船体垂荡、纵摇及纵荡运动的附加质量与阻尼,并与三维势流理论计算结果进行了比较,两者吻合良好。此方法能准确给出浮式结构物的水动力系数,细致描述船舶周围的流场,可广泛应用于船舶与海洋工程浮式结构物的水动力性能研究。     

浅窄航道非线性兴波阻力数值计算

应用Rankine面元法计算船舶在浅窄航道中作匀速直航运动时的非线性兴波阻力．基于势流理论，在船体表面、自由面以及岸壁分布Rankine源来表达流场速度势，利用迭代方法满足非线性自由面边界条件，采用映像法满足水底边界条件，采用网格错位法满足辐射条件．以Wigley数学船型为例进行计算，通过改变航速和航道宽度参数，得到相应的兴波阻力．通过分析不同航速下航道宽度变化对兴波阻力的影响，得到了一些有意义的结论．     

数值波浪水池及顶浪中船舶水动力计算

基于粘流理论建立了三维数值波浪水池,模拟了非线性波浪,并对规则波顶浪中前进的拘束船模的水动力进行了计算.数值模拟中,控制方程-RANS方程和连续性方程使用有限体积法离散,非线性自由面采用VOF方法处理;在入口边界模拟柔性造波板运动产生入射波,使用位于波浪水池尾部的人工阻尼区消波.给出了非线性规则波的模拟结果以及规则波顶浪中前进的拘束船模的水动力计算结果,并与理论解及DUT(Delfi University of Technology)的试验数据进行了比较,二者吻合良好.     

规则迎浪中船舶参数横摇的三维时域预报方法研究

参数横摇是船舶因复原特性改变而引起的典型非线性现象。文章采用三维时域方法预报规则迎浪中船舶的参数横摇运动。该方法引入匹配面将流域分为内域和外域,内域中采用Rankine源来满足物面条件和线性自由面条件,而外域中应用时域格林函数来满足线性自由面条件和远场辐射条件。数值方法中,Froude-Krylov力和恢复力是通过对船舶瞬时湿表面积积分获得,同时考虑了横摇、垂荡和纵摇三自由度之间的耦合作用,以及非线性横摇阻尼的影响。数值结果与试验结果吻合很好,说明该方法可以有效地预报参数横摇。     

Time-domain TEBEM method for mean drift force and moment of ships with forward speed under the oblique seas

A numerical method for solving 3D unsteady potential flow problem of ship advancing in waves is put forward. The flow field is divided into an inner and an outer domain by introducing an artificial matching surface. The inner domain is surrounded by ship wetted surface and matching surface as well as part of the free surface. The free surface condition for the inner domain is formulated by perturbation about the double-body flow or uniform incoming flow assumption. The outer domain is surrounded by matching surface and the rest free surface as well as infinite far-field radiation boundary. The free surface condition for the outer domain is formulated by perturbation about uniform incoming flow. The simple Green function and transient free surface Green function are used to form the boundary integral equation (BIE) for the inner and outer domains, respectively. Taylor Expansion Boundary Element Method (TEBEM) is utilized to solve the double-body flow and inner domain and outer domain unsteady flow BIE. Matching conditions for the inner domain flow and outer domain flow are enforced by the continuity of velocity potential and normal velocity on the matching surface. Direct pressure integration on ship wetted surface is used to obtain the first-order and second-order wave forces (moments). The numerical predictions on the displacement, added resistance, sway mean drift force and yaw mean drift moment of the modified KVLCC2 ship at different forward speeds are investigated by the proposed TEBEM method. It is also compared with the other numerical results. The physical tank experiment results are also developed to validate the accuracy of numerical tank results. Compared with the experiment solutions, a good agreement can be obtained by TEBEM method.

     

采用THOBEM方法研究波流与结构物相互作用的数值模拟(英文)

To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method (HOBEM). By assuming small flow velocities, the velocity potential could be expressed for linear and higher order components by perturbation expansion. A 4th-order Runge-Kutta method was applied for time marching. An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves. Validation of the numerical method was carried out on run-up, wave exciting forces, and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder. The results were in close agreement with the results of a frequency-domain method and a published time-domain method. The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.     

Time-domain simulation for water wave radiation by floating structures （Part A）

Direct time-domain simulation of floating structures has advantages： it can calculate wave pressure fields and forces directly; and it is useful for coupled analysis of floating structures with a mooring system. A time-domain boundary integral equation method is presented to simulate three-dimensional water wave radiation problems. A stable form of the integration free-surface boundary condition （IFBC） is used to update velocity potentials on the free surface. A multi-transmitting formula （MTF） method with an artificial speed is introduced to the artificial radiation boundary （ARB）. The method was applied to simulate a semi-spherical liquefied natural gas （LNG） carrier and a semi-submersible undergoing specified harmonic motion. Numerical parameters such as the form of the ARB, and the time and space discretization related to this method are discussed. It was found that a good agreement can be obtained when artificial speed is between 0.6 and 1.6 times the phase velocity of water waves in the MTF method. A simulation can be done for a long period of time by this method without problems of instability, and the method is also accurate and computationally efficient.     

基于CFD的船舶不同浪向下运动及增阻数值分析

为分析波浪对船舶快速性和耐波性的影响,必须对波浪中航行的船舶阻力增值进行准确预报。本文基于计算流体力学软件FINE/Marine建立了Wigley船模的数值模型,对不同规则波波长下的船体运动和波浪增阻进行了计算,并与试验结果进行对比,验证了数值模型的可行性与准确性。同时计算分析了船舶在规则波中航行时的波浪增阻与浪向之间的变化关系。通过研究发现：随着浪向角的增大船舶波浪增阻逐渐增加,在60°浪向角时波浪增阻达到最大值,浪向角对波浪增阻的影响较大。本文的研究方法可用于船舶有航速下的不同浪向波浪增阻的数值预报。     

Analysis of hydrodynamic characteristics for arbitrary multihull ships advancing in waves

Multihull vessels have emerged as popular alternatives to conventional monohull ships for high-speed crafts. However, the bridging structures connecting the hulls are vulnerable to various wave actions and the wave impact on the bottom of them is the most serious problems associated with multihulled vessels. In this study, prediction of relative wave elevations under the bridging structures is investigated for multihull ships traveling with forward speed in waves. A computer code YNU-SEA using the three-dimensional (3D) Green function method with forward speed has been developed and used to analyze the hydrodynamic radiation and diffraction forces and motion responses for high-speed catamarans in waves. The results of the present calculations are compared with those of previous calculations as well as with experimental results. The numerical results reveal that the present computer code can be used as a powerful tool for the accurate numerical computation of seakeeping problems for multihull ships advancing in waves. Numerical calculations of wave pattern are also carried out including wave interactions between the hulls to analyze the effects of hull form on the free surface flow around catamarans advancing in waves. The analysis of the wave pattern allows the determination of relative wave height including radiation and diffraction waves. Finally, some discussions are included based on these numerical results which may be helpful for the accurate prediction of relative wave height and wave breaking load on the deck associated with multihull ships.     

Three-dimensional numerical wave tank simulations on fully nonlinear wave–current–body interactions

A finite-difference scheme and a marker-and-cell (MAC) method are used for numerical wave tank (NWT) simulations to investigate the characteristics of nonlinear wave motions and their interactions with a stationary three-dimensional body in the presence of steady uniform currents. The Navier–Stokes (NS) equation is solved in the computational domain, and the boundary values are updated at each time-step by a finite-difference time-marching scheme in the frame of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is satisfied by the marker–density function technique developed for two fluid layers. The incident waves are generated from the inflow boundary by prescribing a velocity profile resembling the motions of a flexible flap wavemaker, and the outgoing waves are numerically dissipated inside an artificial damping zone located at the end of the tank. Using the NS–MAC NWT, nonlinear wave and current interactions around a stationary vertical truncated circular cylinder are studied, and the results are compared with the experimental results of Mercier and Niedzwecki, a time-domain NWT based on linear potential theory, a fully nonlinear NWT, and a second-order diffraction computation. Received: July 3, 2001 / Accepted: September 25, 2001     

波浪中母船干扰下潜水员梯的运动及水动力特性分析

潜水员梯是通过销轴连接在母船舷侧为潜水员执行潜水作业时上下船的平台,为保证在波浪中使用时的安全性,对其波浪作用力、运动响应及销轴处的动作用载荷设计有较高要求。为研究母船干扰下潜水员梯的动态响应特性,采用基于深水格林函数为积分内核的边界元方法,先计算母船及潜水梯这一整体的运动响应及其平均湿表面的辐射和绕射压力分布,依据刚体运动合成原理求解得到潜水梯重心处的运动。再以潜水员梯为研究对象,采用压力积分方法获得潜水员梯的波浪力作用,并将其运动、波浪力和销轴处的作用力联立构成动力学平衡方程,逆向求解出销轴处的动作用力。基于上述思路,系统研究了波长、航速及潜水员梯与母船相对位置变化时潜水员梯的运动、波浪作用力及销轴处作用力的变化规律,为潜水员梯的布置、销轴的选型及潜水作业时母船的操纵可提供指导意义。     

Linear seakeeping and added resistance analysis by means of body-fixed coordinate system

This paper presents an alternative formulation of the boundary value problem for linear seakeeping and added resistance analysis based on a body-fixed coordinate system. The formulation does not involve higher-order derivatives of the steady velocity potential on the right-hand side of the body-boundary condition, i.e., the so-called m _j -terms in the traditional formulation when an inertial coordinate system is applied. Numerical studies are made for a modified Wigley I hull, a Series 60 ship with block coefficient 0.7, and the S175 container ship for moderate forward speeds where it is thought appropriate to use the double-body flow as basis flow. The presented results for the forced heave and pitch oscillations, motion responses, and added resistance in head-sea waves show good agreement with experiments and some other numerical studies. A Neumann–Kelvin formulation is shown to give less satisfactory results, in particular for coupled heave and pitch added mass and damping coefficients.     

利用直接边界元法求解兴波问题

本文描述了利用直接边界元法求解定常兴波问题的理论和数值方法，本方法利用线性化自由面条件，辐射条件基于Dawson的四点差分格式进行处理。本文给出了Wigley船的兴波阻力和波高的数值结果，并与Dawson方法和试验结果进行了比较，结果显示本文的结果合理可靠。     

On nonlinear ship waves and wave resistance calculations

A finite difference method based on the Euler equations is developed for computing ship waves and wave resistances. Time marching is carried out using a time-splitting fractional-step method. The second-order central difference is used to discretize the spatial differentials, while the convection terms are discretized by the QUICK scheme. A body- and free-surface-fitted grid system with a cell-centered stencil is used. A Poisson equation for the pressure increment at each time step is solved to enforce mass conservation. The method is validated by comparing the numerical results with the experimental data for a Wigley parabolic hull. The characteristics of ship waves, such as the wave profile along the hull, the wave pattern on the free surface, the pressure distribution on the hull surface, and the wave-making resistance are reasonably predicted. The calculated results are in good agreement with the experimental data. Received for publication on Oct. 13, 1998; accepted on Jan. 20, 1999     

基于MFS无网格方法的船舶耐波性预报

文章提出了一种新的基于MFS法(Method of Fundamental Solutions,即基本解法)的求解船舶势流问题的无网格数值方法。该方法采用势流控制方程的基本解表达流场中的速度势,可以得到任意阶连续的速度势导数;另外,通过引入流场虚边界的概念避免了传统边界元方法中存在奇点的问题。相比传统的边界元方法,MFS法具有无奇性、不依赖网格、不涉及数值积分、数学原理简单和算法易于编程实现的优势。文中对于MFS法在求解浅水二维剖面水动力系数中的应用进行了探索,提出了一套实用的数值算法,并进一步结合频域切片法计算了Wigley船体的垂荡和纵摇响应,通过和频域面元法及试验结果的对比验证了方法的正确性和实用性;最后应用开发的程序分析了不同水深和航速对于耐波性的影响。该研究开创性地将MFS方法应用到了船舶水动力研究领域,为之后更深入广泛地应用这一特性优良的数值方法提供了有益的参考。     

基于N-S方程的航行船舶辐射问题数值模拟

基于N-S(Navier-Stokes)方程,进行了航行船舶辐射问题的数值模拟,包括强制垂荡和强制纵摇的模拟.计算了船舶垂荡和纵摇的附加质量和阻尼.数值模拟中,控制方程--BANS(Reynolds Averaged Navier-Stokes)方程和连续性方程使用有限体积法离散,非线性自由面采用VOF方法处理.文中给出了以不同航速前进的船舶强制垂荡和纵摇的力与力矩,以及船舶垂荡和纵摇的附加质量系数和阻尼系数,并与DUT(Delft University of Technology)的试验数据进行了比较,二者吻合良好.