Curgere potențială
Î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]
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Î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 viteze 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ărei curbe simple conexe este zero. Acest lucru poate fi demonstrat folosind teorema lui Stokes,
unde este elementul de linie pe curbă și este elementul de suprafață al oricărei suprafețe delimitate de curbă. În spațiul multiplu conectat (să zicem, în jurul unei curbe care înconjoară un corp solid bidimensional sau în jurul unei curbe 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 curbă. În jurul unei curbe care înconjoară un cilindru solid infinit cu care curba 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ă
modificareCurgere staționară
modificareTeoria 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ă
modificareEcuaț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
modificareCurgerea 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 stratului 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ă
modificareCurgerea 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Ψ∂x ∂y = ∂2Ψ∂y ∂x ne oferă atunci constrângerea de incompresibilitate ∇ · v = 0.
Exemple de curgeri incompresibile bidimensionale
modificareOrice 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În cazul în care se aplică următoarea transformare conformă cu o lege de putere, de la z = x + iy la w = φ + iψ:[14]
atunci, scriind z în coordonate polare ca z = x + iy = reiθ, 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 = 12: corespunde cu curgerea în jurul unei plăci semi-infinite,
- n = 23: 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ă
modificareDacă 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
modificareDacă 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 θ = yr și cos θ = xr 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
modificareDacă n = 3, curgerea rezultantă este o versiune hexagonală a cazului n = 2 analizat anterior. Liniile de curent sunt date de ψ = 3x2y − y3, iar curgerea în acest caz poate fi interpretată ca o curgere într-un colț de 60°.
Legile puterii cu n = −1: dublet
modificareDacă 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
modificareDacă n = −2, liniile de curent sunt date de:
Acest câmp de curgere este asociat cu un cvadripolar(d).[15]
Sursa și scurgerea liniară
modificareO 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
modificareUn 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
modificarePentru curgeri tridimensionale, potențialul complex nu poate fi obținut.
Sursă și scurgerea punctuală
modificarePotenț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:
Note
modificare- ^ „Fluide Ideale” (PDF).
- ^ a b c Batchelor (1973) pp. 378–380.
- ^ 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
- ^ a b c Batchelor (1973) pp. 99–101.
- ^ „Turbulenţa | PDF” (în engleză). Scribd. Accesat în .
- ^ 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.
- ^ Anderson, J. D. (). Modern compressible flow. McGraw-Hill. pp. 358–359. ISBN 0-07-242443-5.
- ^ 1942, Landau, L.D. "On shock waves" J. Phys. USSR 6 229-230
- ^ Thompson, P. A. (1971). A fundamental derivative in gasdynamics. The Physics of Fluids, 14(9), 1843-1849.
- ^ Lamb (1994) §287, pp. 492–495.
- ^ 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.
- ^ Batchelor (1973) pp. 404–405.
- ^ a b c d e f g h i Batchelor (1973) pp. 106–108.
- ^ a b c Batchelor (1973) pp. 409–413.
- ^ Kyrala, A. (). Applied Functions of a Complex Variable. Wiley-Interscience. pp. 116–117. ISBN 9780471511298.
Bibliografie
modificare- Batchelor, G.K. (), An introduction to fluid dynamics, Cambridge University Press, ISBN 0-521-09817-3
- Chanson, H. (), Applied Hydrodynamics: An Introduction to Ideal and Real Fluid Flows, CRC Press, Taylor & Francis Group, Leiden, The Netherlands, 478 pages, ISBN 978-0-415-49271-3
- Lamb, H. () [1932], Hydrodynamics (ed. 6th), Cambridge University Press, ISBN 978-0-521-45868-9
- Milne-Thomson, L.M. () [1968], Theoretical hydrodynamics (ed. 5th), Dover, ISBN 0-486-68970-0
Lectură suplimentară
modificare- Chanson, H. (), „Le potentiel de vitesse pour les écoulements de fluides réels: la contribution de Joseph-Louis Lagrange [Velocity potential in real fluid flows: Joseph-Louis Lagrange's contribution]”, La Houille Blanche (în franceză), 93 (5), pp. 127–131, doi:10.1051/lhb:2007072
- Wehausen, J.V.; Laitone, E.V. (), „Surface waves”, În Flügge, S.; Truesdell, C., Encyclopedia of Physics, IX, Springer Verlag, pp. 446–778, arhivat din original la , accesat în
Vezi și
modificare- Curgerea potențială în jurul unui cilindru circular(d)
- Cod de curgere potențială aerodinamică(d)
- Transformare conformă
- Deriva Darwin(d)
- Rețea de curgere(d)
- Câmp laplacian
- Ecuația Laplace pentru curgerea irotațională(d)
- Teoria potențialului(d)
- Funcție de curent(d)
- Potențialul vitezei
- Descompunerea Helmholtz(d)
Legături externe
modificare- „Irrotational flow of an inviscid fluid”. University of Genoa(d), Faculty of Engineering. Accesat în .
- „Conformal Maps Gallery”. 3D-XplorMath. Accesat în . — Java applets for exploring conformal maps
- Potential Flow Visualizations - Interactive WebApps