Interaction of two three-dimensional explosion bubbles

The interaction of two underwater explosion bubbles was mathematically analyzed in this paper. Based on the assumption of potential flow, high-order curved elements were used to discretize the boundary integral equation and solve it. Assuming that gas inside the bubble follows the isentropic rule, the Euler-Lagrange method was used to trace the evolution of the bubble, and when calculating the singular integral, the singularity of the double-layer singular integral was eliminated by reconstructing a principal-value integral of double-layer potential so that a more precise result could be obtained. Elastic mesh technique (EMT) was also used when tracing the evolution of the bubble interface, and numerical smoothing wasn't needed. A comparison of calculations using this three-dimensional model with results of the Reyleigh-Plesset bubble model shows that the three-dimensional model and calculation method in this paper is practical. This three-dimensional model was applied to simulate the interaction of two bubbles under the action of gravity, and the dynamic characteristics of two bubbles near the surface was also analyzed. Bubbles influenced by surface effects and gravity present severe non-linearity. This paper provides a reference for research into the dynamics of multi-bubbles.     

船体周围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.     

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

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.     

水下爆炸空化效应作用下全船响应的并行计算(英文)

As well as shock wave and bubble pulse loading, cavitation also has very significant influences on the dynamic response of surface ships and other near-surface marine structures to underwater explosive loadings. In this paper, the acoustic-structure coupling method embedded in ABAQUS is adopted to do numerical analysis of underwater explosion considering cavitation. Both the shape of bulk cavitation region and local cavitation region are obtained, and they are in good agreement with analytical results. The duration of reloading is several times longer than that of a shock wave. In the end, both the single computation and parallel computation of the cavitation effect on the dynamic responses of a full-scale ship are presented, which proved that reloading caused by cavitation is non-ignorable. All these results are helpful in understanding underwater explosion cavitation effects.     

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

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.     

基于混合迭代法的二维水翼空泡性能分析(英文)

In order to study cavitation characteristics of a 2-D hydrofoil, the method that combines nonlinear cavitation model and mixed-iteration is used to predict and analyze the cavitation performance of hydrofoils. The cavitation elements are nonlinearly disposed based on the Green formula and perturbation potential panel method. At the same time, the method that combines cavity shape for fixed cavity length (CSCL) iteration and cavity shape for fixed cavitation number (CSCN) iteration is used to work out the thickness and length of hydrofoil cavitations. Through analysis of calculation results, it can be concluded that the jump of pressure and velocity potentially exist between cavitation end area and non-cavitations area on suction surface when cavitation occurs on hydrofoil. In certain angles of attack, the cavitation number has a negative impact on the length of cavitations. And under the same angle of attack and cavitation number, the bigger the thickness of the hydrofoil, the shorter the cavitations length.     

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

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.     

Numerical Investigation on Two-dimensional Boundary Layer Flow with Transition

As a basic problem in many engineering applications, transition from laminar to turbulence still remains a difficult problem in computational fluid dynamics （CFD）. A numerical study of one transitional flow in two-dimensional is conducted by Reynolds averaged numerical simulation （RANS） in this paper. Turbulence model plays a significant role in the complex flows＇ simulation, and four advanced turbulence models are evaluated. Numerical solution of frictional resistance coefficient is compared with the measured one in the transitional zone, which indicates that Wilcox （2006） k-ω model with correction is the best candidate. Comparisons of numerical and analytical solutions for dimensionless velocity show that averaged streamwise dimensionless velocity profiles correct the shape rapidly in transitional region. Furthermore, turbulence quantities such as turbulence kinetic energy, eddy viscosity, and Reynolds stress are also studied, which are helpful to learn the transition＇s behavior.     

二维边界层流动转捩过程的数值计算(英文)

As a basic problem in many engineering applications, transition from laminar to turbulence still remains a difficult problem in computational fluid dynamics(CFD). A numerical study of one transitional flow in two-dimensional is conducted by Reynolds averaged numerical simulation(RANS) in this paper. Turbulence model plays a significant role in the complex flows' simulation, and four advanced turbulence models are evaluated. Numerical solution of frictional resistance coefficient is compared with the measured one in the transitional zone, which indicates that Wilcox(2006) k-ω model with correction is the best candidate. Comparisons of numerical and analytical solutions for dimensionless velocity show that averaged streamwise dimensionless velocity profiles correct the shape rapidly in transitional region. Furthermore, turbulence quantities such as turbulence kinetic energy, eddy viscosity, and Reynolds stress are also studied, which are helpful to learn the transition's behavior.     

基于改进的Logvinovich法对带横摇运动非对称楔形体入水砰击力的计算(英文)

The water entry problem of an asymmetric wedge with roll motion was analyzed by the method of a modified Logvinovich model (MLM). The MLM is a kind of analytical model based on the Wagner method, which linearizes the free surface condition and body boundary condition. The difference is that the MLM applies a nonlinear Bernoulli equation to obtain pressure distribution, which has been proven to be helpful to enhance the accuracy of hydrodynamic loads. The Wagner condition in this paper was generalized to solve the problem of the water entry of a wedge body with rotational velocity. The comparison of wet width between the MLM and a fully nonlinear numerical approach was given, and they agree well with each other. The effect of angular velocity on the hydrodynamic loads of a wedge body was investigated.     

Characteristics of a bubble jet near a vertical wall

A numerical model of a coupled bubble jet and wall was built on the assumption of potential flow and calculated by the boundary integral method. A three-dimensional computing program was then developed. Starting with the basic phenomenon of the interaction between a bubble and a wall, the dynamics of bubbles near rigid walls were studied systematically with the program. Calculated results agreed well with experimental results. The relationship between the Bjerknes effect of a wall and characteristic parameters was then studied and the calculated results of various cases were compared and discussed with the Blake criterion based on the Kelvin-impulse theory. Our analyses show that the angle of the jet＇s direction and the pressure on the rigid wall have a close relationship with collapse force and the bubble＇s characteristic parameters. From this, the application range of Blake criterion can be determined. This paper aims to provide a basis for future research on the dynamics of bubbles near a wall.     

Design of underwater robot lines based on a hybrid automatic optimization strategy

In this paper, a hybrid automatic optimization strategy is proposed for the design of underwater robot lines. Isight is introduced as an integration platform. The construction of this platform is based on the user programming and several commercial software including UG6.0, GAMBIT2.4.6 and FLUENT12.0. An intelligent parameter optimization method, the particle swarm optimization, is incorporated into the platform. To verify the strategy proposed, a simulation is conducted on the underwater robot model 5470, which originates from the DTRC SUBOFF project. With the automatic optimization platform, the minimal resistance is taken as the optimization goal;the wet surface area as the constraint condition; the length of the fore-body, maximum body radius and after-body＇s minimum radius as the design variables. With the CFD calculation, the RANS equations and the standard turbulence model are used for direct numerical simulation. By analyses of the simulation results, it is concluded that the platform is of high efficiency and feasibility. Through the platform, a variety of schemes for the design of the lines are generated and the optimal solution is achieved. The combination of the intelligent optimization algorithm and the numerical simulation ensures a global optimal solution and improves the efficiency of the searching solutions.     

Fully nonlinear simulation for fluid/structure impact: A review

This paper presents a review of the work on fluid/structure impact based on inviscid and imcompressible liquid and irrotational flow. The focus is on the velocity potential theory together with boundary element method （BEM）. Fully nonlinear boundary conditions are imposed on the unknown free surface and the wetted surface of the moving body. The review includes （1） vertical and oblique water entry of a body at constant or a prescribed varying speed, as well as free fall motion, （2） liquid droplets or column impact as well as wave impact on a body, （3） similarity solution of an expanding body. It covers two dimensional （2D）, axisymmetric and three dimensional （3D） cases. Key techniques used in the numerical simulation are outlined, including mesh generation on the multivalued free surface, the stretched coordinate system for expanding domain, the auxiliary function method for decoupling the mutual dependence of the pressure and the body motion, and treatment for the jet or the thin liquid film developed during impact.     

Time domain simulation of a one line failure for a DP-assisted mooring system

This paper focuses on the research of a semi-submersible platform equipped with a DP-assisted mooring system. Based on the working principles of the DP-assisted mooring system and the model of the platform motion, a time domain simulation program is applied to analyze the impact, in the case of one line failure, on the platform motion, power consumption of the thrusters and the tension of the mooring lines. The results show that, under the 10-year wind dominant, a one line failure will have little impact on the tension of the mooring lines. When the failure line is windward, the power consumption will increase greatly with a weakened position of accuracy. However when the failure line is leeward, the power consumption will be reduced with a partly strengthened oosition of accuracy.     

Numerical simulation of the flow around two-dimensional partially cavitating hydrofoils

In the present study, a new approach is applied to the cavity prediction for two-dimensional （2D） hydrofoils by the potential based boundary element method （BEM）. The boundary element method is treated with the source and doublet distributions on the panel surface and cavity surface by usethe of the Dirichlet type boundary conditions. An iterative solution approach is used to determine the cavity shape on partially cavitating hydrofoils. In the case of a specified cavitation number and cavity length, the iterative solution method proceeds by addition or subtraction of a displacement thickness on the cavity surface of the hydrofoil. The appropriate cavity shape is obtained by the dynamic boundary condition of the cavity surface and the kinematic boundary condition of the whole foil surface including the cavity. For a given cavitation number the cavity length of the 2D hydrofoil is determined according to the minimum error criterion among different cavity lengths, which satisfies the dynamic boundary condition on the cavity surface. The NACA 16006, NACA 16012 and NACA 16015 hydrofoil sections are investigated for two angles of attack. The results are compared with other potential based boundary element codes, the PCPAN and a commercial CFD code （FLUENT）. Consequently, it has been shown that the results obtained from the two dimensional approach are consistent with those obtained from the others.     

The solution of minimum reliability arbitrary curved beam index for modified weak point of under pure bending

This paper is based on the fundamental loading model of pure bending and the analytical model of a circular beam element with arbitrary initial deflection. The L.W. Guo solution is modified and generalized according to the elastic theory, and the analytical solution for the stress of the beam element with arbitrary initial deflection under pure bending is derived. Using yield theory of edge strength, an expression for the safety margin of one point in the arbitrary curved beam under pure bending （ACPB） is built. This paper modifies the model for weak points of service structures and establishes a foundation for safe design and inspection of imperfect structures. Also, according to the theory of the method of advanced first-order second-moment（AFOSM） , this paper derives an expression for the reliability index of one point in ACPB. Lastly, it modifies the solution for weak points by solving the minimal reliability index.     

采用微气泡缩减高速船拖曳阻力(英文)

Ship hull form of the underwater area strongly influences the resistance of the ship.The major factor in ship resistance is skin friction resistance.Bulbous bows,polymer paint,water repellent paint(highly water-repellent wall),air injection,and specific roughness have been used by researchers as an attempt to obtain the resistance reduction and operation efficiency of ships.Micro-bubble injection is a promising technique for lowering frictional resistance.The injected air bubbles are supposed to somehow modify the energy inside the turbulent boundary layer and thereby lower the skin friction.The purpose of this study was to identify the effect of injected micro bubbles on a navy fast patrol boat(FPB) 57 m type model with the following main dimensions:L=2 450 mm,B=400 mm,and T=190 mm.The influence of the location of micro bubble injection and bubble velocity was also investigated.The ship model was pulled by an electric motor whose speed could be varied and adjusted.The ship model resistance was precisely measured by a load cell transducer.Comparison of ship resistance with and without micro-bubble injection was shown on a graph as a function of the drag coefficient and Froude number.It was shown that micro bubble injection behind the mid-ship is the best location to achieve the most effective drag reduction,and the drag reduction caused by the micro-bubbles can reach 6%-9%.     

流线型潜体首型对纵倾斜流涡结构的影响数值模拟（英文）

A flow field around a streamlined body at an intermediate angle of incidence is dominated by cross-flow separation and vortical flow fields. The separated flow leads to a pair of vortices on the leeside of the body; therefore, it is essential to accurately determine this pair and estimate its size and location. This study utilizes the element-based finite volume method based on RANS equations to compute a 3D axisymmetric flow around a SUBOFF bare submarined hull. Cross-flow vortex structures are then numerically simulated and compared for a submarine with SUBOFF and DRDC STR bows. Computed results of pressure and shear stress distribution on the hull surface and the strength and locations of the vortex structures are presented at an intermediate incidence angle of 20°. A wind tunnel experiment is also conducted to experimentally visualize the vortex structures and measure their core locations. These experimental results are compared with the numerical data, and a good agreement is found.     

Dynamic coupled analysis of the floating platform using the asynchronous coupling algorithm

This paper discusses the numerical modeling of the dynamic coupled analysis of the floating platform and mooring/risers using the asynchronous coupling algorithm with the purpose to improve the computational efficiency when multiple lines are connected to the platform. The numerical model of the platform motion simulation in wave is presented. Additionally, how the asynchronous coupling algorithm is implemented during the dynamic coupling analysis is introduced. Through a comparison of the numerical results of our developed model with commercial software for a SPAR platform, the developed numerical model is checked and validated.