自学脚手架——《热学》 by 李椿(第一,二,三,四,五章)
文章目录
- 绪论
- 第一章 温度
- 第二章 气体分子动理论的基本概念
- 第三章 气体分子热运动速率和能量的统计分布律
- 第四章 气体内的输运过程
- 第五章 热力学第一定律(一)
绪论
第一章 温度
- 平衡态需要同时达成力学平衡,热学平衡,化学平衡。
第二章 气体分子动理论的基本概念
- 雷德利克-邝方程(Redlich-Kwong equation)
在现有的含有两个参量的非理想气体物态方程中,精确度最高的是雷德利克-邝方程:
{p+a/[T2ν(ν+b)]}⋅(ν−b)=RT\{p+a/[T^{2}\nu(\nu+b)]\}\cdot (\nu-b)=RT {p+a/[T2ν(ν+b)]}⋅(ν−b)=RT
此方程中的ppp、TTT、ν\nuν分别是气体的压强、热力学温度、摩尔体积,RRR是普适气体常量,aaa和bbb是两个雷德里克-邝常量。
常量可以利用临界温度和临界压强计算。
非常短的时间内的碰撞,动量总是守恒的,因为分子之间的作用力非常大,且作用过程时间很短;传热的过程中,传热几乎没有发生,此时可以看做绝热过程(计算声速的时候会用到)。
经验表明:由量纲分析所得结果中的无量纲量AAA的数量级往往为1.
量纲分析方法一般有瑞利法,Π\PiΠ定理,白金汉法。
注意区分相对质心的无规则运动和质心的定向运动,基于前者才能获得系统的热力学统计规律。
平常的分子运动中有非常高能量的分子存在,但是其告诉的运动(很高的能量)并不能使得电离、核反应等现象发生,就说明了能量的量子化的部分特性。
分子直径与热力学温度的关系
两个分子之间的弹性碰撞是一个散射过程。在此过程中,两个分子中心之间的距离rrr所达到的最小值,可以视为分子直径。因此,分子直径ddd应该与分子的平均平动动能ε\varepsilonε和分子间斥力常量λ\lambdaλ有关。设热力学温度为TTT,分子间斥力为fff,则由量纲分析可知(由于平衡位置的尺度处的斥力远大于引力,所以只使用伦纳德-琼斯势的第一项):
dimλ=dim(frs)=Lsdimfdimε=Ldimfdimd=L=[(Lsdimf)/(Ldimf)]1/(s−1)=(dimλ/dimε)1/(s−1)=dim{[λ/(kT)]1/(s−1)}\begin{aligned} \mathrm{dim}\ \lambda &=\mathrm{dim}\ (fr^{s})=L^{s}\mathrm{dim}\ f\\ \mathrm{dim}\ \varepsilon &=L\mathrm{dim}\ f\\ \mathrm{dim}\ d &=L\\ &=[(L^{s}\mathrm{dim}\ f)/(L\mathrm{dim}\ f)]^{1/(s-1)}\\ &=(\mathrm{dim}\ \lambda/\mathrm{dim}\ \varepsilon)^{1/(s-1)}\\ &=\mathrm{dim}\ \{[\lambda/(kT)]^{1/(s-1)}\} \end{aligned} dim λdim εdim d=dim (frs)=Lsdim f=Ldim f=L=[(Lsdim f)/(Ldim f)]1/(s−1)=(dim λ/dim ε)1/(s−1)=dim {[λ/(kT)]1/(s−1)}
由此可得:
d∝1/T1/(s−1)d\propto 1/T^{1/(s-1)} d∝1/T1/(s−1)
因s=9∼15s=9\sim15s=9∼15,s−1=8∼14s-1=8\sim 14s−1=8∼14,故1/(s−1)=(1/8)∼(1/14)1/(s-1)=(1/8)\sim(1/14)1/(s−1)=(1/8)∼(1/14),其平均值约为1/10=0.11/10=0.11/10=0.1。于是分子直径ddd与热力学温度TTT之间的关系可近似为:
d∝1/T1/10=1/T0.1d\propto 1/T^{1/10}=1/T^{0.1} d∝1/T1/10=1/T0.1
由此可见,分子直径ddd随着热力学温度TTT的增大而略微减小。通常的温度变化并不足以引起分子直径的很大变化,所以平常使用中可以把分子直径看做常量。
其中用到了伦纳德-琼斯势:
f(r)=λrs−μrtf(r)=\frac{\lambda}{r^{s}}-\frac{\mu}{r^{t}} f(r)=rsλ−rtμ
其中sss取值在9∼159\sim 159∼15,ttt取值在4∼74\sim 74∼7,前者通常取131313,后者通常取777。第一项为斥力,第二项为引力。其最初是米氏模型中提出的,后来由伦纳德-琼斯发展。
- 势能曲线与热膨胀
截取分子之间势能曲线的下半截,然后画出分子随温度变化的振动能级,不同能级下的振动中心位置(平衡位置)连成一条线,则这条线代表着物体随着温度的变化造成体积膨胀或缩小的过程,而这个过程在温度小范围内是一个线性的,线性代表了中心位置连线是一条近似直线,且其是有一定斜率的。如果平衡位置连线垂直于xxx轴,即曲线为对称的,则不会发生热膨胀。所以温度趋近于0的时候,热膨胀将会消失,实验中也验证了这一点。
范德瓦尔斯方程中的aaa能正能负?因为粒子碰壁之前即受到内部分子的吸引力,也受到壁上分子的吸引(当然壁上的分子数不确定多少),则如何?
范德瓦尔斯方程算出的压强可能大于理想气体方程算出的压强,也可能小于。
计算压强:
将相同速度的分子分为一组,即vi‾(vxi,vxi,vxi)\overline{v_{i}}(v_{xi},v_{xi},v_{xi})vi(vxi,vxi,vxi),然后撞在器壁上的单位面积为dA\mathrm{d}AdA,则单位时间内打到dA\mathrm{d}AdA的相同速度的分子形成一个斜的柱体,如下图所示:
如图所示,柱体体积为vixdtdAv_{ix}\mathrm{d}t\mathrm{d}AvixdtdA,柱体中速度为vi‾\overline{v_{i}}vi的分子数密度为ni,vx>0=12nin_{i,v_{x}>0}=\frac{1}{2}n_{i}ni,vx>0=21ni,单个分子传递的动量为2mvix2mv_{ix}2mvix,则速度为vi‾\overline{v_{i}}vi的分子传递的总动量为:
2mvix⋅12ni⋅vixdtdA2mv_{ix}\cdot \frac{1}{2}n_{i}\cdot v_{ix}\mathrm{d}t\mathrm{d}A 2mvix⋅21ni⋅vixdtdA
则不同速度的分子传递的总的动量为:
∑i2mvix⋅12ni⋅vixdtdA=∑imvix2nidtdA\sum_{i}2mv_{ix}\cdot \frac{1}{2}n_{i}\cdot v_{ix}\mathrm{d}t\mathrm{d}A=\sum_{i}mv_{ix}^{2}n_{i} \mathrm{d}t\mathrm{d}A i∑2mvix⋅21ni⋅vixdtdA=i∑mvix2nidtdA
则压强ppp为:
p=∑imvix2nidtdAdtdA=∑imvix2ni=mnvx2‾=13mnv2‾p=\frac{\sum_{i}mv_{ix}^{2}n_{i}\mathrm{d}t\mathrm{d}A}{\mathrm{d}t\mathrm{d}A} =\sum_{i}mv_{ix}^{2}n_{i}=mn\overline{v_{x}^{2}}=\frac{1}{3}mn\overline{v^{2}} p=dtdA∑imvix2nidtdA=i∑mvix2ni=mnvx2=31mnv2
- 估算碰壁数
与画圆柱不同的是,当我们直接用vˉ\bar{v}vˉ做量纲分析得出碰壁数的时候,不能再使用圆柱了,因为不同方向相同速率vvv的到达时间是不同的,由于分子是均匀分布的,所以最终打到底面dA\mathrm{d}AdA的分子可以看做是来自不同球面的分子。
根据量纲分析可得:
J∝nvˉJ\propto n\bar{v} J∝nvˉ
其中速度在各个方向上都有分布,在球面dS\mathrm{d}SdS上的某一个速度下的分子均匀分布,此时其中分子对底面dA\mathrm{d}AdA的立体角的面积为dA′=cosθdA\mathrm{d}A'=\cos \theta \mathrm{d}AdA′=cosθdA,由于需要在dS\mathrm{d}SdS上积分,所以可以将系数cosθ\cos\thetacosθ乘到dS\mathrm{d}SdS上,即化为平行于器壁的面积cosθdS\cos \theta\mathrm{d}ScosθdS,其就是dS\mathrm{d}SdS在器壁上的投影。如果球面叠加为柱体,则结果就是上式,但是此时应该为球面,整个球面的投影面积为πr2\pi r^2πr2,而球面的面积为2πr22\pi r^22πr2,所以实际球面时候粒子数应该为柱体时的12\frac{1}{2}21。再考虑到有一半的分子朝着与器壁方向的反方向飞走,所以占比为12\frac{1}{2}21,所以碰壁数可以重新整理为:
J=14nvˉJ=\frac{1}{4}n\bar{v} J=41nvˉ
第三章 气体分子热运动速率和能量的统计分布律
- 麦克斯韦速率分布的约化形式
令x=v/vpx=v/v_{p}x=v/vp,其中vp=2kTmv_{p}=\sqrt{\frac{2kT}{m}}vp=m2kT,将新变量xxx化为无量纲的形式,则麦克斯韦速率分布可以重新表示为:
f(v)dv=4π−1/2[m/(2kT)]3/2⋅e−mv2/(2kT)v2dv=4π−1/2vp−3⋅e−v2/vp2v2dv=4π−1/2x2⋅e−x2dx=F(x)dx\begin{aligned} f(v)\mathrm{d}v &=4\pi^{-1/2}[m/(2kT)]^{3/2}\cdot e^{-mv^{2}/(2kT)}v^{2}\mathrm{d}v\\ &=4\pi^{-1/2}v_{p}^{-3}\cdot e^{-v^{2}/v_{p}^{2}}v^{2}\mathrm{d}v\\ &=4\pi^{-1/2}x^{2}\cdot e^{-x^{2}}\mathrm{d}x\\ &=F(x)\mathrm{d}x \end{aligned} f(v)dv=4π−1/2[m/(2kT)]3/2⋅e−mv2/(2kT)v2dv=4π−1/2vp−3⋅e−v2/vp2v2dv=4π−1/2x2⋅e−x2dx=F(x)dx
约化形式处理和记忆起来都比较方便。
- 麦克斯韦速度分布的约化形式
令x=vx/vpx=v_{x}/v_{p}x=vx/vp,其中vp=2kTmv_{p}=\sqrt{\frac{2kT}{m}}vp=m2kT,则麦克斯韦速度分布可以重新表示为:
Fvx(x)dx=1πe−x2dxF_{v_{x}}(x)\mathrm{d}x=\frac{1}{\sqrt{\pi}}e^{-x^{2}}\mathrm{d}x Fvx(x)dx=π1e−x2dx
- 平均速率与分布的质心之间的关系
如图所示,求平均速率可以转化为求分布的质心:
vˉ=∫vdN∫dN\bar{v}=\frac{\int v \mathrm{d}N}{\int \mathrm{d}N} vˉ=∫dN∫vdN
xc=∫xdm∫dmx_{c}=\frac{\int x\mathrm{d}m}{\int \mathrm{d}m} xc=∫dm∫xdm
使用麦克斯韦速度分布律推倒碰壁数的时候,需要注意,其中的以dA\mathrm{d}AdA为底的柱体,其实是所有柱体的总和,因为此时只限定了柱体的高度vxdtv_{x}\mathrm{d}tvxdt,和柱体的体积vxdtdAv_{x}\mathrm{d}t\mathrm{d}AvxdtdA,还有速度vxv_{x}vx,但是其他两个方向的速度是未知的,所以这个柱体有很多斜着的柱体,不管斜着的柱体倾斜角度多少,最终的体积还是vxdtdAv_{x}\mathrm{d}t\mathrm{d}AvxdtdA,由于所有不同vxv_{x}vx和vyv_{y}vy的分子数密度为nin_{i}ni,当我们把所有不同vxv_{x}vx和vyv_{y}vy的体积折合成同一个(平移对称性)体积vxdtdAv_{x}\mathrm{d}t\mathrm{d}AvxdtdA,最终的分子数密度即为n=∑nin=\sum n_{i}n=∑ni。
麦克斯韦发射速率分布(分子射线分布)
当壁上开一小孔,考虑腔体内的分子射出小孔的情形(分子束发射等),则需要考虑单位时间单位面积下碰壁数J=14nvˉJ=\frac{1}{4}n\bar{v}J=41nvˉ按照速率vvv的分布:
J(v)dvdtdA=J0∼∞(ΔN0∼vN)=14nvˉ∫0vvf(v)dvvˉdtdA=14nvf(v)dvdtdA\begin{aligned} J(v)\mathrm{d}v\mathrm{d}t\mathrm{d}A &=J_{0\sim\infty}\left(\frac{\Delta N_{0\sim v}}{N}\right)\\ &=\frac{1}{4}n\bar{v}\int_{0}^{v}\frac{vf(v)\mathrm{d}v}{\bar{v}}\mathrm{d}t\mathrm{d}A \\ &=\frac{1}{4}nvf(v)\mathrm{d}v\mathrm{d}t\mathrm{d}A \end{aligned} J(v)dvdtdA=J0∼∞(NΔN0∼v)=41nvˉ∫0vvˉvf(v)dvdtdA=41nvf(v)dvdtdA
即将vˉ=∫vf(v)dv\bar{v}=\int vf(v)\mathrm{d}vvˉ=∫vf(v)dv,此分布可以用于核元素的扩散法分离,即利用铀同位素的分子质量不同导致扩散(穿越小孔)速度不同。
如果此时考虑的不是碰壁数,而是单个分子在某个速率区间的概率,则可将上式除以总的碰壁数即可:
fB=J(v)dvdtdAJdtdA=vf(v)vˉf_{B}=\frac{J(v)\mathrm{d}v\mathrm{d}t\mathrm{d}A}{J\mathrm{d}t\mathrm{d}A}=\frac{vf(v)}{\bar{v}} fB=JdtdAJ(v)dvdtdA=vˉvf(v)
根据这个新的分布可以计算出分子出射的平均速率vB‾\overline{v_{B}}vB为:
vB‾=∫vfBdv=9πkT8m\overline{v_{B}}=\int vf_{B}\mathrm{d}v=\sqrt{\frac{9\pi kT}{8m}} vB=∫vfBdv=8m9πkT
而方均根速率vB2‾\overline{v_{B}^{2}}vB2为:
vB2‾=4kTm\overline{v_{B}^{2}}=\sqrt{\frac{4kT}{m}} vB2=m4kT
课本描述验证麦克斯韦速度分布的密勒-库士实验中,需要验证的分布应该为v4e−mv22kTv^{4}e^{-\frac{mv^{2}}{2kT}}v4e−2kTmv2,书中的表述不严谨,需要把v=ωφv=\frac{\omega}{\varphi}v=φω进行微分。
考普定律
考普定律为经验定律,晶体化合物在温度为298K298K298K时的定压摩尔分子热容Cp,298C_{p,298}Cp,298约等于构成此化合物的所有原子的定压摩尔原子热容CiC_{i}Ci的总和,即:
Cp,298=∑NiCiC_{p,298}=\sum N_{i}C_{i} Cp,298=∑NiCi
式中的NiN_{i}Ni为此化合物中所含有的第iii种原子的个数。
- 碰撞数的应用
由于不同质量的气体分子的碰撞数是不同的,利用这一点,可以分离浓缩铀,这种方法被称为气体扩散(Gaseous diffusion)。
气体扩散是一种通过迫使气态六氟化铀(UF 6 ) 通过半透膜来生产浓缩铀的技术。这会在含有铀 235 ( 235 U) 和铀 238 ( 238 U)的分子之间产生轻微的分离。通过使用多级的大级联,可以实现高分离度。这是第一个能够生产工业上有用数量的浓缩铀的工艺,但现在被认为已经过时,已被更有效的气体离心工艺所取代。
图 气体扩散使用半透膜分离浓缩铀
从扩散的角度说其相应的定律为格雷厄姆定律,该定律指出气体的扩散速率与其分子质量的平方根成反比。例如,在带有包含两种气体混合物的半透膜的盒子中,较轻的分子会比较重的分子更快地从容器中排出。离开容器的气体在一定程度上富含较轻的分子,而残余气体则在一定程度上耗尽。通过气体扩散进行富集过程的单个容器称为扩散器。
- 利用麦克斯韦速度分布律求每秒碰到单位面积器壁上的气体分子数
针对书中的推导过程的常见问题:
推导过程中,怎么说明 vvv 方向气柱内的气体一定可以到达面积 dA\mathrm{d} AdA 上,怎么保证中途不会与其他速度的分子发生碰撞。如果发生碰撞,除非是沿速度方向的碰撞才不会有影响,但是这种可能性并不高啊,为什么就一定能到达气壁呢?
为何只考虑一个正气柱内的气体?
下面是对上述问题的分析:
书中气柱的物理图像与真实发生的情况是有不同的,前者是后者在统计上的等效处理。考虑柱体内的碰撞,首先就要注意自由程的问题,注意到推导过程中柱体的高度是 vxdtv_{x}\mathrm{d}tvxdt,此处 dt\mathrm{d}tdt 时间非常短,趋于 0,这会导致柱体的高度非常小,也就是在表面非常薄的一层,之所以会发生这种情况,是因为我们求的是强度量,并且其在时间上保持不变,也就是说非常小的时间 dt\mathrm{d}tdt 和很长的时间内求得的压强应该是相同的。所以说在此非常矮的柱体内,其柱体高度远远小于自由程的尺度,也就是说碰撞发生的概率极低。
另外推导中还有一个假设,就是碰撞到 dA\mathrm{d}AdA 上的分子并不一定只是一个气柱内的,还有可能是其他方向气柱内的,也就是按照相同 vxv_{x}vx,但不同的 vyv_{y}vy,vzv_{z}vz 去画气柱的时候,会出现无数个斜气柱,现在我们只考虑了一个正气柱(与斜气体柱)的原因在于,我们把所有气柱的水平分量忽略了,只考虑 vxv_{x}vx,这样做的依据是我们只关心 xxx 方向上的气体压强,水平方向上的情况与此无关。这个等效的过程中,我们把所有符合vx∼vx+dvxv_{x}\sim v_{x}+\mathrm{d}v_{x}vx∼vx+dvx 的分子都挪到了一个正气柱中,根据整个平面的平移不变性,正气柱的气体分子数密度也应该为 nnn,所以说这个图像也是与真实情况有差别的。
!!(待补充图片)
第四章 气体内的输运过程
- 分子之间的相对平均速率
在热学中,通过建立单个或者两个对象之间的关系,然后拓展到平均大量求和量,是一个化难为简的重要过程。
首先我们把大量气体分子遇到的问题化简为两个气体分子的问题。两个气体分子的相对速率为:
vr12‾=∣v1−v2∣\overline{v_{r12}}=|\bm{v_{1}}-\bm{v_{2}}| vr12=∣v1−v2∣
大量气体分子的平均相对速率vr‾\overline{v_{r}}vr定义为:
vr‾=1N(N−1)∑i,j,i≠jN∣vi−vj∣\overline{v_{r}}=\frac{1}{N(N-1)}\sum_{i,j,i\neq j}^{N}|\bm{v_{i}}-\bm{v_{j}}| vr=N(N−1)1i,j,i=j∑N∣vi−vj∣
注意iii和jjj是分子的编号,而非速度的区间。当分子的个数N→∞N\rightarrow\inftyN→∞时,1N(N−1)≈1N2\frac{1}{N(N-1)}\approx \frac{1}{N^{2}}N(N−1)1≈N21,上式可以约等为:
vr‾=1N2∑i,jN∣vi−vj∣\overline{v_{r}}=\frac{1}{N^{2}}\sum_{i,j}^{N}|\bm{v_{i}}-\bm{v_{j}}| vr=N21i,j∑N∣vi−vj∣
此时这个公式的意思就是,在一团气体系统中随机抽取一个粒子,然后再抽取一个粒子,本来第二次不能再抽到自己,但是由这个公式,我们可以假想到此时还能抽到自己,然后如果抽到自己的时候,后方的相对速度为0,但是这一项并不影响N→∞N\rightarrow\inftyN→∞时vrˉ\bar{v_{r}}vrˉ的值。所以前后两次的可能抽取次数都是NNN,即每个粒子被抽取到的概率是1N\frac{1}{N}N1。
其实根据上述应该得出这样的结论:两个分子的抽取可以看做是在两个相同的独立的系统之中抽取。这个概念是极为重要的。注意到上述的1/N1/N1/N为分子编号下的概率,有了这样一个独立的前提条件,我们就可以把这个编号换到速度区间上,即给相同一类速度的分子进行编号,即把在vi∼vi+dv\bm{v_{i}}\sim \bm{v_{i}}+\mathrm{d}\bm{v}vi∼vi+dv之间的分子编号为vi\bm{v_{i}}vi(将速度空间离散化,最后再连续化,这样有助于我们理解这个物理过程,这是极其重要的!),然后把在vj∼vj+dv\bm{v_{j}}\sim \bm{v_{j}}+\mathrm{d}vvj∼vj+dv之间的分子编号为vj\bm{v_{j}}vj。然后根据这个我们就可以建立一个采样过程,即现在一个独立系统中取一个分子,其速度为vi\bm{v_{i}}vi,然后概率即为f(vi)dvf(\bm{v_{i}})\mathrm{d}\bm{v}f(vi)dv,然后再从另一个独立系统中取一个分子,其速度为vj\bm{v_{j}}vj,然后概率即为f(vj)dvf(\bm{v_{j}})\mathrm{d}\bm{v}f(vj)dv,然后我们可以得出相对速率的表达式:
vri,j=f(vi)dvf(vj)dv∣vi−vj∣\bm{v_{r_{i,j}}}=f(\bm{v_{i}})\mathrm{d}\bm{v}f(\bm{v_{j}})\mathrm{d}\bm{v}|\bm{v_{i}}-\bm{v_{j}}| vri,j=f(vi)dvf(vj)dv∣vi−vj∣
进而可以写出用概率分布表示平均相对速率的表达式:
vr‾=∫−∞+∞∣v1−v2∣f(v1)f(v2)dv1dv2\overline{v_{r}}=\int_{-\infty}^{+\infty}|\bm{v_{1}}-\bm{v_{2}}|f(\bm{v_{1}})f(\bm{v_{2}})\mathrm{d}\bm{v_{1}}\mathrm{d}\bm{v_{2}} vr=∫−∞+∞∣v1−v2∣f(v1)f(v2)dv1dv2
下面具体说明和计算上述内容:
为了计算简单起见,假定每个分子都是质量为mmm,直径为ddd的小球,再假定除一个分子外其他分子都静止不动,只有那个一分子以平均相对速率vrv_{r}vr运动。
设处在平衡态的气体为理想气体。按麦克斯韦速度分布律,气体中的任一分子1处在速度区间v1∼v1+dv1\bm{v_{1}}\sim \bm{v_{1}}+\mathrm{d}\bm{v_{1}}v1∼v1+dv1内的概率dω1\mathrm{d}\omega_{1}dω1为:
dω1=dN1N=f(v1)dv1=(m2πkT)32e−mv122kTdv1\mathrm{d}\omega_{1}=\frac{\mathrm{d}N_{1}}{N}=f(\bm{v_{1}})\mathrm{d}\bm{v_{1}}=\left(\frac{m}{2\pi kT}\right)^{\frac{3}{2}}e^{-\frac{m\bm{v_{1}}^{2}}{2kT}}\mathrm{d}\bm{v_{1}} dω1=NdN1=f(v1)dv1=(2πkTm)23e−2kTmv12dv1
气体中的任一分子2处在速度区间v2∼v2+dv2\bm{v_{2}}\sim \bm{v_{2}}+\mathrm{d}\bm{v_{2}}v2∼v2+dv2内的概率dω2\mathrm{d}\omega_{2}dω2为:
dω2=dN2N=f(v2)dv2=(m2πkT)32e−mv222kTdv2\mathrm{d}\omega_{2}=\frac{\mathrm{d}N_{2}}{N}=f(\bm{v_{2}})\mathrm{d}\bm{v_{2}}=\left(\frac{m}{2\pi kT}\right)^{\frac{3}{2}}e^{-\frac{m\bm{v_{2}}^{2}}{2kT}}\mathrm{d}\bm{v_{2}} dω2=NdN2=f(v2)dv2=(2πkTm)23e−2kTmv22dv2
由于分子1处在速度区间v1∼v1+dv1\bm{v_{1}}\sim \bm{v_{1}}+\mathrm{d}\bm{v_{1}}v1∼v1+dv1,分子2处在速度区间v2∼v2+dv2\bm{v_{2}}\sim \bm{v_{2}}+\mathrm{d}\bm{v_{2}}v2∼v2+dv2是两个独立事件,所以分子1和分子2同事存在于上述两个区间的概率dω\mathrm{d}\omegadω为:
dω=dω1⋅dω2=f(v1)dv1⋅f(v2)dv2=(m2πkT)32e−mv122kTdv1⋅(m2πkT)32e−mv222kTdv2\begin{aligned} \mathrm{d}\omega &=\mathrm{d}\omega_{1}\cdot \mathrm{d}\omega_{2}\\ &=f(\bm{v_{1}})\mathrm{d}\bm{v_{1}}\cdot f(\bm{v_{2}})\mathrm{d}\bm{v_{2}}\\ &=\left(\frac{m}{2\pi kT}\right)^{\frac{3}{2}}e^{-\frac{mv_{1}^{2}}{2kT}}\mathrm{d}\bm{v_{1}}\cdot \left(\frac{m}{2\pi kT}\right)^{\frac{3}{2}}e^{-\frac{mv_{2}^{2}}{2kT}}\mathrm{d}\bm{v_{2}} \end{aligned} dω=dω1⋅dω2=f(v1)dv1⋅f(v2)dv2=(2πkTm)23e−2kTmv12dv1⋅(2πkTm)23e−2kTmv22dv2
假设分子1相对于分子2的速度为vr\bm{v_{r}}vr,则有:
vr=v1−v2\bm{v_{r}}=\bm{v_{1}}-\bm{v_{2}} vr=v1−v2
设两个分子的质心速度为uc\bm{u_{c}}uc,则有:
uc=mv1+mv2m+m=v1+v22\bm{u_{c}}=\frac{m\bm{v_{1}}+m\bm{v_{2}}}{m+m}=\frac{\bm{v_{1}}+\bm{v_{2}}}{2} uc=m+mmv1+mv2=2v1+v2
则v1\bm{v_{1}}v1和v2\bm{v_{2}}v2可以分别表示为:
v1=uc+vr2,v2=uc−vr2\bm{v_{1}}=\bm{u_{c}}+\frac{\bm{v_{r}}}{2},\bm{v_{2}}=\bm{u_{c}}-\frac{\bm{v_{r}}}{2} v1=uc+2vr,v2=uc−2vr
则有:
mv12+mv22=m(uc+vr2)2+m(uc−vr2)2=m(2uc2+vr22)=2muc2+mvr22\begin{aligned} mv_{1}^{2}+mv_{2}^{2} &=m\left(\bm{u_{c}+\frac{\bm{v_{r}}}{2}}\right)^{2}+m\left(\bm{u_{c}-\frac{\bm{v_{r}}}{2}}\right)^{2}\\ &=m\left(2u_{c}^{2}+\frac{v_{r}^2}{2}\right)\\ &=2mu_{c}^{2}+\frac{mv_{r}^2}{2}\\ \end{aligned} mv12+mv22=m(uc+2vr)2+m(uc−2vr)2=m(2uc2+2vr2)=2muc2+2mvr2
根据积分变量替换法则:
dv1⋅dv2=J⋅duc⋅dvr\mathrm{d}\bm{v_{1}}\cdot \mathrm{d}\bm{v_{2}}=J\cdot \mathrm{d}\bm{u_{c}}\cdot\mathrm{d}\bm{v_{r}} dv1⋅dv2=J⋅duc⋅dvr
其中JJJ为雅可比行列式:
J=∣∂v1∂vr∂v1∂uc∂v2∂vr∂v2∂uc∣=∣1/21−1/21∣=1J=\begin{vmatrix} \frac{\partial \bm{v_{1}}}{\partial \bm{v_{r}}} & \frac{\partial \bm{v_{1}}}{\partial \bm{u_{c}}} \\ \frac{\partial \bm{v_{2}}}{\partial \bm{v_{r}}} & \frac{\partial \bm{v_{2}}}{\partial \bm{u_{c}}} \end{vmatrix} =\begin{vmatrix} 1/2 & 1 \\ -1/ 2& 1 \end{vmatrix}=1 J=∂vr∂v1∂vr∂v2∂uc∂v1∂uc∂v2=1/2−1/211=1
将上述三个公式代入到分子1和分子2的联合分布中得:
dω=f(v1)dv1⋅f(v2)dv2=(m2πkT)32e−mv122kTdv1⋅(m2πkT)32e−mv222kTdv2=(2m2πkT)32e−2muc22kTduc⋅(m/22πkT)32e−m/2⋅vr22kTdvr=f(uc)duc⋅f(vr)dvr\begin{aligned} \mathrm{d}\omega &=f(\bm{v_{1}})\mathrm{d}\bm{v_{1}}\cdot f(\bm{v_{2}})\mathrm{d}\bm{v_{2}}\\ &=\left(\frac{m}{2\pi kT}\right)^{\frac{3}{2}}e^{-\frac{mv_{1}^{2}}{2kT}}\mathrm{d}\bm{v_{1}}\cdot \left(\frac{m}{2\pi kT}\right)^{\frac{3}{2}}e^{-\frac{mv_{2}^{2}}{2kT}}\mathrm{d}\bm{v_{2}}\\ &=\left(\frac{2m}{2\pi kT}\right)^{\frac{3}{2}}e^{-\frac{2mu_{c}^{2}}{2kT}}\mathrm{d}\bm{u_{c}}\cdot \left(\frac{m/2}{2\pi kT}\right)^{\frac{3}{2}}e^{-\frac{m/2\cdot v_{r}^{2}}{2kT}}\mathrm{d}\bm{v_{r}}\\ &=f(\bm{u_{c}})\mathrm{d}\bm{u_{c}}\cdot f(\bm{v_{r}})\mathrm{d}\bm{v_{r}} \end{aligned} dω=f(v1)dv1⋅f(v2)dv2=(2πkTm)23e−2kTmv12dv1⋅(2πkTm)23e−2kTmv22dv2=(2πkT2m)23e−2kT2muc2duc⋅(2πkTm/2)23e−2kTm/2⋅vr2dvr=f(uc)duc⋅f(vr)dvr
其中:
f(uc)duc=(2m2πkT)32e−2muc22kTducf(vr)dvr=(m/22πkT)32e−m/2⋅vr22kTdvr\begin{aligned} f(\bm{u_{c}})\mathrm{d}\bm{u_{c}} &=\left(\frac{2m}{2\pi kT}\right)^{\frac{3}{2}}e^{-\frac{2mu_{c}^{2}}{2kT}}\mathrm{d}\bm{u_{c}}\\ f(\bm{v_{r}})\mathrm{d}\bm{v_{r}} &=\left(\frac{m/2}{2\pi kT}\right)^{\frac{3}{2}}e^{-\frac{m/2\cdot v_{r}^{2}}{2kT}}\mathrm{d}\bm{v_{r}} \end{aligned} f(uc)ducf(vr)dvr=(2πkT2m)23e−2kT2muc2duc=(2πkTm/2)23e−2kTm/2⋅vr2dvr
可以证明上面两个分布都满足归一化。后者即为相对速度分布函数f(vr)f(\bm{v_{r}})f(vr),由其可以得到相应的相对速率分布函数:
f(vr)=4π(m/22πkT)32e−m/2⋅vr22kTvr2dvrf(v_{r})=4\pi\left(\frac{m/2}{2\pi kT}\right)^{\frac{3}{2}}e^{-\frac{m/2\cdot v_{r}^{2}}{2kT}}v_{r}^{2}\mathrm{d}v_{r} f(vr)=4π(2πkTm/2)23e−2kTm/2⋅vr2vr2dvr
所以分子间平均相对速率为:
vr‾=∫−∞+∞f(vr)dvr=∫−∞+∞vr4π(m/22πkT)32e−m/2⋅vr22kTvr2dvr=2⋅8kTπm=2vˉ\begin{aligned} \overline{v_{r}} &=\int_{-\infty}^{+\infty}f(v_{r})\mathrm{d}v_{r}\\ &=\int_{-\infty}^{+\infty}v_{r}4\pi\left(\frac{m/2}{2\pi kT}\right)^{\frac{3}{2}}e^{-\frac{m/2\cdot v_{r}^{2}}{2kT}}v_{r}^{2}\mathrm{d}v_{r}\\ &=\sqrt{2}\cdot\sqrt{\frac{8kT}{\pi m}}\\ &=\sqrt{2}\bar{v} \end{aligned} vr=∫−∞+∞f(vr)dvr=∫−∞+∞vr4π(2πkTm/2)23e−2kTm/2⋅vr2vr2dvr=2⋅πm8kT=2vˉ
由此vr‾=2vˉ\overline{v_{r}}=\sqrt{2}\bar{v}vr=2vˉ,证明完毕。
值得注意的是:
1)在上述推导中,假设了两个分子之间速度分布是独立的,这在分子相距很远的时候可以成立,但是相距很近的时候不一定成立,相距较近的时候可能存在关联。但是由于最终取得是整体分子的平均值,所以最终这种关联被平均掉了,如何理解这一点,是需要后续注意的问题。之所以提到这一点是因为我们要使用这个平均相对速率来求碰撞频率,而碰撞频率就是两个分子相距较近的情况,所以这里的计算是否存在误差,以及多少,都是需要再进一步思考的;
2)除了相对速度的分布之外,我们还得到了质心速度的分布函数f(uc)ducf(\bm{u_{c}})\mathrm{d}\bm{u_{c}}f(uc)duc,考虑到质心为两个分子的质心,如果我们把整个体系的分子两两组合成质量为m′=2mm'=2mm′=2m的分子体系,则这个分布即为新分子的速度分布函数。当然得到这个质心速度分布的过程是取尽了所有的两两组合,而此处只取其中的一种就可以了
单位合在一起写的时候注意保留点“⋅\cdot⋅”,否则很容易出现错误,比如mKmKmK与m⋅Km\cdot Km⋅K是代表不同含义的。
注意气体和液体,与固体界面的摩擦或者黏滞是不同的,固体界面的黏滞不符合牛顿黏滞定律。
热导率
传热公式为:
ˉdQ=−K(dTdz)z0dSdt\newcommand \dbar {{\; \bar{} \hspace{-0.3em} \mathrm d}}\dbar Q=-\mathcal{K}\left(\frac{\mathrm{d}T}{\mathrm{d}z}\right)_{z_{0}}\mathrm{d}S\mathrm{d}t ˉdQ=−K(dzdT)z0dSdt
上式可以转化为:
ˉdQdt=−dT1K(dzdS)z0\newcommand \dbar {{\; \bar{} \hspace{-0.3em} \mathrm d}}\frac{\dbar Q}{\mathrm{d}t}=-\frac{\mathrm{d}T}{\frac{1}{\mathcal{K}}\left(\frac{\mathrm{d}z}{\mathrm{d}S}\right)_{z_{0}}} dtˉdQ=−K1(dSdz)z0dT
上式左侧为单位时间内传输的热量,对比电流的定义:
dqdt=UR=dV1σ(dLdS)\frac{\mathrm{d}q}{\mathrm{d}t}=\frac{U}{R}=\frac{\mathrm{d}V}{\frac{1}{\mathcal{\sigma}}\left(\frac{\mathrm{d}L}{\mathrm{d}S}\right)} dtdq=RU=σ1(dSdL)dV
上下两式对比,电导率与热导率的位置相同,这可以作为热导率这个名字的原因,而非叫做传热系数。
另外导热好的材料,通常导电也好。
黏度系数、热导率、扩散系数与压强的关系,在低温下会呈相关性
黏度系数、热导率、扩散系数与温度的关系
黏度系数η\etaη
1)若采用光滑弹性钢球模型的初级理论,使用分子受一次碰撞就完全“同化”的假设,得到η=(1/3)ρvpλˉ\eta=(1/3)\rho v_{p}\bar{\lambda}η=(1/3)ρvpλˉ;
2)金斯对速度住留做过粗略的估计,得到
η=0.461ρvpλˉ\eta=0.461\rho v_{p}\bar{\lambda}η=0.461ρvpλˉ;
3)查普曼和恩斯考克曾用不同的高深理论得到相同结果,它是个级数,由前几项可得:η≈0.499ρvpλˉ\eta\approx0.499\rho v_{p}\bar{\lambda}η≈0.499ρvpλˉ;
在充分考虑速度住留后,得到的黏度的数字系数就不只是一次项1/3,还应增加作用更微弱的二次项(1/3)2(1/3)^{2}(1/3)2,依次类推,可以得到唯象的级数表达式:
η=[(1/3)+(1/3)2+......]ρvpλˉ=(1/2)ρvpλˉ\begin{aligned} \eta &=[(1/3)+(1/3)^{2}+......]\rho v_{p}\bar{\lambda}\\ &=(1/2)\rho v_{p}\bar{\lambda} \end{aligned} η=[(1/3)+(1/3)2+......]ρvpλˉ=(1/2)ρvpλˉ
- 穿过指定截面的分子的自由程分布函数与分子束中的分子的速率分布函数
首先与自由程的分布与其有关f(λ)=1λˉe−1λf(\lambda)=\frac{1}{\bar{\lambda}}e^{-\frac{1}{\lambda}}f(λ)=λˉ1e−λ1,另外还和自由程本身有关,自由程长的就容易穿过去,则穿过指定截面的分子的自由程分布函数为:
f截面(λ)=λf(λ)λˉ=λλˉ2e−1λf_{截面}(\lambda)=\frac{\lambda f(\lambda)}{\bar{\lambda}}=\frac{\lambda}{\bar{\lambda}^{2}}e^{-\frac{1}{\lambda}} f截面(λ)=λˉλf(λ)=λˉ2λe−λ1
通过上式可以计算通过指定截面的平均自由程:
λˉ截面=∫0+∞λ2f(λ)λˉ=λ2‾λˉ=2λˉ\bar{\lambda}_{截面}=\int_{0}^{+\infty}\frac{\lambda^{2} f(\lambda)}{\bar{\lambda}}=\frac{\overline{\lambda^{2}}}{\bar{\lambda}}=2\bar{\lambda} λˉ截面=∫0+∞λˉλ2f(λ)=λˉλ2=2λˉ
注意黏性,传热,扩散公式都是在非平衡态的物理情景下,但是在推导的一些过程中还是使用了平衡态的一些结论,这是一种近似。
注意焦耳定律和焦耳效应的区别
前者为:
U=U(T)U=U(T) U=U(T)
即内能只是温度的函数。后者为焦耳实验中不符合理想气体的现象,即:
U≠U(T)U\neq U(T) U=U(T)
即内能不仅是温度的函数,还有可能是体积的函数。
值得注意的是理想气体下,焓 HHH 也只是温度的函数,这也是为什么焦耳-汤姆孙实验可以验证实际气体是否为理想气体的依据。
- C′C'C′ 为总热容,CCC 为 1mol1\mathrm{mol}1mol 的热容,ccc 为单位质量的热容。三者之间的关系为:
C′=νC=Mc,C=μc=ˉdqdT\newcommand \dbar {{\; \bar{} \hspace{-0.3em} \mathrm d}}C'=\nu C=Mc,C=\mu c=\frac{\dbar q}{\mathrm{d}T} C′=νC=Mc,C=μc=dTˉdq
其中 ˉdq\newcommand \dbar {{\; \bar{} \hspace{-0.3em} \mathrm d}}\dbar qˉdq 为 1mol1\mathrm{mol}1mol 物质吸收的热量。注意热容并不是一个不变量,在绝热过程中,其值为 0。在等温过程中,分母为 0,其值变为无限大。在压缩过程中,气体放热,温度升高,则值为负数。可见根据上述热容的定义,热容值与具体过程有关。
定压热容和定体热容在固体和液体中差不多,一般不做区分,而在气体中其值不做近似。
第五章 热力学第一定律(一)
- p−V−Tp-V-Tp−V−T 函数图
理想气体其他的二维图都可以经上图投影得到。
通常把严格遵守理想气体物态方程的气体称为理想气体。因而也可以说,凡严格服从波意尔定律和阿伏伽德罗定律的气体是理想气体。实际上,理想气体还严格服从焦耳定律——气体的内能只是温度的函数。因此有人认为:只有严格服从波意耳定律、阿伏伽德罗定律和焦耳定律的气体才是理想气体。
但是以后在理论中又可以证明:凡严格遵守理想气体物态方程的气体必定也严格服从焦耳定律。因而,最终可以说:凡严格遵守理想气体物态方程的气体就是理想气体。当然也可以说:凡严格服从波意耳定律和阿伏伽德罗定律的气体就是理想气体。理想气体物态方程的微分形式
理想气体的物态方程为:
pV=νRT=(M/μ)RTpV=\nu RT=(M/\mu)RT pV=νRT=(M/μ)RT
将上式左右取对数:
lnp+lnV=ln(νR)+lnT=ln(MR/μ)+lnT\begin{aligned} \ln p+\ln V &=\ln(\nu R)+\ln T\\ &=\ln(MR/\mu)+\ln T \end{aligned} lnp+lnV=ln(νR)+lnT=ln(MR/μ)+lnT
将上式取全微分:
(dp/p)+(dV/V)=dT/T(\mathrm{d}p/p)+(\mathrm{d}V/V)=\mathrm{d}T/T (dp/p)+(dV/V)=dT/T
上式便是物态方程的微分形式。
在传热的定义中,定义的不是传递热量,因为这是热质说的比喻,而应该从功的角度定义,即一个系统在不受机械功和电力作功的情况下,内能增加的原因来自于热量的输入。
ΔU=A+Q\Delta U=A+QΔU=A+Q
注意上式外界作功为正,从外界吸热为正。但是在工程技术领域,由于研究的对象往往是机械,需要对外作功,所以作功的方向正好相反,以对外作功为正。
- dU=ˉdA+ˉdQ\mathrm{d}U=\newcommand \dbar {{\; \bar{} \hspace{-0.3em} \mathrm d}}\dbar A+\dbar QdU=ˉdA+ˉdQ
注意上式两个不完全微分得到一个完全微分的态函数,后者与过程无关。
注意上图并不是一一对应的,同一个 ppp 可以对应多个 VVV 的值,但这是可以允许的,只需将理想气体状态方程中的温度 TTT 进行变化,脱离等温过程即可。
另外功是一个过程量:ˉdA=−poutdV=−pindV\newcommand \dbar {{\; \bar{} \hspace{-0.3em} \mathrm d}}\dbar A=-p_{out}\mathrm{d}V=-p_{in}\mathrm{d}VˉdA=−poutdV=−pindV,即在准静态作功的情况下外部压强poutp_{out}pout作功可以转换为内部压强pinp_{in}pin。另外,功只在过程中体现,不能代表系统的某一瞬间的状态。而功 AAA 这个量之所以对应不完整微分,也就是因为其不能代表系统的某一时刻的状态,过程结束之后,其随之结束,而之所以不能代表系统某一时刻的状态,就是因为其与路径相关,对应相同的系统,其可能对应不同的值,所以其不具有作为系统态函数的基本属性。
另外注意之所以没有把作功写成ˉdA=ˉdAp+ˉdAV=−pindV+Vdp\newcommand \dbar {{\; \bar{} \hspace{-0.3em} \mathrm d}}\dbar A=\dbar A_{p}+\dbar A_{V}=-p_{in}\mathrm{d}V+V\mathrm{d}pˉdA=ˉdAp+ˉdAV=−pindV+Vdp,是因为体积不变的情况下,力的作用位移为0,所以外界无法对其作功(或者说外界作功的过程都可以表示为 −pdV-p\mathrm{d}V−pdV,这个在理想气体物态图中可以看出,即黎曼积分与勒贝格积分差别)。
上述的内容前提都是在准静态下的过程,当系统处于非静态下的话,在等压过程中,则系统在恒定的外界压强下体积由VA→VBV_{A}\rightarrow V_{B}VA→VB的时候,外界所作的功为:
W=∫VAVB−poutdV=−pout(VB−VA)=−pin(VB−VA)=−pinΔVW=\int_{V_{A}}^{V_{B}}-p_{out}\mathrm{d}V=-p_{out}(V_{B}-V_{A})=-p_{in}(V_{B}-V_{A})=-p_{in}\Delta V W=∫VAVB−poutdV=−pout(VB−VA)=−pin(VB−VA)=−pinΔV
在上述等压过程中虽然系统内部可能发生剧烈的变化,其压强可能并不维持固定,甚至内部各部分的压强也可能并不相等,但是外界的压强一直保持不变,所以可以先计算外界压强的作功情况,然后再把最终结果对应到内部压强中去。得到外部压强的最终表达式之后,可见最终可以只取系统的初态和终态对应的内部压强,而其都是平衡态,初态和终态的压强必定等于外界的压强,所以上式中的ppp仍然为描述系统平衡态的参量。
实验中一般用到定压热容,主要是因为定压的条件比较容易实现。
焦耳定律
由于在焦耳实验中,气体膨胀到真空中,并无对外作功,所以气体系统内能不变,利用内能UUU的全微分公式可得:
dU=(∂U∂T)VdT+(∂U∂V)TdV=0\mathrm{d}U=\left(\frac{\partial U}{\partial T}\right)_{V}\mathrm{d}T+\left(\frac{\partial U}{\partial V}\right)_{T}\mathrm{d}V=0 dU=(∂T∂U)VdT+(∂V∂U)TdV=0
由上式可以得到:
(∂U∂T)V(∂T∂V)U=−(∂U∂V)T\left(\frac{\partial U}{\partial T}\right)_{V}\left(\frac{\partial T}{\partial V}\right)_{U}=-\left(\frac{\partial U}{\partial V}\right)_{T} (∂T∂U)V(∂V∂T)U=−(∂V∂U)T
注意到在一定条件下,U=U(T,V)U=U(T,V)U=U(T,V)可以写成隐函数的形式F(U,T,V)=0F(U,T,V)=0F(U,T,V)=0,之所以写成U=U(T,V)U=U(T,V)U=U(T,V)这个形式,是因为外界不做功,其与ppp无关。
当T,VT,VT,V的数值给定后,UUU的数值必须满足F(U,T,V)=0F(U,T,V)=0F(U,T,V)=0,因而是T,VT,VT,V的函数。不过在这个式子中,T,V,UT,V,UT,V,U三个变量的地位是平等的,因此也可将TTT看做V,UV,UV,U的函数,或者将VVV看做T,UT,UT,U的函数。
由于F(U,T,V)=0F(U,T,V)=0F(U,T,V)=0,所以T,V,UT,V,UT,V,U三个变量的增量dT\mathrm{d}TdT、dV\mathrm{d}VdV、dU\mathrm{d}UdU不是任意的,必须满足一定条件,也就是可以写出F(T,V,U)F(T,V,U)F(T,V,U)的全微分方程:
dF=∂F∂TdT+∂F∂VdV+∂F∂UdU=0\mathrm{d}F=\frac{\partial F}{\partial T}\mathrm{d}T+\frac{\partial F}{\partial V}\mathrm{d}V+\frac{\partial F}{\partial U}\mathrm{d}U=0 dF=∂T∂FdT+∂V∂FdV+∂U∂FdU=0
如果令VVV保持不变,即在上式中令dV=0\mathrm{d}V=0dV=0,得:
(∂U∂T)V=−(∂F∂T)V,U(∂F∂U)V,T,(∂T∂U)V=−(∂F∂U)T,V(∂F∂T)V,U\left(\frac{\partial U}{\partial T}\right)_{V}=-\frac{\left(\frac{\partial F}{\partial T}\right)_{V,U}}{\left(\frac{\partial F}{\partial U}\right)_{V,T}},\left(\frac{\partial T}{\partial U}\right)_{V}=-\frac{\left(\frac{\partial F}{\partial U}\right)_{T,V}}{\left(\frac{\partial F}{\partial T}\right)_{V,U}} (∂T∂U)V=−(∂U∂F)V,T(∂T∂F)V,U,(∂U∂T)V=−(∂T∂F)V,U(∂U∂F)T,V
上述两式比较可得:
(∂U∂T)V=1/(∂T∂U)V\left(\frac{\partial U}{\partial T}\right)_{V}=1/\left(\frac{\partial T}{\partial U}\right)_{V} (∂T∂U)V=1/(∂U∂T)V
同理可以得到:
(∂T∂V)U=1/(∂V∂T)U\left(\frac{\partial T}{\partial V}\right)_{U}=1/\left(\frac{\partial V}{\partial T}\right)_{U} (∂V∂T)U=1/(∂T∂V)U
(∂U∂V)T=1/(∂V∂U)T\left(\frac{\partial U}{\partial V}\right)_{T}=1/\left(\frac{\partial V}{\partial U}\right)_{T} (∂V∂U)T=1/(∂U∂V)T
所以(∂U∂T)V(∂T∂V)U=−(∂U∂V)T\left(\frac{\partial U}{\partial T}\right)_{V}\left(\frac{\partial T}{\partial V}\right)_{U}=-\left(\frac{\partial U}{\partial V}\right)_{T}(∂T∂U)V(∂V∂T)U=−(∂V∂U)T可化为:
(∂U∂T)V(∂T∂V)U(∂V∂U)T=−1\left(\frac{\partial U}{\partial T}\right)_{V}\left(\frac{\partial T}{\partial V}\right)_{U}\left(\frac{\partial V}{\partial U}\right)_{T}=-1 (∂T∂U)V(∂V∂T)U(∂U∂V)T=−1
上述隐函数部分来自汪志诚《热力学·统计物理》附录A中的推导,值得注意的是,dF=0\mathrm{d}F=0dF=0,但是∂F∂T\frac{\partial F}{\partial T}∂T∂F、∂F∂V\frac{\partial F}{\partial V}∂V∂F、∂F∂U\frac{\partial F}{\partial U}∂U∂F不一定为0(待确认!!!
)。
上述证明用到了U=const.U=\mathrm{const}.U=const.这个条件,所以并不普适,更普遍的证明过程可以参考汪志诚《热力学·统计物理》附录A中的内容。
- 考虑范德瓦尔斯模型的焦耳系数
焦耳系数为:
λ=(∂T∂V)U=−(∂U∂V)T(∂U∂T)V=−(∂U∂V)TCV\lambda=\left(\frac{\partial T}{\partial V}\right)_{U}=-\frac{\left(\frac{\partial U}{\partial V}\right)_{T}}{\left(\frac{\partial U}{\partial T}\right)_{V}}=-\frac{\left(\frac{\partial U}{\partial V}\right)_{T}}{C_{V}} λ=(∂V∂T)U=−(∂T∂U)V(∂V∂U)T=−CV(∂V∂U)T
考虑范德瓦尔斯模型的情况下,内能可以表示为:
U=U0+CVT+Ep(V)U=U_{0}+C_{V}T+E_{p}(V) U=U0+CVT+Ep(V)
其中最后一项Ep(V)E_{p}(V)Ep(V)为随体积变化的势能项,可以由内压强pinp_{in}pin来计算:
Ep(V)=∫−∞V−pindV=∫∞V−(−aV2)dV=−aV∣∞V=−aVE_{p}(V)=\int_{-\infty}^{V}-p_{in}\mathrm{d}V=\int_{\infty}^{V}-\left(-\frac{a}{V^{2}}\right)\mathrm{d}V=-\frac{a}{V}|_{\infty}^{V}=-\frac{a}{V} Ep(V)=∫−∞V−pindV=∫∞V−(−V2a)dV=−Va∣∞V=−Va
注意积分限在∞\infty∞即无穷远的时候,势能为0,在体积缩小的过程中势能不断的减小(在大部分范围内为引力),所以最终的得数为负值。
所以在范德瓦尔斯方程下的内能可以重新表达为:
U=U0+CVT−aVU=U_{0}+C_{V}T-\frac{a}{V} U=U0+CVT−Va
将上式带入到焦耳系数表达式,可以得到:
λ=−(∂U∂V)TCV=−aV2CV\lambda=-\frac{\left(\frac{\partial U}{\partial V}\right)_{T}}{C_{V}}=-\frac{a}{V^{2}C_{V}} λ=−CV(∂V∂U)T=−V2CVa
根据上式可以看到焦耳定律在理想气体下虽然为0,但是在范德瓦尔斯方程描述下的气体系数为负值,即体积增大的时候,内能下降。所以焦耳定律仅仅是一个理想气体的标准,即忽略势能的理想气体系统适用,而非实际情况。
注意:作功与传热都是过程量;作功与传热的等当性(焦耳的热功当量实验确定)
焦耳-汤姆孙实验
焦耳实验为气体的绝热自由膨胀,而焦耳-汤姆孙实验为气体的绝热节流膨胀,两者有两字之差。
由于整个系统分为节流前和节流后两部分,等整个系统达到稳态的时候,两部分的能量公式可以写为(考虑绝热过程,所以 ˉdQ=0\newcommand \dbar {{\; \bar{} \hspace{-0.3em} \mathrm d}}\dbar Q=0ˉdQ=0):
dU=ˉdQ−pdV=−pdV⇓U1−U2=−pV1+pV2⇓U1+pV1=U2+pV2⇓H2=H1\newcommand \dbar {{\; \bar{} \hspace{-0.3em} \mathrm d}} \mathrm{d}U =\dbar Q-p\mathrm{d}V=-p\mathrm{d}V\\ \Downarrow\\ U_{1}-U_{2}=-pV_{1}+pV_{2}\\ \Downarrow\\ U_{1}+pV_{1}=U_{2}+pV_{2}\\ \Downarrow\\ H_{2}=H_{1} dU=ˉdQ−pdV=−pdV⇓U1−U2=−pV1+pV2⇓U1+pV1=U2+pV2⇓H2=H1
所以焦耳-汤姆孙实验是一个等焓过程。
在理想气体下焓只为温度的函数:
H(T)=U+pV=U(T)+νRTH(T)=U+pV=U(T)+\nu RT H(T)=U+pV=U(T)+νRT
α=(∂T∂p)h=−(∂h∂p)TCp\alpha=\left(\frac{\partial T}{\partial p}\right)_{h}=-\frac{\left(\frac{\partial h}{\partial p}\right)_{T}}{C_{p}} α=(∂p∂T)h=−Cp(∂p∂h)T
由于焦耳-汤姆孙为一个绝热膨胀的实验,所以压强只会减小,所以当 α>0\alpha>0α>0 时,ΔT<0\Delta T<0ΔT<0 把这个过程称为致冷过程;反之,当 α<0\alpha<0α<0 时,ΔT>0\Delta T>0ΔT>0 把这个过程称为致温过程。
由于在理想气体情况下,焓只为温度的函数。所以如果真实气体的焦耳-汤姆孙实验的结果为焦耳-汤姆孙系数不为 0,相当于pVpVpV 的变化导致内能的变化,过程中克服内聚力做了功,此时真实气体不满足焦耳定律,也不满足理想气体的物态方程。
- 范德瓦尔斯气体的焦耳-汤姆孙效应
由一摩尔气体的范德瓦尔斯方程可知:
RT=[p+(a/ν2)](ν−b)=pν+(a/ν)−b[p+(a/ν2)]\begin{aligned} RT&=[p+(a/\nu^{2})](\nu-b)\\ &=p\nu+(a/\nu)-b[p+(a/\nu^{2})] \end{aligned} RT=[p+(a/ν2)](ν−b)=pν+(a/ν)−b[p+(a/ν2)]
而其摩尔内能 uuu 为:
u=u0+CνT−(a/ν)u=u_{0}+C_{\nu}T-(a/\nu) u=u0+CνT−(a/ν)
因此其摩尔焓 hhh 为:
h=u+pν=u0+CνT−(a/ν)+RT−(a/ν)+b[p+(a/ν2)]=u0+(Cν+R)T−2a(1/ν)+b[p+(a/ν2)]\begin{aligned} h&=u+p\nu\\ &=u_{0}+C_{\nu}T-(a/\nu)+RT-(a/\nu)+b[p+(a/\nu^{2})]\\ &=u_{0}+(C_{\nu}+R)T-2a(1/\nu)+b[p+(a/\nu^{2})] \end{aligned} h=u+pν=u0+CνT−(a/ν)+RT−(a/ν)+b[p+(a/ν2)]=u0+(Cν+R)T−2a(1/ν)+b[p+(a/ν2)]
由于焦耳-汤姆孙实验是一个绝热节流膨胀过程,是等焓过程,故:
Δh=(Cν+R)ΔT−2aΔ(1/ν)+bΔ[p+(a/ν2)]=0\Delta h=(C_{\nu}+R)\Delta T-2a\Delta (1/\nu)+b\Delta[p+(a/\nu^{2})]=0 Δh=(Cν+R)ΔT−2aΔ(1/ν)+bΔ[p+(a/ν2)]=0
由此可求得范德瓦尔斯气体在焦耳-汤姆孙实验中温度的改变 ΔT\Delta TΔT 为:
ΔT=(Cν+R)−1{2aΔ(1/ν)−bΔ[p+(a/ν2)]}\Delta T=(C_{\nu}+R)^{-1}\{2a\Delta(1/\nu)-b\Delta[p+(a/\nu^{2})]\} ΔT=(Cν+R)−1{2aΔ(1/ν)−bΔ[p+(a/ν2)]}
因为 aaa,bbb,CνC_{\nu}Cν,RRR 都是正值常量,所以 Cν+RC_{\nu}+RCν+R 也是正值常量。因此,范德瓦尔斯气体经过绝热节流膨胀过程后其温度的增减也就完全由 2aΔ(1/ν)−bΔ[p+(a/ν2)]2a\Delta (1/\nu)-b\Delta[p+(a/\nu^{2})]2aΔ(1/ν)−bΔ[p+(a/ν2)] 决定。
由于 Δν>0\Delta \nu>0Δν>0,所以 Δ(1/ν2)<0\Delta(1/\nu^{2})<0Δ(1/ν2)<0,并且 Δp<0\Delta p<0Δp<0,因而 −bΔ[p+(a/ν2)]>0-b\Delta[p+(a/\nu^{2})]>0−bΔ[p+(a/ν2)]>0,这就造成 Δ>0\Delta>0Δ>0,出现致温效应。由此可见,范德瓦尔斯改正量 bbb 将导致产生负效应。
由于 Δν>0\Delta \nu>0Δν>0,所以 Δ(1/ν)<0\Delta(1/\nu)<0Δ(1/ν)<0,因而 2aΔ(1/ν)<02a\Delta(1/\nu)<02aΔ(1/ν)<0,这就造成 ΔT<0\Delta T<0ΔT<0,出现致冷效应。由此可见,范德瓦尔斯改正量 aaa 主要将导致产生正效应。
综上,范德瓦尔斯气体对理想气体物态方程的偏离可用 2aΔ(1/ν)−bΔ[p+(a/ν2)]2a\Delta (1/\nu)-b\Delta[p+(a/\nu^{2})]2aΔ(1/ν)−bΔ[p+(a/ν2)] 代表,它既能导致产生正效应,又能导致产生负效应,还能导致产生零效应。
范德瓦尔斯气体对焦耳定律的偏离,可用 aΔ(1/ν)a\Delta(1/\nu)aΔ(1/ν) 代表,它只能导致产生正效应。到底最后出现怎样的效应,还要由 a/ba/ba/b 的具体数值的大小来决定。
从教科书中所列出的 10 种气体的 aaa 和 bbb 的数值来看,诸气体的 bbb 相当接近,而它们的 aaa 却相差较大,因此关键要看 aaa。如果以 atm⋅L2⋅mol2\mathrm{atm}\cdot\mathrm{L^{2}}\cdot\mathrm{mol^{2}}atm⋅L2⋅mol2 作为 aaa 的单位,He\mathrm{He}He 的数量级为 10−210^{-2}10−2,H2\mathrm{H_{2}}H2 和 Ne\mathrm{Ne}Ne 的数量级为 10−110^{-1}10−1,其他气体的数量级为 10010^{0}100。
因此,一般气体在焦耳-汤姆孙实验中往往将出现正效应,温度会降低。而 He\mathrm{He}He 在焦耳-汤姆孙实验中却常常出现负效应,温度会升高。
当我们让真实气体经过一个无限小的绝热节流膨胀过程,前后的温度和压强尽管有一些变化,肯定也是很小的,这时候可以近似看成两个相同的值,然后就可以在 T−pT-pT−p 图中点一个点,表示这个无限小的绝热节流膨胀过程大致是在哪一个态下进行的,进而可以将正效应和负效应表示在 T−pT-pT−p 图中,可以表示为:
其中曲线可以证明为抛物线,称为反转温度曲线,上方为上反转温度曲线,下方为下反转温度曲线。其中曲线与纵轴组成的封闭区域内为正效应,外部为正效应,另外横轴 ppp 只取正值。抛物线的最大值为 T0=2aRbT^{0}=\frac{2a}{Rb}T0=Rb2a,在抛物线上既不会出现正效应也不会出现负效应。
HHH,HeHeHe 常温下都在反转温度曲线的上方,即负效应区域,此时的膨胀过程为致温的,所以无法制冷,所以为了制冷,需要先将气体降温至正效应区域内,再行绝热膨胀,则会致冷/制冷。
- 焦耳实验与焦耳-汤姆孙实验对比
- 焦耳实验是绝热自由膨胀过程,是等内能过程:
ΔU=0,Δν>0\Delta U=0,\Delta \nu>0 ΔU=0,Δν>0
而焦耳-汤姆孙实验是绝热节流膨胀过程,是等焓过程:
ΔH=0,Δp<0\Delta H=0,\Delta p<0 ΔH=0,Δp<0
- 焦耳实验可以用焦耳系数 λ\lambdaλ 描述:
λ=(∂T/∂ν)u=−(∂u/∂ν)T/Cν\begin{aligned} \lambda&=(\partial T/\partial \nu)_{u} &=-(\partial u/\partial \nu)_{T}/C_{\nu} \end{aligned} λ=(∂T/∂ν)u=−(∂u/∂ν)T/Cν
焦耳-汤姆孙实验可用焦耳-汤姆孙系数 α\alphaα 描述:
α=(∂T/∂p)h=−(∂h/∂p)T/Cp\begin{aligned} \alpha&=(\partial T/\partial p)_{h} &=-(\partial h/\partial p)_{T}/C_{p} \end{aligned} α=(∂T/∂p)h=−(∂h/∂p)T/Cp
- 理想气体在焦耳实验中温度不变:
ΔT=0,λ=0\Delta T=0,\lambda=0 ΔT=0,λ=0
这表明理想气体的内能只是温度的函数。
理想气体在焦耳-汤姆孙实验中温度不变:
ΔT=0,α=0\Delta T=0,\alpha=0 ΔT=0,α=0
这表明理想气体的焓只是温度的函数。
- 焦耳实验表明实际气体不服从焦耳定律,但出入不大;有焦耳效应,却不显著。
焦耳-汤姆孙实验表明:实际气体既不服从焦耳定律,又不服从理想气体物态方程;通常不出现零效应。
- 范德瓦尔斯气体在焦耳实验中温度将下降:
ΔT<0,λ<0\Delta T<0,\lambda<0 ΔT<0,λ<0
范德瓦尔斯气体在焦耳-汤姆孙实验中温度既可能降低,又可能升高、还可能保持不变。
- 焦耳-汤姆孙实验相比焦耳实验好在哪里?
- 焦耳-汤姆孙实验是在稳态下进行测量的,即使温度只有 0.1K0.1\mathrm{K}0.1K 的改变,也可以准确地测出来;而焦耳实验则不然。
- 焦耳-汤姆孙实验中温度的改变常常比焦耳实验中温度的改变要大一些,更容易被准确地测量。
- 对于极稀薄的实际气体,焦耳-汤姆孙实验往往总能得到与理想气体不同的结果;而焦耳实验却会得到与理想气体相同的结果,此点特别重要。
如果把真实气体用范德瓦尔斯模型来处理,它们在焦耳实验和焦耳-汤姆孙实验中将会得到明显不同的结果。
范德瓦尔斯气体的摩尔内能为:
u=u0+CνT−(a/ν)u=u_{0}+C_{\nu}T-(a/\nu) u=u0+CνT−(a/ν)
范德瓦尔斯气体的焦耳系数为:
λ=(∂T/∂ν)u=−(∂u/∂ν)T/Cν=−a/(Cνν2)\begin{aligned} \lambda&=(\partial T/\partial \nu)_{u}\\ &=-(\partial u/\partial \nu)_{T}/C_{\nu}\\ &=-a/(C_{\nu}\nu^{2}) \end{aligned} λ=(∂T/∂ν)u=−(∂u/∂ν)T/Cν=−a/(Cνν2)
当 ν→∞\nu\rightarrow\inftyν→∞ 时,λ→0\lambda\rightarrow0λ→0,由此可见,极稀薄的实际气体在焦耳实验中将得到和理想气体相同的结果。
范德瓦尔斯气体的摩尔焓为:
h=u0+(Cν+R)T−(2a/ν)+b[p+(a/ν2)]h=u_{0}+(C_{\nu}+R)T-(2a/\nu)+b[p+(a/\nu^{2})] h=u0+(Cν+R)T−(2a/ν)+b[p+(a/ν2)]
(注意与高崇伊老师的 ppt 中有所区别:h=u0+(Cν+R)T−(2a/ν)+[RTb/(ν−b)]h=u_{0}+(C_{\nu}+R)T-(2a/\nu)+[RTb/(\nu-b)]h=u0+(Cν+R)T−(2a/ν)+[RTb/(ν−b)])
然后可以求得范德瓦尔斯气体的焦耳-汤姆孙系数为:
α=(∂T/∂p)h=−(∂h/∂p)T/Cp=−(b/Cp){(T0/T)[1−(b/ν)]2−1}/{[(T0/T)[1−(b/ν)]2(b/ν)]−1}\begin{aligned} \alpha&=(\partial T/\partial p)_{h}\\ &=-(\partial h/\partial p)_{T}/C_{p}\\ &=-(b/C_{p})\{(T^{0}/T)[1-(b/\nu)]^{2}-1\}/\{[(T^{0}/T)[1-(b/\nu)]^{2}(b/\nu)]-1\} \end{aligned} α=(∂T/∂p)h=−(∂h/∂p)T/Cp=−(b/Cp){(T0/T)[1−(b/ν)]2−1}/{[(T0/T)[1−(b/ν)]2(b/ν)]−1}
其中 T0=2aRbT^{0}=\frac{2a}{Rb}T0=Rb2a。所以当上式中 ν→∞\nu\rightarrow\inftyν→∞ 时,b/ν→0b/\nu\rightarrow0b/ν→0,因而 α→b[(T0/T)−1]/Cp\alpha\rightarrow b[(T^{0}/T)-1]/C_{p}α→b[(T0/T)−1]/Cp。
由此可见,除非 T=T0=2a/(Rb)T=T^{0}=2a/(Rb)T=T0=2a/(Rb),否则极稀薄的实际气体在焦耳-汤姆孙实验中将得到和理想气体不同的,非零的结果。
- 迈耶关系
H=U+pV=U(T)+νRT=H(T)H=U+pV=U(T)+\nu RT=H(T) H=U+pV=U(T)+νRT=H(T)
则在 1mol1\mathrm{mol}1mol 理想气体下:
Cp−Cν=(∂h∂T)p−(∂u∂T)ν=dhdT−dudT=RC_{p}-C_{\nu}=\left(\frac{\partial h}{\partial T}\right)_{p}-\left(\frac{\partial u}{\partial T}\right)_{\nu}=\frac{\mathrm{d}h}{\mathrm{d}T}-\frac{\mathrm{d}u}{dT}=R Cp−Cν=(∂T∂h)p−(∂T∂u)ν=dTdh−dTdu=R
其中 hhh 为 1mol1\mathrm{mol}1mol 气体的焓。注意上述公式是理想气体下的结果。
另外在理想气体中,如果假设 CpCν=γ\frac{C_{p}}{C_{\nu}}=\gammaCνCp=γ,则通过迈耶关系:
Cν=(∂u/∂T)ν=du/dT=R/(γ−1)Cp=(∂h/∂T)p=dh/dT=γR/(γ−1)\begin{aligned} C_{\nu}&=(\partial u/\partial T)_{\nu}=\mathrm{d}u/\mathrm{d}T\\ &=R/(\gamma-1)\\ C_{p}&=(\partial h/\partial T)_{p}=\mathrm{d}h/\mathrm{d}T\\ &=\gamma R/(\gamma-1) \end{aligned} CνCp=(∂u/∂T)ν=du/dT=R/(γ−1)=(∂h/∂T)p=dh/dT=γR/(γ−1)
值得注意的是在理想气体中,CpC_{p}Cp,CνC_{\nu}Cν 一般为温度的函数,所以 γ\gammaγ 也为温度的函数。如果后者不为温度的函数,则前两者也不为温度的函数,反之亦然。
- 内能和焓的改变
ΔU=Qν=ν∫T0TCνTΔH=Qp=ν∫T0TCpT\Delta U=Q_{\nu}=\nu \int_{T_{0}}^{T}C_{\nu}\mathrm{T}\\ \Delta H=Q_{p}=\nu \int_{T_{0}}^{T}C_{p}\mathrm{T} ΔU=Qν=ν∫T0TCνTΔH=Qp=ν∫T0TCpT
其中下标 ν\nuν 为单位体积。
其中若 γ=const.\gamma=\mathrm{const}.γ=const.,即 Cν=const.C_{\nu}=\mathrm{const}.Cν=const.,Cp=const.C_{p}=\mathrm{const}.Cp=const.,则:
ΔU=Qν=νCνΔT=νRΔT/(γ−1)=Δ(νRT)/(γ−1)=(p2V2−p1V1)/(γ−1)\begin{aligned} \Delta U&=Q_{\nu}\\ &=\nu C_{\nu}\Delta T\\ &=\nu R\Delta T/(\gamma-1)\\ &=\Delta (\nu RT)/(\gamma-1)\\ &=(p_{2}V_{2}-p_{1}V_{1})/(\gamma-1)\\ \end{aligned} ΔU=Qν=νCνΔT=νRΔT/(γ−1)=Δ(νRT)/(γ−1)=(p2V2−p1V1)/(γ−1)
ΔH=Qp=νCpΔT=νγRΔT/(γ−1)=γΔ(νRT)/(γ−1)=γ(p2V2−p1V1)/(γ−1)=γΔU\begin{aligned} \Delta H&=Q_{p}\\ &=\nu C_{p}\Delta T\\ &=\nu \gamma R\Delta T/(\gamma-1)\\ &=\gamma\Delta (\nu RT)/(\gamma-1)\\ &=\gamma(p_{2}V_{2}-p_{1}V_{1})/(\gamma-1)\\ &=\gamma\Delta U \end{aligned} ΔH=Qp=νCpΔT=νγRΔT/(γ−1)=γΔ(νRT)/(γ−1)=γ(p2V2−p1V1)/(γ−1)=γΔU
注意上述计算结果只要初末态为平衡态,则对于非等体和非等压的平衡/非平衡过程同样适用,因为中间使用的变换是普适的。
由于内能和焓都是态函数,所以上述计算内能和焓的改变的公式,对于初末态为平衡态的任何过程都是适用的。
- 参考文献:
【热学】兰州大学 高崇伊授课 李椿版
On the Thermodynamics of Solutions. V. An Equation of State. Fugacities of Gaseous Solutions by Otto. Redlich and J. N. S. Kwong
Redlich–Kwong equation of state
王家慧,申兵辉,王卫,韩萍.气体分子平均相对速率与平均速率[J].广西物理,2013,34(03):46-48.
王明美.理想气体分子平均相对速率公式的推导[J].合肥师范学院学报,2009,27(06):27-28.
郑勇.理想气体分子相对速度分布的推导[J].广西物理,2015,36(03):33-36.
朱旭宁.混合理想气体平均相对速率的求解[J].商洛师范专科学校学报,2005(03):94-95+99.
“热力学·统计物理” by 汪志诚
wiki: Gaseous diffusion
待将高崇伊P13节斥力引力有效半径整理一下,以及有效半径与温度的关系。
自学脚手架——《热学》 by 李椿(第一,二,三,四,五章)相关推荐
- 6月29-7月5日成都uber优步司机第一/二/三组奖励政策明细
成都优步司机第一/二/三组奖励更新了,在写下文之前,我先吐槽一下:靠优步uber发财致富已成往事. 滴滴快车单单2.5倍,注册地址:http://www.udache.com/如何注册Uber司机(全 ...
- 第一卷 第一百三十六章 好白菜
第一卷 第一百三十六章 好白菜 其实我最恨两种人,一种是始乱终弃的小白脸,一种就是生孩子不养活的臭娘们儿.** 提供本书txt电子书下载 **既然你能生就要能养啊,虽然我不知道井里的那个死孩子的母亲 ...
- 雷士照明股权连环局(第一 二 三季完整版)
分享一下我老师大神的人工智能教程!零基础,通俗易懂!http://blog.csdn.net/jiangjunshow 也欢迎大家转载本篇文章.分享知识,造福人民,实现我们中华民族伟大复兴! 转自:h ...
- 【数分书单】业务学习《从零开始做运营入门篇》第一/二/三章小结
<从零开始做运营入门篇>by 张亮 第一章 关于运营的一些实话实说 1 苦逼的运营产品两兄弟 运营需要积累,实力的积累,经验的积累,对用户了解程度的积累,对数据敏感程度的积累:长期磨练出 ...
- 用简单的方法做整套UI(教程第一/二/三弹合集)
http://bbs.66rpg.com/thread-329530-1-1.html http://v.tieba.baidu.com/p/2985559487 首先要准备两个工具,"美图 ...
- 第三十五章 二游南斋(恶魔通缉令)
在空中划过一道长长的弧线后,来不及尖叫,身体开始下降. 瞬间的超重让十趾骤然弯曲,白羽转成掣电. 下方是个箭靶么,一圈一圈的,来个十环吧! 呀!半空出现什么,这么滑,像踩在蛋壳上,壳面还附着了大量蛋清 ...
- 李永乐复习全书高等数学 第五章 多元函数微分学
目录 5.1 多元函数的极限.连续.偏导数与全微分(概念) 例2 试求下列二次极限: (1)limy→0x→0x2+y2∣x∣+∣y∣;\lim\limits_{y\to0 \atop x\to ...
- 第一百三十五期:如何模拟一次阿里双11秒杀场景的实现?程序员必看
秒杀活动可以说在互联网上随处可见,从12306抢票,到聚划算抢购,我们生活的方方面面都可以看到秒杀的身影. 作者:IT技术管理那些事儿 秒杀活动可以说在互联网上随处可见,从12306抢票,到聚划算抢购 ...
- 《java编程思想》读书笔记(二)第五章(2)
成员初始化 Java尽力保证:所有变量在使用前都能得到恰当的初始化. 对于方法的局部变量,Java会以编译时报错的形式贯彻这种保证,比如: 1 void f() 2 { 3 int i; //这里编译 ...
- 李永乐复习全书线性代数 第五章 特征值、特征向量、相似矩阵
目录 例题五 例6 设A,B\bm{A},\bm{B}A,B是nnn阶矩阵. (3)A,B\bm{A},\bm{B}A,B均是nnn阶矩阵,证明AB,BA\bm{AB},\bm{BA}AB,BA有相 ...
最新文章
- Pythonic版冒泡排序和快速排序(附:直接插入排序)
- 奇点、技术失控与技术启示录
- python3 赋值 浅拷贝 深拷贝 简介
- pragma指令简介
- 基于stm32f103zet6的DS1302学习
- 7.18 Shell 管道 重定向 链接
- python支持list类型吗_Python-不支持的操作数类型为%:“list”和“int”
- L1-007 念数字 (10 分)—团体程序设计天梯赛
- java.util.ArrayList#add探索
- XDUOJ 1125 Judgement of Orz Pandas
- NOIP2013TGDay1
- 自动阅卷系统/自动阅卷机/网上阅卷
- 微信Native支付申请接入流程-避免踩坑指南
- python解析mht文件_[原创].mht文件图片解析工具
- 那些年你用过最好的键盘
- QQ空间FLASH代码及其使用方法
- 用笔记本做路由器共享4G流量
- STM32内部参考电压的使用
- Mac OSX 升级python six模块版本
- 红光光浴只有流汗的效果?#红光光浴#种光光学