Difference between revisions of "Category:Time-Dependent Linear Water Waves"

From WikiWaves
Jump to navigationJump to search
 
(One intermediate revision by the same user not shown)
Line 9: Line 9:
 
== Introduction ==
 
== Introduction ==
  
Analytical solutions of the linearised time-domain equations for a coupled-motion or forced-motion problems are very rare. In fact, only a small number of analytic time-domain solutions have been obtained (see the papers by [[Kennard 1949]] and [[McIver 1994]]) for structures with simple geometries.  Therefore, a variety of numerical methods have been developed to solve initial-value time-domain wave-structure interaction problems.  The Fourier transform enables the time and frequency domain quantities to be related and the frequency domain forces and potentials have important roles in the development of the numerical time-domain solution methods. This can be partly attributed to the fact that, in the past, analytical or semi-analytical solutions to linearised water-wave problems were more easily obtained in the frequency domain given that the assumption of time-harmonic motion results in a significant simplification of the interaction equations.  Therefore, it was natural to develop time-domain solutions based on frequency domain results and the Fourier transform relation between the domains.
+
Analytical solutions of the linearised time-domain equations for a coupled-motion or forced-motion problems are very rare. In fact, only a small number of analytic time-domain solutions have been obtained (see the papers by [[Kennard 1949]] and [[McIver 1994]]) for structures with simple geometries.  Therefore, a variety of numerical methods have been developed to solve initial-value time-domain wave-structure interaction problems.  The Fourier transform enables the time and frequency domain quantities to be related and the frequency domain forces and potentials have important roles in the development of the numerical time-domain solution methods. This can be partly attributed to the fact that, in the past, analytical or semi-analytical solutions to linearised water-wave problems were more easily obtained in the frequency domain given that the assumption of time-harmonic motion results in a significant simplification of the interaction equations.  Therefore, it was natural to develop time-domain solutions based on frequency domain results and the Fourier transform relation between the domains.
  
 +
[[Fourier Transform Solution for an Incident Wave Packet]]
  
 
== Equations of motion in the Time Domain ==
 
== Equations of motion in the Time Domain ==
Line 22: Line 23:
 
{{equations of motion time domain without body condition}}
 
{{equations of motion time domain without body condition}}
  
The body boundary condition for a fixed body is given in terms
+
{{two dimensional floating body time domain}}
of the 3 rigid body motions, namely surge, heave and pitch which are indexed as <math>\mu=1,3,5</math> in order to be consistent with the three-dimensional problem. We have a kinematic condition
 
<center><math>
 
\partial_{n}\Phi=\sum\partial_t \xi_{\nu}\mathbf{n}_{\nu},\ \mathbf{x}\in\partial\Omega,
 
</math></center>
 
where <math>\xi_{\nu}</math> is the motion of the <math>\mu</math>th mode and
 
<math>\mathbf{n}_{\nu}</math> is the normal associated with this mode.
 
The dynamic condition is the equation of motion for the structure:
 
<center><math>
 
M_{\mu\nu}\partial_t^2 \xi_{\nu}=-\rho\iint_{\partial\Omega}\partial_t\Phi n_{\mu}\, dS- c_{\mu\nu}\xi_{\nu},\quad \textrm{for} \qquad \mu=1,3,5,
 
</math></center>
 
where summation over repeated indices is assumed. In this equation, <math>M_{\mu\nu}</math> are the elements of the mass matrix
 
<center><math>
 
\mathbf{M}=\left[
 
\begin{matrix}
 
M        &      0      &  M(z^c-Z^R)    \\
 
0        &      M      & -M(x^c-X^R)    \\
 
M(z^c-Z^R)& -M(x^c-X^R)  & I^b_{11}+I^b_{33}
 
\end{matrix}
 
\right] , 
 
</math></center>
 
for the structure and <math>c_{\mu\nu}</math> are the elements of the buoyancy matrix
 
<center><math>
 
\mathbf{C}=\left[
 
\begin{matrix}
 
0        &      0        &  0    \\
 
0        &      \rho g W  &  -\rho g I^A_{1}  \\
 
0        &      -\rho g I^A_{1}  & \begin{matrix}\rho g (I^A_{11}+I^V_3)-\\
 
                                                    Mg(z^c-Z^R) \end{matrix}
 
\end{matrix}
 
\right].
 
</math></center>
 
The terms <math>I^b_{11}</math>, <math>I^b_{33}</math> are the moments of inertia of the body about the <math>x</math> and <math>z</math> axes and the terms <math>I_1^{A}</math>, <math>I^{A}_{11}</math> are the first and second moments of the waterplane (the waterplane area is denoted <math>W</math>) about the <math>x</math>-axis (see Chapter 7, [[Mei 1983]]). In addition, <math>(x^c,z^c)</math> and <math>(X^R,Z^R)</math> are the positions of the centre of mass and centre of rotation of the body and <math>I^{V}_{3}</math> is <math>z</math>-component centre of buoyancy of the structure. Thus, the coupled equations of motion for a floating structure have been derived. (N.B. if is assumed that the centre of rotation and the centre of mass of the structure coincide, i.e. if it is assumed that the body is semi-submerged, the mass and buoyancy matrices become diagonal).  Any wave incidence is assumed to be propagating in the positive <math>x</math> direction.) The scattering and radiation problems are simpler than the coupled problem because the motion of the the structure is then prescribed.
 
  
 
The initial conditions are  
 
The initial conditions are  

Latest revision as of 02:24, 7 April 2010



Generally the focus of research is on the Frequency Domain Problem. The time-domain problem can be solved by Generalized Eigenfunction Expansion for Water Waves or by using Memory Effect Function or by the Laplace Transform for Water Waves

Introduction

Analytical solutions of the linearised time-domain equations for a coupled-motion or forced-motion problems are very rare. In fact, only a small number of analytic time-domain solutions have been obtained (see the papers by Kennard 1949 and McIver 1994) for structures with simple geometries. Therefore, a variety of numerical methods have been developed to solve initial-value time-domain wave-structure interaction problems. The Fourier transform enables the time and frequency domain quantities to be related and the frequency domain forces and potentials have important roles in the development of the numerical time-domain solution methods. This can be partly attributed to the fact that, in the past, analytical or semi-analytical solutions to linearised water-wave problems were more easily obtained in the frequency domain given that the assumption of time-harmonic motion results in a significant simplification of the interaction equations. Therefore, it was natural to develop time-domain solutions based on frequency domain results and the Fourier transform relation between the domains.

Fourier Transform Solution for an Incident Wave Packet

Equations of motion in the Time Domain

Two Dimensional Equations for fixed bodies in the time domain

We consider a two-dimensional fluid domain of constant depth, which contains a finite number of fixed bodies of arbitrary geometry. We denote the fluid domain by [math]\displaystyle{ \Omega }[/math], the boundary of the fluid domain which touches the fixed bodies by [math]\displaystyle{ \partial\Omega }[/math], and the free surface by [math]\displaystyle{ F. }[/math] The [math]\displaystyle{ x }[/math] and [math]\displaystyle{ z }[/math] coordinates are such that [math]\displaystyle{ x }[/math] is pointing in the horizontal direction and [math]\displaystyle{ z }[/math] is pointing in the vertical upwards direction (we denote [math]\displaystyle{ \mathbf{x}=\left( x,z\right) ). }[/math] The free surface is at [math]\displaystyle{ z=0 }[/math] and the sea floor is at [math]\displaystyle{ z=-h }[/math]. The fluid motion is described by a velocity potential [math]\displaystyle{ \Phi }[/math] and free surface by [math]\displaystyle{ \zeta }[/math].

The equations of motion in the time domain are Laplace's equation through out the fluid

[math]\displaystyle{ \Delta\Phi\left( \mathbf{x,}t\right) =0,\ \ \mathbf{x}\in\Omega. }[/math]

At the bottom surface we have no flow

[math]\displaystyle{ \partial_{n}\Phi=0,\ \ z=-h. }[/math]

At the free surface we have the kinematic condition

[math]\displaystyle{ \partial_{t}\zeta=\partial_{n}\Phi,\ \ z=0,\ x\in \partial\Omega_{F}, }[/math]

and the dynamic condition (the linearized Bernoulli equation)

[math]\displaystyle{ \partial_{t}\Phi = -g\zeta ,\ \ z=0,\ x\in \partial\Omega_{F}. }[/math]

The body boundary condition for a fixed body is

[math]\displaystyle{ \partial_{n}\Phi=0,\ \ \mathbf{x}\in\partial\Omega, }[/math]


The initial conditions are

[math]\displaystyle{ \left.\zeta\right|_{t=0} = \zeta_0(x)\,\,\,\mathrm{and}\,\,\, \left.\partial_t\zeta\right|_{t=0} = \partial_t\zeta_0(x). }[/math]

Two dimensional equations for a floating body

We now consider the equations for a floating two-dimensional structure. The equations of motion in the time domain are Laplace's equation through out the fluid

[math]\displaystyle{ \Delta\Phi\left( \mathbf{x,}t\right) =0,\ \ \mathbf{x}\in\Omega. }[/math]

At the bottom surface we have no flow

[math]\displaystyle{ \partial_{n}\Phi=0,\ \ z=-h. }[/math]

At the free surface we have the kinematic condition

[math]\displaystyle{ \partial_{t}\zeta=\partial_{n}\Phi,\ \ z=0,\ x\in \partial\Omega_{F}, }[/math]

and the dynamic condition (the linearized Bernoulli equation)

[math]\displaystyle{ \partial_{t}\Phi = -g\zeta ,\ \ z=0,\ x\in \partial\Omega_{F}. }[/math]

The body boundary condition for a floating body is given in terms of the 3 rigid body motions, namely surge, heave and pitch which are indexed as [math]\displaystyle{ \nu=1,3,5 }[/math] in order to be consistent with the three-dimensional problem. We have a kinematic condition

[math]\displaystyle{ \partial_{n}\Phi=\sum_{\nu}\partial_t \xi_{\nu}\mathbf{n}_{\nu},\ \mathbf{x}\in\partial\Omega_{B}, }[/math]

where [math]\displaystyle{ \xi_{\nu} }[/math] is the motion of the [math]\displaystyle{ \mu }[/math]th mode and [math]\displaystyle{ \mathbf{n}_{\nu} }[/math] is the normal associated with this mode. Note that we define all normal derivatives to point out of the fluid. The dynamic condition is the equation of motion for the structure:

[math]\displaystyle{ \sum_{\nu} M_{\mu\nu}\partial_t^2 \xi_{\nu}=-\rho\iint_{\partial\Omega_{B}}\partial_t\Phi n_{\mu}\, dS - \sum_{\nu} C_{\mu\nu}\xi_{\nu},\quad \textrm{for} \qquad \mu=1,3,5, }[/math]

In this equation, [math]\displaystyle{ M_{\mu\nu} }[/math] are the elements of the mass matrix

[math]\displaystyle{ \mathbf{M}=\left[ \begin{matrix} M & 0 & M(z^c-Z^R) \\ 0 & M & -M(x^c-X^R) \\ M(z^c-Z^R)& -M(x^c-X^R) & I^b_{11}+I^b_{33} \end{matrix} \right] , }[/math]

for the structure and [math]\displaystyle{ c_{\mu\nu} }[/math] are the elements of the buoyancy matrix

[math]\displaystyle{ \mathbf{C}=\left[ \begin{matrix} 0 & 0 & 0 \\ 0 & \rho g W & -\rho g I^A_{1} \\ 0 & -\rho g I^A_{1} & \rho g (I^A_{11}+I^V_3)-Mg(z^c-Z^R) \end{matrix} \right]. }[/math]

The terms [math]\displaystyle{ I^b_{11} }[/math], [math]\displaystyle{ I^b_{33} }[/math] are the moments of inertia of the body about the [math]\displaystyle{ x }[/math] and [math]\displaystyle{ z }[/math] axes and the terms [math]\displaystyle{ I_1^{A} }[/math], [math]\displaystyle{ I^{A}_{11} }[/math] are the first and second moments of the waterplane (the waterplane area is denoted [math]\displaystyle{ W }[/math]) about the [math]\displaystyle{ x }[/math]-axis (see Chapter 7, Mei 1983). In addition, [math]\displaystyle{ (x^c,z^c) }[/math] and [math]\displaystyle{ (X^R,Z^R) }[/math] are the positions of the centre of mass and centre of rotation of the body and [math]\displaystyle{ I^{V}_{3} }[/math] is [math]\displaystyle{ z }[/math]-component centre of buoyancy of the structure. Thus, the coupled equations of motion for a floating structure have been derived. (N.B. if is assumed that the centre of rotation and the centre of mass of the structure coincide, i.e. if it is assumed that the body is semi-submerged, the mass and buoyancy matrices become diagonal). Any wave incidence is assumed to be propagating in the positive [math]\displaystyle{ x }[/math] direction.) The scattering and radiation problems are simpler than the coupled problem because the motion of the the structure is then prescribed.

The initial conditions are

[math]\displaystyle{ \left.\zeta\right|_{t=0} = \zeta_0(x)\,\,\,\mathrm{and}\,\,\, \left.\partial_t\zeta\right|_{t=0} = \partial_t\zeta_0(x). }[/math]

and The initial generalised displacements [math]\displaystyle{ \xi_{\mu} }[/math] and velocities [math]\displaystyle{ \dot{\xi}_{\mu} }[/math] of the body must be specified for all modes [math]\displaystyle{ \mu }[/math]

Floating body constrained to move in heave

The simplest type of floating body problem concerns the motion of a body constrained to move in heave in two-dimensions. Apart from the boundary condition on the structure surface which becomes

[math]\displaystyle{ \partial_n \Phi = \partial_t\xi_{3}\mathbf{n}_{3}, }[/math]

where [math]\displaystyle{ \xi_{3} }[/math] is the heave displacement of the structure, the equations governing the motion of the fluid remain the same. The equation of motion for a body constrained to move in heave only is

[math]\displaystyle{ M\partial_t^2 \xi_{3}=-\rho\iint_{\partial\Omega_{\mathrm{B}}}\partial_t\Phi n_{3}\, dS- \rho g W\xi_{3}. }[/math]

which is a significant simplification compared to the general case - the mass matrix and buoyancy matrix are replaced simply by the mass [math]\displaystyle{ M }[/math] and the hydrostatic term [math]\displaystyle{ \rho g W }[/math]. The initial conditions for the fluid and the structure ([math]\displaystyle{ \xi_{3}(0) }[/math], [math]\displaystyle{ \partial_t{\xi}_{3}(0) }[/math]) must also be prescribed to complete the problem specification.

Pages in category "Time-Dependent Linear Water Waves"

The following 3 pages are in this category, out of 3 total.