SJMaths Air Safety Officer
Day 12 Fluid Dynamics Potential Flow Theory Consolidated Master Guide

Potential Flow Past a Cylinder (Non-Lifting & Lifting)

Complete unified master guide covering classical potential flow past a circular cylinder: Uniform flow, doublet superposition, surface velocity and pressure distributions, D'Alembert's paradox, vortex circulation ($Gamma$), stagnation point migration, Magnus effect, Kutta-Joukowski lift theorem, and real-flow vortex shedding for UPSC Air Safety Officer aspirants.

How This Topic Is Useful for Air Safety Officers & Aircraft Operations
DGCA CAR & ICAO Mandate

Cylindrical geometry represents one of the most common aerodynamic forms encountered on aircraft, from pitot-static sensor probes to landing gear oleo struts and external antennas:

Pillar 1: Formal Technical Definition & Nomenclature

Potential Flow Past a Circular Cylinder: An exact analytical solution of Laplace's equation ($ abla^2phi = 0$ or $ abla^2psi = 0$) obtained by the linear superposition of elementary ideal, inviscid, incompressible, and irrotational planar flowfields.

Fundamental Superposition Elements

  • Non-Lifting Cylinder: Uniform Free Stream (velocity $U_infty$) + Doublet (strength $kappa = 2pi U_infty R^2$) located at the coordinate origin $(0,0)$.
  • Lifting Cylinder: Uniform Free Stream + Doublet + Irrotational Point Vortex of circulation strength $Gamma$ located at the origin.

Key Symbols & Nomenclature

Symbol Physical Quantity SI Units Description / Role
$U_infty$ Free-stream velocity $ ext{m/s}$ Undisturbed approach velocity parallel to positive x-axis
$R$ Cylinder radius $ ext{m}$ Radius of the circular boundary streamline ($psi = 0$)
$kappa$ Doublet strength $ ext{m}^3/ ext{s}$ Magnitude of the sink-source dipole ($kappa = 2pi U_infty R^2$)
$Gamma$ Vortex Circulation $ ext{m}^2/ ext{s}$ Line integral of velocity along a closed curve enclosing the cylinder
$C_p$ Pressure Coefficient Dimensionless $C_p = rac{p - p_infty}{ rac{1}{2} ho U_infty^2} = 1 - 4sin^2 heta$ (non-lifting)
$L'$ Sectional Lift Force $ ext{N/m}$ Aerodynamic force perpendicular to free stream per unit span

Pillar 2: Microscopic Physics & Continuum Mechanics

1. Streamline Superposition & Boundary Formation

In potential flow, the stream function $psi$ satisfies the linear Laplace equation $ abla^2psi = 0$. By linearity, individual solutions can be algebraically added:

$$psi_{ ext{total}} = psi_{ ext{uniform}} + psi_{ ext{doublet}} + psi_{ ext{vortex}}$$ $$psi(r, heta) = U_infty r sin heta - rac{kappa}{2pi r}sin heta + rac{Gamma}{2pi}lnleft( rac{r}{R} ight)$$

Setting doublet strength $kappa = 2pi U_infty R^2$ gives:

$$psi(r, heta) = U_infty left(r - rac{R^2}{r} ight)sin heta + rac{Gamma}{2pi}lnleft( rac{r}{R} ight)$$

At $r = R$, the first term identically vanishes for all $ heta$, and when $Gamma=0$, $psi(R, heta) = 0$. Thus, the circular cylinder surface $r = R$ is an exact solid streamline boundary across which no fluid can penetrate.

2. D'Alembert's Paradox vs. Viscous Reality

Potential Flow (Inviscid Ideal Fluid): The surface velocity $v_ heta = -2U_infty sin heta$ is completely symmetric between front ($ heta = 180^circ$) and back ($ heta = 0^circ$). The pressure coefficient $C_p = 1 - 4sin^2 heta$ recovers to $+1.0$ at the rear stagnation point. Integration yields: $$D' = oint -p cos heta , R , d heta = 0, quad L' = oint -p sin heta , R , d heta = 0$$ This is D'Alembert's Paradox: zero net drag in steady potential flow.
Real Fluid (Viscous Navier-Stokes Reality): In real fluids with finite viscosity ($mu > 0$), fluid particles in the boundary layer cannot overcome the strong adverse pressure gradient from $ heta = 90^circ$ ($C_p = -3$) to $ heta = 0^circ$ ($C_p = +1$). The boundary layer separates near $ heta approx 82^circ$ (laminar) or $ heta approx 120^circ$ (turbulent), forming a wide low-pressure wake that generates massive pressure/form drag ($C_D approx 1.2$).

Pillar 3: Mathematical Derivations & Governing Laws

1. Surface Velocity & Pressure Distribution (Non-Lifting Cylinder)

The radial and tangential velocity components in polar coordinates are:

$$v_r = rac{1}{r} rac{partialpsi}{partial heta} = U_inftyleft(1 - rac{R^2}{r^2} ight)cos heta$$ $$v_ heta = - rac{partialpsi}{partial r} = -U_inftyleft(1 + rac{R^2}{r^2} ight)sin heta$$

On the cylinder surface ($r = R$):

$$v_r(R, heta) = 0 quad ( ext{Impermeable wall condition})$$ $$v_ heta(R, heta) = -2U_infty sin heta$$

Applying Bernoulli's equation between free stream and the surface:

$$p_infty + rac{1}{2} ho U_infty^2 = p + rac{1}{2} ho v_ heta^2 = p + rac{1}{2} ho (4U_infty^2 sin^2 heta)$$ $$C_p = rac{p - p_infty}{ rac{1}{2} ho U_infty^2} = 1 - rac{v_ heta^2}{U_infty^2} = 1 - 4sin^2 heta$$

Key Values of $C_p$ and Velocity

  • At Stagnation Points ($ heta = 0^circ, 180^circ$): $sin heta = 0 implies v_ heta = 0 implies C_p = +1.0$ (Maximum Pressure).
  • At Top & Bottom Crests ($ heta = 90^circ, 270^circ$): $sin heta = pm 1 implies |v_ heta| = 2U_infty implies C_p = 1 - 4(1)^2 = -3.0$ (Peak Suction).
  • Zero Pressure Points ($C_p = 0$): $1 - 4sin^2 heta = 0 implies sin heta = pm rac{1}{2} implies heta = 30^circ, 150^circ, 210^circ, 330^circ$.

2. Addition of Circulation $Gamma$ (Lifting Cylinder & Magnus Effect)

Adding a clockwise vortex of circulation $Gamma$ gives tangential velocity on the surface:

$$v_ heta(R, heta) = -2U_infty sin heta - rac{Gamma}{2pi R}$$

Setting $v_ heta = 0$ yields the stagnation point coordinates:

$$sin heta_s = - rac{Gamma}{4pi U_infty R}$$

Three Regimes of Stagnation Point Behavior

Case 1: $Gamma < 4pi U_infty R$ $sin heta_s > -1$. Two distinct stagnation points exist on the lower surface of the cylinder, symmetric about the vertical axis at $ heta_1 = sin^{-1}(- rac{Gamma}{4pi U_infty R})$ and $ heta_2 = pi - heta_1$.
Case 2: $Gamma = 4pi U_infty R$ $sin heta_s = -1$. The two stagnation points coalesce into a single point at the very bottom of the cylinder ($ heta = -pi/2$ or $270^circ$).
Case 3: $Gamma > 4pi U_infty R$ $|sin heta_s| > 1$. No real physical angle on the surface can satisfy the condition. The single stagnation point detaches from the surface and hangs in the fluid below the cylinder at $r_s = rac{Gamma + sqrt{Gamma^2 - 16pi^2 U_infty^2 R^2}}{4pi U_infty}$.

3. Proof of Kutta-Joukowski Lift Theorem

Integrating surface pressure around the cylinder circumference:

$$L' = -int_0^{2pi} p sin heta , R , d heta$$

Using $p = p_0 - rac{1}{2} ho v_ heta^2 = p_0 - rac{1}{2} ho left(2U_inftysin heta + rac{Gamma}{2pi R} ight)^2$:

$$v_ heta^2 = 4U_infty^2sin^2 heta + rac{2U_inftyGamma}{pi R}sin heta + rac{Gamma^2}{4pi^2 R^2}$$

When multiplied by $sin heta$, only the cross-term $sin^2 heta$ survives integration from $0$ to $2pi$ ($int_0^{2pi}sin^2 heta,d heta = pi$):

$$L' = rac{1}{2} ho rac{2U_infty Gamma}{pi R} cdot pi R = ho_infty U_infty Gamma$$ $$D' = -int_0^{2pi} p cos heta , R , d heta = 0$$

Universal Law: Any cylindrical or airfoil body with circulation $Gamma$ generates lift $L' = ho_infty U_infty Gamma$ and zero wave/induced drag in 2D potential flow.

Pillar 4: Dimensional Analysis, SI/CGS Units & Conversion Factors

Parameter Dimensional Formula SI Unit Physical Meaning
Circulation ($Gamma$) $[M^0 L^2 T^{-1}]$ $ ext{m}^2/ ext{s}$ Vorticity flux enclosed by curve
Doublet Strength ($kappa$) $[M^0 L^3 T^{-1}]$ $ ext{m}^3/ ext{s}$ Dipole source-sink moment
Sectional Lift ($L'$) $[M^1 L^0 T^{-2}]$ $ ext{N/m}$ ($ ext{kg/s}^2$) Force per unit span length
Strouhal Number ($St$) $[M^0 L^0 T^0]$ Dimensionless $St = rac{f D}{U_infty} approx 0.21$

Pillar 5: Aviation Engineering Reference Data

Flow Regimes for Crossflow Past a Circular Cylinder

Reynolds Number ($Re_D$) Flow Pattern / Phenomenon Separation Point ($ heta_{sep}$) Drag Coeff ($C_D$) Aviation Application
$Re_D < 5$ Creeping Stokes flow (unseparated) No separation $sim 24/Re_D$ Aerosol droplet sampling
$40 < Re_D < 300$ Laminar Kármán vortex street $82^circ$ $sim 1.0$ Sensor wire / pitot probe singing
$300 < Re_D < 3 imes 10^5$ Subcritical regime; laminar boundary layer, turbulent wake $82^circ$ $1.2$ Landing gear struts, antennas
$3 imes 10^5 < Re_D < 3.5 imes 10^6$ Supercritical regime; turbulent boundary layer ("Drag Crisis") $120^circ$ $0.30 ext{ to } 0.35$ High-speed transport landing gear crossflow

Pillar 6: Aircraft Airworthiness, Systems & DGCA CAR / ICAO Compliance

  • Pitot-Static Probe Position Error (DGCA CAR Section 2, Series I): Pitot-static masts are exposed to crossflow. The potential flow solution explains why static pressure ports must be located at exact angular positions ($ heta approx 30^circ$ or $39.2^circ$) where $C_p = 0$, or calibrated meticulously across angles of attack.
  • Vortex-Induced Vibration (VIV) & Antenna Failure: In accordance with ICAO Annex 8 airworthiness design requirements, all external cylindrical probes, masts, and blade antennas must be dynamically analyzed for resonance. Shedding frequency $f = rac{St cdot U_infty}{D}$ must never coincide with the structural natural frequency of the mast.
  • Landing Gear Aerodynamic Fairings: Because an unfaired circular cylinder has $C_D approx 1.2$ compared to an airfoil's $C_D approx 0.06$, DGCA airworthiness directives enforce streamlined fairings on fixed gear to reduce drag by over 90%.

    • Pillar 7: High-Frequency UPSC Examiner Traps & Misconceptions

    Trap 1: Peak Surface Velocity Location

    ❌ Examiner Trap: Thinking maximum velocity occurs at the front stagnation point ($ heta = 180^circ$).

    ✅ UPSC Truth: At stagnation points, velocity is zero. Maximum velocity occurs at top and bottom crests ($ heta = 90^circ, 270^circ$) where $v = 2 U_infty$.

    Trap 2: Minimum Pressure Coefficient Value

    ❌ Examiner Trap: Assuming the minimum pressure coefficient on a cylinder is $-1.0$ or $-4.0$.

    ✅ UPSC Truth: $C_p = 1 - 4sin^2 heta$. When $sin heta = pm 1$, $C_{p,min} = 1 - 4 = -3.0$.

    Trap 3: Points Where Pressure Equals Free-Stream Pressure ($C_p = 0$)

    ❌ Examiner Trap: Assuming $p = p_infty$ at $ heta = 45^circ$.

    ✅ UPSC Truth: $1 - 4sin^2 heta = 0 implies sin heta = 0.5 implies heta = 30^circ, 150^circ, 210^circ, 330^circ$.

    Trap 4: Circulation Required for Single Stagnation Point

    ❌ Examiner Trap: Believing single stagnation point occurs at $Gamma = 2pi U_infty R$.

    ✅ UPSC Truth: Single stagnation point at the bottom occurs when $sin heta = -1 implies Gamma = 4pi U_infty R$.

    Pillar 8: Step-by-Step Worked Aircraft Numerical Problem

    Realistic Flight Test Scenario

    A cylindrical flight test sensor mast of diameter $D = 0.4 ext{ m}$ (radius $R = 0.2 ext{ m}$) is tested in a wind tunnel at an airspeed $U_infty = 50 ext{ m/s}$ in sea-level standard air ($ ho_infty = 1.225 ext{ kg/m}^3$). A spinning rotor mechanism imparts a circulation of $Gamma = 25 ext{ m}^2/ ext{s}$ to the cylinder. Calculate: (a) the angular locations of the two stagnation points on the cylinder surface, and (b) the lift force per unit span generated by the cylinder.

    Given Parameters:
    Radius $R = 0.2 ext{ m}$
    Velocity $U_infty = 50 ext{ m/s}$
    Density $ ho_infty = 1.225 ext{ kg/m}^3$
    Circulation $Gamma = 25.0 ext{ m}^2/ ext{s}$
    Governing Equations:
    $$sin heta_s = - rac{Gamma}{4pi U_infty R}$$ $$L' = ho_infty U_infty Gamma$$
    Step 1: Calculate Stagnation Point Sine Value $$sin heta_s = - rac{25}{4 imes pi imes 50 imes 0.2} = - rac{25}{40pi} = - rac{25}{125.6637} approx -0.1989$$
    Step 2: Determine Stagnation Angles $ heta_1$ and $ heta_2$ $$ heta_s = sin^{-1}(-0.1989) approx -11.47^circ ext{ or } 348.53^circ$$ $$ heta_2 = 180^circ - (-11.47^circ) = 191.47^circ$$

    Both stagnation points have shifted downward below the horizontal centerline ($191.47^circ$ and $348.53^circ$).

    Step 3: Compute Sectional Lift Force (Kutta-Joukowski) $$L' = ho_infty U_infty Gamma = 1.225 imes 50 imes 25 = 1531.25 ext{ N/m}$$
    Final Conclusion & Answer:

    (a) The stagnation points are located at $191.5^circ$ and $348.5^circ$ on the lower half of the cylinder.
    (b) The cylinder produces a massive upward aerodynamic lift of $1,531.3 ext{ N per meter of span}$ (Magnus effect).

    High-Yield Practice MCQs: Potential Flow Past a Cylinder

    Curated directly from UPSC Engineering Services, GATE Aerospace, and DGCA Air Safety Officer exams.

    Q1. Flow past a non-lifting circular cylinder of radius $R$ in a uniform stream $U_infty$ is obtained by superposing:
    A A uniform flow and a source
    B A uniform flow and a sink
    C A uniform flow and a doublet of strength $kappa = 2pi U_infty R^2$
    D A source and a sink of equal strength
    Correct Answer: C
    Superposition of a uniform flow $U_infty$ and a doublet of strength $kappa = 2pi U_infty R^2$ forms a circular streamline of radius $R$.
    Q2. What is the maximum local velocity on the surface of a circular cylinder in non-lifting potential flow of free-stream velocity $U_infty$?
    A $U_infty$
    B $1.5 U_infty$
    C $2 U_infty$
    D $4 U_infty$
    Correct Answer: C
    Surface velocity is $v_ heta = -2U_inftysin heta$. At $ heta = 90^circ$ and $270^circ$, $|v_ heta| = 2U_infty$.
    Q3. The minimum pressure coefficient $C_{p,min}$ on the surface of a circular cylinder in ideal flow is:
    A 0
    B -1.0
    C -3.0
    D -4.0
    Correct Answer: C
    $C_p = 1 - 4sin^2 heta$. When $sin heta = pm 1$, $C_p = 1 - 4(1) = -3.0$.
    Q4. The condition for the two stagnation points to merge into a single point on the bottom of a lifting circular cylinder of radius $R$ in uniform flow $U_infty$ is:
    A $Gamma = 2pi U_infty R$
    B $Gamma = 4pi U_infty R$
    C $Gamma = pi U_infty R$
    D $Gamma = 8pi U_infty R$
    Correct Answer: B
    Stagnation angle is given by $sin heta_s = - rac{Gamma}{4pi U_infty R}$. Setting $sin heta_s = -1$ gives $Gamma = 4pi U_infty R$.
    Q5. D'Alembert's paradox states that:
    A Real viscous fluids cannot generate lift
    B A 2D body moving steadily through an inviscid fluid experiences zero net drag
    C Induced drag is zero for finite wings
    D Boundary layers never separate on circular cylinders
    Correct Answer: B
    D'Alembert's paradox proved mathematically that steady potential flow past any closed 2D body produces zero net aerodynamic drag force ($D' = 0$).
    ASO Speed Evaluation 10:00
    Question 1 of 3 2 Marks
    At which angular positions $ heta$ on the surface of a non-lifting circular cylinder does the local static pressure equal the free-stream static pressure ($C_p = 0$)?
    Question 2 of 3 2 Marks
    By the Kutta-Joukowski theorem, what is the lift per unit span generated by a rotating cylinder in air ($ ho = 1.2 ext{ kg/m}^3$) moving at $U_infty = 40 ext{ m/s}$ with circulation $Gamma = 10 ext{ m}^2/ ext{s}$?
    Question 3 of 3 2 Marks
    The phenomenon of alternate vortex shedding behind a cylinder that causes vibration of aircraft antennas is governed by which dimensionless number?

    Active Recall Flashcards & Governing Formulas

    Formula 1: Non-Lifting Surface Velocity & Pressure Distribution $$v_ heta = -2U_infty sin heta, quad C_p = 1 - 4sin^2 heta$$

    Max velocity $2U_infty$ and min $C_p = -3.0$ at $ heta = 90^circ, 270^circ$. Stagnation $C_p = 1.0$ at $0^circ, 180^circ$.

    Formula 2: Kutta-Joukowski Lift Theorem $$L' = ho_infty U_infty Gamma, quad D' = 0$$

    Lift is directly proportional to circulation $Gamma$ and perpendicular to relative wind.

    Formula 3: Stagnation Point Migration Equation $$sin heta_s = - rac{Gamma}{4pi U_infty R}$$

    $Gamma = 4pi U_infty R$ gives single stagnation point at $ heta = -90^circ$.

    Formula 4: Strouhal Number for Vortex Shedding $$St = rac{f D}{U_infty} approx 0.21$$

    Calculates vortex shedding frequency $f$ to prevent antenna resonance.