# Generalized Eigenfunction Expansion for Water Waves for a Fixed Body

## Contents

## Introduction

The generalized eigenfunction method goes back to the work of Povzner 1953 and Ikebe 1960. The generalized eigenfunction method has been applied to water-wave problems by Friedman and Shinbrot 1967, Hazard and Lenoir 2002 (for the case of a rigid body in water of infinite depth) and Meylan 2002 (for a thin plate on water of shallow draft). We will present here the theory for a rigid body in water of finite depth. rms of this solutions (which we call the generalised eigenfunctions) because they solve for

The generalized eigenfunction method is based on an inner product in which the evolution operator is self-adjoint. It follows from the self-adjointness that we can expand the solution in the eigenfunctions of the operator. These eigenfunctions are nothing more than the single-frequency solutions. The main difficulty is that the eigenfunctions are associated with a continuous spectrum, and this requires that they be carefully normalized. Once this is done, we can derive simple expressions which allow the time-domain problem to be solved in terms of the single-frequency solutions. The mathematical ideas are discussed in detail in Hazard and Loret 2007.

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]\Omega[/math], the boundary of the fluid domain
which touches the fixed bodies by [math]\partial\Omega[/math], and the free
surface by [math]F.[/math] The [math]x[/math] and [math]z[/math] coordinates are such that [math]x[/math] is
pointing in the horizontal direction and [math]z[/math] is pointing in the
vertical upwards direction (we denote [math]\mathbf{x}=\left( x,z\right) ).[/math] The
free surface is at [math]z=0[/math] and the sea floor is at [math]z=-h[/math]. The
fluid motion is described by a velocity potential [math]\Phi[/math] and free surface by
[math]\zeta[/math].

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

At the bottom surface we have no flow

At the free surface we have the kinematic condition

and the dynamic condition (the linearized Bernoulli equation)

The body boundary condition for a fixed body is

The initial conditions are

## Equations in the Frequency Domain

We also assume that Frequency Domain Problem with frequency [math]\omega[/math] and we assume that all variables are proportional to [math]\exp(-\mathrm{i}\omega t)\,[/math]

The water motion is represented by a velocity potential which is denoted by [math]\phi\,[/math] so that

[math] \Phi(\mathbf{x},t) = \mathrm{Re} \left\{\phi(\mathbf{x},\omega)e^{-\mathrm{i} \omega t}\right\}. [/math]

The equations are the following

(note that the last expression can be obtained from combining the expressions:

where [math]\alpha = \omega^2/g \,[/math]) The body boundary condition for a rigid body is just

The equation is subject to some radiation conditions at infinity. We assume the following.
[math]\phi^{\mathrm{I}}\,[/math]
is a plane wave travelling in the [math]x[/math] direction,

where [math]A [/math] is the wave amplitude (in potential) [math]\mathrm{i} k [/math] is the positive imaginary solution of the Dispersion Relation for a Free Surface (note we are assuming that the time dependence is of the form [math]\exp(-\mathrm{i}\omega t) [/math]) and

In two-dimensions the Sommerfeld Radiation Condition is

where [math]\phi^{\mathrm{{I}}}[/math] is the incident potential.

## Time domain calculations

The solution in the frequency domain can be used to construct the solution in the time domain. This is well-known for the case of a plane incident wave which is initially far from the body, and in this case the solution can be calculated straightforwardly using the standard Fourier transform. However, when we consider an initial displacement which is non-zero around the bodies, we cannot express the time-dependent solution in terms of the single frequency solutions by a standard Fourier transform.

We begin with the equations in the time domain subject to the initial conditions given by Denoting the potential at the surface by

we introduce the operator [math]\mathbf{G}[/math] which maps the surface potential to the potential throughout the fluid domain. The operator [math]\mathbf{G}\psi[/math] is found by solving

and is defined by [math]\mathbf{G}\psi=\Psi.[/math] The operator [math]\partial_{n}\mathbf{G}[/math], which maps the surface potential to the normal derivative of potential at the surface (called the Dirichlet-to-Neumann map) is given by

Therefore the equations in the time domain can be written as

where we can recover the potential using the operator [math]\mathbf{G}[/math]. The evolution operator [math]\partial_{n}\mathbf{G}[/math] is symmetric in the Hilbert space given by the following inner product

where [math] *[/math] denotes complex conjugate, and we assume that this symmetry implies that the operator is self-adjoint. We can prove the symmetry by using Green's second identity

where [math]\phi[/math] and [math]\psi[/math] are the surface potentials associated with [math]\zeta[/math] and [math]\eta[/math] respectively.

## Eigenfunctions of [math]\partial_{n}\mathbf{G}[/math]

The eigenfunctions of [math] \partial_{n}\mathbf{G}[/math] satisfy

To solve for the eigenfunctions of [math]\partial_{n}\mathbf{G}[/math] we need to solve the frequency-domain equations, and the radian frequency [math]\omega[/math] is exactly the eigenvalue. To actually calculate the eigenfunctions of [math]\partial_{n}\mathbf{G}[/math] we need to specify the incident wave potential, and for each frequency we have two eigenfunctions (waves incident from the left and from the right). It is possible for there to exist point spectra for this operator which correspond to the existence of a trapped mode McIver 1996 and the presence of a trapped mode requires that the generalized eigenfunction expansion we derive must be modified.

## Normalization of the Eigenfunctions

The eigenfunctions of [math]\partial_{n}\mathbf{G}[/math] (with eigenvalue [math]\omega[/math]) are denoted by [math] \zeta_{\kappa}(x,k\left( \omega\right) ) [/math], where [math]\kappa = \pm 1[/math] for waves incident from the left or right respectively. We must also specify radiations conditions to normalize the eigenfunctions and they are given by

Note that [math]R_{\kappa}[/math] and [math]T_{\kappa}[/math] are the reflection and transmission coefficients respectively and that we have normalized so that the amplitude (in displacement) is unity. Also [math]k = - \mathrm{i} k_0[/math].

and we will consider both [math]k(\omega)[/math] and [math]\omega(k)[/math]
in what follows as required. The solution of the single-frequency
equation may be computationally challenging and for the generalized
eigenfunction expansion the major numerical work is to determine the
single-frequency solutions.
As mentioned previously, determining [math]\zeta_\kappa[/math] is the major
computation of the generalized eigenfunction method, but we simply
assume that they are known. We know that the eigenfunctions are
orthogonal for different [math]\omega[/math] (from the self-adjointness of
[math]\partial_n\mathbf{G}[/math]), and that the waves incident from the left and right
with the same [math]\omega[/math] are orthogonal from the identity

Mei 1983. It therefore follows that

but we need to determine the normalizing function [math]\Lambda_{n}\left( \omega_{n}\right)[/math]. This is achieved by using the result that the eigenfunctions satisfy the same normalizing condition with and without the scatterers present. This result, the proof of which is quite technical, is well-known and has been shown for many different situations. The original proof was for Schr\"odinger's equation and was due to Povzner 1953 and Ikebe 1960. A proof for the case of Helmholtz equation was given by Wilcox 1975. Recently the proof was given for water waves by Hazard and Lenoir 2002 and Hazard and Loret 2007.

Since the eigenfunctions satisfy the same normalizing condition with and without the scatterers, we normalize with the scatterers absent. This means that the eigenfunctions are simply the incident waves, and the free surface [math]F[/math] is the entire axis. This allows us to derive

This result allows us to calculate the time-dependent solution in the eigenfunctions (or single-frequency solutions).

## Expansion in Eigenfunctions

We expand the solution for the displacement in the time domain as

where [math]f_\kappa[/math] and [math]g_\kappa[/math] will be determined from the initial conditions. Note that here, and in subsequent equations, we are assuming that [math]k[/math] is a function of [math]\omega[/math] or that [math]\omega[/math] is a function of [math]k[/math] as required. If we take the inner product with respect to the eigenfunctions [math]\zeta_\kappa[/math] we obtain

and

We can therefore write, changing the variable of integration to [math]k[/math], as

If we take the case when [math]\partial_t\zeta_0( x) =0[/math] and write the integral given by the inner product explicitly, we obtain

## Case of a trapped mode

If the scattering structure supports a trapped mode at a particular frequency [math]\omega_{0}[/math] then the expression for the free-surface elevation becomes

where [math]\tilde\zeta(x)[/math] is the trapped mode free-surface elevation. This formula can be extend to the case when there is more than one trapped mode.

## An identity linking waves from the left and right

A consequence of the requirement that the displacement be real, if the initial displacement and initial derivative of displacement is real, is that

must be purely real. This can only be true if

## Matlab Code

- A program to solve two small docks gem_two_small_dock.m

with output shown on the right

### Additional code

This program requires

- Calculate Fourier integral k_integrate.m
- A program to calculate the coefficients for the finite dock problems can be found here two_finite_docks_symmetry.m
- dispersion_free_surface.m
- gif_add_frame.m