2.2.18 · D3 · Physics › Fluid Mechanics › Navier-Stokes equations — derivation from Newton's second la
Yeh page parent Navier–Stokes note ka drill ground hai. Wahan humne equation ko F = m a se banaya tha. Yahan hum ise har us situation mein use karte hain jo exam mein aa sakti hai — steady aur unsteady, viscous aur inviscid, rest mein aur gravity-driven, plus ek word problem aur ek exam twist.
Shuru karne se pehle, ek reminder us master sentence ka jisme hum hamesha plug karte hain:
Har symbol, taaki koi line one se hi na bhatak jaaye:
Definition Is page ka alphabet
ρ (rho) — fluid ki density , har cubic metre mein kitna mass packed hai, units kg / m 3 .
u — velocity field : space aur time ke har point par yeh batata hai ki wahan fluid kis direction mein aur kitni tez move kar raha hai. u x is ka x -axis ke along component hai.
P — pressure , kisi bhi surface par inward push per unit area, units Pa = N / m 2 .
μ (mu) — dynamic viscosity , fluid kitna "sticky/thick" hai, units Pa ⋅ s (dekho Viscosity & Newton's law of viscosity ).
g — gravity as a vector (iska ek direction hai, seedha neeche), magnitude likha jaata hai g = 9.8 m/s 2 . Jab bhi hum koi axis choose karte hain hume g ko us par project karna padta hai: axis x ke along component hai g x = g ⋅ x ^ . Toh g (arrow) aur g (number) alag objects hain — hamesha check karo ki formula kaunsa maang raha hai.
∇ (nabla) — "har direction mein slope" operator; ∇ P woh vector hai jo fastest pressure rise ki taraf point karta hai. ∇ 2 ise do baar apply karta hai (curvature).
D t D — Material derivative , acceleration jo chalte hue fluid ke saath measure ki jaati hai: D t D u = ∂ t ∂ u + ( u ⋅ ∇ ) u . Yeh explicitly Ex 5 (pure convection) aur Ex 6 (pure local term) mein aata hai.
∂ aur d ke baare mein ek baat (ek baar padho, har jagah use karo)
Hum ∂ / ∂ y (partial derivative) tab likhte hain jab koi quantity kai variables par depend karti ho aur hum sirf ek vary kar rahe hon, baaki ko fixed rakh ke. Hum d / d y (ordinary derivative) tabhi likhte hain jab problem ek variable ki function mein collapse ho jaaye. Toh ek full field u x ( x , y , z , t ) mein hum ∂ use karte hain; jab ek baar "steady + only-y -dependence" ise u x ( y ) tak reduce kar de, wahi object ek ODE ban jaata hai aur hum honestly d / d y par switch kar lete hain. Ex 2, 3, 6 mein yeh switch dekho — yeh visible proof hai ki variables freeze kiye gaye the.
Har problem jo tum meet kar sakte ho woh in cells mein se ek hai (ya blend). Neeche ke examples par cell label laga hai.
Cell
Case class
Kya switch on / off hai
Example
C1
At rest (u = 0 )
sirf pressure + gravity
Ex 1
C2
Steady viscous, pressure-driven
∂ t = 0 , μ = 0 , gradient drives
Ex 2
C3
Steady viscous, gravity-driven
μ = 0 , ρ g drives (falling film)
Ex 3
C4
Inviscid limit (μ = 0 )
viscous term drop karo → Euler
Ex 4
C5
Convective term dominant
( u ⋅ ∇ ) u = 0 , narrowing pipe
Ex 5
C6
Unsteady (transient start-up)
∂ t = 0
Ex 6
C7
Sign / direction check
∂ x P ka sign kis direction mein push karta hai?
Ex 7
C8
Degenerate / limiting numbers
μ → 0 , h → 0 , G → 0
Ex 8
C9
Real-world word problem
Poiseuille flow rate mein numbers lagao
Ex 9
C10
Exam twist (dimensionless)
NS terms se Reynolds number banao
Ex 10
C11
All terms live at once
unsteady + convective + viscous + free surface
Ex 11
Worked example Paani ka ek tank, 3 m gehraai. Talne par gauge pressure nikalo. (
ρ = 1000 kg/m 3 , g = 9.8 .)
Forecast: padhne se pehle number guess karo. Zyada gehraai → upar paani ka zyada weight → zyada pressure. Roughly 1 0 4 Pa per metre. 3 m ke liye kya milega?
Step 1. u = 0 har jagah set karo. Yeh step kyun? "At rest" exactly u = 0 hai; yeh poora left side aur poora viscous term kill kar deta hai (dono ko moving field chahiye).
0 = − ∇ P + ρ g
Step 2. Vertical component lo. z ko upar point karne do, toh g = ( 0 , 0 , − g ) , yaani iska z -projection g z = − g hai. Yeh step kyun? g ek vector hai; ek scalar equation padhne ke liye hum ise z par project karte hain. Horizontally kuch bhi vary nahi karta, toh ∂ x P = ∂ y P = 0 aur P sirf z par depend karta hai — partial ek ordinary derivative ban jaata hai.
0 = − ∂ z ∂ P − ρ g ⇒ d z d P = − ρ g
Step 3. Surface se neeche integrate karo. Yeh step kyun? d P / d z constant hai, toh pressure depth h ke saath linearly badhta hai.
P = P surface + ρ g h
Verify: P gauge = ρ g h = 1000 × 9.8 × 3 = 29400 Pa . Units: m 3 kg ⋅ s 2 m ⋅ m = m ⋅ s 2 kg = Pa ✓. Yeh exactly Hydrostatic pressure hai.
y = 0 aur y = h = 2 mm par do fixed plates ke beech paani, pressure gradient G = − ∂ P / ∂ x = 500 Pa/m se push ho raha hai, viscosity μ = 1 0 − 3 Pa⋅s . Centreline velocity nikalo.
Forecast: profile ek parabola hai, dono walls par zero (no-slip), middle mein fastest. Max speed millimetres/s hogi ya metres/s?
Figure mein kya dekhna hai: do moti horizontal slate lines fixed plates hain. Lavender curve jo daayein bulge kar rahi hai woh velocity profile u x ( y ) hai — notice karo ki yeh dono plates par zero se pinned hai (no-slip) aur middle mein sabse moti. Har chhota coral arrow us height par local flow speed hai: walls par kuch nahi se badhte hain centre par sabse lambe tak, parabola trace karte hue. Butter arrow neeche pressure push G = − ∂ P / ∂ x hai jo sab kuch drive karta hai, aur mint dot y = h /2 par u m a x mark karta hai. Picture hi derivation hai: constant curvature (Step 1) jo no-slip walls se forced hai (Step 2) sirf ek symmetric parabola ho sakti hai.
Step 1. Steady (∂ t = 0 ), flow sirf x ke along, sirf y par depend karta hai. Kyunki field ab single-variable function u x ( y ) tak reduce ho gayi hai, hum ∂ 2 / ∂ y 2 se ordinary d 2 / d y 2 par switch karte hain. Yeh step kyun? har dropped term literally zero hai — koi time change nahi, koi cross-flow nahi, toh convection ( u ⋅ ∇ ) u = 0 . Saath hi, continuity ∂ u x / ∂ x = 0 guarantee karta hai ki u x x par depend nahi kar sakta — yahi woh hai jo hume sirf u x ( y ) likhne deta hai.
0 = G + μ d y 2 d 2 u x
Step 2. ODE ko no-slip u x ( 0 ) = u x ( h ) = 0 ke saath solve karo. Yeh step kyun? Second derivative constant hai ⇒ do baar integrate karo; walls constants ko pin karte hain. Yeh Poiseuille flow hai.
u x ( y ) = 2 μ G y ( h − y )
Step 3. Max y = h /2 par. Yeh step kyun? Symmetric parabola midpoint par peak karti hai.
u m a x = 8 μ G h 2
Verify: u m a x = 8 × 1 0 − 3 500 × ( 0.002 ) 2 = 8 × 1 0 − 3 500 × 4 × 1 0 − 6 = 0.25 m/s . Units: Pa⋅s ( Pa/m ) ⋅ m 2 = s m ✓.
Worked example Tel ki ek thin film (
μ = 0.1 Pa⋅s , ρ = 900 ) ek vertical wall se neeche flow karti hai, film thickness δ = 1 mm , free surface par air hai (zero shear). Surface velocity nikalo.
Forecast: yahan gravity driver hai, pressure gradient nahi. Free surface sabse fast hai (koi cheez ise drag nahi karti), wall ruki hui hai. Half-parabola.
Step 1 — axes choose karo aur g project karo. x ko wall ke saath seedha neeche point karne do (flow direction) aur y wall se horizontally bahar measure karo. Kyunki x neeche ki vertical ke saath aligned choose kiya gaya hai, gravity vector fully x par project hota hai: g x = g ⋅ x ^ = + g aur g y = 0 . Yeh step kyun? g ek vector hai; x -momentum equation sirf iska x -component chahti hai. x ko g ke saath align karke hume clean g x = g milta hai, aur y -equation sirf thin film ke across (constant) hydrostatic pressure deta hai, jise hum ignore karte hain. Steady, sirf u x ( y ) (continuity phir se x -dependence forbid karta hai, toh ∂ 2 → d 2 ):
0 = ρ g x + μ d y 2 d 2 u x = ρ g + μ d y 2 d 2 u x ⇒ d y 2 d 2 u x = − μ ρ g
Step 2 — free-surface boundary condition. Wall par no-slip u x ( 0 ) = 0 ; free surface par zero shear d y d u x y = δ = 0 . Yeh step kyun? Air oil ko grip nahi kar sakti, toh viscous stress μ d u / d y open top par vanish hona chahiye — yeh free-surface pressure/stress condition hai: ek bahut halke fluid ke saath interface par, tangential stress → 0 aur surface pressure = atmospheric.
Step 3. Integrate karo: u x ( y ) = μ ρ g ( δ y − 2 y 2 ) . Yeh step kyun? Pehla integration + zero-shear slope constant fix karta hai; doosra + no-slip offset fix karta hai. y = δ par surface speed:
u surf = 2 μ ρ g δ 2
Verify: u surf = 2 × 0.1 900 × 9.8 × ( 0.001 ) 2 = 0.2 900 × 9.8 × 1 0 − 6 = 0.0441 m/s ✓.
Worked example Ek frictionless ideal fluid (
μ = 0 ) steadily flow karta hai. Dikhao ki NS Bernoulli-friendly Euler ban jaata hai, aur pressure drop evaluate karo jab speed 2 se 4 m/s tak double ho jaati hai (ρ = 1000 , horizontal, streamline ke along koi gravity term nahi).
Forecast: koi stickiness nahi, toh sirf pressure aur gravity bachte hain. Tez flow → lower pressure. Kitna lower?
Step 1. μ = 0 set karo. Yeh step kyun? Poora viscous term mar jaata hai, Euler equation (inviscid flow) bachta hai:
ρ D t D u = − ∇ P + ρ g
Step 2. Steady, horizontal streamline ke along, integrate karo → Bernoulli: P + 2 1 ρ u 2 = const . Yeh step kyun? ρ ( u ⋅ ∇ ) u kinetic term mein integrate hota hai.
Step 3. Δ P = − 2 1 ρ ( u 2 2 − u 1 2 ) . Yeh step kyun? Total constant hai: extra speed ke liye pay karne ke liye pressure girna padega.
Verify: Δ P = − 2 1 × 1000 × ( 4 2 − 2 2 ) = − 2 1 × 1000 × 12 = − 6000 Pa . Pressure 6 kPa drop karta hai ✓.
Worked example Ek narrowing channel mein steady flow:
x ke along velocity hai u x ( x ) = u 0 ( 1 + α x ) jahan u 0 = 1 m/s , α = 2 m − 1 . Material derivative use karke x = 0.5 m par fluid particle ka acceleration nikalo.
Forecast: bhaale ki ∂ u / ∂ t = 0 ho (steady!), ek particle tez paani mein move karke speed up karta hai. Toh acceleration zero nahi hai. Positive hai ya negative?
Step 1. Material derivative likho aur local term drop karo. Yeh step kyun? Material derivative true particle acceleration hai; steady flow ∂ u / ∂ t ko zero karta hai, sirf convection bachta hai:
D t D u x = = 0 ∂ t ∂ u x + ( u ⋅ ∇ ) u x = u x d x d u x
Step 2. d x d u x = u 0 α , toh D t D u x = u x ( x ) ⋅ u 0 α = u 0 ( 1 + α x ) ⋅ u 0 α = u 0 2 α ( 1 + α x ) . Yeh step kyun? Diye gaye profile par chain rule — yeh convective acceleration concrete form mein hai.
Verify: D t D u x = 1 2 × 2 × ( 1 + 2 × 0.5 ) = 2 × 2 = 4 m/s 2 . Positive → narrows mein accelerate ho raha hai, leaf-in-a-river picture se match karta hai ✓.
Worked example Stokes' first problem: ek fluid plate ke upar rest mein hai.
t = 0 par plate suddenly move karti hai. Plate ke paas equation hai ∂ t ∂ u = ν ∂ y 2 ∂ 2 u jahan ν = μ / ρ = 1 0 − 6 m 2 / s . t = 4 s baad motion kitni gehraai tak "diffuse" ho gayi hai?
Forecast: viscosity momentum ko heat ki tarah spread karti hai. Reach ki gayi depth ν t ki tarah badhti hai — linearly nahi. Penetration depth δ ∼ ν t estimate karo.
Step 1. Sirf ∂ t aur viscous term rakho. Yeh step kyun? Plate ke along koi pressure gradient nahi, koi gravity nahi; yahan flow dono y aur t par depend karta hai, toh hum genuinely partial derivatives rakhte hain (unlike Ex 2/3 jo ek variable tak collapse ho gaye the). Material derivative apne local piece ∂ u / ∂ t tak reduce ho jaata hai kyunki u ka x -dependence nahi hai. Yeh diffusion equation hai.
Step 2. Dimensional balance: t u ∼ ν δ 2 u ⇒ δ ∼ ν t . Yeh step kyun? Do surviving terms ki sizes match karne se growth law reveal hota hai.
Verify: δ = ν t = 1 0 − 6 × 4 = 4 × 1 0 − 6 = 2 × 1 0 − 3 m = 2 mm ✓. Slow diffusion — isliye viscous startup sluggish hota hai.
+ x ke along badhta hai: ∂ P / ∂ x = + 300 Pa/m . Pressure force per volume kis direction mein hai, aur ek fluid parcel ke saath kya hoga (dusre forces ignore karo)?
Forecast: sign trap. Common instinct: "high pressure aage push karta hai." Careful — force hai − ∇ P .
Step 1. Pressure force per volume f x = − ∂ P / ∂ x . Yeh step kyun? Humne parent mein derive kiya tha ki sirf difference matter karta hai, aur minus sign kehta hai fluid high se low pressure ki taraf push hota hai.
Step 2. f x = − ( + 300 ) = − 300 Pa/m = − 300 N/m 3 . Yeh step kyun? Diya gaya gradient plug karo.
Verify: negative → force − x direction mein point karta hai, yaani backward , + x mein move kar rahe parcel ko decelerate karta hai. Intuition se match karta hai: fluid rising pressure ke against climb karta hai ✓. Magnitude ∣ f x ∣ = 300 N/m 3 .
Worked example Ex 2 se Poiseuille result lo
u m a x = 8 μ G h 2 . Teen limits examine karo: (a) μ → 0 , (b) h → 0 , (c) G → 0 . Har ek physically kya matlab rakhta hai, aur u m a x evaluate karo agar h half karke 1 mm kar diya jaaye (Ex 2 ke G , μ rakhte hue)?
Forecast: thinner gap → dramatically slower flow kyunki h squared appear karta hai. h halving speed 4 factor se cut karegi, 2 se nahi.
Step 1 (a). μ → 0 : u m a x → ∞ . Yeh step kyun? Koi friction nahi matlab push resist karne ke liye kuch nahi — unphysical, humein bata raha hai real fluids hamesha kuch μ rakhte hain.
Step 2 (b). h → 0 : u m a x → 0 . Yeh step kyun? Gap band kar do aur dono walls fluid ko har jagah rok dete hain — flow vanish ho jaata hai.
Step 3 (c). G → 0 : u m a x → 0 . Yeh step kyun? Koi driving gradient nahi, koi flow nahi — Cell C1 (at rest) ke saath consistent.
Verify: h half karke 0.001 m karo: u m a x = 8 × 1 0 − 3 500 × ( 0.001 ) 2 = 8 × 1 0 − 3 500 × 1 0 − 6 = 0.0625 m/s . Yeh 0.25/4 hai — exactly Ex 2 ke answer ka quarter ✓.
μ = 3.5 × 1 0 − 3 Pa⋅s maano) ek capillary-jaisi round tube se flow karta hai, radius R = 1 mm , length L = 0.1 m , pressure difference Δ P = 200 Pa . Volume flow rate Q nikalo.
Forecast: round-pipe Poiseuille law hai Q = 8 μL π R 4 Δ P . Notice karo R 4 — flow radius ke liye extremely sensitive hai. m 3 / s mein ek tiny number expect karo.
Step 1. Cylindrical coordinates mein same NS reduction jaisa Ex 2 mein hai u ( r ) = 4 μL Δ P ( R 2 − r 2 ) deta hai. Yeh step kyun? Round geometry flat-plate parabola ko ek paraboloid se replace karti hai, phir bhi μ ∇ 2 u = − ∂ x P se.
Step 2. Speed ko cross-section par integrate karo: Q = ∫ 0 R u ( r ) 2 π r d r = 8 μL π R 4 Δ P . Yeh step kyun? Flow rate speed ko area par sum karna hai.
Verify: Q = 8 × 3.5 × 1 0 − 3 × 0.1 π ( 1 0 − 3 ) 4 × 200 = 2.8 × 1 0 − 3 π × 1 0 − 12 × 200 ≈ 2.244 × 1 0 − 7 m 3 / s ✓ (lagbhag 0.22 mL/s ). Dekho Poiseuille flow .
Worked example NS mein convective term aur viscous term ki sizes compare karke, woh dimensionless group derive karo jo decide karta hai "inertia vs friction." Water in a pipe ke liye evaluate karo:
U = 1 m/s , L = 0.05 m , ν = 1 0 − 6 m 2 / s .
Forecast: jo ratio milega woh Reynolds number hai R e = U L / ν . Bada R e = inertia wins (turbulent-prone); chhota R e = viscosity wins (smooth).
Step 1 — convective (inertia) term ki size nikalo. Har derivative ko (typical change)/(typical length) se replace karo: scale U ki velocity jo distance L par change ho rahi hai. Yeh step kyun? Hume sirf order of magnitude chahiye, toh exact derivatives unnecessary hain.
ρ ( u ⋅ ∇ ) u ∼ ρ U ⋅ L U = L ρ U 2
Step 2 — viscous term ki size nikalo. Viscous term second derivative hai, toh isme do length divisions hain: μ ∇ 2 u ∼ μ U / L 2 . Yeh step kyun? Same order-of-magnitude accounting, ∇ 2 par apply kiya (curvature ⇒ L se do baar divide karo).
μ ∇ 2 u ∼ μ L 2 U
Step 3 — ratio banao. Inertia ko viscosity se divide karo; common piece U / L cancel ho jaata hai aur μ / ρ = ν (kinematic viscosity) appear karta hai. Yeh step kyun? Do forces ka pure ratio dimensionless hona chahiye — woh number Reynolds number hai.
R e = μU / L 2 ρ U 2 / L = μ ρ U L = ν U L
Verify: R e = 1 0 − 6 1 × 0.05 = 50000 . Bada → inertia dominate karta hai, flow likely turbulent ✓. Dimensionless check: m 2 / s ( m/s ) ( m ) = pure number ✓.
Worked example Impulsively started pipe:
t = 0 par ek lamba round pipe (radius R ) pressure gradient G par switch hota hai. Parabola settle hone se pehle, saare NS terms matter kar sakte hain. Confirm karo kaunse survive karte hain, aur steady state ke baad mean speed nikalo agar R = 1 mm , G = 500 Pa/m , μ = 1 0 − 3 Pa⋅s .
Forecast: start-up par flow unsteady hai (∂ t = 0 ); agar pipe narrow bhi ho toh parcel accelerate karta hai (( u ⋅ ∇ ) u = 0 ); friction hamesha present hai (μ ∇ 2 u ); aur pipe outlet air se marta hai (free-surface / atmospheric pressure condition). Final settled mean speed centreline ki half hai.
Step 1 — kaunse terms alive hain? Transient ke dauran, x -momentum equation har term rakhti hai:
ρ ( unsteady ∂ t ∂ u x + convective ( u ⋅ ∇ ) u x ) = G + viscous μ ∇ 2 u x
Yeh step kyun? Yeh ek hi example hai jahan teen left-block effects mein se koi bhi switch off nahi hai — yeh genuine "combined" cell hai jo matrix mein promise ki gayi thi. Straight uniform pipe mein continuity ∂ u x / ∂ x = 0 convection kill karta hai; narrowing pipe mein yeh survive karta.
Step 2 — steady state tak settle karo. Jaise t → ∞ , ∂ t → 0 ; straight pipe mein convection off hai, toh hum steady Poiseuille u ( r ) = 4 μ G ( R 2 − r 2 ) recover karte hain centreline ke saath u m a x = 4 μ G R 2 . Yeh step kyun? Transient die away karta hai aur balance Cell C2/C9 tak reduce ho jaata hai — cells ke across ek healthy consistency check.
Step 3 — mean speed. Round pipe ke liye ek paraboloid ka area-average exactly peak ka half hota hai: u ˉ = 2 1 u m a x = 8 μ G R 2 . Yeh step kyun? Flow rate ÷ area; paraboloid ka mean uska max ka half hai.
Verify: u ˉ = 8 × 1 0 − 3 500 × ( 0.001 ) 2 = 8 × 1 0 − 3 500 × 1 0 − 6 = 0.0625 m/s ✓. (Atmospheric pressure par free-surface exit ke saath consistent jo downstream boundary set karta hai; continuity ne yahan convection ko zero fix kiya.)
Recall Tumhare liye kaunsa cell sabse mushkil tha?
Cell C7 sign trap (force hai − ∇ P ) ::: fluid high se low pressure ki taraf push hota hai, toh positive ∂ x P ek negative x -force deta hai.
Ex 5 mein steady flow ke bawajood acceleration kyun hai? ::: convective term ( u ⋅ ∇ ) u material derivative mein — parcel tez fluid mein move karta hai.
Ex 8 mein h half karne par u m a x kitne factor se change hota hai? ::: ek quarter, kyunki u m a x ∝ h 2 .
Pipe mein convective term kab survive karta hai (Ex 11)? ::: sirf tab jab pipe narrow ho; continuity ise straight uniform pipe mein kill kar deta hai.
Mnemonic Kisi bhi problem ko ek cell mein turn karo
Teen sawaal poocho: Moving? (C1 vs rest) → Time-changing? (C6/C11) → Friction on? (C4 vs viscous C2/C3). Surviving terms tumhara cell batate hain.