Forhold mellem Rankine og Hugoniot
Den Forholdet Rankine-Hugoniot udtrykker diskontinuitet af varierende mængder gennem en chokbølge eller en brudlinje i en gas. Det blev navngivet til ære for Pierre-Henri Hugoniot
og William Rankine .
Den generelle sag
Vi er interesserede i delvise differentialligninger i dimension 1 af typen:
∂w∂t+∂f(w)∂x=0{\ displaystyle {\ frac {\ partial w} {\ partial t}} + {\ frac {\ partial f (w)} {\ partial x}} = 0}
Disse ligninger siges at være hyperbolske, når den Jacobianske matrix af f er diagonaliserbar og har positive reelle værdier, hvilket antages at være verificeret her. Sådanne ligninger tillader regelmæssige løsninger og diskontinuerlige løsninger, som vi er interesserede i.
Vi integrerer konserveringsligningen ovenfor i nærheden af abscissen x c af diskontinuiteten:
∫xvs.-εxvs.+ε∂w∂tdx=f[w(xvs.-ε)]-f[w(xvs.+ε)]{\ displaystyle \ int _ {x_ {c} - \ varepsilon} ^ {x_ {c} + \ varepsilon} {\ frac {\ partial w} {\ partial t}} \, \ mathrm {d} x = f [ w (x_ {c} - \ varepsilon)] - f [w (x_ {c} + \ varepsilon)]}![{\ displaystyle \ int _ {x_ {c} - \ varepsilon} ^ {x_ {c} + \ varepsilon} {\ frac {\ partial w} {\ partial t}} \, \ mathrm {d} x = f [ w (x_ {c} - \ varepsilon)] - f [w (x_ {c} + \ varepsilon)]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bc553ac539a0b39439a628e59298454ac7adf6fd)
Ved hjælp af Leibniz integrationsregel kommer det:
vvs.[w(xvs.-ε)-w(xvs.+ε)]+∂∂t∫xvs.-εxvs.wdx+∂∂t∫xvs.xvs.+εwdx=f[w(xvs.-ε)]-f[w(xvs.+ε)]{\ displaystyle v_ {c} \, [w (x_ {c} - \ varepsilon) -w (x_ {c} + \ varepsilon)] + {\ frac {\ partial} {\ partial t}} \ int _ { x_ {c} - \ varepsilon} ^ {x_ {c}} w \ mathrm {d} x + {\ frac {\ partial} {\ partial t}} \ int _ {x_ {c}} ^ {x_ {c } + \ varepsilon} w \ mathrm {d} x = f [w (x_ {c} - \ varepsilon)] - f [w (x_ {c} + \ varepsilon)]}![{\ displaystyle v_ {c} \, [w (x_ {c} - \ varepsilon) -w (x_ {c} + \ varepsilon)] + {\ frac {\ partial} {\ partial t}} \ int _ { x_ {c} - \ varepsilon} ^ {x_ {c}} w \ mathrm {d} x + {\ frac {\ partial} {\ partial t}} \ int _ {x_ {c}} ^ {x_ {c } + \ varepsilon} w \ mathrm {d} x = f [w (x_ {c} - \ varepsilon)] - f [w (x_ {c} + \ varepsilon)]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3b8581a5df76c5628b9470b5510ada91a8c4fd4e)
Man bemærkede udbredelseshastigheden af diskontinuiteten.
vvs.=dxvs.dt{\ displaystyle v_ {c} = {\ frac {\ mathrm {d} x_ {c}} {\ mathrm {d} t}}}
Ved at gøre dette opnår vi springforholdet mellem w :
ε→0{\ displaystyle \ varepsilon \ rightarrow 0}
vvs.[w(xvs.-ε)-w(xvs.+ε)]=f[w(xvs.-ε)]-f[w(xvs.+ε)]{\ displaystyle v_ {c} \, [w (x_ {c} - \ varepsilon) -w (x_ {c} + \ varepsilon)] = f [w (x_ {c} - \ varepsilon)] - f [w (x_ {c} + \ varepsilon)]}![{\ displaystyle v_ {c} \, [w (x_ {c} - \ varepsilon) -w (x_ {c} + \ varepsilon)] = f [w (x_ {c} - \ varepsilon)] - f [w (x_ {c} + \ varepsilon)]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/083010e89ff28431decd10c50acf6aa0ef7b814c)
For at forenkle notationerne skriver vi
vvs.(wg-wd)=f(wg)-f(wd){\ displaystyle v_ {c} \, (w_ {g} -w_ {d}) = f (w_ {g}) - f (w_ {d})}
Burgers ligning
Et simpelt eksempel er Burgers ligning, der matcher definitionen ovenfor med og .
w=u{\ displaystyle w = u}
f(w)=u22{\ displaystyle f (w) = {\ frac {u ^ {2}} {2}}}
I dette tilfælde skrives springligningen:
vvs.(ud-ug)=ud2-ug22{\ displaystyle v_ {c} \, (u_ {d} -u_ {g}) = {\ frac {u_ {d} ^ {2} -u_ {g} ^ {2}} {2}}}
er
vvs.=ud+ug2{\ displaystyle v_ {c} = {\ frac {u_ {d} + u_ {g}} {2}}}
Et stationært chok indebærer derfor nødvendigvis det .
ud=-ug{\ displaystyle u_ {d} = - u_ {g}}
Eulers ligninger
Ustabilt problem
Vi anvender springforholdet for hver af Eulers ligninger :
Kontinuitet |
|
w=ρ{\ displaystyle w = \ rho} |
|
f(w)=ρV{\ displaystyle f (w) = \ rho V} |
|
vvs.(ρd-ρg)=ρdVd-ρgVg{\ displaystyle v_ {c} \, (\ rho _ {d} - \ rho _ {g}) = \ rho _ {d} V_ {d} - \ rho _ {g} V_ {g}}
|
Mængde af bevægelse |
|
w=ρV{\ displaystyle w = \ rho V} |
|
f(w)=ρV2+s{\ displaystyle f (w) = \ rho V ^ {2} + p} |
|
vvs.(ρdVd-ρgVg)=(ρdVd2+sd)-(ρgVg2+sg){\ displaystyle v_ {c} \, (\ rho _ {d} V_ {d} - \ rho _ {g} V_ {g}) = (\ rho _ {d} V_ {d} ^ {2} + p_ {d}) - (\ rho _ {g} V_ {g} ^ {2} + p_ {g})}
|
Energi |
|
w=ρE{\ displaystyle w = \ rho E} |
|
f(w)=(ρE+s)V{\ displaystyle f (w) = (\ rho E + p) V} |
|
vvs.(ρdEd-ρgEg)=(ρdEd+sd)Vd-(ρgEg+sg)Vg{\ displaystyle v_ {c} \, (\ rho _ {d} E_ {d} - \ rho _ {g} E_ {g}) = (\ rho _ {d} E_ {d} + p_ {d}) V_ {d} - (\ rho _ {g} E_ {g} + p_ {g}) V_ {g}}
|
Vi bemærkede:
-
ρ{\ displaystyle \ rho}
tætheden,
-
V{\ displaystyle V}
fart,
-
s{\ displaystyle p}
tryk,
-
E=e+v22{\ displaystyle E = e + {\ frac {v ^ {2}} {2}}}
samlet energi pr. masseenhed,
-
e{\ displaystyle e}
den interne energi pr. masseenhed.
Diskontinuiteten kan være af to typer:
- chokket, hvor alle mængder er diskontinuerlige,
- kontaktdiskontinuiteten, hvor hastigheden og trykket er kontinuerligt. Dette svarer til to strømlinjer, der glider ved siden af hinanden uden at trænge ind i hinanden og have det samme tryk (bevarelse af momentum).(vvs.=Vd=Vg){\ displaystyle \ left (v_ {c} = V_ {d} = V_ {g} \ right)}

Stationært højre chok
I tilfælde af et stationært chok som ved aerodynamik bliver springforholdene:
vvs.=0{\ displaystyle v_ {c} = 0}
ρdVd=ρgVgρdVd2+sd=ρgVg2+sg(ρdEd+sd)Vd=(ρgEg+sg)Vg{\ displaystyle {\ begin {array} {rcl} \ rho _ {d} V_ {d} & = & \ rho _ {g} V_ {g} \\ [0.6em] \ rho _ {d} V_ {d } ^ {2} + p_ {d} & = & \ rho _ {g} V_ {g} ^ {2} + p_ {g} \\ [0.6em] (\ rho _ {d} E_ {d} + p_ {d}) V_ {d} & = & (\ rho _ {g} E_ {g} + p_ {g}) V_ {g} \ end {array}}}![{\ displaystyle {\ begin {array} {rcl} \ rho _ {d} V_ {d} & = & \ rho _ {g} V_ {g} \\ [0.6em] \ rho _ {d} V_ {d } ^ {2} + p_ {d} & = & \ rho _ {g} V_ {g} ^ {2} + p_ {g} \\ [0.6em] (\ rho _ {d} E_ {d} + p_ {d}) V_ {d} & = & (\ rho _ {g} E_ {g} + p_ {g}) V_ {g} \ end {array}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cca07d839977b85559a6a254c7ba2f498dab4b6b)
Ved at forenkle energihopligningen ved springforholdet på massestrømningshastigheden opnår vi bevarelsen af den totale entalpi:
Hg=Hd{\ displaystyle H_ {g} = H_ {d}}
eller
H=E+sρ=e+V22+sρ{\ displaystyle H = E + {\ frac {p} {\ rho}} = e + {\ frac {V ^ {2}} {2}} + {\ frac {p} {\ rho}}}
I tilfælde af en ideel gas udtrækker vi forholdet, der forbinder "springene", det vil sige forholdet mellem værdierne nedstrøms og opstrøms for chokket ved at indføre Mach-nummeret opstrøms, antages at være en data om problemet
s=(γ-1)ρe{\ displaystyle p = (\ gamma -1) \ rho e}
Mg=VgγrTg{\ displaystyle M_ {g} = {\ frac {V_ {g}} {\ sqrt {\ gamma rT_ {g}}}}}
hvor r er den specifikke gaskonstant .
Disse forhold er som følger
R=ρdρg=VgVd=(γ+1)Mg2(γ-1)Mg2+2{\ displaystyle R = {\ frac {\ rho _ {d}} {\ rho _ {g}}} = {\ frac {V_ {g}} {V_ {d}}} = {\ frac {(\ gamma +1) M_ {g} ^ {2}} {(\ gamma -1) M_ {g} ^ {2} +2}}}
sdsg=1+2γγ+1(Mg2-1){\ displaystyle {\ frac {p_ {d}} {p_ {g}}} = 1 + {\ frac {2 \ gamma} {\ gamma +1}} (M_ {g} ^ {2} -1)}
TdTg=sdsgρgρd{\ displaystyle {\ frac {T_ {d}} {T_ {g}}} = {\ frac {p_ {d}} {p_ {g}}} {\ frac {\ rho _ {g}} {\ rho _ {d}}}}
MdMg=[R2+γ-12Mg2(R2-1)]-12{\ displaystyle {\ frac {M_ {d}} {M_ {g}}} = \ venstre [R ^ {2} + {\ frac {\ gamma -1} {2}} M_ {g} ^ {2} (R ^ {2} -1) \ højre] ^ {- {\ frac {1} {2}}}}
Densitetsforholdet er begrænset når Mg→∞{\ displaystyle M_ {g} \ to \ infty}
ρdρg→γ+1γ-1{\ displaystyle {\ frac {\ rho _ {d}} {\ rho _ {g}}} \ til {\ frac {\ gamma +1} {\ gamma -1}}}
For eksempel til luft dannet af diatomiske molekyler, hvor grænsen for densitetsforholdet er lig med 6.
γ=1,4{\ displaystyle \ gamma = 1.4}
Vi kan også beregne den reducerede entropi-variation , hvor C V er volumenvarmekapacitetenSVSV=lnsργ{\ displaystyle {\ frac {S} {C_ {V}}} = \ ln {\ frac {p} {\ rho ^ {\ gamma}}}}
Sd-SgVSV=ln(sdsg)-γln(ρdρg){\ displaystyle {\ frac {S_ {d} -S_ {g}} {C_ {V}}} = \ ln \ left ({\ frac {p_ {d}} {p_ {g}}} \ højre) - \ gamma \ ln \ left ({\ frac {\ rho _ {d}} {\ rho _ {g}}} \ right)}
Denne værdi er nul for M g = 1 og stiger med Mach-tallet. Værdierne af M g <1 fører til negative værdier, der er forbudt af termodynamik. Der er ingen chok for .
Mg≤1{\ displaystyle M_ {g} \ leq 1}
Stationært skråt stød
Et stationært skråt stød kan være til stede i en geometrisk konfiguration af den dihedrale type (se figur). Inden for rammerne af Eulers ligninger ignorerer man de komplekse fænomener med viskøs interaktion i umiddelbar nærhed af muren.
Skriv bevarelsen af momentumdensitetstensoren gennem chokket
(ρVV+s)⋅ikke=ρ(V⊥2V⊥V∥V⊥V∥V∥2)(10)+s(10)=(ρV⊥2+sρV⊥V∥){\ displaystyle (\ rho \ mathbf {V} \ mathbf {V} + p) \ cdot \ mathbf {n} = \ rho {\ begin {pmatrix} V _ {\ perp} ^ {2} & V _ {\ perp} V_ {\ parallel} \\ V _ {\ perp} V _ {\ parallel} & V _ {\ parallel} ^ {2} \ end {pmatrix}} {\ start {pmatrix} 1 \\ 0 \ end {pmatrix}} + p {\ begin {pmatrix} 1 \\ 0 \ end {pmatrix}} = {\ begin {pmatrix} \ rho V _ {\ perp} ^ {2} + p \\\ rho V _ { \ perp} V _ {\ parallel} \ end {pmatrix}}}
hvor er stødet normalt, og komponenterne i henholdsvis V normale og parallelle. Derfor bevarelsesforholdene
ikke=(10){\ displaystyle \ mathbf {n} = {\ begin {pmatrix} 1 \\ 0 \ end {pmatrix}}}
V⊥{\ displaystyle V _ {\ perp}}
V∥{\ displaystyle V _ {\ parallel}}
ρV⊥d2+sd=ρV⊥g2+sg{\ displaystyle \ rho V _ {\ perp d} ^ {2} + p_ {d} = \ rho V _ {\ perp g} ^ {2} + p_ {g}}
ρV⊥dV∥d=ρV⊥gV∥g{\ displaystyle \ rho V _ {\ perp d} V _ {\ parallel d} = \ rho V _ {\ perp g} V _ {\ parallel g}}
Bevarelse af masse er skrevet
ρdV⊥d=ρgV⊥g{\ displaystyle \ rho _ {d} V _ {\ perp d} = \ rho _ {g} V _ {\ perp g}}
Fra de sidste to ligninger udleder man bevarelsen af den parallelle hastighed , som man således vil bemærke V∥d=V∥g{\ displaystyle V _ {\ parallel d} = V _ {\ parallel g}}
V∥{\ displaystyle V _ {\ parallel}}
Systemet er derfor identisk med det med den rette indvirkning for de normale hastigheder og de parallelle hastigheder og Mach-numrene på disse komponenter, og den sidstnævnte mængde er derfor mindre end enhed. Dette foregriber ikke på nogen måde, om strømmen nedstrøms for chokket er subsonisk eller supersonisk. Fra forholdet til det rigtige chok kan vi derfor give
Vgsyndθ{\ displaystyle V_ {g} \ sin \ theta}
Vdsynd(θ-β){\ displaystyle V_ {d} \ sin (\ theta - \ beta)}
Vgcosθ{\ displaystyle V_ {g} \ cos \ theta}
Vdcos(θ-β){\ displaystyle V_ {d} \ cos (\ theta - \ beta)}
Mgsyndθ{\ displaystyle M_ {g} \ sin \ theta}
Mdsynd(θ-β){\ displaystyle M_ {d} \ sin (\ theta - \ beta)}
- udtrykket, der relaterer den ukendte vinkel θ til Mach-nummeret og vinklen på dihedronen β i implicit form
V⊥2V⊥1=V⊥2V∥V∥V⊥1=tan(θ-β)tanθ=(γ-1)Mg2synd2θ+2(γ+1)Mg2synd2θ{\ displaystyle {\ frac {V _ {\ perp 2}} {V _ {\ perp 1}}} = {\ frac {V _ {\ perp 2}} {V _ {\ parallel}}} {\ frac {V _ {\ parallel}} {V _ {\ perp 1}}} = {\ frac {\ tan (\ theta - \ beta)} {\ tan \ theta}} = {\ frac {(\ gamma -1 ) M_ {g} ^ {2} \ sin ^ {2} \ theta +2} {(\ gamma +1) M_ {g} ^ {2} \ sin ^ {2} \ theta}}}
som er forenklet i
tanβ=2tanθMg2synd2θ-1Mg2(γ+cos2θ)+2{\ displaystyle \ tan \ beta = {\ frac {2} {\ tan \ theta}} \, {\ frac {M_ {g} ^ {2} \ sin ^ {2} \ theta -1} {M_ {g } ^ {2} (\ gamma + \ cos 2 \ theta) +2}}}
- De tilsvarende kurver (se figur) viser, at:
- to mulige stødvinkler kan svare til en given tovinklet vinkel og derfor stød med forskellig styrke ("kraften" måles f.eks. ved tryk spring); for et svagt stød (generelt set) fører en stigning i Mach-antallet til et fald i stødvinklen;
- for givet M er der et maksimum på β, der tillader en ret stød, ud over det vises en anden konfiguration: en buet stød foran på dihedronfoden .
- tryk springet
sdsg=1+2γγ+1(Mg2synd2θ-1){\ displaystyle {\ frac {p_ {d}} {p_ {g}}} = 1 + {\ frac {2 \ gamma} {\ gamma +1}} (M_ {g} ^ {2} \ sin ^ { 2} \ theta -1)}
ρ2ρ1=(γ+1)Mg2synd2θ(γ-1)Mg2synd2θ+2{\ displaystyle {\ frac {\ rho _ {2}} {\ rho _ {1}}} = {\ frac {(\ gamma +1) M_ {g} ^ {2} \ sin ^ {2} \ theta } {(\ gamma -1) M_ {g} ^ {2} \ sin ^ {2} \ theta +2}}}
M2=1synd(θ-β)1+γ-12Mg2synd2θγMg2synd2θ-γ-12{\ displaystyle M_ {2} = {\ frac {1} {\ sin (\ theta - \ beta)}} {\ sqrt {\ frac {1 + {\ frac {\ gamma -1} {2}} M_ { g} ^ {2} \ sin ^ {2} \ theta} {\ gamma M_ {g} ^ {2} \ sin ^ {2} \ theta - {\ frac {\ gamma -1} {2}}}} }}
Referencer
-
(in) ® Salas , " The Curious Events Leading to theory of Shock Waves " , 17. Shock Interaction Symposium , Rom,2006( læs online )
-
Pierre-Henri Hugoniot , " Memoir om forplantning af bevægelser i kroppe og især i ideelle gasser (anden del) ", Journal de l'École Polytechnique , vol. 58,1889, s. 1–125 ( læs online )
-
(in) William Rankine , " Om den termodynamiske teori om bølger af endelige langsgående forstyrrelser " , Philosophical Transactions of the Royal Society of London , vol. 160,1870, s. 277–288 ( læs online )
-
Marc Buffat, " relationer gennem en ret chok "
-
(i) John D. Anderson, Jr., Fundamentals of Aerodynamik , New York / St. Louis / Paris osv., McGraw-Hill Education ,1991, 772 s. ( ISBN 0-07-001679-8 )
eksterne links
<img src="https://fr.wikipedia.org/wiki/Special:CentralAutoLogin/start?type=1x1" alt="" title="" width="1" height="1" style="border: none; position: absolute;">