PIV measurements of hull wake behind a container ship model with varying loading condition

Flow characteristics of the hull wake behind a container ship model were investigated under different loading conditions (design and ballast loadings) by employing the particle image velocimetry (PIV) technique. Measurements were made at four transverse locations and two longitudinal planes for three Reynolds numbers (Re) (=U₀Lpp/ν, where U₀ is the freestream velocity, Lpp is the length between two perpendiculars of the ship model and ν is the kinematic viscosity) of 5.08×10⁵, 7.60×10⁵, and 1.01×10⁶. It was observed that symmetric, large-scale, longitudinal counter-rotating vortices (with respect to centerline) of nearly the same strength were formed in the near wake. For the ballast-loading condition, the vortices appear at propeller plane below the propeller-boss. The vortex center exhibits a significant upward shift near the propeller-boss as the Reynolds number increase, and as the flow moves downstream. Under the design-loading condition, the vortices first appear at a further downstream location than that for the ballast-loading condition above the propeller-boss. This difference in the flow structure can significantly change the inflow conditions to the propeller blades, such as the streamwise mean velocity profiles and turbulence intensity distributions at the propeller plane. In particular, under the ballast-loading condition, asymmetric inflow may weaken the propulsion and cavitation performance of the marine propeller.     

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.     

Comparison of flow patterns in the near wake of a circular cylinder and a square cylinder placed near a plane wall

The flow patterns in the near wake of a cylinder (either circular or square in shape, D=25 mm) placed in the proximity of a fully developed turbulent boundary layer (thickness δ=0.4D) are investigated experimentally using particle image velocimetry (PIV). The effects of changing the gap height (S) between the cylinder bottom and the wall surface, over the gap ratio range S/D=0.1–1.0, have been investigated. The results show that both the ensemble-averaged and instantaneous flow fields are strongly dependent on S/D. The flow patterns for the two types of cylinders share many similarities with respect to the change in S/D, such as the reduced recirculation length and increased velocity fluctuation in the near wake with increasing S/D, as well as the trend of suppression of vortex shedding at small S/D and onset of vortex shedding at large S/D. However, developments of the shear layers, in terms of wake width, flow curvature, etc., are considerably different for these two types of cylinders. In general, the wake development and momentum exchange for the square cylinder are slower those for the circular cylinder at the same gap ratio. Correspondingly, it is shown that the periodic vortex shedding is delayed and weakened in the case of square cylinder, as compared to that of the circular cylinder at the same S/D.     

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

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

Numerical investigation on the hydrodynamic characteristics of a marine propeller operating in oblique inflow

The hydrodynamic characteristics of a marine propeller operating in oblique inflow are investigated by using CFD method. Two propellers with different geometries are selected as the study subjects. RANS simulation is carried out for the propellers working at a wide range of advance coefficients and incidence angles. The effects of axial inflow and lateral inflow are demonstrated with the hydrodynamic force on the propeller under different working conditions. Based on the obtained flow field details, the hydrodynamic mechanism of propeller operating in oblique inflow is analyzed further. The trailing vortex wake of propeller is highly affected by the lateral inflow, resulting in the deflected development path and the circumferentially non-uniform structure, as well as the enhanced axial velocity in slipstream. Different flow patterns are observed on the propeller blade with the variation of circumferential position. Combined with the computed hydrodynamic forces and pressure distribution on propeller, the mechanism resulting in the increase of propulsive loads and the generation of propeller side force is explored. Finally, a systematic analysis is carried out for the propulsive loads and propeller side force as a function of axial and lateral advance coefficients. The major terms that play a dominant role in the modeling of propulsive loads and propeller side force are determined through the sensitivity analysis. This study provides a deeper insight into the hydrodynamic characteristics of propeller operating in oblique inflow, which is useful to the investigation of propeller performance during ship maneuvers.     

A three-dimensional wake impingement model and applications on tandem oscillating foils

Potential flow based vortex numerical methods have been widely used in aerodynamics and hydrodynamics. In these methods, vortices shed from lifting bodies are traced by using vortex filaments or dipole panels. When the wake elements encounter a downstream body, such as a rudder behind a propeller or a stator behind a rotor, a treatment is necessary to divert the wake elements to pass by the body. This treatment is vital to make wake simulations realistic and to satisfy the non-penetration condition during wake body interaction. It also helps to avoid pure numerical disturbances such as when a vortex filament or an edge of a dipole panel passes through the collection point of a body element; this is a singularity for induced velocity and it will introduce a large numerical disturbance. This necessary treatment for three-dimensional problems with geometrical complexity has not been found to date. In this study, a wake impingement model was developed to divert wake elements to slip over the body surface, model the vortex/body interaction, and predict forces on fluctuating components. The model was also tested on configurations of oscillating foils in tandem with an existing panel method code. Simulation results with the wake impingement model are shown to be in closer agreement with limited published experimental data than those without the model. With the established wake impingement model, force fluctuations on the after body due to the wake vortex impingement were investigated based on a series of simulations. The series varied several parameters including distance between two foils, oscillating frequency, span, rear foil pitch angle, swap angle and vertical position.     

Savonius型转轮三角阵列减流性能研究

为了助力海洋牧场减流防护工程, 研究Savonius型转轮阵列减流性能。作者建立Savonius型转轮三角阵列尾流场数值模型, 并通过水池实验验证准确性, 基于可靠数值模型探究转轮阵列尾涡减流机理, 研究三角阵列结构参数L_X、L_Y, 以及动力参数TSR、初始流速、旋向对整体减流性能的影响规律。结果表明,下游转轮产生的涡流呈现非对称分布, 并且产生更多涡流的转轮拥有更好的减流效果。另外, L_X为3D和L_Y为2D时减流性能最佳。最后对比发现, 在叶尖速比为0.9~1.1减流效果更好; 初始流速大小不影响减流效果; 下游转子对称分布时, 随着上游转子改变旋转方向, 减流效果出现明显差异。     

The virtual mass of a rectangular flat plate of finite aspect ratio

The skin friction of a two-dimensional planing flat plate is made up of two opposing components; a drag force from the flow aft of the stagnation line and an opposing thrust force from the jet flow. This paper is concerned only with the drag term and the wake velocity defect which it causes in the water behind the transom.It is concluded that the skin friction is less than would be expected from a flat plate at ambient static pressure (D_fo say) and is approximately equal to D_fo (1 − C_R), where C_R is the normal force coefficient based on wetted area. The wake velocity decrement due to this drag is found to be significant, particularly for surface piercing propellers.     

Numerical simulation of viscous cavitating flow around a ship propeller

In the present study, cavitation and a ship propeller wake are reported by computed fluid dynamics based on viscous multiphase flow theory. Some recent validation results with a hybrid grid based on unsteady Navier-Stokes (N-S) and bubble dynamics equations are presented to predict velocity, pressure and vapor volume fraction in propeller wake in a uniform inflow. Numerical predictions of sheet cavitation, tip vortex cavitation and hub vortex cavitation are in agreement with the experimental data, same as numerical predictions of longitudinal and transversal evolution of the axial velocity. Blade and shaft rate frequency of propeller is well predicted by the computed results of pressure, and tip vortex is the most important to generate the pressure field within the near wake. The overall results indicate that the present approach is reliable for prediction of cavitation and propeller wake on the condition of uniform inflow.     

船舶螺旋桨尾流场的数值分析

利用基于速度势的低阶面元法计算船舶螺旋桨的尾流场。采用计算较为简捷的关于扰动速度势的基本积分微分方程，并采用双曲面形状的面凶以消除面元间的缝隙。Newton-Raphson迭代过程被用来在桨叶随边满足压力Kutta条件，使桨叶面上表面的压力在随边有良好的一致性。在计算面元的影响系数时，应用了Morino导出的解析计算公式，加快了数值计算的速度。从解面元法的基本积分方程得到的偶极强度和源汇强度，直接求得尾流场的速度分布。     

On the fluid structure interaction of a marine cycloidal propeller

Marine cycloidal propulsion system is efficient in maneuvering ships like tugs, ferries, etc. It is capable of vectoring thrust in all direction in a horizontal plane. When used in pair, the system enables a vessel to perform maneuvers like moving sideways, perform rotation about a point, i.e. turning diameter of its own length, etc. In this system, the propeller blades have to change their angle of attack at different angular position of the disc. Due to this reason, the inflow velocity vector to propeller blades changes continuously. The propeller blade oscillates about a vertical axis passing through its body and at the same time rotates about a point. Superposed on these motions is the dynamics of the ship on which the propulsion system is installed. This results in a formidable and challenging hydrodynamics problem. Each of the propeller blade sections could be considered as an aerofoil operating in combined heave and pitch oscillation mode. Due to the constantly varying inflow velocity, the hydrodynamic flow is unsteady. The unsteady hydrodynamic flow is simulated by incorporating the effect of shed vortices at different time instant behind the trailing edge. Due to the kinematics of the problem, the blade is subjected to higher structural deformation and vibration load. The structural deformation and vibration when coupled with the hydrodynamic loading add another level of complexity to the problem. In this paper, the variation of hydrodynamic load on the propeller blade due to steady and unsteady flow is compared. We also model the structural dynamics of the blade and study its effect on the hydrodynamic loading. Finally, we couple the structural dynamics with hydrodynamics loading and study its influence on the propeller blade for different operating regimes.     

Flow induced vibrations of a sharp edge square cylinder in the wake of a circular cylinder

The response of an oscillating circular cylinder at the wake of an upstream fixed circular cylinder was classified by different researchers as galloping, wake induced galloping or wake induced vibration. Furthermore it is already known that a sharp edge square cylinder would undergo galloping if it is subjected to uniform flow. In this study the influence of the wake of a fixed circular cylinder on the response of a downstream square cylinder at different spacing ratios (S/D = 4, 8, 11) is experimentally investigated. The subject appears not to have received previous attention. The lateral displacements, lift forces and the pressure data from gauges mounted in the wake of the oscillating cylinder are recorded and analyzed. The single degree of freedom vibrating system has a low mass-damping parameter and the Reynolds number ranges from 7.7 × 10² to 3.7 × 10⁴.In contrast to that for two circular cylinders in tandem arrangement, the freely mounted downstream square cylinder displays a VIV type of response at all spacing ratios tested. There is no sign of galloping or wake induced galloping with the square cylinder. With increase at the spacing ratio the cross-flow oscillations decrease. It is shown that the vortices arriving from the upstream fixed circular cylinder play a major role on the shedding mechanism behind the downstream square cylinder and cause the square cylinder to shed vortices with frequencies above Strouhal frequency of the fixed square cylinder (St = 0.13). The VIV type of oscillations in the downstream square cylinder is most probably caused by the vortices newly generated behind the square cylinder.     

The distribution and population dynamics of O-group plaice (Pleuronectes platessa L.) on nursery grounds in the Firth of Forth

Midsummer (1 August) population estimates of about 2 million O-group plaice (Pleuronectes platessa L.) were derived for sandy bays around the Firth of Forth in 1979–1980. This is an order of magnitude less than similar estimates made for the Clyde Sea Area in 1973–1974. Autumn population estimates of 0·4–1·0 million fish were comparable to estimates by the Ministry of Agriculture, Fisheries and Food for the area between the Scottish border and Flamborough Head (2·3 million for 1970 and 1973) which represented 4·8% (1973) to 5·3% (1970) of the total number of O-group fish on the English east coast.Largo Bay was the most important nursery area holding 25% of the total population. It is particularly well situated to receive newly metamorphosed plaice carried in water currents along the north side of the Forth from the spawning ground off Fife Ness. Plaice in the Forth are mainly distributed on fine to medium sandy beaches (186–480 μm), the mean number per haul in midsummer (D) being correlated with the median diameter (m.d. in μm) of the low water sediments by the equation: D=−45·7666+0·2327 m.d. (n=11,r=0·68,P<0·02 but>0·01).The shallow inshore water in sandy bays in the outer Firth was well mixed and more marine than estuarine (27·7–35·0‰). The correlation coefficient between fish density and water temperature was low, while that with salinity (S‰) was: D=6·1618+0·2238S (n=23,r=0·62,P<0·005).Regression analysis demonstrated that the relationship between the instantaneous mortality rate (Z) and the initial population density (D_p) was: Z×100=0·7480+0·0546d_p (n=12,r=0·87,P<0·001).The mean mortality rate for the O-group plaice in the Forth nursery areas was 53% month⁻¹.     

The effect of shroud on vortex shedding mechanism of cylinder

In the present study, flow characteristics were investigated experimentally using particle image velocimetry technique (PIV) in a gap between a solid cylinder and a shroud to reveal the effect of shroud diameter (D_s) and porosity (β) on the vortex shedding mechanism of the cylinder. Porosity (varied from β = 0.3 to 0.7) and diameter ratio (D/D_s = 0.4, 0.5 and 0.6) were main parameters examined at a Reynolds number of Re = 5000. For the porosity values of β ≤ 0.5, it is observed that vortex formation of the cylinder occurs only in the gap and shroud produces its own wake flow patterns. Penetrating flow through the shroud extends the shear layers on the both sides of the shroud through the downstream direction and prevents the interaction of shear layers in the near wake region. The diameter ratio and the porosity are impactful on the wake flow patterns in outer region of the shroud since they are determinant of the penetrating flow rate. Force measurements were also performed in the air tunnel in order to reveal the effect of shroud on the drag coefficient of cylinder. It is found that the drag coefficient of the cylinders are reduced significantly by shrouds when compared with that obtained from the bare cylinder case. However, the drag coefficient of the cylinder together with the shroud is higher than the bare cylinder for all cases since the shrouds enlarge the area exposed to the flow.     

Effects of mass and damping ratios on VIV of a circular cylinder

In this study the basic characteristics of the dynamic response and vortex shedding from an elastically mounted circular cylinder in laminar flow is numerically investigated. The Reynolds number ranges from 80 to 160, a regime that is fully laminar. The governing equations of fluid flow are cast in terms of vorticity. The two-dimensional vorticity transport equation is solved using a vortex method. Effects of important parameters on the system response and vortex shedding are investigated; these include: mass ratio, damping ratio, Reynolds number and reduced velocity. The numerical results show that a decrease in either the mass ratio or damping ratio of the system can lead to an increase in both the oscillation amplitude and the reduced velocity range over which lock-in occurs. The results also suggest that the mass-damping parameter may characterize the system response adequately, although the effect of changing mass ratio appears to be a little more profound compared to damping ratio. Vorticity contour plots suggest that the vortex shedding occurs in the 2S mode, although a wake structure similar to the C(2S) mode appears at distances 15–20 diameters downstream in the lock-in region. The simulation results are in good agreement with previously published data.     

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.     

Experimental analysis on the risk of vortex ventilation and the free surface ventilation of marine propellers

The paper presents a discussion of the ventilation inception and air drawing prediction of ships propellers, aiming to predict under what conditions ventilation will happen, and the actual physical mechanism of the ventilation.Three different types of ventilation inception mechanisms are included in our discussion: free surface vortex ventilation, ventilation by sucking down the free surface without forming a vortex as well as ventilation by propeller coming out of the water. Ventilation prediction is based on a series of model tests, where the propeller is tested in different levels of intermittent ventilation. The use of underwater video gives a visual understanding of the ventilation phenomena.Ventilation by vortex formation has analogies with other phenomena, such as the inlet vortex in pump sumps, ground vortex at the inlet of the aircraft engines and the Propeller Hull Vortex Cavitation (PHVC). The paper includes comparison between Propeller Hull Vortex Cavitation (PHVC) and Propeller Free Surface Vortex Ventilation (PFSVV) as well as comparison between PFSVV and vortex formations of aero engines during high power operation near a solid surface.Experimental data based on several different model tests shows the boundary between the vortex forming, non-vortex forming and free surface ventilation flow regimes. For comparison the following parameters, which determined the intensity of the hydrodynamic interaction between the propeller and free surface have been used: propeller load coefficient c_T, tip clearance ratio c/D, propeller submergence ratio h/R, ambient velocity V_i and flow cavitation/ventilation number σ_cav/σ_vent.     

Experimental study on effect of a new vortex control baffler and its influencing factor

The horseshoe vortex generated around the sail-body junction of submarine has an important influence on the non-uniformity of submarine wake at propeller disc.The flow characteristics in the horseshoe vortex generated area are analyzed,and a new method of vortex control baffler is presented.The influence of vortex control baffler on the flow field around submarine main body with sail is numerically simulated.The wind tunnel experiment on submarine model is carried out,and it is proved that the vortex control baffler can weaken the horseshoe vortex and decrease the non-uniformity of the wake at propeller disc.It is shown from the experiment results that the effect of vortex control baffler depends on its installation position;with a proper installation position,the non-uniform coefficient of submarine wake would be declined by about 50%;the Reynolds number of submarine model has an influence on the effect of vortex control baffler too,and the higher the Reynolds number is,the better the effect of the vortex control baffler is.     

The impact of groundwater discharges on mercury partitioning, speciation and bioavailability to mussels in a coastal zone

The Mussel Watch program conducted along the French coasts for the last 20 years indicates that the highest mercury concentrations in the soft tissue of the blue mussel (Mytilus edulis) occur in animals from the eastern part of Seine Bay on the south coast of the English Channel, the “Pays de Caux”. This region is characterized by the presence of intertidal and submarine groundwater discharges, and no particular mercury effluent has been reported in its vicinity. Two groundwater emergence systems in the karstic coastal zone of the Pays de Caux (Etretat and Yport with slow and fast water percolation pathways respectively) were seasonally sampled to study mercury distribution, partitioning and speciation in water. Samples were also collected in the freshwater–seawater mixing zones in order to compare mercury concentrations and speciation between these “subterranean” or “groundwater” estuaries and the adjacent macrotidal Seine estuary, characterized by a high turbidity zone (HTZ). The mercury concentrations in the soft tissue of mussels from the same areas were monitored at the same time.The means of the “dissolved” (< 0.45 μm) mercury concentrations (HgT_D) in the groundwater springs were 0.99 ± 0.15 ng l^− 1 (n = 18) and 0.44 ± 0.17 ng l^− 1 (n = 17) at Etretat and Yport respectively. High HgT_D concentrations were associated with strong runoff over short water pathways during storm periods, while low concentrations were associated with long groundwater pathways. Mean particulate mercury concentrations were 0.22 ± 0.05 ng mg^− 1 (n = 16) and 0.16 ± 0.10 ng mg^− 1 (n = 17) at Etretat and Yport respectively, and decreased with increasing particle concentration probably as a result of dilution by particles from soil erosion. Groundwater mercury speciation was characterized by high reactive-to-total mercury ratios in the dissolved phase (HgR_D/HgT_D: 44–95%), and very low total monomethylmercury concentrations (MMHg < 8 pg l^− 1). The HgT_D distributions in the Yport and Etretat mixing zones were similar (overall mean concentration of 0.73 ± 0.21 ng l^− 1, n = 43), but higher than those measured in the adjacent industrialized Seine estuary (mean: 0.31 ± 0.11 ng l^− 1, n = 67). In the coastal waters along the Pays de Caux dissolved monomethylmercury (MMHg_D) concentrations varied from 9.5 to 13.5 pg l^− 1 (2 to 8% of the HgT_D). Comparable levels were measured in the Seine estuary (range: 12.2– 21.1 pg l⁻¹; 6–12% of the HgT_D). These groundwater karstic estuaries seem to be mostly characterized by the higher HgT_D and HgR_D concentrations than in the adjacent HTZ Seine estuary. While the HTZ of the Seine estuary acts as a dissolved mercury removal system, the low turbid mixing zone of the Pays de Caux receives the dissolved mercury inputs from the groundwater seepage with an apparent Hg transfer from the particulate phase to the “dissolved” phase (< 0.45 μm). In parallel, the soft tissue of mussels collected near the groundwater discharges, at Etretat and Yport, exhibited significantly higher values than those found in the mussel from the mouth of the Seine estuary. We observe that this difference mimics the differences found in the mercury distribution in the water, and argue that the dissolved phase of the groundwater estuaries and coastal particles are significant sources of bioavailable mercury for mussels.     

On a translating and transversely oscillating cylinder: Part 2: Effect of the velocity ratio on the hydrodynamic forces and the near-wake structure

The unsteady flow past a circular cylinder which starts translating and transversely oscillating from rest in a viscous fluid is investigated at a Reynolds numbers of R=10³ and at a Strouhal number of π/4 and for the maximum oscillatory to translational velocity ratios between 0.1 and 1.0. This study is based on numerical solutions of the two-dimensional unsteady Navier–Stokes equations. The object of the study is to examine the effect of increase of velocity ratio on the near-wake structure as well as the hydrodynamic forces acting on the cylinder. For all velocity rates a periodic structure of vortex evaluation and shedding develops which is repeated exactly as time advances. Vortex dynamics close behind the body are affected by changing acceleration of the cylinder and a changeover from one mode to a different mode of vortex formation is observed with increase in velocity ratio. A comparison of the present results with the impulsively started translating case has been included to illustrate the effect of velocity ratio on drag at small values of velocity ratio.