Zero Speed Rankine-Kelvin Hybrid Method With a Cylinder Control Surface

Author(s):  
Hui Li ◽  
Hao Lizhu ◽  
Huilong Ren ◽  
Xiaobo Chen

The solution of hydrodynamic problem with forward speed still has some well-known problems such as high oscillation and slow convergence of the wave term when using a moving and oscillating source as the Green function. Recently, Ten and Chen (2010) has come up with a new method to benefit the merits of both the Rankine source and moving and oscillating source by taking a hemisphere as the control surface which separates the fluid region into two domains, but some troubles have been induced in the process of solution. Therefore, in this paper, a cylindrical surface instead of a hemisphere is selected to be the control surface to make the solution easy, and in this method, the control surface isn’t divided into panels. In the interior domain near the ship, the Rankin Green function is used to simplify the calculation. In the exterior domain some distance from the ship, there is no panels representing the free surface by using the Green function which satisfy the free surface boundary condition. The whole fluid region matches by the condition that the velocity potentials and their normal derivatives in the interior domain and exterior domain are equal on the control surface separately. In this paper, we have validated the Rankine-Kelvin hybrid method is applicable by adopting it to solve the zero speed problem in this work.

Author(s):  
Hui Li ◽  
Baoli Deng ◽  
Chunlei Liu ◽  
Jian Zou ◽  
Huilong Ren

Abstract A novel matching method has been developed to solve the wave-induced motions and loads of ships with forward speed. The fluid domain is divided into two subdomains by a cylindrical control surface: an interior domain and an exterior domain. Unlike the conventional domain decomposition strategy, the control surface is meshless in present method, on which the physical quantities are expanded into Fourier-Laguerre series. Based on forward speed Green function, the source distribution method is adopted to solve the exterior domain. The calculations of boundary integral equation about forward speed Green function over the control surface are performed analytically, and the solution of exterior domain provides a Dirichlet-to-Neumann (DN) relation on the control surface. In the interior domain, the boundary value problem is solved by Rankine source method. In order to be consistent with exterior solution, the control surface is kept meshless. The ship hull is discretized into constant panels. The free-surface is discretized into cubic B-splines to represent the high-order derivatives of velocity potential precisely. Then, the DN relation is used to close the equation system established in the interior domain. Comparisons with known experimental measurements show that the calculations achieve good accuracy. Furthermore, the influences of numerical method used in the exterior domain are described.


Author(s):  
Zhiyuan Pan ◽  
Torgeir Vada ◽  
Kaijia Han

A time domain Rankine source solver is extended to compute the wave added resistance of ships. The proposed approach applies the momentum conservation principle on the near field fluid volume enclosed by the wet surface of a floating body, the free surface and a control surface. The wave added resistance is then calculated by the integration over the control surface of the fluid velocities and free surface elevations. To be able to incorporate the proposed method with the Rankine source code, an interpolation scheme has been developed to compute the kinematics for the off-body points close to (or on) the free surface. Two Wigley ship models, a containership model S175 and a tanker model KVLCC2 are used to validate the present method. In general good agreement is found comparing with the model test data. The convergence behavior is examined for the proposed method including the selection of the time step and location of the control surface. Both Neumann-Kelvin and double body linearization methods are evaluated with the proposed method. It is found that the Neumann-Kelvin linearization can only be applied for slender ship hull, whereas double body method fits also for blunt ships. It is suggested to apply the proposed method with double body linearization to evaluate the wave added resistance of ships with a control surface close to the ship hull.


2009 ◽  
Vol 53 (03) ◽  
pp. 137-150
Author(s):  
Francis Noblesse ◽  
Gérard Delhommeau ◽  
Chi Yang

The linearized potential flow resulting from a distribution of pressure that advances at constant speed along a straight path at the free surface of calm water, of effectively infinite depth and lateral extent, is considered. A practical method for evaluating the free-surface elevation caused by the moving free-surface pressure patch—which can be used to model steady flows of air-cushion vehicles, high-speed planing boats, surface-effect ships, and some types of hybrid ships—is given. The key ingredient of this method is a highly simplified analytical approximation to the local-flow component in the expression for the Green function associated with the classic Michell-Kelvin linearized free-surface boundary condition.


1993 ◽  
Author(s):  
Bruce S. Rosen ◽  
Joseph P. Laiosa ◽  
Warren H. Davis ◽  
David Stavetski

A unique free-surface flow methodology and its application to design and analysis of IACC yachts are discussed. Numerical aspects of the inviscid panel code and details of the free-surface boundary condition are included, along with enhancements developed specifically for the '92 America's Cup defense. Extensive code validation using wind tunnel and towing tank experimental data address several areas of interest to the yacht designer. Lift and induced drag at zero Froude number are studied via a series of isolated fin/bulb/winglet appendages. An isolated surface piercing foil is used to evaluate simple lift/free­surface interactions. For complete IACC yacht models, upright wave resistance is investigated, as well as lift and induced drag at heel and yaw. The excellent correlation obtained for these cases demonstrates the value of this linear free-surface methodology for use in designing high performance sailing yachts.


2018 ◽  
Vol 845 ◽  
pp. 682-712 ◽  
Author(s):  
Zhi Fu Li ◽  
Guo Xiong Wu ◽  
Chun Yan Ji

Wave radiation and diffraction by a circular cylinder submerged below an ice sheet with a crack are considered based on the linearized velocity potential theory together with multipole expansion. The solution starts from the potential due to a single source, or the Green function satisfying both the ice sheet condition and the crack condition, as well as all other conditions apart from that on the body surface. This is obtained in an integral form through Fourier transform, in contrast to what has been obtained previously in which the Green function is in the series form based on the method of matched eigenfunction expansion in each domain on both sides of the crack. The multipole expansion is then constructed through direct differentiation of the Green function with respect to the source position, rather than treating each multipole as a separate problem. The use of the Green function enables the problem of wave diffraction by the crack in the absence of the body to be solved directly. For the circular cylinder, wave radiation and diffraction problems are solved by applying the body surface boundary condition to the multipole expansion, through which the unknown coefficients are obtained. Extensive results are provided for the added mass and damping coefficient as well as the exciting force. When the cylinder is away from the crack, a wide spacing approximation method is used, which is found to provide accurate results apart from when the cylinder is quite close to the crack.


Author(s):  
Kevin Markeng ◽  
Torgeir Vada ◽  
Zhi Yuan Pan

In this paper two methods for modelling the damping in a narrow gap are investigated. The first method is called the Pressure Damping Model. This method has been used in studies of wave energy devices. An attractive feature of this model is that the modified input is directly related to the energy dissipation in the gap, which means that if this dissipation is estimated the input to the model can be obtained directly. The idea of the method is to add a pressure input in the gap to suppress the resonant motion. A challenge with the method is that it contains a non-linear term. The second method is the Newtonian Cooling damping model. The method is based on introducing a dissipation term in the free surface boundary condition. This dissipation term contains a coefficient which is not directly related to the energy dissipation. Hence this method is not so easy to relate directly to the estimated energy dissipation. An advantage with this method is that it is linear and hence can be expected to be more robust. In the first part of the paper a 2-dimensional problem is investigated using both methods. In addition to the numerical performance and robustness, much focus is put on investigation of the energy balance in the solution, and we attempt to relate both models to the energy dissipation in the gap. In the second part the Newtonian cooling method is implemented in a 3-dimensional potential flow solver and it is shown that the method provides a robust way to handle the resonance problem. The method will give rise to a modified set of equations which are described. Two different problems are investigated with the 3D solver. First we look at a side-by-side problem, where the 3D results are also compared with 2D results. Finally, the moonpool problem is investigated by two different 3D solvers, a classical Green’s function based method and a Rankine solver. It is also shown how this damping model can be combined with a similar model on the internal waterplane to remove irregular frequencies.


Author(s):  
M. D. Yang ◽  
B. Teng

A time-domain simulation method is developed for the coupled dynamic analysis of a spar platform with mooring lines. For the hydrodynamic loads, a time domain second order method is developed. In this approach, Taylor series expansions are applied to the body surface boundary condition and the free surface boundary condition, and Stokes perturbation procedure is then used to establish corresponding boundary value problems with time-independent boundaries. A higher order boundary element method is developed to calculate the velocity potential of the resulting flow field at each time step. The free-surface boundary condition is satisfied to the second order by 4th order Adams-Bashforth-Moultn method. An artificial damping layer is adopted on the free surface to avoid the wave reflection. For the mooring-line dynamics, a geometrically nonlinear finite element method using isoparametric cable element based on the total Lagrangian formulation is developed. In the coupled dynamic analysis, the motion equation for the hull and dynamic equations for mooring lines are solved simultaneously using Newmark method. Numerical results including motions and tensions in the mooring lines are presented.


Author(s):  
Yingyi Liu ◽  
Ying Gou ◽  
Bin Teng

The present study aims to develop an efficient numerical method for computing the diffraction and radiation of water waves with horizontal long cylindrical structures, such as floating breakwaters. A higher-order scheme is used to discretize geometry of the structure as well as the relevant physical quantities. As the kernel of this method, Wehausen’s free-surface Green function is calculated by a newly-developed Gauss-Kronrod adaptive quadrature algorithm after elimination of its Cauchy-type singularities. To improve computational efficiency, a Chebyshev approximation approach is applied to a fast calculation of the Green function that needs evaluation thousands of times. In addition, OpenMP parallel technique is used to the formation of influence coefficient matrix, which significantly reduces CPU time. Finally, computations are performed on wave exciting forces and hydrodynamic coefficients for the long cylindrical structures, either floating or submerged. Comparison with other numerical and analytical methods demonstrates good performance of the present method.


2019 ◽  
Vol 63 (4) ◽  
pp. 251-267 ◽  
Author(s):  
Zhi-Ming Yuan ◽  
Liang Li ◽  
Ronald W. Yeung

Ships often have to pass each other in proximity in harbor areas and waterways in dense shipping-traffic environment. Hydrodynamic interaction occurs when a ship is overtaking (or being overtaken) or encountering other ships. Such an interactive effect could be magnified in confined waterways, e.g., shallow and narrow rivers. Since Yeung published his initial work on ship interaction in shallow water, progress on unsteady interaction among multiple ships has been slow, though steady, over the following decades. With some exceptions, nearly all the published studies on ship-to-ship problem neglected free-surface effects, and a rigid-wall condition has often been applied on the water surface as the boundary condition. When the speed of the ships is low, this assumption is reasonably accurate as the hydrodynamic interaction is mainly induced by near-field disturbances. However, in many maneuvering operations, the encountering or overtaking speeds are actually moderately high (Froude number Fn > 0.2, where <inline-graphic xlink:href="josr10180089inf1.tif"/>, U is ship speed, g is the gravitational acceleration, and L is the ship length), especially when the lateral separation between ships is the order of ship length. Here, the far-field effects arising from ship waves can be important. The hydrodynamic interaction model must take into account the surface-wave effects. Classical potential-flow formulation is only able to deal with the boundary value problem when there is only one speed involved in the free-surface boundary condition. For multiple ships traveling with different speeds, it is not possible to express the free-surface boundary condition by a single velocity potential. Instead, a superposition method can be applied to account for the velocity field induced by each vessel with its own and unique speed. The main objective of the present article is to propose a rational superposition method to handle the unsteady free-surface boundary condition containing two or more speed terms, and validate its feasibility in predicting the hydrodynamic behavior in ship encountering. The methodology used in the present article is a three-dimensional boundary-element method based on a Rankine-type (infinite-space) source function, initially introduced by Bai and Yeung. The numerical simulations are conducted by using an in-house‐developed multibody hydrodynamic interaction program “MHydro.” Waves generated and forces (or moments) are calculated when ships are encountering or passing each other. Published model-test results are used to validate our calculations, and very good agreement has been observed. The numerical results show that free-surface effects need to be taken into account for Fn > 0.2.


Sign in / Sign up

Export Citation Format

Share Document