4.7.8 · Maths › Partial Differential Equations
Intuition Ek-sentence mein picture
Heat, hot se cold ki taraf flow karti hai, aur ek chhoti si slab kitni tezi se heat hoti hai yeh depend karta hai ki uska temperature profile kitna curved hai. Ek bump (concave-down) cool hota hai; ek dip (concave-up) warm hota hai. Yeh "curvature drives change in time" wali baat exactly u t = α u xx hai.
Ek patli insulated rod jo x -axis ke saath rakhi hai. Maano
u ( x , t ) = position x par, time t par temperature (units K),
ρ = density (kg/m³), c = specific heat capacity (J/kg·K), k = thermal conductivity (W/m·K).
Hum u ( x , t ) ko govern karne wali PDE chahte hain. Yeh parabolic hai kyunki A u xx + 2 B u x t + C u tt ka discriminant B 2 − A C = 0 deta hai (yahan C = 0 , B = 0 ).
x aur x + Δ x ke beech ka ek slab lo, cross-section area A ke saath.
Worked example Step-by-step energy balance
Step 1 — Slab mein stored energy.
E ( t ) = ∫ x x + Δ x ρ c u ( s , t ) A d s .
Kyun? Energy density ρ c u times volume element A d s , slab ke across sum kiya.
Step 2 — Stored energy ke change ki rate.
d t d E = ∫ x x + Δ x ρ c u t A d s .
Kyun? Integral ke andar differentiate karo (limits space mein fixed hain).
Step 3 — Net heat jo ANDAR aa rahi hai. Flux left face se enter karta hai, right face se leave karta hai:
net in = A q ( x , t ) − A q ( x + Δ x , t ) .
Kyun? q ( x ) flow + x direction mein hai; energy entering left = + A q ( x ) , leaving right = + A q ( x + Δ x ) , toh net gain = in − out.
Step 4 — Equate karo (conservation).
∫ x x + Δ x ρ c u t A d s = A ( q ( x , t ) − q ( x + Δ x , t ) ) .
Step 5 — A Δ x se divide karo, Δ x → 0 lo.
Left side mean value → ρ c u t . Right side:
Δ x q ( x ) − q ( x + Δ x ) Δ x → 0 − ∂ x ∂ q .
Kyun? Yeh q ke derivative ki definition ka minus hai. Toh:
ρ c u t = − ∂ x ∂ q .
Step 6 — Fourier's law q = − k u x daalo:
ρ c u t = − ∂ x ∂ ( − k u x ) = k u xx ( k constant ) .
Step 7 — Diffusivity define karo α = ρ c k (m²/s):
u t = α u xx
Intuition "Parabolic" / infinite speed kyun
Heat eq ko 0 ⋅ u tt + 0 + ( − α ) u xx + u t = 0 likhkar A u tt + 2 B u x t + C u xx classify karne par discriminant B 2 − A C = 0 milta hai ⇒ parabolic . Physically iska matlab hai smoothing/diffusion ke saath infinite signal speed (ek disturbance ek dam se har jagah feel hoti hai, wave equation ke unlike).
Worked example Kya yeh point heat hoga ya cool?
Ek rod ka ek instant par u ( x , t ) = 20 + 5 sin ( π x ) K hai (jahan α = 2 ). x = 1/2 par u t nikalo.
Step 1: u xx = ∂ x 2 [ 20 + 5 sin π x ] = − 5 π 2 sin π x .
Kyun? Do baar differentiate karo; constants vanish ho jaate hain.
Step 2: x = 1/2 par, sin ( π /2 ) = 1 , toh u xx = − 5 π 2 .
Step 3: u t = α u xx = 2 ( − 5 π 2 ) = − 10 π 2 ≈ − 98.7 K/s.
Interpretation: bump ki peak tezi se cool hoti hai — intuition ke consistent hai.
Worked example Ek known solution check karo
Dikhao ki u = e − α π 2 t sin ( π x ) equation u t = α u xx ko solve karta hai.
Step 1: u t = − α π 2 e − α π 2 t sin π x . Kyun? Exponential par chain rule.
Step 2: u x = π e − α π 2 t cos π x , phir u xx = − π 2 e − α π 2 t sin π x .
Step 3: α u xx = − α π 2 e − α π 2 t sin π x = u t . ✓ Yeh exponentially decay karta hai — heat khatam ho jaati hai.
Recall Check karne se pehle predict karo
Agar α double kar do, toh rod jaldi equilibrate hogi ya dheere? Forecast…
Ek linear initial profile u = a + b x : har jagah u t kya hai? Forecast…
Verify: (1) Jaldi — u t , α ke saath scale karta hai, bada diffusivity = faster smoothing. (2) u xx = 0 ⇒ u t = 0 : already steady hai, kuch nahi badlega.
Common mistake "Heat equation mein waves jaisa
u tt hona chahiye."
Kyun sahi lagta hai: wave equation u tt = c 2 u xx famous hai aur space/time mein symmetric hai. Fix: energy storage ρ c u use karta hai (t mein first order), aur Fourier's law flux mein first order hai ⇒ sirf ek time derivative. Yeh time-reversal symmetry tod deta hai (heat spread hoti hai, kabhi un-spread nahi hoti).
Common mistake Fourier's law mein minus sign bhool jana.
Kyun sahi lagta hai: "flux proportional to gradient" clean lagta hai. Fix: minus ke bina, heat hotter regions ki taraf flow karegi, jo thermodynamics violate karta hai. Minus final equation ko + α u xx banata hai (diffusive, stable). Galat sign se u t = − α u xx milega (backward heat equation — bilkul unstable).
α = k treat karna.
Kyun sahi lagta hai: k "hi" conduction constant hai. Fix: PDE ka time scale diffusivity α = k / ( ρ c ) se set hota hai. Ek material achha conduct kar sakta hai (k bada) lekin bahut energy store kar sakta hai (ρ c bada), jisse moderate α milta hai.
Recall Ek 12-saal ke bachche ko samjhao
Socho ek metal ruler, beech mein hot, ends par cold. Heat ek bheed ki tarah hai jo hamesha ek bheed waali (hot) jagah se khaali (cold) jagah ki taraf push karti hai. Ek tiny spot ka temperature kitni tezi se badlega yeh depend karta hai ki woh temperature map par ek "hill" hai ya "valley": hills apne neeche wale neighbours ko heat lose karti hain aur cool hoti hain; valleys heat gain karti hain aur warm hoti hain. Equation u t = α u xx bas yahi kehti hai ki "time mein temperature kitni tezi se badle = (ek constant) × temperature kitni curved hai." Badi curve → fast change; flat → koi change nahi.
"Time-change equals C urvature × D iffusivity" → u t = α u xx .
Aur Fourier: "flux flows down-hill" → q = − k u x (down = minus).
Fourier's law of heat conduction (1D) batao aur minus sign explain karo. q = − k u x ; minus isliye kyunki heat hot se cold ki taraf flow karti hai, yaani temperature gradient ke opposite.
Kaun se do physical principles heat equation derive karte hain? Fourier's law (flux ∝ −gradient) aur control slab par conservation of energy.
1D heat equation likho aur α define karo. u t = α u xx jahan thermal diffusivity α = k / ( ρ c ) , units m²/s.
Heat equation ko parabolic kyun kaha jaata hai? Iska second-order classification discriminant B 2 − A C = 0 deta hai; physically yeh diffusive/smoothing hai.
Derivation mein ρ c u kya represent karta hai? Thermal energy density (J/m³): density × specific heat × temperature.
Sirf ek time derivative kyun hai (u tt nahi)? Energy storage aur Fourier flux dono time mein first order hain, toh balance sirf ek u t deta hai.
Kisi point par u xx > 0 ka physical matlab? Profile concave up hai (ek dip); u t > 0 toh point warm hota hai.
1D mein steady-state temperature profile ki shape (source ke bina)? Linear, kyunki u xx = 0 .
Kya α double karne se equilibration speed up hoti hai ya slow down? Speed up hoti hai; u t , α ke saath scale karta hai.
Internal source f ke saath heat eq? ρ c u t = k u xx + f .
Stored energy E over slab
dE/dt = integral rho c u_t
Net heat in = A q x minus q x+dx
Heat equation u_t = alpha u_xx