2.2.22 · D3 · Physics › Fluid Mechanics › Blasius solution — exact laminar boundary layer solution
Yeh page Blasius solution ke liye ek worked-example gym hai. Parent note ne theory build ki; yahan hum ise har us situation ke against drill karte hain jo exam (ya reality) de sakta hai. Agar koi symbol naya lage, woh parent mein define tha — lekin hum zaroori cheezein yahan bhi re-anchor karte hain taaki aap kabhi na rukein.
Intuition In saare problems mein kya common hai
Har Blasius question secretly EK hi sawaal hai jo chaar tarike se pucha jaata hai: "slow zone kitna mota hai, drag kitna hai, kis speed par, kis jagah par?" Chaaron ek hi dimensionless number par tike hain, Reynolds number R e x = U ∞ x / ν (dekho Reynolds Number ). Pehle R e x nikalo, aur baaki sab multiplication se mil jaata hai. R e x ko master dial samjho.
Kuch bhi compute karne se pehle, parent ke chaar workhorses pin karo, unke plain meaning ke saath:
Definition Do viscosities, aur unka relation
Do "stickiness" numbers hain, aur problems inhe freely mix karti hain:
ν — kinematic viscosity (m²/s), R e ke andar use hoti hai.
μ — dynamic viscosity (Pa·s = kg/(m·s)), jo velocity gradient se multiply hoke stress deti hai.
Ye density se jude hain: μ = ρ ν . Hamare house air ke liye, μ = 1.2 × 1.5 × 1 0 − 5 = 1.8 × 1 0 − 5 Pa⋅s . Jab bhi neeche μ dikhe, iska matlab bas ρ ν hai.
Example 7 (jo layer ke andar jhaankta hai) se pehle parent ka ek aur idea chahiye. Parent ne dikhaya ki velocity profile har x par identical dikhti hai jab aap height ko clever tarike se measure karo:
Definition Similarity coordinate
η aur profile function f
η = y U ∞ / ( ν x ) — stretched wall-distance . Yeh real height y ko rescale karta hai taaki ek universal curve har x par fit ho. Wall par η = 0 ; layer effectively η ≈ 5 tak khatam ho jaati hai.
f ( η ) — ek dimensionless helper (scaled stream function) jiska slope velocity carry karta hai.
f ′ ( η ) — iski slope, aur f ′ ( η ) = u / U ∞ , velocity ratio . Toh wall par f ′ = 0 (no-slip) aur door f ′ → 1 (free stream).
Ye numerically-solved Blasius equation (parent note) se aate hain. Hume us table se sirf do read-offs chahiye: η ≈ 1.6 par f ′ ≈ 0.5 aur f ≈ 0.65 ; aur f ′ = 0.99 at η ≈ 5.0 . Kisi η ko real height mein wapas convert karne ke liye, definition invert karo: y = η ν x / U ∞ .
Har Blasius problem inhi cells mein se kisi ek mein aata hai. Neeche ke examples cell(s) ke label ke saath hain — milke yeh poora grid cover karte hain.
Cell
Kya vary karta hai / degenerate input
Question type
Example
A
Standard interior point, laminar
δ nikalo
Ex 1
B
Same point
τ w nikalo (local force)
Ex 2
C
Poori plate
Total drag F D nikalo
Ex 3
D
Limit x → 0 (leading edge)
Degenerate: δ → 0 , τ w → ∞
Ex 4
E
Ratio / scaling — x ya U ∞ badlo, numbers nahi
δ , τ w ka response kaisa?
Ex 5
F
Validity ki boundary : R e x near 5 × 1 0 5
Kya Blasius allowed bhi hai? Transition x nikalo
Ex 6
G
Layer ke andar — koi specific y , ya v nikalo
Profile / vertical velocity ≠ 0
Ex 7
H
Real-world word problem (alag fluid: water)
Full pipeline
Ex 8
I
Exam twist: measured drag se backwards kaam karo
Inverse problem
Ex 9
Hum ek recurring "house" flow use karenge taaki numbers familiar rahein: air , ν = 1.5 × 1 0 − 5 m 2 / s , ρ = 1.2 kg/m 3 , jab tak problem kuch aur na kahe.
Air (ν = 1.5 × 1 0 − 5 m 2 / s ) U ∞ = 10 m/s par flat plate ke upar. x = 0.5 m par δ nikalo.
Forecast: Layer thickness guess karo — millimetres, centimetres, ya metres? (Boundary layers famously thin hoti hain.)
Step 1 — Reynolds number. R e x = ν U ∞ x = 1.5 × 1 0 − 5 10 × 0.5 = 3.33 × 1 0 5 .
Yeh step kyun? R e x master dial hai; iske bina kuch bhi compute nahi ho sakta, aur yeh bhi batata hai ki Blasius legal bhi hai ya nahi (chahiye R e x ≲ 5 × 1 0 5 ). Yahan 3.33 × 1 0 5 < 5 × 1 0 5 ✓ laminar.
Step 2 — δ mein plug in karo. δ = R e x 5.0 x = 3.33 × 1 0 5 5.0 × 0.5 = 577 2.5 = 4.3 × 1 0 − 3 m .
Yeh step kyun? Boxed formula master dial ko physical length mein convert karta hai.
Verify: Aadhe metre plate par 4.3 mm — slow-zone, travel ki distance ka 116 mein se ek hissa. Thin, as forecast ✓. Units: dimensionless m = m ✓.
Same air, same U ∞ = 10 m/s . x = 0.5 m par wall shear stress τ w nikalo.
Forecast: Air halka hai. Kya drag pressure atmospheric (1 0 5 Pa) ka bada fraction hoga ya tiny?
Step 1 — local friction coefficient. C f = R e x 0.664 = 577 0.664 = 1.15 × 1 0 − 3 .
Yeh step kyun? C f dimensionless hai aur sirf master dial par depend karta hai — pehle ise compute karo, baad mein un-normalize karo. Yeh "physics ki shape" ko "numbers ki size" se alag karta hai.
Step 2 — dynamic pressure se un-normalize karo. τ w = C f ⋅ 2 1 ρ U ∞ 2 = 1.15 × 1 0 − 3 × 2 1 × 1.2 × 1 0 2 = 0.069 Pa .
Yeh step kyun? 2 1 ρ U ∞ 2 flow ka natural pressure scale hai; multiply karne se force-per-area ke real units wapas aate hain.
Verify: 0.069 Pa vs atmospheric ∼ 1 0 5 Pa — bilkul tiny, as forecast (air drag gentle hota hai) ✓. Units: dimensionless × kg/m 3 × ( m/s ) 2 = kg / ( m⋅s 2 ) = Pa ✓.
Plate length L = 0.5 m , width b = 0.2 m , same air, U ∞ = 10 m/s . Ek side par drag nikalo.
Forecast: Tiny stress, chhota area — force newtons mein hoga ya newton ka chhota fraction?
Step 1 — averaged drag coefficient. R e L = 1.5 × 1 0 − 5 10 × 0.5 = 3.33 × 1 0 5 , toh C D = R e L 1.328 = 577 1.328 = 2.30 × 1 0 − 3 .
Yeh step kyun? Local τ w leading edge par sabse zyada hota hai aur 1/ x ki tarah decay karta hai, isliye hum ek point ka stress area se multiply nahi kar sakte. C D plate par C f ka average hai. Factor 2 kahan se aaya? Average C f matlab L 1 ∫ 0 L C f d x , aur kyunki C f ∝ x − 1/2 , integral ∫ 0 L x − 1/2 d x = 2 L hai, toh average L 2 L × ( value at L ) / L − 1 hai — constants carry through karne par exactly C ˉ f = 2 C f ( L ) milta hai. Isliye C D = 2 C f ( L ) .
Step 2 — wetted area par un-normalize karo. F D = C D ⋅ 2 1 ρ U ∞ 2 ⋅ ( b L ) = 2.30 × 1 0 − 3 × 60 × 0.1 = 0.0138 N .
Yeh step kyun? Pehle jaisa hi pressure scale, ab plate area b L = 0.1 m 2 se multiply hua jo fluid actually rub karta hai.
Verify: 0.0138 N ≈ 1.4 grams-force — feather-light, tiny τ w se match karta hai ✓. Cross-check: C D = 2 C f = 2 × 1.15 × 1 0 − 3 = 2.30 × 1 0 − 3 ✓.
Same air aur speed. x → 0 par (leading edge par) δ aur τ w ka kya hota hai? x = 0.5 m , x = 0.05 m , x = 0.005 m par trend evaluate karo. (Yaad karo μ = ρ ν = 1.8 × 1 0 − 5 Pa⋅s , upar define ki gayi hai.)
Forecast: Plate ne abhi drag karna shuru kiya. Layer yahan sabse moti hai ya pateeli — aur drag sabse halka hai ya tezz?
Step 1 — dono ko pure powers of x mein likho. Boxed δ = 5.0 x / R e x se shuru karo aur R e x = U ∞ x / ν substitute karo:
δ = U ∞ x / ν 5.0 x = 5.0 x U ∞ x ν = 5.0 U ∞ x ν x 2 = 5.0 U ∞ ν x ∝ x 1/2 .
Isi tarah, upar derive kiye gaye μ -form se, τ w = 0.332 μ U ∞ U ∞ / ( ν x ) ∝ x − 1/2 .
Yeh step kyun? Akele x ko square-root ke andar push karne se degenerate behaviour expose hota hai: ek zero se grow karta hai, doosra blow up karta hai.
Step 2 — trend tabulate karo.
x (m)
δ = 5.0 ν x / U ∞ (mm)
τ w = 0.332 μ U ∞ U ∞ / ν x (Pa)
0.5
4.33
0.069
0.05
1.37
0.218
0.005
0.433
0.690
Yeh step kyun? x ko tenfold shrink karne se δ , 10 ≈ 3.16 se shrink hota hai aur τ w usi factor se grow karta hai — dono ek doosre ke mirror hain.
Step 3 — limit lo. Jab x → 0 + : δ → 0 aur τ w → + ∞ .
Yeh step kyun? Yeh leading-edge singularity hai: model x = 0 par infinite shear predict karta hai. Yeh physical nahi hai — thin-layer assumption δ ≪ x tab break hoti hai jab δ ∼ x — lekin yeh sahi tarike se warn karta hai ki drag leading edge par concentrate hota hai.
Verify: δ , x ki tarah scale karta hai (x mein ÷10 ⇒ δ mein ÷3.16: 4.33/3.16 = 1.37 ✓), τ w , 1/ x ki tarah (x mein ÷10 ⇒ τ w mein ×3.16: 0.069 × 3.16 = 0.218 ✓). Singularity boxed formulas ki genuine feature hai.
Neeche ki figure exactly is table ko do curves ke roop mein plot karti hai taaki aap mirror-image trends aur origin par blow-up dekh sako.
Worked example Figure — Cell D trends
Blue (left axis): δ ∝ x zero se rise karta hai — layer leading edge par kuch nahi hoti aur phir tezi se thickens. Orange (right axis): τ w ∝ 1/ x vertical axis ki taraf dive karta hai — shear leading edge par sabse tezz hoti hai aur downstream fade hoti hai. Teeno blue dots tabulated x -values hain. Dekho ki dono curves reflections hain: jahan ek shrink kare, doosra grow kare.
Ek plate U ∞ par test hoti hai. Phir speed chaar guna hoke 4 U ∞ ho jaati hai, same fluid, same point x . Kis factor se (a) δ change hoga, (b) τ w change hoga, (c) total drag F D change hoga?
Forecast: Faster flow — thinner ya thicker layer? Zyada ya kam drag, aur kitna?
Step 1 — exponents padho. δ ∝ U ∞ − 1/2 (from δ = 5.0 ν x / U ∞ ), toh δ , 4 − 1/2 = 2 1 se change hota hai.
Yeh step kyun? Hume sirf chahiye ki har quantity U ∞ par kaise depend karti hai; constants ratio mein cancel ho jaate hain. Faster free stream ⇒ inertia zyada jeette hai ⇒ thinner layer.
Step 2 — shear stress. τ w = 0.332 μ U ∞ U ∞ / ( ν x ) ∝ U ∞ 3/2 , toh yeh 4 3/2 = 8 se change hota hai.
Yeh step kyun? U ∞ ka ek power momentum flux se, steeper wall gradient se aadha power — milke U ∞ 3/2 .
Step 3 — total drag. Pehle note karo R e L = U ∞ L / ν ∝ U ∞ (definition se, L , ν fixed). Toh C D = 1.328/ R e L ∝ R e L − 1/2 ∝ U ∞ − 1/2 . Phir F D ∝ C D U ∞ 2 ∝ U ∞ − 1/2 U ∞ 2 = U ∞ 3/2 , toh F D , 4 3/2 = 8 se change hota hai.
Yeh step kyun? Hidden speed-dependence C D ke andar R e L ∝ U ∞ mein rehti hai; usse expose karne ke baad hi powers sahi add hoti hain. F D aur τ w same U ∞ 3/2 law share karte hain.
Verify: δ : × 0.5 ; τ w : × 8 ; F D : × 8 . Sanity: δ down matlab gradient steeper ⇒ zyada shear ⇒ consistent hai ki τ w bada ✓. Numerically 4 − 0.5 = 0.5 , 4 1.5 = 8 ✓.
Same air, U ∞ = 10 m/s . Kisi distance x tr par flow transition Reynolds number R e x , crit = 5 × 1 0 5 tak pahunch jaata hai. Uske aage, kya δ = 5.0 x / R e x valid hai?
Forecast: Guess karo: kitne centimetres, tens of centimetres, ya metres mein turbulence set in hogi?
Step 1 — Reynolds definition invert karo. R e x = U ∞ x / ν ⇒ x tr = U ∞ R e x , crit ν = 10 5 × 1 0 5 × 1.5 × 1 0 − 5 = 0.75 m .
Yeh step kyun? Validity boundary ek Reynolds number ke roop mein di gayi hai, toh ise distance mein convert karo jo plate par measure ki ja sake.
Step 2 — formula par verdict. x > 0.75 m ke liye flow turbulent hai; Blasius δ = 5.0 x / R e x ab apply nahi hota . Turbulent estimate δ / x ≈ 0.37 R e x − 1/5 use karo (dekho Turbulent Boundary Layer ).
Yeh step kyun? Blasius ek laminar solution hai — ise transition ke baad apply karne par wrong power law (x vs x 4/5 ) over-predict hoti hai.
Verify: x = 0.5 m par (Ex 1) hamaare paas R e x = 3.33 × 1 0 5 < 5 × 1 0 5 tha, safely laminar — consistent, kyunki 0.5 < 0.75 m ✓. Recompute: 0.75 × 10/1.5 × 1 0 − 5 = 5 × 1 0 5 ✓.
Same air, U ∞ = 10 m/s , x = 0.5 m par. (a) Kisi height y par velocity free stream ki aadhi hai, u = 5 m/s ? (b) Dikhao ki wahan vertical velocity v zero nahi hai. Upar define kiye Blasius read-offs use karo: f ′ = 0.5 at η ≈ 1.6 , aur wahan f ≈ 0.65 .
Forecast: Half-speed height — kya yeh 4.3 mm layer ke andar wall ke paas, mid-layer, ya top ke paas hai?
Step 1 — η ko y mein wapas convert karo. Yaad karo (definition block se) η = y U ∞ / ( ν x ) , toh y = η ν x / U ∞ . Scale ν x / U ∞ = 1.5 × 1 0 − 5 × 0.5/10 = 8.66 × 1 0 − 4 m . Phir y = 1.6 × 8.66 × 1 0 − 4 = 1.39 × 1 0 − 3 m .
Yeh step kyun? η "clever stretched height" hai jo profile ko universal banata hai; real millimetre reading paane ke liye ise un-stretch karna hoga. Kyunki f ′ ( η ) = u / U ∞ hai, f ′ = 0.5 exactly wahan mark karta hai jahan u = 5 m/s hai.
Step 2 — layer mein locate karo. δ = 4.33 mm (Ex 1), aur y = 1.39 mm , toh half-speed y / δ = 0.32 par baith ta hai — roughly ek tihai raaste upar.
Yeh step kyun? Confirm karta hai ki profile wall ke paas steep hai: sirf ek tihai height ke baad aadhi speed aa jaati hai, phir 0.99 U ∞ tak crawl karta hai.
Step 3 — vertical velocity nonzero hai. Parent ka formula hai v = 2 1 ν U ∞ / x ( η f ′ − f ) . Plug in karo: root hai 1.5 × 1 0 − 5 × 10/0.5 = 3 × 1 0 − 4 = 1.73 × 1 0 − 2 , aur ( η f ′ − f ) = ( 1.6 × 0.5 − 0.65 ) = ( 0.8 − 0.65 ) = 0.15 . Toh v = 0.5 × 1.73 × 1 0 − 2 × 0.15 = 1.30 × 1 0 − 3 m/s .
Yeh step kyun? Continuity fluid ko upar force karti hai jaise layer thickens — ek real, chahe tiny, v > 0 . Ise ignore karna (common mistake) mass conservation todta hai.
Verify: y = 1.39 mm , δ = 4.33 mm ke andar hai ✓ (hona hi chahiye). v = 1.30 × 1 0 − 3 m/s yani U ∞ ka ∼ 1 0 − 4 — small lekin strictly positive ✓. v ke units: m 2 / s 2 = m/s ✓.
Figure poora universal profile dikhata hai taaki aap dekh sako ki half-speed itna neeche kyun hai.
Worked example Figure — Cell G universal profile
Blue curve : Blasius profile f ′ ( η ) = u / U ∞ , real height y (mm) ke against plotted, x = 0.5 m par. Red arrow origin par: no-slip, wall par u = 0 . Orange dot : half-speed point, sirf y ≈ 1.4 mm par — roughly ek tihai raaste upar. Green dashed line : layer ka edge, δ (jahan u = 0.99 U ∞ ). Notice karo ki curve wall ke paas kitni steep hai aur free stream ke paas kaise flatten hoti hai — woh steepness hi wall shear hai.
Ek ship model panel water (ν = 1.0 × 1 0 − 6 m 2 / s , ρ = 1000 kg/m 3 ) mein U ∞ = 2 m/s par run karta hai. Panel: L = 0.3 m , width b = 0.5 m . (a) Trailing edge par δ nikalo, (b) ek side par drag nikalo. Laminar validity check karo.
Forecast: Water kinematically hamare air se 15× kam viscous hai aur bahut denser. Air se thinner ya thicker layer? Bada ya chhota drag?
Step 1 — Reynolds & validity. R e L = 1.0 × 1 0 − 6 2 × 0.3 = 6.0 × 1 0 5 .
Yeh step kyun? Pehle master dial; aur 6.0 × 1 0 5 > 5 × 1 0 5 hai, toh strictly yeh just past laminar transition hai — Blasius sirf estimate deta hai. Hum pehli approximation ke roop mein proceed karte hain aur flag karte hain.
Step 2 — thickness. δ = 6.0 × 1 0 5 5.0 × 0.3 = 775 1.5 = 1.94 × 1 0 − 3 m .
Yeh step kyun? Same boxed formula; low ν ⇒ high R e ⇒ thin layer (≈1.9 mm, length ke relative air ke 4.3 mm se thinner).
Step 3 — drag. C D = 6.0 × 1 0 5 1.328 = 1.71 × 1 0 − 3 . Phir F D = C D ⋅ 2 1 ρ U ∞ 2 ⋅ b L = 1.71 × 1 0 − 3 × 2 1 × 1000 × 4 × 0.15 = 0.514 N .
Yeh step kyun? Water ki density air se ~800× zyada hai, toh small C D ke bawajood force ab sizeable hai (aadha newton, milli-newtons nahi).
Verify: 2 1 ρ U ∞ 2 = 2000 Pa ; b L = 0.15 m 2 ; 2000 × 0.15 = 300 N × 1.71 × 1 0 − 3 = 0.514 N ✓. Dense fluid ⇒ bada drag, as forecast ✓. Caveat: R e L slightly above 5 × 1 0 5 hai, toh upper-ballpark laminar estimate samjho.
Ek wind tunnel mein, ek plate (L = 1.0 m , b = 0.4 m , air ν = 1.5 × 1 0 − 5 , ρ = 1.2 ) ek-sided drag F D = 0.02 N record karta hai. Flow speed U ∞ nikalo. (Laminar Blasius assume karo.)
Forecast: Inverse problem — force pata hai, cause chahiye. Kya U ∞ kuch m/s aayega ya tens of m/s?
Step 1 — F D ko U ∞ ke terms mein poora likho. F D = C D ⋅ 2 1 ρ U ∞ 2 ⋅ b L = U ∞ L / ν 1.328 ⋅ 2 1 ρ U ∞ 2 b L .
Yeh step kyun? Hum U ∞ ko isolate nahi kar sakte jab tak har hidden U ∞ (andar C D mein bhi, via R e L = U ∞ L / ν !) page par na ho — classic trap yeh hai ki C D speed par depend karta hai yeh bhool jana.
Step 2 — U ∞ ka power collect karo. C D ∝ U ∞ − 1/2 aur dynamic pressure ∝ U ∞ 2 , toh F D ∝ U ∞ 3/2 . Explicitly:
F D = 1.328 L ν ⋅ 2 1 ρ b L U ∞ 3/2 = 0.664 ρ b ν L U ∞ 3/2 .
Yeh step kyun? Single power U ∞ 3/2 tak reduce karne se inversion ek clean root-type solve ban jaati hai.
Step 3 — U ∞ solve karo. Prefactor K = 0.664 × 1.2 × 0.4 × 1.5 × 1 0 − 5 × 1.0 = 0.664 × 1.2 × 0.4 × 3.873 × 1 0 − 3 = 1.234 × 1 0 − 3 . Phir U ∞ 3/2 = F D / K = 0.02/1.234 × 1 0 − 3 = 16.21 , toh U ∞ = 16.2 1 2/3 = 6.42 m/s .
Yeh step kyun? 3/2 exponent undo karne ke liye dono sides ko 2/3 power tak raise karo — physics ka algebraic inverse.
Verify (forward plug karo): R e L = 6.42 × 1.0/1.5 × 1 0 − 5 = 4.28 × 1 0 5 < 5 × 1 0 5 laminar ✓. C D = 1.328/ 4.28 × 1 0 5 = 2.03 × 1 0 − 3 ; F D = 2.03 × 1 0 − 3 × 2 1 × 1.2 × 6.4 2 2 × 0.4 = 2.03 × 1 0 − 3 × 2 1 × 1.2 × 41.2 × 0.4 = 0.0201 N ≈ 0.02 N ✓. Round-trip close ho jaata hai.
Recall Quick self-test
Kaun sa cell leading edge par blow up karta hai, aur kaunsi quantity? ::: Cell D — τ w → ∞ jab x → 0 , jabki δ → 0 .
Agar U ∞ double karo, toh total drag kis factor se change hoga? ::: 2 3/2 ≈ 2.83 (drag ∝ U ∞ 3/2 ).
Kissi bhi Blasius answer par trust karne se pehle R e recompute kyun karna chahiye? ::: Confirm karne ke liye ki R e < 5 × 1 0 5 (laminar); uske upar x laws fail ho jaate hain — dekho Turbulent Boundary Layer .
Inverse problem (Cell I) mein kya trap hai? ::: Yeh bhool jaana ki C D khud U ∞ par depend karta hai (R e L ∝ U ∞ ke through), jisse F D ∝ U ∞ 3/2 milta hai, U ∞ 2 nahi.
μ aur ν ka relation kya hai? ::: μ = ρ ν (dynamic = density × kinematic).
Speed 4 guna karo ⇒ layer aadhi, shear ×8, drag ×8. Kyunki δ ∝ U − 1/2 , τ w ∝ U 3/2 , F D ∝ U 3/2 . Distance: δ ∝ x (grow karta hai), τ w ∝ 1/ x (fade hota hai).
Dekho bhi: Prandtl Boundary Layer Theory · Navier–Stokes Equations · Skin Friction Drag · Similarity Solutions in PDEs · Stream Function and Vorticity .