SISUKORDA
Molekulaarfüüsika: erisoojused ja ülekandenähtused Olekuvõrrandi mehhanistlik esitus Gaasi soojusmahtuvus Tasakaaluline kiiruste jaotus - Maxwelli jaotus Baromeetriline valem Boltzmanni jaotus Maxwelli jaotus Järeldused Maxwelli valemist Maxwelli jaotuse eksperimentaalsest kontrollist Ülekandenähtused


Loeng 9. Molekulaarfüüsika: erisoojused ja ülekandenähtused


Olekuvõrrandi mehhanistlik esitus. 1811. a. formuleeris itaalia füüsik Amedeo Avogadro moolruumalade võrdsuse seaduse, mida tänapäeval tunneme Avogadro seadusena:

Samadel füüsikalistel tingimustel on kõigi gaaside moolruumalad võrdsed.

Sisuliselt on see olekuvõrrandi üks võimalik sõnastus ja vaevalt oleks ta erilist tähelepanu esile kutsunud, kui mitte Avogadro ise poleks andnud talle huvitava tõlgenduse: et kõigi ainete gramm-molekulid sisaldavad võrdse arvu osakesi (molekule), siis sõltub molekulide ruumtihedus üksnes rõhust ja temperatuurist. Asendades olekuvõrrandis ainehulga $z$ molekulide arvuga $N=nV=zN_A$ (suurust $N_A$, mis väljendab molekulide arvu ühes moolis aines, nim. Avogadro arvuks), saame
Avogadro seadus on samaväärne Clapeyroni võrrandiga.

\begin{displaymath}pV=zRT={{nV}\over{N_A}}RT~~\Longrightarrow~~p=nkT,\end{displaymath}

kus $k=R/N_A=1.38\cdot 10^{-23}$J/K kannab nimetust Boltzmanni konstant.

Boltzmanni konstanti võib vaadelda kui "vaba molekuli universaalkonstanti".

Olekuvõrrandi kaks kuju erinevad selle poolest, et esimene neist (Clapeyroni-Mendelejevi võrrand) on puht-empiiriline, sidudes otseselt mõõdetavaid suurusi. Teine seevastu lähtub aine diskreetsuse hüpoteesist (kõik ained koosnevad molekulidest). Hüpoteesi kinnituseks pidi olema Avogadro arvu katseline määramine, mis õnnestus alles sadakond aastat hiljem (Jean Perrin, 1909, lähtudes Browni liikumisest). Küll aga võimaldas teine kuju üles ehitada termodünaamika mehhanistliku esituse - gaaside molekulaar-kineetilise teooria.

Molekulaar-kineetilise teooria põhivõrrand tuletatakse eeldusest, et gaasi rõhk tekib molekulide põrgetel vastu anuma seina (gaasi või vedeliku sisse paigutatud keha välispinda). Et impulssesituses avaldub Newtoni II seadus kujul $\Delta \vec
p=\vec F\Delta t$, võime rõhu avaldise kirjutada

\begin{displaymath}p={F\over S}={{\sum\Delta \vec p}\over{S\Delta t}},\end{displaymath}

kus $\sum\Delta\vec p$ on kõigi ajavahemiku $\Delta t$ jooksul vastu pinnatükki $S$ põrganud molekulide summaarne impulsimuutus. (Õnnetuseks tähistatakse nii impulssi kui rõhku $p$-ga, seetõttu eristab neid meie valemites vaid vektorimärk.) Impulsi muutust on lihtne arvutada, kui teame impulssvektori suunda ja suurust enne põrget; paraku pole mingit lootust seda kõigi molekulide kohta teada saada.

Siin tulebki mängu statistika. Me teeme oletuse, et gaas koosneb molekulidest; kujutame neid punktmassidena ja rakendame neile Newtoni seadusi. Saadud tulemusi töötleme statistika valemite abil ning leiame meid huvitava suuruse kui üksikute molekulide parameetrite keskväärtuse.

Teeme näiteks oletuse, et kõik molekulid liiguvad võrdsete kiirustega ja pind $S$, millele mõjuvat rõhku arvutame, asub $yz$-tasandis. Põrgete arvu rehkendame molekulide ruumtiheduse $n$ (molekulide arv kuupmeetris) ja liikumiskiiruse $v$ abil. Kui molekulid liiguksid kõik $x$-telje ehk baasivektori $\vec i$ suunas, jõuaks meie pinnatükini ajavahemiku $\Delta t$ jooksul

$\Delta t$

Ühe molekuli põrkel seinale antud (jõu)impulss.
Kuna me ei tea põrke kestust, siis mingit "seinale mõjuvat jõudu" me sellest arvutada ei saa.

\begin{displaymath}N=nV=n(v_x\Delta t\cdot S)\end{displaymath}

molekuli. Kui lugeda põrge seinaga absoluutselt elastseks, hakkaksid molekulid pärast põrget liikuma vastassuunas ja sama kiirusega -- $\vec v_2=-\vec v_1$, seega muutuks nende impulss $\vec p_2-\vec
p_1=2m\vec v$ võrra. Pinnatükile $S$ mõjuva rõhu valemis asendame ${\sum\Delta \vec p}=N\cdot 2mv$ ja saame rõhu $p$ väärtuseks

\begin{displaymath}p={{nv\Delta tS\cdot 2mv}\over{S\Delta t}}=2nmv^2.\end{displaymath}

Paraku ei anna see valem meile otsitavat rõhku. Molekulaarfüüsika üheks eelduseks, mida kinnitab toosama Browni liikumine, on see, et molekulid liiguvad kaootiliselt, st kõikvõimalikes suundades ja kõikvõimalike kiirustega. Et leida õiget rõhku, tuleb meil arvutada keskväärtus.

Teeme asja pisut lihtsamaks, oletades, et gaas koosneb ühesugustest molekulidest. Ja nüüd statistika:

  1. Pinnatüki $S$ poole liiguvad pooled molekulid - jagame leitud rõhu kahega;
  2. Iga kiiruse ruut avaldub kolme komponendi $v_x,~v_y,~v_z$ ruutude summana. Meile läheb korda vaid $x$-komponent, mis on üks kolmandik kolme komponendi ruutude summast - jagame tulemust veel kolmega;
  3. Asendame molekulide kiiruste ruutude summa nende ruutkeskmise kiirusega

    \begin{displaymath}\bar v^2={1\over N}\sum v^2\end{displaymath}

    ja saame tulemuseks

    \begin{displaymath}p={1\over 3}nm\bar v^2.\end{displaymath}

$molekuli moju$
Molekulide "kollektiivne mõju".
Ajavahemiku $\Delta t$ jooksul põrkuvad vastu silindri põhja (pindala S) kõik selles silindris olevad molekulid
.

Gaasi rõhk on võrdeline molekulide ruutkeskmise kiirusega.

Et valemis on rõhk võrdeline molekulide ruumtihedusega $n$, saame kohe lahti veel ühest piirangust. Me ei pea enam jälgima, et molekulid oleks ühesugused. Kui meil on näiteks segu kolmest gaasist tihedustega $n_1,~n_2,~n_3$, saame kogurõhuks

\begin{displaymath}p=n_1m_1\bar v_1^2+n_2m_2\bar v_2^2+n_3m_3\bar v_3^2=p_1+p_2+p_3.\end{displaymath}

Jälle hea valem: gaaside segu rõhk on võrdne komponentide osarõhkude summaga. Pidage see meeles, võib tarvis minna.

Veelgi huvitavama seose saame, kui asendame $mv^2=2E_k$. Tuleb välja, et $p={2\over 3}nE_k$, ehk rõhk võrdub kahe kolmandikuga molekulide kineetilisest energiast ruumiühiku kohta. Tõepoolest -- ka dimensioon klapib: $J/m^3={{Nm}\over{m^3}}={N\over{m^2}}=Pa$. Asendame siia rõhu avaldise molekulide ruumtiheduse ja temperatuuri kaudu:

\begin{displaymath}p=nkT={2\over 3}nE_k~~~\Longrightarrow~~~E_k={3\over 2}kT\end{displaymath}

Boltzmanni konstant seob molekulide keskmise kineetilise energia gaasi temperatuuriga.

-- molekuli kineetiline energia on võrdeline temperatuuriga!

Seega on molekulaar-kineetilises teoorias rõhk, energia ja temperatuur ühe ja sellesama suuruse -- molekulide kiiruse -- avaldumisvormid. Siit ka teooria nimetus. Ilus teooria ja õige kah - niivõrd, kuivõrd Newtoni mehaanikat saab rakendada nii väikestele osakestele.


 
Gaasi soojusmahtuvus. Soojusmahtuvuse mõiste tuleb mängu, kui uurime kehade soojenemist. Et tõsta mingi keha temperatuuri ühe kraadi võrra, kulub kindel hulk soojust. See ongi soojusmahtuvus.

Soojuse mõiste tuli füüsikasse koos temperatuuriga. Sedagi tuli mõõta ning selleks defineeriti uus suurus -- soojushulk, mille mõõtmiseks oli vaja ka ühikut. Algselt defineeriti soojushulga ühik etaloonina: üks kalor ($cal$) oli soojushulk, mis tõstis ühe kuupsentimeetri vee temperatuuri ühe kraadi võrra. Hiljem, kui inglise füüsik James Joule 1843. a. tõestas töö ja soojuse ekvivalentsuse ($Q=c\Delta E$, $c=4.18J/cal$ on soojuse mehaaniline ekvivalent), hakati lihtsuse mõttes soojushulki mõõtma energiaühikutes. Nii toimitakse ka SI-süsteemis; "mehaanilist ekvivalenti" asendab mõiste "vee erisoojus": $c_{H_2O}=4180J/kg$.

Anname nüüd definitsioonid:

  1.   Soojusmahtuvus = soojushulk dzaulides, mis tõstab keha temperatuuri ühe kelvini võrra

  2.   Erisoojus = soojushulk, mis tõstab antud aine massiühiku (kilogrammi) temperatuuri 1 K võrra

  3.   Moolsoojus = soojushulk, mis tõstab antud aine ühe mooli ...

Soojushulk ja energia on kattuvad mõisted; mõlemat võib mõõta sma ühikuga.


Keha soojusmahtuvus näitab, kui palju soojust (energiat) on vaja tema temperatuuri tõstmiseks ühe kelvini võrra.


Erisoojus ja moolsoojus iseloomustavad mitte kehi, vaid aineid.

Kui vedelikel ja tahketel ainetel piisab ühest erisoojusest, siis gaaside soojusmahtuvus sõltub termodünaamilisest protsessist. Isohoorilisel protsessil (konstantsel ruumalal) kulub temperatuuri tõstmiseks vähem energiat, kui konstantse rõhu, seega muutuva ruumala korral. Seetõttu eristatakse kaht mool- ja erisoojust:
  • $C_v$ - isohooriline moolsoojus - kui soojendamine toimub konstantse ruumala juures;
  • $C_p$ - isobaariline moolsoojus - kui soojendamine toimub konstantsel rõhul.

Kui muutuvad nii rõhk kui ruumala, ei saa temperatuuri muutust moolsoojuste abil arvutada. Siis lähtutakse olekuvõrrandist ja protsessi valemitest.

Seose moolsoojuste väärtuste vahel leidis saksa arst Julius Mayer 1842. a., lähtudes eeldusest, et soojuste vahe kulutatakse paisumisel tehtud tööks. Et paisumise töö isobaarilisel juhul on $A=p\Delta V$ ($A=F\Delta s=pS\Delta s=p\Delta V$), saame ühe mooli gaasi jaoks olekuvõrrandist

\begin{displaymath}A=p\Delta V=R\Delta T,\end{displaymath}

millest järeldub

\begin{displaymath}C_p=C_V+R,\end{displaymath}

seos on tuntud Mayeri valemina.

Moolsoojused ise on erinevatel gaasidel üsna sarnased, sõltudes vaid molekuli ehitusest.
  • Kõigil üheaatomilistel gaasidel (inertgaasid näiteks) on $C_v={3\over 2}R$,
  • kaheaatomilistel ${5\over 2}R$,
  • kolme- ja enamaatomilistel $3R$.

    Kavala tõestuse on välja pakkunud molekulaarfüüsikud: et nende teoorias ühe aatomi energia on $E=E_k={3\over 2}kT$, saavad nad ühe mooli koguenergiaks $E={3\over 2}N_AkT={3\over 2}RT$, kust $\Delta
E={3\over 2}R\Delta T$, millest $C_v={3\over 2}R$ - nagu peabki olema üheaatomilistel gaasidel.

    Et kuidagi seletada mitmeaatomiliste gaaside moolsoojusi, tuldi mõttele sisse tuua energia jagunemine vabadusastmete peale. Üheaatomilise molekuli liikumisel on kolm vabadusastet (kiirusvektori kolm komponenti), mitmeaatomilistel lisandub veel pöörlev liikumine (nurkkiirusvektoril samuti kolm komponenti) - seega kokku kuus vabadusastet.

    Kuna energia salvestub kineetilise energiana ($E_k={{mv^2}\over 2}$ kulg- ja $E_k={{I\omega^2}\over 2}$ pöördliikumisel), sõltub salvestuva energia hulk lisaks kiiruse ruudule ka massist (inertsimomendist). Energia valemisse kuuluv mass on kõigi liikumissuundade jaoks ühesugune, inertsimoment seevastu sõltub molekuli ehitusest. Kaheaatomilisel molekulil, kui aatomid lugeda punktmassideks, on molekuli teljega (aatomeid ühendav sirglõik) ühtiva telje suhtes inertsimoment null ja seetõttu ei kõlba üks vabadusaste kuuest energia salvestamiseks - nii jääb vabadusastmeid mitte kuus, vaid viis.

    Moolsoojuse valemid võime seega anda kujul:

    \begin{displaymath}C_V={i\over 2}R;~~~~~C_p={{i+2}\over 2}R.\end{displaymath}

    - hea lihtne meeles pidada.

    Aga kuidas arvutada erisoojusi?

  • Gaaside soojusmahtuvuse teoorias eristatakse kaht moolsoojust: isobaarilist ja isohoorilist. Neid seob Mayeri valem.



    Molekuli kulgliikumise kineetilise energia seos vabadusastmetega: mv2=mvx2+mvy2+mvz2. Kui liita kokku kõigi molekulide energiad, tuleb iga vabadusastme kohta kolmandik summaarsest energiast.



    Kaheaatomilise molekuli pöörlemisenergia.
    Pöörlemine ümber z -telje energiat ei lisa, kuna inertsimoment on null.


     
    Tasakaaluline kiiruste jaotus -- Maxwelli jaotus. Kuna rõhk ja siseenergia avaldusid keskmise kiiruse kaudu, peaks molekulaarfüüsika seisukohalt olema lõpmatu arv kiiruste jaotusi, mis vastavad ühele ja samale olekule (samale temperatuurile). Molekulaarfüüsikas väidetakse, et kõigist neist võimalikest jaotustest on üks suurima tõenäosusega jaotus, mida nimetatakse tasakaaluliseks jaotuseks. Molekulide omavaheliste põrgete käigus toimub pidev energiavahetus, mille tõttu muutuvad molekulide kiirused ja lõppkokkuvõttes ka kiiruste jaotuse iseloom;

    tasakaaluline jaotus on selline, kus energiavahetused kiiruste jaotust ei muuda - mingist kiiruste grupist $\{v_1, v_1+\Delta v\}$ põrgete tagajärjel lahkuvate molekulide arv on võrdne sinna teistest gruppidest tulevate molekulide arvuga.

    Molekulaarfüüsika seisukohalt on olek tasakaaluline siis, kui molekulide kiiruste jaotus enam ei muutu. Seda viimast nimetataksegi tasakaaluliseks jaotuseks.


     
    Baromeetriline valem. Kõige lihtsam viis molekulide kiiruste jaotuse määramiseks on uurida nende tiheduse jaotust raskusjõu väljas. Siin on meil kasutada kaks valemit:

    • $p=nkT$   --   rõhu valem molekulaarfüüsikast,
    • $p=\rho gh=mngh$   --   rõhu sõltuvus sügavusest hüdrodünaamikas.

    Viimases valemis tuleb aga arvestada, et erinevalt vedelikest on gaas kokkusurutav. Seetõttu sõltub tihedus rõhust ning meie valemid on täpsed vaid lõpmata õhukeses kihis paksusega $dh$.
    Diferentseerime meie võrrandeid, saades

    \begin{displaymath}dp=kT\cdot dn\end{displaymath}

    \begin{displaymath}dp=-mngdh,\end{displaymath}

    (miinusmärk näitab, et kõrguse suurenedes rõhk kahaneb). Lahutame esimesest võrrandist teise, eraldame muutujad ning integreerime:

    \begin{displaymath}kTdn+mngdh=0~~\Longrightarrow~~{{dn}\over
n}=-{{mg}\over{kT}}dh\end{displaymath}

    \begin{displaymath}\ln n=-{{mg}\over{kT}}h+{\rm const}~~\Longrightarrow~~
n=n_0e^{-{{mgh}\over{kT}}},\end{displaymath}

    kus $n_0=n(h=0)$.

    Kui asendada selles valemis $n=p/kT$ ja $n_0=p_0/kT$, saame baromeetrilise valemi, mis kirjeldabki rõhu sõltuvust kõrgusest konstantse temperatuuri korral (isotermilises atmosfääris).


    Baromeetrilise valemi tuletus: igas õhukeses kihis paksusega dh loeme rõhu ja tiheduse konstantseks; rõhk mingil kõrgusel h sisaldab kõigi kõrgemal olevate kihtide kaalude summat (integraali).

    Küsimus: Kas suudate seda valemit siinsamas ja kohe kirja panna?


     

    Boltzmanni jaotus. Kui vaadata hoolega eksponendi $e$ astendajat baromeetrilises valemis, näeme, et lugejas olev suurus vastab massiga $m$ molekuli potentsiaalsele energiale raskusjõu väljas. Kui lasta sellel molekulil kõrguselt $h$ "alla kukkuda", läheb tema potentsiaalne energia üle kineetiliseks: $E_k(h=0)=E_p(h)$. Et molekulide jaotus atmosfääris peaks olema tasakaaluline, saame kiiruste jaotuse kõrgusel $h=0$, kui asendame baromeetrilises valemis potentsiaalse energia kineetilisega ($mgh={{mv^2}\over 2}$):

    Boltzmanni jaotuse saame, kui asendame baromeetrilises valemis potensiaalse energia mgh kineetilise energiaga mv2/2

    \begin{displaymath}n=n_0e^{-{{mv^2}\over{2kT}}}.\end{displaymath}

    Mida aga tähistab sel juhul suurus n? Baromeetrilises valemis näitas see molekulide arvu ruumalaühiku kohta kõrgusel $h$, kuid ei öelnud midagi nende molekulide kiiruste kohta. Võiksime väita, et see tihedus sisaldab kõiki neid molekule, mis võiksid tõusta kõrgemale kõrgusest $h$. Molekulide koguarv $n_0$ vastaks siis neile molekulidele, mis suudavad tõusta kõrgemale kõrgusest $h=0$.

    Kiiruste jaotuses väljendaks sel juhul arv $n(v)$ selliste molekulide ruumtihedust, mille kiirus on suurem kiirusest $v$. Konstant $n_0$ on nullist suurema kiirusega molekulide arv ruumalaühikus ehk teiste sõnadega - molekulide kogutihedus.

    Seda valemit nimetataksegi Boltzmanni jaotuseks. Mikrofüüsikas väljendab ta kõige sagedamini osakeste jaotust energiate järgi. Meie uurime praegu küll kiirusi - tegelikult on see aga ükspuha.

    Boltzmanni jaotus kuulub nn. integraalsete jaotusfunktsioonide hulka (statistikas nimetatakse neid ka kumulatiivseteks - ld. cumulatio, kuhjumine).


     
    Maxwelli jaotus seevastu on diferentsiaalne jaotusfunktsioon, mis väljendab mingi kiirusega osakeste suhtelist hulka. Et täpselt antud kiirusega liikuvate molekulide suhteline hulk on lõpmata väike, defineeritakse diferentsiaalne jaotus valemiga

    \begin{displaymath}f(v)dv={{n(v\in\{v,v+dv\})}\over{N=n(v\in\{-\infty,+\infty\})}}.\end{displaymath}

    Meie juhul, kus kiiruse moodul ei saa olla negatiivne, on nimetajas vahemik nullist lõpmatuseni. Integraalne jaotusfunktsioon saadakse diferentsiaalse jaotusfunktsiooni integreerimisel $v$-st lõpmatuseni; aga võib ka integreerida suvalises vahemikus, saades teatud kiiruste vahemikku kuuluvate molekulide suhtelise hulga.

    Näiteks peab Maxwelli jaotusest tulenema Boltzmanni jaotus:

    \begin{displaymath}n(u>v)=n_0\int\limits _v^\infty f(u)du=n_0e^{-{{mu^2}\over{2kT}}}.\end{displaymath}

    Diferentsiaalse jaotusfunktsiooni leidmiseks tuleb lahendada pöördülesanne: diferentseerides integraalset (Boltzmanni) jaotust jõuda diferentsiaalse (Maxwelli) jaotuseni.

    Lihtne tuletise võtmine annab

    \begin{displaymath}f(u)={{mu}\over{kT}}e^{-{{mu^2}\over{2kT}}},\end{displaymath}

    mis pole paraku õige, kuna ei arvesta kiiruste jaotust suundade järgi ega kogutiheduse sõltuvust temperatuurist.

    Et ülesandega hakkama saada, tuleb välja kannatada lühike geomeetriline ekskursioon faasiruumi. Kui kujutame kiiruste jaotust faasiruumis, saame kolmemõõtmelise sfääri, kus kiiruste vahemikule $\{u,u+du\}$ vastab kiht pindalaga $4\pi u^2$ ning paksusega $du$. Selle kihi ruumala vastabki vaadeldavasse vahemikku langevate kiirustega osakeste suhtelisele hulgale.

    Seega, kui kõik kiirused oleks võrdtõenäosed (faasiruum oleks kiiruste osas ühtlaselt täidetud), saaksime

    \begin{displaymath}f(v)={\rm const}\cdot v^2,\end{displaymath}

    mis on täiesti absurdne oletus. Suure kiirusega osakeste arv peab seetõttu kahanema kiiremini, kui $n\sim v^{-2}$, ja see saavutatakse eksponentteguri $e^{-kx}$ lisamisega valemisse.


    Kõrgusel h = 0 on kiirusest v suurema kiirusega molekule niisama palju, kui neid on kõrgemal kui h3= v32/2g .


    Boltzmanni jaotus kiiruste järgi.

    Integraalne (kumulatiivne) jaotusfunktsioon näitab, kui palju on antud energiast suurema energiaga molekule.

    Diferentsiaalne jaotusfunktsioon näitab, kui suur on mingisse kiiruste (kõrguste) vahemikku kuuluvate molekulide osakaal.

    Mehaanikas nimetatakse faasiruumiks 6-mõõtmelist ruumi, mille koordinaatideks on lisaks tavalistele ruumikoordinaatidele osakeste kiiruse või impulsi komponendid.

    Tulemuseks on laialtkasutatav valem

    \begin{displaymath}f(v)=c_1v^2e^{-c_2}v^n,\end{displaymath}

    kus $c_1$ ja $c_2$ on konstandid ning $n$ funktsiooni asümptoodile $v\to\infty$ vastav astendaja. Selle, astme- ja eksponentfunktsiooni korrutisena saadava kõvera kujunemist illustreerib kõrvaltoodud graafik.

    Meie ülesanne on konkreetsem: tuletada jaotusfunktsioon, mis on ruumis isotroopne ja vastab Boltzmanni jaotusele. Kui liikumine oleks ühemõõtmeline, võiksime joonisel kujutatud loogikast lähtudes kirjutada diferentsiaalse jaotusfunktsiooni tuletisena:

    \begin{displaymath}n(v\ge u)=\int\limits _u^\infty uf(u)du=n_0e^{-{{mu^2}\over{2...
...er{2kT}}}\right)=n_0{{2mu}
\over{2kT}}e^{-{{mu^2}\over{2kT}}}.\end{displaymath}

    \begin{displaymath}n(v\ge u)=\int\limits _u^\infty uf(u)du=n_0e^{-{{mu^2}\over{2...
...er{2kT}}}\right)=n_0{{2mu}
\over{2kT}}e^{-{{mu^2}\over{2kT}}}.\end{displaymath}

    Teguri $u$ lisamine integraali on vajalik, kuna ühemõõtmelisel juhul sõltub oskeste voog läbi "kontrollpinna" lisaks viimase pindalale ka osakeste kiirusest: sama tiheduse korral jõuab ajaühikus pinda läbida seda rohkem osakesi, mida kiiremini nad liiguvad. Et molekulide tihedus $n_0$ on praeguse probleemi seisukohalt ebaoluline, tähistame seda sisaldava konstandi lihtsalt $C$-ga ja leiame ta (pärast võrrandi mõlemate poolte $u$-ga läbi jagamist) normeerimistingimuse

    \begin{displaymath}\int\limits _{-\infty}^\infty f(u)du=1\end{displaymath}

    abil. Kui nii talitada, saame osakeste suhtelise tõenäosuse jaotuse, mis näitab, kui suure osa molekulide koguarvust moodustavad antud kiirusega (kiiruste vahemikus) antud suunas liikuvad osakesed.

    Maxwelli jaotus ja selle "piirajad" - funktsioonid v2 ning
    Integraalse jaotusfunktsiooni saame diferentsiaalse jaotusfunktsiooni integreerimisel.

    NB! Maxwelli jaotuse tuletamine on lisatekst ning ei kuulu kohustusliku (eksamil kontrollitava) materjali hulka. Jaotust ennast ning sellele vastavat kõverat tuleb tunda ning osata seletada.

    Kõik diferentsiaalsed jaotusfunktsioonid normeeritakse, st. jagatakse läbi vastava funktsiooni integraaliga üle kogu määramispiirkonna.

    Integreerime nüüd saadud funktsiooni. Selleks teeme muutuja vahetuse:

    \begin{displaymath}{{mu^2}\over{2kT}}=x^2~~\Longrightarrow~~{{2mudu}\over{2kT}}
...
...u^2}\over{kT}}dx~~
\Longrightarrow~~du=\sqrt{{2kT}\over m}dx,\end{displaymath}

    \begin{displaymath}{{mu^2}\over{2kT}}=x^2~~\Longrightarrow~~{{2mudu}\over{2kT}}
...
...u^2}\over{kT}}dx~~
\Longrightarrow~~du=\sqrt{{2kT}\over m}dx,\end{displaymath}

    saades integraali

    \begin{displaymath}\int\limits _{-\infty}^\infty f(u)du=C\sqrt{{2k}\over
m}\int...
...s _{-\infty}^\infty e^{-x^2}dx=C\sqrt{{2kT}\over
m}\sqrt\pi=1,\end{displaymath}

    kust

    \begin{displaymath}C={1\over{\sqrt\pi}}\sqrt{m\over{2kT}}=\sqrt{m\over{2\pi
kT}...
...m ja}~~f(u)du=\sqrt{m\over{2\pi
kT}}e^{-{{mu^2}\over{2kT}}}du.\end{displaymath}

    Faasiruum: iga osakest iseloomustab lisaks kohavektorile ka kiirusvektor . Kokku kuus koordinaati.

    See oleks siis tõenäosus, et osake liigub $x$-telje suunas kiirusega $u\leq v\leq u+du$. Tõenäosuse, et osake liigub mingis suvalises suunas kiirusega $\vec v=v_x\vec i+v_y\vec j+v_z\vec k$, saame, kui korrutame omavahel kolm tõenäosust:
    1. et $v_x$ on vahemikus $u_x\leq v_x\leq
u_x+du_x$ -- $=f(u_x)$;
    2. et $v_y$ on vahemikus $u_y\leq v_y\leq
u_y+du_y$ -- $=f(u_y)$;
    3. et $v_z$ on vahemikus $u_z\leq v_z\leq
u_z+du_z$ -- $=f(u_z)$.

    Seega

    \begin{displaymath}f(u)du=f(u_x)f(u_y)f(u_z)du_xdu_ydu_z=\left(\sqrt{m\over{2\pi...
...{-{{mu_y^2}\over{2kT}}}e^{-{{mu_z^2}\over{2kT}}}
du_xdu_ydu_z.\end{displaymath}

    \begin{displaymath}f(u)du=f(u_x)f(u_y)f(u_z)du_xdu_ydu_z=\left(\sqrt{m\over{2\pi...
...{-{{mu_y^2}\over{2kT}}}e^{-{{mu_z^2}\over{2kT}}}
du_xdu_ydu_z.\end{displaymath}

    Nüüd teeme kaks matemaatilist nippi:
  • esiteks viime eksponendid ühisele astendajale, saades

    \begin{displaymath}e^{-{{mu_x^2}\over{2kT}}}e^{-{{mu_y^2}\over{2kT}}}e^{-{{mu_z^...
...^{-{{m(u_x^2+u_y^2+u_z^2)}\over{2kT}}}=e^{-{{mu^2}\over{2kT}}},\end{displaymath}

  • teiseks asendame integreerimisel korrutise $du_xdu_ydu_z$ sfäärilise kihiga faasiruumis $4\pi u^2du$.

    Tulemuseks saame:

    \begin{displaymath}f(u)du=\left(\sqrt{m\over{2\pi
kT}}\right)^3e^{-{{mu^2}\over{2kT}}}\cdot 4\pi u^2du.\end{displaymath}

    Kuidas siit jõuda Maxwelli jaotuse traditsioonilise kujuni, jäägu teie endi leida.

    Mida me näeme? Nii jaotuse enda kui ka eksponentliikme astendaja kordajad sõltuvad suurusest $m/T$. Mida väiksem on see suurus (mida kõrgem temperatuur või mida väiksem osakese mass), seda kiiremini muutub jaotusfunktsioon. Kui tegur on suur (temperatuur madal või osakeste mass suur), on kõver laugjam ja madalam.

    Maxwelli jaotus, nagu teisedki shoti füüsiku J. C. Maxwelli valemid on saadud puhtmatemaatilistest kaalutlustest, tuginedes lisaks katseandmetele veel teatud puhtfilosoofilist laadi eeldustele (molekulid on punktikujulised massiga osakesed, elektri- ja magnetväli ei sõltu keskkonnast jne.). See on uut tüüpi lähenemine füüsikalisele maailmatunnetusele - lähtumine formaal-loogilisest vaatepunktist. Oluline on, et erinevalt Newtoni, Ampere'i, Ohm'i jt. empiirilistest seadustest tuletatakse siin matemaatiliselt ka valemite konstandid. See on meetod, mida tänapäeval tuntakse teoreetilise füüsika nime all.
  • Maxwelli jaotus erinevatel temperatuuridel.


    Maxwelli jaotus pole leitud katsest, vaid tuletatud matemaatiliselt.


     
    Järeldused Maxwelli valemist. Lisaks juba tuntud molekulide ruutkeskmisele kiirusele

    \begin{displaymath}\bar{v^2}=\int\limits _0^\infty v^2f(v)dv={3kT}\over m\end{displaymath}

    saame leida veel keskmise kiiruse

    \begin{displaymath}\bar{v}=\int\limits _0^\infty vf(v)dv=\sqrt{{8kT}\over{\pi m}}\end{displaymath}

    ning, lähtudes funktsiooni ekstreemumi tingimusest

    \begin{displaymath}{{df(v)}\over{dv}}={\rm const}\cdot
e{-{{mv^2}\over{2kT}}}v\left(2-{{mv^2}\over{kT}}\right)=0\end{displaymath}

    molekulide tõenäoseima kiiruse

    \begin{displaymath}v_t=\sqrt{{2kT}\over m}.\end{displaymath}

    Maxwelli jaotuse integreerimise teel on võimalik leida nii molekulide keskmist, ruutkeskmist kui tõenäoseimat kiirust.
    Kulgliikumise jaoks tuletatud jaotust kiiruste järgi saab üldistada suvalise kineetilise energia kohta käivaks, kui asendada Maxwelli jaotuses $v^2={{2E_k}\over m}$.
    Siis
    $dv=1/\sqrt{2mE_k}dE_k$
    ning

    \begin{displaymath}f(E_k)=2\pi^{-{1\over 2}}(kT)^{-{3\over
2}}e^{-{{E_k}\over{kT}}}\sqrt{E_k}.\end{displaymath}

    Saadud valem on oluline erisoojuste klassikalises teoorias ning faasiüleminekute uurimisel.


     
    Maxwelli jaotuse eksperimentaalsest kontrollist. Kiiruste jaotust on suhteliselt lihtne uurida: tuleb tekitada osakeste voog ning juhtida sellest eraldatud lõik liikuvale ekraanile. Praktikas on piisava täpsusega katse teostamine omajagu keeruline ning see õnnestus alles 1920. a.

    Osakeste allikat (hõbetatud plaatinaniit) ümbritsevad kaks sama nurkkiirusega pöörlevat koaksiaalset (samateljelist. ld. co+axis, telg) silindrit, milledest sisemise seina on lõigatud teljega paralleelne pilu. Sellest väljuv molekul jõuab välimise silindrini aja $\Delta
t=(R_v-R_s)/v$ pärast, mistõttu tema langemiskoht ei asu mitte sisemise pilu vastas, vaid on nihkunud $\Delta s=\omega
R_v\Delta t$ võrra. Välimisele silindrile sadestunud hõbedakihi paksus on seetõttu võrdeline vastava kiirusega osakeste esinemissagedusega, ehk, teiste sõnadega - osakeste diferentsiaalse jaotusfunktsiooni $f(v)$ väärtusega antud kiiruse juures.

    Kirjeldatud katset on korratud erinevate ainete ning ka eri tüüpi katseseadmete abil. Tulemused on piisavalt heas kooskõlas Maxwelli valemiga.

    Sellegipoolest loetakse Maxwelli valemi peamiseks väärtuseks tema heuristilist väärtust (kr. heurisko - leian). Teadmine, et mingi süsteemi komponentide energiad jaotuvad kindla seaduspärasuse järgi, võimaldab modelleerida ("panna valemeisse") paljusid juhuslikkusel põhinevaid nähtusi. Üht neist kirjeldame kevadsemestri viimase loengu juurde antud lisatekstis.


    Maxwelli jaotuse kontroll.
    Kiired molekulid (kõvera ots B) jõuavad välimise silindrini varem kui aeglased (ots A). Tekkiva metallikihi paksus vastab Maxwelli jaotusele.


     
    Ülekandenähtused. Maxwelli jaotuse abil saame arvutada ka gaasi mittetasakaalulisi olekuid, kus mõni olekuparameeter (näiteks temperatuur või keemiline koostis) on vaadeldava ruumi eri osades erinev. Selline mittetasakaalulisus toob kaasa energia, massi või liikumishulga ülekandumise ühest ruumiosast teise, siit ka nähtuse nimetus.

    Tehnilistes ülesannetes puutume kõige sagedamini kokku just loetletud kolme suuruse ülekandumisega. Insenerifüüsikas räägitakse vastavalt soojusjuhtivusest (energia ülekanne), difusioonist (massi ülekanne) ja sisehõõrdest (impulsi ehk liikumishulga ülekanne). Kõigil kolmel juhul iseloomustab lähteolekut mingi parameetri (temperatuuri $T$, kontsentratsiooni $c_1$, voolukiiruse $\vec u$) ebaühtlane jaotus, mida saab kirjeldada väljafunktsiooniga (parameetri väärtus sõltuvalt asukohast). Ülekande kiirus on alati võrdeline vaadeldava suuruse gradiendiga (vt. loeng 5), võrdetegurit aga saab leida keskmise kiiruse kaudu.


    Ülekandenähtuste all mõistame soojuse, keemilise koostise, impulsi või muu parameetri levikut gaasis.

    Näitena vaatame lihtsamat juhtu, kus parameetri muutus toimub kindlas suunas, näiteks $z$-telje sihis. Gradient asendub sel juhul lihtsa tuletisega ning valemiteks saame:

    \begin{displaymath}{\rm difusioon:}~~M_1=-D{{dc_1}\over{dz}}S;~~D~-~{\rm
difusioonitegur;}\end{displaymath}

    \begin{displaymath}{\rm soojusjuhtivus:}~~Q=-\kappa{{dT}\over{dz}}S;~~\kappa~-~{\rm
soojusjuhtivustegur;}\end{displaymath}

    \begin{displaymath}{\rm siseh\~o\~ore:}~~F=\eta{{du}\over{dz}}S;~~\eta~-~{\rm
siseh\~o\~ordetegur;}\end{displaymath}

    Viimane valem väljendab meile juba hüdrodünaamikast tuntud hõõrdejõudu, teised aga ajaühikus pinda $S$ läbinud gaasi massi $M$ või soojushulka $Q$. Kõik konstandid on avaldatavad keskmise kiiruse $\bar v$ ning molekuli vaba tee pikkuse $\lambda=1/\sqrt 2\sigma n$ kaudu:

    \begin{displaymath}D={1\over 3}\bar v\lambda;\end{displaymath}

    \begin{displaymath}\kappa={1\over 3}\rho c_V\bar v\lambda;\end{displaymath}

    \begin{displaymath}\eta={1\over 3}\bar v\lambda.\end{displaymath}

    Keskmine kiirus on leitav Maxwelli jaotusest; vaba tee pikkus sõltub aga gaasi tihedusest $n$ ning ühest senini empiirilisest suurusest, nn molekuli efektiivsest ristlõikepindalast (väljendab molekulaarjõudude mõjuraadiust) $\sigma$:

    \begin{displaymath}\lambda={1\over {\sqrt 2\sigma n}}.\end{displaymath}


    Soojusülekanne ja difusioon.
    Temperatuuride (kontsentratsioonide) erinevus T1 - T2 (või konsentratsioonide erinevus c1 - c2) toob kaasa energia (või massi) voo läbi pinna S.

    Sisehõõrdumine.
    Aeglasemad molekulid liiguvad läbi pinna S, aeglustades gaasi voolamist ülemistes kihtides.
    Algusesse
    Sisukorda