玻爾茲曼方程

玻爾茲曼方程或玻爾茲曼輸運方程的變化的動力學方程。波爾茲曼方程並不去確定流體中每個粒子的位置和動量，而是求出具有特定位置和動量的粒子的機率分佈。具體而言，考慮某一瞬間，以位置向量 formula_1 末端為中心的無窮小區域內，動量無限接近動量向量 formula_2的粒子的機率分佈。波爾茲曼方程可用於確定物理量是如何變化的，例如流體在輸運過程中的熱能和動量；還可由此推匯出其他的流體特徵性質，例如黏度，熱導率，以及電阻率，詳見對流擴散方程式。波爾茲曼方程是一個非線性的。方程中的未知函式是一個包含了粒子空間位置和動量的六維機率密度函式。方程解的存在性和唯一性問題仍然沒有完全解決，但就最近發表的一些工作而言，對於解決這一問題還是有一定希望的。概述.相空間與密度函式.系統中所有可能的位置 formula_1 和動量 formula_2 組成的集合被稱作此係統的相空間，其中位置座標記為 formula_6，動量座標記為 formula_7。整個空間是六維的：空間中某一點的座標可表示為 formula_8，每個座標均透過時間 formula_9 引數化。微元可寫作：formula_10波爾茲曼方程的核心是“formula_11”函式，它表示的是在一段極短的時間內，每一相空間單位體積中的formula_12個分子在微元 formula_13 中，位置都為 formula_1 且動量都為 formula_2 的機率。透過定義，我們可使機率密度函式 formula_16 滿足以下條件：formula_17formula_18 被定義為在時間 formula_9，位於 formula_20 的空間元 formula_13 中的粒子總數:61-62。對座標空間與動量空間的一個區域積分即可得該區域內所有具有對應位置和動量的粒子的總數：formula_22雖然formula_11是和一群粒子相關的，但此相空間是對於單一粒子的。這裡不使用 r1, p1 表示粒子1，r2, p2 表示粒子2，……，rN, pN 表示粒子N。系統中的粒子被假定是相同的。對於具有超過一種化學組分的混合物，每一種成分都需要有一個分佈函式，見下文。一般形式.方程的一般形式可以寫作：formula_25這裡“force”一詞指的是外部對粒子施加的力，“diff”表示粒子的擴散，“coll”表示粒子的碰撞，指的是碰撞中粒子間相互的作用力。上述三項的具體形式將會在下文給出。注意，一些作者會使用粒子的速度 formula_26，來代替上文的 formula_2；這兩個物理量可以透過定義formula_28來聯絡。“force”項和“diff”項.考慮一群以formula_11分佈的粒子。每個粒子均受到外力formula_30的作用（不包括粒子間作用力。粒子間的作用見後面對“coll”項的處理)。假設在時間 formula_9，一定數量的粒子都有位置 formula_1。如果此時有一個力formula_30在這一瞬作用在每個顆粒上，那麼在時間 formula_37，它們的位置將會是formula_38，動量將變成 formula_39。在沒有碰撞的情況下，formula_11必須滿足formula_41這裡，注意到相空間元 formula_13 是恆定的這個事實可以從哈密頓方程得知。然而，由於存在碰撞，相空間元 formula_13 中的粒子密度是可變的，所以- f(\mathbf{r},\mathbf{p},t)d^3\mathbf{r}d^3\mathbf{p} \\&amp; = \Delta f d^3\mathbf{r}d^3\mathbf{p}\end{align}其中 formula_44 指的是formula_11的總變化量。式除以 formula_46 並取極限 formula_47 和 formula_48 可得formula_11的全微分：\cdot d\mathbf{p} \\&amp; = \frac{\partial f}{\partial t}dt +\nabla f \cdot \frac{\mathbf{p}dt}{m} + \frac{\partial f}{\partial \mathbf{p}}\cdot \mathbf{F}dt\end{align}其中 ∇ 為梯度算符，· 為點積，formula_50是∇的動量類比的一個簡寫，ê"x", ê"y", ê"z" 為笛卡爾座標系下的單位向量。最終形式.對可得：formula_51這裡，formula_52 為流體中作用在粒子上的力場，formula_24為粒子質量。 右邊的一項用於描述粒子間相互碰撞產生的影響；如果此項為零，則說明粒子之間沒有碰撞。無碰撞情況下的玻爾茲曼方程常被稱為。這個方程比上一節“主要論述”中的一般形式更加有用。然而這個方程依舊是不完整的：除非已知formula_11中的碰撞項，否則formula_11是解不出來的。這一項並不像其他項一樣可以簡單地或一般地得到——這一項是表示粒子的碰撞的統計項，需要知道粒子遵守怎樣的統計規律，例如麥克斯韋-玻爾茲曼分佈，費米-狄拉克分佈或玻色–愛因斯坦分佈。碰撞項和分子混沌.玻爾茲曼的一個關鍵見解就是對碰撞項的確定。他假設的碰撞項完全是由假定在碰撞前不相關的兩個粒子的相互碰撞得到的。這個假設被波爾茲曼稱為“Stosszahlansatz”，也叫做“”。根據這一假設，碰撞項可以被寫作單粒子分佈函式的乘積在動量空間上的積分：formula_56其中 formula_57 和 formula_58 表示碰撞前任意兩個粒子的動量， formula_61 和 formula_62 表示碰撞後的動量formula_63指對應動量的大小，formula_64 是碰撞的微分散射截面。對碰撞項的簡化.求解波爾茲曼方程時，許多挑戰都來自於其複雜的碰撞項；因此我們會做一些對碰撞項“建模”和簡化的嘗試。現知最好的模型是由Bhatnagar，Gross和Krook作出的。BGK近似中假設分子的碰撞會迫使一個物理空間中的某一點的非平衡分佈函式回到麥克斯韋平衡分佈函式，且其發生率正比於分子碰撞頻率。於是，波爾茲曼方程可被寫作以下的BGK形式：(也叫做“馳豫時間近似”，relaxation time approximation)formula_65其中 formula_66 是分子碰撞頻率，和馳豫時間 formula_67 具有倒數關係：formula_68。formula_69是此處局域的麥克斯韋分佈函式，由空間中這一點的氣體溫度給定。普適方程.對於具有多種化學組分的混合物，我們以 i =1，2，3，……，n 標記各種成分。則對於組分i的方程是：formula_70其中 formula_71。碰撞項為formula_72其中 formula_73，相對動量的大小是formula_74Iij 是粒子i和粒子j之間的微分散射截面。此積分的和描述的是某一相空間元中，組分i粒子的進出。應用和推廣.守恆方程.玻爾茲曼方程可用於推導流體動力學中的質量守恆，電量守恆，動量守恆，以及能量守恆定律:p 163。對於只含有一種粒子的流體，粒子數密度 formula_75 為：formula_76算符 A 的期望值由下式給出：formula_77由於守恆方程中包含張量，以下使用愛因斯坦求和約定簡化標記，即 formula_78 且 formula_79，其中 formula_80 為粒子速度向量。定義某函式 formula_81，使得其唯一的自變數為動量 formula_82。假設力 formula_83 為位置的函式，且對於 formula_84，formula_11 為0。對玻爾茲曼方程兩邊同乘 formula_86 ，並對動量積分可得如下四項：formula_87formula_88formula_89formula_90因為 formula_86 在碰撞中守恆，所以最後一項為零。令 formula_92，即粒子質量，積分後的玻爾茲曼方程化為質量守恆方程:pp 12,168：formula_93formula_94 為質量密度，formula_95 為平均流體速度。令 formula_96，即粒子動量，積分後的玻爾茲曼方程化為動量守恆方程:pp 15,169：formula_97formula_98 為壓強張量。令 formula_99，即粒子動能，積分後的玻爾茲曼方程化為能量守恆方程:pp 19,169：formula_100formula_101 為動力熱能密度，formula_102 熱通量向量。哈密頓力學.在哈密頓力學中, 玻爾茲曼方程通常寫作formula_103其中 L 是劉維爾運算元(這裡定義的劉維爾運算元和連結文章中的定義不一致)，它描述了相空間體積的演化；C 是碰撞運算元。非相對論下的L 寫作formula_104波爾茲曼方程的解.直到2010年，波爾茲曼方程的準確解才在數學上被證明是的。這意味著，如果對服從波爾茲曼方程的系統施加一個微擾，此係統最終將回到平衡狀態，而不是發散到無窮，或表現出其他的行為。然而，這種存在性證明是無助於我們在現實問題中求解該等式的。 事實上，這個結論只告訴我們某種特定條件下的解是否存在，而不是如何找到他們。在實踐中，數值計算方法被用於尋找各種形式的波爾茲曼方程的近似解，應用範圍從稀薄氣流中的高超音速空氣動力學，到等離子體的流動中都可以見到。