Difference between revisions of "Wave Energy Density and Flux"

From WikiWaves
Jump to navigationJump to search
(the context indicates that energy density should be lower case epsilon)
 
(68 intermediate revisions by 4 users not shown)
Line 5: Line 5:
 
}}
 
}}
  
 +
{{complete pages}}
  
= Introduction =
+
== Introduction ==
  
We are interested in the transport of energy by ocean waves. It is important to realise that under the assumptions of linear theory, there is not net motion of particles, but there is a transport of energy (as would be expected). The energy consists of two part, one kinetic due to the motion of the fluid and the other potential due to the the variation in the fluid height. It is the resonance between these two energies which gives rise to the wave motion. The situation is analagous to Simple Harmonic Motion but more complicated.  
+
We are interested in the transport of energy by ocean waves. It is important to realise that under the assumptions of linear theory, there is no net motion of particles, but there is a transport of energy (as would be expected). The energy consists of two parts, one kinetic due to the motion of the fluid and the other potential due to the variation in the fluid height. It is the resonance between these two energies which gives rise to the wave motion. The situation is analogous to  
 +
[http://en.wikipedia.org/wiki/Simple_harmonic_motion Simple Harmonic Motion] but more complicated.
  
= Energy per volume =
+
== Energy per volume ==
  
 
[[Image:Energy_volume.png|thumb|thumb|right|600px|Energy Volume]]
 
[[Image:Energy_volume.png|thumb|thumb|right|600px|Energy Volume]]
  
 
We begin by defining  
 
We begin by defining  
<math> \varepsilon(t) </math> as the energy in control volume <math> \gamma(t) </math> given by
+
<math> \mathcal{E}(t) </math> as the energy in control volume <math> \Omega(t) </math> given by
<center><math> \varepsilon (t) = \rho \iiint_\Omega \left( \frac{1}{2} |\mathbf{v}|^2 + gz \right) dV </math></center>
+
<center><math> \mathcal{E}(t) = \rho \iiint_\Omega \left( \frac{1}{2} |\mathbf{v}|^2 + gz \right) \mathrm{d}V </math></center>
where <math>\rho</math> is the fluid density, <math>g</math> is the acceleration due to gravity and <math>\mathbf{v}</math> is the vector
+
where <math>\rho</math> is the fluid density, <math>g</math> is the acceleration due to gravity and <math>\mathbf{v}</math> is the vector of fluid velocity.  
of fluid velocity.  
 
 
 
 
The mean energy over a unit horizontal surface area <math> S </math> is given by  
 
The mean energy over a unit horizontal surface area <math> S </math> is given by  
<center><math> \overline{\varepsilon} = \overline{\frac{\varepsilon(t)}{S}} = \rho \overline{ \int_{-h}^{\zeta(t)} \left( \frac{1}{2} |\mathbf{v}|^2 + gz \right) dz} = \frac{1}{2} \rho \overline{ \int_{-h}^{\zeta(t)} |\mathbf{v}|^2 dZ} + \overline{ \frac{1}{2} \rho g ( \zeta^2 - h^2 ) } </math></center>
+
<center><math> \overline{\mathcal{E}} = \overline{\frac{\mathcal{E}(t)}{S}} = \rho \overline{ \int_{-h}^{\zeta(t)} \left( \frac{1}{2} |\mathbf{v}|^2 + gz \right) \mathrm{d}z} = \frac{1}{2} \rho \overline{ \int_{-h}^{\zeta(t)} |\mathbf{v}|^2 \mathrm{d}z} + \overline{ \frac{1}{2} \rho g ( \zeta^2 - h^2 ) } </math></center>
 
where <math> \zeta(t) \, </math> is free surface elevation and the overbar denotes average (which will be important when we consider waves).  
 
where <math> \zeta(t) \, </math> is free surface elevation and the overbar denotes average (which will be important when we consider waves).  
 
Note that we are considering water of constant [[Finite Depth]].
 
Note that we are considering water of constant [[Finite Depth]].
Line 27: Line 27:
  
 
The remaining perturbation component is the sum of the kinetic and potential energy components, that is
 
The remaining perturbation component is the sum of the kinetic and potential energy components, that is
<center><math> \overline{\varepsilon} = \overline{\varepsilon_{kin}} + \overline{\varepsilon_{pot}} </math></center>
+
<center><math> \overline{\mathcal{E}} = \overline{\mathcal{E}_{kin}} + \overline{\mathcal{E}_{pot}} </math></center>
 
where
 
where
<center><math> \overline{\varepsilon_{kin}} = \frac{1}{2} \rho \overline{\int_{-h}^{\zeta(t)} |\mathbf{v}|^2 dz}, \qquad |\mathbf{v}|^2  
+
<center><math> \overline{\mathcal{E}_{kin}} = \frac{1}{2} \rho \overline{\int_{-h}^{\zeta(t)} |\mathbf{v}|^2 \mathrm{d}z}, \qquad |\mathbf{v}|^2  
 
= \nabla\Phi \cdot \nabla \Phi = \partial_x^2\Phi + \partial_z^2\Phi </math></center>
 
= \nabla\Phi \cdot \nabla \Phi = \partial_x^2\Phi + \partial_z^2\Phi </math></center>
 
and
 
and
<center><math> \overline{\varepsilon_{pot}} = \overline{\frac{1}{2} \rho g \zeta^2 (t)} </math></center>
+
<center><math> \overline{\mathcal{E}_{pot}} = \overline{\frac{1}{2} \rho g \zeta^2 (t)} </math></center>
 
Note that we are assuming only two dimensions <math>x</math> and <math>z</math>.
 
Note that we are assuming only two dimensions <math>x</math> and <math>z</math>.
  
= Energy in Linear Waves =
+
=== Energy in [[Linear Plane Progressive Regular Waves]] ===
  
 
Consider now as a special case of [[Linear Plane Progressive Regular Waves]] by the velocity potential in [[Infinite Depth]] water (for simplicity). The velocity potential throughout the fluid domain is then given by  
 
Consider now as a special case of [[Linear Plane Progressive Regular Waves]] by the velocity potential in [[Infinite Depth]] water (for simplicity). The velocity potential throughout the fluid domain is then given by  
<center><math> \Phi = \mathbf{Re} \{ \frac{igA}{\omega} e^{kz-ikx+i\omega t} \} </math></center>
+
<center><math> \Phi = \mathrm{Re} \{ \frac{igA}{\omega} e^{kz-ikx+i\omega t} \} </math></center>
 
The components is the <math>x</math> and <math>z</math> directions are given by  
 
The components is the <math>x</math> and <math>z</math> directions are given by  
<center><math> \Phi_x = \mathbf{Re} \{ \frac{igA}{\omega} (-ik) e^{kz-ikx+i\omega t} \} = A \mathbf{Re} \{ \omega e^{kz-ikx+i\omega t} \} </math></center>
+
<center><math> \partial_x\Phi = \mathrm{Re} \{ \frac{igA}{\omega} (-ik) e^{kz-ikx+i\omega t} \} = A \mathrm{Re} \{ \omega e^{kz-ikx+i\omega t} \} </math></center>
 
and
 
and
<center><math> \Phi_z = \mathbf{Re} \{ \frac{igA}{\omega} k e^{kz-ikx+i\omega t} \}  
+
<center><math> \partial_z\Phi = \mathrm{Re} \{ \frac{igA}{\omega} k e^{kz-ikx+i\omega t} \}  
  = A \mathbf{Re} \{ i \omega e^{kz-ikx+i\omega t} \} </math></center>
+
  = A \mathrm{Re} \{ i \omega e^{kz-ikx+i\omega t} \} </math></center>
 
respectively.  
 
respectively.  
 
 
We require the following <u>Lemma</u> which is easily proved. If
 
We require the following <u>Lemma</u> which is easily proved. If
<center><math> \mathbf{Re} \{ A e^{i\omega t} \} = A(t) </math></center>
+
<center><math> \mathrm{Re} \{ A e^{i\omega t} \} = A(t) \,\!</math></center>
 
and
 
and
<center><math> \mathbf{Re} \{ B e^{i\omega t} \} = B(t) </math></center>
+
<center><math> \mathrm{Re} \{ B e^{i\omega t} \} = B(t) \,\!</math></center>
 
then it follows that
 
then it follows that
<center><math> \overline{A(t)B(t)} = \frac{1}{2} \mathbf{Re} \{ A B^* \} </math></center>.
+
<center><math> \overline{A(t)B(t)} = \frac{1}{2} \mathrm{Re} \{ A B^* \}. </math></center>
  
 
This allows us to write the following expression  
 
This allows us to write the following expression  
<center><math> \overline{\epsilon_{kin}} = \frac{1}{2} \rho \overline{ ( \int_{-\infty}^0 + \int_0^\zeta ) \left( \Phi_x^2 + \Phi_z^2 \right) } dz </math></center>
+
<center><math>
 
+
\begin{align}
<center><math> = \frac{1}{2} \rho \int_{-\infty}^0 \left( \Phi_x^2 + \Phi_z^2 \right) dz + O (A^3) </math></center>
+
\overline{\mathcal{E}_{kin}} &= \frac{1}{2} \rho \overline{ ( \int_{-\infty}^0 + \int_0^\zeta )  
 
+
\left( (\partial_x\Phi)^2 + (\partial_z\Phi)^2 \right) } \mathrm{d}z \\
<center><math> = \rho \frac{\omega^2 A^2}{4k} = \frac{1}{4} \rho g A^2 , \qquad \mbox{for} \ k=\omega^2/g </math></center>
+
&= \frac{1}{2} \rho \int_{-\infty}^0 \overline{ \left( (\partial_x\Phi)^2 + (\partial_z\Phi)^2 \right) } \mathrm{d}z + O (A^3) \\
<center>
+
&= \rho \frac{\omega^2 A^2}{4k} = \frac{1}{4} \rho g A^2 , \qquad \mbox{for} \ k=\omega^2/g \\
<math> \overline{\varepsilon_{pot}} = \frac{1}{2} \rho g {\overline{\zeta(t)}}^2 = \frac{1}{4} \rho g A^2 </math>
+
\overline{\mathcal{E}_{pot}} &= \frac{1}{2} \rho g {\overline{\zeta(t)}}^2 = \frac{1}{4} \rho g A^2 .
</center>
+
\end{align}
 +
</math></center>
 
Note that it is a standard feature of linear oscillations that the average potential and kinetic energies are equal. Hence
 
Note that it is a standard feature of linear oscillations that the average potential and kinetic energies are equal. Hence
<center><math> \overline{\varepsilon} = \overline{\varepsilon_{kin}} + \overline{\varepsilon_{pot}} = \frac{1}{2} \rho g A^2 </math></center>
+
<center><math> \overline{\mathcal{E}} = \overline{\mathcal{E}_{kin}} + \overline{\mathcal{E}_{pot}} = \frac{1}{2} \rho g A^2 </math></center>
  
 
== Energy Flux ==
 
== Energy Flux ==
  
[[Image:Moving_volume.jpg|thumb|right|600px|Moving Volume]]
+
[[Image:volume_normal.png|thumb|right|600px|Moving Volume]]
  
<u> Energy flux </u> <math>E(t)</math> is the rate of change of energy density <math> \varepsilon(t) </math>. It is the flux of energy which is critical to ocean waves. While the individual fluid particles do not move the waves carry energy. We begin by deriving the energy flux in general conditions.
+
''Energy flux'' <math>\mathcal{P}(t)</math> is the rate of change of energy density <math> \epsilon(t) </math>. It is the flux of energy which is critical to ocean waves. While the individual fluid particles do not move the waves carry energy. We begin by deriving the energy flux in general conditions.
  
<center><math> E(t) \equiv \frac{d\varepsilon(t)}{dt} , \ \varepsilon = \iiint_{\Omega(t)} (\frac{1}{2} \rho |\mathbf{v}|^2 +gz ) \mathrm{d}V </math></center>
+
<center><math> \mathcal{P}(t) \equiv \frac{\mathrm{d}\mathcal{E}(t)}{\mathrm{d}t} = \frac{\mathrm{d}}{\mathrm{d}t} \iiint_{\Omega(t)} \epsilon(t) \mathrm{d}V  
 
 
<center><math> E(t) = \frac{\mathrm{d} \varepsilon(t)}{\mathrm{d}t} = \frac{\mathrm{d}}{\mathrm{d}t} \iiint_{\Omega(t)} \epsilon(t) \mathrm{d}V  
 
 
= \iint_{\partial\Omega(t)}  
 
= \iint_{\partial\Omega(t)}  
\frac{\partial \epsilon(t)}{\partial t} \mathrm{d}V + \iint_{\partial\Omega(t)} \epsilon(t) \mathbf{u}_n \mathrm{d} S </math></center>  
+
\frac{\partial \epsilon(t)}{\partial t} \mathrm{d}V + \iint_{\partial\Omega(t)} \epsilon(t) U_n \mathrm{d} S </math></center>  
 
+
(the last result following from  the transport theorem) where <math> U_n</math> is the normal velocity of surface <math> \partial\Omega(t) </math> outwards of the enclosed volume <math> \Omega(t)</math>.
Transport theorem where <math> \mathbf{u}_n</math> is the normal velocity of surface <math> \partial\Omega(t) </math> outwards of the enclosed volume <math> \Omega(t)</math>.
+
We know that
 
+
<center><math>  
<center><math> \frac{\partial \epsilon}{\partial t} = \frac{\partial}{\partial t} \{ \frac{1}{2} \rho |\mathbf{v}|^2 + \rho g z \}  
+
\begin{align}
= \frac{1}{2} \rho \frac{\partial}{\partial t} ( \nabla\Phi \cdot \nabla\Phi) </math></center>
+
\frac{\partial \epsilon}{\partial t} &= \frac{\partial}{\partial t} \{ \frac{1}{2} \rho |\mathbf{v}|^2 + \rho g z \} \\
 
+
&= \frac{1}{2} \rho \frac{\partial}{\partial t} ( \nabla\Phi \cdot \nabla\Phi) \\
<center><math> = \rho \nabla \cdot \left( \frac{\partial\Phi}{\partial t} \nabla\Phi \right) - \rho \frac{\partial\Phi}{\partial t} \nabla^2 \Phi </math></center>
+
&= \rho \nabla \cdot \left( \frac{\partial\Phi}{\partial t} \nabla\Phi \right) - \rho \frac{\partial\Phi}{\partial t} \nabla^2 \Phi  
 
+
\end{align}
<center><math> E(t) = \frac{d \varepsilon(t)}{dt} = \rho \iiint_{\Omega(t)} \nabla \cdot \left( \frac{\partial \Phi}{\partial t} \nabla \Phi \right) \mathrm{d}V  
+
</math></center>
+ \rho \iint_{\partial\Omega(t)} \left( \frac{1}{2} |\mathbf{v}|^2 + gz \right) \mathbf{u}_n \mathrm{d}S </math></center>
+
so that
 
+
<center><math> \mathcal{P}(t) = \rho \iiint_{\Omega(t)} \nabla \cdot \left( \frac{\partial \Phi}{\partial t} \nabla \Phi \right) \mathrm{d}V  
 +
+ \rho \iint_{\partial\Omega(t)} \left( \frac{1}{2} |\mathbf{v}|^2 + gz \right) U_n \mathrm{d}S </math></center>
 
Invoking the scalar form of Gauss's theorem in the first term, we obtain:
 
Invoking the scalar form of Gauss's theorem in the first term, we obtain:
 +
<center><math> \mathcal{P}(t) = \rho \iint_{\partial\Omega(t)} \frac{\partial\Phi}{\partial t} \nabla \Phi \cdot \mathbf{n} \mathrm{d}S
 +
+ \rho \iint_{\partial\Omega(t)} \left( \frac{1}{2} |\mathbf{v}|^2 + gz \right) U_n \mathrm{d}S </math></center>
 +
where <math>\mathbf{n}</math> is the unit normal.
  
<center><math> E(t) = \rho \iint_{\partial\Omega(t)} \frac{\partial\Phi}{\partial t} \nabla \Phi \cdot \mathbf{n} \mathrm{d}S
+
An alternative form for the energy flux <math> \mathcal{P}(t) \, </math> crossing the closed control surface <math> \partial\Omega(t) \, </math>  
+ \rho \iint_{\partial\Omega(t)} \left( \frac{1}{2} |\mathbf{v}|^2 + gz \right) \mathbf{u}_n \mathrm{d}S </math></center>
 
 
 
An alternative form for the energy flux <math> E(t) \, </math> crossing the closed control surface <math> \partial\Omega(t) \, </math>  
 
 
is obtained by invoking Bernoulli's equation in the second term. Recall that:
 
is obtained by invoking Bernoulli's equation in the second term. Recall that:
 
 
<center><math> \frac{P-P_a}{\rho} + \frac{\partial\Phi}{\partial t} + \frac{1}{2} \nabla\Phi \dot \nabla\Phi + gz = 0 </math></center>
 
<center><math> \frac{P-P_a}{\rho} + \frac{\partial\Phi}{\partial t} + \frac{1}{2} \nabla\Phi \dot \nabla\Phi + gz = 0 </math></center>
 
at any point in the fluid domain and on the boundary.
 
at any point in the fluid domain and on the boundary.
 
+
Here we allowed <math> \ P_a </math> the atmospheric pressure to be non-zero for the sake of physical clarity. Upon substitution in  
Here we allowed <math> \ P_a \equiv \mbox{Atmospheric pressure} \ </math> to be non-zero for the sake of physical clarity. Upon substitution in  
+
the equation above for <math> \mathcal{P}(t) </math> we obtain the alternate form:
the equation above for <math> E(t) </math> we obtain the alternate form:
+
<center><math> \mathcal{P}(t) = \rho \iint_{\partial\Omega(t)} \frac{\partial\Phi}{\partial t} \frac{\partial\Phi}{\partial n} \mathrm{d}s  
 
+
- \rho \iint_{\partial\Omega(t)} \left( \frac{P-P_a}{\rho} + \frac{\partial\Phi}{\partial t} \right) U_n \mathrm{d}s </math></center>
<center><math> E(t) = \rho \iint_{\partial\Omega(t)} \frac{\partial\Phi}{\partial t} \frac{\partial\Phi}{\partial n} \mathrm{d}s  
 
- \rho \iint_{\partial\Omega(t)} \left( \frac{P-P_a}{\rho} + \frac{\partial\Phi}{\partial t} \right) \mathbf{u}_n \mathrm{d}s </math></center>
 
 
where <math>\frac{\partial\phi}{\partial n}</math> is <math>\nabla\phi\cdot\mathbf{n}</math>  
 
where <math>\frac{\partial\phi}{\partial n}</math> is <math>\nabla\phi\cdot\mathbf{n}</math>  
 
 
So the energy flux across <math> \partial\Omega(t)\, </math> is given by the terms under the integral sign. They can be collected in the more compact form:
 
So the energy flux across <math> \partial\Omega(t)\, </math> is given by the terms under the integral sign. They can be collected in the more compact form:
<center><math> E(t) = \iint \left\{ \rho \frac{\partial\Phi}{\partial t} \left( \frac{\partial\Phi}{\partial n} - \mathbf{u}_n \right)  
+
<center><math> \mathcal{P}(t) = \iint \left\{ \rho \frac{\partial\Phi}{\partial t} \left( \frac{\partial\Phi}{\partial n} - U_n \right)  
- ( P - P_a) \mathbf{u}_n \right\} \mathrm{d}s </math></center>
+
- ( P - P_a) U_n \right\} \mathrm{d}s </math></center>
 +
Note that <math> \mathcal{P}(t) \, </math> measures the energy flux into the volume <math> \Omega(t) \, </math> or the rate of growth of the energy density <math> \mathcal{E}(t)\, </math>.
  
Note that <math> E(t) \, </math> measures the energy flux into the volume <math> \Omega(t) \, </math> or the rate of growth of the energy density <math> \varepsilon(t)\, </math>.
+
=== Energy transfer for each boundary ===
  
We are ready now to apply the above formulas to the surface wave propagation problem.
+
Break <math> \partial\Omega(t) \, </math> into its components and derive specialized forms of <math> \mathcal{P}(t) \, </math> pertinent to each.
  
Break <math> \partial\Omega(t) \, </math> into its components and derive specialized forms of <math> E(t) \, </math> pertinent to each.
+
* <math> \partial\Omega_F </math> nonlinear position of the free surface. On this <math> \partial_n\Phi= U_n; </math> and <math> P = P_a </math>, so the fluid pressure is equal to the atmospheric pressure. Therefore over <math> \partial\Omega_F </math>  <math>\mathcal{P}(t)=0</math>, as expected, i.e. there is no energy flow into the atmosphere.
 +
*<math> \partial\Omega_B</math> is the non-moving solid boundary, <math> U_n = </math> and <math> \frac{\partial\Phi}{\partial n} = U_n=0</math> which is the no-normal flux condition.
 +
*<math> \partial\Omega^\pm</math> which are the fluid boundaries fixed in space relative to an earth frame <math> U_n = 0 </math> and <math>, \frac{\partial\Phi}{\partial n} \ne 0</math>
 +
*<math>  \partial\Omega_U </math> the fluid boundaries moving with velocity <math>\mathbf{u}</math> relative to an earth frame.
 +
*<math> U_n = \mathbf{u} \cdot \mathbf{n}, \quad \frac{\partial\Phi}{\partial n} \ne 0 </math>. This case will be of interest for ships moving with constant velocity <math> U </math>.
  
<center><math> S_F : \ \mbox{nonlinear position of the free surface} </math></center>
+
The formula derived above are very general for potential flows with a free surface and solid boundaries. We are now ready to apply them to plane progressive waves.
  
<center><math> \frac{\partial\Phi}{\partial n} = U_n; \ \mbox{normal flow velocity} \equiv \mbox{normal velocity of free surface boundary; kinematic condition}. </math></center>
+
== Surface Wave Problem ==
  
<center><math> P = P_a; \ \mbox{fluid pressure} \equiv \mbox{atmospheric} </math></center>
+
We are ready now to apply the above formulas to the surface wave problem.
  
Therefore over <math> S_F; \ P(t) \equiv 0 </math> as expected. No energy can flow into the atmosphere!
+
=== Energy flux across a vertical fluid boundary fixed in space ===
  
<center><math> S_B: \ \mbox{non-moving solid boundary} </math></center>
+
<center><math>  
 
+
\begin{align}
<center><math> U_n = 0, \ \frac{\partial\Phi}{\partial n} = U_n; \ \mbox{no-normal flux condition} </math></center>
+
\mathcal{P}(t) &= - \rho \int_{-\infty}^{\zeta(t)} \frac{\partial\Phi}{\partial t} \Phi_n \mathrm{d}z \\
 
+
&= - \rho \left( \int_{-\infty}^0 + \int_0^\zeta \right) \frac{\partial\Phi}{\partial t} \Phi_n \mathrm{d}z \\
<center><math> S^\pm: \ \mbox{fluid boundaries fixed in space relative to an earth frame} </math></center>
+
&= - \rho \int_{-\infty}^0 \frac{\partial\Phi}{\partial t} \frac{\partial\Phi}{\partial n} \mathrm{d}z + O(A^3)
 
+
\end{align}
<center><math> U_n = 0, \ \frac{\partial\Phi}{\partial n} \ne 0 </math></center>
+
</math></center>
 
 
<center><math> S_U: \ \mbox{fluid boundaries moving w. velocity} \ \vec{U} \ \mbox{relative to an earth frame} </math></center>
 
 
 
<center><math> U_n = \vec U \cdot \vec n, \quad \frac{\partial\Phi}{\partial n} \ne 0 </math></center>
 
 
 
This case will be of interest for ships moving with constant velocity <math> U </math>.
 
 
 
The formulae derived above are very general for potential flows with a free surface and solid boundaries. We are now ready to apply them to plane progressive waves.
 
 
 
= Energy flux across a vertical fluid boundary fixed in space. =
 
 
 
<center><math> P(t) = - \rho \int_{-\infty}^{\zeta(t)} \frac{\partial\Phi}{\partial t} \Phi_n dZ = - \rho \left( \int_{-\infty}^0 + \int_0^\zeta \right) \frac{\partial\Phi}{\partial t} \Phi_n dZ = - \rho \int_{-\infty}^0 \frac{\partial\Phi}{\partial t} \frac{\partial\Phi}{\partial n} dZ + O(A^3) </math></center>
 
 
 
Mean energy flux for a plane progressive wave follows upon substitution of the regular wave velocity potential and taking mean values:
 
 
 
<center><math> \bar P = - \rho \int_{-\infty}^0 \frac{\partial\Phi}{\partial t} \frac{\partial\Phi}{\partial x} dZ = \frac{1}{2} \rho g A^2 \left( \frac{1}{2} \frac{g}{\omega} \right) </math></center>
 
  
 +
Mean energy flux for a [[Linear Plane Progressive Regular Waves]] follows upon substitution of the regular wave velocity potential and taking mean values:
 +
<center><math> \overline{\mathcal{P}} = - \rho \int_{-\infty}^0 \frac{\partial\Phi}{\partial t} \frac{\partial\Phi}{\partial x} \mathrm{d}z = \frac{1}{2} \rho g A^2 \left( \frac{1}{2} \frac{g}{\omega} \right) </math></center>
 +
where <math>A</math> is the wave amplitude.
 
or
 
or
 
+
<center><math> \overline{\mathcal{P}} = \overline{\mathcal{E}} c_g, \quad c_g = \mbox{group velocity} = \frac{1}{2} c_p </math></center>.
<center><math> \bar P = \bar \varepsilon V_g, \quad V_g = \mbox{group velocity} = \frac{1}{2} V_p \equiv \frac{1}{2} C </math></center>
+
It follows from this that the mean energy flux of a plane progressive wave is the product of its mean energy density times  
 
+
the [http://en.wikipedia.org/wiki/Group_velocity group velocity] of deep water waves.
It follows from this that the mean energy flux of a plane progressive wave is the product of its mean energy density times a velocity which equals <math> \frac{1}{2} </math> the phase velocity in deep water. We call this the [http://en.wikipedia.org/wiki/Group_velocity group velocity] of deep water waves and it is defined as
 
 
 
<center><math> V_g = \frac{1}{2} V_P = \frac{1}{2} \frac{g}{\omega} </math></center>
 
  
 
A more formal proof that this is the velocity with which the energy flux of plane progressive waves propagates is to consider what needs to be the horizontal velocity <math>U_n \equiv U </math> of a fluid boundary so that the mean energy flux across it vanishes?
 
A more formal proof that this is the velocity with which the energy flux of plane progressive waves propagates is to consider what needs to be the horizontal velocity <math>U_n \equiv U </math> of a fluid boundary so that the mean energy flux across it vanishes?
Line 159: Line 143:
 
This can be found from the solution of the following equation:
 
This can be found from the solution of the following equation:
  
<center><math> \overline{P(t)} = 0 = \rho \ {\overline{\int_{-\infty}^0 \frac{\partial\Phi}{\partial t} \frac{\partial\Phi}{\partial x} dZ}}^t - U \ {\overline{\int_{-\infty}^0 \left(\frac{P}{\rho} + \frac{\partial\Phi}{\partial t} \right) dZ}}^t = 0 </math></center>
+
<center><math> \overline{\mathcal{P}(t)} = \rho \ {\overline{\int_{-\infty}^0 \partial_t\Phi \partial_x\Phi \mathrm{d}z}} - U \ {\overline{\int_{-\infty}^0 \left(\frac{P}{\rho} + \frac{\partial\Phi}{\partial t} \right) \mathrm{d}z}} = 0 </math></center>
  
 
Where terms of <math> O(A^3) </math> have been neglected. Note that within linear theory, energy density and energy flux are quantities of <math>O(A^2)</math>. If higher-order terms are kept then we need to consider the treatment of second-order surface wave theory, at least.
 
Where terms of <math> O(A^3) </math> have been neglected. Note that within linear theory, energy density and energy flux are quantities of <math>O(A^2)</math>. If higher-order terms are kept then we need to consider the treatment of second-order surface wave theory, at least.
 
 
Solving the above equation for <math> U </math> we obtain:
 
Solving the above equation for <math> U </math> we obtain:
 
+
<center><math> U = \dfrac{\rho \ {\overline{\int_{-\infty}^0 \dfrac{\partial\Phi}{\partial t} \dfrac{\partial\Phi}{\partial x} \mathrm{d}z}}}
<center><math> U = \frac{\rho \ {\overline{\int_{-\infty}^0 \frac{\partial\Phi}{\partial t} \frac{\partial\Phi}{\partial x} dZ}}^t}{{\overline{\int_{-\infty}^0 \left( \frac{P}{\rho} + \frac{\partial\Phi}{\partial t} \right) dZ}}^t}  </math></center>
+
{{\overline{\int_{-\infty}^0 \left( \dfrac{P}{\rho} + \dfrac{\partial\Phi}{\partial t} \right) \mathrm{d}z}}}  </math></center>
  
 
Upon substitution of the plane progressive wave velocity potential and definition of pressure from Bernoulli's equation we obtain:
 
Upon substitution of the plane progressive wave velocity potential and definition of pressure from Bernoulli's equation we obtain:
 +
<center><math> U \equiv c_g = \frac{1}{2} \frac{g}{\omega} = \frac{1}{2} c_p </math></center>
 +
Note that <math> U \equiv c_g </math> by definition. If the above exercise is repeated in water of finite depth the solution for <math> U </math> after some algebra is:
  
<center><math> U \equiv V_g = \frac{1}{2} \frac{g}{\omega} = \frac{1}{2} V_P </math></center>
+
<center><math> U = c_g = \left( \frac{1}{2} + \frac{kh}{\sinh 2kh} \right) c_p </math></center>
 
 
Note that <math> U \equiv V_g </math> by definition. If the above exercise is repeated in water of finite depth the solution for <math> U </math> after some algebra is:
 
 
 
<center><math> U = V_g = \left( \frac{1}{2} + \frac{KH}{\sinh 2KH} \right) V_P </math></center>
 
  
 
with
 
with
  
<center><math> \omega^2 = gK \tanh KH \,</math></center>
+
<center><math> \omega^2 = gk \tanh kh\,</math></center>
  
It may be shown that the group velocity <math> V_g </math> is given in terms of <math> \omega \ne k \, </math> by the relation
+
It may be shown that the group velocity <math> c_g </math> is given by the relation
  
<center><math> V_g = \frac{d\omega}{d K} </math></center>
+
<center><math> c_g = \frac{\mathrm{d}\omega}{\mathrm{d} k} </math></center>
  
= Rayleigh's proof of the group velocity formula =
+
=== Rayleigh's proof of the group velocity formula ===
  
 
This relation follows from the very elegant "device" due to [http://en.wikipedia.org/wiki/John_Strutt%2C_3rd_Baron_Rayleigh Rayleigh] which applies to any wave form:
 
This relation follows from the very elegant "device" due to [http://en.wikipedia.org/wiki/John_Strutt%2C_3rd_Baron_Rayleigh Rayleigh] which applies to any wave form:
 
Consider two plane progressive waves of nearly equal frequencies and hence wavenumbers. Their joint wave elevation is given by
 
Consider two plane progressive waves of nearly equal frequencies and hence wavenumbers. Their joint wave elevation is given by
  
<center><math> \zeta(x,t) = A cos ( \omega_1 t - K_1 x) + A cos ( \omega_2 t - K_2 x) \, </math></center>
+
<center><math> \zeta(x,t) = A \cos ( \omega_1 t - k_1 x) + A \cos ( \omega_2 t - k_2 x) \, </math></center>
  
where the amplitude is assumed to be common and:
+
where the amplitude is assumed to be common and
  
<center><math> \omega_2 = \omega_1 + \Delta \omega, \quad | \Delta\omega | \ll \omega_1 , \omega_2 </math></center>
+
<center><math>  
 
+
\begin{align}
<center><math> K_2 = K_1 + \Delta K, \quad | \Delta K | \ll K_1 , K_2 </math></center>
+
\omega_2 &= \omega_1 + \Delta \omega, \quad | \Delta\omega | \ll \omega_1 , \omega_2 \\
 +
k_2 &= k_1 + \Delta k, \quad | \Delta k | \ll k_1 , k_2
 +
\end{align}
 +
</math></center>
  
 
Converting into complex notation:
 
Converting into complex notation:
  
<center><math> \zeta(x,t) = A \mathbf{Re} \{ e^{i\omega_1 t - i K_1 x} + e^{i\omega_2 t - i K_2 x} \} = A \mathbf{Re} \{ e^{i\omega_1 t - i K_1 x} + e^{i\omega_1 t - i K_1 x + i \Delta\omega t - i \Delta K x} \} </math></center> <br>
+
<center><math>  
<center><math> = A \mathbf{Re} \{ e^{i\omega_1 t - i K_1 x} \left( 1 + e^{i\Delta\omega t - i \Delta K x} \right) \} </math></center>
+
\begin{align}
 
+
\zeta(x,t) &= A \mathrm{Re} \{ e^{i\omega_1 t - i k_1 x} + e^{i\omega_2 t - i k_2 x} \} \\
The combined wave elevation <math> \zeta \,</math> vanishes identically where <math> F \equiv \left( 1 + e^{i\Delta\omega t - i \Delta K x} \right) = 0 \, </math>.
+
&= A \mathrm{Re} \{ e^{i\omega_1 t - i k_1 x} + e^{i\omega_1 t - i k_1 x + i \Delta\omega t - i \Delta k x} \} \\
 
+
&= A \mathrm{Re} \{ e^{i\omega_1 t - i k_1 x} \left( 1 + e^{i\Delta\omega t - i \Delta k x} \right) \}  
<math> F = 0 \, \ </math> when:
+
\end{align}
 
+
</math></center>
<center><math> e^{i(\Delta t - \Delta K x)} = -1 \,</math></center>
 
 
 
or when:
 
 
 
<center><math> \Delta \omega t - \Delta K x = ( 2 n + 1 ) \pi, \qquad n = 0, 1, 2, \cdots </math></center>
 
  
Solving for <math> x \, </math> we obtain:
+
The combined wave elevation <math> \zeta \,</math> vanishes identically where <math> \left( 1 + e^{i\Delta\omega t - i \Delta k x} \right) = 0 \, </math>,
 +
i.e.
 +
<center><math> e^{i(\Delta\omega t - \Delta k x)} = -1 \,</math></center>
 +
or when
 +
<center><math> \Delta \omega t - \Delta k x = ( 2 n + 1 ) \pi, \qquad n = 0, 1, 2, \cdots </math></center>
  
<center><math> x = \frac{1}{\Delta K} \{ (2n+1)\pi + t \Delta\omega \} \equiv X(t) </math></center>
+
Solving for <math> x(t) \, </math> we obtain:
  
For values of <math> X(t)\, </math> given above, <math> \zeta \equiv 0 \, </math>. These are the nodes of the bi-chromatic wave train where at all times the elevation vanishes and hence the evergy density <math> \equiv 0 </math>. The wave group has the form of consequtive packets separated by nodes.
+
<center><math> x(t) = \frac{1}{\Delta k} \{ (2n+1)\pi + t \Delta\omega \} </math></center>
 
 
The speed of the nodes is <math> \frac{dx}{dt} = \frac{\Delta\omega}{\Delta K} \to \frac{d\omega}{dK} \,  </math> and the energy trapped within two consecutive nodes cannot escape so it must travel at the group velocity: <math> V_g = \frac{dx}{dt} = \frac{d\omega}{dK} \, </math> !
 
 
 
Note that Rayleigh's proof applies equally to waves in finite depth or deep water and in principle to any propagating wave form!
 
  
 +
For values of <math> x(t)\, </math> given above, <math> \zeta = 0 \, </math>. These are the nodes of the bi-chromatic wave train where at all times the elevation vanishes and hence the evergy density is zero. The wave group has the form of consecutive packets separated by nodes.
 +
The speed of the nodes is <math> \frac{\mathrm{d}x}{\mathrm{d}t} = \frac{\Delta\omega}{\Delta k} \to \frac{\mathrm{d}\omega}{\mathrm{d}k} \,  </math>
 +
and the energy trapped within two consecutive nodes cannot escape so it must travel at the group velocity: <math> c_g = \frac{\mathrm{d}x}{\mathrm{d}t} = \frac{\mathrm{d}\omega}{\mathrm{d}k} \, </math>.
 +
Note that Rayleigh's proof applies equally to waves in finite depth or deep water and in principle to any propagating wave form.
 
In finite depth it can be shown after some algebra that
 
In finite depth it can be shown after some algebra that
 +
<center><math> c_g = \frac{\mathrm{d}\omega}{\mathrm{d}k} = \left( \frac{1}{2} + \frac{kh}{\sinh 2kh} \right) \frac{\omega}{k} </math></center>
  
<center><math> V_g = \frac{d\omega}{dK} = \left( \frac{1}{2} + \frac{KH}{\tanh KH} \right) \frac{\omega}{K} </math></center>
+
== Summary ==
 
 
= Summary =
 
  
 
The formulae for the energy flux derived above are very general and for potential flow nonlinear surface waves that are not breaking constitute the energy conservation principle.
 
The formulae for the energy flux derived above are very general and for potential flow nonlinear surface waves that are not breaking constitute the energy conservation principle.
Energy flux (power) input into the fluid domain by any mechanism, wavemaker wind (in a conservative manner), a ship or any floating body must be "retreived" at some distance away. Deriving expressions of the energy flux retreived at "infinity" is a powerful method for estimating the wave resistance of ships (more on this later), the wave damping of floating bodies etc.
+
Energy flux (power) input into the fluid domain by any mechanism, wavemaker, wind (in a conservative manner), a ship or any floating body must be "retrieved" at some distance away. Deriving expressions of the energy flux retrieved at "infinity" is a powerful method for estimating the wave resistance of ships (more on this later), the wave damping of floating bodies, etc.
 
 
 
Yet, the only general way of evaluating wave forces on floating bodies (moving or not) or solid boundaries is by applying the [[Wave Momentum Flux| Wave Momentum]].
 
Yet, the only general way of evaluating wave forces on floating bodies (moving or not) or solid boundaries is by applying the [[Wave Momentum Flux| Wave Momentum]].
  
 
-----
 
-----
 
This article is based on the MIT open course notes and the original article can be found
 
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/C2653597-0443-4705-A4E5-07A93AFCA3F5/0/lecture4.pdf here]
+
[http://ocw.mit.edu/courses/mechanical-engineering/2-24-ocean-wave-interaction-with-ships-and-offshore-energy-systems-13-022-spring-2002/lecture-notes/lecture4.pdf here]

Latest revision as of 17:44, 23 February 2011



Introduction

We are interested in the transport of energy by ocean waves. It is important to realise that under the assumptions of linear theory, there is no net motion of particles, but there is a transport of energy (as would be expected). The energy consists of two parts, one kinetic due to the motion of the fluid and the other potential due to the variation in the fluid height. It is the resonance between these two energies which gives rise to the wave motion. The situation is analogous to Simple Harmonic Motion but more complicated.

Energy per volume

Energy Volume

We begin by defining [math]\displaystyle{ \mathcal{E}(t) }[/math] as the energy in control volume [math]\displaystyle{ \Omega(t) }[/math] given by

[math]\displaystyle{ \mathcal{E}(t) = \rho \iiint_\Omega \left( \frac{1}{2} |\mathbf{v}|^2 + gz \right) \mathrm{d}V }[/math]

where [math]\displaystyle{ \rho }[/math] is the fluid density, [math]\displaystyle{ g }[/math] is the acceleration due to gravity and [math]\displaystyle{ \mathbf{v} }[/math] is the vector of fluid velocity. The mean energy over a unit horizontal surface area [math]\displaystyle{ S }[/math] is given by

[math]\displaystyle{ \overline{\mathcal{E}} = \overline{\frac{\mathcal{E}(t)}{S}} = \rho \overline{ \int_{-h}^{\zeta(t)} \left( \frac{1}{2} |\mathbf{v}|^2 + gz \right) \mathrm{d}z} = \frac{1}{2} \rho \overline{ \int_{-h}^{\zeta(t)} |\mathbf{v}|^2 \mathrm{d}z} + \overline{ \frac{1}{2} \rho g ( \zeta^2 - h^2 ) } }[/math]

where [math]\displaystyle{ \zeta(t) \, }[/math] is free surface elevation and the overbar denotes average (which will be important when we consider waves). Note that we are considering water of constant Finite Depth. We can ignore the term [math]\displaystyle{ -\frac{1}{2} \rho g h^2 \, }[/math] which represents the potential energy of the ocean at rest.

The remaining perturbation component is the sum of the kinetic and potential energy components, that is

[math]\displaystyle{ \overline{\mathcal{E}} = \overline{\mathcal{E}_{kin}} + \overline{\mathcal{E}_{pot}} }[/math]

where

[math]\displaystyle{ \overline{\mathcal{E}_{kin}} = \frac{1}{2} \rho \overline{\int_{-h}^{\zeta(t)} |\mathbf{v}|^2 \mathrm{d}z}, \qquad |\mathbf{v}|^2 = \nabla\Phi \cdot \nabla \Phi = \partial_x^2\Phi + \partial_z^2\Phi }[/math]

and

[math]\displaystyle{ \overline{\mathcal{E}_{pot}} = \overline{\frac{1}{2} \rho g \zeta^2 (t)} }[/math]

Note that we are assuming only two dimensions [math]\displaystyle{ x }[/math] and [math]\displaystyle{ z }[/math].

Energy in Linear Plane Progressive Regular Waves

Consider now as a special case of Linear Plane Progressive Regular Waves by the velocity potential in Infinite Depth water (for simplicity). The velocity potential throughout the fluid domain is then given by

[math]\displaystyle{ \Phi = \mathrm{Re} \{ \frac{igA}{\omega} e^{kz-ikx+i\omega t} \} }[/math]

The components is the [math]\displaystyle{ x }[/math] and [math]\displaystyle{ z }[/math] directions are given by

[math]\displaystyle{ \partial_x\Phi = \mathrm{Re} \{ \frac{igA}{\omega} (-ik) e^{kz-ikx+i\omega t} \} = A \mathrm{Re} \{ \omega e^{kz-ikx+i\omega t} \} }[/math]

and

[math]\displaystyle{ \partial_z\Phi = \mathrm{Re} \{ \frac{igA}{\omega} k e^{kz-ikx+i\omega t} \} = A \mathrm{Re} \{ i \omega e^{kz-ikx+i\omega t} \} }[/math]

respectively. We require the following Lemma which is easily proved. If

[math]\displaystyle{ \mathrm{Re} \{ A e^{i\omega t} \} = A(t) \,\! }[/math]

and

[math]\displaystyle{ \mathrm{Re} \{ B e^{i\omega t} \} = B(t) \,\! }[/math]

then it follows that

[math]\displaystyle{ \overline{A(t)B(t)} = \frac{1}{2} \mathrm{Re} \{ A B^* \}. }[/math]

This allows us to write the following expression

[math]\displaystyle{ \begin{align} \overline{\mathcal{E}_{kin}} &= \frac{1}{2} \rho \overline{ ( \int_{-\infty}^0 + \int_0^\zeta ) \left( (\partial_x\Phi)^2 + (\partial_z\Phi)^2 \right) } \mathrm{d}z \\ &= \frac{1}{2} \rho \int_{-\infty}^0 \overline{ \left( (\partial_x\Phi)^2 + (\partial_z\Phi)^2 \right) } \mathrm{d}z + O (A^3) \\ &= \rho \frac{\omega^2 A^2}{4k} = \frac{1}{4} \rho g A^2 , \qquad \mbox{for} \ k=\omega^2/g \\ \overline{\mathcal{E}_{pot}} &= \frac{1}{2} \rho g {\overline{\zeta(t)}}^2 = \frac{1}{4} \rho g A^2 . \end{align} }[/math]

Note that it is a standard feature of linear oscillations that the average potential and kinetic energies are equal. Hence

[math]\displaystyle{ \overline{\mathcal{E}} = \overline{\mathcal{E}_{kin}} + \overline{\mathcal{E}_{pot}} = \frac{1}{2} \rho g A^2 }[/math]

Energy Flux

Moving Volume

Energy flux [math]\displaystyle{ \mathcal{P}(t) }[/math] is the rate of change of energy density [math]\displaystyle{ \epsilon(t) }[/math]. It is the flux of energy which is critical to ocean waves. While the individual fluid particles do not move the waves carry energy. We begin by deriving the energy flux in general conditions.

[math]\displaystyle{ \mathcal{P}(t) \equiv \frac{\mathrm{d}\mathcal{E}(t)}{\mathrm{d}t} = \frac{\mathrm{d}}{\mathrm{d}t} \iiint_{\Omega(t)} \epsilon(t) \mathrm{d}V = \iint_{\partial\Omega(t)} \frac{\partial \epsilon(t)}{\partial t} \mathrm{d}V + \iint_{\partial\Omega(t)} \epsilon(t) U_n \mathrm{d} S }[/math]

(the last result following from the transport theorem) where [math]\displaystyle{ U_n }[/math] is the normal velocity of surface [math]\displaystyle{ \partial\Omega(t) }[/math] outwards of the enclosed volume [math]\displaystyle{ \Omega(t) }[/math]. We know that

[math]\displaystyle{ \begin{align} \frac{\partial \epsilon}{\partial t} &= \frac{\partial}{\partial t} \{ \frac{1}{2} \rho |\mathbf{v}|^2 + \rho g z \} \\ &= \frac{1}{2} \rho \frac{\partial}{\partial t} ( \nabla\Phi \cdot \nabla\Phi) \\ &= \rho \nabla \cdot \left( \frac{\partial\Phi}{\partial t} \nabla\Phi \right) - \rho \frac{\partial\Phi}{\partial t} \nabla^2 \Phi \end{align} }[/math]

so that

[math]\displaystyle{ \mathcal{P}(t) = \rho \iiint_{\Omega(t)} \nabla \cdot \left( \frac{\partial \Phi}{\partial t} \nabla \Phi \right) \mathrm{d}V + \rho \iint_{\partial\Omega(t)} \left( \frac{1}{2} |\mathbf{v}|^2 + gz \right) U_n \mathrm{d}S }[/math]

Invoking the scalar form of Gauss's theorem in the first term, we obtain:

[math]\displaystyle{ \mathcal{P}(t) = \rho \iint_{\partial\Omega(t)} \frac{\partial\Phi}{\partial t} \nabla \Phi \cdot \mathbf{n} \mathrm{d}S + \rho \iint_{\partial\Omega(t)} \left( \frac{1}{2} |\mathbf{v}|^2 + gz \right) U_n \mathrm{d}S }[/math]

where [math]\displaystyle{ \mathbf{n} }[/math] is the unit normal.

An alternative form for the energy flux [math]\displaystyle{ \mathcal{P}(t) \, }[/math] crossing the closed control surface [math]\displaystyle{ \partial\Omega(t) \, }[/math] is obtained by invoking Bernoulli's equation in the second term. Recall that:

[math]\displaystyle{ \frac{P-P_a}{\rho} + \frac{\partial\Phi}{\partial t} + \frac{1}{2} \nabla\Phi \dot \nabla\Phi + gz = 0 }[/math]

at any point in the fluid domain and on the boundary. Here we allowed [math]\displaystyle{ \ P_a }[/math] the atmospheric pressure to be non-zero for the sake of physical clarity. Upon substitution in the equation above for [math]\displaystyle{ \mathcal{P}(t) }[/math] we obtain the alternate form:

[math]\displaystyle{ \mathcal{P}(t) = \rho \iint_{\partial\Omega(t)} \frac{\partial\Phi}{\partial t} \frac{\partial\Phi}{\partial n} \mathrm{d}s - \rho \iint_{\partial\Omega(t)} \left( \frac{P-P_a}{\rho} + \frac{\partial\Phi}{\partial t} \right) U_n \mathrm{d}s }[/math]

where [math]\displaystyle{ \frac{\partial\phi}{\partial n} }[/math] is [math]\displaystyle{ \nabla\phi\cdot\mathbf{n} }[/math] So the energy flux across [math]\displaystyle{ \partial\Omega(t)\, }[/math] is given by the terms under the integral sign. They can be collected in the more compact form:

[math]\displaystyle{ \mathcal{P}(t) = \iint \left\{ \rho \frac{\partial\Phi}{\partial t} \left( \frac{\partial\Phi}{\partial n} - U_n \right) - ( P - P_a) U_n \right\} \mathrm{d}s }[/math]

Note that [math]\displaystyle{ \mathcal{P}(t) \, }[/math] measures the energy flux into the volume [math]\displaystyle{ \Omega(t) \, }[/math] or the rate of growth of the energy density [math]\displaystyle{ \mathcal{E}(t)\, }[/math].

Energy transfer for each boundary

Break [math]\displaystyle{ \partial\Omega(t) \, }[/math] into its components and derive specialized forms of [math]\displaystyle{ \mathcal{P}(t) \, }[/math] pertinent to each.

  • [math]\displaystyle{ \partial\Omega_F }[/math] nonlinear position of the free surface. On this [math]\displaystyle{ \partial_n\Phi= U_n; }[/math] and [math]\displaystyle{ P = P_a }[/math], so the fluid pressure is equal to the atmospheric pressure. Therefore over [math]\displaystyle{ \partial\Omega_F }[/math] [math]\displaystyle{ \mathcal{P}(t)=0 }[/math], as expected, i.e. there is no energy flow into the atmosphere.
  • [math]\displaystyle{ \partial\Omega_B }[/math] is the non-moving solid boundary, [math]\displaystyle{ U_n = }[/math] and [math]\displaystyle{ \frac{\partial\Phi}{\partial n} = U_n=0 }[/math] which is the no-normal flux condition.
  • [math]\displaystyle{ \partial\Omega^\pm }[/math] which are the fluid boundaries fixed in space relative to an earth frame [math]\displaystyle{ U_n = 0 }[/math] and [math]\displaystyle{ , \frac{\partial\Phi}{\partial n} \ne 0 }[/math]
  • [math]\displaystyle{ \partial\Omega_U }[/math] the fluid boundaries moving with velocity [math]\displaystyle{ \mathbf{u} }[/math] relative to an earth frame.
  • [math]\displaystyle{ U_n = \mathbf{u} \cdot \mathbf{n}, \quad \frac{\partial\Phi}{\partial n} \ne 0 }[/math]. This case will be of interest for ships moving with constant velocity [math]\displaystyle{ U }[/math].

The formula derived above are very general for potential flows with a free surface and solid boundaries. We are now ready to apply them to plane progressive waves.

Surface Wave Problem

We are ready now to apply the above formulas to the surface wave problem.

Energy flux across a vertical fluid boundary fixed in space

[math]\displaystyle{ \begin{align} \mathcal{P}(t) &= - \rho \int_{-\infty}^{\zeta(t)} \frac{\partial\Phi}{\partial t} \Phi_n \mathrm{d}z \\ &= - \rho \left( \int_{-\infty}^0 + \int_0^\zeta \right) \frac{\partial\Phi}{\partial t} \Phi_n \mathrm{d}z \\ &= - \rho \int_{-\infty}^0 \frac{\partial\Phi}{\partial t} \frac{\partial\Phi}{\partial n} \mathrm{d}z + O(A^3) \end{align} }[/math]

Mean energy flux for a Linear Plane Progressive Regular Waves follows upon substitution of the regular wave velocity potential and taking mean values:

[math]\displaystyle{ \overline{\mathcal{P}} = - \rho \int_{-\infty}^0 \frac{\partial\Phi}{\partial t} \frac{\partial\Phi}{\partial x} \mathrm{d}z = \frac{1}{2} \rho g A^2 \left( \frac{1}{2} \frac{g}{\omega} \right) }[/math]

where [math]\displaystyle{ A }[/math] is the wave amplitude. or

[math]\displaystyle{ \overline{\mathcal{P}} = \overline{\mathcal{E}} c_g, \quad c_g = \mbox{group velocity} = \frac{1}{2} c_p }[/math]

.

It follows from this that the mean energy flux of a plane progressive wave is the product of its mean energy density times the group velocity of deep water waves.

A more formal proof that this is the velocity with which the energy flux of plane progressive waves propagates is to consider what needs to be the horizontal velocity [math]\displaystyle{ U_n \equiv U }[/math] of a fluid boundary so that the mean energy flux across it vanishes?

This can be found from the solution of the following equation:

[math]\displaystyle{ \overline{\mathcal{P}(t)} = \rho \ {\overline{\int_{-\infty}^0 \partial_t\Phi \partial_x\Phi \mathrm{d}z}} - U \ {\overline{\int_{-\infty}^0 \left(\frac{P}{\rho} + \frac{\partial\Phi}{\partial t} \right) \mathrm{d}z}} = 0 }[/math]

Where terms of [math]\displaystyle{ O(A^3) }[/math] have been neglected. Note that within linear theory, energy density and energy flux are quantities of [math]\displaystyle{ O(A^2) }[/math]. If higher-order terms are kept then we need to consider the treatment of second-order surface wave theory, at least. Solving the above equation for [math]\displaystyle{ U }[/math] we obtain:

[math]\displaystyle{ U = \dfrac{\rho \ {\overline{\int_{-\infty}^0 \dfrac{\partial\Phi}{\partial t} \dfrac{\partial\Phi}{\partial x} \mathrm{d}z}}} {{\overline{\int_{-\infty}^0 \left( \dfrac{P}{\rho} + \dfrac{\partial\Phi}{\partial t} \right) \mathrm{d}z}}} }[/math]

Upon substitution of the plane progressive wave velocity potential and definition of pressure from Bernoulli's equation we obtain:

[math]\displaystyle{ U \equiv c_g = \frac{1}{2} \frac{g}{\omega} = \frac{1}{2} c_p }[/math]

Note that [math]\displaystyle{ U \equiv c_g }[/math] by definition. If the above exercise is repeated in water of finite depth the solution for [math]\displaystyle{ U }[/math] after some algebra is:

[math]\displaystyle{ U = c_g = \left( \frac{1}{2} + \frac{kh}{\sinh 2kh} \right) c_p }[/math]

with

[math]\displaystyle{ \omega^2 = gk \tanh kh\, }[/math]

It may be shown that the group velocity [math]\displaystyle{ c_g }[/math] is given by the relation

[math]\displaystyle{ c_g = \frac{\mathrm{d}\omega}{\mathrm{d} k} }[/math]

Rayleigh's proof of the group velocity formula

This relation follows from the very elegant "device" due to Rayleigh which applies to any wave form: Consider two plane progressive waves of nearly equal frequencies and hence wavenumbers. Their joint wave elevation is given by

[math]\displaystyle{ \zeta(x,t) = A \cos ( \omega_1 t - k_1 x) + A \cos ( \omega_2 t - k_2 x) \, }[/math]

where the amplitude is assumed to be common and

[math]\displaystyle{ \begin{align} \omega_2 &= \omega_1 + \Delta \omega, \quad | \Delta\omega | \ll \omega_1 , \omega_2 \\ k_2 &= k_1 + \Delta k, \quad | \Delta k | \ll k_1 , k_2 \end{align} }[/math]

Converting into complex notation:

[math]\displaystyle{ \begin{align} \zeta(x,t) &= A \mathrm{Re} \{ e^{i\omega_1 t - i k_1 x} + e^{i\omega_2 t - i k_2 x} \} \\ &= A \mathrm{Re} \{ e^{i\omega_1 t - i k_1 x} + e^{i\omega_1 t - i k_1 x + i \Delta\omega t - i \Delta k x} \} \\ &= A \mathrm{Re} \{ e^{i\omega_1 t - i k_1 x} \left( 1 + e^{i\Delta\omega t - i \Delta k x} \right) \} \end{align} }[/math]

The combined wave elevation [math]\displaystyle{ \zeta \, }[/math] vanishes identically where [math]\displaystyle{ \left( 1 + e^{i\Delta\omega t - i \Delta k x} \right) = 0 \, }[/math], i.e.

[math]\displaystyle{ e^{i(\Delta\omega t - \Delta k x)} = -1 \, }[/math]

or when

[math]\displaystyle{ \Delta \omega t - \Delta k x = ( 2 n + 1 ) \pi, \qquad n = 0, 1, 2, \cdots }[/math]

Solving for [math]\displaystyle{ x(t) \, }[/math] we obtain:

[math]\displaystyle{ x(t) = \frac{1}{\Delta k} \{ (2n+1)\pi + t \Delta\omega \} }[/math]

For values of [math]\displaystyle{ x(t)\, }[/math] given above, [math]\displaystyle{ \zeta = 0 \, }[/math]. These are the nodes of the bi-chromatic wave train where at all times the elevation vanishes and hence the evergy density is zero. The wave group has the form of consecutive packets separated by nodes. The speed of the nodes is [math]\displaystyle{ \frac{\mathrm{d}x}{\mathrm{d}t} = \frac{\Delta\omega}{\Delta k} \to \frac{\mathrm{d}\omega}{\mathrm{d}k} \, }[/math] and the energy trapped within two consecutive nodes cannot escape so it must travel at the group velocity: [math]\displaystyle{ c_g = \frac{\mathrm{d}x}{\mathrm{d}t} = \frac{\mathrm{d}\omega}{\mathrm{d}k} \, }[/math]. Note that Rayleigh's proof applies equally to waves in finite depth or deep water and in principle to any propagating wave form. In finite depth it can be shown after some algebra that

[math]\displaystyle{ c_g = \frac{\mathrm{d}\omega}{\mathrm{d}k} = \left( \frac{1}{2} + \frac{kh}{\sinh 2kh} \right) \frac{\omega}{k} }[/math]

Summary

The formulae for the energy flux derived above are very general and for potential flow nonlinear surface waves that are not breaking constitute the energy conservation principle. Energy flux (power) input into the fluid domain by any mechanism, wavemaker, wind (in a conservative manner), a ship or any floating body must be "retrieved" at some distance away. Deriving expressions of the energy flux retrieved at "infinity" is a powerful method for estimating the wave resistance of ships (more on this later), the wave damping of floating bodies, etc. Yet, the only general way of evaluating wave forces on floating bodies (moving or not) or solid boundaries is by applying the Wave Momentum.


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