Visual walkthrough — Computational aerodynamics — panel method (intro), CFD overview

3.1.30 · D2 · Physics › Compressible Flow & Aerodynamics › Computational aerodynamics — panel method (intro), CFD overv

---

Step 1 — "Flow" ka matlab picture mein kya hota hai

KYA karte hain hum: "hawa kisi tarah move kar rahi hai" ko "har point par ek arrow $\vec{V}$ hai" se replace karo.

KYUN: ek wing ko forces sirf isliye feel hote hain kyunki us ke paas ke hawa ke arrows mur jaate hain aur speed up ho jaate hain. Agar hum arrows jaante hain, hum sab kuch jaante hain — pressure, lift, sab kuch.

PICTURE: neeche, door ka uniform wind identical arrows ka ek field hai jo sab ek hi direction mein $V_{\infty }$ speed se point kar rahe hain (padho "$V$-infinity" = wing se door hawa ki speed, jahaan body ka abhi koi influence nahi hai).

$\vec{V}=(u,v)$

• $u$ ::: arrow ka sideways (horizontal, $x$) waala hissa
• $v$ ::: arrow ka up-down (vertical, $y$) waala hissa
• $\vec{V}$ ::: poora arrow, in dono hisson se bana hua

---

Step 2 — Woh ek rule jo wing impose karta hai

"Seedha surface ke andar" ki baat karne ke liye humein ek aur picture-tool chahiye.

Ab, "the part of $\vec{V}$ pointing along $\hat{n}$" nikalne ke liye hum dot product use karte hain.

KYA karte hain hum: solid-wall rule ko ek clean equation mein likhna.

KYUN: yeh single condition shape ka poora fingerprint hai. Aage jo kuch bhi hai woh sirf isi ko satisfy karne ke liye exist karta hai.

PICTURE: wind arrow $\vec{V}$ do hisson mein toot ta hai — wall ke saath waala hissa (allowed, blue) aur wall ke andar waala hissa (forbidden, pink). Hum pink piece ko khatam karna chahte hain.

$\boxed{\vec{V}\cdot \hat{n}=0}\text{on the body}$

• $\hat{n}$ ::: skin ke us point par outward perpendicular unit arrow
• $\vec{V}\cdot \hat{n}$ ::: wall ke andar ja rahi flow ki matra — yeh zero honi chahiye
• $=0$ ::: koi tunnelling nahi; yeh flow-tangency condition hai

Lekin akeli raw wind isko satisfy nahi karti — pink piece non-zero hai. Humein kuch aisa add karna hoga jo ise cancel kare. Enter the sources.

---

Step 3 — Spray nozzle (ek source) aur yeh kyun help karta hai

KYA karte hain hum: in nozzles ko wall ke paas rakhte hain taaki unka outward push wind ke forbidden into-the-wall component ko cancel kar sake.

YEH TOOL KYUN, KUCH AUR NAHI: ek source ka field, wind ke field mein add karne par, phir bhi ek legal flow hai (dono same governing equation follow karte hain — dekho Laplace's Equation & Potential Flow). Aur ek source jo baahaar push karta hai woh bilkul wahi cheez hai jo wind ke andar push karne ko neutralise kar sakta hai. Vortices flow ko spin karte hain; sinks kheenchte hain; source ek natural "anti-tunnel" tool hai.

PICTURE: ek akela source; uske arrows baahaar radiate karte hain, beech mein lambe, door chhote.

Strength $\lambda$ ka source distance $r$ par jo velocity banata hai: $V_{\text{source}}=\frac{\lambda }{2\pi r}$

• $\lambda$ ::: strength — kitni zor se spray karta hai (bada = zyada push)
• $r$ ::: source se us jagah tak ki distance jahaan hum measure kar rahe hain
• $2\pi$ ::: woh poora circle jis par yeh flow spread karta hai; isse divide karne par flow evenly around share hoti hai

---

Step 4 — Outline ko panels mein kaato

KYA karte hain hum: smooth curve ko $N$ flat panels se replace karo; unknowns $N$ numbers ${\lambda }_1,\dots ,{\lambda }_N$ ban jaate hain.

KYUN: $N$ unknown numbers → ek finite, solvable problem. Zyada panels = better shape, lekin (parent ke mistake box ka foreshadowing) cost tezi se badhti hai.

PICTURE: neeche ek smooth airfoil, phir wahi airfoil straight panels ke roop mein. Har panel par: beech mein ek red dot — control point — jahaan hum Step 2 ka rule enforce karenge. Har panel ka apna outward ${\hat{n}}_i$ hai.

• $N$ ::: humne kitne panels choose kiye
• ${\lambda }_j$ ::: panel $j$ par constant source strength (ek unknown per panel)
• control point ::: panel ka midpoint, jahaan flow-tangency test ki jaati hai
• ${\hat{n}}_i$ ::: panel $i$ ka outward normal

---

Step 5 — Sab kuch add karo (superposition)

KYA karte hain hum: kisi bhi point $P$ par total flow = wind + har panel ke source ka push.

KYUN: hum exactly wahi flow build kar sakte hain jo humein chahiye tuned ingredients add karke, jaise paint mix karna.

PICTURE: teen stacked layers — wind arrows, source arrows, aur unka sum — yeh dikhate hue ki combined arrows surface ko hug karne ke liye kaise curve karte hain.

$\vec{V}(P)=\underset{\text{wind}}{\underbrace{{\vec{V}}_{\infty }}}+\underset{\text{all panels}}{\underbrace{\sum _{j=1}^N{\lambda }_j{\vec{v}}_j(P)}}$

• ${\vec{V}}_{\infty }$ ::: uniform far-field wind (fixed, known)
• ${\sum }_{j=1}^N$ ::: "har panel $j$ ko $1$ se $N$ tak add karo"
• ${\lambda }_j$ ::: woh unknown strength jiske liye hum solve kar rahe hain
• ${\vec{v}}_j(P)$ ::: woh flow jo panel $j$ par ek unit-strength source point $P$ par banata hai (pure geometry — ek baar compute ho sakta hai)

---

Step 6 — Rule ko equations mein badlo (influence coefficients)

Ab Step 2 ka rule apply karo — zero into-wall flow — har control point $i$ par. Poore sum ko us panel ke ${\hat{n}}_i$ ke saath dot karo:

${\vec{V}}_{\infty }\cdot {\hat{n}}_i+{\sum }_{j=1}^N{\lambda }_j({\vec{v}}_j\cdot {\hat{n}}_i)=0$

KYA karte hain hum: har control point par ek flow-tangency equation → $N$ equations.

KYUN: $N$ unknowns ${\lambda }_j$ ke liye $N$ equations — ab genuinely solvable.

PICTURE: source $j$ se control point $i$ tak jaata ek arrow, aur sirf uska ${\hat{n}}_i$ ke along wala component (woh hissa jo hum rakhte hain) highlighted. Woh rakhaa gaya hissa hi $A_{ij}$ hai.

$b_i=-{\vec{V}}_{\infty }\cdot {\hat{n}}_i$ likhne par (woh into-wall wind jo humein cancel karni hai), har equation yeh hai:

${\sum }_{j=1}^N{A}_{ij}{\lambda }_j=b_i$

• $A_{ij}$ ::: panel $j$ ka panel $i$ par geometry-only influence
• ${\lambda }_j$ ::: unknown strengths
• $b_i$ ::: panel $i$ par into-wall wind, minus sign ke saath (hume ise undo karna hai)

---

Step 7 — Sabko ek matrix equation mein stack karo

KYA karte hain hum: $N$ equations ko matrix ki rows mein pack karo.

KYUN: "$A\mathbf{\lambda }=\mathbf{b}$ solve karo" computing ka sabse standard kaam hai — machine $A$ ko invert karti hai aur ek saath har ${\lambda }_j$ padh leti hai.

PICTURE: grid $A$, lamba $\mathbf{\lambda }$, aur lamba $\mathbf{b}$, panel $i$ ki row highlighted taaki tum dekh sako ki yeh exactly Step 6 ki equation hai.

Solve hone ke baad, har panel par along-wall (tangential) speed $V_{t,i}$ nikalo aur pressure Bernoulli's Equation se padho: $C_{p,i}=1-{\left(\frac{V_{t,i}}{V_{\infty }}\right)}^2$

• $C_{p,i}$ ::: panel $i$ par pressure coefficient (dimensionless pressure)
• $V_{t,i}/V_{\infty }$ ::: wahan hawa free wind se kitni tezi/dheemi slide karti hai
• fast flow ($V_t>V_{\infty }$) ::: low pressure (lift); slow flow ::: high pressure

---

Step 8 — Degenerate case: sirf sources se ZERO lift milti hai

KYA karte hain hum: notice karo ki sources body ko shape dete hain lekin lift nahi karte — aur dekho ki yeh ek special (degenerate) outcome hai, koi bug nahi.

KYUN: yeh real-world fix ko force karta hai — vortex panels add karo aur Kutta condition impose karo (flow ko sharp trailing edge se smoothly nikalna chahiye), jo ek non-zero $\Gamma$ set karta hai, isliye lift milti hai. Yeh zero-drag-aur-yahan-zero-lift ideal d'Alembert's Paradox in action hai; real drag ke liye Boundary Layer Theory & Skin Friction Drag chahiye.

PICTURE: same airfoil ke upar do flows — top (sirf source) symmetric hai, streamlines front to back mirror karte hain, $\Gamma =0$; bottom (source + vortex + Kutta) tilted hai, hawa trailing edge se cleanly nikalti hai, $\Gamma \ne 0$, lift upar.

• $\Gamma$ ::: circulation — body ke around hawa ka net swirl
• $L^{\prime }=\rho V_{\infty }\Gamma$ ::: lift per unit span; koi swirl nahi ⇒ koi lift nahi
• Kutta condition ::: $\Gamma$ ko pin karta hai taaki flow trailing edge ke around wrap kiye bina nikle

---

Ek-picture summary

Poora safar ek board par: wind shape se takraati hai → outline ko tunable sources se dhako → har midpoint par zero into-wall flow demand karo → yeh $A\mathbf{\lambda }=\mathbf{b}$ ban jaata hai → strengths solve karo → pressure padho → (lift ke liye vortices add karo).

Recall Feynman retelling — seedhe words mein bolo

Socho ki hawa ek solid shape par blow kar rahi hai. Hawa ko metal ke andar jaana allowed nahi hai, sirf us ke saath-saath. Ise force karne ke liye, main outline ko bahut saare chhote spray nozzles se line karta hoon aur har ek ko ek aisi strength par spray karne deta hoon jo main choose kar sakta hoon. Agar main sprays ko sahi aim karoon, unka outward push exactly us hawa ke hisse ko cancel kar deta hai jo wall ko pierce karne ki koshish kar raha tha — toh total flow bas surface ko skim karta hai. "Sahi aim karna" ka matlab hai: har chhote segment ke beech mein, sideways-into-wall flow zero tak add honi chahiye. Yeh ek equation hai per segment. Har nozzle ka har midpoint par push ek fixed geometry number hai (ek influence coefficient), isliye saari equations milke bas ek grid ki rows hain: $A$ times unknown strengths ki list, winds ki list ke barabar jo mujhe cancel karni hain. Computer grid ko invert karta hai, mujhe har strength deta hai, aur along-wall speeds se main Bernoulli ke zariye har jagah pressure padh leta hoon. Ek catch: pure sprays symmetric hain aur lift nahi kar sakti — lift ke liye main chhote swirls (vortices) bhi add karta hoon aur require karta hoon ki hawa sharp tail se smoothly peel off kare. Sprays shape banate hain; swirls lift banate hain.

Recall Quick self-test

• Source strengths ko pin karne wala ek physical rule kya hai? → Har control point par zero normal (into-wall) velocity.
• "Into-wall flow" nikalne ke liye kaun sa math tool use hota hai aur kyun? → Dot product $\vec{V}\cdot \hat{n}$; yeh $\hat{n}$ ke along component ko isolate karta hai.
• Hum wind + sources freely kyun add kar sakte hain? → Superposition: governing equation linear hai.
• $A_{ij}$ kya hai? → Panel $j$ par unit source panel $i$ ke midpoint par kitni into-wall velocity banata hai (pure geometry).
• Sirf sources se zero lift kyun, aur fix kya hai? → Koi circulation nahi ($\Gamma =0$); vortex panels + Kutta condition add karo.