A Similarity Solution of Rear Stagnation-point Flow over a Flat Plate in Two Dimensions

This paper investigates the development of vortex shedding in two-dimensional unsteady incompressible flow at the rear stagnation point of a flat plate.

The Gist

Imagine you are standing in a river, holding a flat board straight up against the current. The water hits the center of your board and splits, flowing smoothly around the edges. This is a "forward stagnation point," and it's a calm, predictable place.

Now, imagine the water flowing away from the back of that board. This is the Rear Stagnation Point. The paper you shared investigates what happens in this specific, tricky zone behind the object. It's a bit like the "wake" behind a boat, but focused on the exact spot where the water tries to rejoin after being pushed apart.

Here is the story of the paper, broken down into simple concepts and everyday analogies:

1. The Setup: The "Traffic Jam" at the Back

When fluid (like water or air) flows past a flat plate, it creates a specific pattern.

• The Ideal World: In a perfect, frictionless world, the water would just glide away smoothly.
• The Real World: Because water has "stickiness" (viscosity), it clings to the plate. As it tries to flow away from the back, it gets confused. It creates a vortex (a swirling whirlpool) and sometimes even flows backward against the main current.

The authors are trying to solve a math puzzle: Can we write a single, perfect formula that describes this messy, swirling behavior behind the plate?

2. The Magic Key: The "Strouhal Number" (The Rhythm of the Swirl)

To solve this, the authors use a special number called the Strouhal number (let's call it $\kappa$ or "kappa").

• Think of it like a metronome: If you tap your foot to a beat, the Strouhal number is the speed of that beat.
• In this flow, it represents how fast the swirling vortices are shedding (dropping off) from the wall.
• The paper asks: Does the math work for every possible beat speed, or are there some rhythms where the math breaks down?

3. The Big Discovery: The "Impossible" Rhythms

The authors found that the math behaves like a picky musician.

• The "Silent" Rhythm ($\kappa = 0$): They proved that if the flow is perfectly steady (no wobbling, no shedding), the math cannot satisfy the rules of the wall. It's like trying to play a song on a piano where the keys are broken; no matter how you press them, you can't get the right note. The flow must be unsteady (wobbly) to exist.
• The "Chaotic" Rhythms ($\kappa > -1.5$): When the rhythm gets too fast or too specific, the math predicts a "singularity." Imagine a traffic jam so bad that the cars suddenly stop existing or the road disappears. In physics terms, the velocity drops to zero instantly or becomes infinite. This tells us the model breaks down, likely because the flow has turned turbulent (chaotic) and our simple formulas can't handle it anymore.

4. The "Magic" Rhythm ($\kappa = -2$)

The paper found one very special, rare rhythm where the math works perfectly.

• When the Strouhal number is exactly -2, the equations simplify into a beautiful, clean solution.
• The Analogy: Imagine a dancer who moves in a perfect, repeating circle. At this specific rhythm, the swirling water behind the plate moves in a predictable, repeating pattern that matches the laws of physics perfectly.
• This solution shows that the water flows away from the wall, but also has a periodic "wiggle" (oscillation) that creates the vortices.

5. The Computer Simulation: Testing the Theory

Since the math is so hard to solve by hand, the authors used a computer (MATLAB) to simulate the flow.

• They tested different "beat speeds" (different values of $\kappa$).
• What they saw:
  • Slow beats ($\kappa \le -2$): The flow is stable and calm far away from the wall.
  • Medium beats ($-2 < \kappa \le -1.5$): The flow starts to wiggle and oscillate. This is the "sweet spot" where vortices are shed regularly.
  • Fast beats ($\kappa > -1.5$): The flow gets weird. The velocity drops sharply, hinting that the flow is about to break down into chaos (turbulence).

6. Why Does This Matter? (The Real World Connection)

You might wonder, "Who cares about math behind a flat plate?"

• The Cylinder Connection: The authors connect this flat plate math to a cylinder (like a bridge pillar or a smokestack).
• The Vibration Danger: When wind or water flows past a cylinder, it sheds vortices. If these vortices shed at the same frequency as the natural vibration of the structure, resonance occurs.
• The Analogy: Think of a child on a swing. If you push the swing at the exact right moment, it goes higher and higher. If the wind pushes a bridge at the right moment (matching the Strouhal number), the bridge can shake violently and potentially collapse.
• The paper helps engineers predict when this dangerous shaking will happen by linking the math of the flat plate to the real-world Strouhal number of a cylinder.

Summary

This paper is a detective story about fluid flow.

1. The Crime: The flow behind an object is messy and creates swirling vortices.
2. The Clue: The "Strouhal number" (the rhythm of the swirl).
3. The Verdict: The math only works perfectly at specific rhythms. At some rhythms, the flow is impossible; at others, it becomes chaotic.
4. The Lesson: By understanding these mathematical limits, we can better predict when structures (like bridges or chimneys) might start shaking dangerously due to the wind or water flowing past them.

The authors successfully found a "perfect dance" (an exact analytical solution) for a specific rhythm, proving that even in the chaotic world of fluid dynamics, there are moments of perfect mathematical order.

Technical Summary

1. Problem Statement

The paper investigates the unsteady, two-dimensional flow of an incompressible viscous fluid at the rear stagnation point of a flat plate. Unlike the classical forward stagnation point (Hiemenz flow), which yields a steady solution, rear stagnation-point flows are characterized by flow separation, reverse flow, and the formation of vortex sheets.

• Core Challenge: Rear stagnation flows generally lack analytic solutions in two dimensions. Previous models (e.g., Proudman and Johnson) relied on neglecting viscous forces far from the wall to obtain asymptotic solutions.
• Objective: To derive an exact similarity solution for the unsteady Navier-Stokes equations at the rear stagnation point, analyze the conditions under which solutions exist or fail, and establish a physical link between the similarity parameter and the Strouhal number (vortex shedding frequency).

2. Methodology

Governing Equations and Similarity Transformation

The study begins with the conservative form of the Navier-Stokes equations and the continuity equation. A stream function $\psi$ is introduced to satisfy continuity.

• Similarity Variables: The authors adopt unsteady similarity variables proposed by Proudman and Johnson:
  $$\psi = -\sqrt{A(t)\nu x} f(\eta, \tau), \quad \eta = \sqrt{\frac{A(t)}{\nu}} y, \quad \tau = At$$
  where $A(t)$ is a time-dependent flow rate parameter, $\nu$ is kinematic viscosity, and $\eta$ is the similarity variable.
• Derivation: Substituting these into the momentum equations reduces the partial differential equations (PDEs) to a third-order nonlinear ordinary differential equation (ODE):
  $$f''' - ff'' - 1 + (f')^2 = \kappa \left( f' + \frac{1}{2}\eta f'' - 1 \right)$$
  Here, $\kappa = \dot{A}/A^2$ is a constant similarity parameter. The paper identifies $\kappa$ as being equivalent to the Strouhal number ($St$), representing the ratio of local inertial forces to convective inertial forces.

Analytical Approach

• Autonomous Transformation: To solve the ODE analytically, the equation is transformed into an autonomous system by introducing a new variable $F = f + \frac{\kappa}{2}\eta$.
• Case Analysis: The authors analyze specific values of $\kappa$ that allow the equation to become autonomous (independent of the explicit variable $\eta$), specifically $\kappa = -2$ and $\kappa = 2/3$.

Numerical Approach

• Runge-Kutta Method: For general values of $\kappa$, the third-order ODE is reduced to a system of three first-order ODEs.
• Implementation: The system is solved numerically using the MATLAB ode23 solver (an explicit Runge-Kutta method) to generate velocity profiles ($f'$), stream functions ($f$), and shear stress ($f''$) for various $\kappa$ values.

3. Key Contributions and Theoretical Findings

Insolubility at $\kappa = 0$

The paper rigorously proves that no solution exists when $\kappa = 0$ (steady flow case) that satisfies the boundary conditions ($f'(0)=0, f'(\infty)=1$).

• Proof Logic: Using Taylor series expansion and derivative analysis, the authors demonstrate that if a solution approaches 1 at infinity, the derivatives force the function to diverge or violate the no-slip condition at the wall. This confirms that steady rear stagnation flow is mathematically insoluble in 2D.

Exact Analytical Solutions

• Case $\kappa = 2/3$: The derived solution leads to a contradiction where the velocity approaches a negative constant at infinity, violating the physical boundary condition. Thus, no valid solution exists for this case.
• Case $\kappa = -2$: This is the primary analytical breakthrough. The authors derive an exact, closed-form solution:
  $$f'(\eta) = 1 - \cos(\sqrt{3}\eta) \pm i \sin(\sqrt{3}\eta)$$
  This solution is periodic and complex-valued. Physically, it implies that the flow field remains unchanged at large distances (satisfying the potential flow condition) and describes a periodic velocity directed toward the wall. This validates the Proudman-Johnson assumption of neglecting viscosity far from the wall for this specific unsteady regime.

Numerical Regimes

Numerical simulations reveal distinct flow behaviors based on the Strouhal number $\kappa$:

1. $\kappa \le -2$: The similarity solution stabilizes and does not change significantly far from the wall.
2. $-2 < \kappa \le -1.5$: The solution exhibits periodic patterns throughout the flow region. This regime corresponds to pronounced vortex shedding.
3. $\kappa > -1.5$: The solution approaches a linear function ($F(\eta) \approx \eta$) but exhibits a sudden, sharp decrease in velocity near the wall, leading to a singularity or velocity discontinuity. This suggests a breakdown of the laminar similarity model, potentially indicating the onset of turbulence.

4. Physical Significance and Applications

Vortex Shedding and Instability

The paper establishes a direct link between the mathematical parameter $\kappa$ and the physical phenomenon of vortex shedding.

• The range $-2 < \kappa \le -1.5$ corresponds to the unsteady oscillatory solutions where vortices are shed from the wall.
• This shedding creates alternating forces that can induce wall vibrations. If the shedding frequency matches the structural resonance frequency, it can lead to catastrophic failure (e.g., in cylinders or plates).

Explicit Relation to Strouhal Number

The authors derive a novel explicit relation between the similarity parameter $\kappa$ and the classical Strouhal number ($St$):
$$St = -\frac{1.3}{\pi} \kappa$$

• Validation: For the analytical case $\kappa = -2$, the calculated $St \approx 0.8276$. After applying a 3D correction factor (approx. 0.25) to account for cylindrical geometry, the result is $St \approx 0.207$.
• Comparison: This matches experimental data for subcritical cylindrical flow ($Re = 10^3 - 10^5$), where $St \approx 0.20 - 0.21$. This validates the 2D similarity model as a robust predictor for 3D vortex shedding frequencies.

5. Conclusion

The paper successfully demonstrates that:

1. Steady rear stagnation-point flow in 2D is insoluble ($\kappa=0$).
2. Unsteady flow admits an exact analytical solution only for specific parameters (specifically $\kappa = -2$), which describes periodic vortex shedding.
3. The similarity parameter $\kappa$ serves as a proxy for the Strouhal number, allowing the prediction of vortex shedding frequencies and the identification of flow regimes (laminar, oscillatory, or singular/turbulent).
4. The derived relation $St = -1.3\kappa/\pi$ bridges the gap between theoretical similarity solutions and experimental fluid dynamics, confirming the model's utility in predicting flow separation and instability in engineering applications.