Difference between revisions of "Standard Linear Wave Scattering Problem"

From WikiWaves
Jump to navigationJump to search
Line 4: Line 4:
 
The water motion is represented by a velocity potential which is
 
The water motion is represented by a velocity potential which is
 
denoted by <math>\phi</math>.  The coordinate system is the standard Cartesian coordinate system
 
denoted by <math>\phi</math>.  The coordinate system is the standard Cartesian coordinate system
with the <math>z</math> axis pointing vertically up. The water occupies the region
+
with the <math>z</math> axis pointing vertically up. The water surface is at
<math>-h<z<0.</math> We denote the free surface by <math>\Gamma_s</math> (located at <math>z=0</math>)
+
<math>z=0</math> and the region of interest is
and the wetted
+
<math>-h<z<0</math>. There is a body which occupies the region <math>\Omega</math>
surface of the ice floe by <math>\Gamma_w</math>. We denote <math>\mathbf{r}= x</math>
+
and we denoted the wetted surface of the body by <math>\partial\Omega</math>
or <math>\mathbf{r}=(x,y)</math> as the horizontal coordinate in two or three dimensions
+
We denote <math>\mathbf{r}=(x,y)</math> as the horizontal coordinate in two or three dimensions
respectively.  
+
respectively and the cartesian system we denote by <math>\mathbf{x}</math>.
 +
We assume that the bottom surface is of constant depth but [[Variable Bottom Topography]]
 +
can easily be included.  
  
 
The equations are the following
 
The equations are the following
 
+
<center><math>
<math>
+
\nabla^{2}\phi=0, \, -\infty<z<0,\,\,\,\mathbf{x}\notin \Omega
\nabla^{2}\phi=0, \, -\infty<z<0
+
</math></center>
</math>
+
<center><math>
 
 
<math>
 
 
\frac{\partial\phi}{\partial z}=0, \, z=h,
 
\frac{\partial\phi}{\partial z}=0, \, z=h,
</math>
+
</math></center>
 
+
<center><math>
<math>
+
  \frac{\partial\phi}{\partial z} = k_{\infty}\phi,\,\mathbf{x}\in\partial\Omega,
  \frac{\partial\phi}{\partial z} = k_{\infty}\phi,\,z\in\Gamma_s,
+
</math></center>
</math>
+
<center><math>
 
 
[[Category:Linear Water-Wave Theory]]
 
<math>
 
 
  \frac{\partial\phi}{\partial z} = L\phi, \, z\in\Gamma_w.
 
  \frac{\partial\phi}{\partial z} = L\phi, \, z\in\Gamma_w.
</math>
+
</math></center>
  
 
where <math>k_{\infty}</math> is the wavenumber in [[Infinite Depth]] which is given by  
 
where <math>k_{\infty}</math> is the wavenumber in [[Infinite Depth]] which is given by  
 
<math>k_{\infty}=\omega^2/g</math> where <math>g</math> is gravity. <math>L</math> is a linear
 
<math>k_{\infty}=\omega^2/g</math> where <math>g</math> is gravity. <math>L</math> is a linear
 
operator which relates the normal and potential on the body surface through the physics
 
operator which relates the normal and potential on the body surface through the physics
of the body. For a fixed body the operator <math>L=0</math>.
+
of the body. The simplest case is for a fixed body  
 +
where the operator is <math>L=0</math> but more complicated conditions are possible.
  
 
The equation is subject to some radiation conditions at infinity. We usually assume that
 
The equation is subject to some radiation conditions at infinity. We usually assume that
 
there is an incident wave <math>\phi^{\mathrm{{In}}}\,</math>   
 
there is an incident wave <math>\phi^{\mathrm{{In}}}\,</math>   
 
is a plane wave travelling in the <math>x</math> direction  
 
is a plane wave travelling in the <math>x</math> direction  
 
+
<center><math>
<math>
 
 
\phi^{\mathrm{{In}}}({r},z)=Ae^{{\rm i}kx}\cosh k(z+h)\,
 
\phi^{\mathrm{{In}}}({r},z)=Ae^{{\rm i}kx}\cosh k(z+h)\,
</math>
+
</math></center>
 
 
 
where <math>A</math> is the wave amplitude and <math>k</math> is the wavenumber which is
 
where <math>A</math> is the wave amplitude and <math>k</math> is the wavenumber which is
 
the positive real solution of the [[Dispersion Relation for a Free Surface]].
 
the positive real solution of the [[Dispersion Relation for a Free Surface]].
Line 49: Line 45:
  
 
In two-dimensions the condition is  
 
In two-dimensions the condition is  
 
+
<center><math>
<math>
 
 
\left(  \frac{\partial}{\partial|x|}-{i}k\right)
 
\left(  \frac{\partial}{\partial|x|}-{i}k\right)
 
(\phi-\phi^{\mathrm{{In}}})=0,\;\mathrm{{as\;}}|x|\rightarrow\infty\mathrm{.}
 
(\phi-\phi^{\mathrm{{In}}})=0,\;\mathrm{{as\;}}|x|\rightarrow\infty\mathrm{.}
</math>
+
</math></center>
 
 
 
where <math>\phi^{\mathrm{{In}}}</math> is the incident potential and <math>k</math>
 
where <math>\phi^{\mathrm{{In}}}</math> is the incident potential and <math>k</math>
 
is the wave number.
 
is the wave number.
  
 
In three-dimensions the condition is  
 
In three-dimensions the condition is  
 
+
<center><math>
<math>
 
 
\sqrt{|\mathbf{r}|}\left(  \frac{\partial}{\partial|\mathbf{r}|}-{i}k\right)
 
\sqrt{|\mathbf{r}|}\left(  \frac{\partial}{\partial|\mathbf{r}|}-{i}k\right)
 
(\phi-\phi^{\mathrm{{In}}})=0,\;\mathrm{{as\;}}|\mathbf{r}|\rightarrow\infty\mathrm{.}
 
(\phi-\phi^{\mathrm{{In}}})=0,\;\mathrm{{as\;}}|\mathbf{r}|\rightarrow\infty\mathrm{.}
</math>
+
</math></center>
 +
 
 +
[[Category:Linear Water-Wave Theory]]

Revision as of 03:48, 21 August 2006

The standard linear wave scattering problem is based on assuming irrotational and inviscid fluid motion and the wave sufficiently small amplitude so that we can linearise all the equations. We also assume that Frequency Domain Problem with frequency [math]\displaystyle{ \omega }[/math] The water motion is represented by a velocity potential which is denoted by [math]\displaystyle{ \phi }[/math]. The coordinate system is the standard Cartesian coordinate system with the [math]\displaystyle{ z }[/math] axis pointing vertically up. The water surface is at [math]\displaystyle{ z=0 }[/math] and the region of interest is [math]\displaystyle{ -h\lt z\lt 0 }[/math]. There is a body which occupies the region [math]\displaystyle{ \Omega }[/math] and we denoted the wetted surface of the body by [math]\displaystyle{ \partial\Omega }[/math] We denote [math]\displaystyle{ \mathbf{r}=(x,y) }[/math] as the horizontal coordinate in two or three dimensions respectively and the cartesian system we denote by [math]\displaystyle{ \mathbf{x} }[/math]. We assume that the bottom surface is of constant depth but Variable Bottom Topography can easily be included.

The equations are the following

[math]\displaystyle{ \nabla^{2}\phi=0, \, -\infty\lt z\lt 0,\,\,\,\mathbf{x}\notin \Omega }[/math]
[math]\displaystyle{ \frac{\partial\phi}{\partial z}=0, \, z=h, }[/math]
[math]\displaystyle{ \frac{\partial\phi}{\partial z} = k_{\infty}\phi,\,\mathbf{x}\in\partial\Omega, }[/math]
[math]\displaystyle{ \frac{\partial\phi}{\partial z} = L\phi, \, z\in\Gamma_w. }[/math]

where [math]\displaystyle{ k_{\infty} }[/math] is the wavenumber in Infinite Depth which is given by [math]\displaystyle{ k_{\infty}=\omega^2/g }[/math] where [math]\displaystyle{ g }[/math] is gravity. [math]\displaystyle{ L }[/math] is a linear operator which relates the normal and potential on the body surface through the physics of the body. The simplest case is for a fixed body where the operator is [math]\displaystyle{ L=0 }[/math] but more complicated conditions are possible.

The equation is subject to some radiation conditions at infinity. We usually assume that there is an incident wave [math]\displaystyle{ \phi^{\mathrm{{In}}}\, }[/math] is a plane wave travelling in the [math]\displaystyle{ x }[/math] direction

[math]\displaystyle{ \phi^{\mathrm{{In}}}({r},z)=Ae^{{\rm i}kx}\cosh k(z+h)\, }[/math]

where [math]\displaystyle{ A }[/math] is the wave amplitude and [math]\displaystyle{ k }[/math] is the wavenumber which is the positive real solution of the Dispersion Relation for a Free Surface. We also have to apply the Sommerfeld Radiation Condition as [math]\displaystyle{ \left|\mathbf{r}\right|\rightarrow \infty }[/math].

In two-dimensions the condition is

[math]\displaystyle{ \left( \frac{\partial}{\partial|x|}-{i}k\right) (\phi-\phi^{\mathrm{{In}}})=0,\;\mathrm{{as\;}}|x|\rightarrow\infty\mathrm{.} }[/math]

where [math]\displaystyle{ \phi^{\mathrm{{In}}} }[/math] is the incident potential and [math]\displaystyle{ k }[/math] is the wave number.

In three-dimensions the condition is

[math]\displaystyle{ \sqrt{|\mathbf{r}|}\left( \frac{\partial}{\partial|\mathbf{r}|}-{i}k\right) (\phi-\phi^{\mathrm{{In}}})=0,\;\mathrm{{as\;}}|\mathbf{r}|\rightarrow\infty\mathrm{.} }[/math]