Difference between revisions of "Ship Kelvin Wake"

From WikiWaves
Jump to navigationJump to search
 
(4 intermediate revisions by one other user not shown)
Line 28: Line 28:
 
The same potential expressed relative to the ship frame is <math> \bar{\Phi}(\mathbf{\bar{x}},t) \, </math>.
 
The same potential expressed relative to the ship frame is <math> \bar{\Phi}(\mathbf{\bar{x}},t) \, </math>.
 
The relation between the two potentials is given by the identity
 
The relation between the two potentials is given by the identity
<center><math> \Phi(x,y,z,t) = \bar{\phi}(\bar{x}-Ut,Y,Z,t) \,</math></center>
+
<center><math> \Phi(x,y,z,t) = \bar{\Phi}(\bar{x},y,z,t)
 +
= \bar{\Phi}(x-Ut,y,z,t) \,</math></center>
 
where the relation between the coordinates of the two coordinate systems has been introduced.
 
where the relation between the coordinates of the two coordinate systems has been introduced.
 +
Note the time dependence occurs in two places in <math> \bar{\Phi} \,</math> and in one place in <math> \Phi \,</math>. The governing equations are always derived relative to the earth coordinate system and time derivatives are initially taken on <math> \Phi  \, </math>.
 +
Therefore
 +
<center><math> \frac{\mathrm{d}\Phi}{\mathrm{d}t} = \frac{\mathrm{d}}{\mathrm{d}t} \bar{\phi} ( x-Ut,y,z,t)
 +
= \frac{\partial\phi}{\partial t} - U \frac{\partial\phi}{\partial x} </math></center>
  
Note the time dependence in the arguments of the two potentials. It occurs in two places in <math> \phi \,</math> and in one place in <math> \Phi \,</math>. The governing equations are always derived relative to the earth coordinate system and time derivatives are initially taken on <math> \Phi  \, </math>:
+
All time derivatives of the earth fixed velocity potential <math> \Phi \, </math> which appear in the free surface condition and the Bernoulli equation can be expressed in terms of derivatives of <math> \bar{\Phi} \, </math> using the Galilean transformation derived above.
 
 
<center><math> \frac{\mathrm{d}\Phi}{\mathrm{d}t} = \frac{\mathrm{d}}{\mathrm{d}t} \phi ( X-Ut,y,z,t) = \frac{\partial\phi}{\partial x} \frac{\mathrm{d}x}{\mathrm{d}t} + \frac{\partial\phi}{\partial t} = \frac{\partial\phi}{\partial t} - U \frac{\partial\phi}{\partial x} </math></center>
 
 
 
The main result to follow is:
 
 
 
<center><math> \frac{\mathrm{d}\Phi}{\mathrm{d}t} = \frac{\partial\phi}{\partial t} - U \frac{\partial\phi}{\partial x}. </math></center>
 
 
 
All time derivatives of the earth fixed velocity potential <math> \Phi \, </math> which appear in the free surface condition and the Bernoulli equation can be expressed in terms of derivatives of <math> \phi \, </math> using the Galilean transformation derived above.
 
  
 
If the flow is steady relative to the ship fixed coordinate system
 
If the flow is steady relative to the ship fixed coordinate system
<center><math> \frac{\partial\phi}{\partial t} = 0 \, </math></center>
+
<center><math> \frac{\partial\bar{\Phi}}{\partial t} = 0 \, </math></center>
 
but
 
but
<center><math> \frac{\mathrm{d}\Phi}{\mathrm{d}t} = -U \frac{\partial\phi}{\partial x} \,</math></center>
+
<center><math> \frac{\mathrm{d}\Phi}{\mathrm{d}t} = -U \frac{\partial\bar{\Phi}}{\partial x} \,</math></center>
 
or, the ship wake is stationary relative to the ship but not relative to an observed on the beach.
 
or, the ship wake is stationary relative to the ship but not relative to an observed on the beach.
 
=== Linear free surface conditions ===
 
 
The linear free surface conditions are as follows,
 
 
<center><math> \frac{\mathrm{d}^2\Phi}{\mathrm{d}t^2} + g \frac{\mathrm{d}\Phi}{\mathrm{d}Z} = 0, \quad Z=0 </math></center>
 
 
We can calculate the second derivative above as
 
<center><math> \frac{\mathrm{d}\Phi}{\mathrm{d}t} = \frac{\partial\phi}{\partial t} - U \frac{\partial\phi}{\partial x} </math></center>
 
and
 
<center><math> \frac{\mathrm{d}^2\Phi}{\mathrm{d}t^2}  = \frac{\partial^2 \phi}{\partial t^2} - 2 U \frac{\partial^2\phi}{\partial x \partial t} + U^2 \frac{\partial^ \phi}{\partial x^2} </math></center>
 
 
This is the free surface condition governing the forward-speed linear seakeeping problem.
 
 
When no ambient waves are present <math> \frac{\partial\phi}{\partial t} = 0 \, </math> and we obtain the free surface condition for the steady Kelvin ship wave problem
 
 
<center><math> U^2 \frac{\partial^2 \phi}{\partial x^2} + g \frac{\partial\phi}{\partial z} = 0, \qquad z=0 </math></center>
 
 
This is the famous Neumann-Kelvin free-surface condition governing the linear steady wave pattern generated by a translating ship.
 
 
Bernoulli equation:
 
 
<center><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\phi}{\partial t} - U \frac{\partial\phi}{\partial x} \right) + \frac{1}{2} \nabla \phi \cdot \nabla \phi + g Z \right] </math></center>
 
 
Free surface elevation:
 
 
<center><math> \zeta = - \left. \frac{1}{g} \frac{\mathrm{d}\Phi}{\mathrm{d}t} \right|_{z=0} = -\frac{1}{g} \left. \left( \frac{\partial\phi}{\partial t} - U \frac{\partial\phi}{\partial x} \right)\right|_{z=0} </math></center>
 
 
In the Kelvin ship wave problem
 
 
<center><math> \frac{\partial\phi}{\partial t} = 0 </math></center>
 
 
Hence
 
 
<center><math> \zeta = \frac{U}{g} \frac{\partial\phi}{\partial x}, \qquad z=0 </math></center>
 
 
So if the velocity potential <math> \phi(\overrightarrow{x}) \,</math> is available in some form relative to the translating frame, the wave pattern follows from the above expression.
 
  
 
== Kelvin wake ==
 
== Kelvin wake ==
Line 90: Line 50:
 
Local view of Kelvin wake consists approximately of a [[Linear Plane Progressive Regular Waves| plane progressive wave group]] propagating in direction <math> \theta\,\!</math>.
 
Local view of Kelvin wake consists approximately of a [[Linear Plane Progressive Regular Waves| plane progressive wave group]] propagating in direction <math> \theta\,\!</math>.
 
As noted above surface wave systems of general form always consist of combinations of plane progressive waves of different frequencies and directions. The same model will apply to the ship kelvin wake.
 
As noted above surface wave systems of general form always consist of combinations of plane progressive waves of different frequencies and directions. The same model will apply to the ship kelvin wake.
Relative to the earth frame, the local plane progressive wave takes the form:
+
Relative to the earth frame, the local plane wave in [[Infinite Depth]] takes the form
 
 
<center><math> \Phi = \frac{igA}{\omega} e^{KZ-iK(X\cos\theta+Y\sin\theta)+i\omega t} </math></center>
 
 
 
Substitute
 
 
 
<center><math> X = x + Ut, Y=y, Z=z. \,</math></center>
 
  
and use the results from [[Translating Coordinate Systems]] to give us, relative to the ship frame
+
<center><math> \Phi = \frac{\mathrm{i}gA}{\omega} e^{kz-\mathrm{i}k(x\cos\theta+y\sin\theta)+ \mathrm{i}\omega t} </math></center>
  
<center><math> \phi = \frac{igA}{\omega} e^{KZ-iK(x\cos\theta+y\sin\theta)-i(kU\cos\theta-\omega)t} </math></center>
+
Relative to the ship frame
 +
<center><math>bar{\Phi} = \frac{igA}{\omega} e^{kz-\mathrm{i}k(x\cos\theta+y\sin\theta)-\mathrm{i}(kU\cos\theta-\omega)t} </math></center>
  
 
But relative to the ship frame waves are stationary, so we must have:
 
But relative to the ship frame waves are stationary, so we must have:
  
<center><math> KU\cos\theta = \omega \,</math></center>
+
<center><math> kU\cos\theta = \omega \,</math></center>
 
or
 
or
<center><math> \frac{\omega}{K} = V_P = U \cos \theta \, </math></center>
+
<center><math> \frac{\omega}{k} = C_p = U \cos \theta \, </math></center>
 
This implies the following
 
This implies the following
 
* The phase velocity of the waves in the kelvin wake propagating in direction <math>\theta\,</math> must be equal to <math> U\cos\theta\,</math>, otherwise they cannot be stationary relative to the ship.
 
* The phase velocity of the waves in the kelvin wake propagating in direction <math>\theta\,</math> must be equal to <math> U\cos\theta\,</math>, otherwise they cannot be stationary relative to the ship.
* Relative to the earth system the frequency of a local system propagating in direction <math> \theta \, </math> is given by the relation
+
* Relative to the earth system the frequency of a local system propagating in direction <math> \theta \, </math> is given by the relation <math> \omega = kU \cos \theta \, </math>
<center><math> \omega = K U \cos \theta \, </math></center>
 
where <math> K = \frac{2\pi}{\lambda}\, </math>.
 
 
 
Relative to the earth system the deep water dispersion relation states:
 
 
 
<center><math> \omega^2 = g K \, </math></center>
 
 
 
or
 
 
 
<center><math> K^2 U^2 \cos^2 \theta = g K  \,</math></center>
 
 
 
<center><math> \Longrightarrow  K = \frac{g}{U^2 \cos^2 \theta} = \frac{2\pi}{\lambda(\theta)} </math></center>
 
  
<center><math> \Longrightarrow \lambda(\theta) = \frac{2\pi U^2 \cos^2 \theta}{g} \, </math></center>
+
Relative to the earth system the [[Infinite Depth]] [[Dispersion Relation for a Free Surface]] states
 +
<center><math> \omega^2 = g k \, </math></center>
 +
so that
 +
<center><math> \lambda(\theta) = \frac{2\pi U^2 \cos^2 \theta}{g} \, </math></center>
  
This is the wavelength of waves in a kelvin wake propagating in direction <math> \theta \, </math> and being stationary relative to the ship.
+
This is the wavelength of waves in a Kelvin wake propagating in direction <math> \theta \, </math> which
 +
are stationary relative to the ship.
  
 
== Application of the Group velocity ==
 
== Application of the Group velocity ==

Latest revision as of 11:37, 6 November 2010

Wave and Wave Body Interactions
Current Chapter Ship Kelvin Wake
Next Chapter Linear Wave-Body Interaction
Previous Chapter Wavemaker Theory



Introduction

Wake created behind a ship

A ship moving over the surface of undisturbed water sets up waves emanating from the bow and stern of the ship. The waves created by the ship consist of divergent and transverse waves. The divergent wave are observed as the wake of a ship with a series of diagonal or oblique crests moving outwardly from the point of disturbance. These wave were first studied by Lord Kelvin


Translating Coordinate System

We have a the standard fixed coordinate system [math]\displaystyle{ x, y, z }[/math] and a moving coordinate systems which is moving in the [math]\displaystyle{ x }[/math] direction with speed [math]\displaystyle{ U }[/math]. We denote the moving coordinate systems in the [math]\displaystyle{ x }[/math] direction by

[math]\displaystyle{ \bar{x} = x + U t }[/math]

Let [math]\displaystyle{ \Phi(\mathbf{x},t) \, }[/math] be the velocity potential describing the potential flow generated by the ship relative to the earth frame. The same potential expressed relative to the ship frame is [math]\displaystyle{ \bar{\Phi}(\mathbf{\bar{x}},t) \, }[/math]. The relation between the two potentials is given by the identity

[math]\displaystyle{ \Phi(x,y,z,t) = \bar{\Phi}(\bar{x},y,z,t) = \bar{\Phi}(x-Ut,y,z,t) \, }[/math]

where the relation between the coordinates of the two coordinate systems has been introduced. Note the time dependence occurs in two places in [math]\displaystyle{ \bar{\Phi} \, }[/math] and in one place in [math]\displaystyle{ \Phi \, }[/math]. The governing equations are always derived relative to the earth coordinate system and time derivatives are initially taken on [math]\displaystyle{ \Phi \, }[/math]. Therefore

[math]\displaystyle{ \frac{\mathrm{d}\Phi}{\mathrm{d}t} = \frac{\mathrm{d}}{\mathrm{d}t} \bar{\phi} ( x-Ut,y,z,t) = \frac{\partial\phi}{\partial t} - U \frac{\partial\phi}{\partial x} }[/math]

All time derivatives of the earth fixed velocity potential [math]\displaystyle{ \Phi \, }[/math] which appear in the free surface condition and the Bernoulli equation can be expressed in terms of derivatives of [math]\displaystyle{ \bar{\Phi} \, }[/math] using the Galilean transformation derived above.

If the flow is steady relative to the ship fixed coordinate system

[math]\displaystyle{ \frac{\partial\bar{\Phi}}{\partial t} = 0 \, }[/math]

but

[math]\displaystyle{ \frac{\mathrm{d}\Phi}{\mathrm{d}t} = -U \frac{\partial\bar{\Phi}}{\partial x} \, }[/math]

or, the ship wake is stationary relative to the ship but not relative to an observed on the beach.

Kelvin wake

Diagram of the Kelvin Wake

Local view of Kelvin wake consists approximately of a plane progressive wave group propagating in direction [math]\displaystyle{ \theta\,\! }[/math]. As noted above surface wave systems of general form always consist of combinations of plane progressive waves of different frequencies and directions. The same model will apply to the ship kelvin wake. Relative to the earth frame, the local plane wave in Infinite Depth takes the form

[math]\displaystyle{ \Phi = \frac{\mathrm{i}gA}{\omega} e^{kz-\mathrm{i}k(x\cos\theta+y\sin\theta)+ \mathrm{i}\omega t} }[/math]

Relative to the ship frame

[math]\displaystyle{ bar{\Phi} = \frac{igA}{\omega} e^{kz-\mathrm{i}k(x\cos\theta+y\sin\theta)-\mathrm{i}(kU\cos\theta-\omega)t} }[/math]

But relative to the ship frame waves are stationary, so we must have:

[math]\displaystyle{ kU\cos\theta = \omega \, }[/math]

or

[math]\displaystyle{ \frac{\omega}{k} = C_p = U \cos \theta \, }[/math]

This implies the following

  • The phase velocity of the waves in the kelvin wake propagating in direction [math]\displaystyle{ \theta\, }[/math] must be equal to [math]\displaystyle{ U\cos\theta\, }[/math], otherwise they cannot be stationary relative to the ship.
  • Relative to the earth system the frequency of a local system propagating in direction [math]\displaystyle{ \theta \, }[/math] is given by the relation [math]\displaystyle{ \omega = kU \cos \theta \, }[/math]

Relative to the earth system the Infinite Depth Dispersion Relation for a Free Surface states

[math]\displaystyle{ \omega^2 = g k \, }[/math]

so that

[math]\displaystyle{ \lambda(\theta) = \frac{2\pi U^2 \cos^2 \theta}{g} \, }[/math]

This is the wavelength of waves in a Kelvin wake propagating in direction [math]\displaystyle{ \theta \, }[/math] which are stationary relative to the ship.

Application of the Group velocity

An observer sitting on an earth fixed frame observes a local wave system propagating in direction [math]\displaystyle{ \theta\, }[/math] travelling at its group velocity [math]\displaystyle{ \frac{d\omega}{dK}\, }[/math] by virtue of the Rayleigh device which states that we need to focus on the speed of the energy density ([math]\displaystyle{ \sim \, }[/math] wave amplitude) rather than the speed of wave crests. So, relative to the earth fixed inclined coordinate system [math]\displaystyle{ (X', Y') \, }[/math]:

[math]\displaystyle{ \frac{X'}{t} = V_g = \frac{d\omega}{dK} \, }[/math]

Or

[math]\displaystyle{ X' = \frac{d\omega}{dK} t \ \Longrightarrow \ \frac{d}{dK} (K X' - \omega t) = 0 }[/math]
[math]\displaystyle{ X' = X \cos \theta + Y \sin \theta = x \cos\theta + y\sin\theta + Ut \cos\theta \, }[/math]

So:

[math]\displaystyle{ KX' - \omega t = K ( x\cos\theta + y\sin\theta ) + (KU\cos\theta - \omega) t \, }[/math]

However [math]\displaystyle{ (KU\cos\theta - \omega) t =0 }[/math] so that the Rayleigh condition for the velocity of the group takes the form:

[math]\displaystyle{ \frac{d}{dK} [ K(\theta) (x\cos\theta + y\sin\theta) ] = 0 \, }[/math]

By virtue of the dispersion relation derived above:

[math]\displaystyle{ K(\theta) = \frac{g}{U^2 \cos^2 \theta} \, }[/math]

It follows from the chain rule of differentiation that Rayleigh's condition is:

[math]\displaystyle{ \frac{d}{d\theta} \left[ \frac{g}{U^2\cos^2\theta} (x\cos\theta + y\sin\theta) ]\right] = 0 \, }[/math]

At the position of the Kelvin waves which are locally observed by an observer at the beach.

  • So the "visible" waves in the wake of a ship are wave groups which must travel at the local group velocity. These conditions translate into the above equation which will be solved and discussed next. More discussion and a more mathematical derivation based on the principle of stationary phase can be found in Newman 1977.

Solution of the equation for angle

Graphical image of the equations

The solution of the above equation will produce a relation between [math]\displaystyle{ \frac{y}{x} \, }[/math] and [math]\displaystyle{ \theta\, }[/math]. So local waves in a Kelvin wake can only propagate in a certain direction [math]\displaystyle{ \theta\, }[/math], given [math]\displaystyle{ \frac{y}{x} \, }[/math].

Simple algebra leads to:

[math]\displaystyle{ \frac{y}{x} = - \frac{\cos\theta\sin\theta}{1+sin^2\theta} = \frac{y}{x}(\theta) }[/math]

which implies that

  • [math]\displaystyle{ \frac{y}{x}(\theta) \, }[/math] is anti-symmetric about [math]\displaystyle{ \theta = 0 \, }[/math] and each part corresponds to the Kelvin wake in the port and starboard sides of the vessel. The physics on

either side is identical due to symmetry.

  • [math]\displaystyle{ \theta = 0 \, }[/math]: waves propagating in the same direction as the ship. These waves can only exist at [math]\displaystyle{ Y = 0\, }[/math] as seen above.
  • [math]\displaystyle{ \theta = \frac{\pi}{2} \, }[/math]: waves propagating at a [math]\displaystyle{ 90^\circ\, }[/math] angle relative to the ship direction of forward translation.
  • [math]\displaystyle{ \theta = 35^\circ 16' \, }[/math]: (or 35,26°) waves propagating at an angle [math]\displaystyle{ \theta = 35^\circ 16' \, }[/math] relative to the ship axis. These are waves seen at the caustic of the Kelvin wake.

Let the solution of [math]\displaystyle{ \frac{y}{x}(\theta) \, }[/math] be of the form, when inverted:

[math]\displaystyle{ \mbox{Region I}: \qquad \theta = f_1 (\frac{y}{x}) \, }[/math]
[math]\displaystyle{ \mbox{Region II}: \qquad \theta = f_2 (\frac{y}{x}) \, }[/math]

Note that observable waves cannot exist for values of [math]\displaystyle{ \frac{y}{x}\, }[/math] that exceed the value shown in the figure or [math]\displaystyle{ \left. \frac{y}{x} \right|_{Max} = 2^{-3/2} \, }[/math]. This translates into a value for the corresponding angle equal to [math]\displaystyle{ 19^\circ 28' \, }[/math] (or 19,47°) which is the angle of the caustic for any speed [math]\displaystyle{ U\, }[/math].

"transverse" and "divergent" wave systems in the Kelvin wake

The crests of the wave system trailing a ship, the Kelvin wake, are curves of constant phase of:

[math]\displaystyle{ \frac{x\cos\theta+y\sin\theta}{\cos^2\theta} = \mathbf{C} \, }[/math]

In [math]\displaystyle{ \mbox{Region I} \, }[/math]:

[math]\displaystyle{ \mathbf{C} = \frac{x\cos f_1(\frac{y}{x}) + y\sin f_1 (\frac{y}{x})}{\cos^2 f_1 (\frac{y}{x})} \equiv G_1 (\frac{y}{x}) }[/math]

In [math]\displaystyle{ \mbox{Region II} \, }[/math]:

[math]\displaystyle{ \mathbf{C} = \frac{x\cos f_2(\frac{y}{x}) + y\sin f_2 (\frac{y}{x})}{\cos^2 f_2 (\frac{y}{x})} \equiv G_2 (\frac{y}{x}) }[/math]

Plotting these curves we obtain a visual graph of the "transverse" and "divergent" wave systems in the Kelvin wake.


This article is based on the MIT open course notes and the original article can be found here

Ocean Wave Interaction with Ships and Offshore Energy Systems