משוואת ולסוב

מתוך ויקיפדיה, האנציקלופדיה החופשית
קפיצה אל: ניווט, חיפוש

משוואת ולסוב היא משוואה דיפרנציאלית המתארת את ההתפתחות בזמן של פונקציית ההתפלגות של פלזמה המכילה חלקיקים טעונים עם אינטראקציה רחוקת-טווח (למשל, קולונית). המשוואה הוצעה לראשונה לתיאור פלזמה על ידי הפיזיקאי הרוסי אנטולי ולסוב בשנת 1938. [1] [2] [3]

הבעייתיות במודל הקינטי הסטנדרטי[עריכת קוד מקור | עריכה]

ולסוב טען כי המודל הקינטי המבוסס על משוואת בולצמן מתקשה בתיאור פלזמה, מכיוון שהמערכת כוללת אינטראקציות רחוקות-טווח כמו האינטראקציה הקולונית. ולסוב ציין את הבעיות הבאות העולות בניסיון להחיל את המודל הקינטי המבוססת על התנגשות של זוג חלקיקים לדינמיקה של פלזמה:

  1. התאוריה של התנגשויות בין זוג חלקיקים נמצאת בסתירה לגילויים של ריילי, אירווינג לאנגמיור ולואי טונקס על תנודות טבעיות בפלזמת אלקטרונים.
  2. התאוריה של התנגשויות בין זוג חלקיקים לא ניתנת ליישום עבור אינטראקצייה קולונית מכיוון שהאיברים הקינטיים מתבדרים.
  3. התאוריה של התנגשויות בין זוג חלקיקים לא מסוגלת להסביר את הניסויים שבוצעו על ידי האריסון מריל והארולד ווב בנושא פיזור אנומלי של אלקטרונים בפלזמה גזית. [4]

ולסוב העלה השערה כי האופי ארוך-הטווח של האינטראקציה הקולונית הוא המקור לקשיים האלו. ולסוב החל לפתח את משוואת בולצמן עבור מקרה ללא התנגשויות, המוצגת בקואורדינטות מוכללות:

\frac{\text{d}f(\bar{q},\bar{p},t)}{\text{d}t}=0,

ובמפורש, כמשוואה דיפרנציאלית חלקית:

\frac{\partial f}{\partial t}+\frac{\partial \bar{q}}{\partial t}\cdot\frac{\partial f}{\partial \bar{q}}+\frac{\partial \bar{p}}{\partial t} \cdot \frac{\partial f}{\partial \bar{p}}=0,

התאים אותה למקרה של פלזמה, וקיבל את מערכת המשוואות המפורטות להלן.

מערכת המשוואות של ולסוב-מקסוול[עריכת קוד מקור | עריכה]

במקום התאוריה הקינטית, המבוססת על התנגשויות, ולסוב השתמש בשדה קולקטיבי עקבי הנוצר על ידי חלקיקי הפלזמה הטעונים על-מנת לתאר את האינטראקציה בין חלקיקי הפלזמה. בתיאור זה מתבצע שימוש בפונקציות ההתפלגות f_e(\vec{r},\vec{p},t) ו-f_i(\vec{r},\vec{p},t) עבור האלקטרונים והיונים. פונקציית ההתפלגות f_\alpha(\vec{r},\vec{p},t) עבור חלקיק מסוג כללי \alpha, בעל תנע \vec{p} בקרבת מקום \vec{r} בזמן t. בתיאור זה, במקום משוואת בולצמן, מתקבלת המערכת הבאה של משוואות לתיאור הפלזמה (ביחידות cgs):

\frac{\partial f_e}{\partial t}+\vec{v}_e \cdot \nabla f_e - e\left[ \vec{E}+\frac{1}{c}\left(\vec{v} \times \vec{B}\right)\right] \cdot \frac{\partial f_e}{\partial \vec{p}} = 0
\frac{\partial f_i}{\partial t}+\vec{v}_i \cdot \nabla f_i - e\left[ \vec{E}+\frac{1}{c}\left(\vec{v} \times \vec{B}\right)\right] \cdot \frac{\partial f_i}{\partial \vec{p}} = 0
\nabla \times \vec{B} = \frac{4 \pi \vec{j}}{c}+\frac{1}{c}\frac{\partial \vec{E}}{\partial t}, ~~ \nabla \times \vec{E} = -\frac{1}{c}\frac{\partial{\vec{B}}}{\partial{t}}
\nabla \cdot \vec{E} = 4\pi \rho, ~~ \nabla \cdot \vec{B} = 0
\rho = e \int \left( f_i-f_e \right) d^3p, ~~ \vec{j} = e \int \left( f_i - f_e \right) \vec{v} d^3 p, ~~ \vec{v}_\alpha = \frac{\vec{p}/m_\alpha}{\left[1+p^2/\left(m_\alpha c\right)^2\right]^{1/2}}

כאשר e מטען האלקטרון, c מהירות האור, m_\alpha מסת חלקיק מסוג \alpha, \vec{E}(\vec{r},t) ו-\vec{B}(\vec{r},t) השדות החשמלי והמגנטי בנקודה \vec{r} בזמן t, בהתאמה. ההבדל הבסיסי בין מערכת המשוואות האלו ומערכת משוואות המתארת חלקיקים בהשפעת שדה אלקטרומגנטי חיצוני נעוץ בעובדה שהשדות במשוואת ולסוב מקיימים אינטראקציה דו-כיוונית עם החלקיקים מכיוון שהתפלגות החלקיקים הטעונים משפיעה על השדה האלקטרומגנטי, ולהיפך.

משוואת ולסוב-פאוסון[עריכת קוד מקור | עריכה]

משוואות ולסוב-פואסון הן קירוב למשוואות ולסוב-מקסוול עבור מקרה לא-יחסותי ללא שדה מגנטי:

\frac{\partial f_\alpha}{\partial t} + \vec{v} \cdot \frac{\partial f_\alpha}{\partial \vec{x}} + \frac{q_a\vec{E}}{m_\alpha} \cdot \frac{\partial f_\alpha}{\partial \vec{v}} = 0,

ומשוואת פואסון עבור השדה החשמלי:

\epsilon \nabla^2\phi+\rho=0.

כאשר q_\alpha מטען החלקיק, m_\alpha מסת החלקיק, \vec{E}(\vec{x},t) השדה החשמלי, \phi(\vec{x},t) הפוטנציאל החשמלי ו-\rho צפיפות המטען החשמלי. משוואות ולסוב-פואסון משמשות לתיאור תופעות רבות בפלזמה בעלות מאפיינים לא-מקסווליאנים בולטים, למשל דעיכת לנדאו ושכבות כפולות, אשר לא ניתן לתאר באמצעות מודלים קינטיים.

משוואות המומנט[עריכת קוד מקור | עריכה]

במודל נוזלי לתיאור פלזמה, למשל מגנטו-הידרודינמיקה (MHD), התפלגות המהירויות איננה מחושבת. עובדה זו מושגת על ידי החלפת f(\vec{r},\vec{v},t) במומנטים של הפלזמה כמו צפיפות מספרית, n, מהירות ממוצעת, \mathbf{u}, וצפיפות, \mathbf{p}. גדלים אלו מכונים מומנטים מכיוון שהמומנט מסדר n של f מתקבל על ידי אינטגרציה של v^n f לפי המהירות. המשתנים המתקבלים תלויים במיקום ובזמן בלבד, כלומר מידע מסוים אבד. בתאוריה מרובת-נוזלים, כל סוג חלקיקים מתואר כנוזל נפרד עם ערכים של לחץ, צפיפות ומהירות זרימה. המשוואות המתארות את מומנטי הפלזמה מכונות משוואות מומנט או משוואות זרימה. שתי משוואות המומנט הנפוצות ביותר מוצגות בהמשך (במערכת היחידות הבינלאומית). קבלת משוואות המומנט ממשוואת ולסוב לא מצריכה הנחות לגבי פונקציית ההתפלגות.

משוואת הרציפות[עריכת קוד מקור | עריכה]

משוואת הרציפות מתארת את שינוי הצפיפות עם הזמן. ניתן לקבלה על ידי אינטגרציה של משוואת ולסוב על-פני מרחב המירויות כולו.

\int\frac{\text{d}}{\text{d}t}fd^3v = \int\left[ \frac{\partial}{\partial t}f + (\vec{v} \cdot \nabla_r)f+\nabla_v \cdot (\vec{a}f) \right]d^3v = 0

לאחר חישובים, ניתן לקבל

\frac{\partial}{\partial t}n + \nabla \cdot (n \mathbf{u})=0

צפיפות החלקיקים, n, והמהירות הממוצעת, \mathbf{u}, הם המומנטים מסדר אפס וסדר ראשון:

n = \int fd^3v
\mathbf{u} = \frac{\int \vec{v}fd^3v}{\int fd^3v} = \frac{1}{n} \int \vec{v} fd^3 v

משוואת התנע[עריכת קוד מקור | עריכה]

קצב השינוי בתנע של חלקיק ניתן על ידי הצבת כוח לורנץ במשוואת החוק השני של ניוטון:

m\frac{\text{d}\vec{v}}{\text{d}t} = q\left(\vec{E} + \vec{v} \times \vec{B} \right)

על ידי שימוש במשוואה זו ומבשוואת ולסוב, משוואת התנע עבור כל נוזל מתקבלת כ:

mn\frac{\text{d}}{\text{d}t}\textbf{u} = -\nabla \cdot \mathbf{p} + qn\vec{E} + qn \mathbf{u} \times \vec{B},

כאשר \mathbf{p} הוא טנזור הלחץ. הנגזרת המלאה לפי הזמן מתקבלת:

\frac{\text{d}}{\text{d}t} = \frac{\partial}{\partial t} + \mathbf{u} \cdot \nabla.

טנזור הלחץ מוגדר כצפיפות המסה כפול מטריצת הקו-וריאנס של המהירות:

p_{ij} = m \int (v_i-u_i)(v_j-u_j)fd^3v.

קירוב השדה הקפוא[עריכת קוד מקור | עריכה]

בדומה למקרה של MHD-אידאלי, ניתן להתייחס את הפלזמה כאל "קשורה" לקווי השדה המגנטי בתנאים מסוימים. ניתן לומר שקווי השדה המגנטי "קפואים" בתוך הפלזמה. ניתן לפתח את התנאים לקירוב זה ממשוואת ולסוב.

נגדיר את הסקלות \mathbf{L, T} ו-\mathbf{V} עבור זמן, מרחק ומהירות, בהתאמה. סקלות אלו מייצגות מדד לפרמטרים השונים המביאים לשינוי גדול ב-f, כאשר שינוי גדול משמעותו

\frac{\partial f}{\partial t}T \sim f, ~~ \left| \frac{\partial f}{\partial \vec{r}} \right|L \sim f, ~~ \left| \frac{\partial f}{\partial \vec{v}} \right|V \sim f.

ניתן לכתוב:

t' = \frac{t}{T}~~, ~~ \vec{r}~'=\frac{\vec{r}}{L}~~, ~~ \vec{v}~' = \frac{\vec{u}}{V}.

כעת ניתן להציג את משוואת ולסוב כך:

\frac{1}{T}\frac{\partial f}{\partial t'} + \frac{V}{L} \vec{v}~' \cdot \frac{\partial f}{\partial \vec{r}} + \frac{q}{mV} \left( \vec{E} + V\vec{v}~'\times \vec{B} \right) \cdot \frac{\partial f}{\partial \vec{v}~'} = 0

עד כה לא בוצעו קירובים. כעת נקבע V = R \omega_g, כאשר \omega_g = qB/m היא תדירות הגירו ו-R הוא רדיוס הגירו (ראה פרמטרים של פלזמה). על ידי חלוקת המשוואה ב-\omega_g נקבל

\frac{1}{\omega_gT}\frac{\partial f}{\partial t'} + \frac{R}{L} \vec{v}~' \cdot \frac{\partial f}{\partial \vec{r}~'} + \left( \frac{\vec{E}}{VB} + \vec{v}~' \times \frac{\vec{B}}{B} \right) \cdot \frac{\partial f}{\partial \vec{v}~'} = 0

בקירוב בו 1/\omega_g \ll T ו-R \ll L, שני האיברים הראשונים יהיו קטנים הרבה יותר מ-f, מכיוון ש-\partial f / \partial t' \sim f ~,~v' \lesssim 1 ~,~\partial f / \partial \vec{r}~' \sim f. האיבר האחרון הוא מסדר-גודל של f, ולכן ניתן להזניח את שני האיברים הראשונים ולקבל

\left( \frac{\vec{E}}{VB} + \vec{v}~' \times \frac{\vec{B}}{B} \right) \cdot \frac{\partial f}{\partial \vec{v}~'} \approx 0 \Rightarrow \left( \vec{E} + \vec{v} \times {\vec{B}} \right) \cdot \frac{\partial f}{\partial \vec{v}} \approx 0

ניתן לפרק את המשוואה לחלקים מקביל וניצב:

\vec{E}_\parallel \frac{\partial f}{\partial \vec{v}_\parallel} + \left( \vec{E}_\perp + \vec{v} \times \vec{B} \right) \cdot \frac{\partial f}{\partial \vec{v}_\perp} \approx 0

כעת נשים לב לכך ש- \vec{v} = \vec{v}_0 + \Delta\vec{v}, כאשר

\vec{v}_0 \times \vec{B} = -\vec{E}_\perp

וכעת ניתן לקבל

\vec{E}_\parallel \frac{\partial f}{\partial \vec{v}_\parallel} + \left( \Delta\vec{v}_\perp \times \vec{B} \right) \cdot \frac{\partial f}{\partial \vec{v}_\perp} \approx 0

עבור שדה מקביל קטן נקבל

\left( \Delta\vec{v}_\perp \times \vec{B} \right) \cdot \frac{\partial f}{\partial \vec{v}_\perp} \approx 0

המשמעות של משוואה זו היא שההתפלגות היא גירו-טרופית. המהירות הממוצעת של התפלגות גירו-טרופית היא אפס, ולכן \vec{v}_0 שווה זהותית למהירות הממוצעת, \mathbf{u}, ונקבל

\vec{E} + \mathbf{u} \times \vec{B} \approx 0

לסיכום, זמן המחזור של תנועת הגירו ורדיוס הגירו קטנים הרבה יותר, בהתאמה, מסקלות הזמן והמרחק שיגרמו לשינויים גדולים בפונקציית ההתפלגות. יש לחשב את התנאים לקירוב השדה הקפוא לכל סוג חלקיקים בנפרד. מכיוון שלאלקטרונים זמן מחזור גירו ורדיוס גירו קצרים יותר מאשר ליונים, התנאים יתקיימו עבורם לעתים קרובות יותר.

ראו גם[עריכת קוד מקור | עריכה]

הערות שוליים[עריכת קוד מקור | עריכה]

  1. ^ A. A. Vlasov (1938). "On Vibration Properties of Electron Gas" (בRussian). J. Exp. Theor. Phys. 8 (3): 291. 
  2. ^ A. A. Vlasov (1968). "The Vibrational Properties of an Electron Gas". Soviet Physics Uspekhi 10 (6): 721. doi:10.1070/PU1968v010n06ABEH003709. Bibcode1968SvPhU..10..721V. 
  3. ^ A. A. Vlasov. Theory of Vibrational Properties of an Electron Gas and Its Applications, 1945. 
  4. ^ H. J. Merrill and H. W. Webb (1939). "Electron Scattering and Plasma Oscillations". Physical Review 55 (12). doi:10.1103/PhysRev.55.1191. Bibcode1939PhRv...55.1191M.