Vortex-induced vibration effect of extreme sea states over the structural dynamics of a scaled monopile offshore wind turbine

In order to demonstrate the relevance of considering Vortex-Induced Vibrations (VIV) in the structural design of marine structures, this study proposes an alternative experimental and analytical approach in wet conditions to measure the fluid–structure interaction in the near field and quantify the viscous damping with measured structural and 3D hydrodynamic accelerations. It was demonstrated that VIV caused and incremented 5.00% of the structural damping coefficient, and the extreme wind loading increased 74% of the offshore monopile’s structural damping, demonstrating the relevance of the high non-linear hydrodynamics effects during selecting parameters into the structural design in offshore applications.


Introduction
The development of offshore engineering projects for nonconventional and marine energies shows worldwide growth, where offshore wind energy technology receives considerable financial resources for research and technological development (Rueda-Bayona et al. 2019). The offshore energy sector's opportunities bring challenges, such as the mitigation of environmental and economic impacts because of structural failures caused by extreme sea states (ESS) (Xie and Aly 2020). These hazards demand rigorous structural B Juan Gabriel Rueda-Bayona juan.gabriel.rueda@correounivalle.edu.co Andrés Guzmán faguzman@uninorte.edu.co design for developing effective damping systems and a better understanding of non-linear fluid-structure interactions. In this sense, the fluid-structure interaction of the hydrodynamic field of offshore monopile wind turbines has been studied, resulting in viscous damping that affects the structural dynamics, thus, in the structural design, vortex-induced vibrations (VIV) generated during extreme sea states are not commonly measured; its estimation for the structural design is relegated to hydrodynamic coefficients with a secondary contribution to the structural excitation frequencies. Also, offshore structural parameters have been identified in dry conditions through shakers and wind tunnels. The aforementioned common practice could be dangerous because structural amplifications may be ignored. After all, the VIV effects are not considered.
Several studies have analysed the effects of linear and nonlinear wind, waves, and currents loads over fixed structures, estimating the natural and damped period of the response (Aggarwal et al. 2017;Marino et al. 2017; Bajrić et al. 2018). Other studies highlighted that wind and wave loading mitigated vibrations of their offshore Monopile Wind Turbine (MWT) during earthquake excitations (Yang et al. 2020). As evidenced further, there is not a significant quantity of studies devoted to coupling physics considering loads acting simultaneously in the specialised literature. The environmental loads (wind, waves, and currents) have been one of the main causes of accidents on-board and structural failures during the last decades provoking operational impairment of the offshore structures (Kandasamy et al. 2016). This explains why numerical and physical modelling to calculate the structural design parameters under complex non-linear environmental loading conditions have been widely studied in the design of the offshore structures (Jafarabad et al. 2014; Rueda-Bayona et al. 2018b). Bachynski et al. (2019) assessed the structural response of a 5 MW offshore MWT under irregular wave loading. They evidenced that the numerical results overestimated the second-order wave excitation, miscalculated the third-order wave excitation, and incorrectly calculated the low structural damping. Other studies applied wave loads calculated from a synthetic-free surface generated by the JONSWAP spectra. These studies did not show details of the fluid-structure interaction because of the hydrodynamic vortex generated by the waves approaching the structure ( Arany et al. (2017) used a method for designing monopiles of an offshore wind turbine for defining types and sizes for financial viability, who pointed the application of linear wave theory for estimating mild and extreme wave loads. The method did not consider the critical structural responses caused by the fluid-structure non-linear interactions, because the calculation of hydrodynamic forces with the linear theory cannot determinate the complex structural responses which must be controlled with specialised viscous dampers. Suja-Thauvin et al. (et al. 2018) carried out numerical and physical modelling of a 4 MW wind turbine at model scale concluding that model assumptions and simplifications cause inaccuracies because of the overestimation of total structural responses, resulting in overinvestments in the offshore engineering project and cost overruns because of the increased repair and maintenance activities. These studies could be limited when simulating fluid-structure interactions, because CFD solvers might omit frequencies for exciting the first mode of structure oscillations (Paulsen et al. 2013). Although numerical modelling is fast and inexpensive for determining structural dynamic parameters in offshore structures, physical modelling with instrumented scaled models in tanks or wave channels (flumes) is strongly recommended.
Computational Fluid Dynamics (CFD) and numerical approaches are the most used methods to evaluate the nonlinear responses of offshore structures because of wind, waves, and current loads (Lee and Mizutani 2009;Li et al. 2011;Rueda-Bayona et al. 2021). Colwell and Basu (Colwell and Basu 2009) applied the Kaimal and JONSWAP spectra to generate loads of wind and waves to assess the structural response of an offshore wind turbine. Wei et al. (2013) utilised the JONSWAP spectra to generate wave loads to develop a modal analysis of an offshore wind turbine considering the sea states in the Korean sea. Mockute et al. (2019) analysed numerically the exerted forces of combined and misaligned wind-wave loading over a MWT during rouge seas and reported that wave kinematics theoretical models did not capture properly the time evolution of loading force due to the local non-linear hydrodynamics. Although these approaches are relevant to consider non-linear responses, at that moment, these authors did not validate their models through experimental procedures.
Other studies have conducted physical modelling of offshore wind turbines in dry conditions using shakers or vibrating mechanisms in the absence of water (no wave flumes) (O'Kelly-Lynch et al. 2020). Ou et al. (2007) developed a damping isolation system to control vibrations of a jacket platform for determining structural parameters. Also, Mojtahedi et al. (2013) applied forced vibration tests with a shaker in dry conditions to analyse the structural health during several damage conditions. Hosseinlou and Mojtahedi (2016) assessed a physical scaled model of a jacket offshore structure through vibration tests in dry conditions to identify initial numerical parameters for structural health purposes.  2017) analysed experimental results of two offshore MWT scaled models (flexible and stiff), measuring the wave elevations and structural accelerations during several sea states generated by the JONSWAP spectra. Also, the free surface (wave loading) aside from the structure was measured and did not consider wind loading and not reported the damping ratios of the rigid monopile. Maes et al. (2018) analysed structural accelerations of a MWT under breaking wave loads. Wave elevations over the structure were measured and pointed out that they did not consider the hydrodynamic viscous effect because of the complexity of estimating the process. The loads were measured when the wave break over the structure, omitting the wave hydrodynamics (orbital velocities) and the viscous effect of the nearby field. Shirzadeh et al. (2013) recommended an experimental study of the offshore structural behaviour under wind and waves loads and the contribution of wind damping. In addition, Esfeh and Kaynia (2020) recommended using realistic wind and wave experimental conditions for estimating structural dynamic parameters during the effect of liquefaction on the performance of monopile foundations.
Most research focussed on reducing vibrations because of fore-aft environmental loading (waves, wind earthquakes). Still, VIV may significantly affect the offshore wind turbine stability, which requires a comprehensive understanding of this turbulent loading (Jahangiri and Sun 2019). About the hydrodynamic viscous damping because of VIV, several studies have been applying numerical and physical modelling to analyse that fluid-structure interaction (Zdravkovich 1996 Their study pointed out that the soil effects reduce the structural natural frequency, but not considering the hydrodynamic viscous damping that could affect the structural dynamics. Other studies utilised the inertia coefficient (C M ) of the Morison equation and diffraction wave theory for considering the viscous damping effect of VIV during wave loading over monopile turbines (Brekke et al. 2005; Bisoi and Haldar 2014). Despite several methods for estimating the viscous terms generated by VIV, its calculation may be challenging. The coefficient selection from tables and curves may omit the inherent viscous effects that could provoke under-over estimations of the structural responses.
The literature review highlighted the improvement needed for understanding the effect of viscous damping over the structural dynamics of fixed offshore wind turbines. The applied numerical modelling does not properly consider the complex non-linear fluid-structure interactions, which generated over or underestimated structural dynamic parameters. The physical modelling use is widely reported in specialised literature using vibrating mechanisms in dry conditions (e.g., shakers, wind tunnels) for identifying structural dynamic parameters under the wind, waves, and earthquake loading. Although those experimental results showed success in their research objectives, dry mode without water may conceal critical structural frequencies. The experimental results retrieved from wave basins or flumes (wet conditions) used records of water elevation measured aside (wave loading), over the structure, or far from the near hydrodynamic field, which could omit the near-field viscous effect during the fluid-structure interaction.
To the best knowledge of the authors, there is no reported experimental research in wet conditions (wave flumes or basins) focused on the effect of extreme winds over structural dynamics during rough sea states. Despite of some researchers use physical modelling in dry conditions for assessing the extreme wind loads effects over offshore structures, that approach does not consider the real effect of viscous damping seen during experiments in wet conditions, hence, they did not compare the natural non-linear interactions between structural accelerations and near hydrodynamic field accelerations. The VIV studies showed numerical and experimental results of fluid-structure interactions, where CFD solvers with measured in situ velocities of waves and currents depicted the near hydrodynamic field. Most of these were applied for scouring and sediment transport around monopiles and analysed how structural vibrations perturbed the near hydrodynamic field. Then, it was not found studies focussed on the measured viscous effect of extreme wave and wind loading over the structural dynamics of an offshore MWT.
The present study proposes an alternative approach to determine the main structural dynamic parameters under extreme waves and wind loading in wet conditions. This study assessed the near hydrodynamic field upstream of the monopile and structural acceleration of the wind turbine for analysing the non-linear effect of wave-wind loads over the structure. In this sense, a Design of Experiment (DOE) and an Analysis of Variance (ANOVA) was performed to evaluate the response of structural damps and periods under the non-linear effects of VIV and extreme loading. This paper is organised starting with a description of the model setup, continuing with the assessment of natural periods and damping coefficients for the modeled structure. In the next stage, the fluid-structure interaction is evaluated, and the near hydrodynamic field is analysed. In the end, we offered a discussion of the results and main conclusions related to nonlinear interactions of the fluid-structure interactions seen in vortex-induced vibrations.

Methods
In order to identify natural and damped periods of the structure in dry and wet conditions, the research was developed following experimental and analytical phases: (1) Model setup., (2) Design of Experiment (DOE), (3) Natural periods, (4) Damping coefficients, (5) Fluid-structure interaction, and (6) Near Hydrodynamic field analysis. The next sections will describe and depict the application of the proposed methodology.

Model set-up
The wave channel is 37.0 m long, 0.80 m wide, and 1.20 m deep and is at the Coastal and Ports Laboratory of the Universidad Nacional Autónoma de México. A dissipative beach with gravel rocks is at the end of the channel for reducing internal wave reflections (Fig. 1). The water level for the experiment was 0.20 m. Four water level sensors recording to 100 Hz were located every 1.76 m (characteristic wavelength) from the centroid of the wind turbine model. A 3D velocity profiler (Vectrino, www.nortekgroup.com) measured x, y and z wave orbital velocities with an 80 Hz measurement rate, and was set at 0.08 m upstream from the model with a vertical recording profile of 0.07 m from the bottom. The 2D Ultrasonic Velocity Profiler (UVP) measured x and z wave orbital velocities and was set 0.16 m upstream from the wind turbine. Finally, inside the model's hub was installed a triaxial accelerometer recording at 200 Hz (Fig. 2).
This study configured a pulling system (Chakrabarti 2005) integrated by a dynameter and a nylon cable to apply extreme wind loads to the structure. The traction loads were constant during the experiment (Condition 2), and the cables pulled back the wind blades simulating constant extreme winds (Fig. 2). Because of each run of the experiment lasted 30 s approximately, the pulley system is considered proper, because this elapsed time may be related to a wind gust. The implemented pulling system simulates the effect of extreme winds such as Hurricanes or wind-gusts over the blades. Because in the real wind turbine systems, the pitch control stops the blade rotation (feathered blades) during extreme winds, the implemented pulling mechanism allowed simulating these extreme weather conditions. The scale factor (variable) for the parameters depicted in the set-up (Fig. 2) followed Froude's law recommended for the physical modelling of offshore structures (Chakrabarti 2005).
To complement the analysis of horseshoe vortex generation, the parameters D/λ was verified according to Hogben (Hogben 1978), and the Keulegan-Carpenter (KC) was calculated using the expression K C = (H s /D)π (Journée and Massie 2002). These non-dimensional parameters used the wave parameters (Hs, Tp, λ) generated by the paddle, hence drawbacks for obtaining proper parameters may occur when residuals waves remain after experimental run; 15 min were considered between each run to prevent noise in the wave measurements. The geometry and kinematic features of the model and prototype model are listed in Table 1.
The results of KC and D/λ derived from the wave parameters (Table 1) suggest that the flow regime with Hs 1 will be in the inertia regime with no wave diffraction; then, the generation of horseshoe vortex nearby the MWT is highly probable.
The DOE has two experimental conditions, each one with 3 runs with two repetitions (9 runs in total), where Condition 1 gathers the parameters for irregular wave loading (JONSWAP spectra, gamma = 1) with no wind loading, and Condition 2 considers combined irregular waves + wind loading ( Table 2).
The pulley system was disconnected during the runs of Condition 1 to avoid artificial damping The wave parameters such as significant wave height (Hs) and peak period (Tp) are related to extreme wave events.

Results and discussion
The structural response of a dynamic system is expressed in terms of its natural periods, damping coefficients and loads (fluid-structure interaction). This section discusses the findings related to this characterisation.

Natural periods
This study applied vibration tests in dry and wet conditions to identify the effect of hydrodynamic viscous damping. Five free vibration tests in dry conditions were run considering the x-dir degree of freedom (bending of the monopile in the x-direction), which is the main structural response because of the foundation restrictions and loading direction. The test comprised a controlled pulling (using a dynamometer) of the wind turbine hub with a cable and releasing it, for measuring the structural accelerations (a) as seen in Fig. 3.
The measured accelerations were utilised for calculating the natural period (Tn) and damping coefficient using the logarithmic decrement method. We performed five free vibration tests in dry conditions concluding that the structure has a natural period of 0.178 s (f = 5.6 Hz), along a damping of 0.081%.
The forced vibration test in wet conditions of this study consisted in the analysis of structural parameters of the turbine during the wave loading generated by the paddle, and wind effect produced by the pulley system. The test required a Fourier analysis for identifying the significant (mean) accelerations of loads and response accelerations derived from the measured hydrodynamic parameters and structural accelerations (x-direction), respectively. The zero spectral moment of accelerations (significant acceleration) and the associated period were selected to calculate the transfer functions and apply the half bandwidth method. Accordingly, six transfer functions were plotted using the x-dir accelerations and hydrodynamic loading measured during all the 9 runs (Table  2) by the water level sensor (SN -water surface elevation), the 3D velocity sensor (3DVS), and the 2D Ultrasonic velocity profiler (UVP); the acceleration loads were calculated using the free surface displacements measured by the SN, and the velocity records of 3DVS and UVP. The generated transfer functions have in the x-axis the load excitation frequency and in the y-axis the normalised acceleration structural acceleration load acceleration SN, Vectrino, UVP : the load excitation frequencies were calculated using the Fast Fourier Transform in the acceleration and period data.
The transfer function calculated through the SN records of the Condition 1 (Fig. 4) showed that wave load was between 0.9 Hz and 1.4 Hz, the hydrodynamic load measured by the 3DVS were 0.41 Hz and 1.6 Hz, and the load excitation frequency derived from UVP records where 0.27 Hz to 0.41 Hz. Then, the frequency of excitations in the near hydrodynamic field captured by the SN and the 3DVS were similar and doubled the frequencies captured by the UVP in the farthest field. In addition, the highest acceleration amplification was observed in the 3DVS's transfer functions, approximately 10,000 times higher than the transfer functions of SN and UVP with a peak frequency of 1.55 Hz.
The results for Condition 2 (waves + wind loading) evidenced damped curves because of the lower amplification of the normalised accelerations in the transfer curves ( Fig. 4d-f). The transfer function derived from the UVP  1 (a, b, c) and Condition 2 (d, e, f) records of Condition 2 (Fig. 4f) showed a normalised acceleration of 0.00025 with an irregular shape compared to the other transfer curves. However, because there were no reported anomalies in the UVP nor human errors, and the frequency range of the run was similar to the transfer function of Condition 1, this study considers that the high non-linear hydrodynamics generated the two amplification peaks of this transfer function (Fig. 4f). The application of the half bandwidth method through the transfer function plots (Fig. 4) revealed the damped natural periods were associated with the forced vibration test. The dynamic parameters were then calculated using the hydrodynamic loads (HS, 3DVS, UVP), as shown in Table 3.
The damped natural periods (T D ) evidenced that UVP records reported the highest values from 2.61 s for Condition 1 and 3.70 s for Condition 2. The lowest damped period was identified using the 3DVS records during the Condition 1 period (T D ), which increased four times when the wind load was applied (Condition 2). The damped periods from SN were similar to 3DVS's in Condition 1, and it was observed no significant changes in T D (0.99 s) when the wind load was applied. According to the previous results, this study considers that SN and UVP acceleration records are not affected by the near hydrodynamic field because of their higher periods and locations compared to the 3DVS (Fig. 2). It was also observed that the 3DVS's results evidenced the increment of 4-times in the damped period in Condition 2. Then, the applied wind load attenuated the structural vibrations and incremented T D to 2.70 s.

Damping coefficients
The pseudo-acceleration response spectrum calculated with the free surface and flow accelerations (Fig. 5) pointed that acceleration peaks of 3DVS (Fig. 5a-c) were narrower than SN and UVP peaks (Fig. 5d-i). The 3DVS's results evidenced two spectra peaks at 0.05 s and 0.91 s; the UVP's results also showed two peaks with a higher magnitude but less differentiated between them compared to the 3DVS's results.
About the SN's results, they showed a single peak spectra similar to a bell curve. The SN's results did not capture two spectra peaks because its location compared to the other sensors (3DVS, UVP) was not in-line with the wave direction; thus, the accelerations were not affected by wave reflection during the fluid-structure interaction (wave loading).
The spectra curves of 3DVS and UVP evidence two spectra peaks generated by the non-linear interaction of the nearby hydrodynamic field. The 3DVS's results showed lower accelerations peaks because of its closer position to the structure compared to the UVP. The proximity of the 3DVS captured the highest viscous damping, mainly because of the hydrodynamic vortex generated in the near field. Then, the farthest the hydrodynamics is measured at the foundation, the lower the hydrodynamic viscous damping is generated.

Fluid-structure interaction
Physical and hydrodynamic interactions, such as wave reflection-diffraction or changes in the near hydrodynamic field because of the generated vortex, may generate the nonlinearities during the fluid-structure interaction. In this sense, to identify if the acceleration records of SN, 3DVS, and UVP captured the non-linearities mentioned above in the wave orbital velocities, a DOE-ANOVA analysis was performed.
The results of DOE-ANOVA were plotted in standardised effects bars and Pareto charts (Figs. 6 and 7), where each sensor was labelled as follows: A (3DVS), B (SN), C (UVP). A previous inspection of normalising or standardising the input data for modifying the x-axes of Pareto charts revealed that curve shapes did not change significantly; then, this study utilised variable x-axes derived from normalising with the highest record of each dataset (Figs. 6 and 7).
The DOE-ANOVA analysis derived from runs of Condition 1 evidenced the high nonlinearity in the acceleration results. Analysing the standardised effects graphs, only one second-order standardised effect seen in run 8 (Fig. 6o) proved that the factor AA (3DVS) influenced the structural acceleration (response). In this sense, the statistical relation (cause-effect) of fluid accelerations measured by the 3DVS with the structural acceleration was evidenced.
The ANOVA second-order effect of 3DVS's accelerations (AA) pointed out that when fluid accelerated and decelerated induced to the structure vibrating in the same load direction (x-dir). Then, the Pareto chart, in most of the runs (Fig. 6), evidenced the forward-backwards bending movements because of wave orbital velocities. The Pareto chart of run 8 (Fig. 6o) mentioned above supported the standardised effect of run 8, because the three convex curves were drawn. Despite run 7 and 9 are repetitions with the same experimental Conditions of run 8, the high nonlinearity involved allowed capturing the factor effect in one of the three same runs. In addition, the results of Fig. 6 evidenced that increasing the wave energy reduced the capability of DOE-ANOVA in identifying statistical evidence of the effect of factor (loading) over the response (structural accelerations).
The DOE-ANOVA analysis for Condition 2 (Fig. 7) with combined wind and wave loading showed more statistical significance compared to the DOE-ANOVA results of Condition 1 (Fig. 7).
The 3DVS's results revealed statistical significance in 4 of the 9 runs according to the standardised effects. The ANOVA second-order effects of the 3DVS (+ AA and − AA) were identified, which support the performed statistical analysis of Condition 1 (Fig. 6) and confirm that the calculated 3DVS's accelerations influenced the structural dynamics of the structure. The + AA and −AA effects of runs 1, 2, 7, and 8 ( Fig. 7a, c, m, o) showed the effect of acceleration and deceleration of wave orbital velocities over x-dir displacements of the structure. The Pareto chart's concave and convex curves reinforced the statement that wave orbital velocities in the near hydrodynamic field cause a high non-linear fluid-structure interaction.
Comparing the DOE-ANOVA results of Condition 1 and Condition 2 revealed that wind loading damped the structural dynamics because of the constant loading in the x-direction. The tilt of the wind turbine because of the aerodynamic loading controlled the x-dir displacements what reduced the normalised accelerations of the transfer functions (Fig. 4). Then, extreme wind loading revealed additional structural dynamic parameters that must be identified for dampers designing.

Near hydrodynamic field analysis
To verify the warning of generation of horseshoe vortex nearby the MWT (Table 1) and to analyse the non-linear effect of wave loading over the structure, were inspected the velocity records of Run 1-Condition 1 ( Table 2) measured by the 3DVS (Fig. 8).
The measured horizontal and vertical velocities by the 3DVS of each 1 cm layer of a 7 cm profile (Fig. 8) allowed analysing the near hydrodynamic field (Fig. 8). The hydrodynamic field's initial records were close to 0 m/s because the experiment started with still water. From the 1000 records (10 s), the hydrodynamic field is excited, and from 1500 records, the x velocities showed positive and negative values associated with the in-line wave orbital velocities. The velocities from record 1500 to 4500 (30 s of elapsed time) pointed velocities about 0.2 m/s (negative and positive), where the velocities depicted a vertical vortex because of the x and z velocities.
In order to identify the mentioned above vortex because of the fluid-structure interaction, which induces viscous damping to the structure, this study zoomed in to the velocity measurements (Fig. 8) and extracted 0.3 s (26 records) of the time series (Fig. 9). As a result, positive horizontal xvelocities were detected from 16.4 to 16.7 s and negative values from 16.7 to 17.0 s, and positive values of vertical z-velocities in all the elapsed time. Then, horizontal x and z vertical component velocities evidenced an elliptical vortex (horseshoe vortex) with rotation, as seen in Fig. 9, which was warned by the KC and D/ λ results of Table 1.
The spectral analysis of the velocity components measured by the 3DVS sensor (Fig. 10) (Fig. 9), the spectral analysis (Fig. 10) and the calculated St, KC and D/λ, a horseshoe vortex was identified governing the near hydrodynamic field of the MWT in this study.
Comparing the damped natural periods of Condition 1 (3DVS) reported in Table 3 against the velocity profiles measured nearby to the structure (Figs. 8 and 9), the spectral curve of 3D velocities (Fig. 10) and the vortex's period mentioned above (Fig. 9), it is evidenced that wind load increased the  3DVS (a, b, c), UVP (d, e, f) and SN (g, h, i) measurements. Free surface accelerations were calculated using the water level records (SN), and flow accelerations were calculated using the 3DVS and UVP records structural damping and the vortex controlled the structural periods. As a result, the viscous damping generated at the near hydrodynamic field attenuated the hydrodynamic load over the structure. The natural periods of the structure in dry Conditions (Fig. 3) were 0.017 s and damped to 0.62 s (Condition 1-3DVS) during forced vibrations in wet Conditions ( Table 3). The increment of the natural period occurred by the wind load effect, and by the viscous effect of the nearby elliptical vortex (Fig. 9) upstream of the structure, which controlled the structural vibration and damped the structure from 0.08% (dry conditions) to 5.00% (wet conditions) ( Table 3). The increment of almost 5.00% of the damping coefficient of this study agrees with the relative differences between excitation frequencies and the natural frequency known as the over-5.00% criterion reported by the research of Ko Y.Y. (Ko 2020) and the guidelines of DNV-GL (DNV-GL 2016). The calculated 0.05% damping coefficient of Condition 1 (3DVS) The normalised accelerations of the transfer functions (Fig. 4) evidenced that 3DVS's transfer curves reported the highest structural amplification, and the DOE-ANOVA analysis (Figs. 6 and 7) pointed out the statistical significance of the effect of 3DVS's accelerations over the structural accelerations (responses). As a result, it was proved the relevance of measuring at the near hydrodynamic field in front of the structure through the 3DVS. In this sense, the vortex generated by the fluid-structure interaction (Fig. 9) could not be captured by the water level sensor (SN), nor the 2D Ultrasonic velocity profiler (UVP) located farther from the structure (Fig. 2). For preliminary stages of the structure, designing may use Pseudo-acceleration spectra curves (Fig. 5), where the fluid-structure interaction is not required for identifying damping coefficients. However, these spectra curves suggest a range of acceleration responses where the critical period (amplification-resonance) could be wrongly identified. However, these spectra curves may be considered for estimating the critical range of periods where the structure requires effective damping, which can be validated by the frequency ranges of the transfer curves and the calculated structural damped period through the application of the Half Bandwidth method in wet conditions.   Previous research demonstrated the complexity of VIV (Tödter et al. 2021). Hence, the authors highlighted that research results strongly depend on structural parameters such as damping, geometry and structural parameters, which provoke different response amplitudes, frequencies and vortex patterns. In this sense, each OWT configuration will generate different structural parameters results. However, the damping ratio and Strouhal-KC numbers may be considered generalized parameters for structure design, no matter the ratio between models and prototypes. These parameters will keep the same, considering their non-dimensionality and defined limits. In this sense, the damping ratio calculated in this study applies to any scale between model and prototype.
In the common practice of offshore structure designing, it is assumed that structural damping may be determined in dry Conditions (e.g., shakers or vibrating tables) or through the mere utilisation of wave loads or measured velocity profiles Fig. 9 Zoom in to the 3D velocity profiles measured by the 3DVS for the Run 1-Condition 1. The bottom level is pointed by the 0 cm of elevation (the lowest point of the profile), and the 7 cm level is the top of the profile without the near fluid-structure interaction. It is then necessary for the offshore structure design to validate damped periods in wet Conditions and extreme excitations such as hurricane winds. The VIV and strong wind loading significantly modify the damped natural period.

Conclusions
The non-linear analysis of near-field fluid-structure interaction was performed after applying five steps or phases: (1) Model set-up, (2) Natural periods, (3) Damping coefficients, (4) Fluid-structure interaction, and (5) Near Hydrodynamic field analysis. As a result, each step revealed significant findings of the relevant effect of vortex-induced vibrations (VIV) over the structural dynamics of an offshore MWT.
The estimation of structural natural periods in dry conditions showed underestimations because of the calculated lower periods and damping coefficients compared to the dynamic parameters calculated in wet conditions.
The transfer function plots showed that normalised accelerations calculated using 3DVS's measurements were about 10,000 times higher compared to the transfer functions of SN and UVP. During irregular wave-wind loading, the transfer functions from 3DVS's records evidenced the increment of 4-times of the structural period compared to the transfer functions during solo irregular wave loading. It was evidenced that wind loading controlled the x-dir structural acceleration, damped the natural periods, and increased the damping coefficient.
The pseudo-acceleration response spectrums calculated through the water level sensor showed a single peak instead of two peaks. As a result, it was evidenced that utilising water level records measured aside from the structure will not capture the non-linearities of the near-field fluid-structure interactions, which will generate single-peaks spectra curves omitting different structural modes.
The DOE-ANOVA analysis confirms that the 3DVS's measurements captured VIV's non-linear effects over the structure, where a vertical vortex near the monopile foundation was detected. The identified vortex controlled the structural vibration during irregular wave loading and generated a damping coefficient of 0.05 according to the over-5.00% criterion, vortex which got intensified during the wind loading effect. The calculated damping ratios of this research and the identified effect of VIV over the structural dynamics, may be applicable to any scale between models and prototypes due to their non-dimensionality and defined limits.
The model scale used in this research may conceal excitation frequencies seen in real-life environments because of currents, then, the derived results of the study are limited to extreme sea states when the potential energy of waves heights predominates over the kinetic energy of ocean currents. For future research, it is suggested to apply this methodology combining waves-wind-current loading in waves flumes to identify how currents affect the fluid-structure non-linear interactions.
Finally, this study recommends identifying dynamic parameters of offshore monopile foundations in wet conditions and measuring 3D velocities of the near hydrodynamic field to capture non-linear interactions of the fluid-structure interactions seen in vortex-induced vibrations. advice, support, and recommendations during the experimental activities.
Funding Open Access funding provided by Colombia Consortium.
Data availability All data generated or analysed during this study are included in this published article.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.