High Resolution Ship Hydrodynamics Simulations in Open Source Environment

The numerical simulation of wake and flee-surface flow around ships is a complex topic that involves multiple tasks： the generation of an optimal computational grid and the development of numerical algorithms capable to predict the flow field around a hull. In this paper, a numerical framework is developed aimed at high-resolution CFD simulations of turbulent, free-surface flows around ship hulls. The framework consists in the concatenation of ＂tools＂, partly available in the open-source finite volume library OpenFOAM. A novel, flexible mesh-generation algorithm is presented, capable of producing high-quality computational grids for free-surface ship hydrodynamics. The numerical frame work is used to solve some benchmark problems, providing results that are in excellent agreement with the experimental measures.     

CFD技术在船舶与海洋工程复杂粘性流动中的应用进展（英文）

Complex flow around floating structures is a highly nonlinear problem, and it is a typical feature in ship and ocean engineering.Traditional experimental methods and potential flow theory have limitations in predicting complex viscous flows. With the improvement of high-performance computing and the development of numerical techniques, computational fluid dynamics(CFD) has become increasingly powerful in predicting the complex viscous flow around floating structures. This paper reviews the recent progress in CFD techniques for numerical solutions of typical complex viscous flows in ship and ocean engineering.Applications to free-surface flows, breaking bow waves of high-speed ship, ship hull–propeller–rudder interaction, vortexinduced vibration of risers, vortex-induced motions of deep-draft platforms, and floating offshore wind turbines are discussed.Typical techniques, including volume of fluid for sharp interface, dynamic overset grid, detached eddy simulation, and fluid–structure coupling, are reviewed along with their applications. Some novel techniques, such as high-efficiency Cartesian grid method and GPU acceleration technique, are discussed in the last part as the future perspective for further enhancement of accuracy and efficiency for CFD simulations of complex flow in ship and ocean engineering.     

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

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.     

船体周围Stokes波浪流耦合VOF模型的数值模拟（英文）

The surface wave generated by flow around a ship hull moving near free surface of water is simulated numerically in this study. The three-dimensional implicit finite volume method(FVM) is applied to solve Reynolds averaged Navier-Stokes(RANS) equation. The realizable k-ε turbulence model has been implemented to capture turbulent flow around the ship hull in the free surface zone. The volume of fluid(VOF) method coupled with the Stokes wave theory has been used to determine the free surface effect of water. By using is a six degrees of freedom model, the ship hull's movement is numerically solved with the Stokes wave together. Under the action of Stokes waves on the sea, the interface between the air and water waves at the same regular pattern and so does the pressure and the vertical velocity. The ship hull moves in the same way as the wave. The amplitude of the ship hull's heave is less than the wave height because of the viscosity damping. This method could provide an important reference for the study of ships' movement, wave and hydrodynamics.     

Numerical Simulation of the Stokes Wave for the Flow around a Ship Hull Coupled with the VOF Model

The surface wave generated by flow around a ship hull moving near free surface of water is simulated numerically in this study. The three-dimensional implicit finite volume method(FVM) is applied to solve Reynolds averaged Navier-Stokes(RANS) equation. The realizable k-ε turbulence model has been implemented to capture turbulent flow around the ship hull in the free surface zone. The volume of fluid(VOF) method coupled with the Stokes wave theory has been used to determine the free surface effect of water. By using is a six degrees of freedom model, the ship hull’s movement is numerically solved with the Stokes wave together. Under the action of Stokes waves on the sea, the interface between the air and water waves at the same regular pattern and so does the pressure and the vertical velocity. The ship hull moves in the same way as the wave. The amplitude of the ship hull’s heave is less than the wave height because of the viscosity damping. This method could provide an important reference for the study of ships’ movement, wave and hydrodynamics.     

A porous cell method for the simulation of fluid-solid interactions

There are many ways of describing a solid,porous or fluid region of the computational domain when solving the Navier-Stokes equations(NSE)for flow motions.Amongst these the porous cell method is one of the most flexible approaches.In this method,a parameter is defined as a ratio of the volume open to water and air in a calculation cell to its cell volume.In the calculation,the same numerical procedure is applied to every cell and no explicit boundary conditions are needed at solid boundaries.The method is used to simulate flow through porous media,around solid bodies and over a moving seabed.The results compare well with experimental data and other numerical results.In our future work the porous cell method will be applied to more complex fluid-solid interaction situations.     

驱逐舰烟气排放过程数值仿真及同模型试验比较研究(英文)

The exhaust smoke dispersion for a generic frigate is investigated numerically through the numerical solution of the governing fluid flow, energy, species and turbulence equations. The main objective of this work is to obtain the effects of the yaw angle, velocity ratio and buoyancy on the dispersion of the exhaust smoke. The numerical method is based on the fully conserved control-volume representation of the fully elliptic Navier-Stokes equations. Turbulence is modeled using a two-equation(k-ε) model. The flow visualization tests using a 1/100 scale model of the frigate in the wind tunnel were also carried out to determine the exhaust plume path and to validate the computational results. The results show that down wash phenomena occurs for the yaw angles between ψ =10° and 20°. The results with different exhaust gas temperatures show that the buoyancy effect increases with the increasing of the exhaust gas temperature. However, its effect on the plume rise is less significant in comparison with its momentum. A good agreement between the predictions and experiment results is obtained.     

Risk zone of wrack hitting marine structure simulated by 2D hydraulic model

The wrack or the ship out of control will drift with flow.One of the most important factors that drive the ship is flow current which moves circularly in tidal area.The wrack from same place always drifts in different ways if the start time is different.So,during the ship drifting period,the drift trace is also determined by both wave and wind forces.The drift direction is limited by water depth which must be deeper than ship draft. These marine structures that can not afford the hit of wrack or will destroy the wrack must be well considered when they are placed near harbor and waterway or other water area with ship running.The risk zone should be consulted according to tide and weather conditions to protect structures and ships in necessary.A method is presented here to simulate the risk zone by 2D numerical hydraulic model with tidal current,wave,wind and water depth considered.This model can be used to built early-warning and protect system for special marine structure.     

海床对海洋管道涡旋脱落和受力系数的影响（英文）

The effect of rigid bed proximity on flow parameters and hydrodynamic loads in offshore pipelines exposed to turbulent flow is investigated numerically. The Galerkin finite volume method is employed to solve the unsteady incompressible 2 D Navier–Stokes equations. The large eddy simulation turbulence model is solved using the artificial compressibility method and dual time-stepping approach. The proposed algorithm is developed for a wide range of turbulent flows with Reynolds numbers of 9500 to 1.5×104. Evaluation of the developed numerical model shows that the proposed technique is capable of properly predicting hydrodynamic forces and simulating the flow pattern. The obtained results show that the lift and drag coefficients are strongly affected by the gap ratio. The mean drag coefficient slightly increases as the gap ratio increases, although the mean lift coefficient rapidly decreases. The vortex shedding suppression happen at the gap ratio of less than 0.2.     

Numerical modeling of exhaust smoke dispersion for a generic frigate and comparisons with experiments

The exhaust smoke dispersion for a generic frigate is investigated numerically through the numerical solution of the governing fluid flow, energy, species and turbulence equations. The main objective of this work is to obtain the effects of the yaw angle velocity ratio and buoyancy on the dispersion of the exhaust smoke The numerical method is based on the fully conserved control-volume representation of the fully elliptic Navier-Stokes equations. Turbulence is modeled using a two-equation （k-ε） model The flow visualization tests using a 1/100 scale model of the frigate in the wind tunnel were also carried out to determine the exhaust plume path and to validatc the computational results. The results show that down wash phenomena occurs for the yaw angles between ψ=10° and 20°. The results with different exhaust gas temperatures show that the buoyancy effect increases with the increasing of the exhaust gas temperature. However, its effect on the plume rise is less significant in comparison with its momentum. A good agreement between the predictions and experiment results is obtained.     

粘流方法的螺旋桨尾流场数值分析(英文)

The computational fluid dynamics(CFD) method is used to numerically simulate a propeller wake flow field in open water.A sub-domain hybrid mesh method was adopted in this paper.The computation domain was separated into two sub-domains,in which tetrahedral elements were used in the inner domain to match the complicated geometry of the propeller,while hexahedral elements were used in the outer domain.The mesh was locally refined on the propeller surface and near the wake flow field,and a size function was used to control the growth rate of the grid.Sections at different axial location were used to study the spatial evolution of the propeller wake in the region ranging from the disc to one propeller diameter(D) downstream.The numerical results show that the axial velocity fluctuates along the wake flow;radial velocity,which is closely related to vortices,attenuates strongly.The trailing vortices interact with the tip vortex at the blades’ trailing edge and then separate.The strength of the vortex shrinks rapidly,and the radius decreases 20% at one diameter downstream.     

周期性阵风流对通气超空泡尾流结构影响的数值研究（英文）

A computational model is established to investigate the effects of a periodic gust flow on the wake structure of ventilated supercavities. The effectiveness of the computational model is validated by comparing with available experimental data.Benefited from this numerical model, the vertical velocity characteristics in the entire flow field can be easily monitored and analyzed under the action of a gust generator; further, the unsteady evolution of the flow parameters of the closed region of the supercavity can be captured in any location. To avoid the adverse effects of mounting struts in the experiments and to obtain more realistic results, the wake structure of a ventilated supercavity without mounting struts is investigated. Unsteady changes in the wake morphology and vorticity distribution pattern of the ventilated supercavity are determined. The results demonstrate that the periodic swing of the gust generator can generate a gust flow and, therefore, generate a periodic variation of the ventilated cavitation number σ. At the peak σ, a re-entrant jet closure appears in the wake of the ventilated supercavity. At the valley σ, a twin-vortex closure appears in the wake of the ventilated supercavity. For the forward facing model, the twin vortex appears as a pair of centrally rolled-up vortices, due to the closure of vortex is affected by the structure. For the backward facing model,however, the twin vortex appears alternately as a pair of centrally rolled-up vortices and a pair of centrally rolled-down vortices,against the periodic gust flow.     

能量耗散效应的多域边界元法(英文)

The wave diffraction and radiation around a floating body is considered within the framework of the linear potential theory in a fairly perfect fluid.The fluid domain extended infinitely in the horizontal directions but is limited by the sea bed,the body hull,and the part of the free surface excluding the body waterplane,and is subdivided into two subdomains according to the body geometry.The two subdomains are connected by a control surface in fluid.In each subdomain,the velocity potential is described by using the usual boundary integral representation involving Green functions.The boundary integral equations are then established by satisfying the boundary conditions and the continuous condition of the potential and the normal derivation across the control surface.This multi-domain boundary element method(MDBEM) is particularly interesting for bodies with a hull form including moonpools to which the usual BEM presents singularities and slow convergence of numerical results.The application of the MDBEM to study the resonant motion of a water column in moonpools shows that the MDBEM provides an efficient and reliable prediction method.     

Analytical study of magnetohydrodynamic propulsion stability

In this paper an analytical solution for the stability of the fully developed flow drive in a magneto-hydro-dynamic pump with pulsating transverse Eletro-magnetic fields is presented. To do this, a theoretical model of the flow is developed and the analytical results are obtained for both the cylindrical and Cartesian configurations that are proper to use in the propulsion of marine vessels. The governing parabolic momentum PDEs are transformed into an ordinary differential equation using approximate velocity distribution. The numerical results are obtained and asymptotic analyses are built to discover the mathematical behavior of the solutions. The maximum velocity in a magneto-hydro-dynamic pump versus time for various values of the Stuart number, electro-magnetic interaction number, Reynolds number, aspect ratio, as well as the magnetic and electrical angular frequency and the shift of the phase angle is presented. Results show that for a high Stuart number there is a frequency limit for stability of the fluid flow in a certain direction of the flow. This stability frequency is dependent on the geometric parameters of a channel.     

多工况涡流管内部分离过程的三维数值研究（英文）

Air separators provide safe, clean, and appropriate air flow to engines and are widely used in vehicles with large engines such as ships and submarines. In this operational study, the separation process in a Ranque–Hilsch vortex tube cleaning(cooling) system is investigated to analyze the impact of the operating gas type on the vortex tube performance; the operating gases used are air, nitrogen, oxygen, carbon dioxide and nitrogen dioxide. The computational fluid dynamic model used is equipped with a three-dimensional structure, and the steady-state condition is applied during computations. The standard k–ε turbulence model is employed to resolve nonlinear flow equations, and various key parameters, such as hot and cold exhaust thermal drops, and power separation rates, are described numerically. The results show that nitrogen dioxide creates the greatest separation power out of all gases tested, and the numerical results are validated by good agreement with available experimental data. In addition, a comparison is made between the use of two different boundary conditions, the pressure-far-field and the pressure-outlet, when analyzing complex turbulent flows in the air separators. Results present a comprehensive and practical solution for use in future numerical studies.     

船首形状对自由液面水下航行体阻力影响的试验研究（英文）

In this paper, towing tank experiments are conducted to study the behavior of flow on a model of the underwater vehicle with various shapes of bows, i.e. tango and standard bows in free surface motion tests. The total resistances for different Froude numbers are considered experimentally. The towing tank is equipped with a trolley that can operate in through 0.05–6 m/s speed with ±0.02 m/s accuracy. Furthermore, the study is done on hydrodynamic coefficients i.e. total, residual and friction resistance coefficients, and the results are compared. Finally, the study on flow of wave fields around bows is done and wave filed around two bows are compared. The Froude number interval is between 0.099 and 0.349. Blockage fraction for the model is fixed to 0.005 3. The results showed that the residual resistance of the standard bow in 0.19 to 0.3 Froude number is more than the tango bow in surface motion which causes more total resistance for the submarine. Finally, details of wave generated by the bow are depicted and the effects of flow pattern on resistance drag are discussed.     

A free surface frequency domain green function with viscous dissipation and partial reflections from side walls

The free-surface Green function method is widely used in solving the radiation or diffraction problems caused by a ship or ocean structure oscillating on the waves. In the context of inviscid potential flow, hydrodynamic problems such as multi-body interaction and tank side wall effect cannot be properly dealt with based on the traditional free-surface frequency domain Green function method, in which the water viscosity is omitted and the energy dissipation effect is absent. In this paper, an open-sea Green function with viscous dissipation was presented within the theory of visco-potential flow. Then the tank Green function with a partial reflection from the side walls in wave tanks was formulated as a formal sum of open-sea Green functions representing the infinite images between two parallel side walls of the source in the tank. The new far-field characteristics of the tank Green function is vitally important for improving the validity of side-wall effects evaluation, which can be used in supervising the tank model tests.     

长椭球体绕流问题中层流向湍流流动状态转变的数值预测（英文）

In this work, the laminar-to-turbulent transition phenomenon around the two-and three-dimensional ellipsoid at different Reynolds numbers is numerically investigated. In the present paper, Reynolds Averaged Navier Stokes(RANS) equations with the Spalart-Allmaras, SST k-ω, and SST-Trans models are used for numerical simulations. The possibility of laminar-toturbulent boundary layer transition is summarized in phase diagrams in terms of skin friction coefficient and Reynolds number.The numerical results show that SST-Trans method can detect different aspects of flow such as adverse pressure gradient and laminar-to-turbulent transition onset. Our numerical results indicate that the laminar-to-turbulent transition location on the 6:1 prolate spheroid is in a good agreement with the experimental data at high Reynolds numbers.     

基于流固耦合方法的两自由度水翼流激振动特性研究（英文）

The flutter of a hydrofoil can cause structural damage and failure, which is a dangerous situation that must be avoided. In this work, based on computational fluid dynamics and structural finite element methods, a co-simulation framework for the flow-induced vibration of hydrofoil was established to realize fluid-structure interaction. Numerical simulation research was conducted on the flow-induced vibration characteristics of rigid hydrofoil with 2-DOF under uniform flow, and the heave and pitc...     

Numerical analysis of transient fluid-structure interaction of warship impact damage caused by underwater explosion using the FSLAB北大核心CSCD

［Objectives］This paper aims to address the numerical simulation problems of the dynamic response of ships subject to near-, medium- and far-field underwater explosions by establishing several numerical methods and calculation models. ［Methods］First, load and fluid-structure interaction models are established on the basis of the Eulerian finite element method and acoustic finite element method using the field-split technique, and FSLAB fluid-structure interaction software is developed. Next, near-, medium- and far-field underwater explosions are numerically simulated respectively. The shock wave propagation law, bubble shape and load evolution characteristics of near free-surface and near-wall underwater explosions are obtained, and the shock response characteristics of a spherical shell and ship subject to far-field underwater explosions are analyzed. Finally, the FSLAB software results are compared with the analytical solutions, reference solutions and experimental data. ［Results］The results show that the FSLAB fluid-structure interaction software developed in this paper is effective and accurate in simulating the impact damage of underwater explosions on warships. ［Conclusion］This study can provide a basis and support for the power assessment of underwater anti-explosion and shock design of warships. © 2022 Journal of Clinical Hepatology. All rights reserved.