Effect of flexibility on the pitch-heave flutter instability of a flexible foil elastically supported on its leading edge

This paper presents a new analytical tool that extends previous fluid-structure interaction formulations to characterize the parametric regions, frequencies, and growth rates of coupled pitch-heave flutter instabilities in flexible foils, demonstrating how reduced stiffness and spring constants significantly enlarge instability ranges to guide the design of future flexible oscillating foil turbines.

The Gist

Imagine you are holding a long, flexible ruler (like a diving board) out of a car window while driving. As the wind hits it, the ruler starts to wiggle. Sometimes, it just flutters harmlessly. Other times, it starts to shake violently, growing larger and larger until it might snap. In engineering, this dangerous, self-amplifying shaking is called flutter.

This paper is about understanding exactly when and why a flexible wing (or foil) starts to shake itself apart, and how we can use that knowledge to build better machines, like wind turbines that harvest energy from the wind.

Here is a breakdown of the paper's findings using simple analogies:

1. The Setup: The "Springy" Ruler

The author, R. Fernandez-Feria, created a mathematical tool to predict this shaking. He didn't just look at a rigid stick; he looked at a flexible foil (like a piece of thin metal or plastic) that is attached to a pivot point at its front edge (the leading edge).

Think of this pivot point as being held by two invisible springs:

• The Up-Down Spring: Allows the front to bounce up and down (Heave).
• The Twisting Spring: Allows the front to rotate up and down (Pitch).

The paper asks: If the foil is flexible, and these springs are loose or tight, how does the wind make it shake?

2. The Old Way vs. The New Way

The Old Way: Previous math models were like looking at a stiff ruler. They worked great for rigid wings or very stiff flexible ones. But if the foil was very floppy (like a wet noodle), the old math broke down. It couldn't predict the shaking accurately because it only looked at the "first bend" of the ruler.

The New Way: This paper introduces a smarter math model. It's like upgrading from looking at just the first bend of the ruler to watching two distinct waves travel along it.

• By including a second "wave" (a second bending mode), the new model can predict the shaking of much more flexible foils (down to stiffness levels 10 times lower than before).
• It's like the difference between predicting how a stiff diving board moves versus predicting how a wet noodle flaps in the wind. The new model catches the "noodle" behavior that the old one missed.

3. The Three Scenarios of Shaking

The author tested three different ways the foil could move, finding some surprising results:

A. The "Locked" Ruler (Clamped)

Imagine the front of the ruler is glued tight so it can't move up/down or twist. It can only bend.

• Finding: If the ruler is too floppy (low stiffness), it will start to shake violently if the wind is fast enough and the ruler is heavy enough.
• Analogy: Think of a flag in the wind. If the pole is too flimsy, the whole flag flaps wildly. The paper maps out exactly how flimsy the pole can be before it goes crazy.

B. The "One-Way" Ruler (Only Up/Down or Only Twist)

Imagine the front is on a track so it can only bounce up and down, but not twist. Or vice versa.

• The Old Belief: Engineers used to think a rigid wing needs both bouncing and twisting to shake itself apart. If you lock one, it's safe.
• The New Discovery: This is true for rigid wings. But for flexible wings, you can get a dangerous shake even if you only allow bouncing!
• The "Coupling" Effect: When the spring holding the front is loose, the "bouncing" motion of the front gets tangled up with the "bending" of the flexible body. They start dancing together, and the shaking gets much worse, much faster. It's like a child on a swing; if you push at the right time (coupling), the swing goes higher.

C. The "Free" Ruler (Both Up/Down and Twist)

This is the classic scenario where the front can bounce and twist freely.

• Finding: Flexibility makes things worse. A flexible wing starts shaking at lower wind speeds than a rigid one.
• The "Sweet Spot": There is a specific range of stiffness and weight where the shaking is most violent. The paper provides a map to find this "danger zone" so engineers can avoid it.

4. Why Does This Matter? (The "So What?")

You might ask, "Why do we care about shaking rulers?"

1. Safety: In airplanes and bridges, we want to avoid flutter because it causes structural failure (things breaking). This new tool helps engineers design wings that are safe even if they are slightly flexible.
2. Energy Harvesting: This is the exciting part! Some new wind turbines don't have spinning blades; they have flexible flaps that flap in the wind to generate electricity.
   • To make these turbines efficient, you want them to flutter, but you need to control it so they don't break.
   • This paper acts as a recipe book. It tells engineers: "If you use a foil of this thickness, attach it with springs of this strength, and put it in wind of this speed, it will flap perfectly to generate power without snapping."

5. The "Gravity" Factor

The paper also added gravity to the mix. Imagine the foil is heavy and submerged in water. Gravity pulls it down, changing its resting shape before the wind even hits it.

• The math shows that gravity changes the "starting position" of the foil, which slightly alters how it reacts to the wind. It's like setting a diving board lower in the pool; the way it bounces changes slightly.

Summary

This paper is a mathematical map for the chaotic world of flexible wings in the wind.

• Old Map: Good for stiff wings, failed for floppy ones.
• New Map: Accurate for floppy wings, predicts exactly when they will shake, how fast they will shake, and how to tune the springs to either stop the shaking (for safety) or encourage it (for energy).

It turns a complex physics problem into a set of clear rules, helping us build better, safer, and more efficient machines that dance with the wind.

Technical Summary

1. Problem Statement

The paper addresses the flutter instability of a two-dimensional flexible foil immersed in an inviscid fluid current. Unlike traditional studies focusing on rigid foils or foils with a fixed (clamped) leading edge, this work investigates a foil that is elastically mounted on its leading edge via linear (heave) and torsional (pitch) springs and dampers.

The primary objective is to characterize how foil flexibility interacts with the passive degrees of freedom (heave and pitch) to alter the onset, growth rates, and parametric regions of flutter. The study aims to extend existing analytical models to cover lower stiffness regimes where previous approximations failed, specifically targeting the coupling between structural spring modes and flexural deformation modes.

2. Methodology

The author employs an analytical approach based on linear potential flow theory and Euler-Bernoulli beam theory, avoiding the high computational cost of numerical CFD simulations while maintaining physical clarity.

• Governing Equations: The fluid-structure interaction (FSI) is modeled by coupling the linearized inviscid flow equations with the Euler-Bernoulli beam equation. Gravity is explicitly included, affecting the static equilibrium shape.
• Approximation of Deformation: To capture the first two natural bending modes of the flexible foil, the foil's displacement $z_s(x,t)$ is approximated using a fifth-order polynomial in space. This expansion includes:
  • Passive heave ($h$) and pitch ($\alpha$) at the leading edge.
  • Two flexural deformation modes ($d_1$ and $d_2$).
• Moment Method: The governing partial differential equations are converted into a system of four ordinary differential equations by taking the first four moments of the beam equation (multiplying by $1, (x+1), (x+1)^2, (x+1)^3$ and integrating).
• Fluid Forces: Unsteady aerodynamic forces (lift, moment, and flexural moments) are calculated using Theodorsen's function for harmonic motion. Closed-form analytical expressions for these forces are derived for the specific polynomial deformation shape.
• Stability Analysis: The system is linearized around a static equilibrium state. By assuming harmonic perturbations ($e^{i\gamma t}$), the problem is reduced to a complex eigenvalue problem ($\det[\mathbf{A}(\gamma)] = 0$).
  • $\gamma = k + i\sigma$: $k$ is the reduced frequency, and $\sigma$ is the growth rate.
  • Instability occurs when $\sigma < 0$ (growth rate is positive).

3. Key Contributions

1. Extended Analytical Validity: The paper extends a previous analytical formulation (which considered only the first flexural mode) by including a second flexural mode. This significantly extends the validity of the analytical tool from stiffness parameters $S \sim O(10)$ down to $S \sim O(10^{-1})$.
2. Inclusion of Gravity and Equilibrium: The model incorporates gravity to determine the static equilibrium location and shape of the foil in both stagnant and flowing fluids, providing necessary initial conditions for stability analysis.
3. Discovery of Coupled Instability Modes: The study identifies new instability mechanisms where the spring instability mode (heave or pitch) couples with the flexural instability mode as stiffness decreases. This coupling was not captured by single-mode analytical models.
4. Parametric Mapping: The tool provides a comprehensive mapping of flutter regions in terms of mass ratio ($R$), stiffness parameter ($S$), and spring constants ($k_h, k_a$), offering critical flutter velocities and frequencies.

4. Key Results

A. Validation

• The analytical results were validated against numerical data (Alben, Michelin, etc.) for clamped foils.
• The model accurately captures the first two natural flexural modes and their associated flutter instabilities for $S \gtrsim 0.3$.
• It correctly predicts the divergence of growth rates and frequencies as stiffness decreases, a region where single-mode models fail.
• Limitation: The model does not capture "flapping-flag" instabilities occurring at extremely low stiffness ($S \lesssim 10^{-1}$), which require higher-order modes.

B. Clamped Foil ($k_h, k_a \to \infty$)

• Flutter instabilities for a clamped flexible foil exist only above a critical mass ratio $R$ that depends on stiffness $S$.
• Growth rates are negligible for high stiffness but increase dramatically as $S$ drops below $\approx 10$.
• The second flexural mode (higher frequency) becomes unstable at lower stiffness values than the first mode.

C. Passive Heave Only ($k_a \to \infty$) or Pitch Only ($k_h \to \infty$)

• Rigid Foil Limit: A rigid foil with only one degree of freedom (heave or pitch) is stable. Flutter requires coupling between heave and pitch.
• Flexible Foil: Flexibility introduces a new instability mechanism. As the spring constant decreases, the flexural mode couples with the spring mode.
• This coupling creates a new unstable region below a threshold stiffness $S$ and above a threshold mass ratio $R$.
• The growth rate increases substantially in this coupled region compared to the uncoupled flexural mode.

D. Coupled Pitch-Heave ($k_h, k_a$ finite)

• For rigid foils, flutter occurs in a specific region of the $(k_h, k_a)$ plane.
• Effect of Flexibility: Introducing flexibility ($S$ finite) enlarges the instability region and lowers the critical mass ratio required for flutter.
• The coupling between the spring modes and the flexural modes (specifically the second mode $d_2$) significantly increases the growth rate, particularly when $S \lesssim 10$.
• The instability region shrinks as the mass ratio $R$ decreases.

E. Effect of Damping

• Increasing damping constants ($b_h, b_a$) reduces both the parametric region of instability and the growth rates.
• Removing damping ($b=0$) significantly expands the unstable region and increases growth rates, confirming the stabilizing role of structural damping.

5. Significance and Applications

• Design Guide for Energy Harvesting: The results are directly applicable to the design of flapping-foil turbines for flow-induced energy harvesting. The analytical tool allows engineers to quickly estimate critical flutter velocities and optimal stiffness/mass ratios to maximize energy extraction without running expensive numerical simulations.
• Aeroelastic Stability: The study clarifies the transition from rigid-body flutter to flexible-foil flutter, highlighting that flexibility can either stabilize or destabilize a system depending on the coupling with spring constants.
• Computational Efficiency: By providing a closed-form analytical solution valid for a wide range of stiffnesses, the paper offers a fast, reliable alternative to numerical FSI solvers for preliminary design and parametric studies of flexible structures in fluid flows.

In summary, this paper provides a robust analytical framework that bridges the gap between rigid-body aeroelasticity and highly flexible flag-like dynamics, specifically elucidating how the coupling of structural springs and foil flexibility creates new, highly unstable flutter regimes relevant to engineering applications.