Ajutaj convergent-divergent

Ajutaj convergent–divergent

În dinamica fluidelor și termodinamică un ajutaj convergent–divergent sau ajutaj de Laval este un ajutaj, adică un canal, a cărui secțiune este inițial descrescătoare, iar apoi crescătoare în sensul curgerii. Are aplicații în special pentru curgerea gazelor. În funcție de raportul de presiuni de la capetele sale și de viteza gazului la intrare el poate transforma energia potențială a gazului în energie cinetică (crește viteza de curgere și scade presiunea statică, de exemplu în ajutajele reactive ale motoarelor de rachetă) sau invers (scade viteza de curgere și crește presiunea statică, de exemplu la intrarea în compresoarele turboreactoarelor avioanelor supersonice).[1]

Istoric

Secțiune longitudinală prin motorul de rachetă RD-107 (Muzeul de Stat de Istoria Cosmonauticii „K. Țiolkovski”)

Giovanni Battista Venturi a proiectat tuburi convergent–divergente cunoscute sub denumirea de tuburi Venturi pentru experimente privind efectele de reducere a presiunii fluidului atunci când fluidul curge prin îngustări ale secțiunii (efectul Venturi). Inginerul și inventatorul german Ernst Körting se presupune că în 1878 la injectoarele cu jet de abur a trecut de la un ajutaj convergent la unul convergent–divergent, dar aceste ajutaje au rămas secret industrial.[2] Ulterior, în anul 1888, inginerul suedez Gustaf de Laval a folosit propria sa soluție la ajutajele convergent–divergente ale turbinelor cu abur cu acțiune.[3][4][5][6]

Ajutajul convergent–divergent a fost aplicat pentru prima dată la un motor de rachetă de către Robert H. Goddard. Cele mai multe motoare rachete moderne care realizează propulsia prin jeturi de gaz fierbinte folosesc ajutaje convergent–divergente.

Construcție și funcționare

Îmbinarea părților convergentă și divergentă se face la secțiunea minimă din gâtul ajutajului (indicată în prima imagine printr-o linie punctată), care este, evident, identică pentru ambele părți. Forma părții convergente și a celei divergente se aleg astfel încât pierderile prin frecare, turbulențe, unde de șoc sau rezistențe gazodinamice să fie cât mai nici.

Funcționarea sa se bazează pe diferitele proprietăți ale gazelor care curg. Legea conservării masei face ca de-a lungul ajutajului debitul masic de gaz să fie constant. Deoarece pierderile sunt relativ mici, curgerea este aproape izentropică (entropia gazului rămâne aproape constantă).

Curgerea în interior poate fi subsonică, sonică sau supersonică, în funcție de viteza de la intrare și raportul presiunilor la ieșire și intrare.

Curgere incompresibilă

Dacă în interior viteza gazului nu depășește valoarea M = 0,3 , curgerea poate fi considerată incompresibilă și poate fi calculată cu ecuația lui Bernoulli, valabilă pentru lichide:

ρ w 2 2 + ρ g h + p = c o n s t a n t {\displaystyle \rho {\frac {w^{2}}{2}}+\rho \,g\,h+p=constant}

unde ρ {\displaystyle \rho } este densitatea, w {\displaystyle w} este viteza iar p {\displaystyle p} este presiunea în punctul respectiv, g {\displaystyle g} este accelerația gravitațională iar h {\displaystyle h} este diferența de înălțime între intrare și ieșire.

Curgere compresibilă

Dacă în interior vitezele depășesc valoarea M = 0,3 , efectele compresibilității nu mai pot fi ignorate, iar calculul se face pe baza energiei, prin mijloace termodinamice. Pentru calcul se fac următoarele ipoteze simplificatoare, care nu schimbă semnificativ situația:

  • Gazul este considerat ideal.
  • Curgerea este staționară și se face de-a lungul liniilor de curent.
  • Deoarece ajutajul nu are părți în mișcare, nu există schimb de lucru mecanic cu exteriorul. Deoarece timpul de parcurgere a ajutajului este foarte scurt, schimbul de căldură cu exteriorul este neglijabil. Ca urmare, ajutajul poate fi considerat un sistem adiabatic.
  • Se neglijează pierderile, astfel curgerea este un proces reversibil. Fiind și adiabatică, este un proces izentropic.

Cu aceste ipoteze ecuația de continuitate devine:[1]

d A A + d w w = d v v {\displaystyle {\frac {dA}{A}}+{\frac {dw}{w}}={\frac {dv}{v}}}

unde A este aria secțiunii, iar v este volumul masic al gazului în punctul considerat.

Valoarea funcției / p s i {\displaystyle /psi} în funcție de raportul dintre presiunea într-un anumit punct și presiunea de la intrare, pentru exponentul adiabatic γ = 1 , 4 {\displaystyle \gamma =1,4}

Ecuația de continuitate se scrie de obicei sub forma m ˙ = A ψ = constant , {\displaystyle {\dot {m}}=A\,\psi ={\text{constant}},} [7] unde funcția ψ {\displaystyle \psi } are expresia:[1][8]

ψ = ( p p 1 ) 1 γ γ γ 1 [ 1 ( p p 1 ) γ 1 γ ] {\displaystyle \psi =\left({\frac {p}{p_{1}}}\right)^{\frac {1}{\gamma }}{\sqrt {{\frac {\gamma }{\gamma -1}}\left[1-\left({\frac {p}{p_{1}}}\right)^{\frac {\gamma -1}{\gamma }}\right]}}}

unde γ {\displaystyle \gamma } este exponentul adiabatic al gazului, care pentru gazele biatomice (de exemplu aer) se poate lua 1,4, iar pentru cele triatomice (de exemplu CO2, abur) 1,3.

p1 este presiunea gazului în secțiunea de intrare a ajutajului, iar
p este presiunea gazului în punctul considerat.

Funcția ψ {\displaystyle \psi } are un maxim pentru raportul critic de presiuni:[1][9][7]

p c r p 1 = ( 2 γ + 1 ) γ γ 1 {\displaystyle {\frac {p_{cr}}{p_{1}}}=\left({\frac {2}{\gamma +1}}\right)^{\frac {\gamma }{\gamma -1}}}

valoarea maximă a funcției fiind:[1][7]

ψ m a x = ( 2 γ + 1 ) 1 γ 1 γ γ + 1 . {\displaystyle \psi _{max}=\left({\frac {2}{\gamma +1}}\right)^{\frac {1}{\gamma -1}}{\sqrt {\frac {\gamma }{\gamma +1}}}\,.}

Pentru γ = 1 , 4 {\displaystyle \gamma =1,4} raportul critic de presiuni este ϵ c r = p c r / p 1 = 0 , 528 , {\displaystyle \epsilon _{cr}=p_{cr}/p_{1}=0,528,} pentru care ψ = 0 , 483. {\displaystyle \psi =0,483.}
Pentru γ = 1 , 3 {\displaystyle \gamma =1,3} raportul critic de presiuni este ϵ c r = p c r / p 1 = 0 , 546 , {\displaystyle \epsilon _{cr}=p_{cr}/p_{1}=0,546,} pentru care ψ = 0 , 472. {\displaystyle \psi =0,472.}

Pentru o curgere izentropică, viteza gazului la ieșire este:[1]

w 2 = 2 ( h 1 h 2 ) + w 1 2 {\displaystyle w_{2}={\sqrt {2(h_{1}-h_{2})+w_{1}^{2}}}}

unde h 1 , h 2 {\displaystyle h_{1},\,h_{2}} sunt entalpiile gazului la intrare, respectiv la ieșire, iar w 1 {\displaystyle w_{1}} este viteza gazului la intrare.

Exprimată în funcție de presiuni, relația precedentă devine:[10][11][12][13][14]

w 2 = 2 γ γ 1 p 1 v 1 [ 1 ( p 2 p 1 ) γ 1 γ ] + w 1 2 {\displaystyle w_{2}={\sqrt {{\frac {2\gamma }{\gamma -1}}p_{1}v_{1}\left[1-\left({\frac {p_{2}}{p_{1}}}\right)^{\frac {\gamma -1}{\gamma }}\right]+w_{1}^{2}}}}

unde v 1 {\displaystyle v_{1}} este volumul masic al gazului la intrare. Relația este cunoscută sub numele de formula lui Saint-Venant.[11]

Curgerea din ajutaj fiind izentropică, temperatura depinde de presiune sau de volumul specific conform relației:

T T 1 = ( p p 1 ) γ 1 γ = ( v v 1 ) 1 γ 1 . {\displaystyle {\frac {T}{T_{1}}}=\left({\frac {p}{p_{1}}}\right)^{\frac {\gamma -1}{\gamma }}=\left({\frac {v}{v_{1}}}\right)^{\frac {1}{\gamma -1}}\,.}

Cazul unei viteze subsonice la intrare

Este cazul cel mai intuitiv, în care curgerea se face dintr-o zonă cu presiune mai mare spre o zonă cu presiune mai mică.

Cazul raportului de presiuni mai mare decât raportul critic

Dacă raportul presiunilor de la ieșire și intrare este mai mare decât raportul critic de presiuni ( 1 > ϵ = p 2 / p 1 > ϵ c r ) , {\displaystyle (1>\epsilon =p_{2}/p_{1}>\epsilon _{cr}),} punctul de funcționare al ajutajului este pe ramura din dreapta a curbei ψ {\displaystyle \psi } (v. figura, curba roșie continuă), iar curgerea este subsonică peste tot, inclusiv în secțiunea minimă. În porțiunea convergentă a ajutajului presiunea scade, iar viteza și volumul masic ale gazului cresc. În porțiunea divergentă presiunea continuă să scadă, volumul masic crește, dar, datorită creșterii secțiunii, viteza gazului scade. Datorită creșterii secțiunii liniile de curent devin divergente, iar dacă această divergență este prea mare, fluidul se poate desprinde de suprafața difuzorului, jetul îngustându-se. Desprinderea este însoțită de variații de presiune relativ mari, ceea ce produce turbulență, respectiv disipare. Se recomandă ca unghiul de evazare al difuzorului să fie de 10–30°, valorile mari fiind acceptate la rapoarte de presiuni apropiate de 1, cum apar la instalațiile de ventilare. Odată cu scăderea raportului presiunilor (ceea ce determină creșterea vitezei de curgere), unghiul de evazare trebuie micșorat.

Cazul raportului de presiuni mai mic decât raportul critic

Variația vitezei (v), a presiunii (p) și a temperaturii (T) la curgerea printr-un ajutaj convergent–divergent

Dacă raportul presiunilor de la ieșire și intrare este mai mic decât raportul critic de presiuni ( ϵ c r > ϵ = p 2 / p 1 > 0 ) , {\displaystyle (\epsilon _{cr}>\epsilon =p_{2}/p_{1}>0),} punctul de funcționare al ajutajului este pe ramura din stânga a curbei ψ {\displaystyle \psi } (v. figura, curba roșie punctată).

În figura de alături intrarea gazului este în stânga, iar ieșirea în dreapta. În partea convergentă presiunea scade, în secțiunea minimă realizându-se presiunea corespunzătoare raportului critic de presiuni. Viteza gazului, inițial subsonică ( M < 1 {\displaystyle M<1} ), crește până la la cea a sunetului ( M = 1 {\displaystyle M=1} ), deoarece o porțiune convergentă nu poate realiza mai mult. În porțiunea divergentă (difuzor) viteza devine supersonică ( M > 1 {\displaystyle M>1} ), continuă să crească, iar presiunea continuă să scadă conform secțiunii difuzorului și a relației Saint-Venant. La ieșire pot să apară următoarele cazuri: [15]

1 - destindere în afara ajutajului; 2 - destindere optimă în ajutaj; 3 - contracție la ieșirea din ajutaj; 4 - contracție în interiorul difuzorului
  • Presiunea de la ieșire este mai mare decât cea a mediului exterior, p 2 > p e x t , {\displaystyle p_{2}>p_{ext},} caz în care jetul continuă să se destindă turbulent în afara ajutajului, energia sa remanentă disipându-se în frecări, unde de șoc și sonore. Situația este caracteristică pentru ajutajele reactive ale rachetelor în atmosfera superioară, unde presiunea este foarte mică.
  • Presiunea de la ieșire este exact cea a mediului exterior, p 2 = p e x t , {\displaystyle p_{2}=p_{ext},} cazul optim, în care este utilă întreaga energie a jetului în condițiile date. Jetul are formă cilindrică, iar viteza sa este disipată în frecări și unde de șoc și sonore în afara ajutajului.
  • Presiunea de la ieșire este mai mică decât cea a mediului exterior, p 2 < p e x t , {\displaystyle p_{2}<p_{ext},} caz în care jetul suferă o contracție cu o creștere de presiune. Contracția poate să înceapă în secțiunea de ieșire sau mai devreme, într-un punct din interiorul difuzorului, sub forma desprinderii de peretele difuzorului. Desprinderea se face cu creșteri bruște de presiune, ceea ce duce la turbulențe și unde de șoc care reduc caracterul izentropic, adică performanțele ajutajului.

Cazul unei viteze supersonice la intrare

În cazul că viteza gazului la intrarea într-un ajutaj convergent-divergent este supersonică, fenomene din el se petrec oarecum invers, însă curgerea necesită îndeplinirea unor anumite condiții.

Se numește presiune de stagnare p0 suma dintre presiunea statică p și cea dinamică:

p 0 = p + ρ w 2 2 {\displaystyle p_{0}=p+\rho {\frac {w^{2}}{2}}}

unde ρ {\displaystyle \rho } este densitatea, iar w {\displaystyle w} este viteza gazului în punctul respectiv.

Deoarece curgerea se face de la presiune mare la presiune mică, presiunea de stagnare de la intrare, p01, trebuie să fie mai mare decât presiunea de stagnare de la ieșire, p02 adică:

p 01 = p 1 + ρ 1 w 1 2 2 > p 2 + ρ 2 w 2 2 2 = p 02 . {\displaystyle p_{01}=p_{1}+\rho _{1}{\frac {w_{1}^{2}}{2}}>p_{2}+\rho _{2}{\frac {w_{2}^{2}}{2}}=p_{02}.}

Practic densitatea și viteza la intrare trebuie să fie suficient de mari.

Dacă raportul presiunilor de stagnare la ieșire și intrare este mai mare decât raportul critic de presiuni ( 1 > ϵ = p 02 / p 01 > ϵ c r ) , {\displaystyle (1>\epsilon =p_{02}/p_{01}>\epsilon _{cr}),} punctul de funcționare al ajutajului este pe ramura din dreapta a curbei ψ {\displaystyle \psi } (v. figura, curba roșie continuă), iar curgerea este supersonică peste tot, inclusiv în secțiunea minimă. Inițial, în partea convergentă, viteza și volumul masic scad, iar presiunea și temperatura cresc. În secțiunea minimă viteza este încă supersonică, iar după secțiunea minimă presiunea și temperatura scad, iar viteza și volumul masic cresc. Acest regim nu are aplicații tehnice.

Dacă raportul presiunilor de stagnare de la intrare și ieșire este mai mic decât raportul critic de presiuni ( ( ϵ c r > ϵ = p 02 / p 01 > 0 ) , {\displaystyle (\epsilon _{cr}>\epsilon =p_{02}/p_{01}>0),} punctul de funcționare al ajutajului este pe ramura din stânga a curbei ψ {\displaystyle \psi } (v. figura, curba roșie punctată).

În partea convergentă presiunea crește, în secțiunea minimă realizându-se presiunea corespunzătoare raportului critic de presiuni. Viteza gazului, inițial supersonică ( M > 1 {\displaystyle M>1} ), scade până la la cea a sunetului ( M = 1 {\displaystyle M=1} ). În porțiunea divergentă (difuzor) viteza devine subsonică ( M < 1 {\displaystyle M<1} ), continuă să scadă, iar presiunea continuă să crească conform secțiunii difuzorului și a relației Saint-Venant. Alura acestor parametri este inversă față de cea din figura de la începutul secțiunii cazul raportului de presiuni mai mic decât raportul critic (ca atunci când intrarea ar fi în dreapta, iar ieșirea în stânga).

Creșterea presiunii împiedică apariția desprinderilor în difuzor. Însă pot să apară unde de șoc în porțiunea convergentă, efect folosit la comprimarea dinamică a aerului în prizele de aer ale avioanelor supersonice.

Debite și forțe

Cunoscându-se parametrii de stare și viteza (numărul Mach) în secțiunea cu aria A, debitul masic de gaz se poate calcula cu relația:[16]

m ˙ = A p 0 T γ M R M a ( 1 + γ 1 2 M a 2 ) γ + 1 2 ( γ 1 ) {\displaystyle {\dot {m}}={\frac {Ap_{0}}{\sqrt {T}}}\cdot {\sqrt {\frac {\gamma M}{R}}}\cdot \mathrm {Ma} \cdot (1+{\frac {\gamma -1}{2}}\mathrm {Ma} ^{2})^{-{\frac {\gamma +1}{2(\gamma -1)}}}}

unde M este masa molară a gazului, iar R este constanta universală a gazului ideal.

Dacă curgerea în secțiunea minimă este sonică (Ma = 1) relația de mai sus se siplifică:

m ˙ = A p 0 T γ M R ( γ + 1 2 ) γ + 1 2 ( γ 1 ) {\displaystyle {\dot {m}}={\frac {Ap_{0}}{\sqrt {T}}}\cdot {\sqrt {\frac {\gamma M}{R}}}\cdot ({\frac {\gamma +1}{2}})^{-{\frac {\gamma +1}{2(\gamma -1)}}}}

De asemenea, se poate folosi relația dedusă din formula Saint-Venant:[10]

m ˙ = A c r γ p 1 v 1 ( 2 γ + 1 ) γ + 1 γ 1 . {\displaystyle {\dot {m}}=A_{cr}{\sqrt {\gamma {\frac {p_{1}}{v_{1}}}\left({\frac {2}{\gamma +1}}\right)^{\frac {\gamma +1}{\gamma -1}}}}.}

Pentru curgeri reale, neizentropice, vitezele, respectiv debitele, trebuie micșorate cu 2...7 %.[10]

Forța exercitată de fluid asupra ajutajului este:

F = m ˙ w 2 {\displaystyle F={\dot {m}}\cdot w_{2}}

În aerodinamică această forță se numește tracțiune și de obicei se raportează la accelerația gravitațională standard.[17]

Aplicații

Principala aplicație a ajutajelor convergent-divergente este realizarea curgerilor gazelor cu viteze supersonice. Aceste viteze sunt necesare în ajutajele reactive ale avioanelor supersonice și rachetelor, pentru generarea jeturilor care creează forțe de propulsie.

Tot în aviație, dispozitivele de admisie a aerului ale avioanelor supersonice sunt realizate sub formă de ajutaje convergent-divergente, care comprimă dinamic aerul înainte de a intra în compresoarele turbomotoarelor.

O altă aplicație a acestor ajutaje este la turbinele de abur cu acțiune de Laval, la care ajutajele care creează jeturi care acționează paletele funcționează la diferențe mari de entalpie, care corespund vitezelor supersonice.

Altă aplicație era la injectoarele de apă (pompe cu jet de abur) care introduceau apă în cazanele locomotivelor cu abur.

Note

  1. ^ a b c d e f Alexandru Dănescu ș.a., Lexicon de termodinamică și mașini termice, vol. 1 A–E, București, Editura Tehnică, 1985, pp. 38–39
  2. ^ en Krehl, Peter O. K. (). History of Shock Waves, Explosions and Impact: A Chronological and Biographical Reference. Springer. ISBN 9783540304210. Arhivat din original la . Accesat în . 
  3. ^ A se vedea:
    • Patentul belgian no. 83,196 (din 29 septembrie 1888)
    • Patentul englez no. 7143 (din 29 aprilie 1889)
    • de Laval, Carl Gustaf Patrik, "Steam turbine" Arhivat în , la Wayback Machine. U.S. Patent no. 522,066 (filed: 1889 May 1; issued: 1894 June 26)
  4. ^ en Theodore Stevens and Henry M. Hobart (). Steam Turbine Engineering. MacMillan Company. pp. 24–27.  Available on-line here Arhivat în , la Wayback Machine. in Google Books.
  5. ^ en Robert M. Neilson (). The Steam Turbine. Longmans, Green, and Company. pp. 102–103.  Available on-line here in Google Books.
  6. ^ en Garrett Scaife (). From Galaxies to Turbines: Science, Technology, and the Parsons Family. Taylor & Francis Group. p. 197.  Available on-line here Arhivat în , la Wayback Machine. in Google Books.
  7. ^ a b c Vlădea, 1974, p. 379
  8. ^ Vlădea, 1974, p. 378
  9. ^ Popa ș.a., MIT1, p. 59
  10. ^ a b c Popa ș.a., MIT1, p. 60
  11. ^ a b Vlădea, 1974, p. 377–379
  12. ^ en „Richard Nakka's Equation 12”. Arhivat din original la . Accesat în . 
  13. ^ en „Robert Braeuning's Equation 1.22”. Arhivat din original la . Accesat în . 
  14. ^ en George P. Sutton (). Rocket Propulsion Elements: An Introduction to the Engineering of Rockets (ed. 6th). Wiley-Interscience. p. 636. ISBN 0-471-52938-9. 
  15. ^ Vlădea, 1974, pp. 386–389
  16. ^ en Hall, Nancy. „Mass Flow Choking”. NASA. Arhivat din original la . Accesat în . 
  17. ^ Lucian Miclăuș, Glosar aviatic, Timișoara: Editura Marineasa, 2001, ISBN: 978-9485-99-5. „copie arhivă” (PDF). Arhivat din original (PDF) la . Accesat în . 

Bibliografie

Vezi și

Legături externe

Portal icon Portal Fizică
Portal icon Portal Inginerie mecanică
  • en Exhaust gas velocity calculator[nefuncționalăarhivă]
  • en Other applications of nozzle theory Flow of gases and steam through nozzles