Influence of slip boundary on flow separation and drag of flow past bluff body at high Reynolds numbers北大核心CSCD

［Objectives］Flow separation increases the drag and noise of underwater vehicles, and influences the controllability of their control surfaces. Therefore, the influence of slip caused by superhydrophobic surfaces on drag reduction and flow separation is studied. ［Methods］A partial slip boundary condition is developed, and the flow around a circular cylinder and foil with a slip boundary at high Reynolds numbers are numerically simulated. ［Results］The results show that the when the slip length increases, the flow around the cylinder goes through three stages: the turbulent Kármán vortex street, laminar Kármán vortex street and non-separation Stokes flow. The drag coefficient increases first and then decreases, and the vortex shedding frequency increases. For flow around a foil, the separation position moves downstream until the separation region disappears when the slip length increases, and the drag coefficient decreases while the lift coefficient increases. ［Conclusions］The results of this study show that for flow past bluff body at high Reynolds number, the slip boundary can control flow separation and reduce drag effectively, providing technical support for the application of superhydrophobic surfaces for the flow control of underwater vehicle appendages. © 2022 Journal of Clinical Hepatology. All rights reserved.     

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

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.     

开源环境下的高精确度船舶水动力学仿真(英文)

The numerical simulation of wake and free-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 Open FOAM. 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.     

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.     

Numerical study of viscous wave-making resistance of ship navigation in still water

The prediction of a ship＇s resistance especially the viscous wave-making resistance is an important issue in CFD applications. In this paper, the resistances of six ships from hull 1 to hull 6 with different hull forms advancing in still water are numerically studied using the solver naoe-FOAM-SJTU, which was developed based on the open source code package OpenFOAM. Different components of the resistances are computed and compared while considering two speed conditions （12 kn and 16 kn）. The resistance of hull 3 is the smallest while that of hull 5 is the largest at the same speed. The results show hull 3 is a good reference for the design of similar ships, which can provide some valuable guidelines for hull form optimization.     

基于B样条应用的船舶波浪增阻计算（英文）

Making an exact computation of added resistance in sea waves is of high interest due to the economic effects relating to ship design and operation. In this paper, a B-spline based method is developed for computation of added resistance. Based on the potential flow assumption, the velocity potential is computed using Green's formula. The Kochin function is applied to compute added resistance using Maruo's far-field method, the body surface is described by a B-spline curve and potentials and normal derivation of potentials are also described by B-spline basis functions and B-spline derivations. A collocation approach is applied for numerical computation, and integral equations are then evaluated by applying Gauss–Legendre quadrature. Computations are performed for a spheroid and different hull forms; results are validated by a comparison with experimental results. All results obtained with the present method show good agreement with experimental results.     

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

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.     

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

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.     

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.     

船体周围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 analysis of the high skew propeller of an underwater vehicle

A numerical analysis based on the boundary element method (BEM) was presented for the hydrodynamic performance of a high skew propeller (HSP) which is employed by an underwater vehicle (UV). Since UVs operate at two different working conditions (surface and submerged conditions), the design of such a propeller is a cumbersome task. This is primarily due to the fact that the resistance forces as well as the vessel efficiency under these conditions are significantly different. Therefore, some factors are necessary for the design of the optimum propeller to utilize the power at the mentioned conditions. The design objectives of the optimum propeller are to obtain the highest possible thrust, minimum torque, and efficiency. In the current study, a 5-bladed HSP was chosen for running the UV. This propeller operated at the stern of the UV hull where the inflow velocity to the propeller was non-uniform. Some parameters of the propeller were predicted based on the UV geometrical hull and operating conditions. The computed results include the pressure distribution and the hydrodynamic characteristics of the HSP in open water conditions, and comparison of these results with those of the experimental data indicates good agreement. The propeller efficiency for both submerged and surface conditions was found to be 67% and 64%, respectively, which compared to conventional propellers is a significantly higher efficiency.     

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.     

Numerical solutions for a two-dimensional airfoil undergoing unsteady motion

Continuous vorticity panels are used to model general unsteady inviscid, incompressible, and two-dimensional flows. The geometry of the airfoil is approximated by series of short straight segments having endpoints that lie on the actual surface. A piecewise linear, continuous distribution of vorticity over the airfoil surface is used to generate disturbance flow. The no-penetration condition is imposed at the midpoint of each segment and at discrete times. The wake is simulated by a system of point vortices, which move at local fluid velocity. At each time step, a new wake panel with uniform vorticity distribution is attached to the trailing edge, and the condition of constant circulation around the airfoil and wake is imposed. A new expression for Kutta condition is developed to study (i) the effect of thickness on the lift build-up of an impulsively started airfoil, (ii) the effects of reduced frequency and heave amplitude on the thrust production of flapping airfoils, and (iii) the vortex-airfoil interaction. This work presents some hydrodynamic results for tidalstream turbine.     

一种水翼空化的被动控制方法（英文）

This paper introduces a new idea of controlling cavitation around a hydrofoil through a passive cavitation controller called artificial cavitation bubble generator(ACG). Cyclic processes, namely, growth and implosion of bubbles around an immersed body, are the main reasons for the destruction and erosion of the said body. This paper aims to create a condition in which the cavitation bubbles reach a steady-state situation and prevent the occurrence of the cyclic processes. For this purpose, the ACG is placed on the surface of an immersed body, in particular, the suction surface of a 2 D hydrofoil. A simulation was performed with an implicit finite volume scheme based on a SIMPLE algorithm associated with the multiphase and cavitation model. The modified k-ε RNG turbulence model equipped with a modification of the turbulent viscosity was applied to overcome the turbulence closure problem. Numerical simulation of water flow over the hydrofoil equipped with the ACG shows that a low-pressure recirculation area is produced behind the ACG and artificially generates stationary cavitation bubbles. The location, shape, and size of this ACG are the crucial parameters in creating a proper control. Results show that the cavitation bubble is controlled well with a well-designed ACG.     

水动力效率最高的鳗鱼推进器(英文)

It is hypothesized that steady anguilliform swimming motion of aquatic animals is purely reactive such that no net vortex wake is left downstream. This is versus carangiform and tunniform swimming of fish, where vortex streams are shed from tail, fins, and body. But there the animal movements are such to produce partial vortex cancellation downstream in maximizing propulsive efficiency. In anguilliform swimming characteristic of the eel family, it is argued that the swimming motions are configured by the animal such that vortex shedding does not occur at all. However, the propulsive thrust in this case is higher order in the motion amplitude, so that relatively large coils are needed to produce relatively small thrust; the speeds of anguilliform swimmers are less than the carangiform and tunniform, which develop first order thrusts via lifting processes. Results of experimentation on live lamprey are compared to theoretical prediction which assumes the no-wake hypothesis. Two-dimensional analysis is first performed to set the concept. This is followed by three-dimensional analysis using slender-body theory. Slender-body theory has been applied by others in studying anguilliform swimming, as it is ideally suited to the geometry of the lamprey and other eel-like animals. The agreement between this new approach based on the hypothesis of wakeless swimming and the experiments is remarkably good in spite of the physical complexities.     

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

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.     

Experimental Research on Flash Boiling Spray of Dimethyl Ether

The high-speed digital imaging technique is applied to observe the developing process of flash boiling spray of dimethyl ether at low ambient pressure, and the effects of nozzle opening pressure and nozzle hole diameter on the spray shape, spray tip penetration and spray angle during the injection are investigated. The experimental results show that the time when the vortex ring structure of flash boiling spray forms and its developing process are determined by the combined action of the bubble growth and breakup in the spray and the air drag on the leading end of spray; with the enhancement of nozzle opening pressure, the spray tip penetration increases and the spray angle decreases. The influence of nozzle hole diameter on the spray tip penetration is relatively complicated, the spray tip penetration is longer with a smaller nozzle hole diameter at the early stage of injection, while the situation is just opposite at the later stage of injection. This paper establishes that the variation of spray angle is consistent with that of nozzle hole diameter.     

Numerical simulation on truncated oval-nosed projectile penetrating into stiffened plate北大核心CSCD

[Objectives]In this paper, the numerical simulation method is used to study the anti-penetration performance and energy absorption mode of a stiffened plate, as well as the influence of different stiffened bars on the flight attitude of the projectile body.[Methods] Finite element software LS-DYNA is used to simulate the process of a truncated oval-nosed projectile penetrating a stiffened plate, and the results of the numerical simulation are compared with an experiment to verify the reliability of the numerical simulation method. The momentum method and mass equivalence method are used to predict the residual velocity of the projectile, and the applicability of different theoretical methods within different velocity ranges is compared. The deformation energy of different regions of the stiffened plate is then extracted to analyze the influence of the initial velocity of the projectile body on the energy absorption mode of the target plate. Finally, the structure of the stiffeners is changed and the influence of the relative position of the stiffeners on the penetration attitude of the projectile body is analyzed.[Results]The results show that the mass equivalence method is more accurate than the momentum method in predicting the residual velocity of the stiffened plate when the initial velocity of the projectile body is in the range of 300–900 m/s. The ratio of the deformation energy of the stiffened plate to the energy loss of the projectile body decreases with the increase of the initial velocity of the projectile body. The effect of a T-stiffened plate on trajectory is greater than that of a rectangular-stiffened plate.[Conclusions]The related calculation method and research results have certain reference value for research and engineering application surrounding the anti-penetration of stiffened plates. © 2023 Chinese Journal of Ship Research. All rights reserved.     

液舱晃荡与浮体非线性耦合运动分析(英文)

Nonlinear interactions among incident wave, tank-sloshing and floating body coupling motion are investigated. The fully nonlinear sloshing and body-surface nonlinear free surface hydrodynamics is simulated using a Non-Uniform Rational B-Spline (NURBS) higher-order panel method in time domain based on the potential theory. A robust and stable improved iterative procedure (Yan and Ma, 2007) for floating bodies is used for calculating the time derivative of velocity potential and floating body motion. An energy dissipation condition based on linear theory adopted by Huang (2011) is developed to consider flow viscosity effects of sloshing flow in nonlinear model. A two-dimensional tank model test was performed to identify its validity. The present nonlinear coupling sway motion results are subsequently compared with the corresponding Rognebakke and Faltinsen (2003)’s experimental results, showing fair agreement. Thus, the numerical approach presented in this paper is expected to be very efficient and realistic in evaluating the coupling effects of nonlinear sloshing and body motion.     

Mathematical simulation of steep waves at a focus point

Mathematical models simulating steep waves at a focus point are presented in this paper. Simulations of extreme waves in a model basin were used to determine the loads on floating structures induced by the waves. Based on a new wave theory, numerical test results show that the simulation procedure is effective and the induced motion of water particles in the front of waves is an important factor influencing impact loads on floating bodies.