AAM  >> Vol. 6 No. 9 (December 2017)

    当碰撞核为零时Bolzmann方程特殊解的存在性
    Existence of Special Solutions to Boltzmann Equations, When the Collision Kernel Is Zero

  • 全文下载: PDF(406KB) HTML   XML   PP.1115-1118   DOI: 10.12677/AAM.2017.69135  
  • 下载量: 560  浏览量: 1,086  

作者:  

张 飞:云南民族大学数学与计算机科学学院,云南 昆明

关键词:
Boltzmann方程解的存在性碰撞核Boltzmann Equations Existence of solution Collision Kernel

摘要:

本文构造了玻尔兹曼方程的两个特殊解,并证明了该解的存在唯一性。

In this paper, we will give a rigorous proof of existence and uniqueness of special solutions to Boltzmann equations, when the collision kernel is zero.

1. 引言

路德维希·玻尔兹曼(Ludwig Boltzmann,1844年2月20日~1906年9月5日),奥地利物理学家,热力学和统计力学的奠基人之一。

玻尔兹曼1844年出生于奥地利的维也纳,1866年获得维也纳大学博士学位。

玻尔兹曼的贡献主要在热力学和统计物理方面。1869年,他将麦克斯韦速度分布律推广到保守力场作用下的情况,得到了玻尔兹曼分布律。1872年,玻尔兹曼建立了玻尔兹曼方程(又称输运方程),用来描述气体从非平衡态到平衡态过渡的过程。

Boltzmann equation又称为玻尔兹曼输运方程,它就是分布函数法中所采用的一种方程,即是非平衡分布函数 F ( x , t , v ) 所满足的一个方程,求解此方程可得到不同条件下的 F ( x , t , v ) ,然后即可求出电子的各种输运参量。玻尔兹曼方程是经典粒子牛顿力学运动模型和能态跃迁的量子力学模型相糅合的产物。如果忽略所有的相干效应,经过一定的简化,可以从量子输运模型中推导出玻尔兹曼方程。

玻尔兹曼输运方程中考虑到了载流子的速度分布和散射的方向性,因此较为精确。

在有电场或温度梯度等外场的情况下,根据分布函数因电场、磁场、温度梯度等外场而引起的漂移变化以及因散射而引起的变化,即可建立起Boltzmann方程,由于其中的散射项应是一个对散射几率的积分, 所以Boltzmann方程是一个微分–积分方程。该方程的求解很复杂,通常采用近似方法,常用的一种近似方法就是弛豫时间近似 [1] 。

玻尔兹曼方程是一个高维的方程,三维波矢空间(v),三维实空间(x),再加上一维时间(t),难于求解,常用蒙特卡罗方法来模拟。本文讨论当碰撞核为零时玻尔兹曼方程的特殊解构造。

2. 准备知识

在碰撞核理论中,气体分子是由概率分布函数 F F ( t , x , v ) 0 ,测量气体分子在时间 t R + ,空间为 x R 3 和瞬时速度 v R 3 。这种函数,通常称为分布函数或者密度函数。下面为方便讨论们讨论方便,我们定义玻尔兹曼方程:

t F + v x = B ( F , F ) (1.1)

这里 B ( F , F ) 是玻尔兹曼碰撞积。这种碰撞积只有在辐角为的 w 密度函数下,并由下式给出:

B ( F , F ) ( x , t , v ) = R 3 × S 2 ( F F 1 F F 1 ) b ( v v 1 , w ) | cos ( v v 1 , w ) | d w d v 1 (1.2)

F 1 , F F 1 分别由指定的 F ( x , t , v 1 ) , F ( x , t , v ) F ( t , x , v 1 ) ,在 v 1 , v v 1 下给出,且 v 1 R 3 , w S 2

v = v ( v v 1 ) w w , v 1 = v 1 + ( v v 1 ) w w . (1.3)

v + v 1 = v + v 1 , | v | 2 + | v 1 | 2 = | v | 2 + | v 1 | 2 (1.4)

关系(1.4)是每一个气体分子之间碰撞的动量守恒和动能守恒定律 [2] 。

3. 构造Boltzmann方程的两个特殊解

问题一:若 F = exp ( | v | 2 ) 时, F = exp ( | v | 2 ) 总是方程(1.1) (1.2)的一个解。

证明:首先检验 F = exp ( | v | 2 ) 是方程(1.1) (1.2)的一个解。

由(1.4)碰撞前后的速度关系,可以得到如下关系式

F F 1 = exp ( | v | 2 + | v 1 | 2 ) exp ( | v | 2 + | v 1 | 2 ) = 0 (1.5)

所以 B ( F , F ) = 0 恒成立

F = exp ( | v | 2 ) 是方程(1.1)(1.2)的一个解。

然后我们再证明,当 [ ( | v A x c t | 2 + | v 1 A x c t | 2 ) ] exp [ ( | v A x c t | 2 + | v 1 A x c t | 2 ) ]

也是方程(1.1)(1.2)的一个解,且能化成上述形式。

其中, A = ( a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33 ) = ( 0 0 0 0 0 0 0 0 0 )

从而 F = exp ( | v A x c t | 2 ) = exp β

β = ( v ˜ 1 i = 1 3 a 1 i x i + c 1 t ) 2 + ( v ˜ 2 i = 1 3 a 2 i x i + c 2 t ) 2 + ( v ˜ 3 i = 1 3 a 3 i x i + c 3 t ) 2

下面计算,
t F = exp β t β = i = 1 3 a 1 i ( v ˜ 1 i = 1 3 a 1 i x i + c 1 t ) v ˜ 1 i = 1 3 a 2 i ( v ˜ 2 i = 1 3 a 2 i x i + c 2 t ) v ˜ 2 i = 1 3 a 3 i ( v ˜ 3 i = 1 3 a 3 i x i + c 3 t ) v ˜ 3

v x F = ( v ˜ 1 v ˜ 2 v ˜ 3 ) ( 2 i = 1 3 a 1 i x i ( v ˜ 1 i = 1 3 a 1 i x i + c 1 t ) 2 i = 1 3 a 2 i x i ( v ˜ 2 i = 1 3 a 2 i x i + c 2 t ) 2 i = 1 3 a 3 i x i ( v ˜ 3 i = 1 3 a 3 i x i + c 3 t ) ) = 2 i = 1 3 a 1 i x i ( v ˜ 1 i = 1 3 a 1 i x i + c 1 t ) v ˜ 1 2 i = 1 3 a 2 i x i 2 ( v ˜ 2 i = 1 3 a 2 i x i + c 2 t ) v ˜ 2 2 i = 1 3 a 3 i x i ( v ˜ 3 i = 1 3 a 3 i x i + c 3 t ) v ˜ 3

下面考察 B ( F , F ) 这一项,有 B ( F , F ) = 0

F F 1 F F 1 = exp [ ( | v A x c t | 2 + | v 1 A x c t | 2 ) ] exp [ ( | v A x c t | 2 + | v 1 A x c t | 2 ) ] = 0

记得我们曾经学习过这样一条性质 [3] :

a 2 + b 2 = c 2 + d 2 ,如果 a + b = c + d ,则有 ( a + m ) 2 + ( b + n ) 2 = ( c + m ) 2 + ( d + n ) 2

又由(1.4) v + v 1 = v + v 1 , | v | 2 + | v 1 | 2 = | v | 2 + | v 1 | 2

由上述性质可 F F 1 F F 1 = exp [ ( | v A x c t | 2 + | v 1 A x c t | 2 ) ] exp [ ( | v A x c t | 2 + | v 1 A x c t | 2 ) ] = 0

B ( F , F ) ( x , t , v ) = R 3 × S 2 ( F F 1 F F 1 ) b ( v v 1 , w ) | cos ( v v 1 , w ) | d w d v 1 = 0

exp β [ c 1 ( v ˜ 1 i = 1 3 a 1 i x i + c 1 t ) + c 2 ( v ˜ 2 i = 1 3 a 2 i x i + c 2 t ) + c 3 ( v ˜ 3 i = 1 3 a 3 i x i + c 3 t ) ]

i = 1 3 c i v i = 1 3 j = 1 3 c j a j i x i + i = 1 3 c i 2 t = 0

其中 β = ( v ˜ 1 i = 1 3 a 1 i x i + c 1 t ) 2 + ( v ˜ 2 i = 1 3 a 2 i x i + c 2 t ) 2 + ( v ˜ 3 i = 1 3 a 3 i x i + c 3 t ) 2

于是令 c 1 ( v ˜ 1 i = 1 3 a 1 i x i + c 1 t ) + c 2 ( v ˜ 2 i = 1 3 a 2 i x i + c 2 t ) + c 3 ( v ˜ 3 i = 1 3 a 3 i x i + c 3 t ) = 0

而且 i = 1 3 a 1 i ( v ˜ 1 i = 1 3 a 1 i x i + c 1 t ) v ˜ 1 i = 1 3 a 2 i ( v ˜ 2 i = 1 3 a 2 i x i + c 2 t ) v ˜ 2 i = 1 3 a 3 i ( v ˜ 3 i = 1 3 a 3 i x i + c 3 t ) v ˜ 3 = 0

通过整理得 k = 1 3 ( j = 1 3 ) c i v i = 1 3 j = 1 3 c j a j i x i + i = 1 3 c i 2 t = 0

而且 k = 1 3 ( j = 1 3 a j k ) ( v ˜ k 2 i = 1 3 a k i x i v ˜ k + c 3 v ˜ k t ) = 0

所以 c = ( c 1 c 2 c 3 ) = ( 0 0 0 )

A = ( a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33 ) = ( 0 0 0 0 0 0 0 0 0 )

于是形如 F = exp ( | v A x c t | 2 ) 的解只有 F = exp ( | v | 2 ) 这一解。

问题二:当 F = exp ( | a x ( a t + b ) v c | 2 ) = exp ( β ) 也是方程(1.1) (1.2)的一个解

其中 v = ( v ˜ 1 v ˜ 2 v ˜ 3 ) , x = ( x ˜ 1 x ˜ 2 x ˜ 3 ) , c = ( c 1 c 2 c 3 ) , a R

可计算 β = i = 1 3 [ a x i ( a t + b ) v i c i ] 2

F t = 2 a i = 1 3 v ˜ i [ a x i ( a t + b ) v i c i ]

v x F = 2 a i = 1 3 v ˜ i [ a x i ( a t + b ) v i c i ]

再由上述性质可知: F F 1 F F 1 = 0

F t + v x F = 0

F = exp ( | a x ( a t + b ) v c | 2 ) 是玻尔兹曼方程(1.1) (1.2)的一个解。

文章引用:
张飞. 当碰撞核为零时Bolzmann方程特殊解的存在性[J]. 应用数学进展, 2017, 6(9): 1115-1118. https://doi.org/10.12677/AAM.2017.69135

参考文献

[1] 四川师大物理系, 编. 统计物理学进展[M]. 1989.
[2] Boardos, C., Golse, F. and Levermore, C.D. (1993) Fluid Dynamic Limits of Kinetic Equations II: Convergence Proofs for Boltzmann Equation. Communications on Pure and Applied Mathematics, 46, 667-753.
[3] 张恭庆, 林源渠. 泛函分析讲义(上册) [M]. 北京: 北京大学出版社, 2006.