AI-Driven Model Prediction of Motions and Mooring Loads of a Spar Floating Wind Turbine in Waves and Wind

1. Introduction

The search for sustainable and reliable energy solutions has increasingly turned toward the oceans, with Floating Offshore Wind Turbines (FOWTs) emerging as pivotal technologies in harnessing marine wind resources [1]. The complexity and stochastic nature of marine environments, characterized by the dynamic interactions of waves, wind and currents, presents significant challenges in the design and operation of FOWTs. Understanding the hydrodynamic responses of these structures, especially under extreme conditions, is crucial for their safety, efficiency, and longevity. Engineers and designers frequently conduct evaluations and forecasting using various [2] numerical time-domain coupled analyses. In this context, OpenFAST has emerged as an integrated tool, capturing the complete behavior of wind turbine systems –including the turbine, support structure, mooring lines, and the marine environment– to facilitate a comprehensive understanding of system dynamics [3]. While foundational, these studies require large computational resources, and face challenges in accurately predicting FOWTs’ non-linear responses to varied sea states, requiring extensive simulations and tunning for validate them.

To obtain the dynamic response of FOWT at a reasonable time, some researchers have improved the time-domain coupling analysis methods. Many of these approaches employ a decoupled analysis in the time domain, treating the platform’s mooring lines as quasi-static springs while enhancing the model by incorporating the dynamic effects of the lines through equivalent linear damping and inertia coefficients [4,5]. However, these approaches have some inherent limitations, as drag-induced damping is non-linear and thus dependent on the amplitude of the oscillation. More recently, Low and Langley [6] introduced a hybrid method that simulates low-frequency motion in the time domain and wave frequency motion in the frequency domain, thereby reducing the computational load and enhancing the efficiency of time-domain evaluation methods. Nevertheless, the methods’ computational demand, though lesser, still poses challenges, highlighting an ongoing quest for more efficient predictive models.

Parallel to these developments, the field of artificial intelligence (AI), particularly data-driven neural network models (ANN), has shown promise in capturing temporal dependencies and complex patterns in sequential data. Many scholars and researchers have considered selecting the appropriate model to resolve issues related to time series prediction problems [7]. Feedforward Neural Networks (FNNs) are broadly used in many time series forecasting applications [8,9,10]. Nevertheless, while effective for certain tasks, these simple models often fall short when handling the non-linearity and long-term dependencies inherent in complex time series data. This limitation has led to the exploration and adoption of more advanced neural network architectures, such as Recurrent Neural Networks (RNN), proposed in the 80’s, which are capable of storing information for a long time in hidden states [11]. RNNs achieve this by using loops within the network structure, allowing information to be passed from one step of the network to the next. This looping mechanism enables the network to maintain a form of memory, capturing dependencies across different time steps. In particular, one of the self-looping-based RNNs, the Long Short-Term Memory (LSTM) [12,13], has gained significant attention for its ability to understand long-term sequential dependencies, making them exceptionally suitable for forecasting tasks. Consequently, LSTMs present a groundbreaking methodology for predicting seakeeping dynamics, leveraging large response datasets.

The ability of ANN models to learn complex long-term dependencies at low computational cost suggests a potential paradigm shift in how to obtain the dynamic response of FOWTs. It aligns with recent researches Bjørni et al. [14] where explored machine learning’s impact on predicting spar floating wind turbine mooring tensions using motion time series as inputs, finding generally strong correlations except in extreme sea states. Silva and Maki [15] crafted an LSTM neural network for simulating a ship’s 6 degrees of freedom (DoF) motions at 20 knots in irregular seas, evaluating model accuracy against wave inputs and the training data quantity. Most recently, Lee et al. [16] presented novel integrated architectures combining an LSTM encoder and decoder with a convolutional neural network have been developed. These have been utilized to forecast ship motion, taking the encoded initial state vector and anticipated ocean wave field information around a vessel.

In response to these sequence learning challenges in modelling FOWT platforms, this study introduces a novel LSTM-based neural network model, implemented using PyTorch, that transcends traditional methodologies in simulating the dynamic responses of FOWTs. It bridges the existing gap in the literature regarding the application of AI techniques to infer time-domain responses for FOWTs, using as input operational and severe sea states and wind conditions. Leveraging a comprehensive dataset derived from over 300 fully coupled simulations using OpenFAST for Time-Domain (TD) and ANSYS AQWA for Frequency-Domain (FD) of the OC3 reference design [17]. The model developed not only forecasts the motion time series but also extends its predictive power to fairlead tensions under operational and extreme wind and wave conditions. This approach promises to accelerate the prediction process of FOWT dynamics while revolutionising real-time monitoring and control strategies, thereby improving offshore wind energy systems’ design and operational efficiency. Through rigorous validation and performance evaluation, our research sets a new benchmark in seakeeping analyses in marine sciences, paving the way for more resilient and efficient renewable energy solutions.

This paper is organized as follows: first, the definition of the case study; second, the description of the numerical simulation model and the neural network setup applied in the present research; third, the results obtained and their evaluation; and finally, some conclusions and future lines of works are discussed.

2. Case study

2.1. Floating Wind Turbine Model

The OC3-Hywind model [17] is a reference wind turbine that comprises a spar-type floating wind turbine (FWT) of 5 megawatts, anchored by three catenary mooring lines. This concept has been chosen in the current work as proof of concept due to its extensive validation of both experimental [18,19] and numerical data [17,20], providing a robust foundation for comparative analysis of the performance and accuracy of the numerical simulations carried out in the present study.

The draft of the platform is 120 m and, between 4 m and 12 m below the sea water level (SWL), extends a linearly tapering conical transition region that connects two cylinder diameters of 9.4 m below the taper with a slender one of 6.5 m. Figure 1 illustrates the concept. More details on the platform’s inertia, mass and mooring properties can be found in NREL Report [17].

The 5-MW baseline wind turbine conventional controller relies in a generator-torque controller and a full-span rotor-collective blade-pitch controller [21]. Each one can work independently; the generator-torque controller maximizes the power capture below the rated operation point while the blade-pitch controller regulates the generator speed above the rated operation point.

2.2. Database generation for training and validation of the neural network model

This FOWT was simulated under different sea states corresponding to Aberdeenshire’s coast (see Figure 2). This location is significant as it is where the 30 MW Hywind Scotland pilot wind farm is situated [22]. Utilizing data from the wave scatter diagram, a selection of 300 irregular waves paired with turbulent wind realizations formed the Latin Hypercube sampling employed in the neural network analysis. The mean wind speed at the hub height for each wave scenario was determined using the IEC 61400 [23] power law (see Eq. (1)).

$$U\left(z\right)=U_{10}{\left({\displaystyle \frac{z}{10}}\right)}^{\alpha },$$

(1)

with the recommended exponent $\alpha$ of $0.14$[23], and z, representing the vertical distance in meters from SWL ($90\mathrm{m}$ in this study for the hub). A range of wind speeds from $6\mathrm{m}/\mathrm{s}$ and $33\mathrm{m}/\mathrm{s}$ (above the cut-out limit for the present controller) were considered. This range includes both operational and extreme cases. One-hour simulations using OpenFAST v3.5 were conducted with different seeds (randomly uniform distributed) of the selected wind and wave spectra combination, ensuring comprehensive coverage of the FOWT’s dynamic responses under varying environmental conditions.

In the current simulations, the water depth is set at 320 meters, following established literature, with both waves and wind propagating along the x-axis (180-degree heading).

3. Numerical Simulations

OpenFAST is a time-domain open-source simulation tool to analyse wind turbines. In the numerical model, the spar FOWT system has been modelled using the platform properties described in Section 2.2. The platform’s hydrodynamic calculations were obtained by merging potential flow theory with Morison’s equation within OpenFAST’s HydroDyn module. The hydrodynamic coefficients obtained from frequency domain analysis are utilized in the time domain by replacing the integral over all frequency-dependent incident-wave-excitation forces [25] with a convolution integral, which in HydroDyn is implemented using Fourier transforms as outlined by Oglive [26] and Duarte et al. [27]. In the present research, although HydroDyn (OpenFAST) primarily processes data files generated by WAMIT, ANSYS Aqwa, a commercial software for undertaking hydrodynamic and mooring analyses, was employed for the frequency domain analysis.

Then, the model has been calibrated and fine-tuned using data from the available literature [28]. Similar methodology was adopted in [20]. For all DoFs, numerical decay tests were performed to calibrate the numerical model. This fine-tuning consisted adding a quadratic viscous-drag term ($C_D$) equal to 0.6 in horizontal direction, and an additional linear damping matrix in surge, sway, heave and yaw, which aligns with the NREL Report [28]. Figure 3 presents an example of a time series of platform motions for both the experimental decay test and the numerical decay test of the tuned numerical model.

Furthermore, motion RAOs were obtained by running a white noise spectrum to validate the numerical model with the available numerical and experimental literature ([17,29,30]). In Figure 4, a comparison of the motion RAOs obtained with the developed numerical model and the ones obtained in [17] is displayed.

In addition, the mooring system presented in [31] is modelled with MoorDyn (OpenFAST) [32], which takes into account the effect of inertia, hydrostatics, axial elasticity, structural damping, hydrodynamic added mass, and hydrodynamic drag. Bending stiffness is neglected since a chain is considered for the mooring system. The line’s contact with the seabed is modelled with a series of spring dampers.

A three-dimensional turbulent wind field was generated according to the mean wind speed at the hub height, using the power law explained previously. TurbSim tool [33] was employed to simulate a one-hour wind time series for each case with a random realization of the Kaimal spectrum. InflowWind module processes these wind time series over the turbine and loads them in OpenFAST. Each irregular wave spectrum was generated using a random phase by HydroDyn with seeds chosen uniformly.

Also, the aerodynamic loads on both the blades and tower were calculated within OpenFAST’s AeroDyn module. To model the operational rotor, we used dynamic Blade-Element/Momentum (DBEM) with an unsteady airfoil aerodynamic (UA) model; when the mean wind speed overpasses the cut-out limit, the rotor is set as idling allowing it to spin freely at low speed so the aerodynamic loads are neglected [34].

Finally, the structural dynamics of the wind turbine is computed with the ElastoDyn module, which receives the forces and the reactions and obtains the position, velocity and accelerations of key points of the structure. In scenarios where the turbine is not operational, the initial rotor speed is set to zero, and the blade pitch is in feathered condition. Furthermore, for these cases, the startup of the generator was disabled and some speed-independent drag was added to the drivetrain by applying a light brake load in the ServoDyn modules [35].

4. Neural Network Model

Although several data-driven and Machine Learning (ML) models have been designed and developed for floating bodies motion prediction in the past, considerable challenges still remain. Regarding the Artificial Neural Network (ANN) techniques, feedforward Neural Networks models (FNN) have been used in many time series forecasting [9], including FOWTs mooring lines’ tension logs prediction [14]. A FNN model is a network without recurrent links, which means the signals only pass in one direction. Data flows through the input layer, where each element calculates a weighted sum of its inputs. These computed values are then fed into the subsequent layer, and this sequence continues through each layer until the final output is produced. The output is finally optimized with the back propagation method. Nevertheless, while these basic models are effective for some tasks, they frequently struggle with high non-linearity and long-term dependencies, characteristic of complex time series data.

To address these limitations, Recurrent Neural Networks (RNNs) were developed in the 1980s [11]. Their architecture with recurrent connections enables the network to maintain a memory of previous inputs in its internal state, which is essential for capturing temporal dependencies in sequential data. However, RNNs face challenges with vanishing or exploding gradients during training. As a result, when dealing with longer sequences, RNNs may struggle to retain earlier information, complicating the learning of dependencies between distant elements within the sequence.

For this reason, in the present study, the authors employed an LSTM-based Neural Network, which mitigates the vanishing gradient problem, implemented using PyTorch. For the FOWTs dynamics, excited both by waves and wind, the response depends not only on the current state but also on the previous states. This approach, distinct from traditional FNN, ensures more accurate and reliable predictions of the complex, time-dependent behaviour of FOWT platforms.

4.1. Model Architecture and Set-up

The developed neural network includes an input linear layer that works as an encoder and a single linear output layer that acts as a decoder, outputting the time series for each DoF motion. The hidden layers contain five LSTM layers designed to find the complex dependence between the inputs and outputs and the sequential dependence with the previous time step ($h_{t-1}$). These layers, each comprising 128 cells, employ the tanh activation function ($\sigma$). This facilitates the model’s ability to learn complex mappings from inputs to outputs and is particularly effective with the normalization used herein (between $-1$ and $+1$).

The LSTM receives the wave elevation and the wind speed time series at each local axis of the nacelle as input vectors (${\eta }_{t_i}$) for a time vector t. The output (${\xi }_{t_i}$) from the developed model is each DoF motion or fairlead tension for the combination of excitations. During training, the network performs a forward pass (Eq. (2)) using updated weights ($W_{ih}$ and $W_{hh}$) and biases ($b_{ih}$ and $b_{hh}$) and computes gradients for backpropagation.

$${\xi }_{t_i}=\sigma ({\eta }_{t_i}{W}_{ih}^T+b_{ih}+h_{t-1}{W}_{hh}^T+b_{hh}).$$

(2)

The losses measure the discrepancy between the predicted outputs and the reference values, quantifying the performance of the model by indicating how well it is learning and generalizing from the training data. The loss function used herein was Mean Squared Error (MSE) averaged over all datasets (see Eq. (3)).

$$\mathit{L}(\xi ,\widehat{\xi })={\displaystyle \frac{1}{N}}\sum _{j=1}^N\left({\displaystyle \frac{1}{T}}\sum _{t_i=1}^T{({\xi }_{t_i}^j-{\widehat{\xi }}_{t_i}^j)}^2\right),$$

(3)

where N is the total number of cases used in each dataset, T is the total length of each time series, ${\xi }_{t_i}^j$ is reference output, and ${\widehat{\xi }}_{t_i}^j$ is the predicted one.

Optimization of the hyperparameters was carried out using the Ray Tune tool. Sampling rate, sequence length, number of layers, activation function and optimizer were some of the hyperparameters optimized to obtain the best results in the minimum training time. Figure 5 displays the validation loss over the different trials. The test coloured in blue remains for the final hyperparameters employed in the NN. Despite not being the lowest loss, the improvement in loss over the computational training time of applying another hidden size or sampling rate combination is minimal.

The final model configuration employed the Adam optimizer [36] to minimize the loss during backpropagation. The chosen sequence length was 80 time steps. An initial learning rate of $0.005$ was implemented, with a decay schedule that reduced by 0.9 each epoch. The detailed hyperparameters used in the presented model are listed in Table 1.

The time series output from numerical simulations of each corresponding movement and mooring fairlead tensions are used as the training output of the network. The model’s inputs include excitation logs, i.e. sea surface elevation and turbulent wind speeds vectors. Both inputs and outputs are normalized to a range between $-1$ and $+1$ based on their respective minimum and maximum values. Scaling of the features prevents vanishing or exploding gradients issues and enhances numeric stability during training.

The input data, originating from numerical simulations with an initial sampling frequency of $100\mathrm{Hz}$, was down-sampled to $2\mathrm{Hz}$ to prevent redundant information and forcing the model to focus on long-term dependencies. High sampling frequencies can lead to similar consecutive data samples, causing redundancy. The Nyquist theorem [37] was checked when selecting the minimum re-sampling frequency.

Finally, the dataset was randomly divided into training ($65\%$) used by the model to learn patterns, validation ($15\%$) to tune hyperparameters and prevent overfitting, but not used to update the model weights and bias, and testing ($20\%$) to assess the model’s generalization performance on unseen data and final check. This split into subsets ensures robust model evaluation. All the datasets have been preprocessed to check if the statistical data from this split are similar. Figure 6 presents the rotated kernel densities from all the subsets from input features, i.e. wave elevation and wind speed in X, Y and Z directions respectively from left to right. It shows that similar distributions and means values are maintained across the training, validation, and testing datasets, indicating consistent statistical properties and ensuring the reliability of the model’s performance evaluations.

5. Results and discussion

The presented ANN model was trained to predict a 30-minute time series for all DoF motions and fairlead tensions, requiring 100 epochs and over 60 hours of computational time in total, with the model running on a GPU with 8 GB memory. Based on the sea state descriptors, 30 minutes of wind and wave excitation represent the forcing series [38,39,40], ensuring that the system’s dynamics and any longer-term patterns and oscillatory events in the tail are adequately captured and modelled.

This duration represents a significant improvement over similar previous studies, which only achieved accurate predictions for durations of up to 150 seconds [14,15,16]. Moreover, the model takes less than 5 seconds to infer 30 minutes of six DoF motions and the three fairlead tensions for each sea state, while the same case, numerically simulated within equal computation conditions, takes more than 10 minutes (a two orders of magnitude gain).

One important aspect to assess the robustness of this type of the models is to show that they do not present neither overfitting nor underfitting issues[41,42,43]. In the case of overfitting, the learning model can only memorize the training data set, displaying poor results for inputs other than such training data set. In the case of underfitting, the learning model is too simple and cannot capture the underlying data patterns .

To verify that the model had no overfitting issues throughout the training, the losses for both training and validation datasets were monitored until they converged to a low MSE averaged value for all the predicted features. The loss curves for heave presented in Figure 7 motion exhibit similar trends for both training and validation losses, further indicating the absence of overfitting.

In regards to underfitting, from the analysis next it will be made clear that the model is able to accurately predict the responses for the test data set.

5.1. Discussion of Selected Test Cases Through Analysis of Time Histories

As reported in Section 2.2, the simulated cases only include 180-degree wave directionality. For this reason, and considering the symmetry of the configuration with respect the x-axis, the results will focus on the excited DoFs, i.e. surge, heave and pitch motions, and the three fairlead tensions. Four representative cases from the test dataset have been selected for the present discussion:

a)

Irregular wave case $H_S$ 2m and $T_P$ 8.5s and turbulent wind with $U_{90}$ 12m/s represents the operational wave with higher occurrence probability.

b)

On the other hand, case corresponding to irregular wave with $H_S$ 3m and $T_P$ 7s and wind speed of 6m/s for $U_{90}$, corresponds to the lowest peak period from the test dataset.

c)

In addition, case with $H_S$ 2m and $T_P$ 21.5s and $U_{90}$ 12m/s is in the highest peak period limits, with low observed events.

d)

Finally, the irregular wave case with $H_S$ 9.5m and $T_P$ 17s and wind speed for $U_{90}$ of 26 m/s exemplifies the extreme sea state for the utilized scatter. In such condition the turbine is parked/idling.

A selection of thirty minutes of predicted outputs from these cases is presented in Figure 8 for motions, and in Figure 9 for mooring loads.

The results for motions fairly match with the reference values, presented in the grey dashed line in Figure 8. Comparing the results between different DoF motions, the heave response exhibits a less accurate prediction compared to surge and pitch. This discrepancy can be attributed to the relatively smaller amplitude of the heave motion compared to the others of the DoFs, despite normalization, making it more prone to larger relative errors. Additionally, the heave motion’s principal frequency response is significantly larger than the other predicted responses, i.e. surge DoF, suggesting the need for adjustments in the sampling frequency and sequence length to improve prediction accuracy. The comparison across different sea states and wind conditions (that will be latter discussed in more detail) reveals that, in general, the LSTM model demonstrates strong predictive accuracy capturing the overall trends. In low wave and wind excitation, such as lower peak period (Figure 8(a), (b) and (c)), all motion predictions closely follow the actual signals, showing the model’s robustness to varying periods. In extreme conditions (Figure 8(d)), the surge and pitch motion predictions remain reliable, albeit with major discrepancies struggling to capture the higher rapid oscillations accurately. Heave predictions continue to show larger deviations.

Furthermore, Figure 9 displays the predicted tensions for each mooring line, with a good correlation with reference values. Fairleads 2 and 3 correspond to the up wind ones facing the incoming waves and wind and Fairlead 1 the down wind one. Fairleads responses are largely dominated by the second-order surge mean drift response. Nevertheless, the model also captures the wave frequency, as can be seen in the operational case (Figure 9(a)). For the extreme case (Figure 9 (c)) there are no noticeable differences between fairleads. In the next section, a detailed evaluation of the model’s accuracy is carried out.

5.2. Evaluation of Selected Test Cases Using a Probabilistic Approach

Although the overall time series comparison in Figure 8 and Figure 9 suggests a strong correlation between the estimated and the reference data, further analysis has been conducted by obtaining and evaluating the heteroskedasticity ($R^2$ score) and comparing the statistics of the predicted values to those of the reference data. These metrics provide deeper insights into the model’s performance, assessing its precision in a probabilistic framework, which is useful from the industrial point of view.

Figure 10 shows the predicted and reference values for the wave cases displayed previously. Each point in the scatter represents a single time step response amplitude –motions and mooring tensions– prediction where the input wave elevation and wind speed are identical between numerical simulation and LSTM inference. The exact 1:1 correlation would be in the straight black line. $R^2$ values are also presented for each case. In the same figure, the residual values, defined as the difference between the obtained and reference values, are presented with a grey scale legend.

It can be appreciated that the residual distribution for surge response is similar for all values and the different selected cases, with $R^2$ coefficients of around $0.97$. Pitch motion prediction displays a high degree of accuracy, with $R^2$ values above 0.9 for all the samples but lower than for surge predictions.

Regarding heave, the prediction presents a lower accuracy except for the sea state with smallest peak period. For this case (Figure 8 (b)) the response is rather regular which might help the LSTM model to predict more accurately, leading to the largest $R^2$ in Figure 10 (b). It is remarkable that for the largest peak periods (cases c and d) the accuracy in the prediction on heave motions falls, in particular, for the extreme sea state, with $H_S=9.5\mathrm{m}$.

The inference of the tensions at the fairleads outputs $R^2$ scores above $0.95$ in all the cases shown indicating that the predictions are rather accurate.

In addition to comparing values in each time step, the Cumulative Distribution Functions (CDF) and the Probability Density Functions (PDF) of the reference and obtained values have also been compared for a comprehensive understanding of the model’s performance. In the case of the CDF, the Probability of Exceedance ($1-CDF$), PoE, is actually displayed instead, as it is of more practical value. PoEs and PDFs are presented in Figure 11 for all the cases selected for discussion. The PDFs were estimated from the CDFs using the Kernel Density Estimation (KDE) method [44].

Upon the statistical study, the Two-Sample Kolmogorov-Smirnov Test (KS-Test) [45] was performed over all motion predictions to compare whether the two independent samples (reference and inference) were drawn from the same probabilistic distribution. In Figure 11 the K-S Test p-Value results are also provided. In addition, the figure displays the 95% confidence intervals of the CDFs and PDFs, in shaded blue, computed using the Moving Block Bootstrap (MBB) method technique [46].

In the operational case, Figure 11(a), the PDF of the predicted signal (blue dashed line) closely follows the reference signal’s PDF (grey solid line) for surge motion. However, for pitch and heave motions, slight deviations are visible around mean values. Focusing on the PoEs, they fairly match the reference ones, with larger deviations for heave at the upper tail. For pitch, it shows smaller predicted values in the central part of the distribution.

For cases with both lower and upper wave periods, Figure 11 (b) and (c), present, respectively, the predicted PDFs. It can be seen that they align well with the reference PDFs for surge and heave motions, with minor deviations. Some deviations are visible, especially in the central region of the distribution for pitch DoF.

The comparison of PoEs and PDFs for the extreme sea state, Figure 11 (d), both for surge and pitch motions, display a good correlation except with minimum over predictions in the peak part of the surge distribution. Nevertheless, in heave motions, the PoE shows some differences for both tails.

Focusing on the KS-Tests, the resulting p-Values are above $0.05$ for most of the cases, indicating no significant difference between the predicted and reference distributions, thus confirming the reliability and accuracy of the model’s predictions even under extreme conditions.

5.3. Global Evaluation of All Test Cases Using a Probabilistic Approach

Each test case can be analyzed through its $R^2$ score and the p-Value from K-S Test. To represent the accuracy of the model for all the cases in the test dataset split (see Section 4.1), heat maps along the whole Convex Hull of such dataset are displayed in Figure 12 and Figure 13. These figures show the general behaviour of the model depending on the inputs, i.e. significant wave height –or wind speed– and the peak period.

In Figure 12 the $R^2$ score for the three DoF motions and the three mooring tensions is shown. The three motions present similar distributions, with lower correlation with the reference data in the cases with large $H_S$. This general trend was also observed for the selected cases in Section 5.1. The dependence with $T_P$ is instead negligible.

In the same Figure 12, tension predictions display high correlations for all the fairleads, with similar results for both Fairlead 2 and Fairlead 3, as expected.

The average across the whole test dataset of the $R^2$ correlation indicator has been computed for all the DoFs and fairlead tensions. Except for heave, their values are all above $0.9$ ($0.94$, $0.91$ for surge and pitch, and $0.98$, $0.94$ and $0.94$ for the three tensions). However, for heave response, this value is $0.78$. This lower correlations is expected considering the analysis from the selected cases in Section 5.1.

Figure 13 focuses on the p-Values from the K-S Tests presented as pass (1), no-pass (0) binary heatmaps. It can be noticed that for surge and pitch motions, most of the cases pass the K-S Test of equal distribution between the reference and the predicted data. It is remarkable that for low $H_S$ there are some cases that fail the K-S Test. This could be related to the reduced response under these conditions, naturally leading to larger relative deviations.

Regarding heave DoF, the heatmap reveals a less favorable outcome, with several cases not passing the K-S Test. Contrary to surch and pitch, the majority of cases that pass the test correspond to low $H_S$ values. Indeed, it is observed that higher $H_S$ conditions induce greater variability and non-linearity in the heave motion, making it more challenging for the model to capture the full statistical distribution of the data with accuracy. All the motions display no dependence on $T_P$, as for the $R^2$ analysis.

As for the fairlead tensions, one can observe that most of the cases pass the K-S Test –mainly in Fairlead 1 with $90\%$ of the cases –. The majority of the failed cases are located in the range of high significant wave height consistently with their larger dependence on the surge motion.

In brief, the presented model demonstrates high accuracy and reliability for most dynamic responses, particularly in surge and pitch motions and mooring tensions, with strong statistical consistency across various conditions.

6. Conclusions

In this paper, a comprehensive study was conducted on the prediction of the dynamics of the OC3-Hywind Floating Offshore Wind Turbine using an LSTM-based neural network model. The model was trained and validated with extensive datasets from fully coupled simulations using OpenFAST under various sea states and wind conditions, including operational and severe sea states. The following conclusions have been derived:

• The LSTM model effectively predicted time series for surge, heave, and pitch motions, as well as fairlead tensions under both operational and extreme wind and wave conditions. The model demonstrated high predictive accuracy, particularly for surge and pitch motions, with $R^2$ values generally above $0.9$.
• Statistical evaluations, including Probability Density Functions (PDFs), Cumulative Distribution Functions (CDFs), and Kolmogorov-Smirnov (K-S) tests, confirmed the reliability of the model’s predictions. The majority of the cases passed the K-S test, indicating that the predicted and actual distributions are very similar.
• Although the model performed well overall, the heave motion predictions were less accurate compared to surge and pitch. This discrepancy is attributed to the smaller amplitude of heave motions and their higher frequency response, indicating the need for further adjustments in sampling frequency and sequence length.
• The model significantly reduced the computational time required for predicting FOWT dynamics. While traditional numerical simulations could take more than 10 minutes to compute, the LSTM model inferred 30 minutes of time series data in less than 5 seconds. This reduction makes the present proposed approach relevant mainly for fatigue analysis, with aiming to discard preliminary designs.

In summary, the reliability of the developed LSTM model in predicting FOWT dynamics enable this type of models as promising tool for its application in real-time monitoring and control strategies. This approach promises to enhance the design and operational efficiency of offshore wind energy systems.

Future research efforts will address extending this study to different types of FOWTs and incorporating wave and wind directionality and misalignment. Further optimization of the model’s hyperparameters and architecture could also improve prediction accuracy, especially for heave motions.