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

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

(12 intermediate revisions by the same user not shown) | |||

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. | ||

+ | |||

+ | <!--[[Image:Submerged_str.png|600px|right||frame|caption here]] --> | ||

+ | |||

+ | == 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, in tandem with the discussion of thin bodies (above), we can obtain an expression for <math>\phi \,</math> around the outer boundary (<math>\partial \Omega \,</math>) of our domain: | ||

+ | |||

+ | |||

+ | <center><math> | ||

+ | \frac{1}{2}\phi(\textbf{x})=\int_{\Gamma} [ \ \phi \ ] \partial_{n^{_{'}}} G ds^{_{'}} +\int_{\partial \Omega} (\phi \partial_{n^{_{'}}} G - G \partial_{n^{_{'}}} \phi )ds^{_{'}} | ||

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

+ | |||

+ | |||

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

+ | |||

+ | Note carefully that <math>\phi \,</math> relates to points on the outer boundary, <math>\partial \Omega \,</math>, whereas <math>[ \ \phi \ ] \,</math> relates to points on the inner boundary, <math>\Gamma \,</math>. | ||

+ | |||

+ | |||

+ | |||

+ | In order to compute <math>[ \ \phi \ ] \,</math>, we will restrict all '''x''' points in the expression above to be on <math>\Gamma \,</math>, and then take the normal derivative of the expression to obtain | ||

+ | |||

+ | |||

+ | |||

+ | <center><math> | ||

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

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

+ | |||

+ | Note that in our case, <math>\oint \,</math> refers to a Hadamard finite-part integral, as opposed to a contour integral. | ||

+ | |||

+ | ===Hadamard Finite Part Integrals === | ||

+ | A hadamard finite part integral is an example of a hypersingular integral equation (an integral equation whose kernel has high order singularities), and must be treated with caution. | ||

+ | |||

+ | In order to evaluate our integral above, it is necessary to express the jump in <math>\phi \,</math> as | ||

+ | |||

+ | <center><math> | ||

+ | [ \ \phi(x^{_{'}}) \ ] = (1-x^{_{'}2})^{1/2} \sum_{m=0}^{M}b_m U_m(x^{_{'}}) | ||

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

+ | where <math>U_m(x^{_{'}}) \,</math> denotes Chebyshev Polynomials of the Second Kind. | ||

+ | |||

+ | It is also necessary to employ the powerful result discussed in [[Martin and Rizzo 1989]]: | ||

+ | |||

+ | <center><math> | ||

+ | \oint_{-1}^{1} \frac{(1-v^2)^{1/2}U_m(v)}{(u-v)^2}dv = -\pi (m+1) U_m(u) | ||

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

+ | Note that the structure has to be non-dimensionalised accordingly. | ||

+ | |||

+ | |||

+ | From this we can obtain | ||

+ | |||

+ | <center><math> | ||

+ | -\frac{1}{2}\sum_{m=0}^{M} b_m (m+1)U_m(x) + \int_{\partial \Omega} \partial_{n^{_{'}}n} G \phi ds^{_{'}}- \int_{\partial \Omega}\partial_n G \partial_{n^{_{'}}} \phi ds^{_{'}} = 0 | ||

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

+ | |||

+ | |||

+ | Multiplying this entire expression by <math>(1-x^2)^{1/2} U_n(x) \,</math> and integrating allows us to obtain the matrix representation | ||

+ | |||

+ | |||

+ | <center><math> | ||

+ | \tilde{M} \vec{b} + \mathbb{Y} (\mathbb{G}^{(1)}\vec{\phi} - \mathbb{G}^{(2)} \vec{\phi}_{n^{_{'}}} ) = 0 | ||

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

+ | |||

+ | where | ||

+ | |||

+ | <center><math> | ||

+ | \begin{align} | ||

+ | \tilde{M}_{nn} &= -\frac{\pi}{4}n \\ | ||

+ | \mathbb{Y} &= \mathbb{U} \ \mathbb{W}^{(1)}\mathbb{W}^{(2)} \\ | ||

+ | \mathbb{U} &= [U_1(x),U_2(x), ... , U_M(x)]\\ | ||

+ | \mathbb{W}^{(1)}_{nn} &= (1-x_n^2)^{1/2} | ||

+ | \end{align} | ||

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

− | |||

− | |||

− | + | Note that <math>\mathbb{W}^{(2)}</math> is a weights matrix for numerically integrating over <math>\Gamma \,</math>, and the two <math>\mathbb{G}</math> matrices are obtained using boundary element method, which is discussed below. | |

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

− | + | For the purposes of this article, the webpages [[:Category:Boundary Element Method|Boundary Element Method]] and [[Boundary Element Method for a Fixed Body in Finite Depth]] provide an excellent explanation to the Boundary Element method techniques used here. | |

+ | |||

+ | Recall the expression | ||

+ | |||

+ | <center><math> | ||

+ | \frac{1}{2}\phi(\textbf{x})=\int_{\Gamma} [ \ \phi \ ] \partial_{n^{_{'}}} G ds^{_{'}} +\int_{\partial \Omega} (\phi \partial_{n^{_{'}}} G - G \partial_{n^{_{'}}} \phi )ds^{_{'}} | ||

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

+ | |||

+ | Using the Chebyshev series from above, we can represent this as | ||

+ | |||

+ | <center><math> | ||

+ | \mathbb{G}^{(3)} \mathbb{W}^{(1)} \mathbb{U} \vec{b}+\mathbb{G}^{(4)}\vec{\phi} - \mathbb{G}^{(5)} \vec{\phi}_{n^{_{'}}} = \frac{1}{2} \vec{\phi} | ||

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

+ | |||

+ | |||

+ | |||

+ | This leaves us with the following two equations | ||

+ | <center><math> | ||

+ | \begin{align} | ||

+ | \tilde{M} \vec{b} + \mathbb{Y} (\mathbb{G}^{(1)}\vec{\phi} - \mathbb{G}^{(2)} \vec{\phi}_{n^{_{'}}} ) &= 0 \\ | ||

+ | \mathbb{G}^{(3)} \mathbb{W}^{(1)} \mathbb{U} \vec{b}+\mathbb{G}^{(4)}\vec{\phi} - \mathbb{G}^{(5)} \vec{\phi}_{n^{_{'}}} &= \frac{1}{2} \vec{\phi} \\ | ||

+ | \end{align} | ||

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

+ | |||

+ | |||

+ | Using the technique outlined in [[Boundary Element Method for a Fixed Body in Finite Depth]] we express | ||

+ | <center><math> | ||

+ | \vec{\phi}_{n^{_{'}}} = A \vec{\phi} -\vec{f} | ||

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

+ | |||

+ | Consequently we can represent our system in the following form | ||

+ | |||

+ | <center><math> | ||

+ | \begin{bmatrix} | ||

+ | \tilde{M} & \mathbb{Y}(\mathbb{G}^{(1)} - \mathbb{G}^{(2)}A) \\ | ||

+ | \mathbb{G}^{(3)}\mathbb{W}^{(1)}\mathbb{U} & (\mathbb{G}^{(4)}- \mathbb{G}^{(5)}A - \frac{1}{2}\mathrm{I}) | ||

+ | \end{bmatrix} | ||

+ | |||

+ | \begin{bmatrix} | ||

+ | \vec{b} \\ | ||

+ | \vec{\phi} | ||

+ | \end{bmatrix} | ||

+ | |||

+ | = | ||

+ | |||

+ | \begin{bmatrix} | ||

+ | -\mathbb{Y}\mathbb{G}^{(2)}\vec{f} \\ | ||

+ | -\mathbb{G}^{(5)}\vec{f} | ||

+ | \end{bmatrix} | ||

+ | |||

+ | \,\!</math></center> | ||

+ | |||

+ | From this we are able to compute the jump in <math>\phi \,</math>, and the potential around the outer boundary <math>\partial \Omega</math> | ||

== Eigenfunction Matching == | == Eigenfunction Matching == | ||

− | + | As an alternative to the the technique outlined on this webpage, [[Eigenfunction Matching for a Submerged Finite Dock]] can be used to determine the potential at any point in the region, as well as the jump in <math>\phi \,</math> | |

===Matlab Code === | ===Matlab Code === | ||

+ | Will be posted presently | ||

==References == | ==References == | ||

− | + | [[Martin and Rizzo 1989]] | |

− | Linton | + | |

+ | [[Linton and McIver 2001]] | ||

+ | |||

+ | ==Bibtex== | ||

+ | The bibtex reference for this webpage is below: | ||

+ | |||

+ | |||

+ | |||

+ | <pre> | ||

+ | @MISC{referenceID,</pre> | ||

+ | author = {wikiwaves}, | ||

+ | title = { {{PAGENAME}} }, | ||

+ | month = {{CURRENTMONTHNAME}}, | ||

+ | year = { {{CURRENTYEAR}} }, | ||

+ | url = { {{fullurl:{{PAGENAME}}}} } | ||

+ | <pre> | ||

+ | } | ||

+ | </pre> |

## Latest revision as of 22:54, 29 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, in tandem with the discussion of thin bodies (above), we can obtain an expression for [math]\displaystyle{ \phi \, }[/math] around the outer boundary ([math]\displaystyle{ \partial \Omega \, }[/math]) of our domain:

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

Note carefully that [math]\displaystyle{ \phi \, }[/math] relates to points on the outer boundary, [math]\displaystyle{ \partial \Omega \, }[/math], whereas [math]\displaystyle{ [ \ \phi \ ] \, }[/math] relates to points on the inner boundary, [math]\displaystyle{ \Gamma \, }[/math].

In order to compute [math]\displaystyle{ [ \ \phi \ ] \, }[/math], we will restrict all **x** points in the expression above to be on [math]\displaystyle{ \Gamma \, }[/math], and then take the normal derivative of the expression to obtain

Note that in our case, [math]\displaystyle{ \oint \, }[/math] refers to a Hadamard finite-part integral, as opposed to a contour integral.

### Hadamard Finite Part Integrals

A hadamard finite part integral is an example of a hypersingular integral equation (an integral equation whose kernel has high order singularities), and must be treated with caution.

In order to evaluate our integral above, it is necessary to express the jump in [math]\displaystyle{ \phi \, }[/math] as

where [math]\displaystyle{ U_m(x^{_{'}}) \, }[/math] denotes Chebyshev Polynomials of the Second Kind.

It is also necessary to employ the powerful result discussed in Martin and Rizzo 1989:

Note that the structure has to be non-dimensionalised accordingly.

From this we can obtain

Multiplying this entire expression by [math]\displaystyle{ (1-x^2)^{1/2} U_n(x) \, }[/math] and integrating allows us to obtain the matrix representation

where

Note that [math]\displaystyle{ \mathbb{W}^{(2)} }[/math] is a weights matrix for numerically integrating over [math]\displaystyle{ \Gamma \, }[/math], and the two [math]\displaystyle{ \mathbb{G} }[/math] matrices are obtained using boundary element method, which is discussed below.

## Boundary Element Methods

For the purposes of this article, the webpages Boundary Element Method and Boundary Element Method for a Fixed Body in Finite Depth provide an excellent explanation to the Boundary Element method techniques used here.

Recall the expression

Using the Chebyshev series from above, we can represent this as

This leaves us with the following two equations

Using the technique outlined in Boundary Element Method for a Fixed Body in Finite Depth we express

Consequently we can represent our system in the following form

From this we are able to compute the jump in [math]\displaystyle{ \phi \, }[/math], and the potential around the outer boundary [math]\displaystyle{ \partial \Omega }[/math]

## Eigenfunction Matching

As an alternative to the the technique outlined on this webpage, Eigenfunction Matching for a Submerged Finite Dock can be used to determine the potential at any point in the region, as well as the jump in [math]\displaystyle{ \phi \, }[/math]

### Matlab Code

Will be posted presently

## References

## Bibtex

The bibtex reference for this webpage is below:

@MISC{referenceID,

author = {wikiwaves}, title = { Wave Scattering by Submerged Thin Bodies }, month = February, year = { 2023 }, url = { https://wikiwaves.org/Wave_Scattering_by_Submerged_Thin_Bodies }

}