Curgere potențială

curgere în care vitezele sunt gradientul unei funcții scalare

În dinamica fluidelor, curgerea potențială sau curgerea irotațională se referă la o descriere a unei curgeri a unui fluid fără vorticitate. O astfel de descriere apare de obicei în limita viscozității care dispare, adică pentru un fluid ideal și fără vorticitate în curgere.

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

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

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

Printre aplicațiile curgerii potențiale se numără: curgerea exterioară pentru profile aerodinamice, undele de apă, curgerea electroosmotică și curgerea apelor subterane. În cazul curgerilor (sau al unor părți ale acestora) cu efecte puternice de vorticitate, aproximația curgerii potențiale nu este aplicabilă. În regiunile de curgere în care se știe că vorticitatea este importantă, cum ar fi urmele valurilor și straturile limită, teoria curgerii potențiale nu este capabilă să ofere predicții rezonabile ale curgerii.[1] 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: curgerea în jurul aeronavelor, curgerea apelor subterane, acustica, undele de apă și curgerea electroosmotică.[2]

Descriere și caracteristici

modificare
 
O curgere potențială se construiește prin adăugarea unor curgeri elementare simple și observarea rezultatului.
 
Linii de curgere pentru curgerea potențială incompresibilă în jurul unui cilindru circular, într-o curgere uniformă.

În curgerea potențială sau irotațională, câmpul vectorial de vorticitate este zero. Acest lucru poate fi exprimat astfel:

 

unde   este câmpul de viteze și   este câmpul de vorticitate. Deoarece orice câmp vectorial cu curbură zero poate fi exprimat ca gradient al unui scalar, câmpul de viteze poate fi scris ca gradient al unui scalar numit potențial de viteză,  . Deoarece curbura gradientului este întotdeauna zero, avem:[3]

 

Potențialul de viteză nu este definit în mod unic, deoarece la acesta se poate adăuga o funcție arbitrară a timpului, de exemplu  , sau un termen de forma   cu   fiind un câmp vectorial arbitrar, fără a afecta mărimea fizică relevantă, care este  . Neunicitatea este eliminată de obicei prin selectarea adecvată a condițiilor inițiale sau la limită satisfăcute de  , iar ca atare procedura poate varia de la o problemă la alta.

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

 

unde   este elementul de linie pe contur și   este elementul de suprafață al oricărei suprafețe delimitate de contur. În spațiul multiconectat (de exemplu, în jurul unui contur care înconjoară un corp solid în două dimensiuni sau în jurul unui contur care înconjoară un tor în trei dimensiuni) sau în prezența vârtejurilor concentrate (de exemplu, în așa-numitele vârtejuri irotaționale sau vârtejuri punctuale, sau în inele de fum), circulația   nu trebuie să fie zero. În primul caz, teorema Stokes nu poate fi aplicată, iar în ultimul caz,   este diferită de zero în regiunea delimitată de contur. În jurul unui contur care înconjoară un cilindru solid infinit de lung cu care conturul se buclează de   ori, avem:

 

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

Curgere incompresibilă

modificare

În cazul unei curgeri incompresibile, cum ar fi un lichid sau un gaz la numere Mach scăzute (dar nu pentru undele sonore), viteza v are divergență zero:[3]

 

Înlocuind aici  , observăm că   satisface ecuația lui Laplace[3]

 

unde 2 = ∇ ⋅ ∇ este operatorul Laplace (uneori scris și Δ). Deoarece soluțiile ecuației 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 complet determinat de cinematica sa: ipotezele de irotaționalitate și divergență zero a curgerii. Dinamica, în legătură cu ecuațiile de impuls, trebuie aplicată ulterior numai dacă se dorește calcularea câmpului de presiune: de exemplu, pentru curgerea în jurul profilurilor aerodinamice, prin utilizarea principiului lui Bernoulli.

În cazul curgerilor incompresibile, contrar unei concepții 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 zero. Incapacitatea curgerii potențiale de a satisface condițiile de limită necesare, în special în apropierea limitelor solide, face ca aceasta să nu fie valabilă în reprezentarea câmpului de curgere necesar. Dacă curgerea potențială satisface condițiile necesare, atunci aceasta reprezintă soluția necesară a ecuațiilor Navier–Stokes incompresibile.

În două dimensiuni, curgerea potențială incompresibilă se reduce la un sistem foarte simplu care este analizat folosind analiza complexă (a se vedea mai jos).

Curgere compresibilă

modificare

Curgere constantă

modificare

Teoria curgerii potențiale poate fi, de asemenea, utilizată pentru a modela curgerea compresibilă irotațională. Derivarea ecuației de guvernare pentru   din ecuația lui Euler este destul de simplă. Ecuațiile de continuitate și de impuls (de curgere potențială) pentru curgerile 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   din cele două ecuații de guvernare rezultă în:

 

Versiunea incompresibilă apare în limita  . Substituind aici   rezultă în[4][5]

 

unde   este exprimată ca o funcție de mărimea vitezei  . Pentru un gaz politropic,  , unde   este raportul de căldură specifică și   este entalpia de stagnare. În două dimensiuni, ecuația se simplifică la:

 

Validitate: În forma actuală, ecuația este valabilă pentru orice curgere potențială neviscidă, indiferent dacă curgerea este subsonică sau supersonică (de exemplu, curgerea Prandtl-Meyer). Cu toate acestea, în cazul curgerilor supersonice și, de asemenea, în cazul curgerilor transonice, pot apărea unde de șoc care pot introduce entropie și vorticitate în curgere, făcând curgerea să devină rotațională. Cu toate acestea, există două cazuri în care curgerea potențială prevalează chiar și în prezența undelor de șoc, care se explică prin ecuația impulsului (nu neapărat potențială) scrisă sub următoarea formă:

 

unde   este entalpia specifică,   este câmpul de vorticitate,   este temperatura și   este entropia specifică. Deoarece în fața undei de șoc conducătoare avem o curgere potențială, ecuația lui Bernoulli arată că   este constantă, fiind de asemenea constantă de-a lungul undei de șoc (condițiile Rankine-Hugoniot) și, prin urmare, putem scrie[4]

 

1) Atunci 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 vorticitate este zero. Undele de șoc la marginea anterioară ascuțită a unui con bidimensional sau a unui con tridimensional (curgere Taylor-Maccoll) au o intensitate constantă. 2) Pentru undele de șoc slabe, saltul de entropie peste unda de șoc este o mărime de ordinul trei în ceea ce privește 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: Atunci când curgerea este predominant unidirecțională, cu mici abateri, cum ar fi în cazul curgerii pe lângă corpuri subțiri, ecuația completă poate fi simplificată în continuare. Fie   curentul principal și să se ia în considerare mici abateri de la acest câmp de viteze. Potențialul de viteză corespunzător poate fi scris sub forma   unde   caracterizează mica abatere 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. Atunci când   este mică (curgere transonică), avem următoarea ecuație neliniară:[4]

 

Unde   este valoarea critică a derivatei Landau  [6][7] și   este volumul specific. Curgerea transonică este complet caracterizată de un singur parametru  , care pentru gazul politropic ia valoarea   . Sub transformarea hodografică, ecuația transonică în două dimensiuni devine ecuația Euler–Tricomi.

Curgere instabilă

modificare

Ecuațiile de continuitate și de impuls (de curgere potențială) pentru curgerile instabile sunt date de:

 

Prima integrală a ecuației impulsului (de curgere potențială) este dată de:

 

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

Înlocuind aici   rezultă în:

 

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

 

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

 

Unde sonore: În undele sonore, mărimea vitezei   (sau numărul Mach) este foarte mică, deși termenul instabil 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[8][4]

 

care este o ecuație de undă liniară pentru potențialul de viteză φ. Din nou, partea oscilatorie a vectorului viteză v este legată de potențialul viteză prin v = ∇φ, în timp ce, ca și înainte, Δ este operatorul Laplace, iar c este viteza medie a sunetului în mediul omogen. Rețineți că și părțile oscilatoare ale presiunii p și densității ρ satisfac fiecare în parte ecuația de undă, în această aproximație.

Aplicabilitate și limitări

modificare

Curgerea potențială nu include toate caracteristicile curgerilor care se întâlnesc în lumea reală. Teoria curgerii potențiale nu poate fi aplicată în cazul curgerilor interne vâscoase,[1] cu excepția curgerilor între plăci foarte apropiate. Richard Feynman a considerat că curgerea potențială este atât de nefizică, încât singurul fluid care se supune ipotezelor este „apa uscată” (citându-l pe John von Neumann).[9] De asemenea, curgerea potențială incompresibilă face o serie de predicții nevalabile, cum ar fi paradoxul lui d'Alembert, care afirmă că rezistența la înaintare a oricărui obiect care se deplasează printr-un fluid infinit, altfel în repaus, este egală cu zero.[10] Mai exact, curgerea potențială nu poate explica comportamentul curgerilor care includ un strat limită.[1] 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), cum ar fi vortexul liber și sursa punctiformă, dispun de soluții analitice rapide. Aceste soluții pot fi suprapuse pentru a crea curgeri mai complexe care să satisfacă o varietate de condiții la limită. Aceste curgeri corespund îndeaproape curgerilor din viața reală în întreaga mecanică a fluidelor. În plus, multe informații valoroase rezultă atunci când se ia în considerare abaterea (adesea ușoară) dintre o curgere observată și curgerea potențială corespunzătoare. Curgerea potențială își are numeroase 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ă în interiorul stratului limită. Absența efectelor stratului limită înseamnă că orice linie de curgere poate fi înlocuită cu o limită solidă fără nicio modificare a câmpului 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ă în două dimensiuni este simplă de analizat cu ajutorul cartografierii conforme, prin utilizarea transformărilor planului complex. Cu toate acestea, nu este necesară utilizarea numerelor complexe, așa cum se întâmplă, de exemplu, în analiza clasică a curgerii fluidelor pe lângă un cilindru. Nu este posibil să se rezolve o curgere potențială folosind numere complexe în trei dimensiuni.[11]

Ideea de bază este de a utiliza o funcție holomorfă (numită și analitică) sau meromorfă f, care face legătura între domeniul fizic (x, y) și domeniul transformat (φ, ψ). În timp ce x, y, φ și ψ au toate valori reale, este convenabil să se definească mărimile complexe:

 

Acum, dacă scriem corespondența f ca[11]

 

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

 

Componentele vitezei (u, v), respectiv în direcțiile (x, y), pot fi obținute direct din f prin diferențierea față de z. Adică[11]

 

Deci câmpul de viteză v = (u, v) este specificat de[11]

 

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

 

Astfel, φ poate fi identificat ca fiind potențialul de viteză, iar ψ se numește funcția de curgere.[11] Liniile de constantă ψ sunt cunoscute sub numele de linii de curent, iar liniile de constantă φ sunt cunoscute sub numele de linii echipotențiale (a se vedea suprafața echipotențială).

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

 

Astfel, curgerea are loc de-a lungul liniilor de constantă ψ și în unghi drept cu liniile de constantă φ.[11]

Δψ = 0 este, de asemenea, satisfăcută. Această relație este echivalentă cu rotația vectorului viteză v fiind egală cu zero: ∇ × v = 0. Așadar, curgerea este irotațională. Condiția automată 2Ψ/xy = 2Ψ/yx rezultă din constrângerea de incompresibilitate ∇ · v = 0.

Exemple de curgeri incompresibile bidimensionale

modificare

Orice funcție diferențiabilă poate fi utilizată pentru f. Exemplele de mai jos utilizează o varietate de funcții elementare; pot fi utilizate și funcții speciale. Rețineți că pot fi utilizate funcții cu mai multe valori, cum ar fi logaritmul natural, dar atenția trebuie limitată la o singură suprafață Riemann.

Legile puterii

modificare
Exemple de hărți conforme pentru legea de putere w = Azn
Exemple de hărți conforme pentru legea de putere w = Azn, pentru diferite valori ale puterii n. Este prezentat planul z, care arată liniile de potențial constant φ și de funcție de curgere ψ, în timp ce w = φ + .

În cazul în care se aplică următoarea hartă conformă a legii puterii, de la z = x + iy la w = φ + :[12]

 

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

 

În figurile din dreapta sunt prezentate exemple pentru mai multe 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 linii echipotențiale. Câteva valori interesante ale lui n sunt:[12]

  • n = 1/2: aceasta corespunde curgerii î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: curgerea datorată unei surse dublete.

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

Legi de putere cu n = −1: curgere uniformă

modificare

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

 

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

Legi de putere cu n = −2

modificare

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

 

care este un sistem de hiperbole dreptunghiulare. 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 viteză este dat de φ, sau

 

În dinamica fluidelor, câmpul de curgere din apropierea originii corespunde unui punct de stagnare. Rețineți că fluidul de la origine este în repaus (acest lucru rezultă din diferențierea lui f(z) = z2 la z = 0). Linia de curgere ψ = 0 este deosebit de interesantă: are două (sau patru) ramificații, urmând axele de coordonate, adică x = 0 și y = 0. Deoarece niciun fluid nu curge de-a lungul axei x, aceasta (axa x) poate fi tratată ca o frontieră solidă. Astfel, este posibil să se ignore curgerea în semiplanul inferior, unde y < 0, și să se concentreze asupra curgerii în semiplanul superior. Cu această interpretare, curgerea este cea a unui jet dirijat 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.

Legi de putere cu n = 3

modificare

Dacă n = 3, curgerea rezultată este o versiune hexagonală a cazului n = 2 prezentat anterior. Liniile de curent sunt definite de ecuația ψ = 3x2yy3. Curgerea poate fi interpretată ca o curgere într-un colț de 60°.

Legi de putere cu n = −1: dublet

modificare

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

 

Acest lucru este mai ușor de interpretat în termeni de componente reale și imaginare:

 

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

 

sau în coordonate polare:

 

Legi de putere cu n = −2: cvadrupol

modificare

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

 

Acesta este câmpul de curgere asociat cu un cvadrupol.[13]

Linie sursă și chiuvetă

modificare

O sursă de linie sau o chiuvetă de putere   (   pentru sursă și   pentru chiuvetă) este dată de potențialul:

 

Unde   de fapt este curgerea de volum pe unitate de lungime pe o suprafață care înconjoară sursa sau chiuveta. Câmpul de viteză în coordonate polare este dat de:

 

adică o curgere pur radială.

Linie vortex

modificare

Un vârtej de forță   este dat de

 

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

 

adică o curgere pur azimutală.

Analiză pentru curgeri incompresibile tridimensionale

modificare

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

Sursă punctiformă și chiuvetă

modificare

Potențialul de viteză al unei surse punctiforme sau al unei chiuvete de putere   (  pentru sursă și   pentru chiuvetă) în coordonate polare sferice este dat de:

 

unde   este curgerea de volum pe o suprafață închisă care înconjoară sursa sau chiuveta. Câmpul de viteză în coordonate polare sferice este dat de:

 
  1. ^ a b c Batchelor (1973) pp. 378–380.
  2. ^ Kirby, B.J. (), Micro- and Nanoscale Fluid Mechanics: Transport in Microfluidic Devices., Cambridge University Press, ISBN 978-0-521-11903-0 
  3. ^ a b c Batchelor (1973) pp. 99–101.
  4. ^ a b c d e 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.
  5. ^ Anderson, J. D. (). Modern compressible flow. McGraw-Hill. pp. 358–359. ISBN 0-07-242443-5. 
  6. ^ 1942, Landau, L.D. "On shock waves" J. Phys. USSR 6 229-230
  7. ^ Thompson, P. A. (1971). A fundamental derivative in gasdynamics. The Physics of Fluids, 14(9), 1843-1849.
  8. ^ Lamb (1994) §287, pp. 492–495.
  9. ^ 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.
  10. ^ Batchelor (1973) pp. 404–405.
  11. ^ a b c d e f g h i Batchelor (1973) pp. 106–108.
  12. ^ a b c Batchelor (1973) pp. 409–413.
  13. ^ 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