On the numerical accuracy of the prediction of resistance coefficients in ship stern flow calculations

This article presents a study on the accuracy of the numerical determination of the friction and pressure resistance coefficients of ship hulls. The investigation was carried out for the KVLCC2 tanker at model- and full-scale Reynolds numbers. Gravity waves were neglected, i.e., we adopted the so-called double-model flow. Single-block grids with H–O topology were adopted for all the calculations. Three eddy viscosity models were employed: the one-equation eddy viscosity and the two-equation models proposed by Menter and the TNT version of the two-equation k-ω model. Verification exercises were performed in sets of nearly geometrically similar grids with different densities in the streamwise, normal, and girthwise directions. The friction and pressure resistance coefficients were calculated for different levels of the iterative error and for computational domains of different size. The results show that on the level of grid refinement used, it is possible to calculate the viscous resistance coefficients in H–O grids that do not match the ship contour with a numerical uncertainty of less than 1%. The differences between the predictions of different turbulence models were larger than the numerical uncertainty; however, these differences tended to decrease with increases in the Reynolds number. The pressure resistance was remarkably sensitive to domain size and far-field boundary conditions. Either a large domain or the application of a viscous–inviscid interaction procedure is needed for reliable results. This work was presented in part at the International Conference on Computational Methods in Marine Engineering—MARINE 2007, Barcelona, June 3–4, 2007.     

带附体潜艇绕流场数值计算与验证

采用κ-ε湍流模型和有限体积法用求解RANS方程的方法,对全附体的SUBOFF AFF-8潜艇模型三维粘性绕流场进行了数值模拟计算.一种结合增强型壁面函数的两层模型被采用,以便在近壁区域采用精细网格来模拟边界层内的流场流动特征.给出了雷诺数Re=1.2×107下模型表面纵中剖面线上半部分压力系数与摩擦阻力系数的分布,以及指挥台围壳与尾翼表面不同高度处的压力系数分布,同时计算了桨盘面处的尾流速度分布,并与有关文献给出的实测值进行了比较和验证,分布趋势基本一致,较相关文献给出的计算结果有相当大的改进,此外从指挥台围壳与尾翼后的速度矢量图观察到了与有关文献相同的涡旋流场结构,从而验证了该数值方法的有效性,为潜艇概念设计中评价设计方案优劣提供了重要的技术手段.     

On the role played by turbulence closures in hull shape optimization at model and full scale

 The practical use of automated computational fluid dynamics (CFD)-based design tools in the ship-building industry requires powerful flow solvers which are able to take into account realistic geometries as well as complex physical phenomena, such as turbulence. A shape optimization tool is developed in this framework. A derivative-free optimizer, yielding both flexibility and robustness, is preferred to the classical gradient-based method, which is more difficult to implement and is still limited to only moderately complex problems. The flow solver included in the design procedure solves the incompressible Reynolds-averaged Navier–Stokes equations on unstructured grids using a finite-volume formulation involving several near-wall low-Reynolds-number turbulence models. The design tool is used to optimize the stern of a modern hull shape at model and full scale, with different purposes being considered. More precisely, the drag reduction and the homogenization of the flow in the wake are expected by controlling the longitudinal vortex generated. Our interest is particularly focused on the influence of turbulence modeling in the design process. The effects of a two-equation model based on the eddy-viscosity assumption and a second-order closure relying on the Reynolds stress transport equations are compared. Received: September 24, 2002 / Accepted: April 14, 2003 RID="*" Acknowledgment. The authors thank the scientific committee of CINES (project dmn2050) for the attribution of CPU time.     

Test of large-eddy simulation models for homogeneous-shear turbulence in stratification

Numerical tests of various subgrid-scale (SGS) models were conducted for turbulence in thermally stratified homogeneous-shear flow at a relatively low Reynolds number. Compared with a direct numerical simulation (DNS), we found that nondynamic isotropic SGS models are not able to represent the energy spectrum very well because the energy decays considerably during the transition between an initial random stage and a stage of coherent turbulent structures. Dynamic models performed well for simulating the energy spectrum and the change of GS properties with time; anisotropy is not a necessary feature under the present simulation conditions, although one of the special features of stratified turbulence is anisotropy. This may be because the present grids for large-eddy simulation were fine enough to resolve the patches of counter-gradient heat fluxes, which play an important role in the evolution of turbulent energy in stratified turbulence. With respect to the domain-averaged values of SGS stresses, only the dynamic two-parameter mixed (DTM) model produced results of the same order of magnitude as those of filtered DNS. This is because of the terms arising from re-decomposing of the SGS stresses in the DTM model. It was also found that this incompetence in simulating the SGS stress is not necessary to simulate GS energy evolution, as is known for wall turbulence. Updated from the Japanese original (J Soc Nav Archit Jpn 2001; 190:27–39)     

Hydrodynamic derivatives investigation of unconventionally arranged pusher-barge systems

Unconventional arrangements of pusher-barge systems were studied in this paper. Pusher-barge systems consisting of 4, 6, and 8 barges with one pusher were tested in various combinations. Captive model testing was performed on the various combinations at the Hiroshima University towing tank. Hydrodynamic derivatives of the systems were obtained from the model test data by using the least-square analysis method. For asymmetric conditions, the hydrodynamic derivatives and N′_ββ were added to the force and moment equations in order to obtain better fitting of the least-square curves. Motion equations were modified to cover the asymmetric cases of pusher-barge systems with lateral force and yaw moment due to the asymmetry arrangement. Turning simulations (with 20° sudden angle change) were carried out and a comparison of advance distance and tactical diameter made. An erratum to this article can be found at     

Mixing in the surface boundary layer of a tropical freshwater reservoir

A combined observational-modeling study was conducted to investigate turbulence mixing, and the relation to surface forcing, in the surface boundary layer (SBL) of a tropical, high-altitude, freshwater reservoir. A suite of vertical profiles of temperature microstructure, collected at three different stations of one-day duration each, provided estimates of dissipation rates of turbulence kinetic energy, , and temperature variance, χ. Numerical simulations of and χ, using state-of-the-art, public domain, two-equation turbulence closure models, compared favorably with the observations and reproduced the dynamics of daytime wind mixing as well as the vertical and temporal turbulence structure during nighttime convective conditions.Two independent estimates of vertical eddy diffusivities in the stably stratified (daytime) SBL, computed from the microstructure measurements, agreed closely, and the near surface heat and buoyancy fluxes, computed from the diffusivities, were similar to those computed independently from surface meteorology. Model generated eddy diffusivities agreed closely with the observed values, except those generated by K profile parameterization (KPP) model simulations. The good agreement provides confidence that nutrient fluxes in the SBL may be accurately computed from the models when forced with regularly measured surface meteorological parameters. The consequences are important for estimation of daily primary productivity rates in the euphotic zone and the ability to predict algal blooms such as those observed in the present reservoir.     

舰船气泡尾流数值模拟(英文)

Based on a volume of fluid two-phase model imbedded in the general computational fluid dynamics code FLUENT6.3.26, the viscous flow with free surface around a model-scaled KRISO container ship (KCS) was first numerically simulated. Then with a rigid-lid-free-surface method, the underwater flow field was computed based on the mixture multiphase model to simulate the bubbly wake around the KCS hull. The realizable k-ε two-equation turbulence model and Reynolds stress model were used to analyze the effects of turbulence model on the ship bubbly wake. The air entrainment model, which is relative to the normal velocity gradient of the free surface, and the solving method were verified by the qualitatively reasonable computed results.     

Investigation of cavitation inception characteristics of hydrofoil sections via a viscous approach

Unlike traditional approaches based on potential theory, this work investigated the cavitation inception characteristics of hydrofoil sections by solving Reynolds-averaged Navier–Stokes (RANS) equations. The viscous effects at high Reynolds numbers were taken into account through calculation of the eddy viscosity using a k- model. This viscous approach delivers the lift and drag forces in a more self-contained way, which cannot be obtained directly from the potential theory. A comparison with smooth surface measurements validates the proposed numerical method. Two approaches to include the roughness effects are also described. The predictions of cavitation bucket characteristics based on potential theory are more conservative than those reached by the proposed viscous computations. Hydrofoil sections with different camber ratios, foil thicknesses, and Reynolds numbers were studied numerically. A decrease in the Reynolds number tends to increase the cavitation-free region without changing the shape of the steep region. An increase in foil thickness gives a range increase in the angle of attack where no cavitation is expected. An increase in the camber ratio moves the cavitation-free location of the angle of attack, and slightly increases its range.     

A pseudofield model approach to simulate compartment-fire phenomena for marine fire safety design

In order to prevent the spread of marine compartment fires, it is necessary to understand the governing factors or characteristics of fire-spread phenomena. We present a pseudofield model approach to this problem. We first described a field model of turbulent heat convection based on a standard k – turbulence model. Two-dimensional numerical simulations of a two-linked compartment fire were carried out in order to predict the turbulent convection flow induced by the heat released from the fire. Then a more complicated fire-spread problem of multilinked compartment fires was analyzed by means of a zone model, in which the amounts of oxygen consumption and gas generation were solved by a gas-balance equations system. The effect of threshold conditions on fire propagation and the effect of the thickness of the heat insulation were investigated with numerical simulations.     

Influence of turbulence closure models on the vortical flow field around a submarine body undergoing steady drift

Manoeuvring underwater vehicles experience complex three-dimensional flow. Features include stagnation and boundary layer separation along a convex surface. The resulting free vortex sheet rolls up to form a pair of streamwise body vortices. The track and strength of the body vortex pair results in a nonlinear increase in lift as body incidence increases. Consequently, accurate capture of the body vortex pair is essential if the flow field around a manoeuvring submarine and the resulting hydrodynamic loading is to be correctly found. This work highlights the importance of both grid convergence and turbulence closure models (TCMs) to the strength and path of the crossflow-induced body vortices experienced by the DOR submarine model at an incidence angle of 15°. Five TCMs are considered; Spalart–Allmaras, k-ε, k-ω, shear stress transport, and the SSG Reynolds stress model. The SSG Reynolds stress model shows potential improvements in predicting both the path and strength of the body vortex over standard one- and two-equation TCMs based on an eddy viscosity approach.     

喷水推进器进口流道水动力性能分析

[目的]研究喷水推进器进口流道主参数对其性能的影响,为喷水推进器设计提供依据.[方法]基于STAR-CCM+商业软件,通过定常雷诺平均NS方程(RANS),数值模拟分析不同进速比(IVR)工况下喷水推进器进口流道轴线高度和进流角度对其水动力性能的影响,并根据国际拖曳水池会议(ITTC)的不确定度分析规程进行数值不确定度...     

破损船舶进水模拟

本文应用RANS方程进行了破损舰艇的进水模拟。使用了流体体积法、两方程湍流模型及壁函数处理自由面附近湍流流动,给出了后台阶流动、 系列60以及破损Burgundy船的CFD模拟计算结果,捕捉了船舶破损相关的基本流动特征。     

非均衡湍流边界层脉动压力测试研究

通过建立非均衡湍流边界层脉动压力测试技术及可提供较高测试信噪比的低噪声风洞试验装置,采用表面贴装式硅微传声器线列阵,测量获得了不同压力梯度模型的湍流边界层脉动压力特性,压力梯度模型表面湍流边界层脉动压力低频平台区幅值增加3-5 dB;回归得到了非均衡湍流边界层流脉动压力频谱与Strauhal数及Re数和压力梯度相关的拟合模型,具有一定的适用性;压力梯度对湍流边界层脉动压力频率—波数谱的影响主要集中在300 Hz以下的低频段,传输波数附近峰脊区谱级增加3-5 dB。     

熵格子Boltzmann方法的重整化群代数湍流模型

重整化群理论所建立的湍流模型能够最大程度地减小模型经验性,因此文章尝试将重整化群代数湍流模型引入到熵格子Boltzmann方法中,建立新型的计算模型以对高雷诺数湍流进行模拟研究。同时为了进行比较研究,还建立了熵格子Boltzmann方法的标准大涡模拟模型。完成了对高雷诺数湍流绕流场的模拟计算。结果表明:所建立的熵格子Boltzmann方法重整化群代数湍流模型能够有效地模拟高雷诺数湍流流动问题;其对紧贴壁面处较小尺度湍涡的模拟结果趋近于大涡模拟的结果;重整化群代数湍流模型在对高雷诺数湍流的模拟中表现出耗散模型的特征。     

PIV velocity field measurements of flow around a KRISO 3600TEU container ship model

The main purpose of this investigation was to demonstrate a useful application of the particle image velocimetry (PIV) method to analyze the complex flow characteristics around a ship. For a sample illustration, the KRISO 3600TEU container ship model was chosen. The flow structure in the stern and near-wake region of the model has been investigated experimentally in a circulating water channel. Instantaneous velocity fields measured by the PIV velocity field measurement technique have been ensemble-averaged to give details of flow structures such as the spatial distributions of the local mean velocity, vorticity, and turbulent kinetic energy. The free-stream velocity was fixed at U _o = 0.6m/s, and the corresponding Reynolds number based on the length between perpendiculars was about 9.0 × 10⁵. The container ship model shows a complicated three-dimensional flow structure in the stern and near-wake regions. The PIV results clearly revealed the formation of large-scale bilge vortices in the stern region and their effect on the flow in the near-wake. The results shown here provide valuable information for hull form design and the validation of viscous ship flow codes and of turbulence models.     

An experimental/numerical approach for evaluating skin friction on full-scale ships with surface roughness

For an accurate evaluation of the increase in skin friction due to various surface topographies on ships, i.e. plate roughness, coatings, or bio fouling, both experimental measurements and a numerical evaluation of those measurements are necessary. The measurements are necessary as no other practical method exists to evaluate the skin friction coefficient on most surface topographies, and numerical evaluation is required to compute the roughness effects of varying ship types and speeds. Therefore, a method for measuring the skin friction coefficient for bio-fouled and structured surfaces at full-scale friction velocity is presented, and a validation of the design and measuring procedures is given. For the ship frictional resistance calculations, the Computational Fluid Dynamics (CFD) code ShipFlow has been modified to take the added friction into account.     

散货船阻力预报中湍流模型的应用

采用直接求解RANS方程的方法对20艘散货船的绕流场和阻力进行了数值模拟,比较了标准模型、RNG模型及SST模型对数值计算结果的影响.通过与试验结果的详细对比,发现只要选取合适的湍流模型,数值模拟结果能够达到较高的精度,在1% -2%之间.这表明该数值计算方法能够满足工程的实际使用要求,可以在工程中推广应用.     

叶片数对喷水推进低比转速轴流泵叶轮水动力性能的影响

喷水推进低比转速轴流泵是随着喷水推进技术的发展而产生的一种新泵型。针对一型低比转速轴流泵叶轮叶片数分别选取为5、6、7建立三组叶轮模型,选择RNG k-ε两方程涡粘模型进行叶轮内流场与功率、总压扬程、效率的数值模拟计算,得到三组叶轮的相关水动力曲线。在此基础上分析叶片数对低比转速轴流泵叶轮水动力性能的影响,得到叶片数对功率、总压扬程与工况点影响较大,而基本不影响最高效率的结论,为我国该类泵型的深入研究提供参考。     

Numerical analysis of ducted propeller performance under open water test condition

This paper investigates the open water performance of the Ka-series propellers at various pitch and expanded area ratios in combination with the 19A duct by employing the panel method panMARE and the RANSE code ANSYS-CFX. An efficient method, Caly, is developed in order to generate the 3D-geometry and the surface numerical grid of the ducted propellers. Caly can be coupled with ANSYS-TurboGrid to automatically produce 3D-grids for the RANSE solver. The numerical results are compared with published experimental data and the flow details are concluded and compared. The influences that the grid resolution, the panel arrangements of duct and blade, and the flow in gap between inside wall of the duct and blade tip on the numerical results are studied. Grids verification, turbulence model dependency analysis and Reynolds number scantling are also discussed.     

基于CFD法的对转螺旋桨敞水性能参数匹配研究

摘要：采用雷诺平均纳维-斯托克斯（RANS)法对B型螺旋桨和对转螺旋桨周围的粘性流场进行数值处理,运用Pro/E和GAMBIT软件分别进行三维建模和网格划分,并采用ORIGIN和TECPLOT后处理软件对结果进行图像处理。分析研究湍流模型对螺旋桨水动力性能的影响,对对转桨进行数值计算验证,同时对其敞水性能参数的匹配问题进行探讨,分析其尾流和速度场的变化情况。结果表明：SST 模型为较佳的湍流模型;对转桨在适当的桨距比和直径比的情况下,后桨可充分吸收前桨的尾涡能量,其推进效率明显提高。