Fluidodinamica dei moti convettivi
1. Teorema del trasporto di Raynolds
Considero mezzo continuo.
Ipotesi di localizzazione.
Equilibrio termodinamico.
Sistema chiuso con m fissata - passo definire V, E, S.
Considero la generica proprietà Ψ e proprietà specifica ψ con Ψ = ∫R ψ dm = ∫R ρψ dV.
Allora: Ψ(t) = ∫Rt ρ(x,y,z,t) ψ(x,y,z,t) dV e Ψ(t+Δt) = ∫Rt+Δt ρ(x,y,z,t+Δt) ψ(x,y,z,t+Δt) dV.
So che generalmente d/dt Ψ = limΔt→0 Ψ(t+Δt) - Ψ(t) / Δt.
Rt = RA ∪ RB, Rt+Δt = RB ∪ RC = Rt - RA + RC.
Allora Ψ(t+Δt) - Ψ(t) = ∫RC ρt+Δtψt+ΔtdV - ∫RA ρt+Δtψt+ΔtdV- ∫Rt ρtψtdV.
Da cui dΨ/dt = ∫Rt ρt ∂ψ/∂t dV + ∫Rt ρψ u⋅m̂ dV.
Allora: limΔt→0 -1/Δt ∫RA ρt+Δtψt+Δt dV quindi: = ∫RC ρt+Δtψt+Δt dV.
Ricordo RC = RA ∪ RB.
Da cui: dΨ/dt = ∂(ρψ) / ∂t dV + ∫p ρψ û⋅m̂ dV = ∫∂/∂t (∇ (ρψ u) dV.
Fluidodinamica dei moti convettivi
1. Teorema del trasporto di Raynolds
Considero mezzo continuo.
Ipotesi di localizzazione.
Equilibrio termodinamico.
Sistema chiuso con m fissata.
Allora: Ψ(t) = ∫Rt ρ(x,y,z,t) ψ(x,y,z,t) dV e Ψ(t+Δt) = ∫Rt+Δt ρ(x,y,z,t)t+Δt ψ(x,y,z,t)t+Δt dV.
So che generalmente Ψ(t+Δt) - Ψ(t) = ∫Rc ρt+Δt ψt+Δt.
Da cui dΨ/dt = ∫Rt ∂ρψ/∂t dV + ∮∂R ρψ ũṙm dS = ∫Rt ∂/∂t (ρψ) dV = ∫Rt ∂ψ/∂t dV.
2. Bilancio locale della massa
Ψ = ρm, m = mt + mf allora dΨ⁄dt = 0.
dρ⁄dt + ∇(ρū) = 0.
Per moto incomprimibile vale ρ = const allora ∇ū = 0.
In generale dρ⁄dt + ∇(ρū)dv = ∫R ∂ρ⁄∂t + ρ(ū ∙ ∇ ψ̅ + ψ̅ ∇∙ū) dv.
d⁄dt ∫R ρΨ dv = ∫R ρ ∂Φ⁄∂t dt.
3. Bilancio della qt. di moto
Ψ = ρq̅ = m∙ū Vettoriale! Scrivo il bilancio per la i-esima componente.
ūi = qi⁄m = ṁi⁄m = υi.
So F = m∙ū = m d q̅⁄dt allora dv1⁄dt = f̅.
Teo di Reynolds ∫R ∂ ⁄∂t dV = Fi.
Distanza (gravità) Fc(b) = ∫p ρ gidV.
Contatto (pressione & viscosità) Fi(s) = ∫ tσij = -pδij + τij Tensore degli sforzi viscosi.
∇2ū + 1⁄3 ∇ ∇∙ū.
Bilancio di energia
Ψ = ΕΨ = ε + u̅i̅ u̅i̅ / 2.
dE / dt = Q̇ + Ẇ.
Teo del trasporto ∫R ρ (dΕ / dt + u̅i̅ Du̅i̅ / Dt) dV = Q̇ + Ẇ.
Q̇ = ∫R q̅g̅ dV - ∫∂ q̅m̅ dΓ̅.
Gauss-Green = ∫R ∇q̅g̅ dV ∫R (q̅g̅ - ∇q̅g̅) dV.
Ẇ = ∫R ρg̅ u̅i̅ dV + ∫∂ t̅ij̅ m̅ dΓ̅ = ∫R ρg̅i̅ u̅i̅ dV + ∫R ∇(σ̅ij u̅i̅) dV.
= ∫R ρg̅i̅ u̅i̅ dV + = ∫R ρg̅i̅ u̅i̅ dV + ∫R u̅i̅ ∂σij̅ / ∂xj dV + ∫ ∂ui̅ σij̅ / ∂xj dV.
σij̅ Dij̅ = Di̅j̅ p + Δια̅ij̅ = - p∇ · u̅ + Di̅:
Da cui ∫R ρ (dΕ / Dt + µ̇ Db̅i̅ / Dt) dV = ∫Rt q̅q - ∇² · q̅ + ρ dΕ / Dt = - ∇ · q̅ + q̅q - ρ∇ · u̅ + Di̅: Lavoro di compressione (se moto inc. = 0).
5) Approssimazione di Oberbeck - Boussinesq
dvi/dt + ∂(vgvi)/∂xj = 0.
ρ0(dvi/dt + uj ∂vi/∂xj) = -∂p/∂xi + ρgi + ∂τ