Development of hybrid URANS-LES methods for flow simulation in the ship stern area

The paper presents the results of the application of a new hybrid URANS-LES method for the investigations of the ship wake behind the tanker KVLCC2. The switching between URANS and LES models is based on the ratio between the turbulence scale and the cell size of the mesh. Ship resistance, fields of the axial velocity and turbulent kinetic energy in the propeller plane are calculated and compared with measurements. Much attention is paid to the analysis of the unsteady velocities, their PDF distributions and spectra. Numerical analysis shows that the instantaneous velocities deviate substantially from their mean values which are usually used as the estimated velocities in modern engineering methodologies. The thrust variation in the unsteady wake is more than twice as large as that in the time averaged (frozen) wake. The results of the present study point out that the unsteadiness in the wake behind full ships can be very large and should be taken into account when propulsion and unsteady loadings are determined.     

An inverse design approach for minimising wake at propeller plane using CFD

Knowledge of wake characteristics in the stern region is important for ensuring good propeller design and performance. This work examines the utility of CFD in the analysis of flow in the case of full aft beam vessels having characteristic cut stern shape to facilitate propeller aperture. The underwater stern shape may be more complex due to the occurrence of stern appendages such as bossings, strut supports and local shape variations. To this extent, CFD offers an effective tool for both qualitative as well as quantitative assessment of the local geometry. Wake estimate is required for choice of the most favorable propeller geometry. In the present method, the analysis quantifies the effects of small changes in stern rake angles and offers an inverse design approach towards finalising the stern shape. The method consists of solving the standard k-ε turbulent model of RANS equations in cell centered finite volume multi zone grid in the flow domain. This approach has been used in estimating the velocity at the propeller plane. The results have been compared with experimentally obtained values of nominal wake. The approach suggests that CFD can provide a cost effective and quick assessment of flow. It is also an attractive means of pre-empting heterogeneous flow related problems such as vibration and noise due to unfavorable wake in the stern region.     

Numerical Experiment on Two-Dimensional Line Thermal

The time evolution of a two-dimensional line thermal-a turbulent flow produced by an initial element with signifi-cant buoyancy released in a large water body, is numerically studied with the two-equation k - e model for turbulence closure. The numerical results show that the thermal is characterized by a vortex pair flow and a kidney shaped concentra-tion structure with double peak maxima; the computed flow details and scalar mixing characteristics can be described by self-similar relations beyond a dimensionless time around 10. There are two regions in the flow field of a line thermal: a mixing region where the concentration of tracer fluid is high and the flow is turbulent and rotational with a pair of vortex eyes, and an ambient region where the concentration is zero and the flow is potential and well-described by a model of doublet with strength very close to those given by early experimental and analytical studies. The added virtual mass coeffi-cient of the thermal motion is found to be approximat     

Estimates of the characteristics of the horizontal turbulent exchange in the northwestern part of the Pacific Ocean

Estimates of the characteristics of the horizontal turbulent exchange (coefficients of horizontal exchange, scales of deviations of the velocity of geostrophic currents from the background flows, and horizontal scales of perturbations) in the upper ocean layer were obtained on the basis of the data of expeditions held in 1978 and 1980 in the northwestern part of the Pacific Ocean. It is shown that the characteristics of the horizontal turbulent exchange strongly depend on the structure of the background flows, and one of the causes of their variability is related to tropical cyclones. A tendency to a decrease (to different degrees) of the turbulent coefficients in the wake of a tropical cyclone is noted.     

湍能封闭模型在渤海潮流模拟中的应用

将一个三维湍能封闭模型应用于渤海潮流模拟，通过流速分解的初边值方法，计算了渤海全海区潮波系统、潮流结构，结果与观测符合较好。计算同时给出了湍粘性系数和湍动能的时空分布，比较常湍粘性系数和湍的垂直抛物分布及湍能封闭模拟表明，湍粘性系数的不同选取会对流场结构有显著影响。     

The upper ocean response to a moving typhoon

The upper ocean response to the translation speed of typhoons is studied using a three-dimensional primitive equation model. Similar models studied previously have applied stability criteria rather than the diffusion term to simulate the vertical mixing process. This study retains the diffusion term and uses the level-2 turbulence closure scheme to estimate the vertical eddy viscosity. The model results indicate that in the forced period, the mixed-layer temperature decrease is greater for a slow-moving storm due to stronger upwelling caused by the longer residence time. A fast-moving storm can attain a similar cooling intensity in the wake period if its residence time allows the wind to resonate with the current. The significant downward momentum diffusion and advection in the first few inertial periods of these events leads to strong, persistent inertial pumping throughout the upper ocean in the wake period. The mixed layer is further cooled by turbulent mixing supported by vertical current shears. Meanwhile, the upper thermocline exhibits a compensating temperature increase. The vertical transfer magnitude and penetration scale are smaller in the slow-moving case, when the inertial motion decays rapidly. The model results also indicate that the dominant cooling process can be inferred from the non-dimensional storm speed. However, this value may be misleading for rapidly moving storms in which the current response is so distant from the storm that little wind work is performed on the ocean.     

Laboratory, numerical, and theoretical modeling of the flow in a far wake in a stratified fluid

The far-wake flow past a sphere towed in a fluid with high Reynolds and Froude numbers and with a pycnocline-form salt-density stratification is studied in a laboratory experiment based on particle image velocimetry and in numerical and theoretical modeling. In the configuration under consideration, the axis of sphere towing is located under a pycnocline. Flow parameters, the profiles of density and average velocity, and the initial field of velocity fluctuation in numerical modeling are specified from the data of the laboratory experiment. The fields of fluid velocity at different times and the time dependences of integral parameters of wake flow, such as the average velocity at the axis and the transverse width of the flow, are obtained. The results of numerical modeling are in good qualitative and quantitative agreement with the data of the laboratory experiment. The results of the laboratory experiment and numerical modeling are compared to the predictions of a quasi-linear and quasi-two-dimensional theoretical model. The time evolution of both the average velocity at the axis and the transverse width of the wake is obtained with the model and is in good agreement with the experimental data. The results of numerical modeling also show that, under the effect of velocity fluctuation in the wake, internal waves whose spatial period is equal to the characteristic period of the wake’s vortex structure are excited efficiently in the pycnocline.     

两层粘性流体中圆柱体受迫振荡数值模拟

研究两层粘性流体中无限长水平圆柱体的受迫振荡问题。在湍流模式下,采用VOF方法追踪两层流体的内界面,基于动网格技术模拟圆柱体的运动边界,对均匀流中横向振荡圆柱体的绕流场进行了数值模拟。计算受迫振荡圆柱体的升力系数、阻力系数随时间的演化曲线和圆柱体的尾涡分布,以及圆柱体的受迫振荡激发两层流体内界面的扰动,并与均匀流体的情况进行了比较分析。研究表明,流体的两层分层效应对受迫振荡圆柱体的升阻力系数和尾涡分布特性都有显著影响,在水下输油气管道涡激振动特性的工程评估中,应考虑流体的密度分层效应。     

Analysis of wake behind a rotating propeller using PIV technique in a cavitation tunnel

A two-frame particle image velocimetry (PIV) technique is used to investigate the wake characteristics behind a marine propeller with 4 blades at high Reynolds number. For each of 9 different blade phases from 0° to 80°, 150 instantaneous velocity fields are measured. They are ensemble averaged to study the spatial evolution of the propeller wake in the region ranging from the trailing edge to one propeller diameter (D) downstream location. The phase-averaged mean velocity shows that the trailing vorticity is related to radial velocity jump, and the viscous wake is affected by boundary layers developed on the blade surfaces and centrifugal force. Both Galilean decomposition method and vortex identification method using swirling strength calculation are very useful for the study of vortex behaviors in the propeller wake region. The slipstream contraction occurs in the near-wake region up to about X/D=0.53 downstream. Thereafter, unstable oscillation occurs because of the reduction of interaction between the tip vortex and the wake sheet behind the maximum contraction point.     

带前缘结节叶片导管桨尾流不稳定性分析

导管桨的尾流不稳定性在其性能评价中非常重要,不但是其能否提供稳定推力的保证,而且也与螺旋桨的尾流噪声直接相关。为了改善导管桨的尾流,提高尾流稳定性,并优化导管桨的流场脉动,根据座头鲸鳍肢前缘结节的仿生原理,对导管桨叶片的导边进行改进,提出了两种仿生桨型,采用IDDES湍流模型对低进速系数下常规导管桨和仿生叶片导管桨进行数值模拟,探究叶片构型对导管桨性能和尾流不稳定性的影响。计算结果表明,前缘结节可以有效降低叶片受力波动的幅值和叶片所受合力的主频域峰值,具有较大结节的叶片对导管桨尾流有明显的优化作用,在尾流远场中扩大了流动稳定区,延后了尾流处涡破碎的发生,改善了能量谱密度的频域分布。进一步,大前缘结节叶片导管桨应用在低速工况下时,可以大量减少尾流泄涡区域的二次涡产生,这是由于前缘结节提升了相邻涡互感的强度,使得尾流更加稳定,而小结节叶片仿生桨型对导管桨尾流则无明显优化作用。研究方法和成果可为螺旋桨尤其是导管桨尾流不稳定性研究提供参考,不仅验证了前缘结节在导管桨叶片应用的合理性,而且揭示了其优化尾流稳定性的机理。     

基于尾流模型的风场偏航控制优化研究

在风机大尺度化与风场大型化的趋势下,如何通过合适的控制策略以降低尾流损失成为关键问题。以包括30台NREL-5MW风机并采用5行6列平行四边形布置方式的小型风场为研究对象,基于显示尾流模型,以各风机偏航角度为优化参数,风场总功率为目标函数,使用粒子群优化算法对比分析了偏航控制对不同风速、风向、湍流强度下的风场性能提升效果。结果表明,偏航控制优化可在风向与风机行或列方向平行时发挥明显效果,当风机行列间距为4倍风轮直径且湍流强度为5%时,在不同风速下偏航控制可分别将风场总体发电量提升15%~20%,但对于布置间距大于7倍风轮直径或湍流强度高于15%时的风场,其作用十分有限,总体发电量提升在5%以内。     

Line broadening on HF ocean surface radar backscatter spectra

Experimental results from an array of moored current meters and an HF ocean surface radar support the idea that line broadening on the radar spectra is caused by the velocity distribution within the radar target cell. The experiment was done in the wake of a small island where the velocity variations were severe. An estimate is made of the line broadening which can be expected. In a turbulent flow with dissipation rate of the orderepsilon sim 10^{-10}m^{2}s^{-3}and target cell size 1 3000 m, the line broadening isDeltaf sim 10^{-3}Hz. This would be resolved with a radar time series ofsim 20min and indicates that the HF ocean surface radar technique has potential in the observation of surface velocity distributions.     

Mathematical modeling of turbulent jets of deep-water sewage discharge into coastal basins

The basic properties of the dynamic model of a turbulent jet formed by a deep-water sewage discharge into the stratified environment of coastal regions are considered. The model developed was used to estimate the parameters of a floating-up jet of deep wastewater discharge from Sand Island into the basin of Mamala Bay (Hawaii) depending on the season and discharge operation mode. The estimates of the float-up depths of the jet and the initial dilution of the jet were estimated on the basis of model calculations using experimental data on the vertical profiles of the water temperature and salinity under the actual conditions of stratification in the study region. It is shown that the further propagation of the wastewater jet depends on tidal events and internal waves generated by tides. The appearance of turbulent jets at the sea surface was recorded. The model estimates of the parameters of the wastewater discharge were compared with the results of experimental measurements. Good agreement was found, which indicates that the physical mechanisms of the propagation of turbulent jets in a stratified medium are adequately described by the model.     

Simulation model of SAR remote sensing of turbulent wake of semi-elliptical submerged body

1 IntroductionThe study of SAR remote sensing of surface wake of a semi-elliptical submerged body is of great importance actually in the ocean. The dynamics and the SAR remote sensing mechanism of this wake are very complicated. Numerous researches (Bergmann et al., 1958; Miles, 1968a,b; Huppert and Miles, 1969; Bon-neton et al., 1993) indicate that in the ocean the imagery and dynamics of this wake are quite different when a semi-elliptical submerged body moves horizontally at high speed …     

Attenuation Measurements Across Surface-Ship Wakes and Computed Bubble Distributions and Void Fractions

A surface ship's wake is composed of several hydrodynamic phenomena. A large part of that wake contains a mixture of air bubbles of various sizes in turbulent water. Eventually, as the wake ages, the turbulence subsides and bubbles begin to rise at rates that are determined by their sizes. These bubbles of various sizes and concentrations control the propagation of acoustic signals inside and across a wake. To further our understanding of these phenomena, a series of three continuous-wave (CW)-pulsed signals were transmitted across a wake as the wake aged. Each transmission contained a set of four 0.5-ms-long pulses. The 12 pulses ranged over frequencies from 30 to 140 kHz in 10-kHz steps. The acoustic attenuations across wakes that were due to varying bubble-size densities within the wakes were determined experimentally. From those data, estimates of the bubble densities as functions of the speed of the wake-generating ship, the wake's age, and acoustic frequency were calculated. From the bubble-density results, power-law fits and void fractions are calculated. The attenuation measurements were taken at 7.5-m intervals behind the wake-generating ship and continued for about 2 km. The experiment was run for wakes generated at ship speeds of 12- and 15-kn wakes, and the 15-kn run was repeated for consistence determination. The bubble densities were observed to have power-law forms with varying parameters with the strongest, for early ages, having an exponent of ${-}$3.6 and a void fraction of 4 $times$ 10$^{-7}$ , and with both diminishing for older wakes, as might be expected.      

On the scale of the turbulent small masses of sea water

The relation between the initial spatial scale and the life time, of turbulent small masses of sea water produced at the strait or the interface between two water masses, is investigated analytically. In the analysis, their initial shape is assumed to be expressed as the Rankine vortex, and the horizontal eddy viscosity coefficient to obey then-th power-law with respect to scale. Particularly in the case ofn=4/3 (the 4/3 power-law), it is obtained that the life time is proportional to the 2/3 power of the initial scale. The proportional coefficient can be determined in the present analysis except a parameter related to the energy dissipation rate of each sea area.     

Experimental analysis of the flow field around horizontal axis tidal turbines by use of scale mesh disk rotor simulators

Understanding the flow field around horizontal axis marine current turbines is important if this new energy generation technology is to advance. The aim of this work is to identify and provide an understanding of the principal parameters that govern the downstream wake structure and its recovery to the free-stream velocity profile. This will allow large farms or arrays of devices to be installed whilst maximising device and array efficiency. Wake characteristics of small-scale mesh disk rotor simulators have been measured in a 21 m tilting flume at the University of Southampton. The results indicate that wake velocities are reduced in the near wake region (close behind the rotor disk) for increasing levels of disk thrust. Further downstream all normalised wake velocity values converge, enforcing that, as for wind turbines, far wake recovery is a function of the ambient flow turbulence. Varying the disk proximity to the water surface/bed introduces differential mass flow rates above and below the rotor disk that can cause the wake to persist much further downstream. Finally, the introduction of increased sea bed roughness whilst increasing the depth-averaged ambient turbulence actually decreases downstream wake velocities. Results presented demonstrate that there are a number of interdependent variables that affect the rate of wake recovery and will have a significant impact on the spacing of marine current turbines within an array.     

Influence of propeller wake sheet on rudder gap flow and gap cavitation

The characteristics of the flow over the rudder’s pintle gap are investigated by using the particle image velocimetry (PIV) technique. The propeller and rudder models are scaled down to 1/28.5. Highly accelerated leakage outflows are separated at the discontinuities of the gap and generate strong cavitation at the suction side of the rudder. In the rudder and propeller configuration, the propeller wake sheet ahead of the gap entrance region starts to induce leakage flow over the lower pintle gaps of the suction side. The gap flow has a velocity magnitude as high as 0.4U₀ in the high leakage flow condition, where the wake sheet locates over the gap entrance. The cross-flow of the propeller wake sheet interferes the gap entrance region and triggers gap cavitation. As the propeller wake sheet moves downstream and weakens, the gap flow velocity decreases over the gap entrance.     

Giant comet marks

Large longitudinal bedforms (comet marks) occur on the current-swept shelf off southeast Africa. They emanate from tight clusters of low rock outcrops and are as long as 135 m and as wide as 30 m. The low length/width ratio (1:4.5) is explained by lateral coalescence of several narrow tails associated with individual obstacles. They seem to form by removal of sediment from the wake region of the obstacles. Hydrodynamic considerations indicate a high probability for sonic suspension transport, yet no sand shadows occur in the velocity-defect region of the inner wake. It is concluded that shear velocities must exceed the minimal velocity required to sustain sediment transport.     

Noise Control Due to the Stator Wake Blade Interaction via Tail Articulation

The biologically inspired method of tail articulation is investigated as a means of reducing tonal noise due to wake deficit blade interaction in underwater vehicles. Experiments are carried out in a water tunnel under typical operating conditions for underwater vehicles. Tail articulation is implemented using a life scale stator model with a hinged flapping tail operating both in free-stream velocities corresponding to Reynolds number in the range 75000 < Re < 300000 and at frequencies up to 30 Hz to investigate the range of Strouhal number 0.0 < St < 0.35. Velocity measurements of the active stator wake are carried out by laser Doppler velocimetry (LDV) and particle image velocimetry (PIV) to investigate the effects of tail articulation on the stator wake. Time-averaged measurements of the stator wake by LDV show that of the tail articulation has a dominant effect on the time mean stator drag. Instantaneous phase-averaged measurements of the stator wake by PIV show a transition in the unsteady stator wake as is increased, from a deflected vortex sheet to a series of rolled up, discrete vortices. Measurements are made of the wake due to both sinusoidal and nonsinusoidal tail motion profiles, which show that significant wake alteration is achieved with tail articulation. A low-order model describing the creation and convection of vorticity by tail articulation is developed which describes wake phenomena observed in LDV and PIV measurements. Finally, a 3-D unsteady propeller simulation using both experimental wake velocity data by PIV and simulated wake velocity data generated with the reduced-order model are used to predict the effect of sinusoidal tail articulation on radiated noise. Results using simulated data indicate that a significant noise alteration is achieved in all cases, and noise reduction of 5-8 dB is achieved in some cases.