|
Szoláris HD és MHD
A hidrodinamika mozgásegyenletei - az Euler- és a Navier-Stokes-egyenlet
A klasszikus mechanika keretein belül az asztrofizikai, illetve a szoláris hidrodinamika általános egyenletei három alapelvre épülnek: első (i) a
tömegmegmaradás, második (ii) az impulzusmegmaradás, és végül (iii) az energiamegmaradás. A következőkben röviden áttekintjük e három elv matematikai
formulákban történő megfogalmazását.
Az impulzus megmaradásának elve a dinamika alaptörvényeként, azaz Newton II. törvényeként szerepel a fizikában. Ezt az axiómát is Leonhard Euler (1707-1783)
alkalmazta először mozgó, áramló gázokra és folyadékokra, felfedezve azt a tényt, hogy az egyenletekben az ismeretlen függvények és deriváltfüggvényeik szorzata is fellép,
azaz matematikai szakkifejezést használva: az egyenletek nemlineárisak. Ennek tulajdonságnak a káoszelmélet felfedezésekor sikerült kimutatni
jelentőségét az 1960-as években.
A folyadékok mozgásegyenletének továbbfejlesztése Claude Navier (1785-1836) francia mérnök nevéhez fűződik, végül ma is használt alakját,
a napjainkban Navier-Stokes-egyenletként emlegetett alakot George Stokes (1819-1903) ír fizikus adta meg a XIX. század második felében.
Tehát a második alapelv az impulzus megmaradását fejezi ki. Newton II. axiómája szerint a
tömeg és a gyorsulás szorzata a v térfogatban lévő anyagra ható erők eredőjével egyenlő, azaz a térfogati és felületi erők összegével.
(i) A térfogati erőknek nevezzük azokat az erőket, amelyek a test minden részére hatnak (mint pl. a nehézségi erő), és így arányosak annak térfogatával
vagy tömegével (amennyiben megfelelően kis tartományra szorítkozunk). Ha elhanyagoljuk a csillag mágneses terét, akkor a tömegegységre jutó térfogati erő
csak a nehézségi erőből származik, azaz
ahol V(x,t) a gravitációs potenciál, melyet az alábbi módon számíthatunk ki:
ahol G a gravitációs konstans.
(ii) A felületi erők viszont rövid hatótávolságuk miatt csak a tekintett anyagrész
felülete mentén hatnak, és ezért a felülettel arányosak. A térfogatelem felülete mentén ható erők a szomszédos anyagrészektől származnak, pontosabban
azon molekuláktól, amelyek a térfogatelemet határoló felület közvetlen környezetében vannak. Mivel ezek az ún. molekuláris erők nagyon
rövid hatótávolságúak, ezért csak a felület közelében érzékelhetők, és emiatt nagyságuk a felülettel arányos. Ha a v(t) térfogatot határoló
felület s(t), akkor jelölje a t vektor az n normálisú egységnyi felületre
ható erőt, amelyet feszültségvektornak vagy röviden feszültségnek nevezünk. Ez általában nem merőleges a felületre.
Tehát t feszültség értéke függ a helytől, az időtől, illetve a dS felületelem helyzetétől:
A deformálható testek mechanikájából ismeretes a feszültségvektor tenzorral kifejezett megadása, azaz
mely összefüggés szerint a t feszültségvektor az n egységvektornak
lineáris vektorfüggvénye, ezért a Tki kilenc mennyiség egy
kétindexes tenzort alkot. Ezt a tenzort feszültségtenzornak nevezzük. A feszültségtenzor fizikai jelentése a következő: az azonos indexű, tehát a főátlóban álló elemek
a felület normálisa irányába mutató feszültségek, a vegyes indexűek pedig a felület érintője mentén hatnak. Ezért az azonos indexű elemeket húzó- vagy
nyomófeszültségeknek, a vegyes indexűeket nyíró- vagy csúsztató feszültségeknek nevezzük.
Felhasználva a térfogati és a felületi erők egyenleteit az impulzusmegmaradást kifejező egyenletet a következő alakban írhatjuk fel:
Ha felhasználjuk a transzport-egyenlet kontinuitási egyenlet segítségével levezetett alakját, akkor a fenti egyenlet az alábbi formában adható meg:
A felületre vonatkozó integrált a Gauss-Osztrogradovszkij-tétel alkalmazásával térfogati integrállá alakíthatjuk:
amivel
Ez az egyenlőség tetszőleges v(t) térfogat esetén is igaz, így fenáll
összefüggés, melyet először Augustin-Louis Cauchy (1789-1857) vezetett le ebben az általános alakban. Ez tehát a plazma
mozgásegyenlete általános alakban. Ez az egyenlet a bennük szereplő fizikai mennyiségeknek a test egyes pontjaiban felvett
értékei között teremt kapcsolatot. A mozgásegyenletben a felületi erők járuléka formálisan úgy jelenik meg, mintha az is
a térfogati erő sűrűsége lenne, hiszen a Tki feszültségtenzor divergenciájának a tér egyes pontjaiban felvett
értékei fordulnak elő a jobb oldalon. Ez azonban csak formálisan van így, mert ennek csak azokon a helyeken van zérustól különböző
járuléka, ahol a feszültségtenzor változik a hellyel.
Az ideális folyadék mozgásegyenlete - Euler-egyenlet |
|
A folyadékok - a tapasztalat szerint - térfogatváltozással nem járó alakváltozással szemben csak igen kicsi ellenállást mutatnak, ha
a deformáció lassan megy végbe. A szomszédos folyadékrétegek nagyon kis erővel elcsúztathatók egymással, mely azt jelenti, hogyha a
deformáció lassú, akkor a nyírófeszültségek elhanyagolhatóan kicsik, azaz pl. derékszögű koordinátákat használva
Az olyan idealizált folyadékot, melynél a mozgás folyamán
nem lépnek fel nyírófeszültségek súrlódásmentes vagy ideális folyadéknak nevezzük.
Vizsgáljuk meg a feszültségtenzor azonos indexű elemeit, melyeket húzó- vagy nyomófeszültségeknek nevezünk. Szintén derékszögű koordinátákat
használva a Txx az x normálisú lap felületegységére ható x irányú erőt jelenti, hasonlóan a másik két
Tyy és Tzz azonos indexű elemnél az y, illetve a z normálisú lap esetén.
A folyadékok elméletében a felületegységre merőlegesen ható erőt nyomásnak nevezzük, és akkor tekintjük pozitívnak, ha a
tekintett térfogatelembe befelé mutat. Ennek megfelelően
azaz
ahol δik=1 ha i=k és δik=0 ha i≠k.
Ha a térfogati erők esetén is behelyettesítjük a f=-gradV öszefüggést, akkor az ideális
folyadék mozgásegyenletére az alábbi összefüggést kajuk
melyet Euler vezetett le, és tiszteletére Euler-egyenletnek hívunk.
|
A súrlódó folyadék mozgásegyenlete - Navier-Stokes-egyenlet |
|
A valódi folyadékok mozgása során - ellentétben az ideális folyadékokkal - a súrlódás következtében nyírófeszültségek
is fellépnek a folyadékrétegek között. Erre mutat egyébként sok tapasztalat is, amely szerint az ideális folyadékok
dinamikáját leíró elméletek egyes eredményei ellentmondásban vannak a megfigyelésekkel. Az ideális folyadékok elméletének továbbfejlesztése
természetes módon adódik úgy, hogy az ideális folyadék feszültségtenzorát kiegészítjük egy taggal,
amely a súrlódási erőtől származó feszültségeket tartalmazza, azaz
ahol τik az ún. súrlódási vagy viszkózus tenzor. Újra a deformálható testek leírásához visszatérve, tudjuk, hogy a
feszültségtenzornak szimmetrikusnak kell lennie, ezért τik súrlódási tenzor is szimmetrikus. Mivel a
nyugvó folyadékokban súrlódás nem lép fel, a viszkózus tenzor a folyadék sebességterével van kapcsolatban oly módon, hogy
nyugalmi állapotban eltűnik. Továbbá a folyadék transzlációs és forgó mozgásánál az egyes rétegek egymáshoz képest nem
mozdulnak el, így nyírófeszültségek nem lépnek fel, ezért a súrlódási tenzor csak a deformációs sebességektől függ.
George Stokes (1819-1903) ír fizikus gondolatmenetéhez hasonlóan a tapasztalattal jó egyezésben felvethető, hogy
lineáris közelítésben az izotrop folyadékokra a τik súrlódási tenzor arányos a
Dik deformációs tenzorral:
ahol a deformációs tenzor tehát
és μ valamint μϑ a belső súrlódásra jellemző anyagi együtthatók. Ezek a folyadék
anyagi minőségétől és a hőmérséklettől függnek.
Behelyettesítve korábbi eredményeinket a mozgásegyenlet Cauchy által levezetett általános alakjában, a következőt kapjuk
ahol
A fenti vektoregyenletet Navier-Stokes-egyenletnek nevezük. Ezen egyenletet tehát a súrlódó folyadékok mozgásegyenlete. A Navier-Stokes-egyenletek
természetesen a μ=μϑ=0 esetben átmennek az ideális folyadékok mozgását leíró Euler-egyenletbe.
|
Mind az Euler- mind a Navier-Stokes-egyenlet tehét a folyadék áramlását határozza meg ideális, illetve súrlódó folyadék esetén.
Azonban ebben az alakban még túl általánosak ahhoz,
hogy belőlük az elmozdulásvektort a térkoordináták és az idő függvényében meghatározzuk, ugyanis a folyadékok áramlását leíró
mozgásegyenletek parciális differenciálegyenletek. A konkrét fizikai feladatok megoldásához ezeket az egyenleteket ki kell egészíteni az
ún. kezdeti és határfeltételekkel.
|