Intuition The big picture (WHY this exists)
A 2D incompressible flow has a velocity field V ⃗ = ( u , v ) \vec{V} = (u,v) V = ( u , v ) . That's two functions of position. But these two functions are not independent — they are tied together by physics:
Incompressibility (mass conservation, ∇ ⋅ V ⃗ = 0 \nabla\cdot\vec V=0 ∇ ⋅ V = 0 ) ties them → lets us pack them into ONE scalar, the ==stream function ψ \psi ψ ==.
Irrotationality (no spin, ∇ × V ⃗ = 0 \nabla\times\vec V=0 ∇ × V = 0 ) ties them → lets us pack them into ONE scalar, the ==velocity potential ϕ \phi ϕ ==.
So instead of solving for two velocity components, we solve for one scalar field . That's the whole point: fewer unknowns, easier math.
Definition Velocity potential
If a flow is irrotational (∇ × V ⃗ = 0 \nabla\times\vec V = 0 ∇ × V = 0 ), there exists a scalar ϕ \phi ϕ such that
V ⃗ = ∇ ϕ ⇒ u = ∂ ϕ ∂ x , v = ∂ ϕ ∂ y . \vec V = \nabla\phi \qquad\Rightarrow\qquad u=\frac{\partial\phi}{\partial x},\quad v=\frac{\partial\phi}{\partial y}. V = ∇ ϕ ⇒ u = ∂ x ∂ ϕ , v = ∂ y ∂ ϕ .
Such a flow is called potential flow .
WHY does ϕ \phi ϕ exist? From vector calculus: the curl of a gradient is always zero, ∇ × ( ∇ ϕ ) = 0 \nabla\times(\nabla\phi)=0 ∇ × ( ∇ ϕ ) = 0 . So if a field can be written as a gradient, it is automatically irrotational. The converse (a key theorem) says: in a simply-connected region, if ∇ × V ⃗ = 0 \nabla\times\vec V=0 ∇ × V = 0 , then V ⃗ \vec V V can be written as ∇ ϕ \nabla\phi ∇ ϕ .
Intuition Why is it like the electric potential?
Same math as E ⃗ = − ∇ V \vec E=-\nabla V E = − ∇ V . Velocity flows "downhill" along ϕ \phi ϕ . The sign convention in fluids is usually + + + (no minus), so velocity points toward increasing ϕ \phi ϕ .
Definition Stream function (2D incompressible)
If a flow is incompressible (∇ ⋅ V ⃗ = 0 \nabla\cdot\vec V=0 ∇ ⋅ V = 0 ) and 2D, there exists a scalar ψ \psi ψ such that
u = ∂ ψ ∂ y , v = − ∂ ψ ∂ x . u=\frac{\partial\psi}{\partial y},\qquad v=-\frac{\partial\psi}{\partial x}. u = ∂ y ∂ ψ , v = − ∂ x ∂ ψ .
WHY this works (DERIVE it): Continuity in 2D is
∂ u ∂ x + ∂ v ∂ y = 0. \frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0. ∂ x ∂ u + ∂ y ∂ v = 0.
We want u , v u,v u , v to come from one function so this is automatically satisfied. Try u = ∂ ψ / ∂ y u=\partial\psi/\partial y u = ∂ ψ / ∂ y , v = − ∂ ψ / ∂ x v=-\partial\psi/\partial x v = − ∂ ψ / ∂ x . Substitute:
∂ ∂ x ( ∂ ψ ∂ y ) + ∂ ∂ y ( − ∂ ψ ∂ x ) = ψ x y − ψ y x = 0. \frac{\partial}{\partial x}\!\left(\frac{\partial\psi}{\partial y}\right)+\frac{\partial}{\partial y}\!\left(-\frac{\partial\psi}{\partial x}\right)=\psi_{xy}-\psi_{yx}=0. ∂ x ∂ ( ∂ y ∂ ψ ) + ∂ y ∂ ( − ∂ x ∂ ψ ) = ψ x y − ψ y x = 0.
It vanishes identically because mixed partials commute. So by construction, any ψ \psi ψ gives a divergence-free flow. That is why the minus sign sits on v v v — to cancel.
Intuition Why streamlines = lines of constant
ψ \psi ψ
Along a streamline, the velocity is tangent to the line: d y d x = v u \frac{dy}{dx}=\frac{v}{u} d x d y = u v , i.e. u d y − v d x = 0 u\,dy - v\,dx = 0 u d y − v d x = 0 .
Now the change in ψ \psi ψ moving by ( d x , d y ) (dx,dy) ( d x , d y ) is
d ψ = ψ x d x + ψ y d y = − v d x + u d y . d\psi=\psi_x\,dx+\psi_y\,dy = -v\,dx + u\,dy. d ψ = ψ x d x + ψ y d y = − v d x + u d y .
This is exactly u d y − v d x u\,dy - v\,dx u d y − v d x ! So d ψ = 0 d\psi=0 d ψ = 0 along a streamline ⇒ streamlines are curves of constant ψ \psi ψ . Beautiful.
If a flow is both incompressible and irrotational ("ideal flow"), both ϕ \phi ϕ and ψ \psi ψ exist.
Each satisfies Laplace's equation. WHY:
Incompressible + V ⃗ = ∇ ϕ \vec V=\nabla\phi V = ∇ ϕ : ∇ ⋅ ∇ ϕ = 0 ⇒ ∇ 2 ϕ = 0 \nabla\cdot\nabla\phi=0\Rightarrow \nabla^2\phi=0 ∇ ⋅ ∇ ϕ = 0 ⇒ ∇ 2 ϕ = 0 .
Irrotational + ψ \psi ψ defs: ∂ v ∂ x − ∂ u ∂ y = 0 ⇒ − ψ x x − ψ y y = 0 ⇒ ∇ 2 ψ = 0 \dfrac{\partial v}{\partial x}-\dfrac{\partial u}{\partial y}=0 \Rightarrow -\psi_{xx}-\psi_{yy}=0\Rightarrow \nabla^2\psi=0 ∂ x ∂ v − ∂ y ∂ u = 0 ⇒ − ψ xx − ψ y y = 0 ⇒ ∇ 2 ψ = 0 .
∇ 2 ϕ = 0 , ∇ 2 ψ = 0 \boxed{\nabla^2\phi=0,\qquad \nabla^2\psi=0} ∇ 2 ϕ = 0 , ∇ 2 ψ = 0
ϕ \phi ϕ -lines ⟂ ψ \psi ψ -lines
∇ ϕ = ( ϕ x , ϕ y ) = ( u , v ) \nabla\phi=(\phi_x,\phi_y)=(u,v) ∇ ϕ = ( ϕ x , ϕ y ) = ( u , v ) is the velocity — perpendicular to equipotential lines.
∇ ψ = ( ψ x , ψ y ) = ( − v , u ) \nabla\psi=(\psi_x,\psi_y)=(-v,u) ∇ ψ = ( ψ x , ψ y ) = ( − v , u ) — perpendicular to streamlines.
Their dot product: ( u ) ( − v ) + ( v ) ( u ) = 0 (u)(-v)+(v)(u)=0 ( u ) ( − v ) + ( v ) ( u ) = 0 . So equipotential lines and streamlines cross at right angles — they form an orthogonal grid (a "flow net").
Worked example Uniform flow
Take ϕ = U x \phi = U x ϕ = U x . Then u = ϕ x = U u=\phi_x=U u = ϕ x = U , v = ϕ y = 0 v=\phi_y=0 v = ϕ y = 0 . Uniform horizontal stream.
Why this step? Differentiating the potential gives velocity directly.
Stream function: u = ψ y = U ⇒ ψ = U y u=\psi_y=U\Rightarrow\psi=Uy u = ψ y = U ⇒ ψ = U y ; v = − ψ x = 0 v=-\psi_x=0 v = − ψ x = 0 ✓. Streamlines ψ = U y = \psi=Uy= ψ = U y = const are horizontal lines. Why? Constant y y y = horizontal, matching the flow direction.
Worked example Source flow (strength
m m m )
Radial outflow u r = m 2 π r u_r=\dfrac{m}{2\pi r} u r = 2 π r m , u θ = 0 u_\theta=0 u θ = 0 . In polar form, ϕ \phi ϕ satisfies u r = ∂ ϕ / ∂ r u_r=\partial\phi/\partial r u r = ∂ ϕ / ∂ r .
Integrate: ϕ = m 2 π ln r \phi=\dfrac{m}{2\pi}\ln r ϕ = 2 π m ln r . Why? Because ∂ ∂ r m 2 π ln r = m 2 π r = u r \frac{\partial}{\partial r}\frac{m}{2\pi}\ln r=\frac{m}{2\pi r}=u_r ∂ r ∂ 2 π m ln r = 2 π r m = u r ✓.
Stream function: u r = 1 r ∂ ψ ∂ θ = m 2 π r ⇒ ψ = m 2 π θ u_r=\frac1r\frac{\partial\psi}{\partial\theta}=\frac{m}{2\pi r}\Rightarrow \psi=\dfrac{m}{2\pi}\theta u r = r 1 ∂ θ ∂ ψ = 2 π r m ⇒ ψ = 2 π m θ .
Streamlines ψ = \psi= ψ = const ⇒ θ = \theta= θ = const ⇒ radial rays out of the origin. Why? A source pushes fluid straight outward, so streamlines are spokes.
Worked example Check orthogonality for the source
Equipotentials ϕ = m 2 π ln r = \phi=\frac{m}{2\pi}\ln r= ϕ = 2 π m ln r = const ⇒ circles (r = r= r = const). Streamlines = radial rays. Circles meet radial rays at 90 ° 90° 90° . Why this matters: confirms the flow-net theorem visually.
Worked example Flux between two streamlines
Between ψ 1 \psi_1 ψ 1 and ψ 2 \psi_2 ψ 2 , Q = ψ 2 − ψ 1 Q=\psi_2-\psi_1 Q = ψ 2 − ψ 1 . For the source between θ = 0 \theta=0 θ = 0 and θ = 2 π \theta=2\pi θ = 2 π : Q = m 2 π ( 2 π ) − 0 = m Q=\frac{m}{2\pi}(2\pi)-0=m Q = 2 π m ( 2 π ) − 0 = m . Why? Total outflow from a source of strength m m m is exactly m m m — internally consistent.
ϕ \phi ϕ exists for every flow."
Why it feels right: Velocity is a vector field, and we love potentials.
The fix: ϕ \phi ϕ exists only if the flow is irrotational (∇ × V ⃗ = 0 \nabla\times\vec V=0 ∇ × V = 0 ). A spinning (vortical) flow has no single-valued ϕ \phi ϕ .
ψ \psi ψ exists for every flow."
Why it feels right: Streamlines always exist, so surely ψ \psi ψ does too.
The fix: The scalar ψ \psi ψ (with the flux property) is guaranteed for 2D incompressible flow. For compressible flow you need a modified (density-weighted) stream function.
Common mistake Forgetting the minus sign:
v = − ψ x v=-\psi_x v = − ψ x .
Why it feels right: Symmetry tempts you to write v = ψ x v=\psi_x v = ψ x .
The fix: The minus sign is required so that continuity cancels (ψ x y − ψ y x = 0 \psi_{xy}-\psi_{yx}=0 ψ x y − ψ y x = 0 ). Drop it and you no longer guarantee incompressibility.
Common mistake Mixing the roles: "
ϕ \phi ϕ =const are streamlines."
Why it feels right: Both are level curves of a scalar.
The fix: ψ = \psi= ψ = const → streamlines (along the flow). ϕ = \phi= ϕ = const → equipotentials (perpendicular to the flow).
Recall Feynman: explain to a 12-year-old
Imagine a river. Instead of tracking how fast water moves in the left–right and up–down directions separately (two things), we use ONE clever number for each spot:
The stream function ψ \psi ψ is like a "lane number." Water never crosses lanes, so any line where ψ \psi ψ stays the same is a path water actually follows. And the gap between two lane numbers tells you exactly how much water flows in that channel.
The velocity potential ϕ \phi ϕ is like "height on a hill" — water slides toward higher ϕ \phi ϕ . It only works when the water isn't spinning in little whirlpools.
When the water is both un-squishable and non-spinning, both numbers exist, and their maps cross like graph paper at perfect right angles.
Mnemonic Remember the definitions
"ψ \psi ψ has the cross, ϕ \phi ϕ goes straight."
ψ \psi ψ : u = ψ y u=\psi_y u = ψ y , v = − ψ x v=-\psi_x v = − ψ x → indices cross and a minus appears.
ϕ \phi ϕ : u = ϕ x u=\phi_x u = ϕ x , v = ϕ y v=\phi_y v = ϕ y → indices match , no minus.
And "STeam = STreamlines" (ψ \psi ψ ), "Potential = Perpendicular Pressure-like" (ϕ \phi ϕ ).
Condition for velocity potential ϕ \phi ϕ to exist Flow must be irrotational,
∇ × V ⃗ = 0 \nabla\times\vec V=0 ∇ × V = 0 (since
∇ × ∇ ϕ = 0 \nabla\times\nabla\phi=0 ∇ × ∇ ϕ = 0 ).
Condition for stream function ψ \psi ψ to exist 2D incompressible flow,
∇ ⋅ V ⃗ = 0 \nabla\cdot\vec V=0 ∇ ⋅ V = 0 .
Define u , v u,v u , v from ϕ \phi ϕ u = ∂ ϕ / ∂ x , v = ∂ ϕ / ∂ y u=\partial\phi/\partial x,\; v=\partial\phi/\partial y u = ∂ ϕ / ∂ x , v = ∂ ϕ / ∂ y .
Define u , v u,v u , v from ψ \psi ψ u = ∂ ψ / ∂ y , v = − ∂ ψ / ∂ x u=\partial\psi/\partial y,\; v=-\partial\psi/\partial x u = ∂ ψ / ∂ y , v = − ∂ ψ / ∂ x .
Why is the minus sign needed in v = − ψ x v=-\psi_x v = − ψ x ? So continuity
ψ x y − ψ y x = 0 \psi_{xy}-\psi_{yx}=0 ψ x y − ψ y x = 0 is satisfied automatically.
What does ψ = \psi= ψ = const represent? A streamline (curve tangent to the velocity).
What does ϕ = \phi= ϕ = const represent? An equipotential line, perpendicular to streamlines.
Flow rate between two streamlines Q = ψ 2 − ψ 1 Q=\psi_2-\psi_1 Q = ψ 2 − ψ 1 (per unit depth).
Equation satisfied by ϕ \phi ϕ and ψ \psi ψ in ideal flow Laplace's equation,
∇ 2 ϕ = 0 , ∇ 2 ψ = 0 \nabla^2\phi=0,\ \nabla^2\psi=0 ∇ 2 ϕ = 0 , ∇ 2 ψ = 0 .
Why are equipotentials ⟂ streamlines? ∇ ϕ ⋅ ∇ ψ = ( u ) ( − v ) + ( v ) ( u ) = 0 \nabla\phi\cdot\nabla\psi=(u)(-v)+(v)(u)=0 ∇ ϕ ⋅ ∇ ψ = ( u ) ( − v ) + ( v ) ( u ) = 0 .
Cauchy–Riemann (fluid) relations ϕ x = ψ y \phi_x=\psi_y ϕ x = ψ y and
ϕ y = − ψ x \phi_y=-\psi_x ϕ y = − ψ x .
ϕ \phi ϕ and ψ \psi ψ for uniform flow U U U ϕ = U x , ψ = U y \phi=Ux,\ \psi=Uy ϕ = U x , ψ = U y .
ϕ \phi ϕ and ψ \psi ψ for a source of strength m m m ϕ = m 2 π ln r , ψ = m 2 π θ \phi=\frac{m}{2\pi}\ln r,\ \psi=\frac{m}{2\pi}\theta ϕ = 2 π m ln r , ψ = 2 π m θ .
Continuity Equation — source of the stream function (∇ ⋅ V ⃗ = 0 \nabla\cdot\vec V=0 ∇ ⋅ V = 0 ).
Vorticity and Irrotational Flow — source of the velocity potential (∇ × V ⃗ = 0 \nabla\times\vec V=0 ∇ × V = 0 ).
Laplace Equation — both ϕ \phi ϕ and ψ \psi ψ are harmonic.
Bernoulli Equation — used once ϕ \phi ϕ gives the velocity field.
Cauchy-Riemann Equations — complex potential w = ϕ + i ψ w=\phi+i\psi w = ϕ + i ψ .
Electrostatic Potential — direct mathematical analogue (E ⃗ = − ∇ V \vec E=-\nabla V E = − ∇ V ).
Flow Nets — orthogonal grid of ϕ \phi ϕ and ψ \psi ψ lines.
auto-satisfies continuity
combined with incompressible
combined with irrotational
Incompressibility div V=0
Intuition Hinglish mein samjho
Dekho, 2D flow mein velocity ke do components hote hain, u u u aur v v v . Do unknown solve karna mushkil hai. Idea simple hai: physics ki do conditions in dono ko jod deti hain, toh hum unhe ek hi scalar function mein pack kar sakte hain. Agar flow incompressible hai (mass conserve, ∇ ⋅ V ⃗ = 0 \nabla\cdot\vec V=0 ∇ ⋅ V = 0 ), toh stream function ψ \psi ψ banti hai: u = ψ y u=\psi_y u = ψ y , v = − ψ x v=-\psi_x v = − ψ x . Ye minus sign zaroori hai taaki continuity automatically satisfy ho jaaye — mixed partials cancel ho jaate hain.
Stream function ka sabse pyaara fact: ψ = \psi= ψ = constant waali lines hi streamlines hain, yaani fluid actually inhi raston pe behta hai. Aur do streamlines ke beech ka difference ψ 2 − ψ 1 \psi_2-\psi_1 ψ 2 − ψ 1 exactly batata hai ki utne channel mein kitna flow Q Q Q ja raha hai. Toh ψ \psi ψ sirf maths nahi, ekdum physical "lane number" hai.
Dusri taraf, agar flow irrotational hai (koi spinning nahi, ∇ × V ⃗ = 0 \nabla\times\vec V=0 ∇ × V = 0 ), toh velocity potential ϕ \phi ϕ exist karta hai, jaise electric potential. V ⃗ = ∇ ϕ \vec V=\nabla\phi V = ∇ ϕ , matlab u = ϕ x u=\phi_x u = ϕ x , v = ϕ y v=\phi_y v = ϕ y . Yaad rakho: ϕ \phi ϕ sirf tab milta hai jab flow ghoomta nahi.
Jab flow dono ho — incompressible aur irrotational (ideal flow) — tab dono ϕ \phi ϕ aur ψ \psi ψ exist karte hain, dono Laplace equation ∇ 2 ϕ = 0 \nabla^2\phi=0 ∇ 2 ϕ = 0 , ∇ 2 ψ = 0 \nabla^2\psi=0 ∇ 2 ψ = 0 follow karte hain, aur unki lines ek dusre se 90 degree pe milti hain (flow net). Exam mein bas yaad rakhna: ψ \psi ψ = streamlines (cross indices + minus), ϕ \phi ϕ = equipotentials (straight indices, no minus). Ye 80/20 ka core hai.