Curgere potențială

curgere în care vitezele sunt gradientul unei funcții scalare
Pentru curgerea potențială în jurul unui cilindru, vedeți Curgerea potențială în jurul unui cilindru circular.

În dinamica fluidelor, curgerea potențială sau curgerea irotațională se referă la descrierea unei curgeri a unui fluid fără vârtejuri. Această descriere este tipică pentru cazul în care viscozitatea este neglijabilă, adică pentru un fluid ideal și în absența vârtejurilor.[1]

Linii de curent⁠(d) ale curgerii potențiale în jurul unui profil aerodinamic NACA 0012⁠(d) la un unghi de atac⁠(d) de 11°, cu straturile de curgere superioare și inferioare identificate. Curgerea este bidimensională, iar profilul aerodinamic are o lungime infinită.

Curgerea potențială descrie câmpul de viteze ca fiind gradientul unei funcții scalare: potențialul vitezei. Ca rezultat, o curgere potențială este caracterizată de un câmp de viteze irotațional⁠(d), care este o aproximare validă pentru diverse aplicații. Irotaționalitatea curgerii potențiale derivă din faptul că gradientul unui scalar⁠(d) este întotdeauna egal cu zero.

În cazul unei curgeri incompresibile, potențialul vitezei satisface ecuația lui Laplace, iar teoria potențialului⁠(d) este aplicabilă. Cu toate acestea, curgerile potențiale au fost utilizate și pentru a descrie curgerile compresibile⁠(d) și curgerile Hele-Shaw⁠(d). Abordarea curgerii potențiale apare în modelarea atât a curgerilor staționare, cât și a celor nestaționare.

Aplicațiile curgerii potențiale includ: câmpul de curgere exterior pentru profile aerodinamice, valuri, electroosmoză și curgerea apelor subterane⁠(d). Pentru curgeri (sau părți ale acestora) cu efecte puternice de vârtej, aproximarea curgerii potențiale nu este aplicabilă. În regiunile de curgere în care se știe că vârtejul este important, cum ar fi urmele valurilor⁠(d) și straturile limită⁠(d), teoria curgerii potențiale nu poate furniza predicții rezonabile ale curgerii.[2] Din fericire, există adesea regiuni mari ale unei curgeri în care ipoteza de irotaționalitate este valabilă, motiv pentru care curgerea potențială este utilizată pentru diverse aplicații, cum ar fi: aripile avioanelor, curgerea apelor subterane⁠(d), acustică, valuri și electroosmoză.[3]

Descriere și caracteristici

modificare
 
O curgere potențială este construită prin adăugarea de curgeri elementare⁠(d) simple și observarea rezultatului.
 
Linii de curent⁠(d) pentru curgerea potențială incompresibilă în jurul unui cilindru circular⁠(d), într-o curgere uniformă.

În curgerea potențială sau irotațională, câmpul vectorial al vârtejului este zero, adică

 

unde   este câmpul de viteze și   este câmpul de vârtej. Ca orice câmp vectorial cu rotație nulă, câmpul de viteză poate fi exprimat ca gradientul unui anumit scalar, să zicem   care se numește potențialul vitezei, deoarece rotația gradientului este întotdeauna zero. Prin urmare, avem[4][5]

 

Potențialul vitezei nu este definit în mod unic, deoarece i se poate adăuga o funcție arbitrară de timp, să zicem  , fără a afecta cantitatea fizică relevantă care este  . Neunicitatea este de obicei eliminată prin selectarea adecvată a condițiilor inițiale sau la limită satisfăcute de   și, ca atare, procedura poate varia de la o problemă la alta.

În curgerea potențială, circulația   în jurul oricărui contur simplu conex   este zero. Acest lucru poate fi demonstrat folosind teorema lui Stokes,

 

unde   este elementul de linie pe contur și   este elementul de suprafață al oricărei suprafețe delimitate de contur. În spațiul multiplu conectat (să zicem, în jurul unui contur care înconjoară un corp solid bidimensional sau în jurul unui contur care înconjoară un tor tridimensional) sau în prezența vârtejurilor concentrate (să zicem, în așa-numitele vârtejuri irotaționale⁠(d) sau vârtejuri punctuale sau în inele de fum), circulația   nu trebuie să fie zero. În primul caz, teorema lui Stokes nu se poate aplica, iar în cazul ulterior,   este diferit de zero în regiunea delimitată de contur. În jurul unui contur care înconjoară un cilindru solid infinit cu care conturul se învârte de   ori, avem

 

unde   este o constantă ciclică. Acest exemplu aparține unui spațiu dublu conectat. Într-un spațiu tuplu conectat  , există   astfel de constante ciclice, și anume,  

Curgere incompresibilă

modificare

În cazul unei curgeri incompresibile — de exemplu, a unui lichid sau a unui gaz la numere Mach mici; dar nu pentru undele sonore — viteza v are divergență zero:[4]

 

Substituind aici   rezultă că   satisface ecuația lui Laplace[4]

 

unde 2 = ∇ ⋅ ∇ este operatorul Laplace (uneori scris și Δ). Deoarece soluțiile ecuației lui Laplace sunt funcții armonice, fiecare funcție armonică reprezintă o soluție de curgere potențială. După cum este evident, în cazul incompresibilității, câmpul de viteze este determinat complet din cinematica sa: ipotezele de irotaționalitate și divergență zero a curgerii. Dinamica, în legătură cu ecuațiile de impuls, trebuie aplicată doar ulterior, dacă cineva este interesat să calculeze câmpul de presiune: de exemplu, pentru curgerea în jurul profilurilor aerodinamice prin utilizarea legii lui Bernoulli.

În curgerile incompresibile, contrar concepției greșite comune, curgerea potențială satisface într-adevăr ecuațiile complete Navier-Stokes, nu doar ecuațiile lui Euler, deoarece termenul vâscos

 

este identic cu zero. Este incapacitatea curgerii potențiale de a satisface condițiile limită necesare, în special în apropierea frontierelor solide, ceea ce o invalidează în reprezentarea câmpului de curgere necesar. Dacă curgerea potențială satisface condițiile necesare, atunci este soluția necesară a ecuațiilor Navier–Stokes incompresibile.

Când este bidimensională, cu ajutorul funcției armonice 𝜑 și a funcției sale armonice conjugate 𝜓 (funcția de curent), curgerea potențială incompresibilă se reduce la un sistem foarte simplu care este analizat folosind analiză complexă (vezi mai jos).

Curgere compresibilă

modificare

Curgere staționară

modificare

Teoria curgerii potențiale poate fi utilizată și pentru modelarea curgerii compresibile irotaționale. Derivarea ecuației de guvernare pentru   din ecuația lui Euler este destul de simplă. Ecuațiile de continuitate și de impuls (pentru curgere potențială) pentru curgeri staționare sunt date de

 

unde ultima ecuație rezultă din faptul că entropia este constantă pentru un element de fluid și că pătratul vitezei sunetului este  . Eliminarea lui   din cele două ecuații de guvernare rezultă în:

 

Versiunea incompresibilă apare la limita  . Substituind aici   rezultă[6][7]

 

unde   este exprimat ca funcție a magnitudinii vitezei  . Pentru un gaz politropic⁠(d),  , unde   este coeficientul de transformare adiabatică și   este entalpia de stagnare⁠(d). Bidimensional, ecuația se simplifică la

 

Valabilitate: În forma actuală, ecuația este valabilă pentru orice curgere potențială ideală, indiferent dacă curgerea este subsonică sau supersonică (de exemplu, curgerea Prandtl-Meyer⁠(d)). Cu toate acestea, în curgeri supersonice și transonice pot apărea unde de șoc care pot introduce entropie și vârtej în curgere, făcând curgerea rotațională. Cu toate acestea, există două cazuri în care curgerea potențială prevalează chiar și în prezența undelor de șoc, care sunt explicate din ecuația de impuls (nu neapărat potențială) scrisă în următoarea formă

 

unde   este entalpia specifică,   este câmpul de vârtej,   este temperatura și   este entropia specifică. Deoarece în fața undei de șoc de vârf avem o curgere potențială, ecuația lui Bernoulli arată că   este constantă, ceea ce este de asemenea constantă pe întreaga undă de șoc (condițiile Rankine-Hugoniot⁠(d)) și prin urmare putem scrie[6]

 

1) Când unda de șoc este de intensitate constantă, discontinuitatea entropiei de-a lungul undei de șoc este, de asemenea, constantă, adică   și, prin urmare, producția de vârtej este zero. Undele de șoc la marginea ascuțită a unei pene bidimensionale sau a unui con tridimensional (curgerea Taylor-Maccoll⁠(d)) au intensitate constantă. 2) Pentru undele de șoc slabe, saltul de entropie de-a lungul undei de șoc este o cantitate de ordinul trei în termeni de intensitatea undei de șoc și, prin urmare,   poate fi neglijat. Undele de șoc în corpurile subțiri sunt aproape paralele cu corpul și sunt slabe.

Curgeri aproape paralele: Când curgerea este predominant unidirecțională cu abateri mici, cum ar fi în curgerea trecută de corpuri subțiri, ecuația completă poate fi simplificată în continuare. Fie   curentul principal și să luăm în considerare abaterile mici de la acest câmp de viteză. Potențialul vitezei corespunzător poate fi scris ca   unde   caracterizează abaterea mică de la curgerea uniformă și satisface versiunea liniarizată a ecuației complete. Aceasta este dată de

 

unde   este numărul Mach constant care corespunde curgerii uniforme. Această ecuație este valabilă cu condiția ca   să nu fie aproape de unitate. Când   este mic (curgere transonică), avem următoarea ecuație neliniară

 

Unde   este valoarea critică a derivatei Landau⁠(d)  [8][9] și   este volumul specific. Curgerea transonică este complet caracterizată de parametrul unic  , care pentru gazul politropic ia valoarea  . Sub transformarea hodografă⁠(en)[traduceți], ecuația transonică bidimensională devine ecuația Euler–Tricomi⁠(d).

Curgere nestaționară

modificare

Ecuațiile de continuitate și de impuls (pentru curgere potențială) pentru curgeri nestaționare sunt date de

 

Integrala primă a ecuației de impuls (pentru curgere potențială) este dată de

 

unde   este o funcție arbitrară. Fără a pierde generalitatea, putem stabili   deoarece   nu este definit în mod unic. Combinând aceste ecuații, obținem

Substituind aici   rezultă

 

Curgeri aproape paralele: Ca și înainte, pentru curgeri aproape paralele, putem scrie (după introducerea unui timp recalculat  )

 

cu condiția ca numărul Mach constant   să nu fie aproape de unitate. Când   este mic (curgere transonică), avem următoarea ecuație neliniară[6]

 

Unde sonore: În undele sonore, magnitudinea vitezei   (sau numărul Mach) este foarte mică, deși termenul nestaționar este acum comparabil cu ceilalți termeni principali din ecuație. Astfel, neglijând toți termenii pătratici și de ordin superior și observând că în aceeași aproximare,   este o constantă (de exemplu, în gazul politropic  ), avem[10][6]

 

care este o ecuație de undă liniară pentru potențialul vitezei φ. Din nou, partea oscilatorie a vectorului de viteză v este legată de potențialul vitezei prin v = ∇φ, în timp ce ca înainte Δ este operatorul Laplace și c este viteza medie a sunetului în mediul omogen⁠(d). De reținut că și părțile oscilatorii ale presiunii p și densității ρ satisfac individual ecuația de undă, în această aproximare.

Aplicabilitate și limitări

modificare

Curgerea potențială nu include toate caracteristicile curgerilor întâlnite în lumea reală. Teoria curgerii potențiale nu poate fi aplicată pentru curgeri vâscoase⁠(d) interne,[2] cu excepția curgerilor dintre plăci apropiate⁠(d). Richard Feynman considera curgerea potențială atât de nefizică încât singurul fluid care să respecte ipotezele era „apa uscată” (citat după John von Neumann).[11] Curgerea potențială incompresibilă face, de asemenea, o serie de predicții invalide, cum ar fi paradoxul lui d'Alembert⁠(d), care afirmă că rezistența la înaintare a oricărui obiect care se deplasează printr-un fluid infinit altfel în repaus este zero.[12] Mai precis, curgerea potențială nu poate explica comportamentul curgerilor care includ un strat limită⁠(d).[2] Cu toate acestea, înțelegerea curgerii potențiale este importantă în multe ramuri ale mecanicii fluidelor. În special, curgerile potențiale simple (numite curgeri elementare⁠(d)), precum vârtejul liber⁠(d) și sursa punctiformă dispun de soluții analitice rapide. Aceste soluții pot fi suprapuse pentru a crea curgeri mai complexe care satisfac o varietate de condiții de frontieră. Aceste curgeri corespund îndeaproape curgerilor reale din întreaga mecanică a fluidelor; în plus, multe perspective valoroase apar atunci când se consideră abaterea (adesea mică) dintre o curgere observată și curgerea potențială corespunzătoare. Curgerea potențială găsește multe aplicații în domenii precum proiectarea aeronavelor. De exemplu, în dinamica fluidelor numerice, o tehnică constă în cuplarea unei soluții de curgere potențială în afara stratului limită cu o soluție a ecuațiilor stratului limită⁠(d) în interiorul Strat limită⁠(d). Absența efectelor stratului limită înseamnă că orice linie de curent poate fi înlocuită cu o frontieră solidă fără nicio schimbare în câmpul de curgere, o tehnică utilizată în multe abordări de proiectare aerodinamică. O altă tehnică ar fi utilizarea solidelor Riabouchinsky⁠(d).

Analiză pentru curgerea incompresibilă bidimensională

modificare

Curgerea potențială bidimensională este simplă de analizat folosind transformările conforme, prin utilizarea transformărilor planului complex. Totuși, spre deosebire de analiza tradițională a curgerii fluidelor în jurul unui cilindru, utilizarea numerelor complexe nu este necesară. Nu este posibil să se rezolve o curgere potențială tridimensională folosind numere complexe.[13]

Ideea de bază este de a utiliza o funcție olomorfă (numită și analitică) sau meromorfă f, care face o conexiune între domeniul fizic (x, y) la domeniul transformat conform (φ, ψ). În timp ce x, y, φ și ψ sunt toate cu valori reale, este convenabil să definim cantitățile complexe

 

Acum, dacă scriem funcția f ca[13]

 

Atunci, deoarece f este o funcție olomorfă sau meromorfă, trebuie să satisfacă ecuațiile Cauchy-Riemann[13]

 

Componentele vitezelor (u, v), în direcțiile (x, y) respectiv, pot fi obținute direct din f prin diferențiere în raport cu z[13]

 

Deci câmpul de viteze v = (u, v) este specificat de[13]

 

Atât φ cât și ψ satisfac ecuația lui Laplace:[13]

 

Deci φ poate fi identificat ca potențialul vitezei, iar ψ este numită funcția curentului⁠(d).[13] Liniile de constantă ψ sunt cunoscute ca linii de curent⁠(d), iar liniile de constantă φ sunt numite linii echipotențiale (vezi suprafața echipotențială⁠(d)).

Liniile de curent și liniile echipotențiale sunt ortogonale între ele, deoarece:[13]

 

Astfel, curgerea are loc de-a lungul liniilor de constantă ψ și perpendicular pe liniile de constantă φ.[13]

De asemenea, este îndeplinită relația Δψ = 0, care este echivalentă cu ∇ × v = 0. Prin urmare, curgerea este irotațională. Condiția automată 2Ψ/xy = 2Ψ/yx ne oferă atunci constrângerea de incompresibilitate ∇ · v = 0.

Exemple de curgeri incompresibile bidimensionale

modificare

Orice funcție diferențiabilă poate fi utilizată pentru f. Exemplele următoare folosesc o varietate de funcții elementare⁠(d); pot fi utilizate și funcții speciale⁠(d). De reținut că funcțiile cu mai multe valori⁠(d), cum ar fi logaritmul natural, pot fi utilizate, dar atenția trebuie limitată la o singură suprafață Riemann⁠(d).

Legile puterii

modificare
Exemple de hărți conforme pentru legea de putere w = Azn
Exemple de hărți conforme pentru legea puterii w = Azn, pentru diferite valori ale puterii n. Este prezentat planul z, arătând liniile de potențial constant φ și funcția de curent ψ, în timp ce w = φ + .

În cazul în care se aplică următoarea transformare conformă de tip lege de putere, de la z = x + iy la w = φ + :[14]

 

atunci, scriind z în coordonate polare ca z = x + iy = re, avem[14]

 

În figurile din dreapta sunt prezentate exemple pentru câteva valori ale lui n. Linia neagră este limita curgerii, în timp ce liniile albastre mai închise sunt linii de curent, iar liniile albastre mai deschise sunt liniile echipotențiale. Unele puteri interesante n sunt:[14]

  • n = 1/2: corespunde cu curgerea în jurul unei plăci semi-infinite,
  • n = 2/3: curgere în jurul unui colț drept,
  • n = 1: un caz trivial de curgere uniformă,
  • n = 2: curgere printr-un colț sau în apropierea unui punct de stagnare, și
  • n = −1: curgere datorată unui dublet sursă

Constanta A este un parametru de scalare: modulul | A | determină scara, în timp ce argumentul⁠(d) său arg(A) introduce o rotație (dacă este diferită de zero).

Legile puterii cu n = −1: curgere uniformă

modificare

Dacă w = Az1, adică o lege de putere cu n = 1, liniile de curent (adică liniile de constantă ψ) sunt un sistem de linii drepte paralele cu axa  . Acest lucru este cel mai ușor de observat scriind în termeni de componente reale și imaginare:

 

astfel obținând φ = Ax și ψ = Ay. Această curgere poate fi interpretată ca o curgere uniformă paralelă cu axa  .

Legile puterii cu n = −2

modificare

Dacă n = 2, atunci w = Az2 și linia de curent corespunzătoare unei anumite valori a lui ψ sunt acele puncte care satisfac

 

care reprezintă un sistem de hiperbole rectangulare. Acest lucru poate fi demonstrat prin rescrierea din nou în termeni de componente reale și imaginare. Observând că sin 2θ = 2 sin θ cos θ și rescriind sin θ = y/r și cos θ = x/r se observă (prin simplificare) că liniile de curent sunt date de

 

Câmpul de viteze este dat de φ, sau

 

În dinamica fluidelor, câmpul de curgere din apropierea originii corespunde unui punct de stagnare⁠(d). De reținut că fluidul din origine este în repaus (ceea ce rezultă din diferențierea lui f(z) = z2 la z = 0). Linia de curent ψ = 0 este deosebit de interesantă: are două (sau patru) ramuri, urmând axele de coordonate, adică x = 0 și y = 0. Deoarece nu curge fluid prin axa  , aceasta (axa  ) poate fi tratată ca o frontieră solidă. Astfel, este posibil să ignorăm curgerea în semiplanul inferior unde y < 0 și să ne concentrăm asupra curgerii în semiplanul superior. Cu această interpretare, curgerea este cea a unui jet direcționat vertical care lovește o placă plană orizontală. De asemenea, curgerea poate fi interpretată ca fiind o curgere într-un colț de 90 de grade, dacă se ignoră regiunile specificate de (să zicem) x, y < 0.

Legile puterii cu n = 3

modificare

Dacă n = 3, curgerea rezultantă este o versiune hexagonală a cazului n = 2 analizat anterior. Liniile de curent sunt date de ψ = 3x2yy3, iar curgerea în acest caz poate fi interpretată ca o curgere într-un colț de 60°.

Legile puterii cu n = −1: dublet

modificare

Dacă n = −1, liniile de curent sunt date de:

 

Această ecuație este mai ușor de interpretat în termeni de componente reale și imaginare:

 

Astfel, liniile de curent sunt cercuri care sunt tangente la axa Ox în origine. Cercurile din semiplanul superior curg astfel în sensul acelor de ceasornic, iar cele din semiplanul inferior curg în sens invers acelor de ceasornic. De remarcat că componentele vitezei sunt proporționale cu r−2; iar valorile lor în origine sunt infinite. Acest model de curgere este de obicei numit dublet sau dipol și poate fi interpretat ca o combinație a unei perechi sursă-scurgere de intensitate infinită, menținută la o distanță infinitezimală. Câmpul de viteze este dat de:

 

sau în coordonate polare:

 

Legile puterii cu n = −2: cvadripolar

modificare

Dacă n = −2, liniile de curent sunt date de:

 

Acest câmp de curgere este asociat cu un cvadripolar⁠(d).[15]

Sursa și scurgerea liniară

modificare

O sursă sau scurgere liniară de intensitate   (  pentru sursă și   pentru scurgere) este dată de potențialul:

 

unde   reprezintă, de fapt, curgerea volumică pe unitate de lungime printr-o suprafață care înconjoară sursa sau scurgerea. Câmpul de viteză în coordonate polare este:

 

adică o curgere pur radială.

Vârtejul liniar

modificare

Un vârtej liniar de intensitatea   este dat de:

 

unde   este circulația în jurul oricărui contur simplu închis care înconjoară vârtejul. Câmpul de viteză în coordonate polare este:

 

adică o curgere pur azimutală.

Analiză pentru curgeri incompresibile tridimensionale

modificare

Pentru curgeri tridimensionale, potențialul complex nu poate fi obținut.

Sursă și scurgerea punctuală

modificare

Potențialul vitezei al unei surse sau scurgeri punctuale de intensitate   (  pentru sursă și   pentru scurgere) în coordonate polare sferice este dat de:

 

unde   reprezintă, de fapt, curgerea volumică printr-o suprafață închisă care înconjoară sursa sau scurgerea. Câmpul de viteză în coordonate polare sferice este:

 
  1. ^ „Fluide Ideale” (PDF). 
  2. ^ a b c Batchelor (1973) pp. 378–380.
  3. ^ Kirby, B.J. (), Micro- and Nanoscale Fluid Mechanics: Transport in Microfluidic Devices., Cambridge University Press, ISBN 978-0-521-11903-0, arhivat din original la , accesat în  
  4. ^ a b c Batchelor (1973) pp. 99–101.
  5. ^ „Turbulenţa | PDF” (în engleză). Scribd. Accesat în . 
  6. ^ a b c d Landau, L.D., Lifshitz, E.M. (2013). Fluid mechanics: Landau And Lifshitz: course of theoretical physics, Volume 6 (Vol. 6). Elsevier. Section 114, page 436.
  7. ^ Anderson, J. D. (). Modern compressible flow. McGraw-Hill. pp. 358–359. ISBN 0-07-242443-5. 
  8. ^ 1942, Landau, L.D. "On shock waves" J. Phys. USSR 6 229-230
  9. ^ Thompson, P. A. (1971). A fundamental derivative in gasdynamics. The Physics of Fluids, 14(9), 1843-1849.
  10. ^ Lamb (1994) §287, pp. 492–495.
  11. ^ Feynman, R. P.; Leighton, R. B.; Sands, M. (), The Feynman Lectures on Physics, 2, Addison-Wesley , p. 40-3. Chapter 40 has the title: The flow of dry water.
  12. ^ Batchelor (1973) pp. 404–405.
  13. ^ a b c d e f g h i Batchelor (1973) pp. 106–108.
  14. ^ a b c Batchelor (1973) pp. 409–413.
  15. ^ Kyrala, A. (). Applied Functions of a Complex Variable. Wiley-Interscience. pp. 116–117. ISBN 9780471511298. 

Bibliografie

modificare

Lectură suplimentară

modificare

Vezi și

modificare

Legături externe

modificare