A time-domain second-order FEM model for the wave diffraction-radiation problem. Validation with a semisubmersible platform

<!-- metadata commented in wiki content


<big>'''A time-domain second-order FEM model for the wave diffraction-radiation problem. Validation with a semisubmersible platform.'''</big>

Borja Servan-Camas<sup>(a,*)</sup>, Jose E. Gutierrez-Romero<sup>(b)</sup>, Julio Garcia-Espinosa<sup>(a,c)</sup>

(a) Centre Internacional de Metodes Numerics en Enginyeria (CIMNE),[http://www.cimne.upc.edu/ http://www.cimne.upc.edu/]

Edifici C1, Campus Norte, UPC, Gran Capitán s/n, 08034 Barcelona, Spain

(b) Depart. Naval Technology, Universidad Politécnica de Cartagena (UPCT), [http://www.upct.es/ http://www.upct.es/],

Paseo Alfonso XIII 48, 30203 Cartagena, Spain

(c) Universidad Politécnica de Cataluña (UPC), C. Gran Capitan s/n, Campus Nord, 08034 Barcelona, Spain.

(*) Corresponding autor.
-->

===Abstract===

A finite element method for the solution of the up-to-second-order wave diffraction-radiation problem in the time-domain is proposed. The solver has been verified against available analytical solutions, and validated using experimental data available for the HiPRWind semisubmersible platform (designed for floating wind turbines). To validate, the wave diffraction-radiation solver is coupled to body dynamics and mooring solvers in the time-domain. The HiPRWind movements and mooring forces have been compared for a large number of test cases, including decay tests, monochromatic waves, and bichromatic waves obtaining good correlation for both body movements and mooring forces.

===1. Introduction===

There is a growing interest in the industry of Floating Offshore Wind Turbines (FOWT) due to their ability to access the vast quantities of wind resources available over deep waters. Despite the existence of full scale prototypes which are already operating, such as the Hywind in Norway [<span id='cite-1'></span>[[#1|1]]] or the Windfloat in Portugal [<span id='cite-2'></span>[[#2|2]]], the industry still faces design and operational challenges which require the development of new modeling tools. In this work, we propose a model to analyse the up-to-second-order response of floating structures, which have been validated through experimental tests conducted for the HiPRWind semisubmersible FOWT model. A review on floating offshore wind technology can be found in [<span id='cite-3'></span>[[#3|3]],<span id='cite-4'></span>[[#4|4]],<span id='cite-5'></span>[[#5|5]]].

<span id='_Ref467220412'></span><span id='_Ref467219995'></span>The hydrodynamics of the semisubmersible concept for FOWTs has received attention in recent literature. For instance, [<span id='cite-6'></span>[[#6|6]]] and [<span id='cite-7'></span>[[#7|7]]] focused on slow-drift and mean-drift forces of semisubmersible platforms. A comparison of a semisubmersible against the SPAR concept can be found in [<span id='cite-8'></span>[[#8|8]],<span id='cite-9'></span>[[#9|9]],<span id='cite-10'></span>[[#10|10]]], and simulations in the Time-Domain (TD) considering different models for the hydrodynamic loads were carried out in [<span id='cite-11'></span>[[#11|11]]].

One of the main concerns regarding semisubmersible platforms are the slow-drift forces. For semisubmersible platforms with catenary mooring lines these forces are usually within the range of the surge natural frequency, leading to large displacements. This is because even though the wave frequencies are usually larger than the natural frequencies, second-order effects contain low frequency components that might cause the slow-drift.

Second-order forces can increase the surge response of semisubmersible platforms, even becoming larger than the first-order response. To estimate these slow drift forces, Newman’s approximation can be used which only depends on the first order solution. However it has been shown in [<span id='cite-_Ref467219995'></span>[[#_Ref467219995|7]]] that this does not provide the necessary precision for this type of problem. Therefore, the second-order effect must be taken into account to accurately compute slow-drift motions, and the mooring system should be designed accordingly.

The impact of slow-drift forces on the design of the mooring systems and difficulties in estimating the corresponding forces, is a problem that requires substantial further research.  The design of the mooring system is based on loads induced by the excursions. These can be predicted on an extensive set of time-domain simulations comprising different environmental conditions (combinations of waves, current, and wind conditions) which aim to represent the metocean statistics of the specific location. For these simulations, a wave diffraction-radiation solver in the frequency-domain is usually used to evaluate the added masses, potential damping, and wave excitation forces.

Quadratic transfer functions (QTFs) obtained by frequency domain solvers are usually used as inputs for second-order time-domain simulations. These QTFs depend on the linear first-order response amplitude operators (RAOs), and when the mooring lines are not linear (like the catenary lines used for semisubmersible platforms), these must be linearised. Therefore, QTFs obtained in the frequency-domain do not contain information about the nonlinearities contained in the mooring system. To overcome this shortcoming, a time-domain solver could be used. While the former are computationally more efficient for linear problems, the latter accounts for these nonlinearities in a natural way.

Apart from the importance of second order effects on the dynamic response of floating structures, recent works have emphasized the importance of the second-order effects on the surge response of FOWT. Coulling et al. [<span id='cite-12'></span>[[#12|12]]] concluded that particular attention must be paid to the motion and load responses of the platforms associated with the second-order difference-frequency forces of environmental wave loads, as the exclusion of the second-order dynamic analysis leads to a reduction of the platform mean excursions. Other works [<span id='cite-13'></span>[[#13|13]],<span id='cite-14'></span>[[#14|14]]] have recently assessed the effects of second-order hydrodynamics on semisubmersible FOWT which are usually neglected in the dynamic behaviour. These forces lead to large oscillations that strain the mooring system or to vibrations that cause fatigue damage to the moored structure. Luptom and Langley state in [<span id='cite-_Ref467220412'></span>[[#_Ref467220412|6]]] that slow-drift forces might be less important for FOWT than for larger and better known floating structures such as those used in the oil industry. An accurate estimation of these forces is mandatory to assess its impact.

<span id='_Ref482260439'></span>More recently, Lopez-Pavon et al. [<span id='cite-_Ref467219995'></span>[[#_Ref467219995|7]]] and Simos et al. [<span id='cite-15'></span>[[#15|15]]] focused on the estimation and verification of second-order wave induced forces on the HiPRWind semisubmersible platform. This work concludes with the following observations:

:1. Accurate calculation of second-order forces may be difficult to achieve and in some cases, different numerical codes (based on different approximations for these forces) may render divergent results.

:2. The experimental verification of the slow-drift effects is quite difficult. Accurate measurement of low-frequency forces is hard to obtain and indirect verifications based on the resonant motions of the floating body depend on other factors such as the viscous damping of the small-scale model, the geometric characteristics of the mooring system, etc.

In this work, a Finite Element Method (FEM) that solves the up-to-second-order wave diffraction-radiation problem in the time-domain is proposed. In this model, non-linear forces such as those arising from the mooring lines can be introduced into the dynamics of the floater with no need of linearization. First, mathematical and numerical models for the wave diffraction-radiation are presented. Then the model is verified by comparing the results with available second-order analytical solutions. A model of the HiPRWind platform is calibrated using decay tests and analysed in monochromatic, bichromatic, and irregular waves in order to validate the proposed numerical approach. Finally conclusions are drawn regarding the model presented and the results obtained.

===2. Up-to-second-order diffraction-radiation governing equations===

The potential flow governing equations for the up-to-second-order wave problem are obtained by applying Taylor expansion on the boundary surfaces of a time-independent domain. This approach enables an approximation of the free surface on <math display="inline">z=</math><math>\xi</math>  and the mean wetted surface <math display="inline">{S}_{B}</math> of a floating body at time t. Then, a perturbed solution based on the Stokes expansion procedure is applied to the velocity potential, the free surface elevation, and the floater motion. Retaining terms up to second order, the resulting equations are:

{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>\Delta {\phi }^{1+2}=0</math>
|}
|  colspan='2'  style="text-align: center;white-space: nowrap;"|<math display="inline">in\, \Omega</math> ,
|  style="text-align: right;text-align: right;white-space: nowrap;"|(1)
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>\frac{\partial {\eta }^{1+2}}{\partial t}-\frac{\partial {\phi }^{1+2}}{\partial z}=</math><math>-{S}^{1}</math>
|}
|  colspan='2'  style="text-align: center;white-space: nowrap;"|<math display="inline">in\, z=</math><math>0</math>,
|  style="text-align: right;text-align: right;white-space: nowrap;"|<span id='_Ref480536160'></span>(2)
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>\frac{\partial {\phi }^{1+2}}{\partial t}+\frac{{P}_{fs}}{\rho }+g{\eta }^{1+2}=</math><math>-{R}^{1}</math>
|}
|  colspan='2'  style="text-align: center;white-space: nowrap;"|<math display="inline">in\, z=</math><math>0</math>,
|  style="text-align: right;text-align: right;white-space: nowrap;"|<span id='_Ref480536163'></span>(3)
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{v}_{\phi }^{1+2}\cdot {n}_{p}^{0}+{v}_{\phi }^{1}\cdot {n}_{p}^{1}=-\left( {v}_{p}^{1}{+v}_{\psi }^{1}\right) \cdot {n}_{p}^{1}</math>

<math display="inline">-\left( {v}_{p}^{1+2}{+v}_{\psi }^{1+2}+{r}_{p}^{1}\cdot \left( \nabla {v}_{\phi }^{1}+\right. \right. </math><math>\left. \left. \nabla {v}_{\psi }^{1}\right) \right) \cdot {n}_{p}^{0}</math>
|}
|  colspan='2'  style="text-align: center;white-space: nowrap;"|<math display="inline">on\, P\in {S}_{B}^{0}</math>,
|  style="text-align: right;text-align: right;white-space: nowrap;"|(4)
|-
|  colspan='2'  style="text-align: right;vertical-align: top;text-align: center;white-space: nowrap;"|<math>{R}^{1}={\eta }^{1}\frac{\partial }{\partial z}\left( \frac{\partial {\phi }^{1}}{\partial t}\right) +</math><math>{\zeta }^{1}\frac{\partial }{\partial z}\left( \frac{\partial {\phi }^{1}}{\partial t}\right) +</math><math>{\eta }^{1}\frac{\partial }{\partial z}\left( \frac{\partial {\psi }^{1}}{\partial t}\right) +</math><math>\frac{1}{2}\nabla {\phi }^{1}\cdot \nabla {\phi }^{1}+\nabla {\psi }^{1}\cdot \nabla {\phi }^{1},</math>
|  colspan='2'  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref480536394'></span>(5)
|-
|  colspan='2'  style="text-align: right;vertical-align: top;text-align: center;white-space: nowrap;"|<math>{S}^{1}=\frac{\partial {\phi }^{1}}{\partial x}\frac{\partial {\eta }^{1}}{\partial x}+</math><math>\frac{\partial {\phi }^{1}}{\partial y}\frac{\partial {\eta }^{1}}{\partial y}+\frac{\partial {\phi }^{1}}{\partial x}\frac{\partial {\zeta }^{1}}{\partial x}+</math><math>\frac{\partial {\phi }^{1}}{\partial y}\frac{\partial {\zeta }^{1}}{\partial y}+</math><math>\frac{\partial {\psi }^{1}}{\partial x}\frac{\partial {\eta }^{1}}{\partial x}+\frac{\partial {\psi }^{1}}{\partial y}\frac{\partial {\eta }^{1}}{\partial y},</math>
|  colspan='2'  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref480536395'></span>(6)
|}


where superscripts <math display="inline">1</math> and <math display="inline">1+</math><math>2</math> denote the components at the first-order and up-to-second-order solution (first-order plus second-order); <math display="inline">{\psi }^{1+2}</math> is the up-to-second-order incident wave velocity potential; <math display="inline">{\zeta }^{1+2}</math> is the up-to-second-order incident free surface elevation; <math display="inline">{\phi }^{1+2}</math> and <math display="inline">{\eta }^{1+2}</math> are the up-to-second-order diffraction-radiation wave velocity potential and free surface elevation; <math display="inline">{P}_{fs}</math> is the free surface pressure; <math display="inline">{S}_{B}^{0}</math> is the mean wetted body surface; <math display="inline">{v}_{p}^{i}</math> is the local velocity induced at point P by the body i<sup>th</sup> order movements; <math display="inline">{v}_{\phi }^{i}</math> is the i<sup>th</sup> order fluid velocity induced by the diffracted-radiated waves; <math display="inline">{v}_{\psi }^{i}</math> is the i<sup>th</sup> order fluid velocity induced by the incident waves; <math display="inline">{n}_{p}^{i}</math> is the normal vector to the body surface <math display="inline">{S}_{B}^{i}</math> at point P.


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">''' [[Image:Draft_Garcia-Espinosa_291836717-image1.jpeg|600px]] '''</span></div>

<span id='_Ref430769885'></span><div id="_Toc429661855" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Figure 1: First and second-order rigid body movement.</span></div>


The fluid pressure at a point P on the body surface is given by:

{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math display="inline">{P}_{p}^{1+2}={P}_{ph}^{0}+{P}_{ph}^{1+2}{+P}_{p\psi }^{1+2}+</math><math>{P}_{p\phi }^{1+2}</math>, 
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(7)
|}


where <math display="inline">{P}_{ph}^{i}</math>, <math display="inline">{P}_{p\psi }^{i}</math>, and <math display="inline">{P}_{p\phi }^{i}</math> stand for the i<sup>th</sup> order hydrostatic, incident wave induced, and diffracted-radiated waves induced pressures. Each pressure component is further decomposed as:

{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math display="inline">{P}_{ph}^{0}=-\rho g{z}_{p}</math>,

<math display="inline">{P}_{ph}^{1+2}=-\rho g{r}_{pz}^{1+2}</math>,

<math display="inline">{P}_{p\phi }^{1+2}=-\rho \frac{\partial {\phi }^{1+2}}{\partial t}-</math><math>\rho \frac{1}{2}\nabla {\phi }^{1}\cdot \nabla {\phi }^{1}-\rho \nabla {\psi }^{1}\cdot \nabla {\phi }^{1}-</math><math>\rho {r}_{p}^{1}\cdot \nabla \left( \frac{\partial {\phi }^{1}}{\partial t}\right) ,</math> 
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(8)
|}


<span id='_Ref466991951'></span><span id='_Ref466538730'></span>where <math display="inline">{r}_{p}^{i}</math> represents the displacement of point P induced by the i<sup>th</sup> order body movement (see Figure 1). Further details on obtaining the governing equations can be found in Servan-Camas and Garcia-Espinosa [<span id='cite-16'></span>[[#16|16]]], and Servan-Camas [<span id='cite-17'></span>[[#17|17]]].

The body dynamics of the floating body are governed by the equation of motion:

{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math display="inline">\overline{\overline{\boldsymbol{M}}}{\boldsymbol{\, X}}_{tt}+</math><math>\overline{\overline{\boldsymbol{K}}}\boldsymbol{\, X}=\boldsymbol{F}</math> 
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(9)
|}


where <math display="inline">\overline{\overline{\boldsymbol{M}}}</math> is the mass matrix of the body; <math display="inline">\overline{\overline{\boldsymbol{K}}}</math> is the hydrostatic restoring matrix (approximates the integral of the hydrostatic pressure); <math display="inline">\boldsymbol{F}</math> is the vector of the hydrodynamic forces induced by dynamic pressures plus any other external forces; and <math display="inline">\boldsymbol{X}</math> represents the movements of the six degrees of freedom of the body.

Loads acting on the body are obtained by direct pressure integration on the body surface underneath the mean water level, with the exception of hydrostatic forces, which are obtained via the corresponding hydrostatic restoring matrices. Also, the second-order loads ( <math display="inline">{\mathit{\boldsymbol{F}}}_{wl}^{2}</math> and <math display="inline">{\mathit{\boldsymbol{M}}}_{wl}^{2}</math>) due to the change of the wetted surface induced by the first order solution are accounted for:

{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{\mathit{\boldsymbol{F}}}_{wl}^{2}\mathit{\boldsymbol{=-}}\frac{1}{2}\rho g\int_{{\Gamma }_{wl}^{0}}^{}{\left( {\xi }^{1}-{r}_{pz}^{1}\right) }^{2}\frac{{\mathit{\boldsymbol{n}}}_{p}^{0}}{\sqrt{1-{{n}_{pz}^{0}}^{2}}}dl,</math>
|}
|  rowspan='2' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(10)
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{\mathit{\boldsymbol{M}}}_{wl}^{2}=-\frac{1}{2}\rho g\int_{{\Gamma }_{wl}^{0}}^{}{\left( {\xi }^{1}-{r}_{pz}^{1}\right) }^{2}\, \overset{\rightarrow}{{\mathit{\boldsymbol{G}}}^{0}{\mathit{\boldsymbol{P}}}^{0}}\times \frac{{\mathit{\boldsymbol{n}}}_{p}^{0}}{\sqrt{1-{{n}_{pz}^{0}}^{2}}}dl.</math>
|}
|}


<span id='_Toc374089684'></span>where <math display="inline">\overset{\rightarrow}{{\mathit{\boldsymbol{G}}}^{0}{\mathit{\boldsymbol{P}}}^{0}}</math> is the vector from the centre of gravity of the floater G to any point P on the wet surface.

===3. Numerical model===

<span id='_Toc435436078'></span>

====3.1 Finite element formulation====

This section presents the formulation based on the Finite Element Method (FEM) to solve the system of equations governing the wave diffraction-radiation problem. This formulation has been developed to be used in conjunction with unstructured meshes in order to enhance geometry flexibility and speed up the initial modelling time.

Let <math display="inline">{Q}_{h}^{\ast }</math> be the finite element space for interpolating functions, when constructed in the usual manner. From this space, it can be constructed that subspace <math display="inline">{Q}_{h,\phi }^{\ast }</math>, which incorporates the Dirichlet conditions for the potential <math display="inline">\phi</math> . The space of test functions, denoted by <math display="inline">{Q}_{h}</math>, is constructed as <math display="inline">{Q}_{h,\phi }</math>, but with functions vanishing on the Dirichlet boundary. The weak form of the problem can be written as follows:

Find <math display="inline">[{\phi }_{h}]\in {Q}_{h,\phi }^{\ast }</math>, by solving the discrete variational problem:

{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>\int_{\Omega }^{}\nabla {\upsilon }_{h}\cdot \nabla {\phi }_{h}\, d\Omega =\int_{{\Gamma }^{B}}^{}{\upsilon }_{h}\cdot {\hat{\phi }}_{n}^{B}\, d\Gamma +</math><math>\int_{{\Gamma }^{R}}^{}{\upsilon }_{h}\cdot {\hat{\phi }}_{n}^{R}\, d\Gamma</math>

<math>+\int_{{\Gamma }^{{Z}_{0}}}^{}{\upsilon }_{h}\cdot {\hat{\phi }}_{n}^{{Z}_{0}}d\Gamma +\int_{{\Gamma }^{{Z}_{-H}}}^{}{\upsilon }_{h}\cdot {\hat{\phi }}_{n}^{{Z}_{-H}}d\Gamma \quad \, {\forall \upsilon }_{h}\in {Q}_{h},</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref465845025'></span>(11)
|}


where <math display="inline">{\hat{\phi }}_{n}^{B}</math>, <math display="inline">{\hat{\phi }}_{n}^{R}</math>, <math display="inline">{\hat{\phi }}_{n}^{{Z}_{0}}</math> and <math display="inline">{\hat{\phi }}_{n}^{{Z}_{-H}}</math> are the potential normal gradients corresponding to the Neumann boundary conditions on bodies; the radiation boundary; the free surface; and the bottom surface of the domain. At this point, it is useful to introduce the associated matrix form of Eq.<span id='cite-_Ref465845025'></span>[[#_Ref465845025|(11)]]:

{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>\overline{\overline{\boldsymbol{L}}}\boldsymbol{\phi }={\boldsymbol{b}}^{B}+{\boldsymbol{b}}^{R}+</math><math>{\boldsymbol{b}}^{{Z}_{0}}+{\boldsymbol{b}}^{{Z}_{-H}},</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(12)
|}


where <math display="inline">\overline{\overline{\boldsymbol{L}}}</math>''' '''is the standard Laplacian matrix, and <math display="inline">{\boldsymbol{b}}^{B}</math>''', ''' <math display="inline">{\boldsymbol{b}}^{R}</math>''', '''and''' ''' <math display="inline">{\boldsymbol{b}}^{{Z}_{0}}</math>, and''' ''' <math display="inline">{\boldsymbol{b}}^{{Z}_{-H}}</math>''' '''are the vectors resulting from integrating the corresponding boundary condition terms. The seabed boundary for the refracted and radiated potential is imposed by taking <math display="inline">{\boldsymbol{b}}^{{Z}_{-H}}=</math><math>0</math>.

<span id='_Toc435436079'></span>

====3.2 Free surface boundary conditions====

Combining the kinematic (Eq. <span id='cite-_Ref480536160'></span>[[#_Ref480536160|(2)]]) and dynamic (Eq. <span id='cite-_Ref480536163'></span>[[#_Ref480536163|(3)]]) free surface boundary conditions, the free surface condition up to second order reads:

{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>\frac{{\partial }^{2}\phi }{\partial {t}^{2}}+g\frac{\partial \phi }{\partial z}+</math><math>\frac{\partial }{\partial t}\left( \frac{{P}_{fs}}{\rho }\right) +\left\{ {Q}^{1}\right\} =</math><math>0.</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref465845046'></span>(13)
|}


where superscripts 1+2 have been omitted (and will be from this point on), and <math display="inline">{Q}^{1}</math> are the source terms from the first-order solution. This condition is implemented as a Neumann boundary condition that fulfils the flux boundary integral:

{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math display="inline">{\boldsymbol{b}}^{{Z}_{0}}={\overline{\overline{\boldsymbol{M}}}}_{{\Gamma }^{{Z}_{0}}}{\hat{\phi }}_{n}^{{Z}_{0}}</math>,
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(14)
|}


where <math display="inline">{\overline{\overline{\boldsymbol{M}}}}_{{\Gamma }^{{Z}_{0}}}</math> ''' '''is''' '''the corresponding boundary mass. The terms <math display="inline">{R}^{1}</math> and <math display="inline">{S}^{1}</math> are given by Eqs. <span id='cite-_Ref480536394'></span>[[#_Ref480536394|(5)]] and <span id='cite-_Ref480536395'></span>[[#_Ref480536395|(6)]] respectively, then:

{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math display="inline">{Q}^{1}={\partial }_{t}{R}^{1}-{S}^{1}</math>,
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(15)
|}


Eq. <span id='cite-_Ref465845046'></span>[[#_Ref465845046|(13)]] is discretized in time using the following numerical scheme:

{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>\frac{{\phi }^{n+1}-2{\phi }^{n}+{\phi }^{n-1}}{\Delta {t}^{2}}=-g{\phi }_{z}^{n}-</math><math>\frac{1}{12}g\left( {\phi }_{z}^{n+1}+10{\phi }_{z}^{n}+{\phi }_{z}^{n-1}\right) -</math><math>\frac{{{P}_{fs}^{n+1}-P}_{fs}^{n-1}}{\rho 2\Delta t}</math>

<math>\quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad -\left\{ \frac{1}{12}\left( {\left( {Q}^{1}\right) }^{n+1}+\right. \right. </math><math>\left. \left. 10{\left( {Q}^{1}\right) }^{n}+{\left( {Q}^{1}\right) }^{n-1}\right) \right\} ,</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(16)
|}


where for the specific case <math display="inline">{P}_{fs}=0</math>, the above scheme becomes a fourth order compact Padé scheme. Once the velocity potential is solved at the new time step, the free surface elevation is computed using the following fourth-order in time numerical scheme:

{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{\eta }^{n+1}=-\frac{1}{g\Delta t\, }\left( {\frac{25}{12}{\phi }^{n+1}-4\phi }^{n}+\right. </math><math>\left. 3{\phi }^{n-1}-\frac{4}{3}{\phi }^{n-2}+\frac{1}{4}{\phi }^{n-3}\right) -</math><math>\frac{{P}_{fs}^{n+1}}{\rho g}\left\{ -{\left( {S}^{1}\right) }^{n+1}\right\} .</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(17)
|}


<span id='_Toc435436080'></span>

====3.3 Radiation condition and wave absorption====

Waves represented by <math display="inline">\phi</math>  are born at the bodies and propagate in all directions away from them. The waves need to either be dissipated or to leave the domain so they will not bounce back at the outer boundary and interact with the bodies. A Sommerfeld radiation condition at the edge of the computational domain is introduced:

{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{\partial }_{t}\phi +c\nabla \phi \cdot {n}_{R}=0</math> <math display="inline">in\, {\Gamma }_{R}</math>,
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(18)
|}


where <math display="inline">{\Gamma }_{R}</math> is the surface bounding the domain horizontally; <math display="inline">{n}_{R}</math> is the normal vector of <math display="inline">{\Gamma }_{R}</math> pointing outwards from the domain; and <math display="inline">c</math> is a prescribed wave phase velocity. This radiation condition will let waves moving at velocity <math display="inline">c</math> to escape the domain. The numerical scheme used to implement the radiation condition is

{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{\left( {\phi }_{\mathit{\boldsymbol{n}}}^{R}\right) }^{n+1}=\, -\frac{{\phi }^{n-1}-{\phi }^{n}}{c\Delta t}</math> <math display="inline">in\, {\Gamma }_{R}</math>''.''
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(19)
|}


The prescribed phase velocity <math display="inline">c</math> will be set for radiating those waves with the smallest frequency (largest wavelengths) considered in each specific case. The typical value of <math display="inline">c</math> is the phase velocity of longer incident waves. However, waves with higher frequencies (smaller phase velocities) will not leave the domain through <math display="inline">{\Gamma }_{R}</math>, so that they will be reflected. Therefore, wave absorption is introduced into the dynamic free surface boundary condition by varying the pressure such that:

{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{P}_{fs}\left( \mathit{\boldsymbol{x,}}t\right) =\kappa \left( \mathit{\boldsymbol{x}}\right) \rho \frac{\partial \phi }{\partial z}.</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref465845106'></span>(20)
|}


Eq. <span id='cite-_Ref465845106'></span>[[#_Ref465845106|(20)]] increases pressure when the free surface is moving upwards, while decreasing the pressure when the free surface is moving downwards. Then energy is transferred from the waves to the atmosphere and waves are damped. However, the coefficient <math display="inline">\kappa (\mathit{\boldsymbol{x}})</math> will be set to zero in the analysis area (near the bodies), so that damping will have no effect on the solution of the wave-body interaction problem. Further details can be found in Servan-Camas and Garcia-Espinosa [<span id='cite-_Ref466991951'></span>[[#_Ref466991951|16]]] and Servan-Camas [<span id='cite-_Ref466538730'></span>[[#_Ref466538730|17]]].

===4. Mooring models===

<span id='_Ref466539385'></span><span id='_Ref466539367'></span>Two different computational models have been implemented in the seakeeping solver to simulate the mooring lines. The first one is an elastic catenary solver, and the second is a non-linear FEM dynamic cable solver. In the following sections, a brief description of the mathematical model is given. Details of the numerical implementation of the mooring solver, the body dynamics solver, and their coupling can be found in [<span id='cite-_Ref466538730'></span>[[#_Ref466538730|17]],<span id='cite-18'></span>[[#18|18]],<span id='cite-19'></span>[[#19|19]],<span id='cite-'></span>[[#|]]].

====4.1 Elastic Catenary model====

The elastic catenary formulation is based on the model proposed in [<span id='cite-20'></span>[[#20|20]]]. Each mooring line is analysed in a local coordinate system located at the anchor. The local z-axis is oriented vertical and the local x-axis is oriented horizontally from the anchor to the position of the fairlead. When a portion of the mooring line rests on the seabed, the equations for the horizontal and vertical distances between the anchor and a given point on the line, <math display="inline">x\left( s\right)</math>  and <math display="inline">z\left( s\right)</math> , can be written as,

{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>x\left( s\right) =\left\{ \begin{matrix}s&for\, \, \, 0\leq s\leq {L}_{B}-\frac{{H}_{F}}{{C}_{B}\omega },\\s+\frac{{C}_{B}\omega }{2EA}\left[ {s}^{2}-2\left( {L}_{B}-\frac{{H}_{F}}{{C}_{B}\omega }\right) s+\left( {L}_{B}-\frac{{H}_{F}}{{C}_{B}\omega }\right) \mathrm{max}\,\left( \left( {L}_{B}-\frac{{H}_{F}}{{C}_{B}\omega }\right) ,0\right) \right] &for\, {L}_{B}-\frac{{H}_{F}}{{C}_{B}\omega }\leq s\leq {L}_{B},\\\begin{matrix}{L}_{B}+\frac{{H}_{F}}{\omega }\mathrm{ln}\,\left[ \frac{\omega \left( s-{L}_{B}\right) }{{H}_{F}}+\sqrt{1+{\left( \frac{\omega \left( s-{L}_{B}\right) }{{H}_{F}}\right) }^{2}\, }\right] +\frac{{H}_{F}L}{EA}\\+\frac{{C}_{B}\omega }{2EA}\left[ -{L}_{B}^{2}+\left( {L}_{B}-\frac{{H}_{F}}{{C}_{B}\omega }\right) \mathrm{max}\,\left( \left( {L}_{B}-\frac{{H}_{F}}{{C}_{B}\omega }\right) ,0\right) \right] \end{matrix}&for\quad \quad \quad \quad \, {L}_{B}\leq s\leq L,\end{matrix}\right.</math> 
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(21)
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>z\left( s\right) =\left\{ \begin{matrix}0&for\, \, \, 0\leq s\leq {L}_{B}-\frac{{H}_{F}}{{C}_{B}\omega },\\\frac{{H}_{F}}{\omega }\sqrt{1+{\left( \frac{\omega \left( s-{L}_{B}\right) }{{H}_{F}}\right) }^{2}\, }+\frac{\omega {\left( s-{L}_{B}\right) }^{2}}{2EA}&for\quad \quad \quad \quad \, {L}_{B}\leq s\leq L,\end{matrix}\right.</math> 
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(22)
|}


being <math display="inline">s</math> the catenary arc length; <math display="inline">{H}_{F}</math> the horizontal component of the effective tension; <math display="inline">{C}_{B}</math> static friction coefficient; <math display="inline">\omega</math>  catenary weight per unit length in water; E the Young modulus; A the cross section area; and <math display="inline">{L}_{B}</math> portion of the length of cable resting on the seabed.

The equation for the effective tension in the line at any point of the line <math display="inline">T\left( s\right)</math>  is written as follows:

{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math display="inline">T\left( s\right) =\left\{ \begin{matrix}\mathrm{max}\,\left( {H}_{F}+{C}_{B}\omega \left( s-{L}_{B}\right) ,0\right) \\\sqrt{{H}_{F}^{2}+\left( \omega {\left( s-{L}_{B}\right) }^{2}\right) }\quad for\quad {L}_{B}\leq s\leq L\end{matrix}\right.</math> .
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(23)
|}


The above formula calculates the total load on the system from all of mooring lines. The restoring load is found by first translating each fairlead tension from its local coordinate system to the global frame, and then adding the tensions from all lines.

====4.2 Dynamic cable model====

The dynamic equations for a mooring cable with length L, negligible bending and torsional stiffness can be formulated as [<span id='cite-_Ref466539385'></span>[[#_Ref466539385|18]]]:

{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>\left( {\rho }_{w}{C}_{m}{A}_{0}+{\rho }_{0}\right) \frac{{\partial }^{2}{r}_{l}}{\partial {t}^{2}}=</math><math>\frac{\partial }{\partial l}\left( E{A}_{0}+\frac{e}{e+1}\frac{\partial {r}_{l}}{\partial l}\right) +</math><math>f\left( t\right) \left( 1+e\right) ,\,</math> 
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(24)
|}


where <math display="inline">{\rho }_{w}</math> is the water density; <math display="inline">{C}_{m}</math> is the added mass coefficient; <math display="inline">{\rho }_{0}</math> is the mass per unit length of the unstretched cable; <math display="inline">\, {r}_{l}</math> is the position vector; <math display="inline">E</math> is the Young's modulus; <math display="inline">{A}_{0}</math> is the cross-sectional area of the cable; <math display="inline">e</math> is the strain; <math display="inline">f\left( t\right)</math>  are the external loads applied on the cable including the self-weight; hydrostatic loads; drag forces and seabed interaction; and <math display="inline">l</math> is the length along the unstretched cable.

The boundary conditions of the mooring cable are given by

{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>\frac{{\partial }^{2}{r}_{l}}{\partial {t}^{2}}=0,\, at\, l=0,\, \left( anchor\, point\right) ,</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(25)
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>\frac{{\partial }^{2}{r}_{l}}{\partial {t}^{2}}=\ddot{{r}_{b}},\, at\, l=L,\, \left( fairlead\, point\right) ,</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(26)
|}


where <math display="inline">\ddot{{r}_{b}}</math> is the second derivative of the position vector at the fairlead connection point.

The above non-linear equation is solved using the standard Finite Element Method. Details of the mathematical and numerical dynamic model are provided in Gutierrez-Romero et al. [<span id='cite-_Ref466539385'></span>[[#_Ref466539385|18]]].

===5. Numerical model verification===

====5.1 Verification case 1: mean drift forces on a hemisphere====

=====5.1.1 Problem description=====

This case consists of estimating the mean drift forces on a hemisphere. In this work, mean drift forces are obtained by averaging the time series of the corresponding second-order force. The analytical solution for the fixed hemisphere was obtained by Fernandes and Levy [<span id='cite-21'></span>[[#21|21]]], and for the freely floating hemisphere was obtained by Kudou [<span id='cite-22'></span>[[#22|22]]] and reported by Pinkster [<span id='cite-23'></span>[[#23|23]]]. In this section, the numerical results are compared with the analytical solution. The details of the hemisphere particulars are given in <span id='cite-_Ref426102176'></span>[[#_Ref426102176|Table 1]]. Both fixed and free floating cases were analysed. Figure 2 (left) shows the mesh used for the calculations. It can be observed that mesh refinement is required in the area of the waterline in order to obtain accurate results of mean drift forces. Figure 2 (right) shows a snapshot of the wave elevation around the hemisphere in one of the cases run. Finally, figure 3 (above) compares the analytical results with the numerical results in the case of the fixed hemisphere, while figure 3 (below) compares the results in the free floating case. Good correlation is observed for the whole range of waves analysed.

<div id="_Ref426102176" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Table 1: hemisphere particulars.</span></div>

{| style="width: 38%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Depth'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Infinite</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Mass'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Displacement</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Radius'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1 m</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''CG coordinates '''[m]</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(0,0,0)</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Number of tetrahedrons'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">274283</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Number of triangles'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">22082</span>
|}


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
{|
|-
| [[Image:Draft_Garcia-Espinosa_291836717-image2.jpeg|234px]]
| [[Image:Draft_Garcia-Espinosa_291836717-image3.jpeg|center|294px]]
|}
</div>

<span id='_Ref426102204'></span><div id="_Toc429661876" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Figure 2: Hemisphere: mesh refinement close up (left) and wave contours (right).</span></div>


{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="text-align: center;vertical-align: top;width: 100%;"|[[Image:Draft_Garcia-Espinosa_291836717-image4.png|468px]]

[[Image:Draft_Garcia-Espinosa_291836717-image5.png|468px]] 
|-
|  style="text-align: center;vertical-align: top;"|<span id='_Ref426102205'></span><span id='_Toc429661877'></span>

<span style="text-align: center; font-size: 75%;">Figure 3: Mean drift forces on Hemisphere.</span>
|}


=====5.1.2 Verification case 2: diffraction of second-order monochromatic waves by semisubmerged horizontal rectangular cylinder=====

<span id='_Toc435436133'></span>

===1.7.1 Problem description===

<span id='_Ref426107620'></span>This test case deals with the solution to the diffraction problem of second-order monochromatic surface waves by a semisubmerged horizontal cylinder of rectangular cross section. The boundary-value problem is resolved and the results are compared with the analytical solution obtained by the method of matched Eigen function expansions presented in [<span id='cite-24'></span>[[#24|24]]]. Horizontal and vertical forces and the moment about the heel of the prismatic cylinder are analysed for different monochromatic waves. An illustration of this problem is shown in figure 4. Relevant geometric parameters are: depth (h = 1 m), half beam (b = 1 m), and draught (d = 0.2 m).

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Garcia-Espinosa_291836717-image6.jpeg|600px]] </div>

<span id='_Ref426106310'></span><div id="_Toc429661878" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Figure 4: Horizontal semi-submerged cylinder: problem layout.</span></div>

The situation considered for analysis is the diffraction of waves by a fixed horizontal cylinder of rectangular cross section. Analysis is undertaken with the following assumptions: the fluid is inviscid and incompressible, the seabed and the cylinder are impervious, and the excitation is provided by normally incident plane waves of small amplitude and frequency. Several cases are run for different wave periods (T = 0.897, 1.003, 1.07, 1.16, 1.445, 2.299, 4.17, and 6.37 seconds), and the simulation time is about 30 seconds, with an initialization time of 10 seconds. All movement is restrained so that the body is completely fixed, therefore, wave diffraction occurs but not radiation.

===1.7.2 Mesh generation===

Mesh properties for the present analysis are summarized in <span id='cite-_Ref426106512'></span>[[#_Ref426106512|Table 2]]. Figure 5 shows an isometric view of the mesh used for the present analysis at the region close to the surface of the body.

<div id="_Ref426106512" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Table 2: Verification case 2: mesh particulars.</span></div>

{| style="margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Minimum element size'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.01 m</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Maximum element size'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1 m</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Number of elements'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">121687</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Number of nodes'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">22940</span>
|}


<div id="_Toc435436135" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Garcia-Espinosa_291836717-image7.jpeg|600px]] </div>

<span id='_Ref426107505'></span><div id="_Toc429661879" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Figure 5: Generated mesh: Top view (up) and front view (down).</span></div>

===1.7.3 Verification of results===

Figure 6 shows the amplitude of the second-order horizontal force <math display="inline">({F}_{x})</math>, vertical force <math display="inline">({F}_{z})</math> forces, and moment about y axis <math display="inline">{(M}_{y})</math>, for both the analytical results reported in [<span id='cite-_Ref426107620'></span>[[#_Ref426107620|25]]] and the numerical results obtained in this work. The second-order components of the forces and moments (double frequency component) are normalized with <math display="inline">\rho gh{A}^{2}</math>, where <math display="inline">\rho</math>  is the density of the fluid, <math display="inline">g</math> the gravity, <math display="inline">{A}^{2}</math> the square of the wave amplitude, and <math display="inline">h</math> the water depth. Results are plotted against the dimensionless wave number ( <math display="inline">kh</math>). As can be observed, good correlation is obtained for the analysed range of wave numbers.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Garcia-Espinosa_291836717-image8.png|600px]] </div>

<span id='_Ref426107589'></span><div id="_Toc429661880" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Figure 6: Second-order wave forces (r=0) and moment (r=1) on a rectangular cylinder: analytical versus numerical.</span></div>

====1.8 Verification case 3: diffraction of second-order bichromatic waves by bottom mounted circular cylinder====

===1.8.1 Problem description===

<span id='_Ref481314099'></span>This test case deals with the diffraction of a bottom mounted vertical cylinder under monochromatic waves with water depth <math display="inline">h=</math><math>4\, m</math>, radius of the cylinder <math display="inline">R=1\, m</math>. More details can be found on [<span id='cite-25'></span>[[#25|25]]].

===1.8.2 Convergence analysis===

<span id='OLE_LINK53'></span><span id='OLE_LINK54'></span><span id='OLE_LINK55'></span><span id='OLE_LINK72'></span><span id='OLE_LINK73'></span><span id='OLE_LINK74'></span>A convergence analysis has been carried out to assess the convergence of the present numerical approach to the mathematical model. <span id='cite-_Ref481313867'></span>[[#_Ref481313867|Table 3]] provides the unstructured mesh details for each case tested. The wave frequencies chosen for this analysis are <math display="inline">{\omega }_{1}^{2}R/g=</math><math>1.4</math> and  <math display="inline">{\omega }_{2}^{2}R/g=2.0</math>. Time step is calculated such that <math display="inline">g\Delta {t}^{2}/h=</math><math>0.34</math>, being <math display="inline">h</math> the characteristic element size at the floating line.

<span id='cite-_Ref481313868'></span>[[#_Ref481313868|Table 4]] provides the dimensionless force for the sum-frequency along with the wave direction, the difference between two consecutive discretization sizes, and the relative difference respect to the corresponding force. Differences have been calculated instead of errors due to the lack of an exact solution to compare with. The results are given in <span id='cite-_Ref481313868'></span>[[#_Ref481313868|Table 4]] and show that the convergence rate is approximately of second order, although care must be taken since convergence tests on unstructured meshes contain uncertainties due to the irregularities in the shape of the elements.

<div id="_Ref481313867" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Table 3: Convergence test: mesh particulars.</span></div>

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;"|
|  colspan='3'  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Characteristic element size [m]</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|
|  rowspan='2' style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Number of elements</span>
|  rowspan='2' style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Number of nodes</span>
|-
|  style="border: 1pt solid black;text-align: center;"|
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Floating line </span> <math display="inline"></math> <math>h</math>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Body and free surface</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Volume</span>
|  style="border: 1pt solid black;text-align: center;"|<span id='OLE_LINK59'></span><span id='OLE_LINK60'></span><span id='OLE_LINK61'></span> <math display="inline">\Delta t</math><span style="text-align: center; font-size: 75%;"> [s]</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span id='_Hlk481400529'>Mesh 1</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.4</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">0.0587</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">727221</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">125926</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Mesh 2</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.071</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1414</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2828</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">0.0494</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">983719</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">169908</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Mesh 3</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.05</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">0.0415</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1683943</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">289527</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Mesh 4</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.035</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0707</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1414</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">0.0349</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3475047</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">594852</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Mesh 5</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.025</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.05</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">0.0294</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">8203778</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1399855</span>
|}


<span id='_Ref481313868'></span><span id='OLE_LINK29'></span><div id="OLE_LINK31" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Table 4: Convergence analysis: values of the sum-frequency dimensionless force for </span> <math display="inline">{\omega }_{1}^{2}R/g=</math><math>1.4</math><span style="text-align: center; font-size: 75%;"> and </span> <math display="inline">{\omega }_{2}^{2}R/g=</math><math>2.0</math><span style="text-align: center; font-size: 75%;">.</span></div>

{| style="width: 65%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;"|
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Mesh 1</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Mesh 2</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Mesh 3</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Mesh 4</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Mesh 5</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span id='OLE_LINK64'></span><span id='OLE_LINK65'></span> <math>h</math>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.071</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.05</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.035</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.025</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<math>{F}_{x}^{\ast }</math>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.498</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.348</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.229</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.159</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.141</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Difference:</span>

<math display="inline">\epsilon =\left| {\left( {F}_{x}^{\ast }\right) }^{i}-{\left( {F}_{x}^{\ast }\right) }^{i-1}\right|</math> 
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">-</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.150</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.119</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.070</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.018</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Relative difference:</span>

<math display="inline">{\epsilon }^{\%}=100\, \epsilon /\left| {\left( {F}_{x}^{\ast }\right) }^{i}\right| \,</math> 
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">-</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">6.39%</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.34%</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.24%</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.84%</span>
|}

===1.8.3 Verification of results===

<span id='cite-_Ref480788660'></span>[[#_Ref480788660|Table 5]] shows a comparison for the non-dimensional amplitude for the sum-frequency and difference-frequency forces along the wave direction. In the previous analysis, Mesh 4 allows to obtain a relative difference below of 5% (see <span id='cite-_Ref481313868'></span>[[#_Ref481313868|Table 4]]), therefore it has been selected to be used in this verification. Compared to numerical and semi-analytical solutions obtained by other authors, the results in this work are within the range of those reported elsewhere.

<div id="_Ref480788660" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Table 5: Comparison of non-dimensional amplitude of the sum-frequency and difference-frequency forces along the wave direction.</span></div>

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  rowspan='2' style="border: 1pt solid black;text-align: center;"|
|  colspan='2'  style="border: 1pt solid black;text-align: center;"|<math display="inline">\frac{{\omega }_{1}^{2}R}{g}=</math><math>1.0</math> <math display="inline">\frac{{\omega }_{2}^{2}R}{g}=1.6</math>
|  colspan='2'  style="border: 1pt solid black;text-align: center;"|<math display="inline">\frac{{\omega }_{1}^{2}R}{g}=</math><math>1.2</math> <math display="inline">\frac{{\omega }_{2}^{2}R}{g}=1.8</math>
|  colspan='2'  style="border: 1pt solid black;text-align: center;"|<math display="inline">\frac{{\omega }_{1}^{2}R}{g}=</math><math>1.4</math> <math display="inline">\frac{{\omega }_{2}^{2}R}{g}=2.0</math>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">sum-freq</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">diff-freq</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">sum-freq</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">diff-freq</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">sum-freq</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">diff-freq</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Shao and Faltinsen [<span id='cite-26'></span>[[#26|26]]]</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.868</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.861</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.190</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.788</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.088</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.759</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Kim and Yue [<span id='cite-27'></span>[[#27|27]]]</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.853</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.856</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.182</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.788</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.094</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.765</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span id='OLE_LINK11'></span><span id='OLE_LINK12'>Eatock Taylor and Huang [</span><span id='cite-_Ref481314099'></span>[[#_Ref481314099|<span style="text-align: center; font-size: 75%;">26</span>]]<span style="text-align: center; font-size: 75%;">]</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.883</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.849</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.294</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.769</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.114</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.777</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Moubayed and Williams [<span id='cite-28'></span>[[#28|28]]]</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.783</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.840</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.091</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.761</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.998</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.734</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Present work</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.815</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.852</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.248</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.765</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.159</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.740</span>
|}


===2. Validation against the HiPRWind model===

====2.1 Case description====

The floating platform geometry considered in this paper has been provided by the HiPRWind FP7 project (EU 7th RTD FP under grant agreement no. 256812) [<span id='cite-29'></span>[[#29|29]]] and is composed of three buoyant columns connected with bracings. Model tests were carried out at Ecole Centrale Nantes’ facilities. A model built in stainless steel to scale of <math display="inline">\lambda =</math><math>1/19.8</math> was used in the tests (see figure 7). These experiments were devised by Simos et al. [<span id='cite-_Ref482260439'></span>[[#_Ref482260439|15]]] in order to validate an alternative frequency-domain method to obtain the second-order response of the floater from the measured motions. Results related to mooring loads are presented here for the first time.

Figure 8 shows an overview of the HiPRWind CAD model generated. <span id='cite-_Ref466541966'></span>[[#_Ref466541966|Table 6]] provides the platform details in full scale, as well as the water depth used in this study.

<div id="_Ref466541966" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Table 6: HiPRWind platform main particulars (full scale).</span></div>

{| style="width: 42%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Depth'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">100 m</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Operation design draft'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">15.5 m</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Distance between column centres'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">35 m</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Column diameter'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">7 m</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Heave plates diameter'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">20 m</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Displacement'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2332 T</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''XG'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0 m</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''YG'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0 m</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''ZG'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">-4.46 m</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Radius of gyration (pitch)'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">22.38 m</span>
|}


Figure 9 shows a view of the mesh used for this case study. This mesh consists of tetrahedral 643603 elements and 119350 nodes. The cylindrical domain has a radius of 500 meters, a height of 100 m water depth, and the absorption area starts 50 meters from the centre of the platform.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Garcia-Espinosa_291836717-image9.jpeg|468px]] </div>

<div id="_Ref465845190" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Figure 7: HiPRWind platform model used.</span></div>

====2.2 Model calibration====

In order to predict seakeeping in real conditions with a potential flow solver, viscous effects are to be incorporated via external forces. These external forces are modelled by simplified formulas accounting for the overall viscous effects acting on the platform. In this work, the viscous effects have been included in the computational solver by means of linear and quadratic models. This model has been calibrated using the experimental information of the extinction tests for surge, heave and pitch motions.

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|''' [[Image:Draft_Garcia-Espinosa_291836717-image10.jpeg|366px]]  [[Image:Draft_Garcia-Espinosa_291836717-image11.jpeg|366px]] '''
|}


<div id="_Ref466541894" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Figure 8: HiPRWind computational model overview.</span></div>

Figure 10 shows the layout of the experimental set-up for these tests. Three elastic lines were used to keep the model in position during the extinction experiments. These lines have a small linear weight distribution so the catenary effect is negligible. The location of the end points for each line is given in <span id='cite-_Ref466360460'></span>[[#_Ref466360460|Table 7]]. Pretension of 550kN was applied to each line.

<div id="_Ref466360460" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Table 7: Mooring lines end points coordinates and stiffness coefficients.</span></div>

{| style="width: 61%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;"|
|  style="border: 1pt solid black;text-align: center;"|<math display="inline">\mathit{\boldsymbol{(}}{\mathit{\boldsymbol{x}}}_{\mathit{\boldsymbol{i}}}\mathit{\boldsymbol{,}}{\mathit{\boldsymbol{y}}}_{\mathit{\boldsymbol{i}}}\mathit{\boldsymbol{,}}{\mathit{\boldsymbol{z}}}_{\mathit{\boldsymbol{i}}}\mathit{\boldsymbol{)}}</math><span style="text-align: center; font-size: 75%;">''' '''[m]</span>
|  style="border: 1pt solid black;text-align: center;"|<math display="inline">\mathit{\boldsymbol{(}}{\mathit{\boldsymbol{x}}}_{\mathit{\boldsymbol{f}}}\mathit{\boldsymbol{,}}{\mathit{\boldsymbol{y}}}_{\mathit{\boldsymbol{f}}}\mathit{\boldsymbol{,}}{\mathit{\boldsymbol{z}}}_{\mathit{\boldsymbol{f}}}\mathit{\boldsymbol{)}}</math><span style="text-align: center; font-size: 75%;">''' '''[m]</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{K\, }}[KN/m]</math>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Line 1'''</span>
|  style="border: 1pt solid black;text-align: center;"|<math>(-323.75,\, 0,\, 10)</math>
|  style="border: 1pt solid black;text-align: center;"|<math>(-23.73,\, 0,\, 10)</math>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">20.8</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Line 2'''</span>
|  style="border: 1pt solid black;text-align: center;"|<math>(207.7,250.4,10)</math>
|  style="border: 1pt solid black;text-align: center;"|<math>(11.87,\, 20.55,\, 10)</math>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">20.8</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Line 3'''</span>
|  style="border: 1pt solid black;text-align: center;"|<math>(207.7,-250.4,10)</math>
|  style="border: 1pt solid black;text-align: center;"|<math>(11.87,\, -20.55,\, 10)</math>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">20.8</span>
|}


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
''' [[Image:Draft_Garcia-Espinosa_291836717-image12.jpeg|600px]] '''</div>

<div id="_Ref466360426" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Figure 9: FEM model mesh.</span></div>

The linear stiffness matrix corresponding to the mooring system is:
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{K}_{mooring}=\left[ \begin{matrix}36.5KN/m&0\, KN/m&0\, KN/m&0\, KN/rad&36.5KN/rad&0\, KN/rad\\0\, KN/m&23.5KN/m&0\, KN/m&-235KN/rad&0\, KN/rad&126KN/rad\\0\, KN/m&0\, NK/m&0\, KN/m&0\, KN/rad&0\, KN/rad&0\, KN/rad\\0\, KNm/m&-235KNm/m&0\, KNm/m&19600KNm/rad&0\, KNm/rad&-1260KNm/rad\\365KNm/m&0\, KNm/m&0\, KNm/m&0\, KNm/rad&25100KNm/rad&0\, KNm/rad\\0\, KNm/m&126KNm/m&0\, KNm/m&309KNm/rad&0\, KNm/rad&39400KNm/rad\end{matrix}\right]</math> 
|}


The viscous damping forces have been divided into two groups. The first group corresponds to the bracings of the structure. The corresponding forces are applied in the centre of gravity of the platform and depend on the velocity of the platform. The second group corresponds to the heave plates and columns, and calibration forces are applied in the centre of the heave plates using Morison elements that take into account the relative velocity between the platform and the incident waves.  <span id='cite-_Ref419733175'></span>[[#_Ref419733175|Table 8]] provides a summary of the coefficients of the damping terms obtained in the calibration phase. A similar calibration process for the BEM frequency domain solver WADAM [<span id='cite-30'></span>[[#30|30]]] has also been performed by a third party, resulting in reasonably similar calibration values.

<div id="_Ref419733175" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Table 8: Model calibration: added mass, linear damping and quadratic damping coefficients.</span></div>

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;"|
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">FEM</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">WADAM/SIMO</span>
|-
|  rowspan='3' style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Applied at CG'''</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Surge linear damping: '''</span> <math display="inline">{\boldsymbol{B}}_{\boldsymbol{11}}</math><span style="text-align: center; font-size: 75%;">[KN/(m/s)]</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">75 </span>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">70</span>
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Heave added mass: '''</span> <math display="inline">{\boldsymbol{A}}_{\boldsymbol{33}}</math><span style="text-align: center; font-size: 75%;">''' '''[t]</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1200 </span>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">1000</span>
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Heave linear damping: '''</span> <math display="inline">{\boldsymbol{B}}_{\boldsymbol{33}}</math><span style="text-align: center; font-size: 75%;">[KN/(m/s)]</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">110 </span>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">110</span>
|-
|  rowspan='2' style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Applied at the centre of each '''</span>

<span style="text-align: center; font-size: 75%;">'''heave plate base'''</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Heave linear damping: '''</span> <math display="inline">{\boldsymbol{B}}_{\boldsymbol{33}}</math><span style="text-align: center; font-size: 75%;">[KN/(m/s)]</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">76</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">50</span>
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Heave quadratic damping: '''</span> <math display="inline">{\boldsymbol{B}}_{\boldsymbol{33}}^{\boldsymbol{2}}</math><span style="text-align: center; font-size: 75%;">[KN/(m/s)<sup>2</sup>]</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">805</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">600</span>
|}


Figure 11 shows experimental versus numerical results after calibration of the surge, heave, and pitch decay tests. Good correlation has been reached for the three degrees of freedom. The natural periods obtained by the FEM solver are given in <span id='cite-_Ref467179576'></span>[[#_Ref467179576|Table 9]]. Natural periods obtained by the BEM model and experiments have an error margin of 1 second.

<div id="_Ref467179576" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Table 9: Natural periods obtained from decay test with linear restoring.</span></div>

{| style="margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Surge'''</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Heave'''</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Pitch'''</span>
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">70 s</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">19 s</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">26 s</span>
|}


====2.3 Catenary mooring====

Three catenary lines were used as mooring lines for the rest of the experiments. The location of the end points for each line is given in <span id='cite-_Ref419790750'></span>[[#_Ref419790750|Table 10]]. <span id='cite-_Ref466537498'></span>[[#_Ref466537498|Table 11]] provides the mooring line details.

<div id="_Ref419790750" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Table 10: Mooring lines end points coordinates.</span></div>

{| style="width: 58%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;"|
|  style="border: 1pt solid black;text-align: center;"|<math display="inline">\boldsymbol{Fairlead}\boldsymbol{\, }\boldsymbol{point}</math><span style="text-align: center; font-size: 75%;">''' (x,y,z) '''[m]</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Anchor point (x,y,z) '''[m]</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Line 1'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(-325.73, 0.00, -80.00)</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(-23.73, 0.00, 10.00)</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Line 2'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(206.79, -248.00, -80.00)</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(15.24, -18.17, 10.00)</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Line 3'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(206.79, 248.00, -80.00)</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(15.21, 18.17, 10.00)</span>
|}


<div id="_Ref466537498" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Table 11: Mooring lines particulars.</span></div>

{| style="width: 45%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Length'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">330 m</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Section'''</span>
|  style="border: 1pt solid black;text-align: center;"|<math display="inline">1.108E{\times 10}^{-2}</math><span style="text-align: center; font-size: 75%;">m<sup>2</sup></span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Equivalent Young modulus'''</span>
|  style="border: 1pt solid black;text-align: center;"|<math display="inline">5.720{\times 10}^{10}</math><span style="text-align: center; font-size: 75%;"> Pa</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Linear weight in water'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1453 N/m</span>
|}


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Garcia-Espinosa_291836717-image13.jpeg|534px]] </div>

<div id="_Ref420661959" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Figure 10: Model basin setup including elastic mooring.</span></div>

====2.4 Validation test 1: monochromatic waves====

In this section, Response Amplitude Operators (RAOs) of the computational model presented in this work are compared with those obtained in the experiments. It should be noted that the experimental data shows a change in the platform response throughout the experiment, meaning, the results of the spectral analysis depend on the time interval used. If the period of time used for calculating the RAOs is chosen towards the end of the experiment, an increase in the response in the low frequencies is observed which raises the question of whether longer waves are suitable for analysis.

<span id='OLE_LINK32'></span><span id='OLE_LINK36'></span><span id='OLE_LINK38'></span>Experimental RAOs were obtained for full scale wave heights of <math display="inline">H=</math><math>2\, m</math> and <math display="inline">H=5\, m</math>, including the elastic lines used for the calibration. The experimental RAO was obtained using a time interval of ten wave cycles, starting at 500 seconds, and using a Fast Fourier Transform (FFT) to filter low frequency components induced by the model basin. The numerical RAOs were obtained using a time interval of ten wave cycles, starting at 100 seconds, and using 50 seconds of initialization. Figure 12 compares the RAOs in surge, heave, and pitch, respectively, against the numerical results obtained by the FEM and WADAM/SIMO solvers. While sufficient correlation is found for the lower periods, the correlation is not as good for the longer ones. This could be caused by the fact that the distance from the wave generator to the platform is 15.10 meters, while the wavelength ranges from 2.84 m to 26.34 m, and, longer waves may not have time to fully develop.

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">''' [[Image:Draft_Garcia-Espinosa_291836717-image14.png|360px]] '''</span> [[Image:Draft_Garcia-Espinosa_291836717-image15.png|360px]]  [[Image:Draft_Garcia-Espinosa_291836717-image16.png|360px]] 
|-
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Figure 11: Comparison between experimental and numerical decay tests.</span>
|}


<span id='OLE_LINK69'></span><span id='OLE_LINK70'></span><span id='OLE_LINK71'></span><span id='OLE_LINK75'></span>The experimental data shows that differences between the RAOs for the two wave heights can be large. We find two possible causes: due to the non-linearities introduced by the mooring system, and/or due to the experimental uncertainties.

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="text-align: center;vertical-align: top;width: 100%;"|[[Image:Draft_Garcia-Espinosa_291836717-image17.png|366px]]  [[Image:Draft_Garcia-Espinosa_291836717-image18.png|372px]]  [[Image:Draft_Garcia-Espinosa_291836717-image19.png|366px]] 
|-
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Figure 12: Comparison of RAOs obtained from experiments and numerical simulations.</span>
|}


====2.5 Validation test 2: bichromatic waves====

A number of tests were carried out with bichromatic waves in order to analyse the second-order response of the platform. The incident wave periods range between 5.5 and 21 seconds. The wave frequency difference ranges from 0.0128 Hz to 0.0167 Hz. The latter frequency is intentionally close to the surge resonant frequency as the focus of the analysis will be on the surge response to the slow drift.

<span id='cite-_Ref465934909'></span>[[#_Ref465934909|Table 12]] provides the experimental test matrix, including incident wave height ( <math display="inline">H</math>), frequency <math display="inline">(f)</math>, and the frequency difference and sum. All these cases have been simulated in the time-domain using the FEM model proposed in this work for the diffraction-radiation problem. The different cases have been analysed using the quasi-static elastic catenary and the non-linear FEM dynamic cable models. Once the simulations were carried out, the time series were transformed via FFT to the frequency domain in order to enable the comparison of results with the experimental data. No filtering was undertaken on the experimental and numerical data.

When analysing the wave elevation in the experiments, it was found that the wave energy spectrum was not bichromatic, showing energy spread around the frequency pair under analysis. In order to better reproduce the experiment, a set of waves reproducing the free surface elevation of the experiment was used, instead of a pure bichromatic wave.

Figure 13 compares the spectrum of the surge movement obtained in the experiments against those obtained numerically using the FEM solver. Overall, in the range of the incident wave frequencies, the numerical results correlate reasonably well with the experimental ones. However, in the low frequency range, larger differences are noted in some cases. Considering the difficulty of matching the RAOs in the monochromatic wave tests mentioned above, the correlation of the bichromatic wave test is acceptable.

When comparing the quasi-static catenary model with the dynamic cable model, it is observed that both models provide reasonably similar results, although the catenary model predicts larger movements.  This could be due to the lack of energy dissipation of the catenary model itself, while in the dynamic cable model Morison forces in the mooring lines are taken into account, leading to energy dissipation.

{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="vertical-align: top;width: 50%;"|[[Image:Draft_Garcia-Espinosa_291836717-image20.png|288px]]

<span style="text-align: center; font-size: 75%;">(a): Case 1</span>
|  style="vertical-align: top;width: 50%;"|[[Image:Draft_Garcia-Espinosa_291836717-image21.png|288px]]

<span style="text-align: center; font-size: 75%;">(b): Case 2</span>
|-
|  style="vertical-align: top;width: 50%;"|[[Image:Draft_Garcia-Espinosa_291836717-image22.png|288px]]

<span style="text-align: center; font-size: 75%;">(c): Case 3</span>
|  style="vertical-align: top;width: 50%;"|[[Image:Draft_Garcia-Espinosa_291836717-image23.png|288px]]

<span style="text-align: center; font-size: 75%;">(d): Case 4</span>
|-
|  style="vertical-align: top;width: 50%;"|[[Image:Draft_Garcia-Espinosa_291836717-image24.png|288px]]

<span style="text-align: center; font-size: 75%;">(e): Case 5</span>
|  style="vertical-align: top;width: 50%;"|[[Image:Draft_Garcia-Espinosa_291836717-image25.png|288px]]

<span style="text-align: center; font-size: 75%;">(f): Case 6</span>
|-
|  colspan='2'  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Figure 13 (a)-(p): Comparison between the experimental, catenary and cable model for surge response to bichromatic waves.</span>
|}


{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">''' [[Image:Draft_Garcia-Espinosa_291836717-image26.png|288px]] '''</span>

<span style="text-align: center; font-size: 75%;">(g): Case 7</span>
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">''' [[Image:Draft_Garcia-Espinosa_291836717-image27.png|288px]] '''</span>

<span style="text-align: center; font-size: 75%;">(h): Case 8</span>
|-
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">''' [[Image:Draft_Garcia-Espinosa_291836717-image28.png|288px]] '''</span>

<span style="text-align: center; font-size: 75%;">(i): Case 9</span>
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">''' [[Image:Draft_Garcia-Espinosa_291836717-image29.png|288px]] '''</span>

<span style="text-align: center; font-size: 75%;">(j): Case 10</span>
|-
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">''' [[Image:Draft_Garcia-Espinosa_291836717-image30.png|288px]] '''</span>

<span style="text-align: center; font-size: 75%;">(k): Case 11</span>
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">''' [[Image:Draft_Garcia-Espinosa_291836717-image31.png|288px]] '''</span>

<span style="text-align: center; font-size: 75%;">(l): Case 12</span>
|-
|  colspan='2'  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Figure 13: (continued)</span>
|}


<span id='cite-_Ref481163282'></span>[[#_Ref481163282|Table 13]] provides data regarding the CPU time required for simulating the bichromatic wave tests. On average, the ratio of CPU time to real time is in the order of 30.5.

{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">''' [[Image:Draft_Garcia-Espinosa_291836717-image32.png|288px]] '''</span>

<span style="text-align: center; font-size: 75%;">(m): Case 13</span>
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">''' [[Image:Draft_Garcia-Espinosa_291836717-image33.png|288px]] '''</span>

<span style="text-align: center; font-size: 75%;">(n): Case 14</span>
|-
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">''' [[Image:Draft_Garcia-Espinosa_291836717-image34.png|288px]] '''</span>

<span style="text-align: center; font-size: 75%;">(o): Case 15</span>
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">''' [[Image:Draft_Garcia-Espinosa_291836717-image35.png|288px]] '''</span>

<span style="text-align: center; font-size: 75%;">(p): Case 16</span>
|-
|  colspan='2'  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Figure 13: (continued)</span>
|}


<div id="_Ref465934909" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Table 12: Bichromatic test matrix.</span></div>

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  rowspan='2' style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Case'''</span>
|  colspan='2'  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Incident wave 1'''</span>
|  colspan='2'  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Incident wave 2'''</span>
|  colspan='2'  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Freq. difference'''</span>
|  colspan='2'  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Freq. Sum'''</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<math display="inline">{\boldsymbol{H}}_{\boldsymbol{1}}</math><span style="text-align: center; font-size: 75%;">''' '''[m]</span>
|  style="border: 1pt solid black;text-align: center;"|<math display="inline">{\boldsymbol{f}}_{\boldsymbol{1}}</math><span style="text-align: center; font-size: 75%;">''' '''[Hz]</span>
|  style="border: 1pt solid black;text-align: center;"|<math display="inline">{\boldsymbol{H}}_{\boldsymbol{2}}</math><span style="text-align: center; font-size: 75%;">''' '''[m]</span>
|  style="border: 1pt solid black;text-align: center;"|<math display="inline">{\boldsymbol{f}}_{\boldsymbol{2}}</math><span style="text-align: center; font-size: 75%;">''' '''[Hz]</span>
|  style="border: 1pt solid black;text-align: center;"|<math display="inline">{\boldsymbol{f}}_{\boldsymbol{2}}-</math><math>{\boldsymbol{f}}_{\boldsymbol{1}}</math><span style="text-align: center; font-size: 75%;">''' '''[Hz]</span>
|  style="border: 1pt solid black;text-align: center;"|<math display="inline">{\boldsymbol{T}}_{\boldsymbol{diff}}</math><span style="text-align: center; font-size: 75%;">''' '''[s]</span>
|  style="border: 1pt solid black;text-align: center;"|<math display="inline">{\boldsymbol{f}}_{\boldsymbol{1}}+</math><math>{\boldsymbol{f}}_{\boldsymbol{2}}</math><span style="text-align: center; font-size: 75%;">''' '''[Hz]</span>
|  style="border: 1pt solid black;text-align: center;"|<math display="inline">{\boldsymbol{T}}_{\boldsymbol{sum}}</math><span style="text-align: center; font-size: 75%;">''' '''[s]</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''1'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.63</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0582</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">4.32</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0735</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0153</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">65.4</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1318</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">7.59</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''2'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.27</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0667</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.54</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0813</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0146</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">68.3</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1480</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">6.76</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''3'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.80</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1053</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.62</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1220</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0167</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">59.9</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2272</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">4.40</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''4'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.13</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1205</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.27</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1351</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0147</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">68.2</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2556</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.91</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''5'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.88</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1300</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.14</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1429</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0128</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">78.0</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2729</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.66</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''6'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.50</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1449</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.92</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1587</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0138</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">72.4</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.3037</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.29</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''7'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.67</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1370</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.02</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1515</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0145</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">68.8</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2885</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.47</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''8'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.35</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1515</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.84</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1667</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0152</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">66.0</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.3182</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.14</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''9'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.22</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1613</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.77</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1754</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0141</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">70.7</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.3367</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.97</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''10'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.11</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1667</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.70</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1818</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0152</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">66.0</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.3485</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.87</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''11'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.83</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0909</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.10</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1053</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0144</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">69.7</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1962</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.10</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''12'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.37</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0833</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.44</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0997</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0163</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">61.2</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1830</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.46</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''13'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">4.59</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0714</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.16</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0850</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0135</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">73.8</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1564</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">6.39</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''14'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.64</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0526</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.16</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0672</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0145</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">68.9</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1198</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">8.35</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''15'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.82</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0526</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.16</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0672</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0145</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">68.9</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1198</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">8.35</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''16'''</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.44</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0476</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">6.03</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0625</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0149</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">67.2</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.1101</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">9.08</span>
|}


<span id='_Ref481163282'></span><div id="_Ref481163277" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Table 13: CPU time for bichromatic wave tests</span></div>

{| style="width: 89%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  rowspan='2' style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Case'''</span>
|  colspan='3'  style="border-top: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''Catenary mooring'''</span>
|  colspan='3'  style="border-top: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''Cable mooring'''</span>
|-
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Simulation Time '''[s]</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''CPU Time '''[s]</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Ratio'''</span>

<span style="text-align: center; font-size: 75%;">[s/s]</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Simulation Time '''[ST(s)]</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''CPU Time '''[s]</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Ratio'''</span>

<span style="text-align: center; font-size: 75%;">[s/s]</span>
|-
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''1'''</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">775.02</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">36227.18</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">46.74</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1200.03</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">72016.11</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">60.01</span>
|-
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''2'''</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1115.14</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">23947.44</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">21.47</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1115.14</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">19757.59</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">17.72</span>
|-
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''3'''</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1279.19</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">28688.95</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">22.43</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1279.19</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">32642.45</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">25.52</span>
|-
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''4'''</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1615.12</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">36035.17</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">22.31</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1615.12</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">41652.91</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">25.79</span>
|-
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''5'''</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1046.01</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">25104.81</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">24.00</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1046.01</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">27602.47</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">26.39</span>
|-
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''6'''</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">935.09</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">29797.70</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">31.87</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">935.09</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">24044.09</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">25.71</span>
|-
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''7'''</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">982.07</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">25378.13</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">25.84</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">982.07</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">25797.82</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">26.27</span>
|-
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''8'''</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">896.05</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">30768.33</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">34.34</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">896.05</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">30533.34</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">34.08</span>
|-
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''9'''</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">854.08</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">30032.72</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">35.16</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">854.08</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">23306.34</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">27.29</span>
|-
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''10'''</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">830.08</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">31300.66</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">37.71</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">830.08</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">38213.82</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">46.04</span>
|-
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''11'''</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">710.02</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">53935.90</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">37.98</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">710.02</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">37674.61</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">53.06</span>
|-
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''12'''</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1613.07</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">46106.76</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">28.58</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1613.07</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">39654.95</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">24.58</span>
|-
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''13'''</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1813.12</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">28062.74</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">15.48</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1813.12</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">24148.93</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">13.32</span>
|-
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''14'''</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1052.09</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">26884.85</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">25.55</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1052.09</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">25070.24</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">23.83</span>
|-
|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''15'''</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1052.09</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">41201.61</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">39.16</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1052.09</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">37945.54</span>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">36.07</span>
|-
|  style="border-left: 1pt solid black;border-right: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''16'''</span>
|  style="border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">836.06</span>
|  style="border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">29020.15</span>
|  style="border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">34.71</span>
|  style="border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">836.06</span>
|  style="border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">23880.50</span>
|  style="border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">28.56</span>
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''Mean value'''</span>
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">30.21</span>
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">30.89</span>
|}


====2.6 Validation test 3: irregular waves====

Two tests of irregular waves are used to validate the capability of the present FEM solver to predict the seakeeping of the HiPRWind, as well as the behaviour of the mooring. The simulation time was 784.7 seconds, and the time interval used for the analysis range from 500 s to 784.7 s. No filtering was used during the analysis process, to the experimental data or to the numerical results. <span id='cite-_Ref482262152'></span>[[#_Ref482262152|Table 14]] provides the particulars of the target wave spectrums. The model discretization used in this test is the same as for the bichromatic tests, but only the cable model is used to simulate the mooring lines.

{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">''' [[Image:Draft_Garcia-Espinosa_291836717-image36.png|534px]] '''</span>

[[Image:Draft_Garcia-Espinosa_291836717-image37.png|534px]] 
|-
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Figure 14: Comparison between experimental and second order numerical incident wave elevation for the time range analyzed. Up: Irregular 1; down: Irregular 2.</span>
|}


The incident wave field was determined by means of a FFT analysis of the incident wave elevation reported by the experiments within the analysis time interval. Figure 14 compares the resulting second order numerical incident wave elevation with the experimental one, finding good correlation between them.

{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="text-align: center;vertical-align: top;"|''' [[Image:Draft_Garcia-Espinosa_291836717-image38.png|534px]] '''

[[Image:Draft_Garcia-Espinosa_291836717-image39.png|534px]] 
|-
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Figure 15: Comparison between experimental and second order numerical surge movements for the time range analyzed. Up: Irregular 1; down: Irregular 2.</span>


|}


Figure 15 compares the second-order numerical and experimental surge response within the analysis time interval. Good correlation is found, although some deviation in the low frequency can be observed. Figure 16 compares the second-order numerical and experimental heave response within the analysis time interval. Again a fair agreement is found, although the numerical solution shows larger amplitudes in the higher frequencies. Figure 17 compares the second-order numerical and experimental pitch response within the analysis time interval. Just like in surge, a good correlation is found, although some deviation in the low frequency can be observed. Regarding the phase agreement, all movements show a good phase correlation between the numerical and experimental results.

{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="text-align: center;vertical-align: top;"|''' [[Image:Draft_Garcia-Espinosa_291836717-image40.png|534px]] '''

[[Image:Draft_Garcia-Espinosa_291836717-image41.png|534px]] 
|-
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Figure 16: Comparison between experimental and second order numerical heave movements for the time range analyzed. Up: Irregular 1; down: Irregular 2.</span>


|}


Figure 18 and 19 provide a spectral analysis of the loads induced by the mooring lines on the platform. The numerical results obtained follow the trend of the experiments well, although the amplitude at some frequencies do not match the experimental ones.

{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="text-align: center;vertical-align: top;"|''' [[Image:Draft_Garcia-Espinosa_291836717-image42.png|534px]] '''

[[Image:Draft_Garcia-Espinosa_291836717-image43.png|534px]] 
|-
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Figure 17: Comparison between experimental and second order numerical pitch movements for the time range analyzed. Up: Irregular 1; down: Irregular 2.</span>
|}


<div id="_Ref482262152" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Table 14: Irregular wave test matrix.</span></div>

{| style="margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">'''Case'''</span>
|  style="border: 1pt solid black;text-align: center;"|<math>{\mathit{\boldsymbol{H}}}_{\mathit{\boldsymbol{s}}}\mathit{\boldsymbol{\, }}[m]</math>
|  style="border: 1pt solid black;text-align: center;"|<math display="inline">{\mathit{\boldsymbol{T}}}_{\mathit{\boldsymbol{p}}}</math><span style="text-align: center; font-size: 75%;">''' '''[s]</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Irregular 1</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.50</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">16</span>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Irregular 2</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">4.00</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">13</span>
|}


====2.7 Viscous effects on mean wave forces and wave drift damping.====

<span id='_Ref488308824'></span>For semi-submersible platforms, when wave drift forces become small, viscous effects may contribute to drift forces. This effect is third order, meaning that it will become more important as the wave amplitudes increase (see [<span id='cite-31'></span>[[#31|31]]]). In the case of the HiPRWind, we can break down the contribution of the viscous effects to the mean force on two main components. The first due to the viscous forces generated at the heave plates in combination with the pitch movement, which provides a horizontal component with non-zero mean value. The second component comes from the variation of the free surface elevation at the columns waterline, which also leads to a non-zero mean horizontal force. The contribution of these forces has been estimated a posteriori, and it has been compared to the mean wave force for two case studies: bichromatic wave case 1, and irregular waves with <math display="inline">{H}_{s}=</math><math>2.5\, m</math> and <math display="inline">{T}_{p}=16\, s</math>. In both cases we found that the contribution of the viscous effects to the mean drift were in the order of magnitude of the wave mean drift. Therefore, this effect should be taken into account in order to evaluate the mean excursion of the platform. However, the non-zero frequency response should not be affected.

{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="text-align: center;vertical-align: top;"|''' [[Image:Draft_Garcia-Espinosa_291836717-image44.png|534px]] '''

[[Image:Draft_Garcia-Espinosa_291836717-image45.png|534px]] 
|-
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Figure 18: Comparison between experimental and second order numerical line 1 loads in the frequency domain. Up: Irregular 1; down: Irregular 2.</span>
|}


A posteriori estimation of the wave drift damping was calculated for all case studies analysed in this work. It was assumed that the major contribution to this damping is due to the columns. For this estimation, the damping of one column was approximated by the damping of a vertical circular cylinder reported in Fig 5.22 of [<span id='cite-_Ref488308824'></span>[[#_Ref488308824|32]]]. An estimation for the HiPRWind wave drift damping in bichromatic waves provides a value of about 7% of the surge viscous damping at the most, while being completely negligible for the irregular wave cases. Therefore the wave drift damping can be considered negligible compared to the viscous damping, and no significant change in the results should be expected if it was considered in the simulations.

{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="text-align: center;vertical-align: top;"|''' [[Image:Draft_Garcia-Espinosa_291836717-image46.png|534px]] '''

[[Image:Draft_Garcia-Espinosa_291836717-image47.png|600px]] 
|-
|  style="text-align: center;vertical-align: top;"|<span id='_GoBack'></span><span style="text-align: center; font-size: 75%;">Figure 19: Comparison between experimental and second order numerical line 1 loads in the frequency domain. Up: Irregular 1; down: Irregular 2.</span>


|}


===Conclusions===

A FEM model for the second-order wave diffraction-radiation problem in the time-domain has been developed. The model has been verified against analytical solutions, comparing mean drift loads for a floating hemisphere, second-order forces for a semi-submerged horizontal rectangular cylinder, and second-order forces induced on a bottom mounted cylinder. In all cases, the numerical results obtained correlate well with the analytical results.

Furthermore, the computational solver has been validated against experiments for the HiPRWind model carried out at the EC Nantes’ facilities. First, the model viscous damping was calibrated to reproduce decay tests for surge, heave, and pitch. A good match between the experiments and the calibrated model was obtained. Secondly, RAOs were compared using the monochromatic wave tests, the numerical results correlate well with the experiments for shorter waves, but not for longer waves, potentially due to the distance of the platform to the wave generator. Thirdly, the surge response of the HiPRWind subject to bichromatic waves were performed and compared to the experiments. Taking into account the experimental difficulties associated with measuring second-order quantities accurately, it has been found that the computed movements of the platform are in reasonable correlation with the experimental data. In the fourth and final stage, a validation of irregular waves under two different condition were carried out and the up to second-order movements amplitudes predicted numerically were compared with the experiments, showing good phase correlation, and reasonable correlation in amplitude.

With regards to the modelling of the mooring lines, results obtained for bichromatic waves correlate well with the experimental results. Both the quasi-static catenary model and the nonlinear FEM dynamic cable provide reasonably similar results. For irregular waves, the loads induced by the dynamic cable on the structure and the experimental results share a similar pattern for the frequency range analysed.

In conclusion, the second-order time-domain FEM model presented in this work, along with the mooring models used, are a good alternative for solving the second-order seakeeping of floating platforms.

===Acknowledgements===

The authors acknowledge EC Nantes whose facilities (used under the EU MARINET program) made this work possible.

The authors thank Acciona Energía and Fraunhofer Institute, and especially Raul Manzanas, for providing the data regarding the FP7 project HiPRWind, and the “Universidad Politécnica de Madrid” for granting access to EU MARINET SEMISO project experimental results.

Thanks to Carlos Lopez Pavon for providing the numerical results obtained using WADAM/SIMO and shown in this work.

The research leading to these results has received funding from the Spanish Ministry for Economy and Competitiveness under Grants ENE2014-59194-C2-1-R and ENE2014-59194-C2-2-R (X-SHEAKS).

The authors also fully acknowledge the support of NVIDIA Corporation for the donation of a TeslaK20 GPU used for this research.

===References===

<div id="1"></div>
[[#cite-1|[1]]] Hanson, T. D., Skaare, B., Yttervik, R., Nielsen, F. G. and Havmoller, O., 2011, “Comparison of Measured and Simulated Responses at the First Full Scale Floating Wind Turbine HYWIND,” Presented at EWEA 2011, European

Wind Energy Association, Brussels, Belgium.

<span id='OLE_LINK39'></span><div id="2"></div>
[[#cite-2|[2]]] Roddier, D., Cermelli, C., Aubault, A. and Weinstein, A. “Windfloat: A Floating Foundation for Offshore Wind Turbines,” J. Renewable Sustainable Energy 2010, 2(3), 033104. doi: [http://dx.doi.org/10.1063/1.3435339 http://dx.doi.org/10.1063/1.3435339]

<div id="3"></div>
[[#cite-3|[3]]] Maine International Consulting LLC, 2012, “Floating Offshore Wind Foundations: Industry Consortia and Projects in the United States, Europe and Japan,” Overview Report, Bremen ME.

<div id="4"></div>
[[#cite-4|[4]]] Breton S.P. and Moe H. Status, plans and technologies for offshore wind turbines in Europe and North America. Renew Energy 2009; 34:646-54.

<div id="5"></div>
[[#cite-5|[5]]] Strach-Sonsalla M. and Muskulus M. Prospects of Floating Wind Energy, International Journal of Offshore and Polar Engineering, 2016. 26:81-87.

<span id='OLE_LINK41'></span><div id="6"></div>
[[#cite-6|[6]]] Lupton, R.C. and Langley, R.S. Scaling of slow-drift motion with platform size and its importance for floating wind turbines. Renewable Energy 2017; 101: 1013-1020. Doi: [https://doi.org/10.1016/j.renene.2016.09.052 https://doi.org/10.1016/j.renene.2016.09.052]

<div id="7"></div>
[[#cite-7|[7]]] Lopez-Pavon, C., Watai, R. A., Ruggeri, F., Simos, A. N. and Souto-Iglesias, A. Influence of Wave Induced Second-Order Forces in Semisubmersible FOWT Mooring Design. Journal of Offshore Mechanics and Arctic Engineering 2015; 137 / 031602-1.

<div id="8"></div>
[[#cite-8|[8]]] Peiffer, A., Aubault, A. and Weinstein, J., 2011, “A Generic 5 MW Windfloat for Numerical Tool Validation & Comparison Against a Generic Spar,” ASME Paper No. OMAE2011-50278.

<div id="9"></div>
[[#cite-9|[9]]] Liu, Y., Li, S., Yi, Q. and Chen, D. Developments in semi-submersible floating foundations supporting wind turbines: A comprehensive review. Renewable and Sustainable Energy Reviews 2016; 60:433–449.

<div id="10"></div>
[[#cite-10|[10]]] Liua M., Xiaoa, L., Lua, H. and Shia J. Experimental investigation into the influences of pontoon and column configuration on vortex-induced motions of deep-draft semi-submersibles. Ocean Engineering 2016; 123: 262–277.

<span id='OLE_LINK42'></span><span id='OLE_LINK43'></span><div id="11"></div>
<span style="text-align: center; font-size: 75%;">[[#cite-11|[11]]] Philippe, M., Babarit, A. and Ferrant, P. Aero-Hydro-Elastic Simulation of a Semi-Submersible Floating Wind Turbine.</span> <span style="text-align: center; font-size: 75%;">J. Offshore Mech. Arct. Eng 136(2), Paper No: OMAE-13-1006; doi: 10.1115/1.4025031.</span>

<div id="12"></div>
[[#cite-12|[12]]] Coulling, J.A., Goupee, J.A., Robertson A.N. and Jonkman, J.M. Importance of Second-order difference wave diffraction forces in the validation of FAST semi-submersible floating wind turbine model.  ASME 2013 32nd International Conference on Ocean, Offshore and Arctic Engineering Nantes, France June 9-14, 2013.

<span id='OLE_LINK66'></span><span id='OLE_LINK67'></span><div id="13"></div>
[[#cite-13|[13]]] Bayati, I., Joknman, J., Robertson, A. and Platt, A. The effects of second-order hydrodynamics on a semisubmersible floating offshore wind turbine. Journal of Physics: Conference Series 524 (2014) 012094 doi:10.1088/1742-6596/524/1/012094.

<div id="14"></div>
[[#cite-14|[14]]] Matos, V.L.F, Simos, A.N. and Sphaier, S.H. Second-order resonant heave, roll and pitch motions of a deep-draft semi-submersible: Theoretical and experimental results. Ocean Engineering 2011. 38: 2227–2243.

<div id="15"></div>
[[#cite-15|[15]]] Simos, A. N., Ruggeri, F., Watai, R. A., Souto-Iglesias, A., Lopez-Pavon, C. Slow-drift of a floating wind turbine: An assessment of frequency-domain methods based on model tests, Renewable Energy, Volume 116, Part A, February 2018, Pages 133-154, ISSN 0960-1481, [https://doi.org/10.1016/j.renene.2017.09.059 https://doi.org/10.1016/j.renene.2017.09.059]

<div id="16"></div>
[[#cite-16|[16]]] Servan-Camas, B. and García-Espinosa, J. Accelerated 3D multi-body seakeeping simulations using unstructured finite elements. J Comput Phys 2013; 252:382e403.

<span id='OLE_LINK50'></span><span id='OLE_LINK51'></span><span id='OLE_LINK52'></span><div id="17"></div>
[[#cite-17|[17]]] Servan-Camas, B. A time-domain finite element method for seakeeping and wave resistance problems. School of Naval Architecture and Ocean Engineering, Technical University of Madrid; 2016 [Doctoral thesis]. [http://oa.upm.es/39794/1/BORJA_SERVAN_CAMAS.pdf http://oa.upm.es/39794/1/BORJA_SERVAN_CAMAS.pdf]

<span id='OLE_LINK44'></span><span id='OLE_LINK45'></span><span id='OLE_LINK48'></span><div id="18"></div>
[[#cite-18|[18]]] Gutiérrez-Romero, J.E., Serván-Camas, B., García-Espinosa, J. and Zamora-Parra, B. Non-linear dynamic analysis of the response of moored floating structures. Marine Structures 2016; 49:116-137.

<div id="19"></div>
[[#cite-19|[19]]] Compassis. SeaFEM Theory Manual. Retrieved from [http://www.compassis.com/soporte. www.compassis.com/soporte.] 2016.

<div id="20"></div>
[[#cite-20|[20]]] Jonkman, J.M. Dynamic modelling and loads analysis of an offshore floating wind turbine, Technical report NREL/TP-500-41958; November 2007.

<div id="21"></div>
[[#cite-21|[21]]] Fernandes, A. C. and Levy, L. A. P. Cálculo de esforços de onda de primeira e segunda ordem em navios e plataformas flutuantes através de integração azimutal analítica, Congresso da SOBENA (1990). Río de Janeiro, Rj, Brazil.

<div id="22"></div>
[[#cite-22|[22]]] Kudou, K. The drifting force acting on a three-dimensional body in waves, J.S.N.A. Japan 1977; Vol. 141.

<div id="23"></div>
[[#cite-23|[23]]] Pinkster, J. A. Low frequency second order wave exciting forces on floating structures. Wageningen: Netherlands Ship Model Basin 650 (1981).

<span id='OLE_LINK68'></span><div id="24"></div>
[[#cite-24|[24]]] Suliz, W. Diffraction of second-order surface waves by semisubmerged horizontal rectangular cylinder. J. Waterway, Port, Coastal, Ocean Eng. 1993; Vol 119, 160-171.

<div id="25"></div>
[[#cite-25|[25]]] Eatock Taylor, R. and Huang, J. B. Semi-analytical formulation for second-order diffraction by a vertical cylinder in bichromatic waves. Journal of Fluids and Structures 1997; 11, 465 – 484.

<div id="26"></div>
[[#cite-26|[26]]] Shao, Y.L. and Faltinsen, O.M. Numerical Study on the Second-Order Radiation/Diffraction of Floating Bodies with/without Forward Speed. 25th Workshop on Water Waves and Floating Bodies 2010, Harbin, China.

<div id="27"></div>
[[#cite-27|[27]]] Kim, M.H. and Yue, D.K.P. The complete second-order dif fraction solution for an axisymmetric body. Part 2 :

bichromatic incident waves and body motions. J. Fluid Mech 1990, 211: 557-593.

<div id="28"></div>
[[#cite-28|[28]]] Moubayed, W.I. and Williams, A.N. Second-Order Hydrodynamic Interactions in An Array of Vertical Cylinders in Bichromatic Waves. Journal of Fluids and Structures 1995, 9: 61-98.

<div id="29"></div>
[[#cite-29|[29]]] The HiPRWind Project. (2017, June 01). Retrieved from [http://hiprwind.eu/ http://hiprwind.eu/]

<div id="30"></div>
[[#cite-30|[30]]] Sesam user manual: Wadam, Det Norske Veritas, 2010. Retrieved from [https://projects.dnvgl.com/sesam/manuals/Wadam_UM.pdf https://projects.dnvgl.com/sesam/manuals/Wadam_UM.pdf]

<div id="31"></div>
[[#cite-31|[31]]] Faltinsen, O. M. Sea loads on ships and offshore structures. Cambridge University Press, 1990.