BlockWerk Guides

Floating Offshore Wind Turbine: Active Wave Damping

Floating Offshore Wind Turbines (FOWTs) represent the frontier of renewable energy in deep coastal waters. Unlike fixed-bottom turbines anchored directly to the seabed in shallow waters (< 50 m), floating foundations such as the NREL OC3-Hywind spar-buoy operate at depths of 100 to 300+ meters.

While deep waters offer stronger and steadier wind resources, they introduce severe structural dynamics challenges. A floating spar-buoy behaves in pitch (rotation θ\theta) as a massive submerged pendulum. Because radiation damping from surrounding seawater is exceptionally small (ζ≈0.03\zeta \approx 0.03), ocean wave excitation at or near resonant frequencies builds compounding rotational oscillations that threaten platform stability and structural life.

This case study investigates active wave damping using a reduced-order 1-DOF spar-buoy pitch model, demonstrating why naive feedback fails and how a second-order notch filter solves the fundamental frequency conflict.

The Physical Setup

The platform is based on the benchmark NREL OC3-Hywind spar-buoy concept. In pitch rotation (mode 5 in naval architecture), the dynamics are governed by Newton-Euler rotational equilibrium:

Itotθ¨+B55θ˙+C55θ=Mwave−MctrlI_{\text{tot}} \ddot{\theta} + B_{55} \dot{\theta} + C_{55} \theta = M_{\text{wave}} - M_{\text{ctrl}}

ParameterMeaningBenchmark Value
I55I_{55}Structural pitch moment of inertia5.00×1010 kg⋅m25.00 \times 10^{10}\ \text{kg}\cdot\text{m}^2
A55A_{55}Hydrodynamic added mass (entrained water)1.70×1010 kg⋅m21.70 \times 10^{10}\ \text{kg}\cdot\text{m}^2
ItotI_{\text{tot}}Total virtual pitch inertia (I55+A55I_{55} + A_{55})6.70×1010 kg⋅m26.70 \times 10^{10}\ \text{kg}\cdot\text{m}^2
B55B_{55}Hydrodynamic radiation damping8.65×108 N⋅m⋅s/rad8.65 \times 10^8\ \text{N}\cdot\text{m}\cdot\text{s/rad}
C55C_{55}Hydrostatic restoring stiffness3.111×109 N⋅m/rad3.111 \times 10^9\ \text{N}\cdot\text{m/rad}
fnf_nNatural pitch resonance frequency0.0343 Hz (ωn=0.2155 rad/s, Tn≈29.2 s)0.0343\ \text{Hz}\ (\omega_n = 0.2155\ \text{rad/s},\ T_n \approx 29.2\ \text{s})
ζ\zetaOpen-loop damping ratio≈0.03\approx 0.03 (Q≈16.7Q \approx 16.7, +24.5 dB+24.5\ \text{dB} resonance peak)

The Degrees of Freedom ("55" Convention)

In naval architecture and marine hydrodynamics, any floating vessel has 6 standard degrees of freedom:

  1. Surge (longitudinal translation)
  2. Sway (lateral translation)
  3. Heave (vertical translation)
  4. Roll (rotation about longitudinal axis)
  5. Pitch (rotation about transverse axis)
  6. Yaw (rotation about vertical axis)

The index 5555 denotes pitch reaction caused by pitch motion (uncoupled single-degree-of-freedom pitch dynamics).

Transfer Function & Spectral Divide

The open-loop transfer function from total input torque to spar-buoy pitch angle θ\theta is:

G(s)=Θ(s)M(s)=1Itots2+B55s+C55G(s) = \frac{\Theta(s)}{M(s)} = \frac{1}{I_{\text{tot}} s^2 + B_{55} s + C_{55}}

In the frequency domain:

  • The spar-buoy structural pitch resonance occurs at ωn=0.2155 rad/s\omega_n = 0.2155\ \text{rad/s} (fn=0.0343 Hzf_n = 0.0343\ \text{Hz}).
  • Dominant ocean wave excitation occurs around fp≈0.10 Hzf_p \approx 0.10\ \text{Hz} (ωp=0.628 rad/s\omega_p = 0.628\ \text{rad/s}, period Tp≈10 sT_p \approx 10\ \text{s}).

Because the wave frequency is nearly three times the structural pitch resonance (ωp/ωn≈2.91\omega_p / \omega_n \approx 2.91), there exists a distinct frequency divide:

  • Structural resonance is low-frequency (0.034 Hz0.034\ \text{Hz}).
  • Wave swell excitation is higher-frequency (0.10 Hz0.10\ \text{Hz}).

Why Naive PD Control Fails

A conventional PD controller employs derivative action (KdsK_d s). Derivative gain increases at +20 dB/decade+20\ \text{dB/decade} with frequency:

∣Kd⋅jω∣=Kd⋅ω|K_d \cdot j\omega| = K_d \cdot \omega

At the wave frequency (0.10 Hz0.10\ \text{Hz}) the derivative term answers every radian of wave-induced tilt with Kp2+(Kdω)2≈561 MN⋅m/rad\sqrt{K_p^2 + (K_d\omega)^2} \approx 561\ \text{MN}\cdot\text{m}/\text{rad} of counter-torque, although the waves are not what it needs to damp. With the waves of this lesson (20 MN⋅m20\ \text{MN}\cdot\text{m}) the peak torque is only about 1.25 MN⋅m1.25\ \text{MN}\cdot\text{m} — 2.5%2.5\% of the ±50 MN⋅m\pm 50\ \text{MN}\cdot\text{m} limit, which would be reached only for a wave-frequency tilt of about 5.1∘5.1^\circ. Nearly all of that effort is spent at the wave frequency, where it does nothing for the resonance; the PD does raise the resonance damping from ζ=0.030\zeta = 0.030 to about 0.0540.054.

Active Wave Damping with a Notch Filter

The control solution places a second-order notch filter N(s)N(s) centered at the peak wave frequency (f0=0.10 Hzf_0 = 0.10\ \text{Hz}) directly into the feedback path:

N(s)=s2+ω02s2+2ζnω0s+ω02with ω0=2πf0N(s) = \frac{s^2 + \omega_0^2}{s^2 + 2\zeta_n \omega_0 s + \omega_0^2} \quad \text{with } \omega_0 = 2\pi f_0

At the wave excitation frequency (s=jω0s = j\omega_0), the numerator becomes zero:

∣N(jω0)∣=0(−∞ dB rejection)|N(j\omega_0)| = 0 \quad (-\infty\ \text{dB rejection})

The filter makes the controller blind to the first-order wave motion, allowing the vessel to ride the wave crests compliantly: in this lesson the peak torque and the dissipated energy drop by roughly 70%. It is not free: at the structural resonance (0.034 Hz0.034\ \text{Hz}) the notch has ∣N∣≈0.93|N| \approx 0.93 but lags the feedback by 21∘21^\circ, so the damping the PD adds shrinks by about 70% (closed-loop ζ≈0.037\zeta \approx 0.037 with the notch, 0.0540.054 without, 0.0300.030 open loop).

Reduced-Order Scope & Visual Context

This lesson models the platform's primary uncoupled pitch dynamics (1-DOF reduced-order spar-buoy). The turbine rotor, mooring catenary lines, and wave heave visual elements serve as illustrative context to aid physical intuition without burdening the fundamental control concepts with multi-megabyte hydroelastic simulation overhead.

See Also

  • StateSpace: compact 2-state MIMO representation of spar-buoy hydrodynamics.
  • TransferFunction: continuous implementation of the 2nd-order notch rejection filter.
  • Integrator: Newton-Euler angular acceleration to pitch angle state chain.
  • Saturation: realistic actuator torque saturation limits.
  • BodePlot: frequency-domain spectral analysis and frequency divide visualization.
  • uPlotDisplay: real-time telemetry tracing for tilt, rate, and torque signals.

Simulate in BlockWerk

Interactive sandbox directly in your browser — zero installation, WASM-powered.

Simulate in BlockWerk →