A critical CAE analysis of the bottom shape of a multi stepped air cavity planing hull

One of the most important aspects, in the ACS and multi stepped hull design, is the choice of the geometrical shape of the cavity and the steps. In this article a complete experimental and numerical campaign on a multi stepped ACS has been carried out, varying the velocity and the air flow rate under the hull. The experimental tests have been conducted in an ITTC Towing Tank allowing to validate the numerical simulations obtained by means of a CFD U-RANSe (Unesteady Reynolds Averaged Navier-Stokes equations) code. The CFD setup is described in detail. From this campaign a critical analysis of the Froude number influence on the air cushion shape has been argued. The authors identified four different behaviours, from low to very high Froude numbers. The use of CFD has allowed to observe quantities of difficult evaluation by means of traditional experimental test, as e.g. the frictional component of the resistance, the airflow path lines and the volume of fraction in transversal and longitudinal sections. The results have been discussed.     

Numerical simulation of micro-bubble drag reduction using population balance model

The phenomenon of drag reduction by the injection of micro-bubbles into turbulent boundary layer has been investigated using an Eulerian-Eulerian two-fluid model. Multiple-size group (MUSIG) based on population balance models, which resolve a wide range of bubble sizes taking into account the bubble break-up and coalescence have been used for this purpose. The simulated results are compared against the experimental findings of Madavan et al. [1984. Reduction of turbulent skin friction by micro-bubbles. Physics of Fluids 27, 356-363] and also other numerical studies explaining the sophisticated phenomena of drag reduction. For the two Reynolds number cases considered, the buoyancy with the plate on the bottom configuration is investigated, as from the experiments it is seen that buoyancy seem to play a role in the drag reduction. Numerical model employed in the investigation comprises of a micro-bubble laden flow wherein two independent sets of Reynolds averaged Navier-Stokes (RANS) transport equations were used to describe both the phases of the flow. The shear stress transport (SST) turbulence model is used as the turbulent closure for the primary phase and a zero equation turbulence model is used for the micro-bubbles. Change in the mean streamwise velocity profiles, void fraction, turbulence modification and other results are presented and discussed with corresponding change in the gas injection rates. The complex mechanism of drag reduction are scrutinised and explained in context to our numerical findings. Special attentions have been also devoted to divulge the effect of bubble coalescence and break-up caused by random collision and turbulent impact. Numerical results showed good agreement for the skin-friction coefficients against experimental data throughout various air injection rates. The MUSIG model was found to be one of the best candidates to resolve the bubble dynamics in micro-bubble-induced drag reduction problems.     

Robust design of microbubble drag reduction in a channel flow using the Taguchi method

This study attempts to obtain optimum parametric levels for robust design of the microbubble drag reduction in a turbulent channel flow. This work was carried out experimentally by measuring the frictional resistance on the upper wall of the channel to analyze the efficiency of drag reduction. Considering the mean flow speed as an indicative factor, several controllable factors that influence the effect of microbubble drag reduction were investigated in this work by using the Taguchi method. The controllable factors in this study were the amount of air injected, area of air injection, and microbubble size. For the condition of optimum parametric levels, the effect of drag reduced could reach up to 21.6%.     

潜航器微气泡减阻数值模拟方法研究

近年来微气泡减阻技术应用于水面舰船的相关研究取得重要进展,许多舰船已经采用了此类技术。然而该技术在潜航器上的应用整体上仍处于理论分析和试验测试阶段。以回转体潜航器简化模型为研究对象,利用多孔介质来简化喷气小孔,并结合不可压缩水体和可压缩理想气体的流体体积（VOF）方法,建立了Realizable ε-k湍流计算模式。通过拖曳试验验证了多孔介质等效喷气小孔的合理性和数值模式的准确性,结合试验和数值结果探究了微气泡减阻技术对潜航器航行阻力的影响。数值结果显示,气泡对潜航器尾流低速区的改变使尾部压力分布产生变化从而导致压差阻力增高。同时气泡可以显著降低其覆盖区域的黏性阻力。并且随着来流速度的提高,气泡覆盖范围扩大,黏性减阻率持续增加。进一步地,建立了加长改进模型,数值模拟结果表明,微气泡减阻技术不仅能大幅减小黏性阻力,还能有效减小模型总阻力。     

Ship hull characteristics from surface wake synthetic aperture radar (SAR) imagery

Improved data collection and processing technologies along with the use of high resolution spectral techniques soon will make it possible to obtain estimates of the Kelvin wave amplitude function A(θ), ship speed U, and ship heading α from synthetic aperture radar images of ship wakes. This paper presents a series of methods for deriving additional hull characteristics such as the length L, volume V, and offsets ζ(x,z) from this spectral and surface wave information. The first method estimates the ship length by taking the Fourier transform of the slope amplitude function |kA|. The remaining estimates make use of the hull inversion code developed at the University of Michigan by Wu in 1991. The accuracy of the hull offsets predicted by the code is first determined for various options for solving the linear inversion problem. In this case, both the magnitude and phase of A(θ) are known in addition to the hull draft H. Since the draft is not often known a priori, the accuracy of the code is determined next by predicting the volume of the ship for an approximate though plausible input value of H. Finally, the accuracy of the non-linear inversion problem of obtaining offsets along the entire hull is investigated when only the magnitude but not the phase of A(θ) is known.     

The effects of drag reduction by ribbons attached to cylindrical pipes

This study proposes a drag reduction device that uses three ribbons attached 120 degree apart to vertical pipes. Experiments were conducted in a circulating water channel to investigate the effects of the ribbon length and the direction of the flow on various current velocities. Drag on a vertical cylinder was measured by a resistance dynamometer. Flow visualizations were conducted using laser sheet beams. Laser Doppler Velocimetry (LDV) was used to measure the velocity field in the wake. This experiment demonstrates that attached ribbons can be used to reduce the drag force on vertical pipes for various directions of incoming flows. The ribbon-type device is very simple and easy to fabricate for field applications. The results are promising for the application to offshore structures.     

Numerical investigation of wave elevation and bottom pressure generated by a planing hull in finite-depth water

A numerical investigation of the bottom pressure and wave elevation generated by a planing hull in finite-depth water is presented. While the existing literature addresses the free-surface deformation and pressure field at the seafloor independently, this work proposes a direct comparison between the two hydrodynamic quantities. The dependence of the pressure disturbances at the ocean floor from the waves generated at the free-surface by a planing hull is studied for several values of both the depth and hull Froude numbers. The methodology employed is Smoothed Particle Hydrodynamics (SPH), a numerical technique based on the discretization of the continuum fields of hydrodynamics through mesh-less particles. The SPH code herein chosen is initially validated against experimental data for transom-stern flow. Subsequently, numerical simulations are presented for a planing hull in high-speed regimes. The results show a direct correlation between surface wave dynamics and hydrodynamic pressure disturbances at the seafloor as the value of the Froude number is varied. This is assessed by studying the inverse dependence of the low-pressure wake angle with the Froude number and by comparison of SPH results with similar works in the cited literature.     

Investigation of drag crisis phenomenon using CFD methods

When fluid flow passes a cylinder, the drag crisis phenomenon occurs between the sub-critical and the super-critical Reynolds numbers. The focus of the present studies was on the numerical prediction of the drag crisis based on CFD methods. In this work, block structured meshes with refined grids near the cylinder surface and in the downstream were employed. Both 2D and 3D simulations were performed using various turbulence models, including the SST k − ω model, the k − ϵ model, the SST with LCTM, the DES model, and the LES model. In the convergence studies, the effects of the grid size, the time step, the first grid size and the aspect ratio (for 3D simulations) on the solutions were examined. The errors due to spatial and time discretizations were quantified according to a V&V procedure. Validation studies were carried out for various Reynolds numbers between Re = 6.31 × 10⁴ and 7.57 × 10⁵. The averaged drag force, the RMS of lift force and the Strouhal number were compared with experimental data. The studies indicated that standard 2D and 3D RANS methods were inadequate to capture the drag crisis phenomenon. The LES method however has the potential to address the problem.     

An experimental study of interceptors for drag reduction on high-performance sailing yachts

Interceptors have been widely used in recent years in fast ferries and small high-speed leisure and commercial craft for ride and trim control, and steering. In the context of high-performance sailing yachts, they first appeared in 2008 on the yacht Ecover 3 which was dismasted while leading the Vendee Globe Challenge race. However, in spite of their popularity in power craft, few studies have been published investigating the impact of interceptors on vessel performance, and apparently none in the case of sailing yachts. In the current study, interceptors are compared with an aerodynamic device known as a Gurney flap. It is shown that interceptors are generally substantially smaller than Gurney flaps. A comprehensive experiment programme is presented exploring the impact of interceptors on the performance of an Open 60 yacht hull. Results show a marked reduction in calm-water resistance over a wide speed range, with benefits of 10–18% in the speed range between 8 and 20 knots, accompanied by reduced sinkage and trim. The gains observed are much larger than those observed in powercraft, and also substantially greater than those achievable through trim changes by moving ballast longitudinally. The benefits appear to be largely sustained in small waves.     

Analysis of interaction between ship bottom air cavity and boundary layer

The successful designs of hulls for ships employing drag reduction by air bottom cavitation have been based on solutions of inverse problems of the theory of ideal incompressible fluid. However, prediction of the drag reduction ratio, the air demand by ventilated cavities and the cavity impact on the hull–propeller interaction is impossible in the framework of this theory because all mentioned characteristics depend on interaction of air cavities with the ship boundary layers. Because the known CFD tools are not fitted to ventilated cavitation at low Froude numbers, an analysis of this interaction requires a novel flow model. This model includes the incompressible air flow in the ventilated cavity, the compressible flow of a water–air mixture in the boundary layer on cavities and downstream of them and the curl-free incompressible outer water flow. The provided 2D computations employing this model allows for explanations of the earlier observed effects and for prediction of the air demand by ventilated cavities. The computed velocity profiles downstream of cavities are in the accordance with the available experimental data.     

Tsunami detectability using open-ocean bottom pressure fluctuations

Rationale for the measurement of open-ocean tsunami signatures are presented, and available pertinent data are reviewed. Models for tsunami signature and background noise are proposed in order to synthesize an optimum tsunami receiver. Using these models, the minimum tsunami amplitude (in cm) to yield the probability of correct tsunami detectionP_{D} = 0.999and probability of false alarmP_{F} = 10^{-3}is found to be0.718/sqrt{f_{0}}, wheref_{0}is the tsunami dominant frequency (in cycles/h). A realizable receiver is proposed and its performance is evaluated using actual tsunami signatures. It is demonstrated that the detection of a tsunami with an average amplitude as small as 0.7 cm is possible for theP_{D}andP_{F}as above. Simulation results using synthesized background noise are shown. Tidal effects on the receiver performance also are considered and are found negligible for a certain range of the receiver parameters, resulting in a considerable reduction of the signal processing required.     

Model test of the surface and submerged vehicles with the micro-bubble drag reduction

This study is based on the effective experiment observation and measuring technology to discuss the interaction influence between liquid turbulent boundary layer and a crowded group micro-bubbles. It is in order to understand and quantify the micro-bubbles clouds inside the turbulent boundary layer to eliminate the capacity of skin friction drag. Whenever the micro-bubbles are over supplied, pile up effect happened which makes micro-bubbles to integrate to each other as a large-size air film. Although they still have the drag reduction effect, the efficiency of drag reduction slowed down at this transition period. In the experiment of vertical type circulating water tunnel, when 1 μm porous medium is at 7 m/s flow speed, the C_v value at 0.056 has the best drag reduction efficiency of 26%. While 10 μm porous medium is at the same flow speed, the drag reduction efficiency is only around 23%.     

Numerical research on drag reduction by ventilated partial cavity based on two-fluid model

This paper presents a numerical study on the drag reduction mechanism created by a ventilated partial cavity and its associated effects by the downstream dispersed microbubbles. A semi-empirical approach is introduced to model the discrete interface of the ventilated cavity and its complex gas leakage rate induced by the local turbulent shear stress. Based on the Eulerian–Eulerian two-fluid modeling framework, a population balance approach based on MUltiple-SIze-Group (MUSIG) model is incorporated to simulate the dynamical effects of bubbly flow along the test body. Particular attention is also directed to grasp a better understanding of the size evolution of microbubble and its associated effects on drag reduction. Model predictions are validated against three experimental measurements carried out in a high-speed water tunnel by Schauer (2003) and Wosnik et al. (2005). Close examination of the flow structures, gas void fraction distributions and its resultant density ratio provides valuable insights on the complex physical phenomenon, helping to consolidate idea to maximize the drag reduction for ventilated cavitating vehicles.     

定量模拟分析船体与空气枪间的安全距离

海洋地震勘探过程中,水中放置的空气枪与船体之间须保持一定安全距离,否则引起船载设施和仪器剧烈抖动,影响船舶设施和仪器的安全运行,甚至更严重情况下,有可能造成船体的损坏.利用海水压力脉冲与船体受到的冲击压力关系,结合建立的空气枪震源子波模型模拟近场震源子波,可计算出船体与空气枪间的安全距离.研究结果表明,根据震源子波模型模拟的近场子波主峰值,总容量150in3G枪、300 in3相干G枪及760 in3相干G枪的安全距离分别为3.90、6.57 m和9.13 m.由此可知,空气枪与船体间距离只有大于安全距离,才能避免船体大的振动或设备免遭损害,保证海洋地震勘探有效、安全的进行.     

SWAN模型中不同风拖曳力系数对风浪模拟的影响

本文以荷兰哈灵水道海域为实验区域,通过敏感性实验,研究了在14 m/s、31.5 m/s和50 m/s（分别代表一般大风、强热带风暴和强台风的极端条件）定常风速下SWAN模型中不同风拖曳力系数对风浪模拟的影响程度。结果表明,对于近岸浅水区域（水深小于20 m）,风拖曳力系数计算方案的选择对有效波高影响较小,而且当风速增加到一定程度后,波浪破碎成为影响波高值的主要因素;对于深水区域（水深大于30 m）,一般大风条件下风拖曳力系数计算方案的选择对有效波高影响仍然较小,随着风速的继续增大,风拖曳力系数计算方案的选择对有效波高的影响逐渐显著。对于平均周期,风拖曳力系数计算方案的选择和风速的改变对其影响均较小,而由水深变浅导致的波浪破碎对其影响较为显著。根据敏感性实验结果,本文对SWAN模型中风拖曳力系数计算方案的选择做出如下建议:计算近岸浅水区域风浪场或深水区域一般大风条件风浪场时,其风拖曳力系数可以直接采用模型默认选项;而对于深水区域更大风速条件,可首先采用模型默认选项试算,然后结合当地海域实测波浪资料进行修正。     

A metrologic method of anomaly field amplitude bottom reduction in undersampled geomagnetic marine surveys

We show the results obtained by means of a seabed reduction technique on the intensity of geomagnetic anomaly fields applied to a synthetic case and then to the real case of a geomagnetic survey of eastern Ligurian Sea (Italy). The eastern Ligurian Sea has very intense short waves anomaly fields and a sea bed that varies greatly in depth. As a result the geomagnetic space signal is characterized by a very large spectral content; in these conditions it is not possible to obtain a full sampled marine survey and vertical continuation analytic procedures and standard numerical bottom reduction based on a single vertical incremental parameter, whichever is applicable, fails to give accurate results. The present technique, which has been fine-tuned over 4 years of experimentation in environmental researchs, aims to provide a simple and efficient means to reduce the distortion of geomagnetic anomalies field caused by the variation of distance between survey plane and magnetic outcrop source position. The compensation procedure is based on evaluation, by comparison of two measurements carried out at different altitudes, of the mean vertical increment typical of each anomaly field principal frequency component bands. The component anomaly fields are then corrected by application of the corresponding vertical increments and lastly, the anomaly geomagnetic field reduced to the sea-bed is computed as Inverse Fourier Transform of a spectrum built as synthesis of the component anomaly fields' spectra. The results obtained have shown a notable increase in definition of anomaly field intensity without the production of appreciable distortions or false geomagnetic echoes. This revised version was published online in November 2006 with corrections to the Cover Date.     

Turbulent free-surface flow around a Wigley hull using the slightly compressible flow formulation

The slightly compressible flow formulation is applied to the free-surface, three-dimensional turbulent flow around a Wigley hull. Two turbulence models (large eddy simulation and Baldwin–Lomax) are used and compared. The simulation conditions are the ones for which experimental and numerical results exist. The computational grid is built using an algebraic grid generator with the model fixed in space. The codes use the interface-capturing technique for computing the free-surface displacements and the Beam and Warming scheme for marching in time the numerical model. The results compare well with the experimental data available.     

Reflection profiling in the deep ocean using a near bottom hydrophone

A 20 km long high resolution seismic reflection profile was carried out approximately 300 km southwest of Bermuda. The data were collected using a small airgun sound source and a single hydrophone receiver towed 100 m above the seafloor at a depth of 5400 m. Comparisons with nearby conventional seismic reflection profiles show the considerable improvement of resolving power provided by this method, particularly of the basement morphology beneath the 700 m thick sediment column. The data were recorded digitally and selected data examples show the enhancement provided by filtering, stacking, source deconvolution and corrections for hydrophone motion. The precise picture of basement topography that results from this data is compared with deep tow bathymetry profiles of the present day mid-Atlantic Ridge flanks, and is seen to be remarkably similar.     

An inverse method for underwater bottom topography by using SAR imagery

I~IOXThe image technology of the synthetic aperture radar (SAR) is one of most important advances in spaceberne microwave remote sensing (Larson et al., 1976). It has extensive application in the ocean remote sensing, such as inversion of the sea subdue winds and currentS, reconstruction of underwater hatom topography and so on. Owing tO the high conductivity of sea water, an electromagnetic wave actually cannot penetrate into the sea water, with extremely smallPenetration depth (abbot 1 c…     

用合成孔径雷达图像反演浅海水下地形的一种方法

海流与海底地形的相互作用导致了海表面的粗糙起伏,从而引起对海面观测的雷达散射截面的变化。在浅海海流速度垂直分布廓线为均匀的假定下,由流体连续性方程和驰豫时间近似下流体力学弱相互作用的理论,已证明雷达散射截面的空间变化正比于海流方向上海流速度变化的梯度。研究了合成孔径雷达(SAR)对海面的观测图像在无槽道信息和海流方向预知条件下,提出用散射系数的二维相关函数确定海流方向。在海流速度和方向已知条件下,推导了由雷达散射截面的空间变化迭推反演浅海水下地形的公式,并研究了航天飞机SIR-CSAR在我香港特区海域的二维图像反演浅海水下地形的应用。