Viscosity and nonlinearity effects on the forces and waves generated by a floating twin hull under heave oscillation

Nonlinear hydrodynamics of a twin rectangular hull under heave oscillation is analyzed using numerical methods. Two-dimensional nonlinear time-domain solutions to both inviscid and viscous problems are obtained and the results are compared with linear, inviscid frequency-domain results obtained in [26] to quantify nonlinear and viscous effects. Finite-difference methods based on boundary-fitted coordinates are used for solving the governing equations in the time domain [2]. A primitive-variables based projection method [6] is used for the viscous analysis and a mixed Eulerian–Lagrangian formulation [11] for inviscid analysis. The algorithms are validated and the order of accuracy determined by comparing the results obtained from the present algorithm with the experimental results of Vugt [22] for a heaving rectangle in the free surface. The present study on the twin-hull hydrodynamics shows that at large and non-resonant regular frequencies, and small amplitude of body oscillation, the fluid viscosity does not significantly affect the wave motion and the radiation forces. At low frequencies however the viscosity effect is found to be significant even for small amplitude of body oscillation. In particular, the hydrodynamic force obtained from the nonlinear viscous analysis is found to be closer to the linear inviscid force than the nonlinear inviscid force to the linear inviscid force, the reason for which is attributed to the wave dampening effect of viscosity. Since the wave lengths generated at smaller frequencies of oscillation are longer and therefore the waves could have a more significant effect on the dynamic pressure on the bottom of the hulls which contribute to the heave force, the correlation between the heave force and the wave elevation is found to be larger at smaller frequencies. Because of nonlinearity, the wave radiation and wave damping force remained nonzero even at and around the resonant frequencies – with the resonant frequencies as determined in [26] using linear potential flow theory. As to be expected, the nonlinear effect on the wave force is found to be significant at all frequencies for large amplitude of oscillation compared to the hull draft. The effect of viscosity on the force, by flow separation, is also found to be significant for large amplitude of body oscillation.     

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.     

Numerical Simulation of Viscous Flow Around a Rolling Cylinder with Ship-Like Section

Two-dimensional unsteady incompressible viscous flow around a rolling cylinder with ship-like section is numerically simulated by employing the computational scheme previously developed by the authors, in which the continuity and momentum equations are satisfied simultaneously at each time step for oscillating flow. The numerical results show that the motion of vortices around a rolling ship hull is cyclical. It is found that the location of the vortices is very similar to the existing experimental result. Using these simulation results, we can calculate the roll damping of ships including viscous effects.     

Vortex Patterns and Roll Damping on Various Cross Sections of Ships

A numerical study is presented on roll damping of ships by solving Navier-Stokes equation.Two-Dimensional unsteady incompressible viscous flow around the rolling cylinders of various ship-likecross sections are numerically simulated by use of the computational scheme previously developed by theauthors.The numerical results show that the location of the vortices is very similar to the existing experi-mental result.For comparison of vortex patterns and roll damping on various ship-like cross sections,vari-ous distributions of shear stress and pressure on the rolling ship hull surface are presented in this paper.Itis found that there are two vortices around the midship-like section and there is one vortex around the foreor stern section.Based on these simulation results.the roll damping of a ship including viscous effects iscalculated.The contribution of pressure to the roll moment is larger than the contribution of frictionalshear stress.     

Propulsive performance of a tanker hull form in damaged conditions

Many disastrous oil spill accidents from damaged vessels become worse especially when the early treatment is not prompt enough. To properly handle this type of accidents and prevent further disasters, International Maritime Organization establishes and imposes various rules and regulations. Better understanding of the propulsive performance of damaged vessels is important for containing the oil spill while the vessels are being towed or self-propelled. In the present study, both experimental and computational methods were used to investigate the flow phenomena around the hull and the hydrodynamic performances of a VLCC in various damaged conditions. From the resistance and self-propulsion test results, it is found that higher power is required to propel the ship especially with the bow trim. Wake measurement data provide physical insight into the factors to be considered for the propeller operation in damaged conditions.     

Hydrodynamic shape optimization by high fidelity CFD solver and Gaussian process based response surface method

A computational framework for hydrodynamic shape optimization of complex ship hull form is proposed and applied to improve the calm water performance of the KRISO Container Ship (KCS). The framework relies on three key features: a novel shape morphing method based on a combination of subdivision surfaces and free form deformations, a robust three dimensional viscous computational fluid dynamic solver based on the openFOAM open-source libraries and a Gaussian process-response surface method (GP-RSM) based on ordinary Kriging model which has been created to speed-up the evaluation of the quantity of interest (QoI) of the design process.The accuracy of the hydrodynamic solver is proven by comparing the obtained results against available experimental measurements. A preliminary sensitivity analysis on the mesh size has been carried out aiming at reducing the computational burden required by the CFD predictions. Three GP-RSMs have been trained relying on increasing number of hull designs. Each surrogate model has been cross-validated by both leave-one-out and k-fold techniques. The behaviours of these multi-dimensional surfaces have been analyzed in details by sampling the investigated design space with 10⁷ points according to a Full-Factorial algorithm, highlighting the regions of maximum deviation with respect to the resistance of the reference hull. The three optimum designs provided by the corresponding GP-RSM models have been verified by using high-fidelity CFD simulations with a refined mesh configuration. Calm water resistance, wave patterns and pressure distributions over the selected hull surfaces have been discussed in the light of the generated shape variations.     

Practical evaluation of sinkage and trim effects on the drag of a common generic freely floating monohull ship

A practical method to account for the influence of sinkage and trim on the drag of a freely floating (free to sink and trim) common monohull ship at a Froude number F ≤ 0.45 is considered. The sinkage and the trim are estimated via two alternative simple methods, considered previously. The drag is also estimated in a simple way, based on the classical Froude decomposition into viscous and wave components. Specifically, well-known semiempirical expressions for the friction drag, the viscous pressure drag and the drag due to hull roughness are used, and the wave drag is evaluated via a practical linear potential flow method. This simple approach can be used for ship models as well as full-scale ships with smooth or rough hull surfaces, and is well suited for early ship design and optimization. The method considered here to determine the sinkage and the trim, and their influence on the drag, yields theoretical predictions of the drag of the Wigley, S60 and DTMB5415 hulls that are much closer to experimental measurements than the corresponding predictions for the hull surfaces of the ships in equilibrium position at rest. These numerical results suggest that sinkage and trim effects, significant at Froude numbers 0.25 < F, on the drag of a typical freely floating monohull ship can be realistically accounted for in a practical manner that only requires simple potential flow computations without iterative computations for a sequence of hull positions.     

万米级深海陶瓷耐压结构水下内爆流场数值模拟

耐压结构是深海潜器的重要组成部分,但在深海的高压环境中却存在内爆的风险。为研究陶瓷耐压结构水下内爆的流场特性,使用针对可压缩多相流问题开发的开源代码,采用直接数值模拟,应用自适应直角网格,对两种压力条件下的耐压结构水下内爆进行了数值模拟。通过低压模拟结果与理论解和试验值比较,验证了模拟方法的有效性,进而开展万米级深海陶瓷耐压结构水下内爆模拟。分析发现：陶瓷耐压结构发生内爆后,其内部气腔存在多次压缩—反弹现象,深海环境压力越大则反弹越不明显;气腔反弹阶段,在结构外部将产生数倍于深海环境压力的冲击波,且传播速度接近声速;冲击波压力峰值与到球心距离呈负指数幂函数关系;在相同深海环境压力下,耐压结构外部监测点的冲击波压力与球体半径呈正比例关系。     

Numerical calculation of free-surface potential flow around a ship using the modified Rankine source panel method

A modified Rankine source panel method is presented for solving a linearized free-surface flow problem with respect to the double-body potential. The method of solution is based on the distribution of Rankine sources on the hull as well as its image and on the free surface. An iterative algorithm is used for determining the free surface and wave resistance using upstream finite difference operator. A verification of numerical modeling is made using the Wigley hull and the validity of the computer program is examined by comparing the details of wave profiles and wave making resistance with Series 60 model.     

Study of air-ventilated cavity under model hull on water surface

The use of air cavities beneath ship hulls can lead to significant drag reduction. A study of air-ventilated cavities under a simplified hull has been undertaken. Experiments with a 56-cm-long stepped-hull model were carried in an open-surface water channel at flow velocities 28–86 cm/s. The air-cavity parameters were measured at different model positions. Different cavity forms, a strong growth of the cavity length with the flow velocity, and an optimal trim angle for the largest air-cavity area were identified. Numerical studies were conducted using a linear potential-flow method and the finite-volume viscous code Fluent. The computationally inexpensive three-dimensional potential-flow modeling predicted air-cavity shapes and provided qualitative agreement with the measured average length of the air cavity. Two-dimensional viscous modeling reasonably predicted macroscopic features and viscous effects in the air-cavity flow, while exaggerated the mixed-phase flow regions.     

Optimization Design of Minimum Total Resistance Hull Form Based on CFD Method

In order to reduce the resistance and improve the hydrodynamic performance of a ship, two hull form design methods are proposed based on the potential flow theory and viscous flow theory. The flow fields are meshed using body-fitted mesh and structured grids. The parameters of the hull modification function are the design variables. A three-dimensional modeling method is used to alter the geometry. The Non-Linear Programming (NLP) method is utilized to optimize a David Taylor Model Basin (DTMB) model 5415 ship under the constraints, including the displacement constraint. The optimization results show an effective reduction of the resistance. The two hull form design methods developed in this study can provide technical support and theoretical basis for designing green ships.     

Simplified CFD modeling for bilge keel force and hull pressure distribution on a rotating cylinder

A horizontal, circular cylinder fitted with one bilge keel is forced to rotate harmonically around its axis. The bilge keel load and hull pressure distribution are investigated. A fully submerged condition (infinite fluid), and three partly-submerged conditions are considered. A two-dimensional numerical study is performed, and the results are validated against recently published experimental data by van’t Veer et al. [30]. In addition, comparisons for mass and drag coefficients are also made with experimental data for plate in infinite fluid (Keulegan and Carpenter [8]), and wall-mounted plate (Sarpkaya and O’Keefe [9]) in oscillatory flow.A Navier–Stokes solver based on the Finite Volume Method is adopted for solving laminar flow of incompressible water. The free-surface condition is linearized by neglecting the nonlinear free-surface terms and the influence of viscous stresses in the free surface zone, while the body-boundary condition is exact. This simplified modeling of the problem required the mesh to be fine only around the bilge keels, leading to a total number of cells around N ≃ 1 ×10⁴, which reduced computational cost significantly.The influence of draft and amplitude of oscillations on the bilge keel force and hull pressure distribution are considered. The bilge keel force is presented in terms of non-dimensional drag and mass coefficients including higher harmonic components. The numerical results are also compared with the industry standard empirical method for calculation of roll damping proposed by Ikeda et al. [4]. In general, a good agreement between the results of the present numerical method and the experimental data is obtained and the differences with those predicted by the empirical method are addressed.     

Modelling of multi-bodies in close proximity under water waves—Fluid forces on floating bodies

This work presents two-dimensional numerical results of the dependence of wave forces of multiple floating bodies in close proximity on the incident wave frequency, gap width, body draft, body breadth and body number based on both viscous fluid and potential flow models. The numerical models were validated by the available experimental data of fluid oscillation in narrow gaps. Numerical investigations show that the large amplitude responses of horizontal and vertical wave forces appear around the fluid resonant frequencies. The convectional potential flow model is observed to un-physically overestimate the magnitudes of wave forces as the fluid resonance takes place. By introducing artificial damping term with appropriate damping coefficients μ∈[0.4, 0.5], the potential flow model may work as well as the viscous fluid model, which agree with the damping coefficients used in our previous work for the predication of wave height under gap resonance. In addition, the numerical results of viscous fluid model suggest that the horizontal wave force is highly dependent on the water level difference between the opposite sides of an individual body and the overall horizontal wave force on the floating system is generally smaller than the summation of wave force on each body.     

Numerical analysis of control surface effects on AUV manoeuvrability

Computational fluid dynamics, CFD, is becoming an essential tool in the prediction of the hydrodynamic efforts and flow characteristics of underwater vehicles for manoeuvring studies. However, when applied to the manoeuvrability of autonomous underwater vehicles, AUVs, most studies have focused on the determination of static coefficients without considering the effects of the vehicle control surface deflection. This paper analyses the hydrodynamic efforts generated on an AUV considering the combined effects of the control surface deflection and the angle of attack using CFD software based on the Reynolds-averaged Navier–Stokes formulations. The CFD simulations are also independently conducted for the AUV bare hull and control surface to better identify their individual and interference efforts and to validate the simulations by comparing the experimental results obtained in a towing tank. Several simulations of the bare hull case were conducted to select the k–ω SST turbulent model with the viscosity approach that best predicts its hydrodynamic efforts. Mesh sensitivity analyses were conducted for all simulations. For the flow around the control surfaces, the CFD results were analysed according to two different methodologies, standard and nonlinear. The nonlinear regression methodology provides better results than the standard methodology does for predicting the stall at the control surface. The flow simulations have shown that the occurrence of the control surface stall depends on a linear relationship between the angle of attack and the control surface deflection. This type of information can be used in designing the vehicle's autopilot system.     

The effect on resistance of displacement hull forms due to fore-to-aft pressure re-distribution

The results of an experimental study of the resistance reduction from modification of the pressure field around displacement hulls creating fore-aft pressure re-distribution are presented. The fore-aft pressure re-distribution is achieved by hull surface porosity created by orifices through the hull which are connected to each other by fore-aft ducting. This affects the fore-aft pressure gradient when moving. Tests were completed with Model A, which has parabolic shaped waterlines, wall sides, a flat bottom and no parallel midbody. Resistance and surface pressure measurements were obtained with and without the pressure re-distribution orifices. The extent of pressure re-distribution was varied by using three ducting sizes. The resistance measurements and surface pressure data show that in some cases there is a large fore-aft surface pressure gradient and the fore-aft pressure re-distribution results in a significant reduction in the vessel resistance. The test results suggest that fore-aft pressure re-distribution may have benefit in reducing displacement hull and/or SWATH strut resistance.     

A comparative study of numerical simulations for fluid–structure interaction of liquid-filled tank during ship collision

Although International Maritime Organization (IMO) has taken many measures to minimize ship collisions, ships carrying liquid cargo sometimes do get struck by other vessels. The outflow of crude oil causes very serious consequences to the environment. In such cases it is necessary to analyze the response of structure of struck liquid cargo-filled tank to account for fluid–structure interaction accurately. In this paper, numerical simulation of collision between a container ship with double hull very large crude carrier (VLCC) is presented. Three different numerical simulation mothods were adopted to model fluid–structure interaction in liquid-filled cargo tank, namely arbitrary Lagrangian–Eulerian finite element method, Lagrangian finite element method and linear sloshing model. The numerical simulation results reveal that the fluid–structure interaction of liquid cargo-filled tank has a significant effect on the motion and structural response of the struck cargo tank. Compared with the calculation results of ALE FE method, the linear sloshing model underestimates the influence of fluid–structure interaction of liquid cargo tank while the Lagrangian–Eulerian finite element method may be considered as the practical method for engineering applications as it provided more reasonable results with a relatively low central processing unit (CPU) time.     

Source models of flow through and around screens and gauzes

The source panel method is used to investigate the singularity models of flow through and around screens that have been developed by various authors. These studies are extended to complete a framework of models in terms of the level of complexity, the nature of the flow in the wake and the form of pressure drop boundary condition. The models that best represent the limited experimental data that is available are identified and the need for more experimental data, in order to carry out a thorough assessment, is noted.     

Estimation of hydrodynamic coefficients for flow around cylinders in side-by-side arrangement with variation in separation gap

A finite-volume method (FVM) using a multi-block grid technique has been applied to investigate the hydrodynamic coefficients and to categorize the two-dimensional flow pattern around a pair of adjacent cylinders (known here after as double cylinders). Matching grids were used at the interface between blocks (cylinders), and the boundary conditions are imposed around the periphery to exchange physical values across the interface. To validate the developed numerical method, a case with laminar viscous flow around a single cylinder was carried out, the results showed good agreement compared to previously published data. Hydrodynamic coefficients, Strouhal number, and stagnation point change were thoroughly investigated as a function of the separation between the two cylinders. It is found that the changes of hydrodynamic coefficients due to variations in separation are well characterized by the five different vortex flow patterns, which can also be used to explain the flow and vortex pattern around the double cylinders.     

非均匀风影响下风力机三维气动粘性流场的数值模拟

基于面向对象的开源软件OpenFOAM,选择美国国家新能源实验室(NREL)Phase VI风力机为对象,对以往研究较少的非均匀来流风速作用下风力机三维气动粘性流场进行数值模拟。采用较为接近于真实情况的指数型风剖面,计算了轮毂处风速分别为5、10、15和25 m/s四种工况下的叶片表面压力分布、叶片的推力、尾涡等气动力数据,并与均匀来流风速下的风力机气动力学性能进行详细的对比,探讨非均匀风剖面对风力机流场结构和流动特性影响的物理现象和规律。     

Numerical wave tank computations of nonlinear motions of two-dimensional arbitrarily shaped free floating bodies

Time domain simulations of nonlinear motions of two-dimensional floating bodies in waves are presented. The so called `body exact' approach is adopted in a numerical wave tank. A new scheme for pressure evaluation on the wetted hull is developed and systematically used with good results in terms of accuracy and stability. Strongly flared geometries are successfully handled even at very large amplitude motions. The validation of the code is carried out according to the 20th ITTC Seakeeping Committee recommendations through internal checks of consistency and through comparisons with available experimental data.