KCS迎浪中耐波性的CFD数值模拟

采用计算流体动力学(CFD)软件STAR-CCM,对国际集装箱标模KCS在迎浪规则波中的运动响应和波浪增阻等主要耐波性能进行数值计算。采用动态重叠网格处理船体运动,同时利用VOF方法模拟自由面波形演化情况。计算结果显示:动态重叠网格具有较强的适应性,能较好地处理船体的大幅运动,计算得到的船体运动幅值能与试验结果相一致,静水和波浪中的阻力计算结果与试验结果能较好地吻合;从自由面波形方面看,KCS船体兴波与来波存在较强的相互干扰,当波峰撞击船首时,波面起伏较大,数值计算能较好地模拟这些流动现象。     

分段式尾压浪板对高速船阻力性能的影响

基于CFD技术,以排水型高速船Model 5b为模型,寻求改善高速船阻力性能的尾压浪板新形式。首先尾压浪板新形式的确定在静水条件下进行,然后在波浪条件下验证该压浪板的阻力性能。基于CFD软件建立三维数值波浪水池,静水条件下采用切割体网格技术预报船模的阻力性能。波浪条件下数值水池入口采用直接造波方法,尾部采用人工阻尼消波方法,自由面采用VOF方法处理,采用重叠网格技术预报船模的阻力性能以及运动响应。确定一种比常规压浪板阻力性能优良的分段式压浪板,为船舶节能附体的研究提供参考。     

Validation of 6-DOF motion simulations for ship turning in regular waves

In this study, a six degrees of freedom (6-DOF) motion simulation method of a ship steering in regular waves is validated. The proposed simulation model is based on the two-time scale concept where the 6-DOF motions are expressed as the sum of the low-frequency maneuvering motions and high-frequency wave-induced motions. Turning simulations of a KCS container ship model with a rudder angle of \(\pm 35^\circ\) in calm water and regular waves are performed and the obtained results are compared with the results of a free-running model test. The model tests were conducted using a ship model of length 3.057 m in a square tank at the National Research Institute of Fisheries Engineering, Japan. The wave conditions were as follows: the wave height was 3.6 m at full-scale, ratio of wavelength to ship length was 1.0, and the ship approached in the head wave direction before it was steered. The present method can simulate both the turning motion and wave-induced motions in regular waves with practical accuracy.

     

URANS analysis of a broaching event in irregular quartering seas

Ship motions in a high sea state can have adverse effects on controllability, cause loss of stability, and ultimately compromise the survivability of the ship. In a broaching event, the ship losses control, naturally turning broadside to the waves, causing a dangerous situation and possibly capsizing. Classical approaches to study broaching rely on costly experimental programs and/or time-domain potential or system-based simulation codes. In this paper the ability of Reynolds averaged Navier–Stokes (RANS) to simulate a broaching event in irregular waves is demonstrated, and the extensive information available is used to analyze the broaching process. The demonstration nature of this paper is stressed, as opposed to a validated study. Unsteady RANS (URANS) provides a model based on first principles to capture phenomena such as coupling between sway, yaw, and roll, roll damping, effects of complex waves on righting arm, rudders partially out of the water, etc. The computational fluid dynamics (CFD) method uses a single-phase level-set approach to model the free surface, and dynamic overset grids to resolve large-amplitude motions. Before evaluating irregular seas two regular wave cases are demonstrated, one causing broaching and one causing stable surf riding. A sea state 8 is imposed following an irregular Bretschneider spectrum, and an autopilot was implemented to control heading and speed with two different gains for the heading controller. It is concluded that the autopilot causes the ship to be in an adverse dynamic condition at the beginning of the broaching process, and thus is partially responsible for the occurrence of the broaching event.     

Optimization of bow shape for a non ballast water ship

In this research,a commercial CFD code ＂Fluent＂ was applied to optimization of bulbous bow shape for a non ballast water ships（NBS）.The ship was developed at the Laboratory of the authors in Osaka Prefecture University,Japan.At first,accuracy of the CFD code was validated by comparing the CFD results with experimental results at towing tank of Osaka Prefecture University.In the optimizing process,the resistances acting on ships in calm water and in regular head waves were defined as the object function.Following features of bulbous bow shapes were considered as design parameters： volume of bulbous bow,height of its volume center,angle of bow bottom,and length of bulbous bow.When referring to the computed results given by the CFD like resistance,pressure and wave pattern made by ships in calm water and in waves,an optimal bow shape for ships was discovered by comparing the results in the series of bow shapes.In the computation on waves,the ship is in fully captured condition because shorter waves,λ/Lpp 0.6,are assumed.     

A unified seakeeping and maneuvering analysis of ships in regular waves

The behavior of a ship in regular waves during maneuvering was studied by using a two-time scale model. The maneuvering analysis was based on Söding’s (Schiffstechnik 1982; 29:3–29) nonlinear slender-body theory generalized to account for heel. Forces and moments due to rudder, propeller, and viscous cross-flow follow from the state-of-the-art procedures. The developed unified theory of seakeeping and maneuvering was verified and validated for calm water by comparing it with experimental and calculated zigzag and circle maneuvers. Linear wave-induced motions and loads were determined by generalizing the Salvesen-Tuck-Faltinsen (Trans SNAME 1970; 78:250–287) strip theory. The mean second-order wave loads in incident regular deep water waves in oblique sea conditions were estimated by the potential flow theories of Faltinsen et al. (Proc 13th Symp Naval Hydrody 1980), Salvesen (Proc Intl Symp Dynam Mar Vehicl Struct Wave 1974), and Loukakis and Sclavounos (J Ship Res 1978; 22:1–19). The considered theories cover the whole range of important wavelengths. Comparisons between the different mean second-order wave load theories and available experimental data were carried out for different ship hull forms when the ship was advancing forward on a straight course. The mentioned methods have been incorporated into the maneuvering model. Their applicability from the perspective of the maneuvering ability of the selected types of ships was investigated in given wave environments. The wave conditions are valid for realistic maneuvering cases in open coastal areas. It was demonstrated that the incident waves may have an important influence on the maneuvering behavior of a ship. The added resistance, mean second-order transverse force, and yaw moment also play important roles.     

波浪中螺旋桨非定常水动力性能研究（英文）

The speed of a ship sailing in waves always slows down due to the decrease in efficiency of the propeller. So it is necessary and essential to analyze the unsteady hydrodynamic performance of propeller in waves. This paper is based on the numerical simulation and experimental research of hydrodynamics performance when the propeller is under wave conditions. Open-water propeller performance in calm water is calculated by commercial codes and the results are compared to experimental values to evaluate the accuracy of the numerical simulation method. The first-order Volume of Fluid(VOF) wave method in STAR CCM+ is utilized to simulate the three-dimensional numerical wave. According to the above prerequisite, the numerical calculation of hydrodynamic performance of the propeller under wave conditions is conducted, and the results reveal that both thrust and torque of the propeller under wave conditions reveal intense unsteady behavior. With the periodic variation of waves, ventilation, and even an effluent phenomenon appears on the propeller. Calculation results indicate, when ventilation or effluent appears, the numerical calculation model can capture the dynamic characteristics of the propeller accurately, thus providing a significant theory foundation forfurther studying the hydrodynamic performance of a propeller in waves.     

迎浪规则波中波浪增阻和船体垂向运动的数值预报(英文)

The numerical prediction of added resistance and vertical ship motions of one ITTC (International Towing Tank Conference) S-175 containership in regular head waves by our own in-house unsteady RANS solver naoe-FOAM-SL JTU is presented in this paper. The development of the solver naoe-FOAM-SJTU is based on the open source CFD tool, OpenFOAM. Numerical analysis is focused on the added resistance and vertical ship motions (heave and pitch motions) with four very different wavelengths (0.8Lpp≤λ≤1.5Lpp) in regular head waves. Once the wavelength is near the length of the ship model, the responses of the resistance and ship motions become strongly influenced by nonlinear factors, as a result difficulties within simulations occur. In the paper, a comparison of the experimental results and the nonlinear strip theory was reviewed and based on the findings, the RANS simulations by the solver naoe-FOAM-SJTU were considered competent with the prediction of added resistance and vertical ship motions in a wide range of wave lengths.     

A time-domain higher-order boundary element method for 3D forward-speed radiation and diffraction problems

A time-domain higher-order boundary element method for seakeeping analyses in the framework of linear potential theory is newly developed. Ship waves generated by two modified Wigley models advancing at a constant forward speed in calm water or incident waves and the resultant radiation and diffraction forces are computed to validate this code. A rectangular computational domain moving with the same forward speed as the ship is introduced, in which an artificial damping beach is installed at an outer portion of the free surface except the downstream side for satisfying the radiation condition. The velocity potential on the ship hull and the normal velocity on the free surface are calculated directly by solving the boundary integral equation. An explicit time-marching scheme is employed for updating both kinematic and dynamic free-surface boundary conditions, with an embedding of a second-order upwind difference scheme for the derivative in the x-direction to stabilize the calculation. Extensive results including the exciting forces, added mass and damping coefficients, wave profiles, and wave patterns for blunt Wigley and slender Wigley hulls with forward speed are presented to validate the efficiency of the proposed 3D time-domain approach. The corresponding physical tests of the radiation and diffraction problems in a towing tank are also carried out. Computed numerical results show good agreement with the corresponding experimental data and other numerical solutions.     

船舶在波浪中操纵运动仿真

本文研究了双浆双舵船在规则波中的回转运动，首先进行了约束模型试验，得到了操纵运动数学模型中的水动力系数，然后，进行了静水操纵运动数值仿真，并与自航模型试验结果进行了比较。最后，预报船舶在规则波中的回转运动，对一些影响回转运动的因素进行了讨论。     

计算流体力学在船舶线型优化中的应用

以某船舶的计算流体力学(简称CFD)阻力性能计算为例,结合传统的实验流体力学(简称EFD)模型试验计算结果,比较分析了该船原型和改型的性能差异,总结和归纳了CFD在船型优化方面的优势。     

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

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

Evaluation of added resistance in regular incident waves by computational fluid dynamics motion simulation using an overlapping grid system

A computational fluid dynamics simulation method called WISDAM-X was developed to evaluate the added resistance of ships in waves. The Reynolds-averaged Navier–Stokes (RANS) equation was solved by the finite-volume method and a MAC-type solution algorithm. An overlapping grid system was employed to implement rigorous wave generation, the interactions of ships with incident waves, and the resultant ship motions. The motion of the ship is simultaneously solved by combining the solution of the motion of the ship with the solution of the flow about the ship. The free surface is captured by treatment by the density-function method. The accuracy of WISDAM-X is examined by a comparison with experimental data from a container carrier hull form, and shows a fairly good agreement with respect to ship motion and added resistance. Simulations were also conducted for a bow-form series of a medium-speed tanker to examine the effectiveness of the WISDAM-X method as a design tool for a hull form with a smaller resistance in waves. It was confirmed that the WISDAM-X method can evaluate the added resistance with sufficient relative accuracy and can be used as a design tool for ships.     

Added resistance acting on hull of a non ballast water ship

In this paper, added resistances acting on a hull of non ballast water ship(NBS) in high waves is discussed. The non ballast water ships were developed at the laboratory of the authors at Osaka Prefecture University, Japan. In the present paper, the performances of three kinds of bow shapes developed for the NBS were theoretically and experimentally investigated to find the best one in high waves. In previous papers, an optimum bow shape for the NBS was developed in calm water and in moderated waves. For a 2 m model for experiments and computations, the wave height is 0.02 m. This means that the wave height is 15% of the draft of the ship in full load conditions. In this paper, added resistances in high waves up to 0.07 m for a 2 m model or 53% of the full load draft are investigated. In such high waves linear wave theories which have been used in the design stage of a ship for a long time may not work well anymore, and experiments are the only effective tool to predict the added resistance in high waves. With the computations for waves, the ship is in a fully captured condition because shorter waves, λ/Lpp0.6, are assumed.     

基于数值模拟的船艏抗气泡性能分析

对船舶静水中航行问题进行了CFD数值模拟,分析了船艏抗气泡性能。为获得船体周围气泡分布规律,在船前来流中添加气泡,并基于拉格朗日描述法追踪气泡的运动。根据计算得到的气泡分布特点,通过修改船艏线型改善了船舶的抗气泡性能。此外,开展了模型试验,通过船体表面的流线验证了数值模拟方法的可靠性。     

弹性船体波激响应时域过程的数值模拟

本文从入射波的波浪现实出发，用非线性水弹性积分─微分法进行船体运动和结构动响应的时间历程模拟，直接将试验测量结果与理论计算在时域内作比较，这有利于揭示砰击和上浪等非线性瞬态作用导致冲荡响应的力学机理。这种模拟及其验证方法避免了对理论和试验结果的统计特性进行比较时，因子样长度的限制和数据处理过程所带来的不确定性，比峰值比较更能说明波激响应数学模型的精确性和可靠性。     

基于RANS法和边界层理论预报三维船体阻力

基于船体阻力成因和分类,应用二因次RANS法并结合EFD技术为研究手段,研究了船体各阻力成分和兴波波形,考虑边界层第一层网格高度对数值计算结果的影响,通过数值模拟计算和兴波波形比较,分析表明第一层网格无量纲高度y+值充分影响数值结果精度,总阻力在高傅汝德数下其数值误差有增大趋势,同一航速下摩擦阻力值随y+增加略微增大,剩余阻力值却明显减少,且自由面兴波随傅汝德数增大首部波高愈发明显,船中及尾部区域兴波向外扩散,船首尾肩部开尔文波系明显且波峰后移,其趋势与EFD波形一致,具有实用性。     

高速三体船波浪中运动与增阻CFD计算研究

基于数值波浪水池技术,对波浪中高速三体船运动及增阻进行CFD计算研究。控制方程—RANS方程和连续性方程使用有限体积法离散,非线性自由面采用VOF方法处理;在入口边界模拟柔性造波板运动产生入射波,使用位于波浪水池尾部的人工阻尼区消波。首先对规则波顶浪中单个主船体的运动和增阻进行了计算,并与高速细长体理论计算结果进行了比较;随后进行了三体船运动和增阻计算,分析了侧片体对主船体阻力增加的影响。为高速三体船耐波性研究提供了数值工具。     

三维波对船舶运动响应影响研究

通常情况,实际海况的波浪更接近三维波,而在船舶运动预报时,目前主要考虑二维波(长峰波)中的响应,不能真正反映船舶在实际海浪中的运动.文中利用不同方法(频域方法、时域方法以及模型)研究了三维波的方向扩散性对船舶运动响应的影响.三种方法得到的三维波对船舶运动的影响是一致的.结果显示,在研究实际海况中的船舶运动时,考虑三维波的方向扩散性是非常重要的.     

Ship motion-sloshing interaction using a field method

Even in relatively calm waters, low amplitude wave-induced motions of an LNG carrier may induce large amplitude liquid sloshing inside the ship's partially filled tanks, and the interaction between ship motions and sloshing may affect the ship's seekeeping properties. A computational procedure, here referred to as the RANS-RANS method, was developed to account for this interaction, and this method was then employed to predict the free surface flow inside the tanks and the corresponding motions of the ship in regular head and beam waves. This method coupled a compressible VoF technique with a generic wave generation and absorption scheme to obtain wave-induced ship motions with and without considering the effects of sloshing in the ship's tanks. Systematic grid studies were performed to obtain a sufficiently fine grid needed to yield converging predictions. The resulting wave patterns, ship motions, and internal sloshing elevations were compared with results obtained from a computational method, here referred to as the RANS-BEM method, that relied on a boundary element method to obtain ship motions. This RANS-RANS method was validated against model test measurements.