Linear Potential Flow Theory

QBlade is capable of calculating hydrodynamic forces on submerged elements due to floater motion and waves by making use of potential flow theory. This allows the full interaction between the three-dimensional floater geometry and the wave field to be accounted for.

Potential flow and boundary conditions

Here it is assumed that the flow field is irrotational. This allows the flow field to be described in terms of a scalar potential $$\phi$$, which satisfies Laplace’s differential equation $$\nabla^2\phi=0$$. The velocity field $$\vec{v}$$ can be expressed as the gradient of this potential:

$$$\vec{v} = \nabla \phi = \frac{\partial \phi}{\partial x}\vec{e}_x + \frac{\partial \phi}{\partial y}\vec{e}_y + \frac{\partial \phi}{\partial z}\vec{e}_z$$ \textrm{ .}$

Three boundary conditions are require to specify fully the problem 1. The first enforces continuity at the free surface. This is linearised (hence linear potential flow) to give:

$$$\frac{\partial \phi}{\partial z} - \frac{\omega^2}{g}\phi = 0 \hspace{5mm} \textrm{on} \hspace{5mm} z=0$$ \textrm{ ,}$

where $$g$$ is the acceleration due to gravity and $$\omega$$ is the discrete frequency being analysed. The second enforces the boundary condition on the sea bottom:

\begin{split}\begin{align} \nabla\phi \rightarrow 0 \hspace{2mm} \textrm{as} \hspace{2mm} z \rightarrow -\infty & \hspace{10mm} \textrm{Infinite depth} \\ \frac{\partial \phi}{\partial z} = 0 \hspace{2mm} \textrm{on} \hspace{2mm} z = -h & \hspace{10mm} \textrm{Finite depth } h \textrm{ .} \\ \end{align}\end{split}

Finally, the Sommerfeld condition states that wave energy associated with disturbance due to the body is radiated in all directions. A few assumptions are outlined here for the following discussion:

• The floater undergoes negligibly small motions away from the equilibrium position.

• Solutions are assumed to be harmonic $$\hat{\phi} = \phi\,e^{i\omega t}$$.

• The floater is submerged in a fluid with density $$\rho$$ $$\textrm{kg}\,\textrm{m}^{-3}$$.

Equations of Motion

The motion of a generic floating body in dimension $$j$$: $$x_j(t)$$ can be modelled with the equation due to Cummins2:

$$$(M_{ij}+A_{ij}^{\infty})\ddot{x}_j(t) + \int_{-\infty}^{t}K_{ij}(t-\tau)\dot{x}_j(\tau)\,d\tau + C_{ij}x_j(t) = F_j^{w}(t) - F_j^{e}(x,\dot{x},t) \textrm{ .}$$$

In the above expression the indices $$i$$ and $$j$$ represent the degree of freedom of motion and the acting force, respectively. In this equation:

• $$M_{ij}$$ represents the inertia of the floater

• $$A_{ij}^{\infty}$$ is the added mass matrix (see Radiation Forces below)

• $$K_{ij}$$ is the radiation damping matrix (see Radiation Forces below)

• $$C_{ij}$$ is the hydrostatic stiffness matrix (see Hydrostatic Forces)

• $$F_j^{w}$$ are the forces due to waves (see Excitation forces below)

• $$F_j^{m}$$ are external forces due to moorings

These terms shall be described in the following sections.

The motion of the floater in an undisturbed field generates waves which radiate away from the body. Newton’s third law dictates that the force required to set this disturbance in motion gives rise to an equal and opposite force on the body. This acts in the form of a pressure disturbance. An example of this is shown in the following figure for a motion in the surge direction.

The disturbance potential $$\phi_r$$ due to this motion is sought. The body is taken as being impermeable which implies that the local velocity (gradient of the potential) must be equivalent to the local body motion:

$$$\frac{\partial \phi_r}{\partial n} = \nabla\phi\cdot\vec{n} = \vec{u}$$ \textrm{ ,}$

where $$\vec{n}$$ is the normal vector at the surface. The solution for $$\phi_r$$ on the geometry represents the hydrodynamic pressure acting. This is integrated over the floater surface $$S_b$$ to give the total forcing due to motion (note that this is complex as a harmonic solution has been assumed):

$$$A_{ij} - \frac{i}{\omega}B_{ij} = \rho \iint_{S_b}n_i\,\phi_{r,j}\,dS$$ \textrm{ ,}$

where $$A_{ij}$$ and $$B_{ij}$$ are referred to as the added mass and radiation damping matrices. It is important to note that the terms above are calculated in a frequency domain analysis. The added mass matrix can be taken directly from the frequency domain analysis as:

$$$A_{ij}^{\infty} = \lim_{\omega\to\infty} A_{ij}(\omega)$$ \textrm{ ,}$

where it is important that a sufficiently high analysis frequency $$\omega$$ is taken to ensure that $$A_{ij}(\omega)$$ has converged. An impulsive motion in any direction generates a force which is time-varying, this is accounted for with the time convolution in the equations of motion above. The time convolution kernel is referred to as the impulse response function, or IRF and is calculated as:

$$$K_{ij}(t) = \frac{2}{\pi}\int_{0}^{\infty} \omega A_{ij}(\omega)\sin \omega t\,d\omega = \frac{2}{\pi}\int_{0}^{\infty} B_{ij}(\omega)\cos \omega t\,d\omega$$ \textrm{ .}$

In practise the second form is used due to its easier numerical integration.

The arrays for $$A_{ij}(\omega)$$ and $$B_{ij}(\omega)$$ can be imported into QBlade in NEMOH, WAMIT, or BEMUse formats. This integration is carried out numerically with a frequency step size $$\Delta_{\omega}$$.

$$$K_{ij}(t) = \frac{2}{\pi}\int_{0}^{\infty} B_{ij}(\omega)\cos \omega t\,d\omega \approx \frac{2}{\pi} \sum_{n=1}^{n=\omega_{max}} \Delta_{\omega} B_{ij}(n\Delta_{\omega})\cos tn\Delta_{\omega}$$ \textrm{ .}$

The decay of $$K_{ij}$$ implies that the time convolution can be truncated to a finite time $$T$$. The radiation force $$F^{r}$$ can be calculated by carrying out a time convolution numerically with timestep $$\Delta_t$$:

$$$F_j^{r}(t) = \int_{t-T}^{t}K_{ij}(t-\tau)\dot{x}_j(\tau)\,d\tau \approx \sum_{i=1}^{i=T/\Delta_t} \Delta_t K_{ij}(i\Delta_t)\dot{x}_j(t-i\Delta_t)$$ \textrm{ .}$

Excitation forces

The boundary condition on the surface of the floater causes incoming waves to be reflected away. As with the radiation forces, this gives rise to a disturbance potential $$\phi_d$$ and a corresponding force $$X_j$$ which acts on the floater. The Haskind relations allows these forces to be expressed in terms of the radiation potential:

$$$X_j = -i\omega\rho\iint_{S_b} \left( n_i\phi_0- \phi_{r,j}\frac{\partial \phi_0}{\partial n}\right)dS$$ \textrm{ ,}$

where $$\phi_0$$ is the potential of the incoming wave. An IRF for this is calculated as:

$$$H_{ij}(t) = \frac{1}{2\pi}\int_{-\infty}^{\infty}X_j(\omega)e^{i\omega t}\,d\omega$$ \textrm{ .}$

This is numerically integrated with a frequency step size $$\Delta_{\omega}$$. The lower limit of the integral indicates that the excitation IRF is non-causal, which means the chosen input is not the cause for the output. In this context, the non-causality may be explained by the fact that the incident wave hits the body and exerts a wave force before the wave reaches the chosen reference point for the body (usually located at the geometrical center), see Falnes3:.

As with the radiation forces, the time-domain excitation forces $$F^{e}$$ are calculated with a time convolution with the IRF given above:

$$$F_j^{e}(t) = \int_{-\infty}^{\infty} H_{ij}(\tau)\zeta(x_0,y_0,t-\tau)\,d\tau$$ \textrm{ .}$

where $$\zeta(x_0,y_0,t)$$ is the wave elevation at the reference position $$(x_0,y_0)$$ during time $$t$$.

Since the non-causality of the excitation IRF means $$H_{ij}(t) ≠ 0$$ for $$(t) < 0$$, future wave information is required before the waves reach the neutral reference position of the floater, see Falnes3. This integral is again calculated numerically over a truncated time period $$T$$ with the time step $$dT$$.

1

J.N. Newman. Marine Hydrodynamics. MIT Press, 2018. ISBN 9780262534826.

2

W.E. Cummins. The impulse response function and ship motions. Shiffstechnik, pages 101–109, 1961.

3(1,2)

J. Falnes. On non-causal impulse response functions related to propagating water waves. Applied Ocean Research, 17(6):379–389, 1995. doi:https://doi.org/10.1016/S0141-1187(96)00007-7.