DTMB5512 surface combatant in regular head waves

This benchmark corresponds to the simulation of the surface combatant DTMB5512 with constant forward speed submitted to regular head waves. The DTMB5512 is a 3.048m geosym of the DTMB5415 which has been chosen by the International Towing Tank Conference (ITTC) as an international benchmark for CFD validation. Two configurations are considered, captive motion and free motion in pitch and heave (3DOF). Description of this test case and comparison to experimental results are summarized in the following.

Description of the test case

In this simulation, the 1:46.6 scale is considered which corresponds to the laboratory scale of the DTMB 5512 geometry used in the experimental campaigned made by the Iowa Institute of Hydraulic Research (IIHR) [Gui2001] [Irvine2008]. The hull model with no appendage, no propeller and no rudder is considered. The CAD model and data are freely available on the IIHR website. The main characteristics of the vessel are summarized in the following table. The x-axis is pointing forward to the vessel and z-axis is pointing upward as represented in the following figure.

DTMB5512 geometry

Fig. 36 Description of the test case.

Parameter Unit Value
Scale
46.6
Length (Lpp) m 3.048
Beam m 0.405
Draft m 0.132
Wetted surface area \(m^2\) 1.371
LCG m 1.536
VCG m 0.162

Hydrostatic coefficients are summarized in the next table.

Hydrostatic coeff. Value \((x10^3)\)
k33 \(9.68\)
k55 \(5.42\)
k35 \(1.25\)

The DTMB5512 is mounted on a carriage with constant forward speed as represented in next figure. A video is also available for the semi captive test.

DTMB5512 simulation

For the captive test , the DTMB5512 model is fixed on the carriage whereas for pitch and heave motion test the DTMB5512 model is free to move in heave and pitch. Four different speed are considered which correspond to a low (0.19), medium (0.28), mid-high (0.34) and high (0.41) froude numbers.

Integration of steady forces due to forward speed

Resistance force (ITTC57)

Due to the forward speed the vessel is submitted to a resistance force in surge motion and additional heave and pitch steady force. These forces are determined from the experimental results in steady motion with no waves.

Resistance force is computed according to the ITTC57 standard with a hull factor \(k=0.03\). Based on experimental results, a residual coefficient \(C_R\) is estimated for the different vessel speed.

Froude Number Residual coeff.
0.19 \(5,37.10^{-4}\)
0.24 \(9,07.10^{-4}\)
0.34 \(1,68.10^{-3}\)
0.41 \(4,02.10^{-3}\)

User defined forces

Based on the experimental results with forward speed and no waves, steady forces for pitch \(M_{y,user}\) and \(F_{z,user}\) are defined:

\[\begin{split}F_{z,user} = -12.3 u^3 - 2.9 u^2 \\ M_{y,user} = 4.3 u^6 - 9.1 u^5 - 9.7 u^4 + 34.2 u^3 - 22.7 u^2\end{split}\]

where \(u\) is the steady forward speed of the vessel.

Wave drift force

For regular wave, the wave drift force is estimated from the following formula:

\[F_{wd} = A^2 . C_{wd}(\omega)\]

where \(A\) is the wave amplitude, \(\omega\) is the wave frequency and \(C_{wd}\) the mean wave drift coefficient.

In this test case, the mean wave drift coefficient \(C_{wd}\) is estimated from experimental results and added resistance given by :

\[C_{T,ad} = \bar{C_T} - C_{T,st}\]

where \(\bar{C_T}\) is the mean unsteady resistance and \(C_{T,st}\) is the steady resistance in calm water.

Estimation of the mean wave drift coefficient in surge \(C_{wd, x}\) for various wave frequencies is represented in Fig. 38

Mean Wave drift

Fig. 38 Estimation of the mean wave drift coefficient depending on wave frequencies.

Adimentionalization of the forces

To compare forces and moment applied on the vessel to the experimental results, the following adimentionalization is applied:

\[\begin{split}C_t = \frac{F_x}{0.5 \rho U^2 S} \\ C_h = \frac{F_z}{0.5 \rho U^2 S} \\ C_m = \frac{M_y}{0.5 \rho U^2 L S}\end{split}\]

Following [ref], harmonic decomposition of the forces and moment are applied as follows:

\[X_F(t) = \frac{X_0}{2} + \sum_n X_n cos(2 \pi n f_e t)\]

where \(X_0\) is the constant part and \(X_n\) the \(n^{th}\) harmonic coefficient of the function \(X_F\).

Captive test results

The regular wave field propagates in negative x-direction. Four different wave amplitudes are considered corresponding to small (0.025), small-median (0.05, 0.075) and median (0.1) steepness. The wave period is equal to 0.22 seconds.

The zeroth and first harmonic coefficients of \(C_T\), \(C_H\) and \(C_M\) are compared to experimental results [Gui2002] in Fig. 39 and Fig. 40.

Zeroth coefficient

Fig. 39 Comparison of the zeroth coefficients from FRyDoM (continuous lines) with experimental results (discontinuous lines).

First coefficient

Fig. 40 Comparison of the first harmonic coefficient from FRyDoM (continuous lines) with experimental results (discontinuous lines).

Pitch and Heave Motion test results

The Response Amplitude Operator (RAO) of the DTMB5512 in pitch and heave motion are compared to experimental results. Five different forward speed are considered corresponding to different value of the Froude number equal to 0., 0.19, 0.28, 0.34 and 0.41. The numerical results agree well with the experimental results for heave motion. The natural period of the model is around 1 Hz and the peak of amplitude increase with increasing speed. The results in pitch motion agree well with the experimental results in the case of no forward speed. Discrepancy are more pronounced for the case with forward speed. Work are still to be done for this test case in order to characterize these discrepancies.

Heave motion

Fig. 41 Comparison of the Heave RAO results from FRyDoM with experimental data for different regular wave frequency. Results for different Froude number (0.0, 0.19, 0.28, 0.34, 0.41) from left to right and top to bottom.

Pitch motion

Fig. 42 Comparison of the Pitch RAO results from FRyDoM with experimental data for different regular wave frequency. Results for different Froude number (0.0, 0.19, 0.28, 0.34, 0.41) from left to right and top to bottom.

References

[Gui2001]
  1. Gui, J. Longo, B. Metcalf, J. Shao, F. Stern, “Forces, moment, and wave pattern for surface combatant in regular head waves. Part 1 : Measurment systems and uncertainty assessment”, Experiments in Fluids, Vol 31, 2001, pp 674-680.
[Gui2002]
  1. Gui, J. Longo, B. Metcalf, J. Shao, F. Stern, “Forces, moment, and wave pattern for surface combatant in regular head waves. Part 2 : Measurment results and discussions”, Experiments in Fluids, Vol 32, 2002, pp 27-36.
[Irvine2008]
  1. Irvine, J. Longo, F. Stern, “Pitch and Heave Tets Uncertainty assessment for a surface combatant in regular head waves”, Journal Ship Research, Vol 52, No 2, June 2008, pp 146-163.