Monolithic coupling of the pressure and rigid body motion equations in computational marine hydrodynamics

In Fluid Structure Interaction(FSI) problems encountered in marine hydrodynamics, the pressure field and the velocity of the rigid body are tightly coupled. This coupling is traditionally resolved in a partitioned manner by solving the rigid body motion equations once per nonlinear correction loop, updating the position of the body and solving the fluid flow equations in the new configuration. The partitioned approach requires a large number of nonlinear iteration loops per time–step. In order to enhance the coupling, a monolithic approach is proposed in Finite Volume(FV) framework,where the pressure equation and the rigid body motion equations are solved in a single linear system. The coupling is resolved by solving the rigid body motion equations once per linear solver iteration of the pressure equation, where updated pressure field is used to calculate new forces acting on the body, and by introducing the updated rigid body boundary velocity in to the pressure equation. In this paper the monolithic coupling is validated on a simple 2D heave decay case. Additionally, the method is compared to the traditional partitioned approach(i.e. "strongly coupled" approach) in terms of computational efficiency and accuracy. The comparison is performed on a seakeeping case in regular head waves, and it shows that the monolithic approach achieves similar accuracy with fewer nonlinear correctors per time–step. Hence, significant savings in computational time can be achieved while retaining the same level of accuracy.     

Control of Vortex Shedding and Drag Reduction through Dual Splitter Plates Attached to a Square Cylinder

In this paper we have made a numerical study on the control of vortex shedding and drag reduction of a cylinder by attaching thin splitter plates. The wake structure of the cylinder of square cross-section with attached splitter plates is analyzed for a range of Reynolds number, based on the incident stream and height of the cylinder, in the laminar range. The Navier-Stokes equations governing the flow are solved by the control volume method over a staggered grid arrangement. We have used the semi-implicit method for pressure-linked equation(SIMPLE) algorithm for computation. Our results show that the presence of a splitter plate upstream of the cylinder reduces the drag, but it has a small impact on the vortex shedding frequency when the plate length is beyond 1.5 time the height of the cylinder. The presence of a downstream splitter plate dampens the vortex shedding frequency. The entrainment of fluid into the inner side of the separated shear layers is obstructed by the downstream splitter plate. Our results suggest that by attaching in-line splitter plates both upstream and downstream of the cylinder, the vortex shedding can be suppressed, as well as a reduction in drag be obtained. We made a parametric study to determine the optimal length of these splitter plates so as to achieve low drag and low vortex shedding frequency.     

Solving Nonlinear Differential Equation Governing on the Rigid Beams on Viscoelastic Foundation by AGM

In the present paper a vibrational differential equation governing on a rigid beam on viscoelastic foundation has been investigated. The nonlinear differential equation governing on this vibrating syst...     

采用局部的微分求积法求解激波管问题

The localized differential quadrature (LDQ) method is a numerical technique with high accuracy for solving most kinds of nonlinear problems in engineering and can overcome the difficulties of other methods (such as difference method) to numerically evaluate the derivatives of the functions.Its high efficiency and accuracy attract many engineers to apply the method to solve most of the numerical problems in engineering.However,difficulties can still be found in some particular problems.In the following study,the LDQ was applied to solve the Sod shock tube problem.This problem is a very particular kind of problem,which challenges many common numerical methods.Three different examples were given for testing the robustness and accuracy of the LDQ.In the first example,in which common initial conditions and solving methods were given,the numerical oscillations could be found dramatically;in the second example,the initial conditions were adjusted appropriately and the numerical oscillations were less dramatic than that in the first example;in the third example,the momentum equation of the Sod shock tube problem was corrected by adding artificial viscosity,causing the numerical oscillations to nearly disappear in the process of calculation.The numerical results presented demonstrate the detailed difficulties encountered in the calculations,which need to be improved in future work.However,in summary,the localized differential quadrature is shown to be a trustworthy method for solving most of the nonlinear problems in engineering.     

SPSM and its application in cylindrical shells

In naval architectures, the structure of prismatic shell is used widely. But there is no suitable method to analyze this kind of structure. Stiffened prismatic shell method (SPSM) presented in this paper, is one of the harmonic semi-analytic methods. Theoretically, strong stiffened structure can be analyzed economically and accurately. SPSM is based on the analytical solution of the governing differential equations for orthotropic cylindrical shells. In these differential equations, the torsional stiffness, bending stiffness and the exact position of each stiffener are taken into account with the Heaviside singular function. An algorithm is introduced, in which the actions of stiffeners on shells are replaced by external loads at each stiffener position. Stiffened shells can be computed as non-stiffened shells. Eventually, the displacement solution of the equations is acquired by the introduction of Green function. The stresses in a corrugated transverse bulkhead without pier base of an oil tanker are computed by using SPSM.     

Analyzing parabolic profile path for underwater towed-cable

This article discusses the dynamic state analysis of underwater towed-cable when tow-ship changes its speed in a direction making parabolic profile path. A three-dimensional model of underwater towed system is studied. The established governing equations for the system have been solved using the central implicit finite-difference method. The obtained difference non-linear coupled equations are solved by Newton＇s method and satisfactory results were achieved. The solution of this problem has practical importance in the estimation of dynamic loading and motion, and hence it is directly applicable to the enhancement of safety and the effectiveness of the offshore activities.     

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.     

Mathematical modelling of response amplitude operator for roll motion of a floating body: Analysis in frequency domain with numerical validation

This paper investigates mathematical modelling of response amplitude operator （RAO） or transfer function using the frequency-based analysis for uncoupled roll motion of a floating body under the influence of small amplitude regular waves. The hydrodynamic coefficients are computed using strip theory formulation by integrating over the length of the floating body. Considering sinusoidal wave with frequency （ω ） varying between 0.3 rad/s and 1.2 rad/s acts on beam to the floating body for zero forward speed, analytical expressions of RAO in frequency domain is obtained. Using the normalization procedure and frequency based analysis, group based classifications are obtained and accordingly governing equations are formulated for each case. After applying the fourth order Runge-Kutta method numerical solutions are obtained and relative importance of the hydrodynamic coefficients is analyzed. To illustrate the roll amplitude effects numerical experiments have been carried out for a Panamax container ship under the action of sinusoidal wave with a fixed wave height. The effect of viscous damping on RAO is evaluated and the model is validated using convergence, consistency and stability analysis. This modelling approach could be useful to model floating body dynamics for higher degrees of freedom and to validate the result.     

Analysis and experiment of a vessels shell cover in submarine structure

This paper aims to achieve analysis and experiment results that relate to mechanics capability and structural parameter of a special saddle shell of revolution. Theoretically speaking, the saddle shell of revolution consists of a toroidal shell and a spherical shell. The shells simultaneous equations can he solved with harmonious terms. Where, the fundamental equations can be solved by asymptotic exponential perturbation method. The equations of special solution can be solved by Hovozhilovs special solution. This new idea is from a study of some existing solutions of the toroidal shell. The results have heen proved by eompared with some experimental results. The experiments aims to study the effect caused by change of material parameter,or by change of different geometrie dimensions of the saddle shell,which include the ehange of thickness, the change of radius of shell, and the change of ribs. Finally, the accepted product of the saddle shell were reinforeed by a toroidal rib has been submitted.     

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.     

双平面焊接圆管截面DKT节点的中心撑杆鞍点及冠点位置的应力集中系数(英文)

A set of parametric stress analyses was carried out for two-planar tubular DKT-joints under different axial loading conditions.The analysis results were used to present general remarks on the effects of the geometrical parameters on stress concentration factors(SCFs) at the inner saddle,outer saddle,and crown positions on the central brace.Based on results of finite element(FE) analysis and through nonlinear regression analysis,a new set of SCF parametric equations was established for fatigue design purposes.An assessment study of equations was conducted against the experimental data and original SCF database.The satisfaction of acceptance criteria proposed by the UK Department of Energy(UK DoE) was also checked.Results of parametric study showed that highly remarkable differences exist between the SCF values in a multi-planar DKT-joint and the corresponding SCFs in an equivalent uni-planar KT-joint having the same geometrical properties.It can be clearly concluded from this observation that using the equations proposed for uni-planar KT-connections to compute the SCFs in multi-planar DKT-joints will lead to either considerably under-predicting or over-predicting results.Hence,it is necessary to develop SCF formulae specially designed for multi-planar DKT-joints.Good results of equation assessment according to UK DoE acceptance criteria,high values of correlation coefficients,and the satisfactory agreement between the predictions of the proposed equations and the experimental data guarantee the accuracy of the equations.Therefore,the developed equations can be reliably used for fatigue design of offshore structures.     

Coupled motions of two ships in irregular waves in time domain

A three-dimensional time domain approach is used to study the coupled motions of two ships with forward speed in waves. In this approach, the boundary condition is satisfied on the mean wetted hull surface of the moving bodies and the free surface condition is linearized. The problem is solved by using a transient free-surface Green function source distribution on the submerged hulls. After solving the response amplitude operator, the method of spectral analysis is employed to clearly express the motion energy spectrum and significant amplitude of two ships. For verifying the code, two same circular cylinders at beam wave are selected to calculate coupled motions by comparison with the results obtained by 3D frequency method which has been proved to be efficient for solving such problems.Two Wigley ships of different sizes with the same forward speed are chosen for numerical calculation of the interaction effect, and some useful suggestions are obtained for underway replenishment at sea.     

An investigation into the isolation performance of mono-and bi-stable systems

Motivated by the need for improving the isolation performance, many research studies have been performed on isolators with nonlinear characteristics. Based on the shape of their phase portrait, such devices can be configured as either a mono-or bi-stable isolator. This paper focuses on investigating the relative performance of these two classes under the same excitations. Force transmissibility is used to measure the isolation performance, which is defined in terms of the RMS of the ratio of the transmitted force to the excitation force. When the system is subjected to harmonic excitation, it is found that the maximum reduction of the force transmissibility in the isolation range using Quasi-Zero stiffness is achieved. When the system is subjected to random excitation, it has the same effect of Quasi-Zero stiffness. Further, optimum damping can be changed with stiffness and has minimum value.     

Exciting forces for a wave energy device consisting of a pair of coaxial cylinders in water of finite depth

Two coaxial vertical cylinders-one is a riding hollow cylinder and the other a solid cylinder of greater radius at some distance above an impermeable horizontal bottom,were considered.This problem of diffraction by these two cylinders,which were considered as idealization of a buoy and a circular plate,can be considered as a wave energy device.The wave energy that is created and transferred by this device can be appropriately used in many applications in lieu of conventional energy.Method of separation of variables was used to obtain the analytical expressions for the diffracted potentials in four clearly identified regions.By applying the appropriate matching conditions along the three virtual boundaries between the regions,a system of linear equations was obtained,which was solved for the unknown coefficients.The potentials allowed us to obtain the exciting forces acting on both cylinders.Sets of exciting forces were obtained for different radii of the cylinders and for different gaps between the cylinders.It was observed that changes in radius and the gap had significant effect on the forces.It was found that mostly the exciting forces were significant only at lower frequencies.The exciting forces almost vanished at higher frequencies.The problem was also investigated for the base case of no plate arrangement,i.e.,the case having only the floating cylinder tethered to the sea-bed.Comparison of forces for both arrangements was carried out.In order to take care of the radiation of the cylinders due to surge motion,the corresponding added mass and the damping coefficients for both cylinders were also computed.All the results were depicted graphically and compared with available results.     

浮体横摇运动传递函数数学建模的频域分析研究及相关数值验证(英文)

This paper investigates mathematical modelling of response amplitude operator(RAO) or transfer function using the frequency-based analysis for uncoupled roll motion of a floating body under the influence of small amplitude regular waves. The hydrodynamic coefficients are computed using strip theory formulation by integrating over the length of the floating body. Considering sinusoidal wave with frequency( ω) varying between 0.3 rad/s and 1.2 rad/s acts on beam to the floating body for zero forward speed, analytical expressions of RAO in frequency domain is obtained. Using the normalization procedure and frequency based analysis, group based classifications are obtained and accordingly governing equations are formulated for each case. After applying the fourth order Runge-Kutta method numerical solutions are obtained and relative importance of the hydrodynamic coefficients is analyzed. To illustrate the roll amplitude effects numerical experiments have been carried out for a Panamax container ship under the action of sinusoidal wave with a fixed wave height. The effect of viscous damping on RAO is evaluated and the model is validated using convergence, consistency and stability analysis. This modelling approach could be useful to model floating body dynamics for higher degrees of freedom and to validate the result.     

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.     

自由衰减模型试验的横摇阻尼评估方法研究(英文)

The methodology to obtain the non-linear roll damping from decay tests is very old. It has been proposed by Froude in the 19th century and used from then on. Behind it there is a quadratic model [θ│θ│] for the damping and a subsequent equivalent linearization. Probably all model basin in the world follows this approach to assess the damping from a decay test. This is well documented and so is the methods to get the p1-p2 coefficients. This is very general in the sense that in principle,it could be applied t...     

基于改进的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.     

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.