# Difference between revisions of "Wave Scattering by Submerged Thin Bodies"

Mike smith (talk | contribs) (initializing) |
Mike smith (talk | contribs) (thin bodies) |
||

Line 4: | Line 4: | ||

We consider the problem of a thin submerged structure assumed to be rigid, submerged in a semi-infinite domain of constant depth. This submerged structure is then subject to some incident wave forcing, and the resulting scattering problem is investigated. | We consider the problem of a thin submerged structure assumed to be rigid, submerged in a semi-infinite domain of constant depth. This submerged structure is then subject to some incident wave forcing, and the resulting scattering problem is investigated. | ||

+ | |||

+ | == Equations == | ||

+ | The governing equations for our problem are as follows: | ||

+ | |||

+ | <center><math> | ||

+ | \begin{align} | ||

+ | & \Delta \phi=0, &\qquad in \ \Omega\, \\ | ||

+ | & \partial_z \phi=0, &\qquad at\ z=-h\, \\ | ||

+ | & \partial_z \phi=\alpha \phi, &\qquad at\ z=0\, \\ | ||

+ | \end{align} | ||

+ | </math></center> | ||

+ | |||

+ | Plus radiation conditions: | ||

+ | <center><math> | ||

+ | \partial_x \phi \pm i k \phi = 0, \qquad as \ x \rightarrow \pm \infty \, | ||

+ | </math></center> | ||

+ | |||

+ | |||

+ | For the purposes of using Boundary Element Methods later in this article, we restrict our domain to a finite region by imposing artificial boundary conditions at <math> x= a_1 \,</math> and <math> x= a_2 \,</math>. We denote these conditions by | ||

+ | |||

+ | <center><math> | ||

+ | \begin{align} | ||

+ | & \phi=\tilde{\phi}_1(z), &\qquad at \ x=a_1\, \\ | ||

+ | & \phi=\tilde{\phi}_2(z), &\qquad at \ x=a_2\, \\ | ||

+ | \end{align} | ||

+ | </math></center> | ||

+ | |||

+ | and will be considered in more detail later on. | ||

== Thin bodies == | == Thin bodies == | ||

− | have a | + | Additionally, we have a condition for the response of our structure (which is submerged at constant finite depth <math> -d </math>): |

+ | |||

+ | |||

+ | <center><math> | ||

+ | \begin{align} | ||

+ | & \partial_n \phi = f(x,t), &\qquad x \in (-L,L), \ z=-d \, | ||

+ | \end{align} | ||

+ | </math></center> | ||

+ | |||

+ | We denote this segment of the boundary by <math> \Gamma \, </math>. | ||

− | |||

− | |||

− | == | + | There are additional considerations to be made when dealing with thin obstacles - we need to split <math> \Gamma \,</math> into two regions (<math> \Gamma^{\pm} \,</math>) and distinguish between their respective normal derivatives (<math> \partial_n \rightarrow \partial_{n^{\pm}}\,</math>). This allows us to express the boundary condition in the form |

+ | |||

+ | <center><math> | ||

+ | \begin{align} | ||

+ | & \partial_{n^{\pm}} \phi = \pm f(x,t), &\qquad x \in (-L,L), \ z=-d \, | ||

+ | \end{align} | ||

+ | </math></center> | ||

+ | |||

+ | |||

+ | Observe that the normal derivatives are related in the following manner | ||

+ | |||

+ | <center><math> | ||

+ | \begin{align} | ||

+ | & \partial_{n^{-}} = - \partial_{n^{+}} & \, | ||

+ | \end{align} | ||

+ | </math></center> | ||

+ | |||

+ | |||

+ | From this, we can express our problem involving the submerged structure entirely in terms of <math> \Gamma^{+} \,</math>, and omit the notation in the future (ie <math>\partial_{n^{+}} = \partial_n</math>). | ||

+ | |||

+ | |||

+ | For the case of a rigid structure, we define the associated boundary condition to be | ||

+ | |||

+ | |||

+ | <center><math> | ||

+ | \begin{align} | ||

+ | & \partial_n \phi = 0, &\qquad on \ \Gamma \, | ||

+ | \end{align} | ||

+ | </math></center> | ||

+ | |||

+ | == Integral Equation formulation == | ||

+ | Using Green's second identity, it can be shown that | ||

+ | |||

+ | <center><math> | ||

+ | \int_{\Gamma} (\phi \partial_{n^{_{'}}} G - G \partial_{n^{_{'}}} \phi )ds^{_{'}}=\int_{\Gamma} [ \ \phi \ ] \partial_{n^{_{'}}} G ds^{_{'}} | ||

+ | </math></center> | ||

+ | |||

+ | where <math>[ \ \phi \ ]=\phi(x^{+})-\phi(x^{-})</math>, and is typically referred to as the 'jump in <math>\phi \,</math>' | ||

== Boundary Element Methods == | == Boundary Element Methods == | ||

refer to the wikiwaves page etc | refer to the wikiwaves page etc |

## Revision as of 20:59, 28 September 2009

## Introduction

We consider the problem of a thin submerged structure assumed to be rigid, submerged in a semi-infinite domain of constant depth. This submerged structure is then subject to some incident wave forcing, and the resulting scattering problem is investigated.

## Equations

The governing equations for our problem are as follows:

Plus radiation conditions:

For the purposes of using Boundary Element Methods later in this article, we restrict our domain to a finite region by imposing artificial boundary conditions at [math]\displaystyle{ x= a_1 \, }[/math] and [math]\displaystyle{ x= a_2 \, }[/math]. We denote these conditions by

and will be considered in more detail later on.

## Thin bodies

Additionally, we have a condition for the response of our structure (which is submerged at constant finite depth [math]\displaystyle{ -d }[/math]):

We denote this segment of the boundary by [math]\displaystyle{ \Gamma \, }[/math].

There are additional considerations to be made when dealing with thin obstacles - we need to split [math]\displaystyle{ \Gamma \, }[/math] into two regions ([math]\displaystyle{ \Gamma^{\pm} \, }[/math]) and distinguish between their respective normal derivatives ([math]\displaystyle{ \partial_n \rightarrow \partial_{n^{\pm}}\, }[/math]). This allows us to express the boundary condition in the form

Observe that the normal derivatives are related in the following manner

From this, we can express our problem involving the submerged structure entirely in terms of [math]\displaystyle{ \Gamma^{+} \, }[/math], and omit the notation in the future (ie [math]\displaystyle{ \partial_{n^{+}} = \partial_n }[/math]).

For the case of a rigid structure, we define the associated boundary condition to be

## Integral Equation formulation

Using Green's second identity, it can be shown that

where [math]\displaystyle{ [ \ \phi \ ]=\phi(x^{+})-\phi(x^{-}) }[/math], and is typically referred to as the 'jump in [math]\displaystyle{ \phi \, }[/math]'

## Boundary Element Methods

refer to the wikiwaves page etc

## Eigenfunction Matching

refer to the wikiwaves page etc

### Matlab Code

## References

Article by Yang? Linton McIvor?