【等离子体】麦克斯韦速度/速率分布函数

【等离子体】麦克斯韦速度/速率分布函数1 基础知识 概率密度 分布函数的理解

大家好,欢迎来到IT知识分享网。

麦克斯韦速度/速率分布函数

1、基础知识——概率密度/分布函数的理解

  • 离散随机变量:概率函数 P ( x ) P(x) P(x)、概率分布 F ( X ) F(X) F(X)

    P ( x = x i ) = N x i N P(x=x_i)=\frac{N_{x_i}}{N} P(x=xi)=NNxi F ( X ) = ∑ P ( x ≤ X ) F(X)=\sum_{}{}P(x≤X) F(X)=P(xX)

  • 连续随机变量:概率密度函数 f ( x ) f(x) f(x)、概率分布 F ( X ) F(X) F(X)
    • 不能简单将离散随机变量中的概率函数概念简单延伸至连续随机变量
    • 抓住根本定义:概率分布 F ( X ) F(X) F(X)累积概率,为 x ≤ X x≤X xX的累积概率和。
    • 考虑“概率微元” P ( X ≤ x ≤ X + h ) = F ( X + h ) − F ( X ) P(X≤x≤X+h)=F(X+h)-F(X) P(XxX+h)=F(X+h)F(X)。即随机变量 x x x落在 ( X , X + h ) (X,X+h) (X,X+h)区间上的概率。如果令 h → 0 h→0 h0,那么则会导向 lim ⁡ h → 0 P ( X ≤ x ≤ X + h ) = 0 \lim_{h \to 0}P(X≤x≤X+h)=0 limh0P(XxX+h)=0,即单点概率为0。该问题说明:单计连续随机变量每一点的概率 P ( x = x i ) P(x=x_i) P(x=xi)是有“问题”的,就像“单计连续金属杆每一点的质量”(质量密度与质量的关系)。
    • 由此,引入概率密度:随机变量落在无穷小区间内的平均概率 f ( x ) = lim ⁡ h → 0 F ( x + h ) − F ( x ) h = d F d x f(x)=\lim_{h \to 0}{\frac{F(x+h)-F(x)}{h}}=\frac{dF}{dx} f(x)=limh0hF(x+h)F(x)=dxdF
    • 综上: F ( X ) = ∫ − ∞ X f ( x ) d x F(X) = \int_{-\infty}^{X}f(x)dx F(X)=Xf(x)dx!!!
  • 概率密度函数的物理意义
    • 概率密度 = 概率 / “距离”: f ( x ) = P ( X ≤ x ≤ X + h ) d x f(x) = \frac{P(X≤x≤X+h)}{dx} f(x)=dxP(XxX+h)
    • 概率 P ( X ≤ x ≤ X + h ) = d N N P(X≤x≤X+h) = \frac{dN}{N} P(XxX+h)=NdN
    • 所以 f ( x ) = d N N ⋅ d x f(x) = \frac{dN}{N·dx} f(x)=NdxdN
    • 意义: 重要等式,在热学和统计力学中会多次用到, N N N为随机变量总数/粒子总数。

2、麦克斯韦速度分布推导

  • 注意:名为“分布”,实则速度的概率密度函数!!!
  • 推导前提:①分子平均动能方程( 1 2 m v 2 ‾ = 3 2 k T \frac{1}{2} m \overline{v^2} = \frac{3}{2} kT 21mv2=23kT, k k k玻尔兹曼常量)、②三轴速度 v x 、 v y 、 v z v_x、v_y、v_z vxvyvz均匀各向同性且独立同分布
  • 推导:
    • 根据②特性: F ( v x , v y , v z ) = F ( v 2 ) = F ( v x 2 + v y 2 + v z 2 ) = g ( v x 2 ) g ( v y 2 ) g ( v z 2 ) F(v_x,v_y,v_z)=F(v^2)=F(v_x^2+v_y^2+v_z^2)=g(v_x^2)g(v_y^2)g(v_z^2) F(vx,vy,vz)=F(v2)=F(vx2+vy2+vz2)=g(vx2)g(vy2)g(vz2)
    • 两边取对数: ln ⁡ [ F ( v x 2 + v y 2 + v z 2 ) ] = ln ⁡ [ ( ] g ( v x 2 ) ] + ln ⁡ [ g ( v y 2 ) ] + ln ⁡ [ g ( v z 2 ) ] \ln[F(v_x^2+v_y^2+v_z^2)]=\ln[(]g(v_x^2)]+\ln[g(v_y^2)]+\ln[g(v_z^2)] ln[F(vx2+vy2+vz2)]=ln[(]g(vx2)]+ln[g(vy2)]+ln[g(vz2)]
    • 两组求导:
      ∂ [ F ( v 2 ) ] ∂ v x 2 = ∂ [ F ( v 2 ) ] ∂ v 2 ⋅ ∂ v x 2 ∂ v 2 = ∂ [ F ( v 2 ) ] ∂ v 2 \frac{\partial [F(v^2)]}{\partial v_x^2}=\frac{\partial [F(v^2)]}{\partial v^2}·\frac{\partial v_x^2}{\partial v^2} = \frac{\partial [F(v^2)]}{\partial v^2} vx2[F(v2)]=v2[F(v2)]v2vx2=v2[F(v2)]
      ∂ [ F ( v 2 ) ] ∂ v y 2 = ∂ [ F ( v 2 ) ] ∂ v 2 ⋅ ∂ v y 2 ∂ v 2 = ∂ [ F ( v 2 ) ] ∂ v 2 \frac{\partial [F(v^2)]}{\partial v_y^2}=\frac{\partial [F(v^2)]}{\partial v^2}·\frac{\partial v_y^2}{\partial v^2} = \frac{\partial [F(v^2)]}{\partial v^2} vy2[F(v2)]=v2[F(v2)]v2vy2=v2[F(v2)]
      ∂ [ F ( v 2 ) ] ∂ v z 2 = ∂ [ F ( v 2 ) ] ∂ v 2 ⋅ ∂ v z 2 ∂ v 2 = ∂ [ F ( v 2 ) ] ∂ v 2 \frac{\partial [F(v^2)]}{\partial v_z^2}=\frac{\partial [F(v^2)]}{\partial v^2}·\frac{\partial v_z^2}{\partial v^2} = \frac{\partial [F(v^2)]}{\partial v^2} vz2[F(v2)]=v2[F(v2)]v2vz2=v2[F(v2)]
      ∂ [ F ( v 2 ) ] ∂ v x 2 = 1 g ( v x 2 ) ⋅ g ′ ( v x 2 ) \frac{\partial [F(v^2)]}{\partial v_x^2}=\frac{1}{g(v_x^2)}·g'(v_x^2) vx2[F(v2)]=g(vx2)1g(vx2)
      ∂ [ F ( v 2 ) ] ∂ v y 2 = 1 g ( v y 2 ) ⋅ g ′ ( v y 2 ) \frac{\partial [F(v^2)]}{\partial v_y^2}=\frac{1}{g(v_y^2)}·g'(v_y^2) vy2[F(v2)]=g(vy2)1g(vy2)
      ∂ [ F ( v 2 ) ] ∂ v z 2 = 1 g ( v z 2 ) ⋅ g ′ ( v z 2 ) \frac{\partial [F(v^2)]}{\partial v_z^2}=\frac{1}{g(v_z^2)}·g'(v_z^2) vz2[F(v2)]=g(vz2)1g(vz2)
      ⇒ 1 g ( v x 2 ) ⋅ g ′ ( v x 2 ) = 1 g ( v y 2 ) ⋅ g ′ ( v y 2 ) = 1 g ( v z 2 ) ⋅ g ′ ( v z 2 ) \Rightarrow \frac{1}{g(v_x^2)}·g'(v_x^2)=\frac{1}{g(v_y^2)}·g'(v_y^2)=\frac{1}{g(v_z^2)}·g'(v_z^2) g(vx2)1g(vx2)=g(vy2)1g(vy2)=g(vz2)1g(vz2)
    • 根据各向速度相互独立,或者说函数 1 g ( v 2 ) ⋅ g ′ ( v 2 ) \frac{1}{g(v^2)}·g'(v^2) g(v2)1g(v2)的不同取值结果相同,说明:
      1 g ( v 2 ) ⋅ g ′ ( v 2 ) = C o n s t \frac{1}{g(v^2)}·g'(v^2) = Const g(v2)1g(v2)=Const
      ⇒ 1 g ( v i 2 ) ⋅ g ′ ( v i 2 ) = B e − α v i 2 \Rightarrow \frac{1}{g(v_i^2)}·g'(v_i^2) = Be^{-\alpha v_i^2} g(vi2)1g(vi2)=Beαvi2
      ⇒ F ( v 2 ) = B 3 e − α ∑ v i 2 = B 3 e − α v 2 \Rightarrow F(v^2) = B^3e^{-\alpha \sum_{}{} v_i^2} = B^3e^{-\alpha v^2} F(v2)=B3eαvi2=B3eαv2
    • 要求 F ( v 2 ) F(v^2) F(v2)满足归一化条件:
      ∫ − ∞ + ∞ ∫ − ∞ + ∞ ∫ − ∞ + ∞ F ( v x , v y , v z ) d v x d v y d v z = 1 \int_{-\infty}^{+\infty} \int_{-\infty}^{+\infty}\int_{-\infty}^{+\infty} F(v_x,v_y,v_z)dv_xdv_ydv_z = 1 +++F(vx,vy,vz)dvxdvydvz=1
      ⇒ B = α π \Rightarrow B = \sqrt{\frac{\alpha}{\pi}} B=πα
    • 根据①分子的平均动能(需用到泊松积分 I = ∫ − ∞ + ∞ e − A x 2 d x = π A I = \int_{-\infty}^{+\infty}e^{-Ax^2}dx = \sqrt{\frac{\pi}{A}} I=+eAx2dx=Aπ
      ):
      v ‾ 2 = 3 k T m \overline v^2 = \frac{3kT}{m} v2=m3kT
      ⇒ ∫ − ∞ + ∞ ∫ − ∞ + ∞ ∫ − ∞ + ∞ v 2 F ( v x , v y , v z ) d v x d v y d v z = 3 k T m \Rightarrow \int_{-\infty}^{+\infty} \int_{-\infty}^{+\infty}\int_{-\infty}^{+\infty} v^2F(v_x,v_y,v_z)dv_xdv_ydv_z = \frac{3kT}{m} +++v2F(vx,vy,vz)dvxdvydvz=m3kT
      ⇒ ∫ − ∞ + ∞ ∫ − ∞ + ∞ ∫ − ∞ + ∞ v 2 B 3 e − α v 2 d v x d v y d v z = 3 k T m \Rightarrow \int_{-\infty}^{+\infty} \int_{-\infty}^{+\infty}\int_{-\infty}^{+\infty} v^2 B^3e^{-\alpha v^2}dv_xdv_ydv_z = \frac{3kT}{m} +++v2B3eαv2dvxdvydvz=m3kT
      ⇒ B 3 ∫ − ∞ + ∞ ∫ − ∞ + ∞ ∫ − ∞ + ∞ − ∂ e − α v 2 ∂ α d v x d v y d v z = 3 k T m \Rightarrow B^3 \int_{-\infty}^{+\infty} \int_{-\infty}^{+\infty} \int_{-\infty}^{+\infty} -\frac{\partial e^{-\alpha v^2}}{\partial \alpha} dv_xdv_ydv_z = \frac{3kT}{m} B3+++αeαv2dvxdvydvz=m3kT
      ⇒ − B 3 ∂ ∂ α [ ∫ − ∞ + ∞ ∫ − ∞ + ∞ ∫ − ∞ + ∞ − e − α v 2 d v x d v y d v z ] = 3 k T m \Rightarrow -B^3 \frac{\partial}{\partial \alpha}[\int_{-\infty}^{+\infty} \int_{-\infty}^{+\infty} \int_{-\infty}^{+\infty} -e^{-\alpha v^2} dv_xdv_ydv_z] = \frac{3kT}{m} B3α[+++eαv2dvxdvydvz]=m3kT
      ⇒ − B 3 ∂ ∂ α [ ( π α ) 3 2 ] = 3 k T m \Rightarrow -B^3 \frac{\partial}{\partial \alpha}[(\frac{\pi}{\alpha})^{\frac{3}{2}}] = \frac{3kT}{m} B3α[(απ)23]=m3kT
      ⇒ α = m 2 k T \Rightarrow \alpha = \frac{m}{2kT} α=2kTm
    • 麦克斯韦速度分布函数
      F ( v x , v y , v z ) = ( m 2 π k T ) 3 2 e − m 2 k T ( v x 2 + v y 2 + v z 2 ) F(v_x,v_y,v_z) = (\frac{m}{2 \pi kT})^{\frac{3}{2}} e^{-\frac{m}{2kT}(v_x^2+v_y^2+v_z^2)} F(vx,vy,vz)=(2πkTm)23e2kTm(vx2+vy2+vz2)

3、麦克斯韦速率分布推导

  • 逻辑推导:
    • 速度 ( v x , v y , v z ) (v_x,v_y,v_z) (vx,vy,vz)带方向,速率 v v v不带方向( v = ( v x 2 + v y 2 + v z 2 v = \sqrt{(v_x^2+v_y^2+v_z^2} v=(vx2+vy2+vz2
      )。
    • 置于三维坐标轴,考虑速度各向同性且独立同分布,故以原点为中心,以 v v v为半径的球壳为速度概率密度“等值面”,故该球壳的“体积”即为速率 v v v的概率密度。
    • 即: f ( v ) = 4 π v 2 F ( v x , v y , v z ) f(v)=4 \pi v^2 F(v_x,v_y,v_z) f(v)=4πv2F(vx,vy,vz) f ( v ) d v = 4 π v 2 F ( v x , v y , v z ) d v f(v)dv=4 \pi v^2 F(v_x,v_y,v_z)dv f(v)dv=4πv2F(vx,vy,vz)dv(之所以带微分 f ( v ) d v f(v)dv f(v)dv是因为 f ( v ) d v = d N N f(v)dv=\frac{dN}{N} f(v)dv=NdN表示概率)
    • 速率密度函数
      f ( v ) = 4 π v 2 ( m 2 π k T ) 3 2 e − m 2 k T ( v x 2 + v y 2 + v z 2 ) f(v) = 4 \pi v^2 (\frac{m}{2 \pi kT})^{\frac{3}{2}} e^{-\frac{m}{2kT}(v_x^2+v_y^2+v_z^2)} f(v)=4πv2(2πkTm)23e2kTm(vx2+vy2+vz2)
  • 初步分析:影响速率函数的热力学参量有 m 、 T m、T mT,即只有温度分子质量会影响速率函数。定性地,温度越高、分子质量越小,速率分布曲线越平缓(速率很大的概率增加,使得速率分布趋于平均)。

4、麦克斯韦速率分布的引申

  • 三种重要速率
    • 平均速率 v ‾ = 8 k T π m \overline v = \sqrt{\frac{8kT}{\pi m}} v=πm8kT
      v ‾ = ∫ 0 + ∞ v f ( v ) d v \overline v = \int_{0}^{+\infty} vf(v)dv v=0+vf(v)dv):用于计算分子平均自由程、气体碰壁数、气体分子碰撞概率
    • 最概然速率 v p = 2 k T m v_p = \sqrt{\frac{2kT}{m}} vp=m2kT
      ∂ f ( v ) ∂ v = 0 \frac{\partial f(v)}{\partial v} = 0 vf(v)=0):用于讨论不同速度分布曲线
    • 方均根速率 v 2 ‾ = 3 k T m \overline {v^2} = \frac{3kT}{m} v2=m3kT:用于计算分子平均动能
    • 三个速度的比例: v p : v ‾ : v 2 ‾ = 1 : 1.128 : 1.224 v_p:\overline v:\overline{v^2}=1:1.128:1.224 vp:v:v2=1:1.128:1.224

免责声明:本站所有文章内容,图片,视频等均是来源于用户投稿和互联网及文摘转载整编而成,不代表本站观点,不承担相关法律责任。其著作权各归其原作者或其出版社所有。如发现本站有涉嫌抄袭侵权/违法违规的内容,侵犯到您的权益,请在线联系站长,一经查实,本站将立刻删除。 本文来自网络,若有侵权,请联系删除,如若转载,请注明出处:https://haidsoft.com/142981.html

(0)
上一篇 2025-05-06 18:26
下一篇 2025-05-06 18:33

相关推荐

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注微信