应用多方过程玻尔兹曼分布律估算大气层总质量

1、应用多方过程的玻尔兹曼分布律估算大气层总质量（玉溪师范学院理学院物理系05级物理2班 云南 玉溪 653100）摘要： 在大多数的热学教材中只讨论等温情况下的玻尔兹曼分布律，但在研究大气过程时，若把大气过程看成是绝热过程或多方过程，这时等温过程的玻尔兹曼分布律不再适用。本论文通过对玻尔兹曼分布律的研究，把大气过程看成多方过程，在多方过程模型下推导出了多方过程玻尔兹曼分布律，最终应用多方过程玻尔兹曼分布律，对球周围大气层总质量进行精密估算。关键词： 多方过程 玻尔兹曼分布律 大气层总质量1.引言在物理学中，经常对某个物理量的值进行估算。如估算太阳质量、地球质量等。那如何对地球周围大气层的总质量进

2、行估算呢？对于这个问题，我们可以应用玻尔兹曼分布律对地球周围大气层的总质量进行估算1,2。若把大气过程看成绝热过程或多方过程3，这时就需要应用绝热过程或多方过程的玻尔兹曼分布律估算地球周围大气层的总质量。许多文献在估算地球周围大气层的总质量时，往往把大气过程看成等温过程，应用等温过程的玻尔兹曼分布律估算大气层的总质量4。在普通热力学教材中只给出了等温情况下的玻尔兹曼分布律，只能用于等温过程。若把大气过程看成绝热或多方过程，这时等温过程的玻尔兹曼分布律不在适用。必须建立绝热过程或多方过程气体模型5，在这个模型下推导出绝热过程和多方过程的玻尔兹曼分布律，并应用多方过程的玻尔兹曼分布律估算地球周围大

3、气层的总质量。本文包括五部分内容：第一部分分别对等温过程、绝热过程、多方过程的玻尔兹曼分布律进行推导；第二部分将等温过程、绝热过程、多方过程的玻尔兹曼分布律进行比较，分析得出把大气过程看成多方过程的科学性。第三部分应用多方过程的玻尔兹曼分布律对地球周围大气层的总质量进行精密估算。第四部分总结和讨论应用多方过程的玻尔兹曼分布律估算地球周围大气层的总质量的精密性和存在的不足。2.玻尔兹曼分布律的推导推导玻尔兹曼分布律有多种方式，可以从麦克斯韦分布律出发，推导玻尔兹曼分布律6,7,8，也可以应用大气压强与大气层重量之间存在的关系，推导玻尔兹曼分布律。2.1.等温过程玻尔兹曼分布律在等温过程中，由于温

4、度是一个恒定值,即不考虑温度的影响，而麦克斯韦分布律也不考虑温度的影响9,10，则等温过程的玻尔兹曼分布律可由麦克斯韦分布律推导得出。若大气处于平衡状态,当气体分子间的相互作用可以忽略时,麦克斯韦分布律为11： （1）把麦克斯韦分布律推广到分子在保守立场中的运动的情况810，在这种情形下，以总能量代替，这时分子不但限定在一定的速度区间内，而且它们的位置也限定在一定的坐标区间。当系统在力场中处于平衡状态时，坐标介于区间；内，同时速度介于；内的分子数为15： （2）上式对所有可能的速度积分，考虑到麦克斯韦分布函数应满足的归一化条件： （3）则（2）式变为：这里的表示分布在；内具有各种速度的分子总数

5、。再以除上式，则得出分布在坐标区间；内单位体积内的分子数,即分子数密度为： (4) 在应用麦克斯韦分布律推导玻尔兹曼分布律的过程中，可以看出式子中的始终是一个常数，即温度是一个定值。(4)式即为等温过程的玻尔兹曼分布律。 由理想气体状态方程： （5）把（5）代入（4）式得： （6）上式给出了压强随高度的变化关系，称为等气压公式。2.2.绝热过程、多方过程玻尔兹曼分布律在推导玻尔兹曼分布律的过程中，无论把大气模型看成是绝热过程或是多方过程气体模型，地球周围大气层的温度T都不能看成恒定值，温度应该看成一个高度Z的函数12,13，这时就必须考虑温度对玻尔兹曼分布律的影响。在麦克斯韦分布律中并不考虑温

6、度的影响，所以绝热过程和多方过程的玻尔兹曼分布律，不能应用麦克斯韦分布律推导得出，必须考虑用大气压强与大气层重量之间存在的关系推导绝热过程和多方过程玻尔兹曼分布律14,15。根据大气压强与大气层重量之间存在的关系，建立如图（1）所示的模型，我们认为大气过程是一个准静态过程，大气为理想气体，设平衡气体的压强随高度变化的函数关系为，在气体中取一柱气体，其上下端面水平面积为，柱体的高为，气柱上下端面所受压力分别为和，二者之差与气柱所受的重力平衡16。 图（1）在推导绝热过程和多方过程的玻尔兹曼分布律的过程中，要考虑温度的影响，就必须要先求出温度随高度的变关系即温度梯度。2.2.1 温度随高度的变化关

7、系设气体的分子数密度为N，压强为P,高度为z，温度为T，由大气压强与大气层重量之间的关系： （7）由上式得： （8）把空气看为理想气体,大气过程看为是多方过程,设多方指数为（即为绝热，所以可以把绝热看成多方的特殊情况）。对(8)式用复合函数求导可得: （9）由热力学可知,多方过程方程为： （10）对(10)式求偏导数可得: (11)将(11)式代入(9)式可得: (12)将(5)式理想气体气态方程代入(12)式得: (13)上式说明在多方过程中,大气的温度梯度与高度无关,为一个常数17,18。对(13)式进行积分并考虑初始条件（ ）得： (14)上式可以看出，温度是高度的一次函数，(14)式给

8、出了温度随高度的变化关系。2.2.2绝热过程、多方过程玻尔兹曼分布律的推导现在在图（1）建立的模型和温度与高度的变化关系的基础上推导绝热过程和多方过程的玻尔兹曼分布律。由（5）式对Z求导得： (15)把（8）式（13）式（14）式代入（15）式得： (16)对(16)式分离变量并积分得： (17)有Z=0时得 (18)把（18）式代入（17）式得 (19)由（19）式得出分子数密度为 (20)上式即为多方过程的玻尔兹曼分布律.它给出了分子数密度随高度的变化关系。式中，N,P分别为高度z=0处和高度为z处的分子数密度和压强, 为高度z=0处的温度,n为多方指数,为空气的摩尔质量,R为普适气体常数

9、。 由(6)式理想气体状态方程得 (21)把(21)式和(14)式代入(20)式得 (22)上式给出了多方过程（为绝热过程）和高度的变化关系，即多方过程（绝热过程）的压强公式。3.等温模型、绝热模型、多方模型的比较我们可以把大气过程看作等温过程，也可以将大气过程看作绝热过程或多方过程，究竟建立何种大气模型比较接近真实情况，下面将三种模型进行比较。3.1.等温过程与绝热过程、多方过程比较若我们把大气过程看成等温过程即可得到等温过程的玻尔兹曼分布律。若我们把大气过程看成多方过程即可得到多方过程的玻尔兹曼分布律。然而把大气过程看成是等温过程不够精密。因为在实际情况下，大气层的温度与高度有关, 大气层

10、的温度是一个高度的函数，温度随高度增加而降低。温度随高度的变化（见表一）。表 一 高度与温度的关系12 高度Km温度K(实测值)0288.1451281.6512275.1543268.6594262.1665255.6766249.1877242.7008236.2159229.73310223.252由表一可以看出，大气温度随高度的升高而降低，而且变化幅度还较大，若把大气过程看成等温过程显然较为初略。若把多方指取为，即等温的情况。把代入(16)式得： （23）对(23)式分离变量，积分并考虑初始条件（时）得：上式即回到等温情况下的玻尔兹曼分布律。 我们将,代入上式得: （24）分别将，,等

11、代入（24）式得出等温过程分子数密度和高度的关系（见表二）。表 二 分子数密度与高度的关系（等温过程）高度（Km） 实测值（）计算值（）误差02.54672.5467012.31062.2623-2.3%22.09222.0094 -4.4%31.89001.7848-5.5.%41.70311.5854-6.8%51.53171.4083-8.2.%61.37211.2509-9.1%81.09290.9870-10.2%100.85950.7787-10.4%从表二中可以看出，等温情况下计算出的分子数密度随高度的增加误差也逐渐增大。当时，误差达到了10%。等温情况下的玻尔兹曼分布律存在较大

12、误差，主要是由于没有考虑温度的影响。为了更精密的描述大气过程，我们在考虑温度的影响的情况下，推导出了绝热过程和多方过程的玻尔兹曼分布律。3.2.绝热过程与多方过程的比较将(20)式与绝热过程的玻尔兹曼分布律作一比较，当（为绝热系数）时，（20）式变为： （25）代入理想气体的状态方程得： （26）(25)式即为绝热过程时的玻尔兹曼分布律。（26）式为绝热模型压强和高度的关系。因此我们可以将绝热过程看成是多方过程的一种特殊情况。若我们把大气过程看成绝热过程，这也不够精密，因为大气在流动的过程中有一定的热交换，以及太阳及地面对大气热量的供给，大气过程并不是一个严格的绝热系统，因此在绝热模型下推倒出

13、的玻尔兹曼分布律必然也不够精密。若把大气过程看成一个介于等温过程和绝热过程之间的一个多方过程则比较合理，这时多方指数n应该是在1和（）之间的某一个值1。多方指数n可以由温度梯度来确定。由（13）式 可以知道，在多方过程中，大气的温度梯度与高度无关为一个常数。这个常数可以从温度与高度的关系中（表一）求出。从表一可得在0到10Km之间，温度梯度的平均值为6.489K/mol.将这个值代入（13）式，得出n的值为1.235，我们求出的n值的确是介于1和1.4之间。 由得：我们将，代入到（25）中得到绝热过程的分子数密度和高度的关系式： （27）分别将z=1km，2km，等代入（27）式中得到表三。我

14、们将，代入到（20）中得： （28）同样分别将z=1km，2km，等代入（28）式中得到表四。表 三 分子数密度与高度的关系（绝热过程）高度（Km）实测值（）计算值（）误差02.54672.5467012.31062.33691.1%22.09222.13792.2%31.89001.94943.1%41.70311.77124.0%51.53171.60324.7%61.37211.44515.3%81.09291.15805.9%100.85950.90816.2%表 四 分子数密度与高度的关系（多方过程）高度（Km）实测值（）计算值（）误差02.54672.5467012.31062.3

15、1120.02%22.09222.09270.03%31.89001.8897-0.01%41.70311.7034-0.02%51.53171.5309-0.03%61.37211.3708-0.04%81.09291.0924-0.05%100.85950.8587-0.09%通过比较等温过程、绝热过程、多方过程的玻尔兹曼分布律的误差项，可以看出多方过程的误差相对较小，估算出的分子数密度比较接近真实值，所以用多方过程来描述大气过程也比较科学，由多方过程得出的玻尔兹曼分布律也比较精密。4.地球周围大气层总质量估算通过以上的分析可知，多方过程描述大气过程比等温过程、绝热过程更为精密。因此我们应

16、用多方过程的玻尔兹曼分布律来估算地球周围大气层的总质量。在估算地球周围大气层总质量时,我们把大气过程看为多方过程，考虑一个理想化的大气分子模型，把空气视为理想气体。在这个模型中,地球为一半径为的刚性球体,球外引力场中重力加速度为常数g。地球表面以上共有个质量为m的大气分子。根据(20)式多方过程的玻尔兹曼分布律,分子密度在重力场中随高度的分布规律: (29)式中为地球表面（）所在处大气分子数密度。于是地球上空大气层总质量可通过下式计算,于是得到地球周围大气层总质量 （30）有（28）式可以看出N随Z的增大而减小，并不收敛，这是就必须考虑其物理意义，N最小值只能为零，则由（29）式N=0得： （

17、31）,由得：，,，代入（31）式得：,把(29)式代入（30）式并考虑积分上下限，得: （32） 代入(32)式得: （33）计算上式的积分得：通过上式计算得到地球周围大气层的总质量为： （34）上式即为应用多方过程玻尔兹曼分布律估算出的地球周围大气层的总质量。5.结论本论文基于等温模型、绝热模型、多方模型三种大气过程模型，推导出了等温过程、绝热过程、多方过程的玻尔兹曼分布律。通过比较等温过程、绝热过程、多方过程的玻尔兹曼分布律，得出了应用多方过程玻尔兹曼分布律估算地球周围大气层总质量更为精密的结论。这主要是由于在考虑大气过程时，多方过程和绝热过程的玻尔兹曼分布律没有把温度作为恒定值处理，考

18、虑了温度梯度对分子数密度的影响，因此比等温过程玻尔兹曼分布律更精密。然而把大气过程看成多方过程，估算出的分子数密度比较接近测量值，与等温过程、绝热过程相比误差较小，因此用多方过程的玻尔兹曼分布律估算出的地球周围大气层的总质量比用等温过程和绝热过程的玻尔兹曼分布律估算出的地球周围大气层总质量更接近其真实值。应该指出，虽然多方过程玻尔兹曼分布律估算出的地球周围大气层的总质量比用等温过程和绝热过程的玻尔兹曼分布律估算出的地球周围大气层总质量精密。但是应用多方过程玻尔兹曼分布律计算地球周围大气层总质量时也只是近似的。因为:第一,(34)式的最终求得用到了(20)式,而（20）式只是对地球以外实际大气的

19、一种近似。第二,在（32)式的积分区间中，地球引力场中的重力加速度g是个变量,而我们在考虑理想化的大气分子模型中却把它作为一个常数,具体计算时取了。第三,大气实际上是含有多种成分的混合气体,而我们在估算大气总质量时却是以空气平均分子量来计算的。第四,实际的大气并不是理想气体，而我们在推导的过程中却把它视为理想气体分子来处理。第五,从（33）式到（34）式代入数值时做了近似的处理。所以若要想更为精密的估算地球周围大气层的总质量还必须做进一步的努力和研究，采用其他方法或模型对大气层质量进行更精密的估算。致 谢：本文在写作过程中得到了周菁老师的倾力指导和大力支持，在此表示衷心的感谢！参考文献：1 李

20、强.玻耳兹曼分布律的推广和应用J.达县师范高等专科学校学报,2004 14(5):29-302 KEVIN E.TRENBERTH , LESLEY SMITH. The Mass of the Atmosphere A Constraint on Global Analyses. Manuscript received 17 November 2003,in final form 29 June 20043 姚云珍，程俊辉.绝热模型和多方模型下大气压强与高度的关系J. 扬州教育学院学报,2000,3 ：68-70 4 郑玉民,凌瑞良.地球周围大气层总质量的估算J.常熟高专学报,2003,17

21、(4)：39-415 方小兵.大气压与高度关系的再讨论J.安康师专学报,2005,17(2）:95-97 6 Yan Guo .The Vlasov-Maxwell-Boltzmann system near MaxwelliansJ. Invent math,2003,10:593-6107 李椿,章力源,钱尚武.热学M.高等教育出版社,987,91-93 8 程守洙,江之永.普通物理学第一册.M.高等教育出版社,1998,6：81-829 Eric Carlen·Ester Gabetta·Eugenio RegazziniJ. On the Rate of Explo

22、sion for Infinite Energy Solutions of the Spatially Homogeneous Boltzmann Equation, J Stat Phys 2007,129:69972310 张三慧.麦克斯韦速度分布定律的适用范围J.大学物理, 2001,20(8)：10-1211 赵凯华,罗蔚茵.新概念物理教程.热学M.高等教育出版社,1998：92-10412 陈竞余.大气压强与高度的精密公式J.大理师专学报(自然科学版）,1994,2：46-4913 Jian Zheng .Spectrum of transition radiation from h

23、ot electrons generated in ultra-intense laser plasma interactionJ. PHYSICS OF PLASMAS, 2002, 9(8):3610-361614 Viviane Pierrard .The Boltzmann Equation to Study the Escape of LightAtoms from Planetary Atmospheres. J Stat Phys 2007,129:69972315 王长明. 玻耳兹曼分布律的推导J.白城师范高等专科学校学报, 2002,16(3)：65-66 16 王力新.关于

24、气体压强在大气中与“大气重力压强”的讨论J.高师理科学刊, 2001,21(1)：69-7117 吴浩然.大气压强和高度的关系J. 张家口 师专学报（自然科学版), 1997,2：72-7318刘成周.温度对重力场中粒子数按高度分布律的影响J.滨州师专学报， 2003,19(4)：43-45Application of multi-process Boltzmann Distribution Law to Estimate the Total Mass of the AtmosphereMing-Hu Li(Department of physics, Yuxi normal university, Yuxi 653100)Dir