Jatkuvuusyhtälö

Jatkuvuusyhtälö akselisymmetriselle virtaukselle

(a)

Tarkastellaan akselisymmetristä virtauskenttää, joka ilmaistaan lieriökoordinaatistossa (r, φ, z), jossa kaikki virtausmuuttujat eivät riipu atsimutaalisesta kulmasta φ – esimerkiksi aksiaalinen virtaus pyörivän kappaleen yli. Jos nopeuskomponentit (u, w) vastaavat koordinaattisuuntia (r, z), osoitetaan, että jatkuvuusyhtälö saadaan seuraavalla kaavalla:

∂u∂r+ur+∂w∂z=0

b)

Osoitetaan, että jatkuvuusyhtälö voidaan automaattisesti tyydyttää virtausfunktiolla ψ, jonka muoto on sellainen, että

u=1r∂ψ∂z,w=-1r∂ψ∂r

2.2

Jatkuvuusyhtälö kaksiulotteiselle virtaukselle polaarikoordinaatistossa

(a)

Harkitaan kaksiulotteista virtauskenttää, joka ilmaistaan lieriökoordinaatistossa (r, φ, z), jossa kaikki virtausmuuttujat ovat riippumattomia atsimutaalikulmasta φ – esimerkiksi virtaus ympyränmuotoisen sylinterin yli. Jos nopeuskomponentit (u, v) vastaavat koordinaattisuuntia (r, φ), osoita, että jatkuvuusyhtälö saadaan seuraavalla kaavalla:

∂u∂r+ur+1r∂v∂ϕ=0

b)

Osoita, että jatkuvuusyhtälö voidaan automaattisesti tyydyttää virtausfunktiolla ψ, jonka muoto on sellainen, että

u=1r∂ψ∂ϕ,v=-∂ψ∂r

2.3

Kuljetusyhtälö epäpuhtaudelle kaksiulotteisessa virtauskentässä

Monissa teknisissä sovelluksissa ollaan kiinnostuneita epäpuhtauden kulkeutumisesta nestevirrassa. Epäpuhtaus voi olla mitä tahansa saastuttavasta kemikaalista hiukkasiin. Hallitsevan yhtälön johtamiseksi on tunnustettava, että mikäli epäpuhtautta ei synny virtauskentässä, epäpuhtauden massa säilyy. Epäpuhtaudet voivat kulkeutua kahden eri fysikaalisen mekanismin, konvektion ja molekyylidiffuusion, avulla. Olkoon C epäpuhtauden pitoisuus (ts, massa nesteen tilavuusyksikköä kohti); tällöin epäpuhtauden kulkeutumisnopeus pinta-alayksikköä kohti saadaan kaavalla

-D∇C=-Di∂C∂x+j∂C∂y

missä i ja j ovat yksikkövektorit x- ja y-suunnissa ja D on diffuusiokerroin (yksikköinä m2/s, sama kuin kinemaattinen viskositeetti).

Huomaa, että diffuusio kuljettaa epäpuhtautta pitoisuusgradienttia alaspäin (ts, kuljetus tapahtuu suuremmasta pitoisuudesta pienempään pitoisuuteen); tästä johtuu miinusmerkki. Se on analoginen lämmön johtumisen kanssa.

(a)

Harkitaan äärettömän pieni suorakulmainen kontrollitilavuus. Oletetaan, että sen sisällä ei synny epäpuhtauksia ja että epäpuhtaus on riittävän laimea, jotta nestevirtaus säilyy muuttumattomana. Tarkastellaan kontrollitilavuuden massatasetta, osoitetaan, että epäpuhtauden kulkeutumisyhtälö kaksiulotteisessa virtauskentässä saadaan

∂C∂t+u∂C∂x+v∂C∂y-D∂2C∂x2+∂2C∂y2=0

(b)

Miksi on välttämätöntä olettaa epäpuhtauden laimea suspensio? Millaisen muodon kuljetusyhtälö saisi, jos tätä oletusta ei tehtäisi? Lopuksi, miten yhtälöä voitaisiin muuttaa siten, että siinä otettaisiin huomioon, että epäpuhtautta syntyy kemiallisessa reaktiossa nopeudella m˙c tilavuusyksikköä kohti.

2.4

Eulerin yhtälöt akselisymmetriselle virtaukselle

(a)

Virtauskentälle ja koordinaatistolle, joka on esitetty harjoituksessa 2.1, osoita, että Eulerin yhtälöt (invisidiset momenttiyhtälöt) ovat muodossa

ρ∂u∂t+u∂u∂r+w∂u∂z=ρgr-∂p∂rρ∂w∂t+u∂w∂r+w∂w∂z=ρgz-∂p∂z

2.5

Navier-Stokesin yhtälöt kaksiulotteiselle aksisymmetriselle virtaukselle

a)

Osoita, että harjoituksessa 2 kuvatun kaltaisen aksisymmetrisen viskoosivirtauksen muodonmuutosnopeudet ja pyörteisyys.1 kuvatun kaltaista virtausta, saadaan

ε˙rr=∂u∂r;ε˙zz=∂w∂z;ε˙ϕϕ=urγ˙rz=12∂w∂r+∂u∂z;η=∂w∂r-∂u∂z

Vihje: Huomaa, että atsimutaalinen venymisnopeus ei ole nolla. Helpoin tapa määrittää se on tunnistaa, että ε˙rr+ε˙φφ+ε˙zz=0 on vastattava jatkuvuusyhtälöä. (b)

Osoita, että Navier-Stokesin yhtälöt aksisymmetriselle virtaukselle saadaan

ρ∂u∂t+u∂u∂r+w∂u∂z=ρgr-∂p∂r+μ∂2u∂r2+1r∂u∂r-ur2+∂2u∂z2ρ∂w∂t+u∂w∂r+w∂w∂z=ρgz-∂p∂z+μ∂u∂r2w∂r2+1r∂u∂r+w∂u∂r+w∂2w∂z2

2.6

Eulerin yhtälöt kaksiulotteiselle virtaukselle polaarikoordinaatistossa

(a)

Harjoituksessa 2 kuvatulle kaksiulotteiselle virtaukselle.2, osoita, että Eulerin yhtälöt (inviscid momenttiyhtälöt) ovat muotoa

ρ∂u∂t+u∂u∂r+vr∂u∂ϕ-v2r=ρgr-∂p∂rρ∂v∂t+u∂v∂r+vr∂v∂r+vr∂v∂v∂ϕ-uvr=ρgphi-1r∂p∂p∂ϕ↪Ll_9553> Vihjeitä: (i) Momenttikomponenteilla, jotka ovat kohtisuorassa elementtisäätötilavuuden sivupintoihin nähden ja jotka tulevat ja lähtevät niistä, on pienet komponentit radiaalisuunnassa, jotka on otettava huomioon; samoin (ii) näihin pintoihin vaikuttavilla painevoimilla on pienet radiaalikomponentit. 2.7

On osoitettava, että rasitusnopeudet ja pyörteisyys harjoituksen 2 virtaukselle ja koordinaatistolle.6 saadaan

ε˙rr=∂u∂r;ε˙φφ=1r∂v∂φ+urγ˙rφ=12∂v∂r-vr+1r∂u∂φ;ζ=1r∂u∂φ-∂v∂r+vr

Vihje: (i) Sivupinnan vääristymiskulma (β) on määriteltävä suhteessa suoraan, joka yhdistää origon O infinitesimaalisen kontrollitilavuuden keskipisteeseen. 2.8

Virtausta kapeassa raossa (jonka leveys on h) kahden pituudeltaan L:n pituisen konsentrisen sylinterin välissä, joista sisempi, jonka säde on R, pyörii kulmanopeudella ω, voidaan approksimoida Navier-Stokesin yhtälöiden Couette-ratkaisulla. Osoita, että akselin pyörittämiseen pyörimisnopeudella ω rad/s tarvittava vääntömomentti T ja teho P saadaan

T=2πμωR3Lh,P=2πμω2R3Lh

2.9

Akselisymmetrinen pysähdyspistevirtaus

Toteuta samanlainen analyysi kuin jaksossa 2.10.3 käyttäen Navier-Stokesin yhtälöiden aksisymmetristä muotoa, joka on annettu harjoituksessa 2.5 aksisymmetriselle pysähtymispistevirtaukselle, ja osoita, että yhtälön (2.118) ekvivalentti on

.

Vastaa

Sähköpostiosoitettasi ei julkaista.