Computation of Wave Added Resistance by Control Surface Integration

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.

Author(s):  
Heinrich Söding

A 3-dimensional Rankine source panel method for simulating a rigid floating body in steep waves is being developed. The aim is to obtain the same quality as free-surface RANSE methods, which are well suited for this application, but to require only a small fraction of the computing time needed by RANSE methods. The body may have forward speed or perform maneuvering motions. The exact boundary conditions are satisfied at the actual location of the fluid boundaries. The waves are generated not by a material wave maker, but by an approximate wave potential which needs not satisfy the exact free-surface condition. No wave damping regions are required. Whereas for steep waves without a body the method appears satisfactory, it needs further improvements if a body is present.


Author(s):  
B. Padmanabhan ◽  
R. C. Ertekin

This work is motivated by the many instances of intake/discharge flows from openings on floating or submerged ocean vessels and structures that may affect the wave field around them. Damaged vessels may release oil, or water may enter these vessels through openings. In oil skimming operations, for example, a very thin layer of oil must be skimmed off a large surface area, and therefore, oil skimming vessels require large intakes. Floating OTEC plants also require large intake and discharge volumes to sustain their operations. A linear theory is developed to obtain the motions of a 2-dimensional, freely floating body (from which steady intake/discharge flows originate) that encounters incoming waves. The boundary-value problem is formulated within the assumptions of linear potential theory by decomposing the total potential into its oscillatory and steady components. The steady potential is further decomposed into the double-model and perturbation potentials. The time-harmonic potential is coupled with the steady potential through the free-surface condition. The potentials are obtained by use of the quadratic boundary-element method based on the Rankine source. The effect of the steady intake/discharge flows on the diffraction loads, hydrodynamic force coefficients, as well as the motions of a 2-dimensional prismatic body floating on the free surface are presented. It is shown that the exciting wave forces and the hydrodynamic coefficients other than the damping coefficients are not appreciably affected by the intake/discharge flows of low Froude number for a 100MW floating OTEC plant.


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):  
Boris Horel ◽  
Pierre-Emmanuel Guillerm ◽  
Jean-Marc Rousset ◽  
Bertrand Alessandrini

The modeling of ship behavior in astern seas requires a large range of maneuverability and seakeeping knowledge since the understanding of the ship motions returns to solve a fluid structure interactions problem between waves and the ship hull. The broaching phenomenon is known as an abrupt change in motion in the horizontal plane, resulting in a loss of ship’s heading. It is characterized by a sudden divergence of yaw. Thus, there is a transfer of the kinetic energy on the roll axis that increases the risk of ship capsize. In the aim of modeling this phenomenon, the developed model uses the capture of the intersection between the ship hull and the free surface. Thus, we can overcome the hydrostatic stiffness matrix and integrate directly the hydrostatic pressure on the immersed surface. This method has the advantage of taking into account non-linearities of the wave profile into the calculation of the immersed surface, directly by performing a remodeling of the facets near the free surface. In the literature, three main factors are likely to affect the stability: the loading of the vessel, the presence of external disturbance torques and inadequate conditions of navigation, as is the case when a ship is caught in a storm. The first two factors are taken into account in the study of static stability, while the third factor is considered in the study of the instantaneous stability. Hydrostatic behavior of a ship is interesting when one wants to know her intact stability limits in calm seas. However, in the study of the ship behavior in following seas, the ship is no longer in usual conditions of navigation, but in unsuitable conditions requiring the study of the instantaneous stability. In the model formulation, the dynamic torsor comes from the general non-linear maneuverability equations and the time advance is solved by a 4th order Runge Kutta scheme with a constant time step. The torsor of the total applied mechanical action on the ship hull is expressed as the superposition of six torsors (gravity, hydrostatic, Froude Krylov, radiation, hydrodynamics and maneuverability) expressed in the center of gravity of the ship. Thus, we obtain a strong coupling between the maneuverability and seakeeping equations. Validation cases will be conducted and presented. The improvement of the model will require the implementation of test campaigns that will be specific for the study of ship behavior in astern seas. Validation of the model will help to define new stability criteria for ships in wave.


2003 ◽  
Vol 125 (1) ◽  
pp. 41-47 ◽  
Author(s):  
B. Padmanabhan ◽  
R. C. Ertekin

A linear theory is developed to obtain the motions of a two-dimensional, freely floating body (from which steady intake/discharge flows originate) that encounters incoming waves. The boundary-value problem is formulated within the assumptions of linear potential theory by decomposing the total potential into its oscillatory and steady components. The steady potential is further decomposed into the double-model and perturbation potentials. The time-harmonic potential is coupled with the steady potential through the free-surface condition. The potentials are obtained by use of the quadratic boundary-element method based on the Rankine source. The effect of the steady intake/discharge flows on the diffraction loads, hydrodynamic force coefficients, as well as the motions of a two-dimensional prismatic body floating on the free surface are presented. It is shown that the exciting wave forces and the hydrodynamic coefficients other than the damping coefficients are not appreciably affected in the case of low intake/discharge Froude numbers that are estimated, for example, for a 100 MW floating OTEC plant.


2011 ◽  
Vol 55 (03) ◽  
pp. 163-184 ◽  
Author(s):  
Fuat Kara

The prediction of the added resistance of the ships that can be computed from quadratic product of the first-order quantities is presented using the near-field method based on the direct pressure integration over floating body in time domain. The transient wave-body interaction of the first-order radiation and diffraction problems are solved as the impulsive velocity of the floating body by the use of a three dimensional panel method with Neumann-Kelvin method. These radiation and diffraction forces are the input for the solution of the equation of the motion that is solved by the use of the time marching scheme. The exact initial-boundary-value problem is linearized about a uniform flow, and recast as an integral equation using the transient free-surface Green function. A Wigley III hull form with forward speed is used for the numerical prediction of the different parameters. The calculated mean second-order added resistance and unsteady first-order impulse-response functions, hydrodynamics coefficients, exciting forces, and response amplitude operators are compared with experimental results.


2013 ◽  
Vol 57 (01) ◽  
pp. 13-23
Author(s):  
Wei Qiu ◽  
Hongxuan (Heather) Peng

Motions of a floating body in waves are computed in the time domain by solving the body-exact problem with the panel-free method and exact geometry. In the present study, the body boundary condition is imposed on the instantaneous wetted surface exactly at each time step. The free surface boundary is assumed linear so that the time-domain Green function can be applied. The body geometry is represented by NonUniform Rational B-Spline surfaces. At each time step, the instantaneous wetted surface is obtained by trimming the entire body surface. With the panel-free method, the body-exact problems are solved without involving repanelization of the wetted hull surface at each time step. Validation studies have been carried out for a submerged sphere, a flared body, and a Wigley hull. The hydrodynamic forces on the submerged sphere undergoing large-amplitude motion were computed and compared with analytical solutions. For the flared body oscillating in a free surface and the Wigley hull in waves, numerical results were compared with experimental data and solutions by other numerical methods.


2007 ◽  
Author(s):  
Adrian Lungu ◽  
Theodore E. Simos ◽  
George Psihoyios ◽  
Ch. Tsitouras

Author(s):  
Shanti Bhushan ◽  
Pablo Carrica ◽  
Jianming Yang ◽  
Frederick Stern

Scalability studies and computations using the largest grids to date for free-surface flows are performed using message-passing interface (MPI)-based CFDShip-Iowa toolbox curvilinear (V4) and Cartesian (V6) grid solvers on Navy high-performance computing systems. Both solvers show good strong scalability up to 2048 processors, with V6 showing somewhat better performance than V4. V6 also outperforms V4 in terms of the memory requirements and central processing unit (CPU) time per time-step per grid point. The explicit solvers show better scalability than the implicit solvers, but the latter allows larger time-step sizes, resulting in a lower total CPU time. The multi-grid HYPRE solver shows better scalability than the portable, extensible toolkit for scientific computation solver. The main scalability bottleneck is identified to be the pressure Poisson solver. The memory bandwidth test suggests that further scalability improvements could be obtained by using hybrid MPI/open multi-processing (OpenMP) parallelization. V4-detached eddy simulation (DES) on a 300 M grid for the surface combatant model DTMB 5415 in the straight-ahead condition provides a plausible description of the vortical structures and mean flow patterns observed in the experiments. However, the vortex strengths are over predicted and the turbulence is not resolved. V4-DESs on up to 250 M grids for DTMB 5415 at 20° static drift angle significantly improve the forces and moment predictions compared to the coarse grid unsteady Reynolds averaged Navier–Stokes, due to the improved resolved turbulence predictions. The simulations provide detailed resolution of the free-surface and breaking pattern and vortical and turbulent structures, which will guide planned experiments. V6 simulations on up to 276 M grids for DTMB 5415 in the straight-ahead condition predict diffused vortical structures due to poor wall-layer predictions. This could be due to the limitations of the wall-function implementation for the immersed boundary method.


2021 ◽  
Vol 153 (A2) ◽  
Author(s):  
Q Yang ◽  
W Qiu

Slamming forces on 2D and 3D bodies have been computed based on a CIP method. The highly nonlinear water entry problem governed by the Navier-Stokes equations was solved by a CIP based finite difference method on a fixed Cartesian grid. In the computation, a compact upwind scheme was employed for the advection calculations and a pressure-based algorithm was applied to treat the multiple phases. The free surface and the body boundaries were captured using density functions. For the pressure calculation, a Poisson-type equation was solved at each time step by the conjugate gradient iterative method. Validation studies were carried out for 2D wedges with various deadrise angles ranging from 0 to 60 degrees at constant vertical velocity. In the cases of wedges with small deadrise angles, the compressibility of air between the bottom of the wedge and the free surface was modelled. Studies were also extended to 3D bodies, such as a sphere, a cylinder and a catamaran, entering calm water. Computed pressures, free surface elevations and hydrodynamic forces were compared with experimental data and the numerical solutions by other methods.


Sign in / Sign up

Export Citation Format

Share Document