Difference between revisions of "Ship Kelvin Wake"

From WikiWaves
Jump to navigationJump to search
 
(40 intermediate revisions by 5 users not shown)
Line 1: Line 1:
<u>Kelvin wake</u>
+
{{Ocean Wave Interaction with Ships and Offshore Structures
 +
| chapter title = Ship Kelvin Wake
 +
| next chapter = [[Linear Wave-Body Interaction]]
 +
| previous chapter =  [[Wavemaker Theory]]
 +
}}
  
Local view of kelvin wake consists approximately of a plane progressive wave group propagating in direction <math> \theta\,</math>.
+
{{incomplete pages}}
  
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.
+
= Introduction =
  
Relative to the earth frame, the local plane progressive wave takes the form:
+
[[Image:Wake.avon.gorge.arp.750pix.jpg|thumb|right|450px|Wake created behind a ship]]
  
<center><math> \Phi = \frac{igA}{\omega} e^{KZ-iK(X\cos\theta+Y\sin\theta)+i\omega t} </math></center>
+
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 [http://en.wikipedia.org/wiki/Wake 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 [http://en.wikipedia.org/wiki/Lord_Kelvin Lord Kelvin]
  
Substitute:
 
  
<center><math> X = x + Ut, Y=y, Z=z. \,</math></center>
 
  
Relative to the ship frame:
+
== Translating Coordinate System ==
  
<center><math> \phi = \frac{igA}{\omega} e^{KZ-iK(x\cos\theta+y\sin\theta)-i(kU\cos\theta-\omega)t} </math></center>
+
We have a the standard fixed coordinate system <math>x, y, z</math> and a moving coordinate systems
 +
which is moving in the <math>x</math> direction with speed <math>U</math>. We denote the moving coordinate
 +
systems in the <math>x</math> direction by
 +
<center><math>
 +
\bar{x} = x + U t
 +
</math></center>
  
But relative to the ship frame waves are stationary, so we must have:
+
Let <math> \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> \bar{\Phi}(\mathbf{\bar{x}},t) \, </math>.
 +
The relation between the two potentials is given by the identity
 +
<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.
 +
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>
  
<center><math> KU\cos\theta = \omega \,</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> \bar{\Phi} \, </math> using the Galilean transformation derived above.
  
<center><math> \frac{\omega}{K} = V_P = U \cos \theta \, </math></center>
+
If the flow is steady relative to the ship fixed coordinate system
 +
<center><math> \frac{\partial\bar{\Phi}}{\partial t} = 0 \, </math></center>
 +
but
 +
<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.
  
* 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.
+
== Kelvin wake ==
  
* Relative to the earth system the frequency of a local system propagating in direction <math> \theta \, </math> is given by the relation
+
[[Image:Kelvin_wake.jpg|thumb|right|600px|Diagram of the Kelvin Wake]]
  
<center><math> \omega = K U \cos \theta \, </math></center>
+
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.
 +
Relative to the earth frame, the local plane wave in [[Infinite Depth]] takes the form
  
where <math> K = \frac{2\pi}{\lambda}\, </math>.
+
<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>
  
Relative to the earth system the deep water dispersion relation states:
+
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>
  
<center><math> \omega^2 = g K \, </math></center>
+
But relative to the ship frame waves are stationary, so we must have:
  
 +
<center><math> kU\cos\theta = \omega \,</math></center>
 
or
 
or
 +
<center><math> \frac{\omega}{k} = C_p = U \cos \theta \, </math></center>
 +
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.
 +
* 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> K^2 U^2 \cos^2 \theta = g K \quad \Longrightarrow \, </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>
<center><math> K = \frac{g}{U^2 \cos^2 \theta} = \frac{2\pi}{\lambda(\theta)} </math></center>
+
so that
 
+
<center><math> \lambda(\theta) = \frac{2\pi U^2 \cos^2 \theta}{g} \, </math></center>
<center><math> \Longrightarrow \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.
  
An observer sitting on an earth fixed frame observes a local wave system propagating in direcion <math> \theta\,</math> travelling at its <u>group</u> velocity <math>\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> \sim \, </math> wave amplitude) rather than the speed of wave crests!
+
== Application of the Group velocity ==
  
So, relative to the earth fixed inclined coordinate system <math> (x', y') \,</math>:
+
An observer sitting on an earth fixed frame observes a local wave system propagating in direction <math> \theta\,</math> travelling at its [[Wave Energy Density and Flux|group velocity]] <math>\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> \sim \, </math> wave amplitude) rather than the speed of wave crests. So, relative to the earth fixed inclined coordinate system <math> (X', Y') \,</math>:
  
<center><math> \frac{x'}{t} = V_g = \frac{d\omega}{dK} \, </math></center>
+
<center><math> \frac{X'}{t} = V_g = \frac{d\omega}{dK} \, </math></center>
  
 
Or
 
Or
  
<center><math> X' = \frac{d\omega}{dK} t \ \Longrightarrow \ \frac{d}{dK} (K x' - \omega t) = 0 </math></center>
+
<center><math> X' = \frac{d\omega}{dK} t \ \Longrightarrow \ \frac{d}{dK} (K X' - \omega t) = 0 </math></center>
  
 
<center><math> X' = X \cos \theta + Y \sin \theta = x \cos\theta + y\sin\theta + Ut \cos\theta \,</math></center>
 
<center><math> X' = X \cos \theta + Y \sin \theta = x \cos\theta + y\sin\theta + Ut \cos\theta \,</math></center>
Line 59: Line 88:
 
So:
 
So:
  
<center><math> Kx' - \omega t = K ( x\cos\theta + y\sin\theta ) + (KU\cos\theta - \omega) t \,</math></center>
+
<center><math> KX' - \omega t = K ( x\cos\theta + y\sin\theta ) + (KU\cos\theta - \omega) t \,</math></center>
 
+
However <math>(KU\cos\theta - \omega) t =0</math> so that the Rayleigh condition for the velocity of the group takes the form:
And the Rayleigh condition for the velocity of the group takes the form:
 
  
 
<center><math> \frac{d}{dK} [ K(\theta) (x\cos\theta + y\sin\theta) ] = 0 \, </math></center>
 
<center><math> \frac{d}{dK} [ K(\theta) (x\cos\theta + y\sin\theta) ] = 0 \, </math></center>
Line 67: Line 95:
 
By virtue of the dispersion relation derived above:
 
By virtue of the dispersion relation derived above:
  
<center><math> K(\theta) = \frac{g}{U^2 \cos^\theta} \, </math></center>
+
<center><math> K(\theta) = \frac{g}{U^2 \cos^2 \theta} \, </math></center>
  
 
It follows from the chain rule of differentiation that Rayleigh's condition is:
 
It follows from the chain rule of differentiation that Rayleigh's condition is:
Line 75: Line 103:
 
At the position of the Kelvin waves which are locally observed by an observer at the beach.
 
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 <u>wave groups</u> 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 derivaion based on the principle of stationary phase may be found in MH.
+
* So the "visible" waves in the wake of a ship are <u>wave groups</u> 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 [http://en.wikipedia.org/wiki/Principle_of_stationary_phase principle of stationary phase] can be found in [[Newman 1977]].
  
The solution of the above equation will produce a relation between <math> \frac{y}{x} \,</math> and <math> \theta\,</math>. So local waves in a Kelvin wake can only propagate in a certain direcion <math> \theta\,</math>, given <math> \frac{y}{x} \, </math>.
+
== Solution of the equation for angle ==
 +
 
 +
[[Image:Graphical_image.jpg|thumb|right|600px|Graphical image of the equations]]
 +
 
 +
The solution of the above equation will produce a relation between <math> \frac{y}{x} \,</math> and <math> \theta\,</math>. So local waves in a Kelvin wake can only propagate in a certain direction <math> \theta\,</math>, given <math> \frac{y}{x} \, </math>.
  
 
Simple algebra leads to:
 
Simple algebra leads to:
  
 
<center><math> \frac{y}{x} = - \frac{\cos\theta\sin\theta}{1+sin^2\theta} = \frac{y}{x}(\theta)  
 
<center><math> \frac{y}{x} = - \frac{\cos\theta\sin\theta}{1+sin^2\theta} = \frac{y}{x}(\theta)  
 
 
</math></center>
 
</math></center>
 +
which implies that
 +
* <math>\frac{y}{x}(\theta) \,</math> is anti-symmetric about <math> \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>\frac{y}{x}(\theta) \,</math> is anti-symmetric about <math> \theta = 0 \, </math> each
+
* <math> \theta = 0 \,</math>: waves propagating in the same direction as the ship. These waves can only exist at <math> Y = 0\, </math> as seen above.
  
part corresponds to the Kelvin wake in the port and starboard sides of the vessel. The physics on
+
* <math> \theta = \frac{\pi}{2} \, </math>: waves propagating at a <math> 90^\circ\,</math> angle relative to the ship direction of forward translation.
  
either side is identical due to symmetry.
+
* <math> \theta = 35^\circ 16' \, </math>: (or 35,26°) waves propagating at an angle <math> \theta = 35^\circ 16' \, </math> relative to the ship axis. These are waves seen at the caustic of the Kelvin wake.
  
* <math> \theta = 0 \,</math>: waves propagating in the same direction as the ship. These waves
+
Let the solution of <math> \frac{y}{x}(\theta) \, </math> be of the form, when inverted:
  
can only exist at <math> Y = 0\, </math> as seen above.
+
<center><math> \mbox{Region I}: \qquad \theta = f_1 (\frac{y}{x}) \, </math></center>
  
* <math> \theta = \frac{\pi}{2} \, </math>: waves propagating at a <math> 90^\circ\,</math> angle
+
<center><math> \mbox{Region II}: \qquad \theta = f_2 (\frac{y}{x}) \, </math></center>
  
relative to the ship direction of forward translation.
+
Note that observable waves cannot exist for values of <math> \frac{y}{x}\,</math> that exceed the value shown in the figure or <math> \left. \frac{y}{x} \right|_{Max} = 2^{-3/2} \,</math>. This translates into a value for the corresponding angle equal to <math> 19^\circ 28' \, </math> (or 19,47°) which is the angle of the caustic for any speed <math> U\,</math>.
  
* <math> \theta = 35^\circ 16' \, </math>: waves propagating at an angle <math> \theta = 35^\circ
+
[[Image:Kelvin_wave_image.jpg|thumb|right|600px|"transverse" and "divergent" wave systems in the Kelvin wake]]
  
16' \, </math> relative to the ship axis. These are waves seen at the caustic of the Kelvin wake.
+
The crests of the wave system trailing a ship, the Kelvin wake, are curves of constant phase of:
  
Let the solution of the <math> \frac{y}{x})\theta \, </math> be of the form, when inverted:
+
<center><math> \frac{x\cos\theta+y\sin\theta}{\cos^2\theta} = \mathbf{C} \, </math></center>
  
<center><math> \mbox{Region} I: \qquad \theta = f_1 (\frac{y}{x})\, </math></center>
+
In <math> \mbox{Region I} \, </math>:
  
<center><math> \mbox{Region} II: \qquad \theta = f_2 {\frac{y}{x} \, </math></center>
+
<center><math> \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></center>
  
Note that observable waves cannot exist for values of <math> \frac{y}{x}\,</math> that exceed the
+
In <math> \mbox{Region II} \, </math>:
  
value shown in the figure or <math> \left. \frac{y}{x} \right|_Max = 2^{-3/2} \,</math>. This translates into a value for the corresponding angle equal to <math> 19^\circ 28' \, </math> which is the angle of the caustic for any speed <math> U\,</math>.
+
<center><math> \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></center>
  
The crests of the wave system trailing a ship, the Kelvin wake, are curves of constant phase of:
+
Plotting these curves we obtain a visual graph of the "transverse" and "divergent" wave systems in the Kelvin wake.
  
<center><math> \frac{x\cos\theta+y\sin\theta}{\cos^2\theta} = \mathbf{C} \, </math></center>
+
------
  
In <math> \mbox{Region} I \, </math>:
+
This article is based on the MIT open course notes and the original article can be found
 +
[http://ocw.mit.edu/NR/rdonlyres/Mechanical-Engineering/2-24Spring-2002/80D65339-8DD7-4909-A936-117BE9C97F5A/0/lecture8.pdf here]
  
<center><math> \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></center>
+
[[Ocean Wave Interaction with Ships and Offshore Energy Systems]]
  
In <math> \mbox{Region} II \, </math>:
+
[[Category:Linear Water-Wave Theory]]
 
 
<center><math> \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></center>
 
 
 
Plotting these curves we obtain a visual graph of the "transverse" and "divergent" wave systems in the Kelvin wake.
 

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