1. Introduction
Large roll motions in parametric roll and dead ship conditions are serious risks for the safety of a ship in rough sea conditions. To predict the roll motions accurately at resonance condition, estimation of the roll damping is essential. However, the accurate prediction of a ship roll damping is difficult, except by means of high cost experiments. Numerical approaches like CFD are an alternative option to estimate roll damping by considering the viscous effect.
In general, most of the roll damping calculation methods are based on potential flow theory and empirical method. The most common empirical method is Ikeda’s method [
1]. Though this method can be used quite well for conventional ships, the prediction results are sometimes conservative or underestimated for unconventional ships [
2]. Roll damping is strongly nonlinear and is influenced by fluid viscosity and flow characteristics such as flow separation and vortex shedding. In theory, empirical or semi-empirical methods cannot take full consideration of different characteristics of a complex flow. Currently, vulnerability criteria for parametric roll and dead ship conditions are under development by the International Maritime Organization (IMO) as a second-generation intact stability criteria, in which roll damping coefficients are proposed, using Ikeda’s simplified method. The calculation for traditional ships by Ikeda’s simplified method can fit experimental data quite well at small roll angles. However, when the roll angle is large and out of the acceptable range of Ikeda’s method, the accuracy of the damping coefficient is low.
In addition to Ikeda’s simplified method, the correspondence group on Intact Stability regarding the second generation of intact stability criteria also proposed that the roll damping could be calculated by roll decay/forced roll test or CFD simulation. This suggestion can overcome the limitation of the model tests, which can predict roll damping very well but it is costly and time-consuming. Most of the experimental data is limited to a certain frequency range and particular geometry, which makes it impossible for the large-scale expansion of the application [
3].
For the accurate calculation of roll damping, the influence of viscosity must be considered. CFD numerical simulations can consider flow characteristics and also reduce the cost of experiments. By improving CFD technology, it is possible to estimate damping coefficients precisely. Over recent years, numerous research projects based on CFD and experimental simulations have been conducted, for instance, Roddier [
4] investigated two-dimensional simulations. The model was constrained to three degrees of freedom where it was free in roll, heave and sway motion. This numerical simulation used a random vortex method applied on a rectangular box in a beam wave condition. The results were validated against experimental data and showed that having bilge keels decreases roll resonance. They found that considering applied mechanical friction in the code can improve the accuracy of simulations. Na [
5] carried out an investigation on a rectangular box with and without bilge keels to question harmonic force roll motion using experimental simulations. Various bilge keels’ geometries were investigated with different width and angle to observe influences on the damping coefficient. As a result, it was found that bilge keels with larger lengths and horizontal orientation could improve the damping coefficient significantly.
Jung [
6] used a particle imaging velocimetry (PIV) method to perform experimental simulations of a box to analyse vortex and turbulence generation in roll motion under beam waves conditions. It was found that fixing the box increases the intensity of turbulence due to intensification of relative velocity around the box. The flow around the corners was more turbulent because of separation. Yi-Hsiang [
7] simulated the harmonic force roll motion of a floating production storage and offloading (FPSO) hull where the task was heavily involved with 2D CFD analysis. Bilge keels were attached in 45° inclination from the horizontal axis, which produced larger added moments of inertia and damping. The amplitude of added moments of inertia remained equal with horizontal or vertical bilge keels, while with the horizontal bilge keels it produced larger damping amplitude. Later, Kinnas [
8] utilized the incompressible Navier-Stokes solver to analyse a 2D simulation of an FPSO hull in harmonic force roll motion with and without bilge keel appendages. The results showed that in inviscid flow, there is a linear relation between roll moment and roll angle amplitude. However, a non-linear variation exists between the roll moment and roll angle in viscous flow condition and it was observed that by adding the bilge keel, the nonlinearity could increase. Wilson [
9] introduced numerical simulations based on unsteady Reynolds averaged Navier-Stokes (RANS) code to analyse the naval combatant’s motions and wave patterns. The numerical simulations were conducted with and without bilge keels to investigate harmonic excited roll motion. The outcomes showed a good correlation with experimental data especially in case of a hull with bilge keels. However, the numerical approach had difficulties for simulating the free surface in a large roll angle.
Yu [
10] utilized a 2D incompressible Navier-Stokes solver to investigate the roll damping of a rounded bilge box with and without bilge keels, as well as a sharp corner bilge box with and without a step at the keels. In that study, the exciting moment was subject to roll motion and it was observed that bilge keel increases the amplitude of damping and the relation between the roll moment and roll angle was nonlinear due to viscosity effects. Considering the roll damping of a rectangular barge, Bangun [
11] performed a numerical simulation using a 2D incompressible Navier-Stokes solver in order to investigate roll motion. The model was examined with various conditions considering with and without bilge keel, different width sizes and angle of bilge keels. In total, 12 cases were examined. It was concluded that the barge with a smaller bilge keel angle from the horizontal axis produced larger roll damping. Thiagarajan [
12] carried out experimental and numerical studies of FPSO in which the model was scaled 1:350. In numerical simulations, based on a free surface random vortex method (FSRVM), the model was forced to roll and results were in good agreement with the experimental data. Finally, it was concluded that the amplitude of damping is a function of roll angular velocity and width of bilge keels. An equation based on the relation between damping ratio and bilge keel width was proposed with some assumptions.
Avalos [
13] performed a numerical simulation to investigate roll decay. The results gained from the numerical simulations were in an acceptable range which correlates well with experimental results. Through the process, the size of the vortex was a function of roll motion amplitude and width of the bilge keel. The roll decay technique is generally not a preferable method to estimate roll-damping coefficients in large roll motions especially with forward speed. In case of roll decay method, the water is initially at rest and roll dissipation occurs in a transient condition. Instead, damping obtained from harmonic excited roll motion (HERM) technique is based on steady state rolling motions, where the initial transient has already completed and the system is undergoing harmonic periodic motions. Therefore, the uncertainty of results is lower than the decay test where the roll angle magnitude will be decreases quickly over one cycle especially in case with forward speed. Blume [
14] introduced a method to calculate roll damping coefficients effectively called the HERM technique, where this method excites the model in resonance frequency. However, this technique requires a longer time to determine the resonance frequency of the model. Another disadvantage of Blume’s method is the dependency of the roll damping coefficient to maximum roll angle, metacentric height and heel angle, where each of them can be subjected to errors. Handschel [
15] developed the HERM technique to estimate the damping coefficient in a range of frequencies that are very close to the resonance frequency. The technique considers the phase shift between the exciting moment and roll angles other than 90°. Begovic, Day [
16] carried out CFD simulations using STAR CCM+ to calculate the roll damping of the DTMB 5415 trough roll decay technique for both intact and damage conditions. The results are compared against experimental measurement with reasonable accuracy. Mancini, Begovic Mancini, Begovic [
17] conducted roll decay tests using numerical and experimental simulations to extract the roll damping coefficients. They considered the grid convergence index instead of correlation factor method to compute the uncertainty for the numerical simulations, because the solution was not close to asymptotic range. Zhou, Ning Zhou, Ning [
18] conducted numerical and experimental simulations to estimate the roll damping of four different types of ships based on roll decay technique in zero forward speed. The results from experiments and numerical simulations were in good agreement. Somayajula and Falzarano [
19] developed an advanced system of identification to compute frequency dependent roll damping from model test results in irregular waves. The results showed that the method can be used to predict a ship roll motion accurately compared to the potential flow method and empirical methods. Irkal and Nallayarasu Irkal, Nallayarasu [
20] performed experimental and numerical simulations to compute the impact of bilge keels on roll damping. PIV method was used for the experiment to measure velocity field around the model during free oscillation tests. They found that the roll damping coefficient of the model without bilge keels is linear, however, with the bilge keels is strongly non-linear. Wassermann and Feder Wassermann, Feder [
21] carried out model tests based on roll decay and HERM technique to calculate the roll damping of a container ship. They proposed various methods without additional filtering, curve fitting and offset manipulation of the recorded time series. They found that the HERM technique is more reliable in cases with higher forward speed and larger damping values. Oliva-Remola and Bulian Oliva-Remola, Bulian [
22] conducted the HERM technique by shifting a mass harmonically inside the model in the lateral direction to generate excitation. The computed roll damping from HERM technique was smaller than roll decay tests for the same roll angle, because the model reached to a steady state rolling in HERM technique whereas the roll dissipation of roll decay tests occurs in a transient condition. They also proposed a 1 DOF mathematical model to predict the roll motion and calculate the roll damping. It was observed that tuning of dry roll inertia is critical to achieve good results, because the model was considered free in roll and sway motion.
There are a limited number of studies regarding the roll damping coefficients of an entire model and most studies considered a segment of the ship especially the middle section. The focus of numerical simulations was on 2D and overlooked the effect of other motions such as pitch and longitudinal turbulence in the case with forward speed. Therefore, the impact of different degrees of freedom (DOF) on the roll damping coefficient is unknown. In the present study, the numerical simulations of a whole containership model are conducted based on the HERM technique. The impact of forward speed, DOF and excitation frequency at beam sea and oblique sea conditions on roll motion characteristics and roll damping coefficients are investigated.
2. Theoretical Background
The following discussion involves the ‘harmonic excited roll motion’ HERM technique to determine required roll damping coefficients. Regarding Blume’s experimental setup [
14], the model possesses two masses at the centre of gravity that are rotating contrarily around the vertical axis. Specifically, one of the masses is rotating in a clockwise and the other is rotating in an anti-clockwise direction. The rotating mass shares the same frequency in the opposite direction to minimise the yaw motion. During contrary motion, the two masses meet at both sides of the model twice per rotation period, which imposes the maximum roll excitation moment. Various roll amplitudes can be achieved by setting different weights of the rotational masses. From the experiments [
15], the amplitude of the roll exciting moment should ideally be equal to the restoring moment of the heel angle which refers to two masses of rotation on one side of the model. Originally, the equation of roll motion was formulated by using Newton’s second law is balanced between the ships’ motions and external moments as [
23]:
whereas,
| Mass and added mass moment of inertia coefficients |
| Damping moment coefficient |
| Restoring coefficient |
| Roll excitation moment |
The magnitude of damping moment in roll motion is generally less than the total moment of inertia and restoring moment. However, it still is essential to estimate the damping moment because inertia and restoring moments could be shifted in 180 degrees and counteract each other. As a result, only the damping moment limits the roll motion [
24]. The exciting moment changes harmonically:
The induced roll angle and roll velocity are as below:
and
The dissipated energy during a harmonic roll period is:
Substituting the Equations (3) and (4) into Equation (5) with phase angle (𝜗 = 0) gives:
Allowing the general solution of the integral gives:
The solution to define the dissipated damping energy over a cycle is:
With respect to the roll moment, the roll angle is phase shifted by 𝜗. The following formula provides an estimation for the work done by the exciting moment in one roll period:
Substituting Equation (4) into Equation (9) and after solving integral form, the formula gives:
The work done by the exciting moment and the dissipated energy over one roll period should relatively be the same. There are a couple of methods available to determine the roll damping coefficient, however, this paper focuses on the relationship of 𝐸
𝐴 = 𝐸
𝐸 where the roll damping can be calculated by the following formula:
The damping coefficient can now be established to a dimensionless form, which was recommended by the ITTC [
25].