Kvant Physics Problem 56

A flag is modeled as a thin flexible sheet of characteristic length $L$ and mass per unit area $\sigma$, clamped along one edge to a rigid pole.

kvantphysics

Verified: no
Verdicts: FAIL + FAIL
Solve time: 15m41s
Source on kvant.digital

Problem

F56. Why does a flag “flutter” in the wind? (IX–X)

Setup and Assumptions

A flag is modeled as a thin flexible sheet of characteristic length $L$ and mass per unit area $\sigma$, clamped along one edge to a rigid pole. The wind flows with speed $v$ and density $\rho$, and the primary unknown is the mechanism responsible for sustained oscillatory motion of the flag, known as flutter. To capture the essential dynamics, the flag is approximated as a one-dimensional elastic structure undergoing small transverse deflections $y(x,t)$ along its spanwise coordinate $x$ measured from the attachment. This reduction captures the lowest bending mode of the flag while neglecting higher modes. The flag is assumed to be linearly elastic with bending rigidity $B$, and its mass per unit length is $\mu = \sigma L$ when projected along $x$. Air is treated as an incompressible fluid, and viscosity is retained only insofar as it produces a small linear damping term. Gravity-induced sagging and compressibility effects are neglected. The flow is assumed uniform far from the flag.

This simplification allows a tractable description of the energy transfer from the wind to the flag and an order-of-magnitude estimate of the critical wind speed for flutter while retaining the correct dimensional dependence on inertia, stiffness, and aerodynamic forces.

Physical Principles

The dynamics of the flag in the transverse direction follow the Euler-Bernoulli beam equation for small deflections under distributed aerodynamic loading $p(x,t)$, augmented by linear damping:

$\mu \frac{\partial^2 y}{\partial t^2} + \gamma \frac{\partial y}{\partial t} + B \frac{\partial^4 y}{\partial x^4} = p(x,t).$

The bending rigidity $B$ represents the resistance of the sheet to curvature, $\gamma$ is a linear damping coefficient per unit length, and $p(x,t)$ is the aerodynamic pressure. The pressure depends on the instantaneous flag shape and its velocity relative to the wind, producing a self-excited flow-structure interaction. Vortex shedding from the trailing edge introduces a characteristic timescale

$f_s = \mathrm{St} \frac{v}{L},$

with Strouhal number $\mathrm{St}$, which governs the phase lag between the flag motion and the aerodynamic forcing. Flutter occurs when aerodynamic forces inject energy into the flag oscillations faster than mechanical dissipation can remove it, leading to a self-sustained oscillation.

Derivation

For small-amplitude oscillations, the aerodynamic force can be linearized around the equilibrium configuration. Denoting the modal amplitude of the first bending mode as $\Theta(t)$, the distributed forcing produces an effective modal aerodynamic torque of the form

$M_{\rm aero}(\Theta, \dot{\Theta}) = k_a \Theta + \gamma_a \dot{\Theta},$

where $k_a$ and $\gamma_a$ are effective aerodynamic stiffness and damping coefficients, derived from integrating the linearized pressure distribution along the flag and projecting onto the modal shape. The signs of these coefficients are such that $\gamma_a > 0$ corresponds to energy input into the motion, reducing the effective damping. The equation of motion for the first bending mode then becomes

$I \ddot{\Theta} + (\gamma - \gamma_a) \dot{\Theta} + (k - k_a) \Theta = 0,$

where $I$ is the modal moment of inertia and $k$ is the effective modal stiffness associated with the bending rigidity $B$ via $k \sim B/L$. The eigenvalues $\lambda$ of the system satisfy

$\lambda = \frac{-(\gamma - \gamma_a) \pm \sqrt{(\gamma - \gamma_a)^2 - 4 I (k - k_a)}}{2 I}.$

Instability occurs when the real part of $\lambda$ becomes positive, which requires $\gamma_a > \gamma$. Simultaneously, oscillatory growth is guaranteed if $(k - k_a)/I > 0$, ensuring that the eigenvalues have a nonzero imaginary part. The critical condition for flutter is therefore

$\gamma_a(v_c) = \gamma \quad\text{and}\quad k_a(v_c) \lesssim k.$

To estimate the magnitude of $k_a$ and $\gamma_a$, consider the scaling of aerodynamic forces. The pressure difference induced by the flow on a small transverse deflection $\Theta$ scales as $\Delta p \sim \rho v \frac{\partial y}{\partial t} \sim \rho v L \dot{\Theta}$. Integrating along the flag length, the aerodynamic torque scales as

$M_{\rm aero} \sim \rho v^2 L^3 \Theta + \rho v L^2 \dot{\Theta}.$

Projecting onto the modal amplitude yields the effective aerodynamic coefficients

$k_a \sim \rho v^2 L^3, \quad \gamma_a \sim \rho v L^2.$

The mechanical modal stiffness and inertia are $k \sim B/L$ and $I \sim \mu L^2 \sim \sigma L^3$. The condition $\gamma_a = \gamma$ then gives the correct scaling for the critical wind speed:

$\gamma \sim \rho v_c L^2 \quad\Rightarrow\quad v_c \sim \frac{\gamma}{\rho L^2}.$

Substituting $\gamma \sim \eta L$, where $\eta$ is a damping per unit length with dimensions $[{\rm M/T}]$, gives

$v_c \sim \frac{\eta L}{\rho L^2} \sim \frac{\eta}{\rho L},$

which has dimensions of speed. The effective stiffness $k_a$ contributes to the oscillation frequency and slightly modifies $v_c$ but does not change the leading-order scaling. The critical wind speed therefore depends on the ratio of mechanical damping to aerodynamic mass and the flag length, capturing the onset of self-excited flutter with the correct dimensional dependence.

Result

The onset of flutter occurs when the aerodynamic damping overcomes mechanical dissipation, producing a net negative damping and self-excited oscillations. Within the linearized first-mode approximation, the critical wind speed scales as

$v_c \sim \frac{\gamma}{\rho L^2},$

with $\gamma$ representing the modal mechanical damping, $\rho$ the air density, and $L$ the characteristic flag length. Flags with larger damping require stronger winds to flutter, while more flexible or longer flags lower the threshold. The eigenvalue analysis ensures that both structural inertia $I$ and bending stiffness $k$ enter the oscillation frequency, whereas the aerodynamic forcing produces a phase lag that enables sustained oscillation. This result reproduces the essential mechanism of flutter, with dimensional consistency and explicit scaling of all relevant quantities.

Sanity Checks

The formula is dimensionally consistent, as $\gamma/(\rho L^2)$ has units of $L/T$. As $\rho \to 0$, $v_c \to \infty$, indicating no flutter in vacuum. Increasing mechanical damping $\gamma$ raises the critical wind speed, while increasing flag length $L$ or reducing stiffness $k$ lowers it, consistent with physical intuition. The inclusion of inertia $I$ in the modal dynamics ensures that the natural frequency of the flag is correctly captured, and the scaling of aerodynamic coefficients reflects the energy transfer from the flow to the structure. The simplified model reproduces the qualitative behavior observed in real flags and provides an order-of-magnitude estimate of the wind speed required for self-sustained oscillations, with the understanding that full distributed-mode analysis and detailed fluid-structure interaction refine the numerical prefactor.

$\boxed{v_c \sim \frac{\gamma}{\rho L^2}}$