Strip method for a laterally drifting ship in waves

Lateral drift occurs due to the effects of wind forces, wave drifting forces, or both on ships sailing in actual seas. It is important therefore to investigate the influence of lateral drift on seakeeping performance for improved ship operation. The velocity potential was expanded as an asymptotic power series in terms of the lateral speed parameter, τ, defined as ω _e V ₀/g, where ω _e is the frequency of wave encounter; V ₀ denotes the lateral velocity, which is assumed to be sufficiently small; and g is the acceleration due to gravity. By combining this technique with the strip method, two sets of motion equations of all the hydrodynamic force coefficients for ship seakeeping were derived. The first set is for ships without lateral drift and is the same as the equations in the new strip method, and the second set is for the additional motions induced by lateral drift. It was found that all ship motion modes except surge are coupled when a ship drifts laterally in waves.     

采用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.     

二阶波浪力数值计算与试验方法研究

目前准确测量波浪中船舶与海洋结构物二阶波浪力是一个难点,本文提出采用半约束法和弹簧系统约束法两种试验方法测量二阶波浪力.首先介绍了半约束法和弹簧系统约束法两种试验方法,然后针对某散货船采用这两种方法开展了二阶平均波浪力试验测量研究,最后将试验结果与理论计算结果进行了对比,并对两种试验方法的优缺点进行了分析.研究结果表明:两种试验方法均能准确测量二阶平均波浪力.     

Broaching prediction in the light of an enhanced mathematical model,with higher-order terms taken into account

 We have attempted to develop a more consistent mathematical model for capsizing associated with surf-riding in following and quartering waves by taking most of the second-order terms of the waves into account. The wave effects on the hull maneuvring coefficients were estimated, together with the hydrodynamic lift due to wave fluid velocity, and the change in added mass due to relative wave elevations. The wave effects on the hydrodynamic derivatives with respect to rudder angles were estimated by using the Mathematical Modelling Group (MMG) model. Then captive ship model experiments were conducted, and these showed reasonably good agreements between the experiments and the calculations for the wave effects on the hull and the rudder maneuvring forces. It was also found that the wave effects on restoring moments are much smaller than the Froude–Krylov prediction, and the minimum restoring arm appears on a wave downslope but not on a wave crest amidship. Thus, an experimental formula of the lift force due to the heel angle of the ship is provided for numerical modelling. Numerical simulations were then carried out with these second-order terms of waves, and the results were compared with the results of free-running model experiments. An improved prediction accuracy for ship motions in following and quartering seas was demonstrated. Although the boundaries of the ship motion modes were also obtained with both the original model and the present one, the second-order terms for waves are not so crucial for predicting the capsizing boundaries themselves. Received: June 20, 2002 / Accepted: October 10, 2002 Acknowledgments. This research was supported by a Grant-in-Aid for Scientific Research of the Ministry of Education, Culture, Sports, Science and Technology of Japan (No. 13555270). The authors thank Prof. N. Rakhmanin of the Krylov Ship Research Institute for providing the Russian literature, as well as Mr. H. Murata of NHK (Japan Broadcasting Corporation) for translating it into Japanese. Address correspondence to: N. Umeda (e-mail: umeda@naoe.eng.osaka-u.ac.jp)     

Large-amplitude motion responses of a Ro–Ro ship to regular oblique waves in intact and damaged conditions

This article presents a nonlinear time-domain simulation method for the prediction of large-amplitude motions of a Ro–Ro ship in regular oblique waves in an intact and a damaged condition. Numerical computations and model tests have been carried out to investigate the dynamic motion responses of Ro–Ro ship Dextra to various wave amplitudes at three different wave headings. The results of numerical and experimental investigations for stern quartering waves are reviewed. Comparisons between predictions and measurements show good agreement except in the roll-resonant region. Nonlinear effects are significant in horizontal modes of motion, and resonant roll motion, and there is strong coupling between all modes of motion in the roll-resonant region for large-amplitude responses. On the other hand, the time-domain simulation technique suffers from numerical drift in horizontal modes of motion as wave amplitude increases. This is due to nonlinear equations of motion and the lack of a restoring force and moment in horizontal motion. Received: April 30, 2002 / Accepted: August 9, 2002 Acknowledgments. II Programme of the European Community Commission under contract No. BRPR-CT97-0513. Address correspondence to: H.S. Chan (hoi-sang.chan@ncl.ac.uk)     

并靠两浮体上的波浪漂移力(英文)

An innovative hydrodynamic theory and numerical model were developed to help improve the efficiency, accuracy, and convergence of the numerical prediction of wave drift forces on two side-by-side deepwater floating bodies. The wave drift forces were expressed by the double integration of source strength and the corresponding Green function on the body surface, which is consistent with the far field formula based on momentum conservation and sharing the advantage of near field calculations providing the drift force on each body. Numerical results were validated through comparing the general far field model and pressure integral model, as well as the middle field model developed using the software HydroStar.     

FPSO内孤立波载荷实验观测及理论预报模型研究

利用大型密度分层水槽开展了下凹型内孤立波作用在FPSO上的载荷特性系列实验;并依据实验工况,考虑KdV、eKdV和MCC内孤立波理论的适用性条件,数值研究了FPSO内孤立波载荷成分构成;基于实验结果和载荷成分构成,建立了FPSO内孤立波载荷的理论预报模型.研究表明:FPSO内孤立波水平载荷由粘性力和Froude-Krylov力组成,而垂向载荷主要为垂向Froude-Krylov力;Froude-Krylov力可通过动压力沿FPSO浮体湿表面积分得到,粘性力则通过经实验回归的摩擦系数Cf、形状修正因数K乘以内孤立波诱导水质点切向速度沿FPSO浮体湿表面积分得到.系列实验结果得出:摩擦系数Cf和形状修正因数K与雷诺数Re、KC数和流体层深度比h1/h有关;摩擦系数与Re呈自然对数关系;而形状修正因数K与KC数呈幂函数关系.理论预报模型预报的水平载荷、垂向载荷结果均与系列实验和数值结果吻合较好,并且发现随着内孤立波振幅的增加,载荷幅值近乎线性增加,而且上层流体深度对水平力幅值有明显的影响.     

多功能支持船波浪载荷研究

波浪载荷研究是船舶性能预报和结构安全评估的关键内容。文中利用势流理论直接数值求解多功能支持船遭受的波浪载荷(包括一阶波频力和二阶漂移力),利用统计方法得到波浪载荷的长期预报和短期预报极值,包括波浪总纵弯矩和垂向波浪剪力以用于后续的船舶整体结构安全评估,以及纵向漂移力、横向漂移力和首摇漂移力矩用于船舶的动力定位能力分析中。研究结果表明:求解波浪载荷的势流理论方法精确,能为后续的船舶性能预报和结构安全评估提供可靠的载荷输入。     

Numerical simulation of interference effects for a high-speed catamaran

The simulations of the flow around a high-speed vessel in both catamaran and monohull configurations are carried out by the numerical solution of the Reynold averaged Navier–Stokes (RANS) equations. The goal of the analysis is the investigation of the interference phenomena between the two hulls, with focus on its dependence on the Reynolds number (Re). To this aim, numerical simulations are carried out for values of Re ranging from 10⁶ to 10⁸ for two different values of the Froude number (Fr = 0.30, 0.45). Wave patterns, wave profiles, limiting streamlines, surface pressure and velocity fields are analyzed; comparison is made between the catamaran and the monohull configurations. Dependence of the pressure and viscous resistance coefficients, as well as of the interference factor, on the Reynolds number is investigated. Verification and validation for both resistance coefficients and wave cuts is also performed.     

An approach to estimation of near-surface turbulence and CO2 transfer velocity from remote sensing data

The air–sea CO₂ exchange is primarily determined by the boundary-layer processes in the near-surface layer of the ocean since it is a water-side limited gas. As a consequence, the interfacial component of the CO₂ transfer velocity can be linked to parameters of turbulence in the near-surface layer of the ocean. The development of remote sensing techniques provides a possibility to quantify the dissipation of the turbulent kinetic energy in the near-surface layer of the ocean and the air–sea CO₂ transfer velocity on a global scale. In this work, the dissipation rate of the turbulent kinetic energy in the near-surface layer of the ocean and its patchiness has been linked to the air–sea CO₂ transfer velocity with a boundary-layer type model. Field observations of upper ocean turbulence, laboratory studies, and the direct CO₂ flux measurements are used to validate the model. The model is then forced with the TOPEX POSEIDON wind speed and significant wave height to demonstrate its applicability for estimating the distribution of the near-surface turbulence dissipation rate and gas transfer velocity for an extended (decadal) time period. A future version of this remote sensing algorithm will incorporate directional wind/wave data being available from QUIKSCAT, a now-cast wave model, and satellite heat fluxes. The inclusion of microwave imagery from the Special Sensor Microwave Imager (SSM/I) and the Synthetic Aperture Radar (SAR) will provide additional information on the fractional whitecap coverage and sea surface turbulence patchiness.     

Coupled renewal model of ocean viscous sublayer, thermal skin effect and interfacial gas transfer velocity

A previously developed renewal model included parameterizations for the thermal skin effect and interfacial gas transfer velocity. The more readily available cool skin data were used for an adjustment of the gas transfer parameterization. In this work, the renewal concept is extended to include the velocity difference across the viscous sublayer and to account for the stage of surface wave development. As a result, the empirical coefficients that enter the renewal model have been specified more accurately using laboratory data on the surface wind drift current. In addition, the coefficient linking the cool skin and gas transfer parameterization formulas has been determined from the probability distribution function for renewal events. A comparison of the upgraded renewal model with the thermal skin data collected during the COARE and more recent field programs and with gas transfer data collected during GasEx-01 experiment suggests that the renewal model can be a useful tool for producing a physically based parameterization for the interfacial CO₂ transfer velocity. Model uncertainties associated with the effect of surface films are discussed.     

Higher-order boundary element method for the interaction of a floating body with both waves and slow current

A hybrid boundary element method is suggested to solve the problem of the interaction of floating structures with both waves and slow current. A pulsating source and its mirror image referring to the sea bottom are adopted as the Green's function. The velocity potentials are expanded into an eigenfunction expansion in the outer region of the fluid domain while higher order elements are used to discretize the boundary surface surrounding the inner region. The method is validated by comparing calculated results for a circular cylinder with the semi-analytical solutions. The method is then applied to ellipsoids of various breadth and draft to investigate the influence of body shape on the wave drift damping.     

The hydrodynamic forces acting on a cylinder array oscillating in waves and current

The problem of the interaction of multiple cylinders oscillating in waves and slow current is considered. The interaction is represented by waves emitted from adjacent cylinders towards the cylinder under consideration. Wave drift forces and moment in the horizontal plane are calculated by the far-field method based on the conservation of momentum or angular momentum. A semianalytical formula for the calculation of the wave drift damping is then deduced. The conservation of the integrals in these formulae is proved. Special treatments to improve the accuracy of results are discussed. Comparisons between calculated results and experimental measurements are made, showing satisfactory agreement. Effects of various combinations of current direction and incident wave angle on the wave drift damping and damping moment are also examined.     

Hydrodynamic characteristics of slotted breakwaters under regular waves

The performance of a slotted breakwater consisting of one row of vertical slots was investigated theoretically and experimentally under normal regular waves. A simple theoretical model based on an eigenfunction was developed. The wave transmission, reflection, energy loss, and hydrodynamic force exerted on the breakwater were calculated for different values of the wave and structure parameters. The validity of the theoretical model was examined by comparing its results with theoretical and experimental results obtained from different studies. It was found that the transmission coefficient decreases with increasing dimensionless wavenumber (kh), increasing wave steepness (H _i/L), and decreasing breakwater porosity (ε). The reflection coefficient showed the opposite trend to the transmission coefficient. Also, about 20–50% of incident wave energy was lost due to the effect of the breakwater. In addition, the proposed theoretical model can be used for predicting the performance of slotted breakwaters and the hydrodynamic forces exerted on these structures using the friction coefficient f = 1.5.     

波浪对潮流影响的数值模拟及其应用

文章将水动力因素分为潮流特征周期与波浪特征周期两部分,建立了考虑波浪影响的潮流方程。通过考虑波浪引起的辐射应力、表面和底部边界层、Stokes漂流以及Coriolis波浪应力等因素,利用数值方法分析了波浪对潮流场的影响,并将其应用到渤海区域中,分析渤海内波浪对于潮流的影响。     

An improved method for estimating internal solitary waves force on small circular cylindrical structures

Theoretically calculation and laboratory experiments were performed to investigate the horizontal force exerted by internal solitary wave (ISW) on small circular cylindrical structures. Improvement in calculation of drag force and inertia force by introducing two-drag-coefficient method and spatial derivative of velocity into the calculation of force by eigenfunction was made based on systematically study on theoretical foundations. Particle image velocimetry (PIV) method and high precision force sensor were used to measure the velocity and force exerted by ISW, respectively. ISW waveform was traced by analyzing image frames follow grayscale threshold method. For the first time, components of horizontal force were calculated from PIV data, which show good agreement with theoretical one. Moreover, series comparisons between experimental data and calculational results were carried out to examine the proposed method. The results show ideal accuracy of the method and significant improvement comparing to the traditional method.     

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.

     

两船在波浪载荷作用下的运动响应分析

基于波浪辐射/绕射理论,利用水动力分析软件AQWA研究两船在波浪作用下的运动响应特性,给出驳船受浮式生产储油船(FPSO)影响时的附加质量和辐射阻尼,分析了驳船的漂移力和运动响应。计算结果表明,当两船间距为5 m时驳船的附加质量比间距为10 m时的值大2倍多,驳船的附加质量和阻尼随着间距的减小逐渐增大;在间距30 m时,驳船的速度越小所受的漂移力就越大,并且驳船的横摇和首摇运动就越剧烈。     

基于势流理论的船舶操纵水动力分析

以船舶操纵水动力预报为研究背景,以无界流场势流理论为基础,采用边界元法,对目标Wigly船在斜航运动状态下的流场速度和转首力矩进行计算。分析了斜航时船舶周围流场的分布规律,研究了转首力矩随船舶漂角的变化规律。结果表明,该研究对预报斜航时船舶的运动规律具有较好的应用前景。