Extended five-term nonlinear drag model for a wide range of cylinder wakes

This paper proposes a new five-term nonlinear reduced-order drag model that incorporates both linear and quadratic coupling with the lift coefficient to accurately reproduce drag signals in cylinder wakes where traditional two-to-one frequency assumptions fail under harmonic excitation.

The Gist

The Big Picture: The "Dancing" Cylinder

Imagine a long pole (like a telephone pole or a flagpole) standing in a river. As the water flows past it, the water doesn't just slide smoothly; it swirls and peels off the back of the pole in little whirlpools called vortices.

These swirling vortices push and pull on the pole.

• Lift: The force pushing the pole side-to-side (like a leaf blowing in the wind).
• Drag: The force pushing the pole downstream (like the wind trying to knock the pole over).

For a long time, scientists had a simple rule of thumb for how these two forces relate: "The Drag force wiggles twice as fast as the Lift force."

Think of it like a drummer and a bassist. The bassist (Lift) plays a slow, steady beat. The drummer (Drag) plays a fast, double-time beat. This "2-to-1" rhythm worked perfectly for a pole that was either standing still or shaking strictly side-to-side.

The Problem: The "Tilted" Dance

The author of this paper, Osama Marzouk, noticed a problem. What happens if the pole doesn't just shake side-to-side, but wiggles in a tilted, diagonal direction?

Imagine the pole is dancing a tango. It's moving forward, backward, left, and right all at once. When this happens, the old "drummer/bassist" rule breaks. The Drag force starts wiggling at the same speed as the Lift force, not twice as fast. The rhythm gets messy, and the old math models fail to predict what's happening. It's like trying to predict a jazz solo using a strict marching band sheet music; it just doesn't fit.

The Solution: A New "Five-Term" Recipe

Marzouk decided to fix this by creating a new, more flexible mathematical model. He called it an "Extended Five-Term Nonlinear Drag Model."

Here is the recipe for his new model, broken down simply:

1. The Mean (The Average): The baseline weight of the drag.
2. The Quadratic Terms (The Old Rhythm): These are the old "2-to-1" rules that work when the pole shakes straight up and down.
3. The Linear Terms (The New Twist): These are the new ingredients added to the recipe. They account for the fact that when the pole tilts, the Drag force starts listening to the Lift force at the same speed.

The Analogy:
Think of the old model as a two-ingredient smoothie (Banana + Milk). It tastes great if you only have bananas.
The new model is a five-ingredient smoothie (Banana + Milk + Strawberries + Honey + Ice).

• If you only have bananas, the strawberries and honey disappear automatically, and you get the old smoothie.
• But if you have a "tilted" situation (strawberries), the new model adds them in, making the smoothie taste exactly right for that specific mix.

How They Tested It

To prove this new recipe works, the author didn't just guess; he built a virtual wind tunnel on his computer.

1. The Simulation: He used super-computers to simulate water flowing around a cylinder at a specific speed (Reynolds number 300).
2. The Experiment: He made the cylinder vibrate in different directions (straight up, straight down, and at weird angles).
3. The Comparison: He compared the computer's "real" data against his new "five-term" math model.

The Result:
The new model was a perfect match. It could predict the forces on the cylinder whether it was shaking straight up, shaking diagonally, or standing still. The old models failed miserably at the diagonal angles, but the new model handled them with ease.

Why Does This Matter?

This isn't just about math homework; it has real-world applications for safety and engineering.

• Offshore Oil Rigs: These are giant pipes in the ocean. If the currents push them diagonally, the old models might underestimate the stress, leading to a collapse.
• Wind Turbines: The towers are constantly vibrating. If the wind hits them from a weird angle, engineers need accurate models to ensure they don't break.
• Bridges and Skyscrapers: Understanding exactly how wind and water forces interact helps architects build structures that can survive storms without swaying too much.

The Bottom Line

The author discovered that the "old rules" of fluid dynamics are too rigid for the real world, where things rarely move in perfect straight lines. By adding a few extra "ingredients" (linear terms) to the mathematical model, he created a universal tool that works for almost any situation, making our predictions for bridges, ships, and turbines much safer and more accurate.

In short: He upgraded the map from a simple grid to a GPS that works even when you're driving off-road.

Technical Summary

1. Problem Statement

The paper addresses a limitation in existing reduced-order models (ROMs) for fluid-structure interaction (FSI) and vortex-induced vibration (VIV).

• The Assumption: Traditional wake oscillator models for cylinders assume a strict two-to-one frequency relationship between the lift ($C_L$) and drag ($C_D$) coefficients. Specifically, the fundamental frequency of the drag oscillation is assumed to be exactly twice the fundamental frequency of the lift oscillation (quadratic coupling).
• The Failure: This assumption holds for fixed cylinders or cylinders moving strictly in the cross-stream direction. However, it fails when a cylinder undergoes inclined harmonic vibration (vibrating at an angle between the streamwise and cross-stream directions). In these "non-traditional wakes," the drag signal develops significant components at the fundamental lift frequency, causing the standard quadratic-only models to inaccurately reproduce the drag coefficient time history and phase relationships.

2. Methodology

The study employs a combination of Computational Fluid Dynamics (CFD) and nonlinear dynamical systems analysis to develop and validate a new model.

• Numerical Simulation (DNS):
  • Governing Equations: The author solved the two-dimensional incompressible Navier–Stokes equations using Direct Numerical Simulation (DNS) at a Reynolds number ($Re$) of 300.
  • Setup: A cylinder was subjected to sinusoidal harmonic vibrations at various inclination angles (from $40^\circ$ to $90^\circ$ relative to the free stream) and a fixed amplitude of 7.5% of the cylinder diameter ($D$).
  • Technique: The Arbitrary Lagrangian-Eulerian (ALE) method was used to handle the moving mesh. The lift and drag coefficients were computed by integrating the pressure and friction coefficients over the cylinder surface.
• Data Analysis:
  • Spectral Analysis: Power spectra of the lift and drag signals were analyzed to identify frequency components ($a_{1L}, a_{2L}, a_{1D}, a_{2D}$, etc.).
  • Phase Analysis: The phase angles between specific frequency components of lift and drag were calculated.
  • Limit Cycle Projection: The trajectories of the lift-drag relationship were projected onto the $C_L$-$C_D$ plane to observe the shape of the limit cycles (e.g., transition from two-loop to single-loop structures).
• Model Development:
  • Based on the spectral and phase data, the author derived a new algebraic reduced-order model. The parameters of this model were determined using closed-form analytical expressions derived from the training data, rather than trial-and-error tuning.

3. Key Contributions

The primary contribution is the proposal of a new five-term nonlinear drag model that generalizes the relationship between lift and drag.

• The Five-Term Model: The proposed drag coefficient equation ($C_D$) includes:

  1. Mean Component: $C_{D,mean}$ (Time-averaged drag).
  2. Linear Coupling Term 1: $\lambda_1 \frac{a_{2D}}{a_{1L}} C_L$ (Linear dependence on lift).
  3. Linear Coupling Term 2: $\lambda_2 \frac{a_{2D}}{a_{1L}} \dot{C}_L$ (Linear dependence on the derivative of lift).
  4. Quadratic Coupling Term 1: $2q_1 \frac{a_{2D}}{a_{1L}^2} (C_L^2 - \overline{C_L^2})$ (Standard quadratic term).
  5. Quadratic Coupling Term 2: $2q_2 \frac{a_{2D}}{\omega_S a_{1L}^2} C_L \dot{C}_L$ (Standard quadratic derivative term).

• Key Innovations:

  • Inclusion of Linear Terms: The model explicitly adds linear coupling terms to account for the drag component at the fundamental lift frequency ($a_{1D}$), which becomes significant in inclined vibrations.
  • Frequency Adaptation: The model replaces the Strouhal frequency ($\omega_S$) of a fixed cylinder with the vibration frequency ($\Omega$) of the moving cylinder, making it applicable to forced oscillations.
  • Analytical Parameterization: The coefficients ($\lambda_1, \lambda_2, q_1, q_2$) are defined analytically based on the phase angles ($\psi$) between the spectral components of lift and drag, ensuring the model is fully analytical and not purely empirical.

4. Results

The proposed model was validated against the DNS training data across several motion cases:

• Fixed Cylinder & Cross-Stream Motion ($90^\circ$): The model successfully reduced to the traditional form. The linear terms naturally vanished (or became negligible), and the model accurately reproduced the classic two-loop limit cycle and the two-to-one frequency relationship.
• Inclined Motion ($40^\circ - 80^\circ$):
  • Spectral Shift: As the inclination angle decreased, the amplitude of the drag component at the lift frequency ($a_{1D}$) increased significantly, sometimes exceeding the component at twice the frequency ($a_{2D}$).
  • Limit Cycle Distortion: The $C_L$-$C_D$ projection transformed from a symmetric two-loop shape to a distorted single-loop shape.
  • Model Accuracy: The five-term model showed excellent agreement with the DNS simulation data for these complex cases.
  • Comparison: When the linear terms were removed (reverting to the old quadratic-only model), the reconstructed drag signal showed significant phase and amplitude errors, particularly in cases where $a_{1D} > a_{2D}$.
• Efficiency: The algebraic model can generate drag signals approximately three orders of magnitude faster than running a full CFD simulation, making it suitable for real-time FSI applications.

5. Significance

• Generalization of Wake Oscillators: This work extends the applicability of wake oscillator models beyond the restrictive "cross-stream only" or "fixed" cylinder scenarios. It provides a robust framework for modeling VIV and FSI in complex 3D environments where structures may vibrate in multiple degrees of freedom or at oblique angles.
• Physical Insight: The study clarifies the mechanism of lift-drag coupling in non-traditional wakes, demonstrating that the "two-to-one" rule is a special case rather than a universal law.
• Engineering Application: The proposed reduced-order model allows for efficient simulation of fluid-structure interactions in offshore risers, wind turbine towers, and power lines without the computational cost of full Navier-Stokes solvers, while maintaining high accuracy in predicting drag forces under complex motion conditions.
• Future Potential: The author notes that while tested at $Re=300$, the coherent vortex structures in 3D flows suggest the model's utility could extend to higher Reynolds numbers and finite-length cylinders, provided the spectral characteristics remain distinct.