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.     

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

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.     

Application of dual reciprocity boundary element method to predict acoustic attenuation characteristics of marine engine exhaust silencers

In marine engine exhaust silencing systems, the presence of exhaust gas flow influences the sound propagation inside the systems and the acoustic attenuation performance of silencers. In order to investigate the effects of three-dimensional gas flow and acoustic damping on the acoustic attenuation characteristics of marine engine exhaust silencers, a dual reciprocity boundary element method （DRBEM） was developed. The acoustic governing equation in three-dimensional potential flow was derived first, and then the DRBEM numerical procedure is given. Compared to the conventional boundary element method （CBEM）, the DRBEM considers the second order terms of flow Mach number in the acoustic governing equation, so it is suitable for the cases with higher Mach number subsonic flow. For complex exhaust silencers, it is difficult to apply the single-domain boundary element method, so a substructure approach based on the dual reciprocity boundary element method is presented. The experiments for measuring transmission loss of silencers are conducted, and the experimental setup and measurements are explained. The transmission loss of a single expansion chamber silencer with extended inlet and outlet were predicted by DRBEM and compared with the measurements. The good agreements between predictions and measurements are observed, which demonstrated that the derived acoustic governing equation and the DRBEM numerical procedure in the present study are correct.     

The choice of speed and clearance for RAS on 3D method

In this paper, a 3D source distribution technique is used to calculate the coupled motions between two ships which advance in the wave with the same speed. The numerical results of coupled motions for a frigate and a supply ship have a good agreement with the experimental results. Based on the 3D coupled motions of two ships, a spectral analysis is employed to clearly observe the effect of speed, clearance and wave heading on the significant relative motion amplitude (SRMA) of two ships.The method presented in this paper will be helpful to select suitable clearance, speed and wave heading for underway replenishment at sea(RAS).     

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

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.     

基于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.     

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.     

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

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.     

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

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.     

Influence of exhaust parameters on mid-infrared radiative characteristics of the plume北大核心CSCD

[Objectives] Infrared radiative (IR) characteristics of the plume are the main part of ship infrared radiation. The influence analysis of exhaust parameters on radiative characteristics is very important for developing ship stealth and anti-stealth technology. [Methods]In this work, the line-by-line method is applied to solve radiative transfer in gaseous media. Infrared radiative images of plumes under different temperatures and flow rates are obtained. Radiance, spectral radiative intensity, and spectrum-integrated radiative intensity over 3−5 µm are calculated and analyzed. [Results]Results show that there are obvious peaks of spectral radiative intensity curves at 4.18µm and 4.43 µm, and an obvious valley at 4.26 µm. [Conclusions] The infrared radiation characteristics of the ships exhaust plume are greatly affected by the exhaust temperature, but less affected by the exhaust flow rate. © 2023 The Author(s).     

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.     

Numerical simulation of three-dimensional combustion flows

A finite-rate method is used to simulate the three-dimensional combustion process in a plasma generator with CH4 as the fuel. The simulation was run with RNG k-ε model to simulate turbulence, with eddy-dissipation-concept （EDC） model to simulate the combustion and with discrete ordinates model to simulate radiation. The numerical results show that the flow field characteristics and the parameter distributions are under the condition of rich fuels, and these results provide valuable information when optimizing the plasma generator design and organizing its flow fields.     

基于模糊集合的船舶燃气轮机临界因素的Dematel评估方法（英文）

In power production, gas turbines are commonly used components that generate high amount of energy depending on size and weight. They function as integral parts of helicopters, aircrafts, trains, ships, electrical generators, and tanks. Notably, many researchers are focusing on the design, operation, and maintenance of gas turbines. The focal point of this paper is a DEMATEL approach based on fuzzy sets, with the attempt to use these fuzzy sets explicitly. Using this approach, the cause–effect diagram of gas turbine failures expressed in the literature is generated and aimed to create a perspective for operators. The results of the study show that, "connecting shaft has been broken between turbine and gear box" selected the most important cause factor and"sufficient pressure fuel does not come for fuel pump" is selected the most important effect factor, according to the experts.     

3D numerical modeling of wave forces on tandem fixed cylinders using the BEM

In this paper a 3D numerical model was developed to study the complicated interaction between waves and a set of tandem fixed cylinders.The fluid was considered to be inviscid and irrotational.Therefore,the Helmholtz equation was used as a governing equation.The boundary element method（BEM） was adopted to discretize the relevant equations.Open boundaries were used in far fields of the study domain.Linear waves were generated and propagated towards tandem fixed cylinders to estimate the forces applied on them.Special attention was paid to consideration of the effect on varying non-dimensional cylinder radius and distance between cylinders,ka and kd on forces and trapped modes.The middle cylinder wave forces and trapped modes in a set of nine tandem cylinders were validated utilizing analytical data.The comparisons confirm the accuracy of the model.The results of the inline wave force estimation on n tandem cylinders show that the critical cylinder in the row is the middle one for odd numbers of cylinders.Furthermore the results show that the critical trapped mode effect occurs for normalized cylinder radiuses close to 0.5 and 1.0.Finally the force estimation for n tandem cylinders confirms that force amplitude of the middle cylinder versus normalized separation distance fluctuates about that of a single cylinder.     

Research on vibration and near-field soundradiation of ring-ribbed cylindrical shell coated deadening and decoupling materials

The research on structural vibration and sound radiation of underwater ring-ribbed cylindrical shell, which is coated with a kind of deadening and decoupling materials, becomes a focus in recent years. This paper analyzes the problem on two aspects: model experiment and numerical calculation. The model experiment is carried out including three cases firstly, in which the structural vibration response and radiating acoustic field are measured respectively, and the results gained in these three cases are analyzed to discuss the effect of reducing structural vibration and radiating noise of the deadening and decoupling materials. The coupling FEM/BEM and the SEA methods are both used in numerical calculation, i.e. the arithmetic of the coupling FEM/BEM method is adopted to calculate the low frequency characteristics and the SEA method is adopted to calculate the medium-high frequencies characteristics of the model. By comparing experimental results with numerical calculation results, it is proved that the algorithm adopted in this paper is reasonable.     

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.     

Free Surface Flow Generated by Submerged Twin-cylinders in Forward Motion Using a Fully Nonlinear Method

The free surface flow generated by twin-cylinders in forced motion submerged beneath the free surface is studied based on the boundary element method. Two relative locations, namely, horizontal and vertical, are examined for the twin cylinders. In both cases, the twin cylinders are starting from rest and ultimately moving with the same constant speed through an accelerating process. Assuming that the fluid is inviscid and incompressible and the flow to be irrotational, the integral Laplace equation can be discretized based on the boundary element method. Fully-nonlinear boundary conditions are satisfied on the unknown free surface and the moving body surface. The free surface is traced by a Lagrangian technique. Regriding and remeshing are applied, which is crucial to quality of the numerical results. Single circular cylinder and elliptical cylinder are calculated by linear method and fully nonlinear method for accuracy checking and then fully nonlinear method is conducted on the twin cylinder cases, respectively. The generated wave elevation and the resultant force are analysed to discuss the influence of the gap between the two cylinders as well as the water depth. It is found that no matter the kind of distribution, when the moving cylinders are close to each other, they suffer hydrodynamic force with large absolute value in the direction of motion. The trend of force varying with the increase of gap can be clearly seen from numerical analysis. The vertically distributed twin cylinders seem to attract with each other while the horizontally distributed twin cylinders are opposite when they are close to each other.     

3-D computational method of wave loads on turret moored FPSO tankers

A three-dimensional method of calculating wave loads of turret moored FPSO (Floating Production Storage and Offloading) tankers is presented. The linearized restoring forces acting on the ship hull by the mooring system are calculated according to the catenary theory, which are expressed as the function of linear stiffness coefficients and the displacements of the upper ends of mooring chains. The hydrodynamic coefficients of the ship are calculated by the three-dimensional potential flow theory of the linear hydrodynamic problem for ships with a low forward speed. The equations of ship motions are established with the effect of the restoring forces from the mooring system included as linear stiffness coefficients. The equations of motions are solved in frequency domain, and the responses of wave-induced motions and loads on the ship can be obtained. A computer program based on this method has been developed, and some calculation examples are illustrated. Analysis results show that the method can give satisfying prediction of wave loads.     

Application of large eddy simulation in the performance study of wave blocker

The configuration and aerodynamic performance of the inlet system are important aspects in the process of installing a gas turbine on a naval vessel. Under the requirements, large eddy simulation （LES） is used to simulate the three-dimensional fluid flow in the wave blocker of a marine inlet filter. The Smagorinsky-Lilly sub-grid model was used to model motions of small-scale structures. During numerical simulation, the SIMPLE algorithm was applied. The central-differencing spatial discretization scheme and the second order accuracy finite difference for the temporal discretization were used. Simulation gives satisfactory distribution of the vorticity fields and turbulent kinetic energy. Compared with the k-ε turbulent model, the results of LES are better for the distribution of parameters. The results of experimental study in a small-scale wind tunnel indicate that numerical calculation has higher accuracy. Therefore, the methods used are worthy of reference and introduction for the design of an inlet system.