Iterative Frequency-Domain Response of Floating Offshore Wind Turbines with Parametric Drag

Methods for coupled aero-hydro-servo-elastic time-domain simulations of Floating Offshore Wind Turbines (FOWTs) have been successfully developed. One of the present challenges is a realistic approximation of the viscous drag of the wetted members of the floating platform. This paper presents a method for an iterative response calculation with a reduced-order frequency-domain model. It has heave plate drag coefficients, which are parameterized functions of literature data. The reduced-order model does not represent more than the most relevant effects on the FOWT system dynamics. It includes first-order and second-order wave forces, coupled with the wind turbine structural dynamics, aerodynamics and control system dynamics. So far, the viscous drag coefficients are usually defined as constants, independent of the load cases. With the computationally efficient frequency-domain model, it is possible to iterate the drag, such that it fits to the obtained amplitudes of oscillation of the different members. The results show that the drag coefficients vary significantly across operational load conditions. The viscous drag coefficients converge quickly and the method is applicable for concept-level design studies of FOWTs with load case-dependent drag.


Introduction
Research on floating instead of fixed-bottom offshore wind energy started more than ten years ago.Numerous demonstration projects with commercial-scale wind turbines have proven the feasibility of harvesting offshore wind energy at water depths beyond the limits of substructures rigidly connected to the seabed [1].

Coupled Aero-Hydro-Servo-Elastic Dynamics
For a prediction of the dynamic behavior of a Floating Offshore Wind Turbine (FOWT), a multidisciplinary analysis is necessary.There are different coupling effects between the submodels of aerodynamics, hydrodynamics, servodynamics and structural dynamics.This results from the soft substructure restoring forces, especially of semi-submersible concepts, leading to a slow fore-aft motion of the rotor and therefore an influence on the aerodynamics.One important example for the dynamic coupling effects is the influence of the wind turbine blade pitch controller on the global fore-aft mode of the FOWT.The controller pitches the blades, due to a change of the relative wind speed seen by the rotor, originated from a fore-aft oscillation of the platform.This pitching reduces the thrust and can eventually yield a fore-aft system instability as shown by [2][3][4], among others.The effect was experimentally shown for the semi-submersible FOWT concept of this work in [5].Several open-source and commercial integrated FOWT models, respecting these couplings, have been developed.The first was the FAST model by National Renewable Energy Laboratory, Boulder, CO, USA (NREL) (see [6]).
Higher-fidelity models are necessary to model the impact the fore-aft motion of the rotor has on the inflow, the wake and the local flow around the airfoils.These unsteady aerodynamic effects have been investigated first by [7,8].Experimental approaches followed by [9,10].Various high-fidelity simulations have been performed, i.e., [11][12][13][14][15][16][17].They show that these unsteady aerodynamic effects of the wake depend largely on the frequency and the amplitude of the fore-aft motion but also on the Tip Speed Ratio (TSR) of the rotor.
The common methodology for hydrodynamic FOWT simulations is a time-domain formulation, which respects the frequency-dependence of the hydrodynamic coefficients, derived from panel codes [18,19].In the time-domain, the frequency-dependence results in a convolution integral, representing fluid memory effects (see Cummins' equation [20]).Viscous effects are usually accounted for through the inclusion of the quadratic drag term of Morison's equation [21].This drag results on the one side in an excitation force, due to the wave particle motion, next to the Froude-Krylov and diffraction pressures from the panel code computation.On the other side, Morison drag results in a viscous damping of the floating body.Simplified models have been published neglecting the excitation part, modeling only the damping [22].

Hydrodynamic Viscous Drag
A common practice for the determination of the damping properties is to identify a linear damping matrix or quadratic generalized platform damping from model tests, such as in [22,23].Another way is to identify the member-based Morison drag coefficients, instead of the lumped, generalized damping coefficients, see [24,25] or [26].Commonly, free-decay tests are used to obtain these drag or damping coefficients from model tests.Computational Fluid Dynamics (CFD) simulations have also been used for the identification of FOWT drag [27].
As for the load case simulations in the seminal Offshore Code Comparison Collaboration (OC3) through Offshore Code Comparison Continuation, Continued, with Correlation (OC5) projects, the damping properties are usually estimated based on a calculation of the Reynolds number for different sea-states and a fixed platform.The obtained drag coefficients are then used equally for all load conditions in the participating tools.
The present study addresses the question of the load-case dependence of the drag coefficients for FOWTs.As we showed experimentally in [28], the viscous damping properties of FOWTs can vary significantly, depending on the load case.This is due to the strong Reynolds dependence (Re = vD/ν, with velocity amplitude v, characteristic length D and kinematic viscosity ν) and also a dependence on the Keulegan-Carpenter number (KC) (KC = vT/D, with period of oscillation T) of the drag coefficients, see ( [29], Chapter 7) and, additionally, the body surface roughness [30].For circular cylinders, the Reynolds and KC-dependence of the drag is given by Sumer and Fredsøe ( [31], p. 144).It can be seen in that reference that the drag is a complex, not monotonous function of both nondimensional numbers.Therefore, generally reasonable values, only dependent on the roughness, are given in ( [30], p. 117).For plates with sharp edges (like heave plates of semi-submersible FOWTs), it is different: here, flow separation is almost always present and the drag is independent of the Reynolds number.Thus, a parameterization is well possible.In [32], a monotonically decaying drag as function of KC is identified.KC is a function of the body diameter, the motion velocity amplitude and the period of oscillation.Thus, the heave plate drag is a function of the dynamic response of the FOWT.The drag values of [32] could be confirmed for FOWT model tests in [28] for the public TripleSpar design [33], which will be used in this work.
This dependence of the drag on the motion response makes it difficult to determine their magnitude a priori.The response would have to be obtained with an initial guess of the drag coefficients and be updated subsequently, based on the obtained response.Simulation times of common coupled FOWT models are approximately half the simulated time, which yields a large computational effort with such an iterative procedure.
Linear frequency-domain methods are therefore more appropriate for an iterative solution of the drag coefficients.Methods for the linearization of the quadratic drag are available, i.e., [34].A simplified approach for FOWTs was presented recently by Pegalajar-Jurado et al. [35].It accounts for the relative velocity between the body and the fluid.A more complex method for an iterative drag linearization for offshore structures in the frequency-domain was published by Langley in [36].It was applied to FOWTs by Lupton [37].The KC-dependence is not a subject of the latter publication as they focus on the linearization itself.The influence of the wind turbine aerodynamics on the response was also not considered.
The motivation of this work is to include all relevant aero-hydro-servo-elastic effects in a coupled dynamic frequency-domain model for an iterative solution of the response with a parameterized heave plate drag as function of KC.For this end, a generalized model for the member-based Morison viscous drag coefficients will be defined in Section 2. Subsequently, the linearized drag is obtained from the quadratic viscous drag.A reduced-order coupled model will be used to obtain the global system response.The results of Section 3 show the iteratively obtained heave plate drag coefficients and the global response predicted by the reduced-order model together with the one predicted by FAST [38], followed by the conclusions in Section 4.

Numerical Simulation Methods
This section starts with a description of the reference FOWT model used in this work.The hydrodynamic viscous drag model and its linearization will be addressed afterwards and subsequently its integration into a coupled low-order simulation model.

Reference Design
The FOWT design "TripleSpar" is shown in Figure 1.The left part of the figure shows an image of the scaled model of [39], which will be used in this paper, next to the full-scale prototype.The right part of Figure 1 shows the coupled simulation model, which will be introduced in Section 2.
The design was initially developed in the project INNWIND.EU [40].It is a concept design for the research community and can be downloaded [41].The main parameters can be found in Table 1.It is a deep-drafted semi-submersible with a large portion of the hydrostatic restoring coming from the low center of mass.While the columns are made from pre-stressed concrete, the thin heave plates at the lower ends of the columns are made out of steel.Although it has slender columns, it has a lower draft than a spar, due to the three columns.Spars and semi-submersibles with heave plates are common FOWT substructures today [1], which makes this concept well suited for the presented research.
The experimental model was tested in a scale of 1:60 in 2016 at the Danish Hydraulic Institute (DHI) with results presented in [5,28,39].The test program included irregular wave with wind Load Cases (LCs) with different significant wave heights H s and peak spectral periods T p .Two of these cases are considered in the present work, LC 7 (H s = 5.5 m and T p = 8.4 s) and LC 9 (H s = 9.5 m and T p = 11.1 s).

Generalized Morison's Equation
Morison's original formulation [21] was extended for the present work to account, firstly, for the relative fluid velocity due to the body motion and secondly, for heave plates with force components in vertical direction (in the body-fixed frame).The substructure is discretized through nodes as shown in Figure 2.For a flexible model setup, each node k has modified Morison coefficients C * D,i for all three directions i with the area A ik , projected on direction i in body-fixed coordinates (attention needs to be paid to vertical drag coefficients for heave plates C D,z = C D,hp .Here, the drag force is calculated as F D,z = C * D,z v 2 z , whereas, in HydroDyn [42], it is assumed that for each heave plate a node is defined for both ends and half the drag force applies to each end).The viscous drag term of Morison's equation is with the modified hydrodynamic coefficients The nodal body velocities v b,ik are a function of the generalized coordinates, shown in Figure 1.The wave particle velocities v w,ik can be calculated in the frequency-domain using deepwater approximation with linear potential flow theory [18], giving v w,x (ω, z) = ζ(ω)e kz with the wave number k, the upward vertical coordinate z at Still Water Level (SWL) and the wave height amplitude spectrum ζ(ω).The substructure has transverse drag forces in horizontal direction only for the slender columns, not for the heave plates.

Linearized Viscous Drag Model
A linearized description of the drag is necessary for the iterative frequency-domain solution.
Morison's drag term, Equation ( 1), contains both damping and excitation forces due to the relative velocity between the undisturbed fluid and the body.The next sections address first the method to obtain the linearized drag coefficients.Thereafter, both the calculation of the excitation part through a transfer function from the wave height to the nodal forces follows as well as the nodal damping of the FOWT platform will be addressed.

Drag Linearization
Possible linearization methods for Equation ( 1) are harmonic or stochastic linearization.In the case of harmonic linearization, the error due to neglected higher-order terms is minimized for a sinusoidal velocity time series.Stochastic linearization methods improve the linearized model quality for stochastic velocity time series.A discussion with a comparison of alternatives can be found in [43].
The stochastic linearization used here follows Borgman [34].It is based on the Standard Deviation (STD) of the relative fluid velocity.In this work, an efficient method for the linearization of the velocity-dependent part of Equation ( 1) is chosen: the nonlinear drag coefficient C D,ik is split into one part dependent on the wave-velocity and a part dependent on the body-velocity.A comparable approach was selected by Pegalajar-Jurado et al. [35].The linearized nodal Morison force FD,ik results with the approximation as in direction i.The part representing the external drag-excitation is ( Cw * D,ik v w,ik ), and the damping part is (− Cb * D,ik v b,ik ).Only the latter part was considered in [35], while the former was neglected.The approach selected here with Equation (3) neglects the cross-correlation between fluid and body velocity.However, this cross-correlation can be considered when calculating the linearized counterpart CD,ik of the nonlinear drag coefficients C D,ik .
The linearized drag coefficients CD,ik for wave and body velocities are, following [34], a function of the STD σ(v ik ) of the velocity in direction i of a node k with characteristic cross-sectional area A ik , normal to direction i as A challenge for the implementation of Equation ( 4) is related to the established coupled simulation method, described in the introduction.Commonly, the first-order wave forces and also the slow-drift forces are assumed to be independent of the body states.They are therefore calculated a priori and the forces are an external input to the simulation.For an iterative solution of Equation ( 4), this a priori calculation is not possible anymore.A modification of Equation (4) will therefore be made.

Drag Excitation
The external Morison excitation force can be written as a transfer function, when linearized.The goal is to obtain these transfer functions for the generalized platform coordinates ξ (for the low-order coupled model, subject of Section 2.4, they are reduced to ξ = [x p , z p , β p ] T ).With the wave particle velocities v w (ω, x, z), the transfer functions from the wave height ζ 0 at the reference position to the three translational (superscript t) forces are For the rotational direction (superscript r) from ζ 0 to the three generalized moments on the platform about the platform reference point, here the Center of Flotation (CF) holds It is the cross product with the position vector r k to each node k.Equations ( 5) and ( 6) can be unified with the transformation tensor for the nodes H k , assuming here a rigid platform, as The cross-product is written with the operator S (a) b = a × b and I denotes the identity matrix.Then, results. Figure 3 shows the obtained linear transfer functions for the TripleSpar.The calculation of the wave kinematics at the longitudinal column position results in the observed interference effect.

Sfrag replacements frequency [Hz]
∠G 55 Fmor (ω) In order to keep the above-mentioned a priori calculation of the exciting forces, the generalized body velocity response spectrum for the calculation of v b,ik , Equation ( 3) is calculated from the panel code Response Amplitude Operator (RAO) ξ(ω)/ζ(ω) and not from the coupled model of Section 2.4.This is convenient because the panel code calculation is performed as pre-processing of the coupled simulation.The RAO is differentiated to obtain the generalized velocity transfer function The cross-spectral density matrix of the generalized platform velocities is with the complex conjugate transpose G * T .The velocity spectra S v k v k (ω) of the nodes k of Figure 2 are obtained with the cross-spectral density matrix S ξ ξ (ω) of the reference point.The nodal velocity cross-spectral density matrices As stated, the response spectra S ξ ξ result from the RAO, independent of the coupled system response.This approach yields a significant simplification of the computational implementation, while considering the relative fluid velocities instead of the absolute ones.The approach mainly respects relatively well the body response from first-order wave forces.It neglects, however, the low-frequency response at the platform eigenfrequencies and all wind turbine-related effects.This simplification can be seen as a separation of frequency ranges: the range of the platform eigenfrequencies and the one of the first-order waves.The discussed drag excitation forces act at the first-order wave frequencies, above the platform eigenfrequencies.Thus, neglecting the response at the platform eigenfrequencies is consistent, when looking only at the first-order wave frequencies.It is noted, however, that this frequency-domain perception is a simplification of the complex nonlinear flow phenomena in the time-domain.
For the TripleSpar design (Figure 1), it was found that the velocity-induced excitation of Equation ( 1) is small, about 20 times smaller than the acceleration-dependent Froude-Krylov forcing.The results of this study will be published in the PhD dissertation related to this work.

Drag Damping
In this section, the drag coefficients of Equation (3), related to the damping Cb * D,ik , will be identified.The damping is part of the system properties: in a unified state-space description of the FOWT, it is part of the system matrix (A) and influences the eigendynamics.These eigendynamics are especially important for the prediction of the system stability and for the controller design, as mentioned in the Introduction.Consequently, this coupling should be taken into account when determining the drag.To do this, an iteration is necessary: the hydrodynamic damping is updated as a function of the entire system response of the coupled model, until the calculation converges.
Figure 4 (left) shows the velocity STD over the 1:60 scaled TripleSpar columns of [39] for two sea states.The body velocity depends on the depth, in a way determined by the inherent wave response properties of the floater.The drag coefficient, Figure 4 (right), calculated with Equation (4), depends in the same way on the body velocity.The linearized nodal drag force FD,x1 for the upper node at SWL for one column is shown in Figure 5, together with the nonlinear drag force resulting from the nodal velocity v b,x1 .The figure shows FD,x1 for the same two environmental conditions as Figure 4.The linearization depends clearly on the response magnitude (STD, marked by dashed lines).
The calculation of the linearized drag coefficients Cb * D,ik with Equation (4) for the damping problem is again simplified.In order to separate the a priori calculation of external forces, which is based on the undisturbed wave kinematics v w,ik , from the simulation, the undisturbed wave kinematics are neglected for the damping problem.This approach is consistent with the above-mentioned separation of frequencies: the low-frequency response, for which the damping is of critical influence, is outside the first-order wave frequencies.The obtained linear nodal drag coefficients Cb * D,ik can be integrated giving a generalized coefficient matrix D ∈ R (6×6) such that the generalized linearized Morison drag force Fb D on the platform CF can be written as This is a very straightforward description, which can be used across models, i.e., to update the RAO with the viscous drag forces.The integration is again written using the transformation tensor H(r k ) of Equation ( 7) with the damping matrix of the nodes The iteration and the obtained damping matrix D are shown in Figure 6 for the scaled TripleSpar.The solution converges already after about four iterations.
In order to assess the quality of the linearization of the damping problem, the overall Morison drag forces F b D (t) were calculated for an LC of the scaled experiment [39].The time series were calculated once with the quadratic drag force (Equation ( 1)) and once with the linearized drag force (Equation ( 3)).For a clear comparison, only the damping force is shown and the drag excitation part is set to zero.The left part of Figure 7 shows the generalized platform velocities.The right side shows the integrated generalized forces in the respective coordinates.The linearized drag is well in line with the nonlinear one, except for the larger peaks of the velocity.Here, the nonlinear model reaches larger values.This completes the linear Morison model.The next section will introduce the low-order coupled simulation model.

Low-Order Aero-Hydro-Servo-Elastic Model
A mechanical sketch of the low-order coupled simulation model, which is used in this work, was shown in Figure 1.The model has been developed with the focus of a computationally efficient representation of the overall dynamics of the FOWT for control design and an overview on critical load situations.It was used in comparable configurations in previous works [4,[44][45][46].Its derivation starts with physical (as opposed to black-box) models, keeping the basic nonlinearities.In a next step, the nonlinear models are linearized.A frequency-domain computation of the response spectra will be used in the present work.The model will be called Simplified Low-Order Wind Turbine (SLOW) in the following.
The structural model is based on a flexible multibody formulation with rigid and flexible bodies according to [47,48].For this work, the coupled equations of motion have been set up for six Degrees of Freedom (DoF) in 2D: the floating platform rigid body in surge (x p ), heave (z p ) and pitch (β p ), a flexible tower with one representative coordinate x t , the rotor speed Ω and the blade pitch actuator θ 1 (see Figure 1).These DoFs are selected in order to reduce the system to the most important effects, which are critical for the main system dynamics.These act usually in the direction of wind and waves (x).The main simplifications of the aerodynamics and hydrodynamics significantly accelerate the model execution, even in frequency-domain.
The hydrodynamic model includes first-order panel code coefficients and second-order slow-drift forces through Newman's approximation.As in most FOWT simulators, wave drift damping, in the sense of ( [30], p. 153), is not included.The radiation problem has been simplified through a frequency-independent added mass as in [49].Radiation damping is completely neglected because it has negligible effects for the TripleSpar.

Results
The previously developed methods for modeling both the linearized wave excitation forces and the linearized damping forces from drag will now be coupled with a parametric function for the heave plate drag.The results are shown for environmental conditions used in the EU-project LIFES50+ for fatigue assessments ([50], Chapter 7).They are listed in Table 2.

Parametric Drag
The parameterization of the column drag coefficients C D,ik is difficult due to their dependence on both Reynolds and KC-number, as described in the Introduction.In [28], the heave plate drag C D,hp was identified for coupled model tests of the TripleSpar, including the wind turbine and controller.The identified drag compares well with the values published by Tao et al. ( [32], p. 1012).The authors of that publication identified the heave plate drag as a function of KC as shown in Figure 8.
For this work, the values of Figure 8 were parameterized for an online interpolation, based on the instantaneous KC.Although the column drag C D,ik is not parameterized, the parameterization of the heave plate drag C D,hp improves the representation of the global FOWT dynamics.This is because the heave plate drag influences strongly the fore-aft dynamics of the rotor, which induce the couplings with aerodynamics and the controller, as described in the Introduction.

Solution Procedure
The computational implementation of the iterative hydrodynamic model is shown, as a summary of the previous sections, for the drag excitation forces and the damping forces in Figure 9.It visualizes the simplified representation of the relative fluid velocity between the wave particles and the body of Section 2.3.The procedure for the excitation problem focuses on the first-order wave frequencies, while the procedure for the damping problem focuses on the lower frequencies, at the platform eigenfrequencies.The linearization is the same for the columns and the heave plates.The nonlinear drag coefficients of the heave plates, however, are parameterized and change, dependent on the member response magnitude.

Iterative Solution
The procedure of Figure 9 will be executed iteratively with C D,hp = C D,hp (KC).In order to obtain the velocity amplitude v for the estimation of KC, the STD of the stochastic response was used with a constant amplification factor.Assuming a Rayleigh distribution of the maximum amplitudes of the nodal velocity response over the columns, it is possible to approximate a representative velocity amplitude v.This requires a stationary Gaussian distribution of the platform velocity response, which is a good assumption for wave-induced motions [51].Here, the mean of the third highest response amplitudes is taken as a reference, the same as for the calculation of the significant wave height H s .The characteristic velocity v is then The iteration was performed with the simulation model of Section 2.4.The computation was done for a representative load case at above rated winds of Table 2 with the full-scale TripleSpar concept.The resulting generalized damping matrix D over the wind speed bins is shown in Figure 10.A clear dependence of the linearized damping on the environmental bins is observable.It is noted that the wind speeds are correlated with higher sea-states.2.
One effect driving the observed characteristic is that the linearized drag Cb * D,ik increases for larger velocity response amplitudes as shown in Figure 5.The other effect is the KC-dependence for heave plates: Figure 8 shows that the drag coefficient decreases for larger KC-numbers, or larger sea states (given the same periods are present).Thus, two opposing effects drive the sensitivity of Cb * D,hp with respect to the velocity amplitude.
The main conclusion drawn from Figure 10 is that the overall generalized damping depends significantly on the sea state, or the load case.This means that neglecting the load case-dependence of the damping might yield an over-conservative estimation of the stability for larger sea-states and likely a too optimistic stability at small sea-states.
Figure 11 shows KC as well as the heave plate drag over all operational wind speeds.Both drag coefficients are shown, the one of the damping problem (blue) and the one of the excitation problem (orange).It can be seen that their qualitative shape is equal, although differences occur for the higher sea-states (see Table 2).The origin of the peak at v = 22 m/s was analyzed.It is not caused by the interaction with the controller.Therefore, it is a result of the system property sensitivities with respect to the operating point and the wind and wave conditions of Table 2.The drag coefficient for the heave plates C D,hp varies in the shown load conditions significantly in the range C D,hp = 11 . . .18.In the Offshore Code Comparison Collaboration, Continued (OC4) code comparison project [23], a value of C D,hp = 9.6 was used for all simulations.

Verification
With the developed iterative procedure, based on the linear frequency-domain SLOW model, the nonlinear drag coefficients of Figure 11 are available as a function of the environmental conditions.These values for C D,hp have been used for the nonlinear FOWT design tool FAST, version 8.16, with all 25 DoFs (including collective blade pitch actuator) enabled.Figure 12 shows the results for an operational case above rated wind.In spite of the reduced number of DoFs of the SLOW model, the main response in surge x p , pitch β p and tower-top bending x t is well reproduced.This holds for the low-frequency response, but also for the first-order wave response.The results show that the low-order model is capable of predicting viable values for the nonlinear heave plate drag C D,hp as a function of the KC-number of the body-response.Nonlinear effects are small enough for a good representation through a frequency-domain model.This does not hold for the nonlinear switching of the wind turbine controller around rated winds and also for the rotor azimuth-dependent blade excitation from wind speeds.The latter effect can be seen in the tower bending x t at a frequency around 0.43 Hz.

Conclusions
A new method has been developed to iterate the heave plate drag of a FOWT, parameterized as a function of KC, using a low-order aero-hydro-servo-elastic frequency-domain simulation model.The heave plate drag is thus calculated taking into account the FOWT coupled response including the wind turbine controller and excitations from wind and waves.The KC-number is calculated based on this coupled response for each heave plate node.
The procedure involves an efficient representation of the relative fluid between water particles and the body.A separation was made between the drag, responsible for damping and the one responsible for excitation.Despite these simplifications, the comparison of the time series of the drag force has shown a reasonable agreement between the nonlinear and the linearized model.
The employed frequency-domain simulation model is capable of representing the platform motion response, compared to the state-of-the-art engineering model FAST.Thus, the computationally expensive multidisciplinary FOWT system can be reduced to allow for an iterative solution of the viscous drag.
The results show that the heave plate drag varies strongly over realistic operational conditions.A constant drag, which is the common practice for load case simulations, can therefore severely over or underpredict the damping.This is especially important for the fore-aft direction, as the interaction with the wind turbine controller can even further reduce the global damping in this direction.Thus, the proposed method is a realistic means to obtain individual heave plate drag coefficients for each load condition.These coefficients can then be used in higher-fidelity models for a more realistic prediction of the extreme and fatigue loads of FOWTs.A significant advantage of the presented approach is that it is computationally efficient and therefore suitable for concept-level simulations.At this stage, experimental tests or expensive CFD-calculations to obtain the drag coefficients are not yet available.
This study is expected to set off a discussion on how to model the excitation drag, the viscous drag at low-frequencies and possibly the wave drift damping for FOWTs in higher-fidelity models realistically but also approximately in simplified models.

Figure 2 .
Figure 2. Discretized platform shape with nodes of horizontal and vertical drag coefficients.In this work, horizontal drag forces are only applied to the column nodes (dark blue) and vertical drag forces are only applied to the heave plate nodes (light orange).

Figure 3 .
Figure 3. Morison excitation transfer functions G Fmor of Equation (8) with magnitude (left) and phase (right) for drag-induced component only with C D = 0.6 and C D,hp = 20 for scaled TripleSpar design of[39].Top: surge (dark blue) and heave (light orange), bottom: pitch.

Figure 6 .
Figure 6.(a) iterative solution of Morison drag coefficient matrix D for all relevant indices i, j for aligned wind and waves (scaled model of [39], LC 7) and (b) obtained generalized damping matrix.

Figure 7 .
Figure 7. Overall Morison drag force time series F b D (t), nonlinear (dark blue) and linearized (light orange) for scaled TripleSpar in LC 9 of [39].

Figure 9 .
Figure 9. Overview on calculation scheme for Morison drag excitation as well as damping forces as input to linear frequency-domain model.The left column shows the excitation problem, the right column shows the damping problem.The presented approach has distinct procedures for the wave kinematics and the body kinematics.Parametric drag is considered for heave plates.

Figure 10 .
Figure 10.Global linearized Morison drag matrix elements for full-scale TripleSpar over environmental conditions of Table2.

Figure 11 .
Figure 11.Heave plate KC-number and resulting quadratic drag coefficient C D,hp from data of Figure 8 for environmental conditions of Table 2. Case (1): Damping (blue) with Closed Loop (CL) system model iterated until convergence.Case (2): Excitation (orange) with body velocity calculated based on panel code RAO and wave kinematics with linear wave theory.

Figure 12 .
Figure 12.Model verification: Power Spectral Density (PSD) of response to LIFES50+ operational LC of Table 2 at wind speed v = 17.9 m/s (H s = 4.3 m, T p = 10 s).Linearized SLOW model (dark blue), FAST (light orange).

Table 2
(1)ase(1): Damping (blue) with Closed Loop (CL) system model iterated until convergence.Case (2): Excitation (orange) with body velocity calculated based on panel code RAO and wave kinematics with linear wave theory.