Forward-Speed Ship Wave Flows

From WikiWaves
Jump to: navigation, search

Linear free surface conditions

The linear free surface conditions are as follows,

[math] \frac{\mathrm{d}^2\Phi}{\mathrm{d}t^2} + g \frac{\mathrm{d}\Phi}{\mathrm{d}z} = 0, \quad z=0 [/math]

We can calculate the second derivative above as

[math] \frac{\mathrm{d}^2\Phi}{\mathrm{d}t^2} = \frac{\partial^2 \bar{\Phi}}{\partial t^2} - 2 U \frac{\partial^2\bar{\Phi}}{\partial x \partial t} + U^2 \frac{\partial \bar{\Phi}}{\partial x^2} [/math]

When no ambient waves are present [math] \frac{\partial\bar{\Phi}}{\partial t} = 0 \, [/math] and we obtain the free surface condition for the steady Kelvin ship wave problem

[math] U^2 \frac{\partial^2 \bar{\Phi}}{\partial x^2} + g \frac{\partial\bar{\Phi}}{\partial z} = 0, \qquad z=0 [/math]

This is the famous Neumann-Kelvin free-surface condition governing the linear steady wave pattern generated by a translating ship.

Bernoulli equation:

[math] P = - \rho \left( \frac{\mathrm{d}\Phi}{\mathrm{d}t} + \frac{1}{2} \nabla\Phi \cdot \nabla\Phi + gz \right) = - \rho \left[ \left( \frac{\partial\bar{\Phi}}{\partial t} - U \frac{\partial\bar{\Phi}}{\partial x} \right) + \frac{1}{2} \nabla \phi \cdot \nabla \phi + g z \right] [/math]

Free surface elevation:

[math] \zeta = - \left. \frac{1}{g} \frac{\mathrm{d}\Phi}{\mathrm{d}t} \right|_{z=0} = -\frac{1}{g} \left. \left( \frac{\partial\bar{\Phi}}{\partial t} - U \frac{\partial\bar{\Phi}}{\partial x} \right)\right|_{z=0} [/math]

In the Kelvin ship wave problem

[math] \frac{\partial\phi}{\partial t} = 0 [/math]


[math] \zeta = \frac{U}{g} \frac{\partial\bar{\Phi}}{\partial x}, \qquad z=0 [/math]

Neumann-Kelvin linearization of [math]U\gt0[/math] ship seakeeping problem

Let [math] \Phi(X_0,Y_0,Z_0,t) \,[/math]be the total potential relative to the inertial coordinate system

[math] X_0 = x + U t \, [/math]

Let [math] \Phi( X,Y,Z,t) \, [/math] be the same potential expressed relative to the translating frame. It was shown before that

[math] \frac{\mathrm{d}\Phi}{\mathrm{d}t} = \left( \frac{\partial}{\partial t} - U \frac{\partial}{\partial X} \right) \phi (\vec{X}, t) \, [/math]

where now the time dependence of [math]\phi\,[/math] w.r.t. time will be the [math] e^{i\omega t} \,[/math] form in the ship seakeeping problem.

The total potential [math] \Phi\,[/math] consists of the sum of two components in a linearized setting

[math] \mathbf{\Phi}_{TOT} = \mathbf{\bar{\Phi}} + \mathbf{\Phi} \,[/math]

where [math] \mathbf{\bar{\Phi}}\,[/math] is the velocity potential due to the vessel forward translation with constant speed [math] U\,[/math] and [math] \Phi\,[/math] is the seakeeping component due to vessel motions in waves.

Relative to the ship frame:

[math] \mathbf{\bar{\Phi}} = \bar{\Phi} (X,Y,Z) \, [/math]
[math] \mathbf{\Phi} = \Phi(X,Y,Z,t) = \mathbf{Re} \left\{ \phi(X,Y,Z)e^{i\omega t} \right\} \,[/math]

where [math]\omega\,[/math] is the encounter frequency and

[math] \phi = \phi_I + \sum_{j=1}^7 \phi_j \,[/math]

with [math] \phi_j, \ j=1,\cdots,6\,[/math] being the radiation and [math] \phi_7\,[/math] being the diffraction potentials.

Boundary-value problem for [math] \bar{\Phi}[/math]

Free surface condition:

[math] U^2 \bar{\Phi}_{XX} + g \bar{\Phi}_Z = 0, \quad Z=0 [/math]

Ship-hull condition

[math] \vec{n} \cdot \nabla\vec{\Phi} = \vec{n} \cdot \vec{U} = U n_1 \, [/math]

where [math] \vec{n} = ( n_1, n_2, n_3 ) \, [/math] is the unit vector pointing inside the ship hull.

Far from the ship [math] \bar{\Phi}\,[/math] represents outgoing waves which are known as the Kelvin ship wake studied earlier.

The solution for [math] \bar{\Phi}\,[/math] by the above formulation known as the Neumann-Kelvin problem and its generalizations discussed in the literature is carried out by panel methods.

The principal output quantities of interest in practice are:

  • Free-surface elevation
[math] \zeta = - \frac{1}{g} \left. \frac{\mathrm{d}\bar{\Phi}}{\mathrm{d}t} \right|_{\mbox{Inertial frame}} = \frac{U}{g} \left. \bar{\Phi}_X \right|_{\mbox{Ship frame}}, \quad Z=0 [/math]
  • Hydrodynamic pressure (linear)
[math] P = - \rho \left. \frac{\mathrm{d}\bar{\mathbf{\Phi}}}{\mathrm{d}t} \right|_{\mbox{Inertial frame}} = \rho U \left. \bar{\Phi}_X \right|_{\mbox{Ship frame}} [/math]
  • Hydrodynamic pressure (total)
[math] P_T = \rho \left( \frac{\mathrm{d}\bar{\mathbf{\Phi}}}{\mathrm{d}t} = \frac{1}{2} \nabla \bar{\mathbf{\Phi}} \cdot \nabla \bar{\mathbf{\Phi}} + g Z \right)_{\mbox{Inertial frame}} [/math]
[math] = - \rho \left( - U \bar{\Phi_X} + \frac{1}{2} \nabla \bar{\Phi} \cdot \nabla \bar{\Phi} + g Z \right)_{\mbox{Ship frame}} \,[/math]

If [math]\bar{S_w}\,[/math] is the ship wetted surface due to its steady forward translation on a free surface and [math]\bar{n}\,[/math] is the unit normal vector pointing out of the fluid domain the ship ideal-fluid force is given by

[math] \overrightarrow{F} = \iint_{\overline{S_W}} P_T \vec{n} \mathrm{d}S \, [/math]

The wave resistance is: [math] R_W = \vec{i} \cdot \vec{F} \, [/math].

We will derive boundary value problems for the potentials [math] \bar{\Phi}\,[/math] and [math] \Phi\,[/math] relative to the ship fixed frame.

The principal assumption underlying the ensuing derivation is that the ship is slender, thin or flat or in general streamlined in the longitudinal direction. More explicitly, if [math] B\,[/math] is the ship beam, [math] T\,[/math] its draft and [math] L\,[/math] its length we will assume that:

[math] \frac{B}{L}, \quad \frac{T}{L} = O (\varepsilon), \quad \varepsilon \ll 1 \,[/math]

where [math] \varepsilon\,[/math] is the slenderness parameter assumed to be small compared to [math] 1 \,[/math].

The ship slenderness justifies the use of the linear free-surface condition in the forward-speed problem for a broad range of speeds and hull shapes.

Boundary-value problem for time-Harmonic velocity potential

From the Galilean transformation:

[math] \left( \frac{\partial}{\partial t} - U \frac{\partial}{\partial X} \right)^2 \Phi + g \Phi_Z = 0, \quad Z=0 \, [/math]

Relative to the ship frame. In terms of the complex potential:

[math] \Phi = \mathbf{Re} \left\{ \phi e^{i\omega t} \right\} \, [/math]
[math] \left( i\omega - U \frac{\partial}{\partial X} \right)^2 \phi + g \phi_Z = 0 , \quad Z=0 \, [/math]

where [math] \omega\, [/math] is the encounter frequency and [math] \phi\,[/math] is any of the [math] \phi_j\,[/math] potentials.

  • The above time harmonic Neumann-Kelvin free surface condition is being treated by state-of-the-art panel methods. An important simplification for slender ships and large values of [math]\omega\,[/math] will lead to the popular strip theory.
  • The solution for [math] \bar{\Phi}\,[/math] is far from simple numerically. A lot of research has been devoted to this effort, in particular towards the evaluation of the ship Kelvin wake and the ship wave resistance.
  • The linearization of the pressure and vessel wetted surface [math]\overline{S_W}\,[/math] about its static value in calm water must be carried out carefully! Nonlinear effects are known to contribute appreciably to the wave resistance.
  • If available, a fully nonlinear solution of the forward-speed steady ship wave problem is preferable. Numerical issues must be carefully treated and are the subject of state-of-the-art research.
  • Coupling with viscous effects is often strong and important for predicting the total resistance of the ship.

Relative to the ship-fixed coordinate system the ambient wave elevation oscillates with frequency [math]\omega\,[/math].


[math] \zeta = - \frac{1}{g} \left( \frac{\partial}{\partial t} - U \frac{\partial}{\partial X} \right) \Phi_I [/math]


[math] \Phi_I = \mathbf{Re} \left\{ \frac{i g A}{\omega_0} e^{KZ-iKX\cos\beta-iKY\sin\beta+i\omega t} \right\} [/math]
[math] \frac{\partial}{\partial t} = i \omega, \qquad \frac{\partial}{\partial X} = - i K \cos \beta [/math]
[math] \frac{\partial}{\partial t} - U \frac{\partial}{\partial X} = i \left( \omega + K U \cos \beta \right) = i \omega_0 \, [/math]


[math] \zeta = \mathbf{Re} \left\{ A e^{-iKX\cos\beta-iKY\sin\beta+i\omega t} \right\} \,[/math]
[math] = \mathbf{Re} \left\{ A e^{i\omega t} \right\}_{X=Y=0} [/math]

where [math]\Pi_j(\omega) \,[/math] is the complex amplitude of the vessel motion in mode [math]-j\,[/math], a function of the frequency of encounter [math] \omega\,[/math], known as the response amplitude operator (RAO).

The ship equations of motion follow as in the [math] U=0\,[/math] case using linear system theory:

[math] \sum_{j=1}^6 \left[ - \omega^2 \left(M_{ij} + A_{ij} \right) + i\omega B_{ij} + C_{ij} \right] \Pi_j (\omega) = \mathbf{X}_i(\omega),\quad i=1,\cdots,6 [/math]

where the hydrodynamic coefficients [math]A_{ij}(\omega), \ B_{ij}(\omega) \, [/math] and exciting forces are now functions of the encounter frequency [math]\omega\,[/math] and other forward-speed effects.

  • Inertia & hydrostatic matrices the same as in the zero-speed case
  • Will derive BVP's governing the coefficients [math]A_{ij}(\omega), \ B_{ij}(\omega) \, [/math] and exciting forces [math]\mathbf{X}_i (\omega) \, [/math].


[math] \omega = \omega_0 - U \frac{\omega_0^2}{g} \cos\beta \, [/math]
  • [math] \omega \gt \ \mbox{or} \ \lt 0 \, [/math]: Both positive and negative values of [math] \omega\,[/math] are possible. In practice will always deal with the absolute value of [math] \omega\,[/math].
  • Given the absolute wave frequency [math] \omega_0 \gt 0 \, [/math] there exists a unique [math] \omega\,[/math].
  • Given a positive absolute encounter frequency [math] |\omega|\,[/math], there exist possibly multiple [math] \omega_0\,[/math]'s satisfying the above relation. More discussion of this topic will follow.
  • Assuming small amplitude motions the ship responses are modeled after linear system theory, input signal [math]\sim e^{i\omega t} \ \longrightarrow\,[/math] output signal [math] \sim e^{i\omega t} \,[/math].

Relative to the earth-fixed frame the ambient wave velocity potential takes the form:

[math] \Phi_I = \mathbf{Re} \left\{ \phi_I \right\} \,[/math]
[math] \phi_I = \frac{i g A}{\omega_0} e^{KZ_0 - iKX_0 \cos\beta - iKY_0 \sin\beta +i\omega_0 t} \, [/math]

where in deep water: [math] K = \frac{\omega_0^2}{g} \, [/math]


[math] X_0 = x+ Ut \qquad Y_0 = y \qquad Z_0 = z \,[/math]

It follows that:

[math] \phi_I = \frac{i g A}{\omega_0} e^{ Kz - iKx\cos\beta - iKy\sin\beta + i (\omega_0 -UK \cos\beta) t} [/math]


[math] \omega = \omega_0 - UK\cos\beta \,[/math]

Be defined to be the encounter frequency which accounts for the Doppler effect included in the second term in the RHS.

Comments on N-K formulation

  • The ship is assumed to be streamlined in order to justify the decomposition of the steady & time harmonic components.
  • The vessel motions are assumed small and of the same order as the ambient wave amplitude. Terms omitted are of [math]O\left(A^2\right)\,[/math].
  • When Taylor expanding the free-surface and body-boundary condition about [math]Z=0\,[/math] and [math]\overrightarrow{S_B}\,[/math] respectively, the steady flow potential [math] \bar{\Phi}\simeq 0 \,[/math].
  • For ships with appreciable thickness a better approximation for [math]\bar\Phi\,[/math] is that of the double-body flow disturbance such that [math] \overrightarrow{\Phi_Z} = 0\,[/math] on [math] Z=0\,[/math] and [math] \overrightarrow{\Phi_n} = U n\,[/math], on [math]\overrightarrow{S_B}\,[/math], This leads to the state-of-the-art linear 3D steady flow and seakeeping formulation discussed later in connection with panel methods.
  • The N-K formulation is the staring point of strip theory.

Ocean Wave Interaction with Ships and Offshore Energy Systems