Airy wave (horizontal velocity): $u(z,t) = \dfrac{\pi H}{T}\dfrac{\cosh k(z+d)}{\sinh kd}\cos(\omega t)$
KC number: $KC = U_{max}T / D$, $KC < 5$: inertia-dominated, $KC > 20$: drag-dominated
Dispersion relation: $\omega^2 = gk\tanh(kd)$
Real-time inertia, drag, and total force calculation using the Morison equation. Airy wave velocity field visualization, overturning moment, and KC number.
Airy wave (horizontal velocity): $u(z,t) = \dfrac{\pi H}{T}\dfrac{\cosh k(z+d)}{\sinh kd}\cos(\omega t)$
KC number: $KC = U_{max}T / D$, $KC < 5$: inertia-dominated, $KC > 20$: drag-dominated
Dispersion relation: $\omega^2 = gk\tanh(kd)$
The total horizontal force per unit length on a cylindrical member is given by the Morison equation, which sums the inertia and drag components.
$$F(z,t) = \underbrace{\rho C_m \frac{\pi D^2}{4}\frac{\partial u}{\partial t}}_{\text{Inertia Force}}+ \underbrace{\frac{1}{2}\rho C_d D \, u|u|}_{\text{Drag Force}}$$Where:
$\rho$ = water density, $C_m$ = inertia coefficient, $C_d$ = drag coefficient, $D$ = cylinder diameter.
$u(z,t)$ = total horizontal fluid velocity (wave + current), $\frac{\partial u}{\partial t}$ = local fluid acceleration.
The $u|u|$ term ensures the drag force always opposes the direction of flow.
The fluid kinematics (velocity $u$ and acceleration) are modeled using linear Airy wave theory combined with a uniform current.
$$u(z,t) = \underbrace{\frac{\pi H}{T}\frac{\cosh(k(z+d))}{\sinh(kd)} \cos(\omega t)}_{\text{Wave Velocity}}+ \underbrace{U_c}_{\text{Current}}$$Where:
$H$ = wave height, $T$ = wave period, $d$ = water depth, $k$ = wave number ($2\pi/L$), $\omega$ = wave angular frequency ($2\pi/T$), $U_c$ = current velocity.
The $\cosh/\sinh$ term describes how the wave motion decays with depth $z$ below the surface. This is why the forces are highest at the water surface in the visualization.
Jacket Platform Design: Engineers use this exact analysis, per standards like DNV-OS-J101, to calculate the extreme wave and current loads on the tubular legs and braces of offshore jackets. The overturning moment calculated by the simulator is a direct input for the foundation design to prevent the platform from toppling.
CFD Model Validation: The coefficients $C_m$ and $C_d$ are not constants—they depend on surface roughness, Reynolds number, and KC number. High-fidelity OpenFOAM CFD simulations of flow around cylinders are run to calibrate these coefficients, which are then used in the simpler, faster Morison equation for full structural analysis.
FEM Load Input: The distributed force profile $F(z,t)$ calculated here is applied as an external load in Finite Element Analysis (FEA) software like ANSYS or ABAQUS. This allows engineers to compute the resulting stresses, deflections, and fatigue life in every weld and joint of the offshore structure.
Mooring Line & Riser Analysis: For floating platforms, the forces on vertical risers (pipes bringing oil to the surface) and mooring lines are assessed using the Morison equation. Understanding the relative magnitude of inertia vs. drag force is crucial for predicting vortex-induced vibrations (VIV) that can cause metal fatigue.
First and foremost, please do not think of the Morison equation as a universal tool that can calculate anything. It is targeted at slender, cylindrical members (slender bodies). For example, calculating the motion of a large floating structure (like a ship or a semi-submersible platform) in waves requires a different theory (potential flow theory). Another crucial point is that the "drag coefficient Cd" and the "inertia coefficient Cm" are not constants. While existing explanations mention the Keulegan–Carpenter (KC) number determining dominance, in reality, these coefficients vary significantly due to surface roughness (e.g., marine growth), flow turbulence, and the arrangement of columns (group effect). For instance, experimental values show that even for the same cylinder, Cd can vary from 0.6 to over 1.2. While tweaking coefficients in a tool is good for building intuition, for actual design, it is essential to refer to experimental values or code-recommended values appropriate for your specific environment.
Another point is to understand the limitations of Airy wave theory. This theory assumes "small-amplitude waves" where wave height is small compared to wavelength. This means it cannot strictly represent huge, steep waves with strong nonlinearity, like those during typhoons. For example, it might still be applicable for a wave height of 5m and wavelength of 50m (height/length = 0.1), but for a wave height of 10m and wavelength of 60m (ratio 0.17), consideration of higher-order theories like Stokes wave theory becomes necessary. Don't take the tool's output at face value; always get into the habit of asking yourself, "Does the input wave satisfy the small-amplitude wave assumption?"