Morison equation

The Morison equation is the sum of two force components: an inertia force in phase with the local flow acceleration and a drag force proportional to the (signed) square of the instantaneous flow velocity. The inertia force is of the functional form as found in potential flow theory, while the drag force has the form as found for a body placed in a steady flow. In the heuristic approach of Morison, O'Brien, Johnson and Schaaf these two force components, inertia and drag, are simply added to describe the inline force in an oscillatory flow. The transverse force—perpendicular to the flow direction, due to vortex shedding—has to be addressed separately.

The Morison equation contains two empirical hydrodynamic coefficients—an inertia coefficient and a drag coefficient—which are determined from experimental data. As shown by dimensional analysis and in experiments by Sarpkaya, these coefficients depend in general on the Keulegan–Carpenter number, Reynolds number and surface roughness.^[4][5]

The descriptions given below of the Morison equation are for uni-directional onflow conditions as well as body motion.

In an oscillatory flow with flow velocity ${\displaystyle u(t)}$, the Morison equation gives the inline force parallel to the flow direction:^[6]

${\displaystyle F\,=\,\underbrace {\rho \,C_{m}\,V\,{\dot {u}}} _{F_{I}}+\underbrace {{\frac {1}{2}}\,\rho \,C_{d}\,A\,u\,|u|} _{F_{D}},}$

where

• ${\displaystyle F(t)}$ is the total inline force on the object,
• ${\displaystyle {\dot {u}}\equiv {\text{d}}u/{\text{d}}t}$ is the flow acceleration, i.e. the time derivative of the flow velocity ${\displaystyle u(t),}$
• the inertia force ${\displaystyle F_{I}\,=\,\rho \,C_{m}\,V\,{\dot {u}}}$, is the sum of the Froude–Krylov force ${\displaystyle \rho \,V\,{\dot {u}}}$ and the hydrodynamic mass force ${\displaystyle \rho \,C_{a}\,V\,{\dot {u}},}$
• the drag force ${\displaystyle F_{D}\,=\,{\scriptstyle {\frac {1}{2}}}\,\rho \,C_{d}\,A\,u\,|u|}$ according to the drag equation,
• ${\displaystyle C_{m}=1+C_{a}}$ is the inertia coefficient, and ${\displaystyle C_{a}}$ the added mass coefficient,
• A is a reference area, e.g. the cross-sectional area of the body perpendicular to the flow direction,
• V is volume of the body.

For instance for a circular cylinder of diameter D in oscillatory flow, the reference area per unit cylinder length is ${\displaystyle A=D}$ and the cylinder volume per unit cylinder length is ${\displaystyle V={\scriptstyle {\frac {1}{4}}}\pi {D^{2}}}$. As a result, ${\displaystyle F(t)}$ is the total force per unit cylinder length:

${\displaystyle F\,=\,C_{m}\,\rho \,{\frac {\pi }{4}}D^{2}\,{\dot {u}}\,+\,C_{d}\,{\frac {1}{2}}\,\rho \,D\,u\,|u|.}$

Besides the inline force, there are also oscillatory lift forces perpendicular to the flow direction, due to vortex shedding. These are not covered by the Morison equation, which is only for the inline forces.

In case the body moves as well, with velocity ${\displaystyle v(t)}$, the Morison equation becomes:^[6]

${\displaystyle F=\underbrace {\rho \,V{\dot {u}}} _{a}+\underbrace {\rho \,C_{a}V\left({\dot {u}}-{\dot {v}}\right)} _{b}+\underbrace {{\frac {1}{2}}\rho \,C_{d}A\left(u-v\right)\left|u-v\right|} _{c}.}$

where the total force contributions are:

• a: Froude–Krylov force, due to the pressure gradient at the body's location induced by the fluid acceleration ${\displaystyle {\dot {u}}}$,
• b: hydrodynamic mass force,
• c: drag force.

Note that the added mass coefficient ${\displaystyle C_{a}}$ is related to the inertia coefficient ${\displaystyle C_{m}}$ as ${\displaystyle C_{m}=1+C_{a}}$.

• The Morison equation is a heuristic formulation of the force fluctuations in an oscillatory flow. The first assumption is that the flow acceleration is more-or-less uniform at the location of the body. For instance, for a vertical cylinder in surface gravity waves this requires that the diameter of the cylinder is much smaller than the wavelength. If the diameter of the body is not small compared to the wavelength, diffraction effects have to be taken into account.^[7]
• Second, it is assumed that the asymptotic forms: the inertia and drag force contributions, valid for very small and very large Keulegan–Carpenter numbers respectively, can just be added to describe the force fluctuations at intermediate Keulegan–Carpenter numbers. However, from experiments it is found that in this intermediate regime—where both drag and inertia are giving significant contributions—the Morison equation is not capable of describing the force history very well. Although the inertia and drag coefficients can be tuned to give the correct extreme values of the force.^[8]
• Third, when extended to orbital flow which is a case of non uni-directional flow, for instance encountered by a horizontal cylinder under waves, the Morison equation does not give a good representation of the forces as a function of time.^[9]