采用基于雷诺平均的N-S方程和势流方法进行散货船阻力和运动的预报研究（英文）

Resistance prediction of ships using computational fluid dynamics has gained popularity over the years because of its high accuracy and low cost. This paper conducts numerical estimations of the ship resistance and motion of a Japan bulk carrier model using SHIP_Motion, a Reynolds-averaged Navier–Stokes(RaNS)-based solver, and HydroSTAR, a commercial potential flow(PF)-based solver. The RaNS solver uses an overset-structured mesh and discretizes the flow field using the finite volume method, while the PF-based solver applies the three-dimensional panel method. In the calm water test, the total drag coefficient,sinkage, and trim were predicted using the RaNS solver following mesh dependency analysis, and the results were compared with the available experimental data. Next, calm water resistance was investigated for a range of Froude numbers. The added resistance in short-wave cases was simulated using both Ra NS and PF solvers, and the results were compared. The PF solver showed better agreement with the RaNS solver for predicting motion responses than for predicting added resistance. While the added resistance results could not be directly validated because of the absence of experimental data, considering the previous accuracy of the RaNS solver in added resistance prediction and general added resistance profile of similar hull forms(bulk carriers), the prediction results could be concluded to be reliable.     

Ship hull plate processing surface fairing with constraints based on B-spline

The problem of ship hull plate processing surface fairing with constraints based on B-spline is solved in this paper. The algorithm for B-spline curve fairing with constraints is one of the most common methods in plane curve fairing. The algorithm can be applied to global and local curve fairing. It can constrain the perturbation range of the control points and the shape variation of the curve, and get a better fairing result in plane curves. In this paper, a new fairing algorithm with constraints for curves and surfaces in space is presented. Then this method is applied to the experiments of ship hull plate processing surface. Finally numerical results are obtained to show the efficiency of this method.     

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

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.     

Determining the hydrodynamic forces on a planing hull in steady motion

A combination of methods was developed that can determine hydrodynamic forces on a planing hull in steady motion. Firstly, a potential-based boundary-element method was used to calculate the hydrodynamic pressure, induced resistance and lift. Then the frictional resistance component was determined by the viscous boundary layer theory. Finally, a particular empirical technique was applied, to determine the region of upwash geometry and determine spray resistance. Case studies involving four models of Series 62 planing craft were run. These showed that the suggested method is efficient and capable, with results that are in good agreement with experimental measurements over a wide range of volumetric Froude numbers.     

Determining the hydrodynamic forces on a planing hull in steady motion

A combination of methods was developed that can determine hydrodynamic forces on a planing hull in steady motion.Firstly,a potential-based boundary-element method was used to calculate the hydrodynamic pressure,induced resistance and lift.Then the frictional resistance component was determined by the viscous boundary layer theory.Finally,a particular empirical technique was applied.to determine the region of upwash geometry and determine spray resistance.Case studies involving four models of Series 62 planing craft were run.These showed that the suggested method is efficient and capable,with results that are in good agreement with experimental measurements over a wide range of volumetric Froude numbers.     

航行于波浪中的肥大型船波浪增阻计算研究(英文)

Under the background of the energy saving and emission reduction, more and more attention has been placed on investigating the energy efficiency of ships. The added resistance has been noted for being crucial in predicting the decrease of speed on a ship operating at sea. Furthermore, it is also significant to investigate the added resistance for a ship functioning in short waves of large modern ships. The researcher presents an estimation formula for the calculation of an added resistance study in short waves derived from the reflection law. An improved method has been proposed to calculate the added resistance due to ship motions, which applies the radiated energy theory along with the strip method. This procedure is based on an extended integral equation (EIE) method, which was used for solving the hydrodynamic coefficients without effects of the irregular frequency. Next, a combined method was recommended for the estimation of added resistance for a ship in the whole wave length range. The comparison data with other experiments indicate the method presented in the paper provides satisfactory results for large blunt ship.     

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.     

Free transverse vibration analysis of an underwater launcher based on fluid-structure interaction

A pneumatic launcher is theoretically investigated to study its natural transverse vibration in water. Considering the mass effect of the sealing cover, the launcher is simplified as a uniform cantilever beam with a top point mass. By introducing the boundary and continuity conditions into the motion equation, the natural frequency equation and the mode shape function are derived. An iterative calculation method for added mass is also presented using the velocity potential function to account for the mass effect of the fluid on the launcher. The first 2 order natural frequencies and mode shapes are discussed in external flow fields and both external and internal flow fields. The results show good agreement with both natural frequencies and mode shapes between the theoretical analysis and the FEM studies. Also, the added mass is found to decrease with the increase of the mode shape orders of the launcher. And because of the larger added mass in both the external and internal flow fields than that in only the external flow field, the corresponding natural frequencies of the former are relatively smaller.     

二维冲箱式造波机的间隙和形状效应数值研究（英文）

The installation of plunger-type wave makers in experimental tanks will generally include a gap between the back of the wedge and the wall of the tank. In this study, we analyze the influence of this gap on the wave making performance of the plunger using two-dimensional(2 D) CFD calculations for a range of nearly linear wave conditions and compare the results with both experimental measurements and linear potential flow theory. Three wedge-shaped profiles, all with the same submerged volume, are considered. Moreover, the generated waves are compared with the predictions of linear potential flow theory. The calculations are made using the commercial ANSYS FLUENT finite-volume code with dynamic meshes to solve the Navier–Stokes equations and the volume of fluid scheme to capture the air–water interface. Furthermore, the linear potential flow solution of Wu(J Hydraul Res 26:481–493, 1988) is extended to consider an arbitrary profile and serve as a reference solution. The amplitude ratios of the generated waves predicted by the CFD calculations compare well with the predictions of linear potential flow theory for a simple wedge, indicating that viscous effects do not influence this ratio for small-amplitude motions in 2 D. By contrast,significant higher harmonic components are produced by larger amplitude motions. Also, the simple wedge is found to produce the smallest spurious higher harmonic content in the far-field wave.     

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.     

逆压边界层流动的湍流模式研究（英文）

In this paper, a numerical study of flow in the turbulence boundary layer with adverse and pressure gradients(APGs) is conducted by using Reynolds-averaged Navier-Stokes(RANS) equations. This research chooses six typical turbulence models, which are critical to the computing precision, and to evaluating the issue of APGs. Local frictional resistance coefficient is compared between numerical and experimental results. The same comparisons of dimensionless averaged velocity profiles are also performed. It is found that results generated by Wilcox(2006) k-ω are most close to the experimental data. Meanwhile, turbulent quantities such as turbulent kinetic energy and Reynolds-stress are also studied.     

Turbulence Model Investigations on the Boundary Layer Flow with Adverse Pressure Gradients

In this paper, a numerical study of flow in the turbulence boundary layer with adverse and pressure gradients(APGs) is conducted by using Reynolds-averaged Navier-Stokes(RANS) equations. This research chooses six typical turbulence models, which are critical to the computing precision, and to evaluating the issue of APGs. Local frictional resistance coefficient is compared between numerical and experimental results. The same comparisons of dimensionless averaged velocity profiles are also performed. It is found that results generated by Wilcox(2006) k-ω are most close to the experimental data. Meanwhile, turbulent quantities such as turbulent kinetic energy and Reynolds-stress are also studied.     

半无限惯性表面产生的波浪散射问题研究的另一种方法(英文)

A new method to solve the boundary value problem arising in the study of scattering of two-dimensional surface water waves by a discontinuity in the surface boundary conditions is presented in this paper. The discontinuity arises due to the floating of two semi-infinite inertial surfaces of different surface densities. Applying Green’s second identity to the potential functions and appropriate Green’s functions, this problem is reduced to solving two coupled Fredholm integral equations with regular kernels. The solutions to these integral equations are used to determine the reflection and the transmission coefficients. The results for the reflection coefficient are presented graphically and are compared to those obtained earlier using other research methods. It is observed from the graphs that the results computed from the present analysis match exactly with the previous results.     

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

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.     

机翼在接近地面的起飞模拟数学模型（英文）

Aircraft flying close to the ground benefit from enhanced efficiency owing to decreased induced drag and increased lift. In this study, a mathematical model is developed to simulate the takeoff of a wing near the ground using an Iterative Boundary Element Method(IBEM) and the finite difference scheme. Two stand-alone sub-codes and a mother code, which enables communication between the sub-codes, are developed to solve for the self-excitation of the Wing-In-Ground(WIG) effect. The aerodynamic force exerted on the wing is calculated by the first sub-code using the IBEM, and the vertical displacement of the wing is calculated by the second sub-code using the finite difference scheme. The mother code commands the two sub-codes and can solve for the aerodynamics of the wing and operating height within seconds. The developed code system is used to solve for the force, velocity, and displacement of an NACA6409 wing at a 4° Angle of Attack(AoA) which has various numerical and experimental studies in the literature. The effects of thickness and AoA are then investigated and conclusions were drawn with respect to generated results. The proposed model provides a practical method for understanding the flight dynamics and it is specifically beneficial at the pre-design stages of a WIG effect craft.     

Optimization of bow shape for a non ballast water ship

In this research,a commercial CFD code ＂Fluent＂ was applied to optimization of bulbous bow shape for a non ballast water ships（NBS）.The ship was developed at the Laboratory of the authors in Osaka Prefecture University,Japan.At first,accuracy of the CFD code was validated by comparing the CFD results with experimental results at towing tank of Osaka Prefecture University.In the optimizing process,the resistances acting on ships in calm water and in regular head waves were defined as the object function.Following features of bulbous bow shapes were considered as design parameters： volume of bulbous bow,height of its volume center,angle of bow bottom,and length of bulbous bow.When referring to the computed results given by the CFD like resistance,pressure and wave pattern made by ships in calm water and in waves,an optimal bow shape for ships was discovered by comparing the results in the series of bow shapes.In the computation on waves,the ship is in fully captured condition because shorter waves,λ/Lpp 0.6,are assumed.     

The method of imbedded Lagrangian element to estimate wave power absorption by some submerged devices

A simple approach is described to estimate the wave power absorption potential of submerged devices known to cause wave focusing and flow enhancement. In particular, the presence of a flow-through power take-off （PTO） system, such as low-head turbines, can be accounted for. The wave radiation characteristics of an appropriately selected Lagrangian element （LE） in the fluid domain are first determined. In the limit of a vanishing mass, the LE reduces to a patch of distributed normal dipoles. The hydrodynamic coefficients of this virtual object are then input in a standard equation of motion where the effect of the PTO can be represented, for example, as a dashpot damping term. The process is illustrated for a class of devices recently proposed by Carter and Ertekin （2011）, although in a simplified form. Favorable wave power absorption is shown for large ratios of the LE wave radiation coefficient over the LE added mass coefficient. Under optimal conditions, the relative flow reduction from the PTO theoretically lies between 0.50 and 1 2 ≈ 0.71, with lower values corresponding to better configurations. Wave power capture widths, the sensitivity of results to PTO damping and sample spectral calculations at a typical site in Hawaiian waters are proposed to further illustrate the versatility of the method.     

Numerical simulation of viscous liquid sloshing by moving-particle semi-implicit method

A meshless numerical simulation method,the moving-particle semi-implicit method(MPS)is presented in this paper to study the sloshing phenomenon in ocean and naval engineering.As a meshless method,MPS uses particles to replace the mesh in traditional methods,the governing equations are discretized by virtue of the relationship of particles,and the Poisson equation of pressure is solved by incomplete Cholesky conjugate gradient method(ICCG),the free surface is tracked by the change of numerical density.A numerical experiment of viscous liquid sloshing tank was presented and compared with the result got by the difference method with the VOF,and an additional modification step was added to make the simulation more stable.The results show that the MPS method is suitable for the simulation of viscous liquid sloshing,with the advantage in arranging the particles easily,especially on some complex curved surface.     

Numerical simulation of viscous liquid sloshing by moving-particle semi-implicit method

A meshless numerical simulation method, the moving-particle semi-implicit method （MPS） is presented in this paper to study the sloshing phenomenon in ocean and naval engineering. As a meshless method, MPS uses particles to replace the mesh in traditional methods, the governing equations are discretized by virtue of the relationship of particles, and the Poisson equation of pressure is solved by incomplete Cholesky conjugate gradient method （ICCG）, the free surface is tracked by the change of numerical density. A numerical experiment of viscous liquid sloshing tank was presented and compared with the result got by the difference method with the VOF, and an additional modification step was added to make the simulation more stable. The results show that the MPS method is suitable for the simulation of viscous liquid sloshing, with the advantage in arranging the particles easily, especially on some complex curved surface.