Potential Flow Theory
Theory and Physics
What is Potential Flow?
Professor, what is potential flow? It's also called irrotational flow, right?
Potential flow is an inviscid, irrotational (zero vorticity) flow. If the vorticity is zero, the velocity field can be expressed as the gradient of a scalar potential $\phi$.
Substituting this into the incompressible continuity equation $\nabla \cdot \mathbf{u} = 0$ yields the Laplace equation.
The Laplace equation appears in many different fields, doesn't it?
Exactly. Many physical phenomena, such as electrostatic fields, steady-state heat conduction, and groundwater flow, reduce to the Laplace equation. The theory of potential flow is mathematically equivalent to these fields.
Basic Flow Elements
Since the Laplace equation is linear, its solutions can be superimposed, right?
Excellent point. That's the greatest strength of potential flow theory. Complex flows can be constructed by superimposing basic flow elements.
| Flow Element | Velocity Potential $\phi$ | Stream Function $\psi$ | Physical Meaning |
|---|---|---|---|
| Uniform Flow (x-direction) | $U_\infty x$ | $U_\infty y$ | Freestream flow far away |
| Source (strength $m$) | $\frac{m}{2\pi}\ln r$ | $\frac{m}{2\pi}\theta$ | Fluid emission from a point |
| Vortex (circulation $\Gamma$) | $\frac{\Gamma}{2\pi}\theta$ | $-\frac{\Gamma}{2\pi}\ln r$ | Rotation around a point |
| Doublet (strength $\mu$) | $-\frac{\mu \cos\theta}{2\pi r}$ | $-\frac{\mu \sin\theta}{2\pi r}$ | Limit of source + sink |
What happens if you combine a source and a uniform flow?
You get a Rankine half-body. Placing a source of strength $m$ in a uniform flow $U_\infty$ creates a stagnation point at $(x,y) = (-m/(2\pi U_\infty), 0)$, and the streamline passing through it forms the body surface.
Flow Around a Cylinder
The most famous example is flow around a cylinder, right?
The potential flow around a cylinder of radius $a$ in a uniform flow $U_\infty$ is obtained by superimposing a doublet and a uniform flow.
The surface velocity (at $r=a$) is $u_\theta = -2U_\infty \sin\theta$, reaching a maximum value of $2U_\infty$ at $\theta = \pi/2$ (the top).
You can also get the pressure distribution using Bernoulli's theorem, right?
The pressure coefficient becomes $C_p = 1 - 4\sin^2\theta$. This leads to the famous d'Alembert's paradox. Since the pressure distribution is symmetric fore and aft, drag becomes zero in inviscid potential flow.
Kutta-Joukowski Theorem and Lift
Drag is zero, but can lift be generated?
Adding circulation $\Gamma$ around the cylinder generates lift. Adding a vortex to the velocity potential:
According to the Kutta-Joukowski theorem, the lift per unit span is
The magnitude of the circulation determines the lift. For airfoils, the Kutta condition (the condition that the flow leaves the trailing edge smoothly) uniquely determines the value of the circulation.
This theorem is the foundation of airfoil design, right?
Exactly. Using the Joukowski transformation, one can analytically find the flow around an airfoil from the cylinder solution. It's a theory that can be considered the starting point of aeronautical engineering.
Scope of Application and Limitations
In what cases is potential flow theory effective?
It provides a good approximation under the following conditions.
- High Reynolds Number: When viscous effects are confined to the boundary layer
- Regions Far from the Body: Vorticity is nearly zero outside the boundary layer
- Flows Without Separation: Such as airfoils at small angles of attack (before stall)
- Steady or Quasi-Steady: Conditions without vortex shedding
Conversely, it is not applicable to flows involving separation, low Reynolds numbers, or strongly unsteady vortex flows. In actual CFD, the Navier-Stokes equations are solved, but potential flow theory still plays an active role in rapid evaluation during the initial design stage and in validating the reasonableness of CFD results.
D'Alembert's Paradox—The Contradiction of "Zero Drag"
In the 18th century, Jean le Rond d'Alembert used potential flow theory to prove that "the drag on a body moving in an ideal fluid is zero." This is called "d'Alembert's paradox." In reality, it is obviously not zero—this contradiction remained unsolved for over 100 years, hindering the development of fluid mechanics. The solution was Prandtl's boundary layer theory (1904). The answer is that "even a tiny amount of viscosity changes the pressure distribution through a thin boundary layer, generating drag." The behavior at the limit (viscosity→0) is completely different from the solution for zero viscosity—this concept of "singular perturbation" remains an important theme in mathematics and physics today.
Physical Meaning of Each Term
- Time Term $\partial(\rho\phi)/\partial t$: Think of the moment you turn on a faucet. At first, the water comes out in an unstable, splashing manner, but after a while, it becomes a steady flow, right? This "period of change" is described by the time term. The pulsation of blood flow due to the heartbeat, the flow fluctuations each time an engine valve opens and closes—all are unsteady phenomena. So what is steady-state analysis? It looks only at "after sufficient time has passed and the flow has settled down"—in other words, it sets this term to zero. Since this drastically reduces computational cost, solving first under steady-state conditions is a basic CFD strategy.
- Convection Term $\nabla \cdot (\rho \mathbf{u} \phi)$: What happens if you drop a leaf into a river? It gets carried downstream by the flow, right? This is "convection"—the effect where fluid motion transports things. The warm air from a heater reaching the far corner of a room is also because the "carrier," air, transports heat via convection. Here's the interesting part—this term contains "velocity × velocity," making it nonlinear. That is, as the flow becomes faster, this term rapidly strengthens, making control difficult. This is the root cause of turbulence. A common misconception: "Convection and conduction are similar things" → They are completely different! Convection is carried by flow, conduction is transmitted by molecules. There is an order of magnitude difference in efficiency.
- Diffusion Term $\nabla \cdot (\Gamma \nabla \phi)$: Have you ever put milk in coffee and left it? Even without stirring, after a while it naturally mixes, right? That's molecular diffusion. Now a question—honey and water, which flows more easily? Obviously water, right? Honey has high viscosity ($\mu$), so it flows poorly. When viscosity is large, the diffusion term becomes strong, and the fluid moves in a "thick" manner. In low Reynolds number flows (slow, viscous), diffusion is dominant. Conversely, in high Re number flows, convection overwhelms, and diffusion plays a supporting role.
- Pressure Term $-\nabla p$: When you push the plunger of a syringe, the liquid shoots out forcefully from the needle tip, right? Why? Because the plunger side is high pressure, the needle tip is low pressure—this pressure difference provides the force that pushes the fluid. Dam discharge works on the same principle. On a weather map, where isobars are densely packed? That's right, strong winds blow. "Where there is a pressure difference, flow is generated"—this is the physical meaning of the pressure term in the Navier-Stokes equations. A point of confusion here: "Pressure" in CFD is often gauge pressure, not absolute pressure. If results become strange immediately after switching to compressible analysis, it might be due to confusion between absolute/gauge pressure.
- Source Term $S_\phi$: Warmed air rises—why? Because it becomes lighter (lower density) than its surroundings, so it is pushed upward by buoyancy. This buoyancy is added to the equation as a source term. Other examples: chemical reaction heat generated by a gas stove flame, Lorentz force acting on molten metal in a factory's electromagnetic pump... These are all actions that "inject energy or force into the fluid from the outside" and are expressed by source terms. What happens if you forget the source term? In natural convection analysis, if you forget to include buoyancy, the fluid doesn't move at all—you get a physically impossible result, like turning on a heater in a winter room but the warm air doesn't rise.
Assumptions and Applicability Limits
- Continuum Assumption: Valid for Knudsen number Kn < 0.01 (mean free path of molecules ≪ characteristic length)
- Newtonian Fluid Assumption: Shear stress and strain rate have a linear relationship (non-Newtonian fluids require viscosity models)
- Incompressibility Assumption (for Ma < 0.3): Treat density as constant. For Mach numbers above 0.3, compressibility effects must be considered
- Boussinesq Approximation (Natural Convection): Consider density changes only in the buoyancy term, using constant density in other terms
- Non-Applicable Cases: Rarefied gases (Kn > 0.1), supersonic/hypersonic flows (shock wave capturing required), free surface flows (requires VOF/Level Set, etc.)
Dimensional Analysis and Unit Systems
| Variable | SI Unit | Notes / Conversion Memo |
|---|---|---|
| Velocity $u$ | m/s | When converting from volumetric flow rate for inlet conditions, pay attention to cross-sectional area units |
| Pressure $p$ | Pa | Distinguish between gauge pressure and absolute pressure. Use absolute pressure for compressible analysis |
| Density $\rho$ | kg/m³ | Air: approx. 1.225 kg/m³ @20°C, Water: approx. 998 kg/m³ @20°C |
| Viscosity Coefficient $\mu$ | Pa·s | Be careful not to confuse with kinematic viscosity coefficient $\nu = \mu/\rho$ [m²/s] |
| Reynolds Number $Re$ | Dimensionless | $Re = \rho u L / \mu$. Indicator for laminar/turbulent transition |
| CFL Number | Dimensionless | $CFL = u \Delta t / \Delta x$. Directly related to time step stability |
Numerical Methods and Implementation
Basics of Panel Method
How do you solve potential flow numerically?
The most widely used method is the Panel Method. The body surface is divided into panels (line segments or surface elements), and singularities (sources, doublets, vortices) are placed on each panel. It's a method that discretizes the boundary integral equation to find the strength of the singularities.
Here $G$ is the Green's function (in 2D, $G = -\frac{1}{2\pi}\ln|\mathbf{x}-\mathbf{x'}|$), $\sigma$ is the source strength, and $\mu$ is the doublet strength.
Compared to solving the N-S equations in 3D space, the computation is completed only on the surface, right?
That's the biggest advantage of the panel method. A 3D problem reduces to a 2D surface problem, drastically reducing computational cost. Mesh generation also only requires a surface mesh.
Hess-Smith Panel Method
Please teach me the most basic panel method algorithm.
The Hess-Smith Panel Method is used for flow around bodies without lift. The algorithm is as follows.
1. Divide the body surface into $N$ panels
2. Place a constant-strength source $\sigma_j$ on each panel
3. Impose the condition that the normal velocity = 0 at each panel's control point (midpoint)
4. Solve the $N \times N$ system of equations $[A]\{\sigma\} = \{b\}$
5. Calculate the tangential velocity on each panel and find the pressure using Bernoulli's equation
The influence coefficient matrix $A_{ij}$ is the normal velocity component induced at panel $i$'s control point by the source on panel $j$.
How do you calculate lift?
Add vortex panels to Hess-Smith. Add a vortex distribution $\gamma$ to each panel and add the Kutta condition (equal velocities on the upper and lower panels at the trailing edge) to the system of equations. This yields the circulation and lift.
Higher-Order Panel Methods
How can accuracy be improved?
The basic Hess-Smith assumes a constant (zeroth-order) singularity distribution on each panel, but there are the following improvements.
| Panel Method Type | Singularity Distribution | Accuracy | Computational Cost |
|---|---|---|---|
| Constant Panel (Hess-Smith) | Constant | 1st order | Low ($O(N^2)$) |
| Linear Panel | Linear distribution | 2nd order | Medium ($O(N^2)$) |
| Quadratic Panel | Quadratic distribution | 3rd order | High ($O(N^2)$) |
| High-Order Panel + FMM | Arbitrary | High order | $O(N \log N)$ or $O(N)$ |
So computational cost can be reduced with FMM.
Combining with FMM (Fast Multipole Method), which calculates the influence of distant panels collectively using multipole expansions, can reduce computational cost to $O(N)$. This is an essential technique for 3D full-aircraft analysis where the number of panels reaches tens to hundreds of thousands.
Representative Panel Method Codes
Related Topics
なった
詳しく
報告