Difference between revisions of "Wave Drift Forces"

From WikiWaves
Jump to navigationJump to search
Line 119: Line 119:
 
====Mean Drift Force on a Floating 2D Body====
 
====Mean Drift Force on a Floating 2D Body====
  
* Ambient regular waves with amplitude <math>A\,</math> are incident on a floating body in 2D which is allowed to oscillate freely in heave, sway and roll.
+
Ambient regular waves with amplitude <math>A\,</math> are incident on a floating body in 2D which is allowed to oscillate freely in heave, sway, and roll.
  
* As a result of the diffraction and radiation wave disturbances, a reflected plane progressive wave with complex amplitude <math>A^-\,</math> appears at <math>X=-\infty\,</math> over a control surface <math> S^-\,</math> and a radiated/diffracted wave with amplitude <math>A^+\,</math> appears at <math> X=+\infty\,</math> over <math>S^+\,</math>.
+
As a result of the diffraction and radiation wave disturbances, a reflected plane progressive wave with complex amplitude <math>A^-\,</math> appears at <math>x=-\infty\,</math> over a control surface <math> S^-\,</math> and a radiated/diffracted wave with amplitude <math>A^+\,</math> appears at <math> x=+\infty\,</math> over <math>S^+\,</math>.
  
* Applying the momentum conservation principle in the x-direction within the volume <math>V(t)\,</math> bounded by <math>S^+ + S^- + S_F(t) + S_B(t) \, </math> we obtain:
+
Applying the momentum conservation principle in the x-direction within the volume <math>V(t)\,</math> bounded by <math>S^+ + S^- + S_F(t) + S_B(t) \, </math> we obtain:
  
<center><math> D \equiv {\overline{F_X}}^t = - {\overline{\iint_{S^++S^-} \left[ P \vec{n}+\rho\vec{V}\left(V_n-U_n\right)\right]dS}}^t </math></center>
+
<center><math> D \equiv {\overline{F_x}}^t = - {\overline{\iint_{S^++S^-} \left[ P \vec{n}+\rho\vec{V}\left(V_n-U_n\right)\right]dS}}^t </math></center>
  
 
where the momentum conservation principle derived earlier has been applied. Mean values in time are taken leading to a vanishing mean value for the momentum rate of change within the volume <math>V(t)\,</math>, for fixed <math>S^\pm\,</math>.
 
where the momentum conservation principle derived earlier has been applied. Mean values in time are taken leading to a vanishing mean value for the momentum rate of change within the volume <math>V(t)\,</math>, for fixed <math>S^\pm\,</math>.
Line 131: Line 131:
 
* The mean momentum flux across <math>S_F(t)\,</math> is zero.
 
* The mean momentum flux across <math>S_F(t)\,</math> is zero.
  
* The mean horizontal momentum flux across the body section is the drift force <math>D\equiv{\overline{F_X}}^t\,</math>.
+
* The mean horizontal momentum flux across the body section is the drift force <math>D\equiv{\overline{F_x}}^t\,</math>.
  
 
* There remains to evaluate the integrals over <math>S^\pm\,</math> which are surfaces fixed in space and therefore, <math>U_n=0\,</math> on <math>S^\pm\,</math>. Recall that <math>\vec{V}=\nabla\phi\,</math> is the total fluid velocity and <math>P\,</math> is given by the Bernoulli equation.
 
* There remains to evaluate the integrals over <math>S^\pm\,</math> which are surfaces fixed in space and therefore, <math>U_n=0\,</math> on <math>S^\pm\,</math>. Recall that <math>\vec{V}=\nabla\phi\,</math> is the total fluid velocity and <math>P\,</math> is given by the Bernoulli equation.
Line 139: Line 139:
 
Define the velocity potentials for the plane progressive waves described above:
 
Define the velocity potentials for the plane progressive waves described above:
  
<center><math> \phi_I = \mathbb{R}\mathbf{e} \left\{ \frac{igA}{\omega} e^{KZ-iKX+i\omega t} \right\} \, </math></center>
+
<center><math> \phi_I = \mathbb{R}\mathbf{e} \left\{ \frac{igA}{\omega} e^{kz-ikx+i\omega t} \right\} \, </math></center>
<center><math> \phi^+ = \mathbb{R}\mathbf{e} \left\{ \frac{igA^+}{\omega} e^{KZ-iKX+i\omega t} \right\} \, </math></center>
+
<center><math> \phi^+ = \mathbb{R}\mathbf{e} \left\{ \frac{igA^+}{\omega} e^{kz-ikx+i\omega t} \right\} \, </math></center>
<center><math> \phi^- = \mathbb{R}\mathbf{e} \left\{ \frac{igA^-}{\omega} e^{KZ+iKX+i\omega t} \right\} \, </math></center>
+
<center><math> \phi^- = \mathbb{R}\mathbf{e} \left\{ \frac{igA^-}{\omega} e^{kz-ikx+i\omega t} \right\} \, </math></center>
  
 
Upon substitution into the momentum flux expression defined above and after some simple algebra that has been illustrated earlier, it follows that:
 
Upon substitution into the momentum flux expression defined above and after some simple algebra that has been illustrated earlier, it follows that:

Revision as of 15:02, 6 August 2010

Wave and Wave Body Interactions
Current Chapter Wave Drift Forces
Next Chapter Second-Order Wave Theory
Previous Chapter Panel Methods




Nonlinear Effects

Some of the most important nonlinear effects arising in connection with wave-body interactions are

  • Drift forces. They are the mean forces exerted on floating or submerged bodies by ambient waves. They may be treated very well by perturbation theory.
  • Slamming. These are highly nonlinear effects arising when a ship section acts upon the water surface or when a steep or breaking wave impinges upon a floating structure. They may be modeled by fully or partially nonlinear potential flow models of analytical or numerical nature
  • Forces due to viscous flow separation around floating structures and their subsystems, e.g. risers, mooring lines etc. Vortex induced vibrations (VIV) is an important example. Such effects can be treated experimentally and computationally by solving the Navier-Stokes equations.
  • Nonlinear ship motions in steep waves. These effects are mostly of potential flow nature and are being treated by nonlinear Rankine panel methods. The primary nonlinearity is the variable wetness of the ship hull, the nonlinearity of the kinematics of ambient waves and the numerical solution of the equations of motion in the time domain.
  • Nonlinear responses of deep water offshore platforms in certain flexural modes of their tethers. These effects are known as springing & ringing and are treated by a combination of perturbation and nonlinear methods and experiments.

Drift Forces

Drift forces will first be considered in regular waves. The main results will then be extended in random waves. One very important property of drift forces other than their practical significance is that they depend only on the linear solution.

Drift Forces by Pressure Integration

Mean Drift Force on a Vertical Wall

The incident and diffracted velocity potentials can be expressed as

[math]\displaystyle{ \phi_I = \mathrm{Re} \left\{ \frac{i g A}{\omega} e^{kz-ikx+i\omega t} \right\} \, }[/math]
[math]\displaystyle{ \phi_D = \mathrm{Re} \left\{ \frac{i g A}{\omega} e^{kz+ikx+i\omega t} \right\} \,. }[/math]

The total velocity potential is the sum of the incident and diffracted velocity potentials:

[math]\displaystyle{ \phi = \phi_I + \phi_D \,. }[/math]

On the wall, we have the following boundary condition as the velocity across the wall must be zero:

[math]\displaystyle{ \frac{\partial\phi}{\partial x} = 0, \; x = 0 \,. }[/math]

The nonlinear hydrodynamic pressure is expressed as

[math]\displaystyle{ P = - \rho \left( \frac{\partial\phi}{\partial t} + \frac{1}{2} \nabla\phi \cdot \nabla\phi + g z \right) \,. }[/math]

The nonlinear horizontal force on the wall is then

[math]\displaystyle{ F_x = \int_{-\infty}^\zeta P \mathrm{d}z, \quad \zeta = - \frac{1}{g} \frac{\partial\phi}{\partial t} \biggr\vert_{z=0} + O \left( z^2 \right) \,. }[/math]

We need to evaluate [math]\displaystyle{ F_x\, }[/math] correct to [math]\displaystyle{ O \left( A^2 \right) \, }[/math]. Note that the mean time value of effects of [math]\displaystyle{ O(A) \, }[/math], which are linear, is zero.

Assume a perturbation expansion for [math]\displaystyle{ \phi\, }[/math]:

[math]\displaystyle{ \phi = \phi_1 + \phi_2 + \phi_3 + \cdots \, }[/math]

where the potential derived earlier is the linear term denoted above by [math]\displaystyle{ \phi_1 \, }[/math].

In evaluating the leading order effect in the mean horizontal force we drop terms with zero mean values or of order [math]\displaystyle{ A^3\, }[/math] and higher. We note without proof that

[math]\displaystyle{ {\overline{\phi_2(t)}}^t = 0 \, }[/math]

The total horizontal force [math]\displaystyle{ F_x \, }[/math] also accepts the expansion:

[math]\displaystyle{ F_x = F_1 + F_2 + F_3 + \cdots \, }[/math]
[math]\displaystyle{ F_1 = \int_{-\infty}^0 P_1 \mathrm{d}z = - \rho \int_{-\infty}^0 \frac{\partial\phi}{\partial t} \mathrm{d}z = O (A) }[/math]
[math]\displaystyle{ {\overline{F_1(t)}}^t = 0 \, }[/math]
[math]\displaystyle{ F_2 = \int_{-\infty}^0 P_2 \mathrm{d}z + \int_0^\zeta P_1 \mathrm{d}z = O \left(A^2\right) }[/math]
[math]\displaystyle{ = - \rho \int_{-\infty}^0 \left( \frac{\partial\phi_2}{\partial t} + \frac{1}{2} \nabla\phi_1 \cdot \nabla\phi_1 \right) \mathrm{d}z }[/math]
[math]\displaystyle{ - \rho \int_0^\zeta \left( \frac{\partial\phi_1}{\partial t} + g z \right) \mathrm{d}z + O \left( A^3 \right) \, }[/math]

with errors of [math]\displaystyle{ O \left( A^3 \right) \, }[/math]. The last integral may be approximated by Taylor expanding about [math]\displaystyle{ z=0\, }[/math] the first term and by direct integration of the second. It follows that:

[math]\displaystyle{ \int_0^\zeta \left( \frac{\partial\phi_1}{\partial t} + g z \right) \mathrm{d}z = \zeta \left. \frac{\partial\phi_1}{\partial t} \right|_{z=0} + \frac{1}{2} g \zeta^2 + O \left( A^3 \right) }[/math]
[math]\displaystyle{ = - g \zeta^2 + \frac{1}{2} g \zeta^2 = - \frac{1}{2} g \zeta^2 \,. }[/math]

Collecting terms, we have:

[math]\displaystyle{ F_2(t) = - \rho \int_{-\infty}^0 \frac{\partial\phi_2}{\partial t} \mathrm{d}z - \frac{1}{2} \rho \int_{-\infty}^0 \nabla\phi_1 \cdot \nabla\phi_1 \mathrm{d}z + \frac{1}{2} \rho g \zeta^2(t) + O \left( A^3 \right) }[/math]

The mean time value [math]\displaystyle{ {\overline{\frac{\partial\phi_2}{\partial t}}}^t = 0 \, }[/math] for any stationary signal [math]\displaystyle{ \phi_2(t)\, }[/math]. So the second-order potential does not contribute to the mean drift force as stated above!

Of the remaining two terms the quadratic Bernoulli term contributes a suction force which is "pulling" the wall into the wave (counterintuitive but true!) while the last term is always positive pushing the wall in the direction of the wave as expected.

It follows that all of the mean drift force arises from the pressure integration over the surf-zone which is more than enough to overcome the suction force.

Upon substitution of the linear velocity potential derived above and use of the familiar identity:

[math]\displaystyle{ \overline{\mathrm{Re} \left(A_1 e^{i\omega t} \right) \mathrm{Re} \left(A_2 e^{i\omega t} \right)} = \frac{1}{2} \mathrm{Re} \left\{ A_1 A_2^* \right\} \,, }[/math]

it is easy to verify that

[math]\displaystyle{ {\overline{F_x}}^t = \frac{1}{2} \rho g A^2 + O \left( A^3 \right) \,. }[/math]

So the mean horizontal force on a vertical wall of infinite draft by a plane progressive wave has a finite mean value [math]\displaystyle{ \sim A^2 \, }[/math].

If the plane regular wave is incident at an angle [math]\displaystyle{ \beta }[/math], the mean horizontal force may be shown to take the value

[math]\displaystyle{ \overline{F_x} = \frac{1}{2} \rho g A^2 \cos^2\beta \,. }[/math]

In spite of their simplicity the above results have a number of useful applications in practice when waves that are sufficiently short interact with floating structures. Some examples follow:

  1. For a barge-like ship with draught [math]\displaystyle{ T }[/math], where [math]\displaystyle{ \frac{\lambda}{2} \lt T \Rightarrow \lambda \lt \frac{T}{2} \, }[/math], the ship with wall sided geometry acts like a wall. For this ship with length [math]\displaystyle{ L\, }[/math] the mean sway drift force in reasonably short waves is approximately given by
    [math]\displaystyle{ \overline{F_Y} \simeq \frac{1}{2} \rho g A^2 L \, }[/math]
  2. The above result can be extended to cases when the ship hull section has flare with slope [math]\displaystyle{ \alpha\, }[/math].
  3. Short waves incident upon a ship at an angle are locally reflected as if they encounter a continuum of vertical walls inclined at varying angles. We can thus apply the result derived above. Only part of the ship waterline will encounter waves. This region is can be determined by a simple geometrical argument. Denote this portion of the ship waterline by [math]\displaystyle{ \Gamma\, }[/math]. It is easy to show that the mean drift force in the [math]\displaystyle{ (x,y)\, }[/math] directions is given by
[math]\displaystyle{ \vec{F} = \left( \begin{matrix} \overline{F_x} \\ \overline{F_y} \end{matrix} \right) = \frac{1}{2} \rho g A^2 \int_{\Gamma} \left| \vec{n} \cdot \left( \begin{matrix} \cos\beta \\ \sin\beta \end{matrix} \right) \right|^2 \vec{n} dl \,. }[/math]

This result is known as Ray Theorem.

  1. Consider a vertical circular cylinder piercing the free surface and short waves incident upon it. The result from ray theorem stated above can be used to show that the mean drift force in the [math]\displaystyle{ x- }[/math]direction is given by
[math]\displaystyle{ \overline{F_x} = \frac{2}{3} \rho g A^2 R \,. }[/math]
  1. It is possible to extend ray theorem with little additional complexity to the case where the vessel advances with a moderate forward speed (see Faltinsen, 1990).

General case

In the more general case of surface waves interacting with floating bodies oscillating due to the ambient waves, a more general expression may be derived by directly integrating the pressure over the instantaneous position of the vessel wetted surface and linearizing about its mean position keeping consistently terms of [math]\displaystyle{ O \left( A^2 \right) }[/math]. The expression is not given here for simplicity.

Drift Forces by Momentum Conservation

It is possible to derive an expression for the mean drift force by applying the momentum conservation principle derived earlier. This approach provides an expression for the drift force in terms of the wave systems in the far field and can be shown to produce theoretically an identical result to that obtained from the considerably more elaborate task of using near-field pressure integration. The full three-dimensional result will not be presented here, as for the pressure integration, and can be found in Faltinsen 1990 and Wehausen and Laitone 1960. It is noted that it is numerically superior to pressure integration and is often preferred in practice.

Below we prove from first principles the result in two dimensions which is itself quite useful in practice and educational.

Mean Drift Force on a Floating 2D Body

Ambient regular waves with amplitude [math]\displaystyle{ A\, }[/math] are incident on a floating body in 2D which is allowed to oscillate freely in heave, sway, and roll.

As a result of the diffraction and radiation wave disturbances, a reflected plane progressive wave with complex amplitude [math]\displaystyle{ A^-\, }[/math] appears at [math]\displaystyle{ x=-\infty\, }[/math] over a control surface [math]\displaystyle{ S^-\, }[/math] and a radiated/diffracted wave with amplitude [math]\displaystyle{ A^+\, }[/math] appears at [math]\displaystyle{ x=+\infty\, }[/math] over [math]\displaystyle{ S^+\, }[/math].

Applying the momentum conservation principle in the x-direction within the volume [math]\displaystyle{ V(t)\, }[/math] bounded by [math]\displaystyle{ S^+ + S^- + S_F(t) + S_B(t) \, }[/math] we obtain:

[math]\displaystyle{ D \equiv {\overline{F_x}}^t = - {\overline{\iint_{S^++S^-} \left[ P \vec{n}+\rho\vec{V}\left(V_n-U_n\right)\right]dS}}^t }[/math]

where the momentum conservation principle derived earlier has been applied. Mean values in time are taken leading to a vanishing mean value for the momentum rate of change within the volume [math]\displaystyle{ V(t)\, }[/math], for fixed [math]\displaystyle{ S^\pm\, }[/math].

  • The mean momentum flux across [math]\displaystyle{ S_F(t)\, }[/math] is zero.
  • The mean horizontal momentum flux across the body section is the drift force [math]\displaystyle{ D\equiv{\overline{F_x}}^t\, }[/math].
  • There remains to evaluate the integrals over [math]\displaystyle{ S^\pm\, }[/math] which are surfaces fixed in space and therefore, [math]\displaystyle{ U_n=0\, }[/math] on [math]\displaystyle{ S^\pm\, }[/math]. Recall that [math]\displaystyle{ \vec{V}=\nabla\phi\, }[/math] is the total fluid velocity and [math]\displaystyle{ P\, }[/math] is given by the Bernoulli equation.
  • The definition of the plane progressive wave forms defined above will be introduced and an expression will be derived for [math]\displaystyle{ D\, }[/math] in terms of [math]\displaystyle{ A,A^+,A^-\, }[/math]. The result is "mildly" surprising!

Define the velocity potentials for the plane progressive waves described above:

[math]\displaystyle{ \phi_I = \mathbb{R}\mathbf{e} \left\{ \frac{igA}{\omega} e^{kz-ikx+i\omega t} \right\} \, }[/math]
[math]\displaystyle{ \phi^+ = \mathbb{R}\mathbf{e} \left\{ \frac{igA^+}{\omega} e^{kz-ikx+i\omega t} \right\} \, }[/math]
[math]\displaystyle{ \phi^- = \mathbb{R}\mathbf{e} \left\{ \frac{igA^-}{\omega} e^{kz-ikx+i\omega t} \right\} \, }[/math]

Upon substitution into the momentum flux expression defined above and after some simple algebra that has been illustrated earlier, it follows that:

[math]\displaystyle{ \mathbb{D} = \frac{\rho g A}{4} \left(A^+ +A^+* \right) - \frac{\rho g}{4} \left\{ \left|A^+\right|^2 - \left|A^-\right|^2 \right\} }[/math]

where [math]\displaystyle{ (*)\, }[/math] denotes the complex conjugate. Note that [math]\displaystyle{ A\, }[/math] is real and [math]\displaystyle{ A^\pm\, }[/math] are complex quantities.

  • Also note that the drift force as stated above does not appear to depend on [math]\displaystyle{ A\, }[/math] since the ambient wave appears at both infinities. The above expression simplifies a lot by invoking energy conservation.
  • The flow around the floating body is conservative if viscous effects or breaking wave effects are ignored. The former are often present when vortices are shed around corners of the body, bildge keels etc. In such cases the expressions below may be approximately valid or invalid.
  • In the absence of energy losses the mean energy flux across [math]\displaystyle{ S^-\, }[/math] must be identical to that across [math]\displaystyle{ S^+\, }[/math] within the limitations of linear theory:

[math]\displaystyle{ \rightarrow\, }[/math] Energy flux across [math]\displaystyle{ S^+\, }[/math]:

[math]\displaystyle{ P^+ = \frac{1}{4} \rho g V_P \left| A + A^+ \right|^2 \, }[/math]

[math]\displaystyle{ \rightarrow\, }[/math] Energy flux across [math]\displaystyle{ S^-\, }[/math]:

[math]\displaystyle{ P^- = \frac{1}{4} \rho g V_P \left( \left| A \right|^2 - \left| A^- \right|^2 \right) \, }[/math]
[math]\displaystyle{ P^+ \equiv P^- \Rightarrow\, }[/math]
[math]\displaystyle{ \left|A\right|^2 = \left|A^-\right|^2 + \left|A+A^+\right|^2 \, }[/math]
[math]\displaystyle{ = \left|A^-\right|^2 + \left|A\right|^2 + \left|A^+\right|^2 + A \left( A^+ + A^+* \right) }[/math]
[math]\displaystyle{ \Rightarrow \ A \left( A^+ + A^+* \right) = - \left( \left|A^+ \right|^2 + \left|A^- \right|^2 \right) }[/math]
  • This relation (restriction) upon the wave amplitudes is the result of energy conservation. Upon substitution in [math]\displaystyle{ \mathbb{D}\, }[/math]:
[math]\displaystyle{ \mathbb{D} = \frac{1}{2} \rho g \left|A^- \right|^2 }[/math]

This interesting result states that the mean drift force on a two-dimensional body freely floating in ambient regular waves (that sheds no vortices) is identically equal, according to linear theory, to the above expression that depends only on the amplitude of the reflected wave.

  • If the amplitude [math]\displaystyle{ \left|A^-\right|=0\, }[/math] then the body will experience zero mean drift force!

Submerged fixed circle

It can be shown that a submerged fixed circle does not reflect any of the ambient wave disturbance. Hence

[math]\displaystyle{ \left|A^-\right| = 0 \ \Rightarrow \ \mathbb{D} = 0 \ , }[/math]
  • This is the only case we know in 2D or 3D where the mean drift force is zero. If the circle is allowed to move then [math]\displaystyle{ \mathbb{D}\gt 0\, }[/math].
  • We know of no other shape in 2D with the same property. It is noteworthy that [math]\displaystyle{ \mathbb{D}\gt 0\, }[/math] in the case of a circle fixed in finite depth.
  • In practice the flow over the circle will separate and the above result is not valid. See OMF for a discussion of viscous effects around bluff bodies.


Ocean Wave Interaction with Ships and Offshore Energy Systems