1. Introduction
Mathematical modelling is crucial for the design, analysis and optimisation of WECs. An essential development trajectory to an economically competitive WEC requires early device experimentation and refinement using numerical tools [
1]. Significant effort has been undertaken to develop and implement models describing the hydrodynamic interaction between the fluid and the WEC [
2]. An accurate wave-to-wire model for a WEC must also include the dynamics of the WEC power take-off (PTO) system [
3]. Equally, for classes of floating WECs, a suitable mathematical model of the WEC mooring system is required.
For floating WECs, a mooring system is required to keep the device on station. A mooring system is any type of cable, chain, rope, or tether that connects a WEC to an anchoring system fixed on the sea floor. Mooring costs represent a significant portion of total costs for a WEC installation [
4], and are therefore a driver of WEC economic competitive ability. Additionally, the mooring system influences the dynamic behaviour of the WEC, contributing nonlinear resistance and reactance to the WEC dynamics and, unlike free-floating structures, moored floating bodies experience enhanced coupling between various modes of motion. MMs are therefore important for the design, analysis and optimisation of the overall WEC performance, as well as the mooring system itself.
Mathematical analysis is particularly important for mooring systems because physical modelling of the moorings is problematic. As an example, during both tank testing and prototype trials for the Wave Dragon device reported by Friis-Madsen et al. [
5], the relatively shallow waters did not allow a proper catenary system to be tested. In the full scale tests performed by Vickers and Johanning [
6] in real sea conditions, the lack of controllability of the open sea experiments led to anomalies being observed that could not necessarily be quantified. One major obstacle in physical testing of mooring systems is due to the dynamic response of mooring lines being scale-dependent, whereby perfect dynamic similitude between full-scale prototypes and small-scale physical models is difficult to achieve [
7]. Difficulties with scaling effects in physical modelling of mooring systems, plus their inherent costs, therefore increase the relevance/importance of MMs for the analysis of full scale systems [
8,
9,
10,
11].
However, despite their importance, a relatively limited amount of analysis on mathematical modelling for WEC mooring systems is available in the literature. A review of WEC mooring systems is given by Harris et al. [
4], but a review of mathematical modelling methods for these systems does not exist. A few reviews of MMs from general offshore engineering fields, were performed at the end of the 20th century. The major contributions came from: Brown et al. [
12], who performed a comparative study on numerical calculation of mooring line dynamic loading by 15 different mooring solvers. Thomas [
13] details a comparison of the predictions of mooring line dynamics using different time integration schemes, indicating which scheme is to be preferred, highlighting the numerical problems which can be expected to affect the quality of the solution, and how best to avoid these numerical problems. Chatjigeorgiou and Mavrakos [
14] comparatively examine two different numerical schemes dealing with mooring dynamics.
However, mooring systems for WECs entail a number of key differences from general offshore applications, as outlined in
Section 2.2. A limited amount of analysis has been performed for the specific case of mooring systems for WECs. For example, the various dynamic analysis options available for WEC mooring system design are discussed by Thomsen et al. [
15]. A method and application for validating numerical models for WEC mooring systems is presented in Harnois et al. [
16]. Yang et al. [
17] compare four simulation procedures for the analysis of fatigue of WEC moorings, while Bhinder et al. [
18] compared the modelling of three different mooring nonlinear effects using three different software packages. Vissio et al. [
19] compare two different mooring models for use in their numerical simulations of the “ISWEC” WEC device, and state that “...One of the main challenges with this system is finding a suitable mooring line model, in terms of accuracy and the corresponding computational expense of the simulation”.
The objective of this paper is provide a review and analysis of mathematical modelling for WEC mooring systems to fill this gap. To fufil this objective, the paper aims to identify the different mathematical modelling methods for WEC mooring systems, review their usage in the literature relating to WEC analysis and discuss the strengths and weaknesses of the different modelling methods for the different required applications in WEC analysis. The paper is organised as follows:
Section 2 provides an overview of mooring systems in wave energy conversion, detailing the different types of mooring systems, the requirements of WEC mooring systems and a comparison with moorings in other offshore engineering fields.
Section 3 outlines the mathematical modelling of mooring systems, identifying the variables and parameters in the mooring system and the physical effects that the models aim to describe.
Section 4 then details the different model types used to analyse mooring systems, such as static, quasi-static and dynamic models. The different methods available to numerically implement the dynamic models are presented, looking at the different spatial and temporal discretisation options, linearisation and frequency domain techniques, and system identification of simpler parametric models. The type of coupling between the mooring and WEC models is then discussed.
Section 5 features the different software packages available for performing mathematical modelling of mooring systems, listing the different commercial, open-source and in-house software packages used within the literature reviewed in this paper.
Section 6 outlines the types applications that MMs are used for in WEC analysis. For each application a literature review is collated, providing examples of when MMs were utilized for the different applications, the type of WEC and mooring system investigated, and what type of model was employed.
2. Mooring Systems in Wave Energy Conversion
Just as there is diversity among the different types of WEC devices, there is similar diversity among the associated mooring systems and their requirements. An overview of different generic types of WECs and discussion of their mooring requirements is provided by [
4]. Other good overviews of WEC mooring systems include Johanning et al. [
20,
21] and Fitzgerald and Bergdahl [
22].
The relevance of a mooring system for a particular type of WEC depends on the operating principle of the WEC. Obviously, for nonfloating WECs, a mooring system is completely irrelevant. However, for floating WECs, the relevance of the mooring system may range from minimal where the mooring is only required for station-keeping, to highly relevant where the mooring is an integral part of the WEC’s PTO system. For some classes of directional floating WECs, such as terminator or attenuator devices, the mooring system must allow the WEC to weathervane into the predominant wave direction, while for axisymmetric point absorber type WECs the wave directionality is not important. Mooring systems have been classified into three categories, depending on their relevance to the WEC operation [
15,
23,
24,
25]; passive, active and reactive:
Passive mooring: when the only purpose of the mooring is station keeping and it does not influence the power extraction of the WEC from the waves.
Active mooring: when, in addition to providing station-keeping, the mooring system also has a significant influence on the dynamic response and power extraction of the WEC.
Reactive mooring: when the mooring system provides reaction forces for the WEC to extract power from the waves. In this case, the mooring is an integral part of the system, perhaps linking the floating part of a WEC to the PTO, and are therefore especially suited when the PTO exploits the relative movements between the WEC and the sea floor.
2.1. Requirements
A good mooring system must maintain the position of the WEC under extreme loading conditions and ensure its survivability, while also allowing efficient conversion of wave energy in operational conditions [
24]. In general, the requirements for WEC mooring systems can be summarised [
22,
26,
27] as follows:
Accommodate tides.
All components must have adequate strength, fatigue life and durability for the operational lifetime. Corrosion must be considered and controlled in the design, and abrasion due to bottom contact or contact with other lines must be minimised or avoided if possible.
Assure device self-orientation with incoming wave direction, if required.
Be easy to monitor and maintain.
Be economical.
Ensure device survivability.
Minimize environmental impacts on the seabed or native flora and fauna.
Not adversely affect the performance of the WEC device by interfering with the device motion used for energy capture.
Occupy as little space as possible on the seabed, in order to allow the devices to be installed close to each other in arrays.
Reduce/remove loads on the electrical cables and connections.
Restrain device excursion within its reserved area.
2.2. Comparison with Moorings in Other Offshore Engineering Fields
While there exist many similarities between moorings in wave energy and in other offshore engineering fields, such as shipping and offshore oil and gas, WEC mooring systems pose their own unique set of challenges. Karimirad et al. [
25] detail the applicability of different offshore mooring and foundation technologies for WECs. A comparison of station-keeping design for WECs compared to offshore oil and gas platforms is performed by Paredes et al. [
28], which concludes that a deeper analysis of mooring systems for WECs is required, compared to what is demanded for offshore oil and gas platforms, in order to properly assess the effect on the PTO, and survivability. Likewise, Vickers [
29] states that directly implementing the experience of the offshore oil and gas industry could be seen as an appropriate first step; however, due to the added complexity of the operation of WEC devices, the design approach for WEC moorings needs deeper consideration.
Design standards of WEC moorings compared to other offshore fields given by Johanning et al. [
20], who discusses how existing standards used in the offshore oil and gas industry might be adapted, and raises pertinent issues related to the station-keeping of WECs. Johanning et al. [
20] suggest that mooring design for WECs must consider new requirements, other than those that are captured in existing offshore standards such as DNV-OS-E301 [
30], a view also shared by [
22,
31,
32].
The key differences between WEC mooring systems compared to those from other offshore applications are outlined below:
Mooring systems for WECs differ from typical offshore structures, since the WEC is typically designed to resonate, rather than to avoid resonance [
33]. The objective is often to maximize the motion and induce resonant responses at the dominant wave frequency (WF). Therefore, WEC mooring systems need to be designed with regard to dynamics in the WF range [
34], whereas, in other offshore industries, the “static” design principles require that the mooring system’s resonance is far from any exciting frequencies, and the ideal system has no motion response to waves or wave groups [
10].
WECs are usually installed in relatively shallow water depths, compared to other floating offshore structures, where the tidal range is large compared to the water depth and may be more difficult to accommodate. Additionally, the tidal flow can be significant at these nearshore locations. Due to the relatively small size of WECs, and their being moored in relatively shallow waters, the effect of waves, tide and current can be of greater significance than for other floating offshore systems [
10,
35].
WECs absorb large quantities of energy from the waves, which affects the characteristics of the drift forces [
34]. As an example, while the horizontal drift forces will, in general, have the direction of the wave propagation, Hong et al. [
36] show that the Backward-Bent Duct Buoy (BBDB) WEC has horizontal mean drift forces opposite to the wave direction. Similarly, Retzler [
37] presents the results of an experimental investigation of the slow drift motions induced by second order wave drift forces on the Pelamis WEC. The drift forces were found to be mostly due to wave power absorption, and Retzler [
37] concludes that, because of the absorbed wave power, the Pelamis is subjected to considerably larger drift forces than a vessel of similar dimensions.
To maximise their economic viability, it is likely that WECs will need to be deployed in densely packed arrays, as developers attempt to exploit the maximum amount of energy for a given sea surface area. This will also help to reduce the length of expensive power cables connecting devices. The seabed footprint area and the allowable excursion of each device within the array will therefore be restricted. The challenge this presents to the mooring system is unique [
6,
11,
22]. Anchor sharing between multiple mooring lines in a WEC array, will cause cyclic loading on the anchor to alternate in many different directions, which is not well understood in offshore geotechnics, since typical floating offshore oil and gas platform facilities only have one mooring line connected to one anchor [
38].
WECs will have power cables/umbilicals connected which may affect the device dynamics [
39,
40,
41]. Although this is in general different to other offshore applications, the case of floating production, storage and offloading (FPSO) units for offshore oil and gas have risers, which have comparable effects and can be modelled using similar methods.
Finally, there is typically a difference in economics between wave energy conversion and other offshore industries. For example, the revenue generated by a WEC is orders of magnitude lower than that of an oil platform [
22].
Although this section has shown a number of differences between mooring systems for WECs and other offshore fields, many of the mathematical modelling techniques developed to analyse mooring systems in the offshore engineering fields are useful for wave energy applications, as discussed in
Section 3.
2.3. Mooring Components and Configurations
Figure 1 depicts that a mooring system consists of different components and materials, arranged in a certain configuration. A good review of the different mooring system types and classifications is given by Harnois [
35], and the different components and configurations are reviewed by Karimirad et al. [
25]. The report in [
26] establishes an industry knowledge base for existing anchoring and mooring techniques as applied to WECs and provides practical guidance on anchors and mooring techniques and sources of reliable design information necessary to the design of WECs. The report contains information on anchors, mooring lines and connecting hardware, and mooring configurations.
2.3.1. Components and Materials
A mooring system is comprised of mooring lines, anchors and connectors. Additionally, buoys and clump weights can be attached to the mooring lines. Different materials can be chosen for a mooring line: for example chains, fibre ropes or wire ropes. A comparison of different mooring line materials is given by Harnois [
35]. The choice of material will influence factors such as the inertia, elastic stiffness and damping of the mooring line. The role of an anchor is to keep a mooring line attached to a fixed point on the seabed. The main requirements for an anchor are to be able to resist high, horizontal and in some cases vertical, loads in a given seabed type (soft to hard), to be easy to install, and to be cost-effective [
35]. Different anchor designs are available, such as: dead weight, drag embedment anchor, pile anchor and plate anchor. Discussions on anchors for WEC systems are given in [
26,
42].
2.3.2. Configurations
Various types of mooring configurations are used for general floating structures, as depicted in
Figure 2. The different configurations are classified as a single point,
Figure 2a,c, or spread mooring,
Figure 2b,d, according to whether there is one or multiple lines connected to the WEC, respectively. Single point mooring allows the WEC to weathervane around the connection point; however, unlike a spread mooring, this arrangement has no redundancy, should a line fail.
Mooring configurations are also classified as taut,
Figure 2a,b, or slack/catenary,
Figure 2c,d, depending on whether the line is stretched tight under pretension, or hangs loose, respectively. A catenary line relies on the weight of chain links or clump weights to provide a horizontal restoring force to the WEC. A taut-line system is composed of wire or synthetic ropes and is normally highly pre-tensioned. As a result, taut-line systems are usually quite stiff, both in horizontal and vertical directions, and can significantly reduce vertical motions of the WEC, if connected to the sea bottom directly [
24]. Compared to a catenary configuration, a taut mooring system usually occupies a smaller footprint area on the seabed and the excursion of the floating structure is reduced.
If the influence of the mooring system on the vertical motion of the WEC needs to be minimised, then the top part of mooring lines, directly connected to the WEC, should remain as horizontal as possible. In such cases, buoys can be attached, as shown in
Figure 2e,f, where the buoy is attached to the sea bottom by a vertical line and connected to the WEC horizontally. If the compliance of the mooring system needs to be increased, then subsurface buoys and clump weights can be attached to the mooring line, as shown in the lazy-S configuration in
Figure 2g, to increase the geometric compliance of the mooring system (see
Section 3.5.4).
A study into the consideration of mooring cables configurations for WECs is given in Fitzgerald and Bergdahl [
22]. Thomsen et al. [
43] consider and compare mooring configurations for large floating, motion-independent WECs, based on traditions from other offshore fields. Harris et al. [
4] analyse conventional mooring systems in terms of suitability for wave energy application, describing and assessing different types of mooring configurations. The configurations rated by Harris et al. [
4] as highly suitable for wave energy applications are listed below.
Spread moorings: Catenary mooring, mutli-catenary mooring, and taut spread mooring.
Single point moorings: Catenary anchor leg mooring (CALM), and single anchor leg mooring (SALM).
3. Mathematical Modelling of Mooring Systems
Mathematical models for mooring systems have been developed for over 50 years, following the early work of Walton and Polachek [
44], Polachek et al. [
45], and Harleman and Shapiro [
46] in the 1960s. Generally speaking, MMs must be able to predict the hydrostatic equilibrium position and the response due to environmentally induced loadings of the moored system [
11]. The MM represents the relationship between the mooring system parameters, variables and environmental inputs. This section identifies the different parameters, variables and environmental inputs relevant to a WEC mooring system. Then, the important physical effects, forces, loads, etc acting on the mooring system are presented, and the mathematical relationships describing these effects detailed.
3.1. Co-Ordinate Systems
In modelling a WEC mooring system, there are three main coordinate frames of references to consider, as depicted in
Figure 3:
The global frame of reference ,
The WEC frame of reference , and
The curvilinear distance along the mooring line, s.
The global frame of reference and the WEC frame of reference are three dimensional, defining the two horizontal and one vertical cartesian co-ordinate. The origin for the WEC frame of reference, , is defined at the centre of mass of the WEC, and for the global frame of reference the origin, O, is defined at any convenient fixed location. The curvilinear distance is a one dimensional co-ordinate, measured along the mooring line, where corresponds to the bottom end of the mooring line connected to the anchor, and to the top end connected to the WEC (L is the length of the mooring line).
Additionally, the following vectors and angles are useful to define, as depicted in
Figure 3:
: is the unit vector tangential to the mooring line,
: is the unit vector normal to the mooring line,
: is the angle of the mooring line top end with the WEC at the connection point, and
: is the angle of the mooring line with the ground at the touchdown point on the sea floor.
: is a vector containing three angles:
- −
corresponds to the angle between x and due to the WEC pitch displacement.
- −
corresponds to the angle between y and due to the WEC roll displacement.
- −
corresponds to the angle between z and , due to the WEC yaw displacement.
3.2. Mooring System Parameters
There are a number of parameters that define the mooring system and whose values must be selected when designing the system. These parameters are listed in
Table 1. A number of the parameters relate to the topology of the mooring line geometry, as depicted in
Figure 4, such as the number of mooring lines, the position of the anchor(s) on the sea floor, the location along the line(s) where any buoys or weights are connected, the initial length of the mooring line(s), the location of the connection point(s) on the WEC body where the mooring line(s) are attached, and the diameter of the mooring line(s). The other parameters relate to the material properties of the mooring line, such as its mass density, elasticity, pre-tension and internal damping co-efficient.
3.3. Mooring System Variables
The mooring system variables are listed in
Table 2 and are explained, in greater detail, in the following subsections.
3.3.1. Length and Strain
Due to elasticity, the length of the mooring line can stretch in response to tensile forces arising from the motion of the WEC. The length of the mooring line can be defined as:
The strain of the mooring line,
, is the amount it has been stretched and is defined as the ratio of the instantaneous length,
L, to the original length,
:
3.3.2. Tension
The mooring line tension is typically the value of most interest calculated by a MM. The tension defines the force that the mooring lines impose on the WEC. Additionally, the tension values indicate the maximum and cyclic loads that the mooring lines are subjected to, and are therefore vital in mooring design and in extreme load and fatigue analysis.
The tension is related to the mooring line strain or deformation,
, via a constitutive tension-strain relation with the general form [
47]:
where the overdot represents the derivative with respect to time. For the case of a linear stress-strain relationship:
where
A and
E are the the cross-sectional area and Young’s modulus of the mooring line.
The tension is a vector,
, pointing in the direction of the tangent to the mooring line,
. The tension can be decomposed into its horizontal,
, and vertical,
, components. The tension at the top end connection point to the WEC is therefore decomposed into its horizontal and vertical components,
and
, respectively, to determine the effects of the mooring forces on the WEC motion:
For a taut mooring system, the tension is constant throughout the line. For a slack mooring system, the horizontal tension is constant throughout the line, whereas the vertical tension varies proportional to the effective weight per meter of the mooring line in water,
W, and the distance along the mooring line from the touchdown point:
The tension at any point along a slack mooring line can therefore be given by:
3.3.3. Position
The position of the mooring line in the global reference frame,
, is largely determined by the positions of its two end points, the anchor and the connection point on the WEC. The position of the top end of the mooring line,
, is determined by the position of the WEC’s centre of mass,
, and the mooring line connection point on the WEC,
:
For a taut mooring system, the profile of the mooring line is a straight line between the two end points. For a slack mooring line, the profile of a mooring line under static equilibrium, is reasonably well described by the catenary equations:
For a slack mooring line, the touchdown point is where the mooring line comes into contact with the sea floor and is locally horizontal. When the WEC drifts away from the anchor, the mooring line will become taut and the touchdown point will move towards the anchor as more line is lifted from the sea floor. Conversely, when the WEC drifts towards the anchor, the touchdown point will move away from the anchor, as more line is laid down on the sea floor. Because the touchdown point depends on the geometry of the sea floor, the local bathymetry may also be considered a parameter of the mooring system. The catenary equations for a mooring line with sea floor contact are described in Garza Rios et al. [
48] and Jonkman and Buhl [
49]. If the WEC is able to drift sufficiently close to the anchor, the mooring line might become fully slack, running vertically from the WEC to the sea floor, in which case the catenary curve is replaced by an “L” shape and the catenary equations cannot be used anymore.
3.3.4. Velocity and Acceleration
The velocity and acceleration of the mooring line are also important variables. Often, it is useful to decompose these variables into the components normal and tangential to the mooring line: , where the subscripts t and n represent the tangential and normal components, respectively.
3.4. Environmental Inputs
The environmental inputs which drive the mooring line dynamics are the waves, currents, wind and the variation in water depth due to the tide, as depicted in
Figure 5. External forces are directly applied to the mooring lines from the fluid, due to the action of the waves and currents. Additionally, the WEC applies a force to the end of the mooring line, due to the WEC motion, which itself is driven by the waves, currents, winds and the tidal varying water depth.
The exciting forces are a combination of three excitation modes, as described in [
50]:
Steady current, mean wind and mean wave drift forces are creating a mean force, or static loading;
Low Frequency (LF) forces, on the general range of 0–0.02 Hz, are induced by slowly varying wave drift forces, unsteady wind forces, and slowly varying tidal forces;
WF forces, on the general range of 0.03–0.3 Hz, are induced by first order wave forces.
Depending on the type of analysis performed, some or all of these inputs can be considered for the MM. For example, recommendations on design inputs for offshore position mooring systems are given in the offshore standard, DNV OS-E301, by Det Norske Veritas (DNV) [
30], and for the specific case of WEC mooring systems by Bergdahl [
32].
3.5. Physical Effects
This section details the different physical effects, present in the mooring system, that the MM must represent. These effects include: the coupling to the WEC, the buoyancy and weight of the cable and attached buoys and clump weights, the interaction with the sea floor, snap loads, and the mechanical impedance of the cable. The mechanical impedance of the cable includes the resistive forces due to damping and the reactive forces due to the inertia and compliance of the cable.
3.5.1. Connection to the WEC
The motion of the mooring line connection point on the WEC, due to the WEC motion, defines the boundary condition for the motion of the upper end of the mooring line. While the WEC motion is driven by environmental forces, it is also restrained by forces from the mooring line at the connection point(s). Therefore, the WEC and mooring systems form a coupled system requiring a coupled analysis. However, as shown in
Section 4.3.4, uncoupled analysis can also be performed to reduce the computational demands, but it gives approximate answers.
3.5.2. PTO Force
Certain types of WECs react their wave-induced motion against the sea floor or submerged damper plates. The PTO reaction force is transmitted across the spanning distance, between the WEC body and the sea floor/damper plate, by a taut mooring cable. The PTO force is therefore applied to the WEC via the tension in the mooring line at the connection point. Examples of analysis of these types of systems include [
52,
53,
54,
55]
3.5.3. Buoyancy and Weight
The net buoyancy force,
equals the difference between the weight of the mooring line acting vertically downwards and the hydrodstatic bouyancy force acting upwards:
where
g is the gravitational constant (9.81 m·s
),
D is mooring line diameter,
is the mass density of the water and
the mass density of the mooring line. The effective weight per metre of the mooring line in water, introduced in
Section 3.3.2, can therefore be expressed as:
Buoys and clump weights can be attached along the line, which will add positive or negative buoyancy to the attachment points. The introduction of hydrostatic stiffness, in the form of surface buoys or submerged buoys, can reduce the overall cable weight, reduce induced loads and reduce the influence of the cable on the motion response of the device, as shown by Fitzgerald and Bergdahl [
22]. The beneficial effects of buoys in reducing the mooring line dynamic tension is shown by Mavrakos et al. [
56], provided that proper selection of the buoy size, number and location is performed. Mavrakos and Chatjigeorgiou [
57] and Yuan et al. [
58] also study the effects and the numerical analysis of attached buoys and clump weights on mooring lines.
3.5.4. Compliance
The motion of the WEC must be accommodated by compliance or flexibility in the mooring system. Compliance is typically achieved elastically by stretching of the mooring lines, or geometrically, whereby any curvature or slack in the mooring lines is straightened. Taut mooring lines therefore provide most of their compliance elastically, for which dynamic tension is dominant, whereas slack moorings are able to provide large amounts of geometrical compliance through the straightening of the catenary curve, for which static tension is dominant. When the WEC moves from its equilibrium position, the elastic stretching and/or geometric deformation of the mooring lines will result in a restoring force applied to the WEC from the mooring lines.
Elastic stretching of the mooring lines results in tensions opposing the direction of motion, increasing in a way dependent of the line material stiffness [
35]. For small line extensions, the increasing tension can be represented by the linear stress/strain relation in Equation (
4). However, certain types of mooring lines may display a nonlinear stress-strain relationship, with material nonlinearity arising from a tension-dependent Young’s modulus of the mooring fibre, as is the case of the polyester mooring lines analysed by Tahar and Kim [
59]. The stretching of fibre mooring lines is investigated in [
60] which concludes that nonlinear springs would better represent the characteristics of some fibre ropes. Therefore, Equation (
4) can be modified to include a nonlinear relationship between the elastic stretching of the mooring line and the restoring force on the WEC. For example, the software package ANSYS AQWA [
61] uses a fifth order polynomial to model nonlinear mooring line stiffness. Alternatively, Tsukrov et al. [
62] present a method whereby the elastic modulus is updated at each time step, dependent on the current value of strain in the mooring line.
Slack moorings mainly produce static restoring forces, resulting from the weight of the line. The restoring force can be calculated from the catenary equations, by taking the derivative of the tension with respect to the displacement. Jain [
63] presents a simplified method of calculating the equivalent stiffnesses of the mooring cables, derived analytically, by using the basic catenary equation for the cable equilibrium. When the WEC displacement increases, tensions in the mooring lines opposing the direction of motion increase nonlinearly, until the mooring line is fully lifted form the sea floor and becomes taut [
11]. Accessory buoys or clump weights can be used to locally modify the weight or buoyancy of the mooring line and consequently the stiffness characteristics of the mooring system [
35]. The advantage of using an accessory buoy as a hydrostatic spring is shown by Fitzgerald and Bergdahl [
22].
Additionally, slack mooring lines may also present dynamic restoring forces in the case of high frequency (HF) WEC motion, whereby small displacements at the top end of the mooring line, can result in large lateral displacements at the midpoint of the mooring line, as depicted in
Figure 6. As the frequency increases, nonlinear drag (which increases with the square of frequency) provides very large lateral resistance, causing the mooring line to ’freeze’ in place and the WEC motion to be then accommodated elastically by the mooring line [
64]. The thesis by Gobat [
65] focuses on the dynamics of geometrically compliant moorings.
3.5.5. Inertia
The inertia of the mooring system consists of the mass,
m, of the mooring line, plus the added mass,
, of the surrounding fluid that the mooring line must also move. The inertia force can be calculated from Morison’s equation [
66]:
where
is the empirically derived added mass coefficient, and
the magnitude of the normal component of the mooring line acceleration.
3.5.6. Resistance
Mooring systems will have a resistive lossy effect, introducing damping on WEC motions, see Johanning et al. [
10,
67], which is important to consider for estimating the energy production of the WEC in a realistic manner [
24]. There are a number a damping mechanisms present in a mooring system, namely: hydrodynamic drag, vortex-induced vibrations (VIV), internal line damping and seabed interaction [
12]. Here, the mooring line damping due to drag, vortex induced vibrations and internal line damping are outlined. Damping effects, due to the interaction with the sea floor, are addressed in
Section 3.5.7.
Calculating mooring line damping has been of the focus of a number of studies, due to its important role in the LF motions of moored vessels and platforms. The thesis by Raaijmakers [
68] focuses entirely on this subject. Historically, predictions of the LF dynamic behaviour of moored offshore structures assumed that the mooring system only influences the restoring force characteristics. However, the work by Huse [
69] indicates that the mooring line provides up to 80% of the total damping for the surge motion of a floating offshore oil and gas platform. Huse [
70] had previously proposed that drag forces, due to the motion of the mooring lines through the water, were the main source of energy dissipation, the reason being that a small offset of the top end can result in large transverse displacements along the mooring line (
Figure 6), inducing considerable drag forces. A number of other authors, such as Liu [
71], Bauduin and Naciri [
72], and Lie et al. [
73], have built upon the early work of Huse [
69,
70], to improve the modelling of mooring line-induced damping.
Energy dissipation associated with catenary mooring line damping is studied by Johanning et al. [
10], through a series of experimental tests, which concluded that mooring line damping increases in a nonlinear manner, as the line becomes increasingly taut. This is due to the aforementioned effect of small fairlead motions causing large transverse motions along the line (
Figure 6), being more pronounced for taut cables than for slack cables. The experiments also show increasing levels of damping with increasing frequency, which becomes more relevant as the pretension increases. Webster [
8] shows that the contributions of drag damping and internal damping compete: drag damping is associated with motions transverse to the mooring lines, while internal damping is associated with motions along the mooring line.
Drag
The drag force,
, on the mooring line is proportional to the square of the mooring line velocity through the water. This relationship is expressed via Morison’s equation:
where
is the drag coefficient and the term
is the relative velocity between the water and the mooring line. However, depending on the level of modelling detail, sometimes the motion of the fluid is not taken into account, in which case
is simply the velocity of the line. Liu and Bergdahl [
74] present a thorough treatment of drag, formulating it as:
where the subscripts
n and
t denote the contributions normal and tangent to the mooring line line, respectively.
is the current velocity and
is the wave excited water particle velocity.
The drag coefficient of the cable is therefore a critical parameter for the calculating the drag force, but is often chosen with some uncertainty [
75]. For example, marine growth on the mooring line can change its drag properties. Recommended practice for selection of drag coefficients is reported in the offshore standard DNV-RP-F205 [
50].
Internal Line Damping
The internal line damping,
, caused by frictional forces within the line, is typically modelled as a linear dashpot, proportional to the strain velocity:
where
is the internal line damping parameter, whose value depends on the mooring line material and its construction. Weller et al. [
76,
77] document experiments on different mooring lines, providing datasets to facilitate the development of rope and mooring system simulation tools. The experiments involve sinusoidal loading of the lines, by anchoring one end and attaching the other end to an oscillating piston in a test rig, from which properties such as the damping rate and stiffness values can be determined. The damping rate is calculated via:
, where
is the energy dissipated over each load–unload cycle,
the angular frequency of the oscillation, and
X the amplitude of piston displacement.
Vortex Induced Vibrations
VIV refers to the dynamic loading from fluctuations in fluid pressure, caused by vortex formation in the wake of the mooring line. Analysis of the forces due to vortex shedding, and calculations of the induced loading, is complex and typically relies on computational fluid dynamics (CFD) or laboratory testing [
78]. However, the frequencies at which VIVs occur can be relatively easily determined, from the Strouhal number.
The Strouhal number relates the vortex shedding frequency,
, to the velocity of the ambient fluid flow,
u, and a characteristic length dimension of the body (in this case the diameter of the mooring line). The Strouhal number is defined as:
, and its value for a cylinder is between 0.17 to 0.20 over a wide range of flow velocities [
68]. The phenomenon of ’lock-in’ occurs when
becomes close to the natural frequency of vibration of the mooring line, causing a resonant response and increased levels of drag and internal line damping. Drag forces for wire lines, in particular, can be amplified by VIVs, whereas for chain the amplification is considered negligible [
12,
51]. An extensive investigation into VIV of cables in ocean currents is given by Vandiver [
79].
3.5.7. Interaction with the Sea Floor
A significant length of catenary mooring lines is generally in contact with the sea floor. The interaction of the mooring line with the sea floor can result in energy dissipation through friction due to horizontal motion, kinetic energy loss from vertical motion, as the mooring line lays down on the sea floor, and then suction and static friction, as the cable is lifted from the sea floor. A popular approach to model the sea floor interaction, proposed by Webster [
8], is to treat the sea floor as a very stiff distributed elastic support, as depicted in
Figure 7. The support provides no local vertical force when the mooring line is more than one radius off the sea floor, but as the mooring line becomes closer to the bottom, the vertical force increases linearly. The spring constant of the elastic support is selected so that, when the mooring line is on the sea floor, the vertical force exactly balances the weight, in water, of the grounded mooring line. The energy dissipation is then modelled as a linear damper. Other studies to adopt versions of this approach include [
47,
74,
80,
81,
82,
83,
84,
85,
86].
3.5.8. Snap Loads
Snap or snatch loads are impulse-type mooring loads, with an amplitude much larger than static and dynamic mooring loads. Snap loads are caused when a mooring line becomes slack and is then suddenly re-tensioned, and are essentially discontinuous shock waves. For wave energy applications, the ratio between the horizontal motions of the floater and the water depth is usually large, thus increasing the probability of snap loads [
34]. Additionally, PTO end-stop collisions for tautly moored WECs may also cause snap loads, as experimentally investigated by Hann et al. [
87]. A good overview of snap loads in mooring systems for wave energy converters is detailed by Palm et al. [
86]. The numerical mooring line model must have adequate spatial and temporal resolution to accurately capture snap loads, with appropriate numerical methods for this purpose reported by Vassalos and Huang [
88] and Palm et al. [
86].
3.5.9. Bending and Torsional Effects
The rotational degrees of freedom of the mooring line may also be subjected to damping and restoring forces, as the lines bend and twist. Bending and torsional effects are often considered for pipe-style moorings and risers found in offshore oil and gas platforms. However, for typical mooring line materials and configurations in WEC mooring systems, the bending and torsional effects are not significant and can be neglected for offshore renewable energy applications [
89]. A sensitivity study was performed by Hall et al. [
90], revealing that bending and torsional elasticity were of negligible importance to the mooring line tension responses for the three different stability classes of floating wind turbines investigated, regardless of load case. The results were the same whether the model used the significant bending stiffness of a wire rope or the nonexistent bending stiffness of a chain. Including bending and torsional effects into WEC MMs would therefore unnecessarily increase the model complexity and computational requirements, for negligible improvement in the accuracy of the results.
3.6. Relevant Effects for Different Mooring Systems
Table 3 compares the relevant effects for taut and slack mooring systems.
4. Model Types
Mathematically describing the behaviour of the mooring system with relationships between the variables, parameters and environmental inputs, can be implemented with different model types, of varying degrees of complexity. These different model types can be classified as static, quasi-static or dynamic.
4.1. Static Models
A static model considers constant loads only, such as gravity, buoyancy, non-time varying current and wind, and mean wave-drift forces. Static analysis determines the equilibrium between the constant or mean environmental loads and the restoring force of the mooring lines on the WEC. Static equilibrium is determined in a series of iterative steps [
91]:
At the start of the calculation, the initial positions of the WEC, and attached buoys and clump weights, are defined; these, in turn, define the initial positions of the ends of any lines connected to them.
The equilibrium configuration for each line is then calculated, with the line ends fixed.
The out-of-balance loads acting on each free body (WEC, buoy, weight, node, etc.) are then calculated and a new position for the bodies is estimated.
The process, steps 1–3, is repeated until the out of balance load on each free body is zero (up to the specified tolerance).
The results of the static analysis provide the equilibrium position of the WEC. The mooring system parameters, which can determine the equilibrium position, are [
20]: (i) the mooring line pre-tension; (ii) the number of mooring lines, their orientation and material properties; and (iii) the length and dimension of these lines. By displacing the WEC through prescribed distances and calculating the static restoring load, the static analysis also provides relationships between the WEC displacement and the static restoring force of the mooring line, as depicted in
Figure 8.
For the majority of systems, the static analysis process is very quick and reliable. As an example, static solutions for multi-segment mooring lines with attached buoys and weights are derived by Oppenheim and Wilson [
92] and Smith and MacFarlane [
93]. A static design is often used at the initial stage for conceptual design purposes, to identify initial components and configurations for further analysis. Additionally, the static analysis is often used to provide a starting configuration for dynamic simulations [
91].
4.2. Quasi-Static Models
The quasi-static approach assumes that the motion of the system is uniform and linear between two static positions, during a given time step for which the loads on the systems are assumed constant [
29]. The dynamic effects on the mooring system are ignored by this method, omitting the motion dependency of mass, damping and fluid acceleration on the system [
20]. Quasi-static models derive the mooring line shape and tension from the catenary formulations, based on the assumption that the line is in static equilibrium in each time step, that inertia effects can be neglected and that the line profile is reasonably well described by the catenary equations [
72]. The mooring force on the WEC therefore consists only of the position dependent static restoring force, shown for example in
Figure 8. Implementation of quasi-static models for offshore mooring applications, and a review of their development, can be found in [
72,
84,
94,
95] and, for the specific case of wave energy applications, in Nava et al. [
96].
The main disadvantage of quasi-static models is that they neglect hydrodynamic and inertial forces on the line, which can affect the WEC response and are especially important for predicting the mooring loads [
89]. Additionally, for point absorber type WECs, which respond quickly due to WF motions, a situation similar to that illustrated in
Figure 9 could occur [
29], where the top of the line follows the WEC motion but further down the line there is a delayed response, violating the catenary profile assumptions of a quasi-static model. An illustration of the large possible discrepancies, between the predicted outputs from a quasi-static and a dynamic MM, can be seen in the various graphs of calculated displacements and tensions in Yang et al. [
97].
4.3. Dynamic Models
Like most mechanical systems, the dynamics of the mooring line are formulated using Newton’s second law (of motion), where the resulting motions are due to the forces acting on the body. The mooring line dynamics are governed by the following equation [
98]:
where
is a vector comprising the external forces acting on the line. Solutions to Equation (
19) can be difficult to obtain, except under simplified specific cases. An example of such a simplified case is given by Aranha and Pinto [
99], where an algebraic approximation is derived for an expression of the dynamic tension under sinusoidal excitation. Otherwise, in general, a solution to these equations must be obtained either by linearising the model or by employing numerical techniques to approximate the governing nonlinear differential equations [
64]. The numerical approximations are obtained via a spatial discretisation of the mooring line to form a set of ordinary differential equations (ODEs), which are subsequently discretised in time and solved by an integration algorithm. A dynamic simulation will often use a static analysis solution as its initial configuration and then evolve forward in time from there.
The following subsections detail the different spatial and temporal discretisation methods commonly used for nonlinear dynamic MMs. A linearisation approach is then outlined, which allows a frequency-domain representation of the mooring system. Obtaining computationally efficient dynamic mooring system models, via system identification techniques, is then presented. Lastly, the coupling between the dynamic models for the WEC and mooring system is discussed.
4.3.1. Discretisation
Walton and Polachek [
44] published the first treatment of discretised mooring line modelling by formulating the equations of motion for discrete elements, with hydrodynamic drag and added mass forces applied uniformly over each element. Centered finite differences were used to discretise the time derivatives and evolve the solution forward in time. Polachek et al. [
45] later added cable extensibility, thereby offering a thorough method for treating the nonlinear mooring line dynamics. Their discretisation method would nowadays be categorized as a lumped mass method, whereby the mass and externally applied forces are lumped at a discrete number of points/nodes, which are connected by massless springs [
65]. The lumped mass method and other spatial discretisation approaches, such as finite differences and finite element methods (FEMs), are outlined in this subsection. A comparison of these different spatial discretisation approaches is given by Mascoila et al. [
100], who explore the features of the various models and weigh the advantages of each.
Lumped Mass
The lumped mass method involves lumping of all the effects of mass, external forces and internal reactions at a finite number of nodes,
N, along the line, as depicted in
Figure 10. This procedure implies that the behaviour of a continuous line is modelled as a set of concentrated masses connected by massless springs [
101]. By applying the equations of dynamic equilibrium and continuity (stress/strain) to each mass, a set of discrete equations of motion is derived. The lumped mass method is shown to have the advantage of a strictly diagonal mass matrix, eliminating the requirement of added computation for matrix inversion [
100]. Detailed descriptions of example usage of lumped mass MMs can be found in [
81,
89,
101,
102,
103,
104] and, for the specific case of WEC moorings, in Vissio et al. [
19].
Finite Differences
The finite difference approach differs from lumped mass methods by using an infinitesimally small differential element, rather than a finite discrete element. Spatial derivatives may then be treated algebraically as:
using forward, backward or centred finite differences, respectively. Detailed descriptions of implementing finite differences for mooring line, cable and/or riser models are given in [
47,
105,
106,
107] and, for the specific case of WEC moorings, in Cerveira et al. [
85].
Finite Element
The FEM also discretises the mooring line into small differential elements; however, whereas the finite differences method uses the differential form of the governing equations, the FEM uses their integral form counterpart. Variables internal to each element are then described via basis/interpolation functions and the integrals are solved using numerical integration techniques. Therefore, unlike the finite difference method, which estimates the spatial gradients in the governing partial differential equation (PDE) with first-order finite difference functions, the FEM eliminates the spatial gradients by approximating the PDE locally by a set of ordinary differential equations (ODEs), similar to the lumped mass method. Whilst lumped mass models must place all masses at discrete nodes, FEMs can derive the governing equations with an integration of the mass over the entire element, leading to a “consistent” mass formulation [
80]. Although the resulting mass matrix achieves an upper and lower bandwidth of two, compared to the more computationally efficient strictly diagonal matrix for the lumped mass formulation [
100], the FEM is able to produce accurate results with a smaller number of elements compared to the lumped mass formulation. Examples of different solution methods include: collocation [
14], Galerkin [
88,
108,
109], and discontinous Galerkin [
110] methods. A detailed description of finite element modelling for mooring lines is given in [
108,
109,
110,
111] and, specifically for the case of mooring systems for WECs, in Palm [
112].
Other Spatial Discretisations
Whilst the lumped mass, finite differences and FEMs are the most established and popular discretisation schemes, a number of other methods have been used for mooring lines and are detailed here. Winget and Huston [
113] and Kamman and Huston [
114] utilise a finite segment scheme, consisting of a series of ball-and-socket connected rigid rods, to discretise a mooring cable. Nichol et al. [
115] implement a finite segment model for the mooring system of floating tidal energy devices. Garrett [
116] utilises a similar scheme consisting of elastic rods. Filipich and Rosales [
117] model a chain mooring line, by discretising the chain into its individual links.
Temporal Discretisation
The spatial discretisation of a mooring line produces a series of ODEs, which are then solved and advanced in time by an integration algorithm. Most temporal integration schemes are a member of the Newmark family [
118], which is characterized by the following equations [
119]:
where
is the time step;
,
and
are stepwise approximations to the unknown nodal values of the displacement, velocity and acceleration, respectively, at time
;
and
are adjustable parameters defining each particular member of the Newmark family. Implicit algorithms, such as the trapezoidal rule (
), are better for inertial problems, where the dynamic response is dominated by mode shapes with higher period values. Explicit algorithms such as the central difference method (
) are more adequate for transient problems [
119]. Explicit methods calculate the state of a system at time
from the state of the system at the time
, while implicit methods use the state of the system at
and
. The explicit scheme is robust, but may be slow because it requires a very short time step, while the implicit scheme is faster but may produce inaccurate results [
35].
The history, development and usage of the different numerical schemes used for the temporal discretisation and integration of mooring line dynamics is outlined in Gobat [
65] and the references therein. Thomas [
13] performs an investigation of time integration schemes for mooring line dynamics and gives a thorough review of the different available schemes and their derivations. Hearn et al. [
120] examines the influence of using different time integration schemes to solve the dynamic equations of motion applicable to a mooring line.
Finite-difference models discretise the ODEs in both space and time, thereby formulating an inherent time-stepping approach. Lumped-mass models, however, only discretise the ODEs in space, and therefore leave continuous time derivatives, requiring an external integration scheme to advance the model in time [
89].
The most popular finite difference scheme is the box method, depicted in
Figure 11, in which the governing equations are discretised on the half-grid point, in both space and time [
121]. By applying the discretisation on the half-grid points, with spatial and temporal averaging of adjacent grid points, the box method is second-order accurate [
121]. Cable dynamics applications that do not employ the box method, often use temporal integration, such as backward difference, which is only first-order accurate, or the generalized trapezoidal rule, which is second-order accurate only in its most unstable form [
121]. The box method is unconditionally stable for linear problems, but has no numerical dissipation and is susceptible to Crank-Nicholson noise [
121]. Therefore, Gobat [
65] demonstrates that the box method is seldom the best choice of temporal discretisation scheme for the cable dynamics problem. Koh et al. [
122] also come to this same conclusion and propose a modified box method that uses backward difference for the temporal discretisation. However, Gobat [
65] combines the box method spatial discretisation with the generalized-
method for temporal integration. Compared to other schemes, commonly employed for the temporal integration of the cable dynamics equations, including the box method, trapezoidal rule, backward difference, and Newmark’s method, the generalized-
algorithm has the advantages of second-order accuracy, controllable numerical dissipation, and improved stability, when applied to the nonlinear problem.
4.3.2. Linearisation and Frequency Domain Models
Linear models are computationally simpler than their nonlinear counterparts and offer the advantage of frequency domain analysis. Therefore, it is often desirable to linearise the governing dynamics to simplify the analysis and reduce the computational requirements of the model. For example, power production estimates for a WEC require simulating the WEC operation across a wide range of sea states, which can be computationally prohibitive in the time domain and is therefore typically calculated in the frequency domain. For example, Cerveira et al. [
85] argue that during power production conditions, the WEC motions are of relatively small amplitude, therefore it is possible to approximate the nonlinear mooring characteristics to the linear order. Cerveira et al. [
85] model the mooring system by adding linear mass, damping and stiffness coefficients to the linear equations of motion for the WEC, which are then solved in the frequency domain.
A frequency domain methodology is presented by Fitzgerald and Bergdahl [
123] to evaluate the effect of moorings on WEC power extraction. The nonlinear mooring dynamics are calculated numerically in the time domain, for different conditions and then approximated to the linear order. The approximation is considered to be valid for small motion amplitudes and load induced tensions that are small relative to the equilibrium pretension. A similar approach is taken by Spanos et al. [
55,
124,
125], who use a statistical linearisation scheme. A time-domain analysis of the nonlinear behavior of the system is performed, and then the mean-squared error between the nonlinear system and the linear one is minimised to determine, iteratively, the effective linear stiffness and damping of the mooring system.
Other examples of linearised MMs for WECs exist. Bachynski et al. [
33,
126] perform a simple frequency domain analysis on a point absorber WEC, model the mooring as additional linear restoring force and then add the mass of the mooring lines to the mass of the WEC mass. In Bachynski et al. [
127] and Zaroudi et al. [
128], linearised stiffness coefficients for the mooring system are computed using the method from Jain [
63], which allows fast calculation of the stiffness coefficients of catenary lines, for small floater motions, as a function of the pretension, cable weight in water per-unit arc length, fairlead and anchor locations, and line angle at the seabed. Vincente et al. [
129] use a frequency domain analysis to investigate the mooring effects on an array of WECs, modelling the interbody and bottom mooring lines as linear springs/dampers.
For general offshore applications, established techniques of statistical linearisation have been employed to linearise the nonlinearities in the restoring forces and drag, based on minimizing the mean-squared error and assuming Gaussian statistics [
130]. A treatment of the current and seabed friction in the frequency domain, obtained via statistical linearisation, is presented by Liu and Bergdahl [
74]. In Sarkar and Taylor [
131,
132], the nonlinear drag is linearised by preserving the energy dissipation over each cycle. Similarly, Lie et al. [
73] identify mooring line damping coefficients from finite element time domain simulations, using a method based on energy dissipation. Other papers which include linearised, or frequency domain treatment, of the dynamic mooring system model include [
57,
81,
106,
107,
116,
133,
134,
135,
136,
137,
138,
139].
However, although linear models and frequency domain techniques offer computationally efficient MMs, the necessary linearisation of nonlinear terms can induce errors in final predictions [
9]. As shown in [
140,
141,
142,
143,
144] , complex nonlinear responses are found in mooring systems with geometric nonlinearities and nonlinear drag, with multiple coexisting solutions and existance of chaotic explosions, which could not be predicted by equivalent linearisation methods. Additionally, Van den Boom [
101] shows that dynamic tension amplification is strongly influenced by nonlinearities, due to catenary effects, elasticity and drag, and therefore would not be captured in frequency domain models. The MM comparative study by Brown et al. [
12] reveals that line tension and damping results, using nonlinear time domain methods, were generally in good agreement with each other. In contrast, results generated using frequency domain techniques, and particularly damping results, showed large scatter and significant differences, in comparison to their nonlinear time domain counterparts.
4.3.3. System Identification
While the linear models in
Section 4.3.2 overcome the computational demands and complexity of the nonlinear, spatially disretised, time domain simulations in
Section 4.3.1, the accuracy of the linear models is limited and is unacceptable for many situations. Another modelling approach, described in this section, is that of system identification, where computationally efficient nonlinear parametric models are identified from input/output data of the system behaviour obtained from high-fidelity numerical simulations. The computationally expensive high-fidelity numerical simulations therefore only need to be performed once, and then the identified nonlinear parametric models can reproduce the system behaviour in future simulations at a fraction of the computational time. This approach has been used, for example, for WEC hydrodynamics, by using computationally expensive CFD experiments to produce realistic data of the WEC behaviour, which is then used to identify nonlinear parametric hydrodynamic models [
145].
System identification has been used to obtain computationally efficient MMs for general offshore mooring applications. For example: Pascoal et al. [
146,
147] represent the mooring force by a Taylor series and identify the parameters using data from a lumped mass time domain simulation. Umar and Datta [
148] and Murhpy et al. [
149] represent nonlinear restoring force terms with polynomials, whose coefficeints are identified by fitting the polynomials to the restoring force versus displacement curves obtained from numerical calculations. The dynamics of the mooring system have also been replicated by Volterra models [
150,
151] and neural networks [
152,
153,
154,
155,
156,
157].
Physical experiments have also been used to provide the identification data, instead of high fidelity numerical simulations. For example, Narayanan et al. [
158,
159] and Raman et al. [
160,
161] identify the polynomial coefficients for nonlinear restoring force curves from physical experiments. Gobat and Grosenbaugh [
162] develop a frequency domain model for dynamic tension identified from experimental data. Other examples of system identification used on experimental data to identify MMs include Grosenbaugh et al. [
163] and Gottlieb et al. [
164]. The work by Weller et al. [
76] aims to provide datasets, on internal line damping properties obtained from physical experiments, to aid in the identification and development of MMs.
An example of the system identification approach, applied to WEC moorings, is presented by Cerveira et al. [
85], where a reduction in model complexity for the mooring system dynamics was required because using the fully nonlinear mooring dynamics code was prohibitively time consuming, and prone to collapse, due to numerical instabilities. The mooring line force on the WEC was approximated by a Taylor series, where the coefficients were identified from data produced by simulations from the fully nonlinear mooring dynamics code.
4.3.4. Coupled or Uncoupled Models
A WEC and its mooring lines are intrinsically a coupled system; the dynamics of the WEC affect the mooring dynamics and vice versa. This requires a dynamic analysis method that can iterate the solution for the WEC and the mooring lines until they converge on a compatible solution [
29]. However, to save on computation, uncoupled analyses were traditionally performed, by solving the WEC and mooring dynamics separately. An early implementation of the coupled approach was performed by Ormberg and Larsen [
165], where motions and mooring line tensions from physical model tests and simulations, using coupled and uncoupled procedures, were compared. The main conclusion was that the traditional uncoupled approach may be severely inaccurate, especially for floating structures operating in deep waters.
The uncoupled analysis is usually undertaken in two steps. The first step involves the computation of the WEC motions, either without any influence from the moooring system included or only with quasi-static restoring forces from the moorings. The second step then involves applying the precalculated WEC motions, as boundary conditions at the top end of the mooring, to perform a dynamic analysis of the mooring system. The uncoupled approach reduces the complexity of the numerical simulation and often offers greater numerical stability, compared to the coupled approach [
166].
In the coupled approach, the equations of motion of the WEC and mooring lines are solved simultaneously. Gao and Moan [
24] prescribe the following procedure for a coupled mooring analysis: a static equilibrium is first found by a nonlinear static analysis, considering only the mean environmental forces. Mean offset and mean mooring line tension are then obtained. Next, the motion of the WEC and the mooring line tension are solved simultaneously, under the dynamic wave and wind forces, using an appropriate time integration scheme. However, two types of coupling schemes can be performed; weak coupling and strong coupling.
Weak coupling: At each time step, for the equations of motions of the WEC, one (or more) steps of a nonlinear dynamic analysis of the mooring lines is performed, using the WEC displacement from the previous time step at the fairlead of each line. The resulting forces and moments, at the top of the lines, are then included in the forces acting on the WEC, at the current time step.
Strong coupling: At each time step, the dynamic analysis of the mooring lines is performed, where the WEC is considered a “node” of this model, alongside the other discretised nodes of the mooring lines.
Jacob et al. [
167] present results of case studies employing the weak and strong coupling formulations, and the behavior of these formulations is compared, in terms of their accuracy, computational performance, and applicability to different types of offshore systems. Rodrigues et al. [
119] show that the weak coupling scheme provides very good results, with considerably better efficiency than the strong coupling scheme.
Other examples of coupled analyses include Tahar and Kim [
168], where the dynamics of hull, mooring lines, and risers were solved simultaneously in a combined matrix at each time step. Idris [
169] performs coupling between a buoy and cable by adopting the buoy equations of motion as boundary conditions, at one end, for the mooring problem. Garrett [
116] details a coupling scheme for moorings and risers for use with offshore floating platforms.
The coupling of the WEC and the mooring lines, discussed in this section, is referred to as Type 1 coupling by Low and Langley [
130]. Another type of coupling to be considered in the mathematical modelling of mooring systems, termed Type 2 coupling by Low and Langley [
130], is the coupling between the system response at different time scales, namely the mean offset, the response at the incident WF, and the LF second order wave drift motions. For example, the dynamic response of the mooring system at the WF is influenced by the catenary profile of the mooring line which changes with the LF drift motion, and the damping of the LF drift motions is influenced by the WF motions. Type 2 coupling is discussed in further detail in
Section 7.3.
5. Software Packages
There exist numerous software packages that implement MMs for offshore applications. The various software packages used within the publications reviewed in this paper are listed in
Table 4, and are grouped into commercial, open-source and in-house packages. The type of models available in each software package are indicated in the table, along with the references reporting on the targetted usage of each software package. Some software packages are able to simulate the WEC hydrodynamics as well as the MMs, which is also indicated in the table. However, Rhinefrank et al. [
170] note that, due to the specific nature of ocean energy harvesting devices, off-the-shelf software packages often need to be supplemented with custom made numerical approaches, an example of which would be a PTO model.
The commercial packages include AQWA, DNV SESAM, Flexcom, OrcaFlex and Proteus-DS. AQWA was developed by ANSYS and provides a toolset for investigating the effects of environmental loads on floating and fixed offshore structures. SESAM is a software suite for structural and hydrodynamic analysis of ships and offshore structures, originally developed by the Norwegian University of Science and Technology and later bought by Det Norske Veritas (DNV). SESAM consists of several modules, such as Deep C, MIMOSA, RIFLEX, SIMA and SIMO. The SESAM modules are compared by Chrolenko [
198] for the dynamic analysis and design of mooring lines. Flexcom is a structural analysis software package, developed by MCS Kenny, whose engineers are immersed in the offshore oil and gas industry. Flexcom offers a broad range of capabilities for modelling and analysing vessel and mooring behaviour due to relevant environmental loading. OrcaFlex is a software package, developed by Orcina, for the dynamic analysis of offshore marine systems. OrcaFlex has the capability to simulate the WEC and mooring coupled dynamics using both nonlinear time domain and linear frequency domain analysis, or can be used as a library, allowing integration into third party software and automation possibilities.
Table 4 shows that OrcaFlex is by far the most popular software package used in the publications reviewed in this paper. Proteus-DS is a time-domain simulation tool, developed by Dynamic Systems Analysis Ltd, used to analyse and design a variety marine, offshore and subsea systems and technologies using a variety of hydrodynamic and finite-element analysis techniques, and has been used to simulate WEC dynamics.
Two open-source software packages were identified, MoorDyn and Mooring Analysis Program (MAP). MoorDyn is an open-source lumped-mass mooring line model, developed at the University of Victoria, Canada, that supports arbitrary line interconnections, clump weights and floats, and different line properties. The model accounts for internal axial stiffness and damping forces, weight and buoyancy forces, hydrodynamic forces from Morison’s equation, and vertical spring-damper forces from contact with the seabed. MAP is developed by the National Renewable Energy Laboratory (NREL) in USA, and is designed to model the static and dynamic forces of a mooring system. The implementation of a multi-segmented, quasi-static mooring model in MAP is investigated by Masciola et al. [
95]. Seabed contact, seabed friction, and externally applied forces can be modeled with this tool, and it allows multi-element mooring systems, with arbitrary connection configurations, to be analysed. NREL, in partnership with Texas A&M University, have also developed an open-source dynamic FEM mooring model, FEAMooring [
199], for the analysis of floating offshore wind turbines, which may also prove useful for WEC applications.
The in-house software packages identified are: AQUA-FE, MooDy, MODEX and WHOI cable. AQUA-FE is a finite-element program developed at the University of New Hampshire to model partially and completely submerged structures and their moorings in the open ocean, and has mostly been applied to mooring systems for aquaculture applications. MooDy is a software application developed by the Department of Shipping and Marine Technology at Chalmers University of Technology, Sweden, for computing cable dynamics. MooDy was created as a tool for in-depth analysis of the tension force propagation in mooring cables and is designed to accurately capture the effect of snap loads. MooDy is based on a FEM with high-order polynomial basis functions and spatial discretisation using the local discontinuous Galerkin method. Chalmers University of Technology also developed MODEX, which is an older program that uses a FEM, with rod elements, to discretise a submerged cable. WHOI cable is proprietary software of Woods Hole Oceanographic Institution, built upon the groundwork from Howell [
106] and Tjavaras [
200], consisting of a collection of programs for nonlinear cable mechanics designed specifically to solve for typical oceanographic mooring systems.
Apart from the available software packages, it is, of course, possible to implement MMs by programming the governing equations in general purpose software environments. For example, Martinelli et al. [
201] present a model for the dynamic simulation of a catenary mooring line, by implementing the equations describing the dynamic movements of a chain in Comsol Multiphysics, and shows that the solution fully agrees with experimental measurements.
6. Applications in WEC Analysis
Mathematical models for WEC mooring systems are important for numerous applications in the design, analysis and optimisation of WEC systems, such as:
Mooring design
Analysis of mooring line fatigue
Analysis of extreme loads
WEC simulation
WEC control design and evaluation
WEC array layout
Each of these applications is presented in more detail within this section. A tabulated literature review is collated for each application, listing the type of WEC, mooring system and models used in each reference. (Note: For brevity, the following abreviations are used in
Table 5,
Table 6,
Table 7,
Table 8,
Table 9 and
Table 10: TD = Time Domain, FD = Frequency Domain, QS = Quasi-Static, PA = Point Absorber, M. = Mooring, NL = Nonlinear, L. = Linearised, RF = Restoring Force).
6.1. Mooring Design
The mooring design must be an integral part of the overall WEC design, because the mooring system not only represents a large fraction of the capital costs, but can also influence the WEC performance in both operational and survival modes. The objective is to find design solutions to mitigate peak mooring loads, while avoiding the costly overdesign or compromising the power production of the WEC. MMs play a vital role in the design process, where the design solution is usually obtained iteratively, since the loads on the lines depend on the characteristics of the lines themselves [
202].
The current state of mooring design for a number of Danish WECs is assessed by Thomsen et al. [
15], providing a good overview of the procedures and analyses involved in mooring design. Fitzgerald and Bergdahl [
22] and Martinelli et al. [
23] also detail many important aspects of mooring design for WECs. For example, Fitzgerald and Bergdahl [
22] uses MMs to compare different alternatives for the design of the mooring cables of a point absorber WEC in 50-m water depth. The MM analysis compares the different mooring designs in terms of the:
Overall weight and cost required to handle the baseline extreme loads
Mechanical impedance of the mooring cable at the attachment point and its likely influence on the motion of the floater
Seabed footprint and horizontal movement at the surface, due to steady or slowly varying loads, and its likely influence on array efficiency
Deployment considerations, required bollard pull for anchor embedding, etc.
From the results of the analysis by Fitzgerald and Bergdahl [
22], the authors conclude that the use of synthetic cables with higher elasticity has the potential to greatly reduce the overall weight of the system.
Discussion of the range of environmental inputs to be considered, when designing a WEC mooring system is given by Bergdahl [
32]. An illustrative example is presented to show the procedure for designing a mooring system, using static, quasi-static and dynamic models in different stages of the design process. A quasi-static design loop is used to initially design the mooring, and is then followed by sophisticated dynamic simulations. Dynamic calculations may first be run in the frequency domain, to allow many cases to be run, from which a few critical cases should be run in the time domain.
The procedure of initially using a static analysis to identify potential mooring designs, to be then further refined and optimised with high fidelity dynamic analysis, is also recommended by Fonesca et al. [
202]. The static formulation is extremely fast, providing estimates in a few seconds, while the dynamic formulation, depending on the complexity of the excitation and system, can take hours to provide stable statistical results [
202]. Likewise, Johanning et al. [
20] recommend using a static analysis, within a preliminary design approach, to estimate the maximum tension and stresses on the WEC mooring components due to extreme environmental conditions, from which the initial dimensions of the material components can be identified for further analysis.
Monarcha and Fonesca [
203] present a static analytical method for the preliminary design, of multiple line mooring systems, for floating offshore structures. The thesis by Wang [
182], presents an evolutionary optimisation study on offshore mooring system design, detailing the use of MMs with numerical optimisation techniques to provide tools to aid in the design and optimisation of offshore mooring systems.
6.2. Fatigue Analysis
Fatigue assessment of the mooring lines is crucial to ensure the integrity of the mooring system over its operational lifetime. A WEC experiences a multitude of different sea states (mild, moderate and severe) during its lifetime, and the fatigue assessment of the mooring system should be based on the aggregate of the load effects from the various sea states [
137]. Using a quasi-static MM significantly underpredicts the mooring loads, especially for fatigue [
89]; therefore, a dynamic model is required. For each sea state condition, a dynamic analysis is required to capture the stresses and loads on the system, encapsulating the HF system responses within a lengthy simulation duration, to ensure sufficient cycles for the LF responses are represented. Therefore, fatigue analysis for WEC mooring systems requires efficient modelling procedures.
The fatigue life of a mooring line is assessed by the number of stress cycles the mooring line sustains before failure. There are a number of techniques used to perform the fatigue assessment, ranging from the more accurate but computationally costly cycle counting methods in the time domain, to the less accurate but computationally faster frequency domain methods, such as the simple summation, the combined spectrum and the dual narrow-band approaches. A detailed description of the different methods can be found in the Position Mooring standard DNV-OS-E301 [
30] and in the thesis by Larsen [
207], which provides an assessment and comparison of the different methods for fatigue analysis of mooring systems.
Low [
137] and Low and Cheung [
208] outline techniques and model requirements for analysing the fatigue of moorings and risers in floating production systems for offshore oil and gas, presenting a hybrid time/frequency domain method to address the requirement of efficient modelling. Christiansen et al. [
153] use system identification techniques to train neural networks on a FEM MM output to allow efficient mooring line fatigue assessment. Thies et al. [
209] present a method to analyse annual fatigue conditions from a limited number of experimental tanks tests, and from field measurements [
210].
The thesis by Yang [
211] focuses on the fatigue characteristics in mooring lines for WECs. The study compares various simulation procedures for the analysis of fatigue of WEC moorings, with the objective of recommending the type of simulation procedure that can be used to make reliable fatigue assessment of WEC mooring systems, for a reasonable computational effort. In Yang et al. [
17], the authors compare using the coupled and uncoupled approaches and conclude that, although the uncoupled approach is believed to predict less accurately, the differences in terms of predicted fatigue damage is minor between the two approaches and, therefore, the uncoupled approach should be favoured due to the reduced amount of model preparation and computational effort required. However, in a later study by Yang et al. [
166], the same authors conclude that the coupling effects are crucial for the fatigue damage analysis of the mooring line.
Snap loads have been seen to cause significant damage in a number of experiments and field tests [
35,
209,
212]. Therefore, capturing the effect of snap loads is important for assessing the damage and fatigue to mooring lines. Palm et al. [
86] highlight that using numerical methods which are not locally conservative, such as standard continuous finite element or centred finite difference methods, will have problems in correctly capturing snap loads. Palm et al. [
86] present a finite element model, based on the discontinuous Galerkin method, that aims to capture the correct amplitude and duration of snap loads in mooring cables.
6.3. Extreme Loads
The occurrence of extreme, or peak, loads on WEC mooring systems must be carefully evaluated, to ensure a robust and efficient mooring design [
213]. One of the main challenges, for a commercially successful WEC, is to ensure survivability in extreme conditions, at a reasonable cost [
214]. Calculating the maximum tension in the mooring system requires a dynamic analysis; however, historically, quasi-static mooring analysis is used, and phenomena affecting the maximum line tension, that are neglected by this modelling procedure, were accounted for by an overall safety factor, with a typical value of three, for operational conditions, and two for survival conditions [
101].
The thesis by Liu [
134], focuses on the dynamic response analysis, and extreme value problems, of moored floating platforms. The usual process to estimate the extreme mooring line tensions is to perform many time domain simulations, which is computationally expensive. However, Liu [
134] shows that only a few simulations, plus the extreme value theory presented therein, are required to predict the extreme mooring line tension. Liu calculates the expected extreme mooring cable tensions, due to WF excitations, using statistical methods, on nonlinear time domain simulations [
215], and then extreme responses for moored floating platforms, under combined LF–WF excitation [
216].
To enable efficient analysis of peak mooring loads, Low [
136] presents a hybrid time domain—frequency domain approach for predicting extreme vessel motion, and line tension, for deep water floating systems. Muliawan et al. [
175] examine reducing the computational requirements for extreme loads on a WEC mooring system, through parsimonious selection of input sea states to be simulated. The environmental conditions to be considered, for extreme load analysis of WEC mooring systems, are investigated experimentally, from field test measurements, in Harnois et al. [
213].
6.4. WEC simulation
Numerical simulation is an essential tool in analysing and optimising the performance of a WEC. The mooring system will effect the WEC performance, adding nonlinear resistance and reactance to the WEC dynamics; therefore, a MM should be included into the WEC simulation to capture these effects. The effect of moorings on the motions of a point absorber WEC is studied in Fitzgerald and Bergdahl [
123], where it is demonstrated that moorings can influence the absorbed wave energy and that the influence depends not just on the type of cable, but also on where it is attached to the device, and in what degree of freedom of the floater that power is extracted.
To properly account for the way in which the mooring system alters the impedance of the WEC dynamics, some form of dynamic analysis should be used [
28]. The fidelity of the MM should match that of the rest of the WEC simulation; for example, if a frequency domain simulation of the WEC hydrodynamics and PTO system is being implemented, then the MM should be linear to allow a frequency domain implementation.
6.5. WEC Control
The effect of the mooring system on the WEC dynamics should be included when formulating the optimal energy maximisation control strategy. However, considering that mooring forces in the control formulation for WEC energy maximisation is relatively new, with none of the reviewed papers in
Table 9 existing pre-2013. Model-based control solutions, where real-time computation is an absolute requirement, is typically fulfilled using simplified models. Including mooring effects into the WEC model would therefore require a very computationally efficient MM to fulfil the real-time computation requirement.
Apart from energy maximisation, one important objective for MMs in WEC control, not currently reported anywhere in the literature, is the minimisation of mooring fatigue or extreme loads, using control strategies. The objective of minimising mooring fatigue using control strategies is similar to the work of Ferri et al. [
228] and Tom et al. [
229], who investigate minimising structural fatigue of a WEC via control. Parmeggiani et al. [
214] investigate using control strategies to reduce the peak mooring loads in extreme wave conditions, and presents a method which is able to reduce the extreme loads in the main mooring line by 20–30% in wave conditions over a 100 year return period; however, this is an exclusively experimental investigation and no MMs are used.
6.6. Arrays
The complexity of WEC mooring systems increases for the case of WEC arrays, with the number of degrees of freedom of the system increasing with the number of bodies and mooring lines within the array which, in turn, increases the complexity of the MM. Additionally, factors such as stricter restrictions on device excursions to maintain optimal spacing between devices, anchor sharing, and number of inter-body lines, add to the complexity of the analysis of mooring systems for WEC arrays.
Integrated mooring systems may offer cost reductions, compared to individual moorings for each WEC in the array. This was investigated by Gao and Moan [
24], who analysed the WEC motions and mooring tensions for an array compared with the motions and tensions for a single WEC in operational, and survival, conditions. Large mooring line tension in horizontal inter-body mooring lines were observed when the dynamic motions of two adjacent WECs were driven by waves with a phase difference of 180 degrees, resulting in a very high demand on mooring line capacity. In this case, the individual mooring system would prove more feasible than the integrated mooring system, however, by increasing the line lengths and by using clump weights, the tensions of the inter-body lines in the integrated mooring system could have been reduced.
8. Conclusions
Mathematical modelling for WEC mooring systems borrows many established techniques from other offshore engineering fields, due to the common need to design, analyse and optimise moorings for floating structures. Various modelling techniques of increasing complexity and fidelity are available, with the most accurate model types being very computationally costly compared to their simpler counterparts. The type of model will depend on the operating principles of the WEC, the type of mooring system used, and the analysis application.
The full range of effects to be captured by MMs spans many time scales, covering very HF effects such as snap loads and dynamic tensions, WF oscillations, LF drift motions, and very LF tidal variations. This requires the simulations to have very small time steps, to capture the HF effects, but long simulation durations to capture the LF effects. Additionally, the WEC and mooring system models should be subjected to a range of sea states and environmental conditions, which need to be analysed by separate simulations. Therefore, pragmatic choices need to be made regarding the type of model to use and which inputs to consider, to achieve the best balance between computational requirements and accuracy of results.
In this review, the literature reporting the usage of MMs is collated and separated into the application areas of design, fatigue analysis, extreme loads, WEC simulation, WEC control and array analysis. Fatigue and extreme load analysis are seen to require the highest fidelity MMs, able to capture the dynamic tensions and effects throughout the mooring line. WEC simulation and control applications mainly require the effect of the mooring system on the WEC dynamics, i.e., the mechanical impedance of the mooring system, which generally allows more simplified MMs to be utilised. Simpler models are useful for the real time computational requirements of WEC control, and for the long simulation durations and numerous sea state simulations required for power production assessments from WEC simulations. Mooring design and array layout optimisation applications are seen to utilise simple models, initially, and then increase the model fidelity, as the design and analysis becomes more refined.