大气热力学(5)——绝热过程

本篇文章源自我在 2021 年暑假自学大气物理相关知识时手写的笔记,现转化为电子版本以作存档。相较于手写笔记,电子版的部分内容有补充和修改。笔记内容大部分为公式的推导过程。

5.1 气块的概念

为了研究大气运动的热力性质,需要研究一个无限小尺度的气块活动,对该气块作如下假设:

  • 气块与周围大气没有热量交换,所以当它作上升或下降运动时,其温度作绝热变化;
  • 气块作升降运动时,其压强与周围大气的压强始终保持相等;
  • 环境空气处于流体静力平衡状态;
  • 气块运动十分缓慢,所以其动能与气块的总能量相比可以忽略。

这个简单、理想化的模型有助于我们理解大气垂直运动的某些物理过程。

5.2 热力学第一定律的几种微分形式

我们知道,热力学第一定律可写为如下微分形式:

d Q = d E + d W = d E + p d V \mathrm{d} Q = \mathrm{d} E + \mathrm{d} W = \mathrm{d} E + p \mathrm{d} V dQ=dE+dW=dE+pdV

现在求 d E \mathrm{d} E dE。我们写出定容比热的定义式(注意定容过程有 d W = 0 \mathrm{d} W = 0 dW=0,于是有 d Q = d E \mathrm{d}Q = \mathrm{d}E dQ=dE):

C V , m = d Q d T = d E d T = i 2 R C_{V, m} = \frac{\mathrm{d}Q}{\mathrm{d}T} = \frac{\mathrm{d}E}{\mathrm{d}T} = \frac{i}{2} R CV,m=dTdQ=dTdE=2iR

整理得:

d E = i 2 R d T = C V , m d T \mathrm{d}E = \frac{i}{2} R \mathrm{d}T = C_{V, m} \mathrm{d}T dE=2iRdT=CV,mdT

代入到上式得:

d Q = C V , m d T + p d V \mathrm{d} Q = C_{V, m} \mathrm{d}T + p \mathrm{d} V dQ=CV,mdT+pdV

这就是热力学第一定律的第二种微分形式

由迈耶公式 C p , m = C V , m + R C_{p,m} = C_{V,m} + R Cp,m=CV,m+R,上式又可以写为:

d Q = C p , m d T + p d V − R d T \mathrm{d} Q = C_{p, m} \mathrm{d}T + p \mathrm{d} V - R \mathrm{d}T dQ=Cp,mdT+pdVRdT

对气体状态方程 p V = R T pV = RT pV=RT 两边微分得:

p d V + V d p = R d T p\mathrm{d}V + V \mathrm{d}p = R \mathrm{d}T pdV+Vdp=RdT

代入上式得:

d Q = C p , m d T − V d p \mathrm{d} Q = C_{p, m} \mathrm{d}T - V \mathrm{d}p dQ=Cp,mdTVdp

这就是热力学第一定律的第三种微分形式,下一节的绝热过程推导需要用到。

对于运动中的单位质量( m = 1 k g m = 1 \mathrm{kg} m=1kg)的气块,注意到流体静力学方程为 d p = − g ρ d z \mathrm{d} p = - g \rho \mathrm{d}z dp=gρdz,如果我们把热力学第一定律与重力位势 d Φ = g d z \mathrm{d} \Phi = g \mathrm{d}z dΦ=gdz 联系起来,那么可以得到如下推导:

d Q = C p , m d T − V d p = C p , m d T + V ⋅ g ρ d z = C p , m d T + ρ V d Φ = C p , m d T + m d Φ = C p , m d T + d Φ \begin{aligned} \mathrm{d} Q &= C_{p, m} \mathrm{d}T - V \mathrm{d}p \\ &= C_{p, m} \mathrm{d}T + V \cdot g \rho \mathrm{d}z \\ &= C_{p, m} \mathrm{d}T + \rho V \mathrm{d} \Phi \\ &= C_{p, m} \mathrm{d}T + m \mathrm{d} \Phi \\ &= C_{p, m} \mathrm{d}T + \mathrm{d} \Phi \end{aligned} dQ=Cp,mdTVdp=Cp,mdT+Vgρdz=Cp,mdT+ρVdΦ=Cp,mdT+mdΦ=Cp,mdT+dΦ

这就是运动中单位气块的热力学表达式。

5.3 干绝热过程

自然界中存在的各种各样的物质,绝大多数都是以固、液、气三种聚集态存在着。为了描述物质的不同聚集态,而用“相”来表示物质的固、液、气三种形态的“相貌”。不同相之间的相互转变,称为相变或称物态变化

大气中进行的物理过程,通常伴有不同形式的能量转换。在能量转换过程中,空气的状态要发生改变。在气象学上,任一气块与外界之间无热量交换时的状态变化过程,叫做绝热过程。在大气中,作垂直运动的气块,其状态变化通常接近于绝热过程。当升、降气块内部既没有发生水相变化,又没有与外界交换热量的过程,称作干绝热过程

对于一团单位质量的气体,其绝热过程的特征为: d Q = 0 \mathrm{d} Q = 0 dQ=0,热力学第一定律可写为:

d Q = C p , m d T − V d p = 0 移项得: C p , m d T = V d p \mathrm{d} Q = C_{p, m} \mathrm{d}T - V \mathrm{d}p = 0 \\ 移项得:C_{p, m} \mathrm{d}T = V \mathrm{d}p dQ=Cp,mdTVdp=0移项得:Cp,mdT=Vdp

再由状态方程 p V = R T pV = RT pV=RT 得到 V = R T p V = \frac{RT}{p} V=pRT,代入上式可得:

C p , m d T = R T p d p 整理得: d T T = R C p , m d p p C_{p,m} \mathrm{d}T = \frac{RT}{p} \mathrm{d}p \\ 整理得:\frac{\mathrm{d}T}{T} = \frac{R}{C_{p,m}} \frac{\mathrm{d}p}{p} Cp,mdT=pRTdp整理得:TdT=Cp,mRpdp

设干绝热过程中,气块的初态为 ( p 0 , T 0 ) (p_0, T_0) (p0,T0),终态为 ( p , T ) (p, T) (p,T),对方程两边进行积分得:

∫ T 0 T d T T = R C p , m ∫ p 0 p d p p ln ⁡ T T 0 = R C p , m ln ⁡ p p 0 T T 0 = ( p p 0 ) R C p , m \begin{aligned} \int_{T_0}^{T} \frac{\mathrm{d}T}{T} &= \frac{R}{C_{p,m}} \int_{p_0}^{p} \frac{\mathrm{d}p}{p} \\ \ln \frac{T}{T_0} &= \frac{R}{C_{p,m}} \ln \frac{p}{p_0} \\ \frac{T}{T_0} &= \bigg( \frac{p}{p_0} \bigg) ^{\frac{R}{C_{p,m}}} \end{aligned} T0TTdTlnT0TT0T=Cp,mRp0ppdp=Cp,mRlnp0p=(p0p)Cp,mR

因为 R C p , m = 0.286 \frac{R}{C_{p,m}} = 0.286 Cp,mR=0.286,所以:

T T 0 = ( p p 0 ) 0.286 \frac{T}{T_0} = \bigg( \frac{p}{p_0} \bigg) ^{0.286} T0T=(p0p)0.286

上式称为干绝热方程,又称为泊松方程

5.4 干绝热递减率(干绝热直减率)

由:

d Q = C p , m d T + d Φ = 0 \mathrm{d} Q = C_{p, m} \mathrm{d}T + \mathrm{d} \Phi = 0 dQ=Cp,mdT+dΦ=0

移项并上下同除以 d z \mathrm{d} z dz 得:

C p , m d T d z = − d Φ d z C_{p, m} \frac{\mathrm{d}T}{\mathrm{d} z} = - \frac{\mathrm{d} \Phi}{\mathrm{d} z} Cp,mdzdT=dzdΦ

整理得到:

− d T d z = g C p , m -\frac{\mathrm{d}T}{\mathrm{d} z} = \frac{g}{C_{p, m}} dzdT=Cp,mg

定义下式为干绝热递减率(用 γ d \gamma_d γd 表示):

γ d = − ( d T d z ) 干空气块 \gamma_d = - \bigg( \frac{\mathrm{d}T}{\mathrm{d} z} \bigg)_{干空气块} γd=(dzdT)干空气块

g = 9.81   m / s 2 g = 9.81 \ \mathrm{m/s^2} g=9.81 m/s2 C p , m = 1004   J / ( ℃ ⋅ k g ) C_{p,m} = 1004 \ \mathrm{J/(℃ \cdot kg)} Cp,m=1004 J/(kg) 代入可求得干绝热递减率为 γ d = 0.98  ℃ / 100 m = 0.98   K / 100 m \gamma_d = 0.98 \ ℃ / 100 \mathrm{m} = 0.98 \ \mathrm{K} / 100 \mathrm{m} γd=0.98 ℃/100m=0.98 K/100m。实际工作中取 γ d = 1  ℃ / 100 m \gamma_d = 1 \ ℃ / 100 \mathrm{m} γd=1 ℃/100m,这就是说,在干绝热过程中,气块每上升 100m,温度约下降 1℃。

【注意】 γ d \gamma_d γd γ \gamma γ气温直减率)的含义是完全不同的。 γ d \gamma_d γd 是干空气在绝热上升过程中气块本身的降温率,它近似于常数;而 γ \gamma γ 是表示周围大气的温度随高度的分布情况。大气中随地-气系统之间热量交换的变化, γ \gamma γ 可有不同数值,即可以大于、小于或等于 γ d \gamma_d γd

干绝热递减率还有另外一种计算方法。先由热力学第一定律的微分形式得:

C p , m d T = R T d p p C_{p, m} \mathrm{d}T = RT \frac{\mathrm{d}p}{p} Cp,mdT=RTpdp

两边同除以 d z \mathrm{d}z dz 得:

C p , m d T d z = R T d p p d z C_{p, m} \frac{\mathrm{d}T}{\mathrm{d} z} = RT \frac{\mathrm{d}p}{p \mathrm{d} z} Cp,mdzdT=RTpdzdp

引入定义式 γ d = − d T d z \gamma_d = - \frac{\mathrm{d}T}{\mathrm{d} z} γd=dzdT 和 静力学方程 d p d z = − g ρ \frac{\mathrm{d}p}{\mathrm{d}z} = -g \rho dzdp=gρ 可得:

C p , m γ d = R T g ρ p C_{p, m} \gamma_d = RT \frac{g \rho}{p} Cp,mγd=RTpgρ

整理得:

γ d = R T g ρ p C p , m \gamma_d = RT \frac{g \rho}{p C_{p, m}} γd=RTpCp,mgρ

利用气体状态方程 p = ρ R T p= \rho RT p=ρRT 得到:

γ d = R T g ρ ρ R T C p , m = g C p , m \gamma_d = RT \frac{g \rho}{\rho RT C_{p, m}} = \frac{g}{C_{p, m}} γd=RTρRTCp,mgρ=Cp,mg

这与我们第一次推导的结果一致。

5.5 湿绝热过程

5.5.1 潜热

在讨论是绝热过程前,需要先引入潜热的定义:指物质在等温等压情况下,从一个相变化到另一个相吸收或放出的热量。潜热分为两种:

  • 蒸发潜热:蒸发过程中,由于具有较大动能的水分子脱出液面,使液面温度降低。如果保持其温度不变,必须从外界吸收热量,这部分热量即为蒸发潜热。
  • 凝结潜热:与蒸发过程相反的是凝结过程。当水汽发生凝结时,这部分潜热又会全部释放出来,这就是凝结潜热。

对于绝热过程中的干空气和湿空气的潜热情况:

  • 干空气:气块没有液态水和固态水,在气块的绝热过程中水不发生相变,所以没有产生潜热;
  • 湿空气:气块中存在液态水和固态水,在气块的绝热过程中水发生相变,产生了潜热。更具体地讲,饱和湿空气绝热上升时,如果只是膨胀降温,亦应每上升 100m 减温 1℃。但是,水汽已经饱和了,就要因冷却而发生凝结,同时释放凝结潜热,加热气块。

显然,湿空气的绝热过程要比干空气复杂许多,因为凝结的水分可能会脱离气块,也可能随气块一起运动,不妨考虑两种极端情况:

  • 可逆湿绝热过程:水汽相变产生的水成物不脱离气块,随气块上升或下降,所释放的潜热全部保留在气块内部。这时无论气块上升或下降,温度都按湿绝热递减率而变化,这是与干绝热过程一样的可逆过程。
  • 假绝热过程:水汽相变产生的水成物全部脱离气块,此时气块下降时已不再是饱和湿空气,于是气块就按干绝热变化,这是一种不可逆过程。注意这时所释放的潜热仍留在气块中。

实际大气的湿绝热过程往往处于两者之间。不过,从热力学角度来看,这两者差别很小(作上升运动时,干绝热递减率与湿绝热递减率非常接近),我们一般用假绝热过程代替湿饱和气块的绝热过程。

5.5.2 假绝热过程

对于湿空气,设 1g 饱和湿空气中含有水汽 q s q_s qsg,绝热上升,凝结了 d q s \mathrm{d} q_s dqsg 水汽,所释放出的潜热被定义为:

d Q = − L d q s \mathrm{d} Q = -L \mathrm{d} q_s dQ=Ldqs

式中 L L L 表示水汽的凝结潜热,单位为 J / k g \mathrm{J/kg} J/kg。上式右边的负号表示当有水汽凝结时得到热量,因为这时水汽减少, d q s < 0 \mathrm{d} q_s < 0 dqs<0,则 d Q > 0 \mathrm{d} Q > 0 dQ>0;当水分蒸发时消耗热量,这时 d q s > 0 \mathrm{d} q_s > 0 dqs>0,则 d Q < 0 \mathrm{d} Q < 0 dQ<0

现在我们把上式代入到热力学第一定律中,可得:

− L d q s = C p , m d T − R T d p p -L \mathrm{d} q_s = C_{p, m} \mathrm{d}T - RT \frac{\mathrm{d}p}{p} Ldqs=Cp,mdTRTpdp

由于这个方程中只包含湿空气相变所产生的热量,而没有考虑其它的热量,所以上式又称为湿绝热方程。饱和湿空气上升时,方程又可写为:

d T = R T C p , m d p p − L C p , m d q s \mathrm{d}T = \frac{RT}{C_{p, m}} \frac{\mathrm{d}p}{p} - \frac{L}{C_{p, m}} \mathrm{d} q_s dT=Cp,mRTpdpCp,mLdqs

上式说明,饱和湿空气上升时,温度随高度的变化是由两种作用引起的:

  • 由气压变化引起,例如上升时气压减小, d p < 0 \mathrm{d} p < 0 dp<0,这使得温度降低;
  • 由水汽凝结时释放潜热引起,上升时水汽凝结, d q s < 0 \mathrm{d} q_s < 0 dqs<0,造成温度升高。

因此,凝结作用可抵消一部分由于气压降低而引起的温度降低。有水汽凝结时,空气上升所引起的降温将比没有水汽凝结时要缓慢。

5.6 湿绝热递减率(湿绝热直减率)

与干绝热递减率类似,定义下式为湿绝热递减率(用 γ m \gamma_m γm 表示):

γ m = − ( d T d z ) 湿空气块 \gamma_m = - \bigg( \frac{\mathrm{d}T}{\mathrm{d} z} \bigg)_{湿空气块} γm=(dzdT)湿空气块

对湿绝热方程两边同除以 d z \mathrm{d}z dz 得:

d T d z = R T C p , m d p p d z − L d q s C p , m d z \frac{\mathrm{d}T}{\mathrm{d}z} = \frac{RT}{C_{p, m}} \frac{\mathrm{d}p}{p \mathrm{d}z} - \frac{L \mathrm{d} q_s}{C_{p, m} \mathrm{d}z} dzdT=Cp,mRTpdzdpCp,mdzLdqs

引入定义式 γ m = − d T d z \gamma_m = - \frac{\mathrm{d}T}{\mathrm{d} z} γm=dzdT 和 静力学方程 d p d z = − g ρ \frac{\mathrm{d}p}{\mathrm{d}z} = -g \rho dzdp=gρ 可得:

γ m = R T C p , m g ρ p + L d q s C p , m d z \gamma_m = \frac{RT}{C_{p, m}} \frac{g \rho}{p} + \frac{L \mathrm{d} q_s}{C_{p, m} \mathrm{d}z} γm=Cp,mRTpgρ+Cp,mdzLdqs

利用气体状态方程 p = ρ R T p= \rho RT p=ρRT 得到:

γ m = g C p , m + L C p , m d q s d z \gamma_m = \frac{g}{C_{p, m}} + \frac{L}{C_{p, m}} \frac{\mathrm{d} q_s}{\mathrm{d}z} γm=Cp,mg+Cp,mLdzdqs

即:

γ m = γ d + L C p , m d q s d z \gamma_m = \gamma_d + \frac{L}{C_{p, m}} \frac{\mathrm{d} q_s}{\mathrm{d}z} γm=γd+Cp,mLdzdqs

当饱和湿空气上升时, d z > 0 , d q s < 0 \mathrm{d}z > 0,\mathrm{d} q_s < 0 dz>0dqs<0,则有 d q s d z < 0 \frac{\mathrm{d} q_s}{\mathrm{d}z} < 0 dzdqs<0;当饱和湿空气下降时, d z < 0 , d q s > 0 \mathrm{d}z < 0,\mathrm{d} q_s > 0 dz<0dqs>0,则有 d q s d z < 0 \frac{\mathrm{d} q_s}{\mathrm{d}z} < 0 dzdqs<0。所以, γ m \gamma_m γm 总是小于 γ d \gamma_d γd

需要注意的是,由于 d q s d z \frac{\mathrm{d} q_s}{\mathrm{d}z} dzdqs 是气压和温度的函数,所以 γ m \gamma_m γm 不是常数,而是气压和温度的函数。平均而言, γ m = 0.5 ℃ / 100 m \gamma_m = 0.5 ℃/100\mathrm{m} γm=0.5℃/100m

5.7 从 T-lnP 图看干绝热线和湿绝热线

这里先简单介绍一下 T-lnP 图是什么东西。

T-lnP 图,即温度-对数压力图,是一种在单站天气预报中使用的热力学图。其横坐标是温度( T ),以摄氏度(℃)为单位,每隔 10 ℃ 标出度数值,从图左端 -85 ℃ (188K)起向右递增,直至 40 ℃ (313K);纵坐标是气压(对数标尺),以百帕(hpa)为单位,从基准线 1000 hpa 向上递减,至 200 hpa。纵坐标可以简单理解为大气高度。注意!因为采用了对数尺度,所以纵坐标每 1 hpa 都不是等间距的。T-lnP 图反映了大气的许多热力学性质。

在这里插入图片描述

上图是干绝热线和湿绝热线的对比,使用了 T-lnP 图绘制。这幅图可以简单理解为:初始温度和气压相同的干空气块和湿空气块,在作上升运动时,气块内部的温度变化。干空气块上升时,温度是线性变化地下降;而湿空气上升时,温度先是非线性变化地下降,到达某处后再线性变化地下降。具体来说,湿空气经历的是如下过程:

  • 湿空气从 1000 hPa 高度抬升,进行湿绝热过程,水汽凝结释放潜热;
  • 由于我们视该绝热过程为假绝热过程,湿空气到达某处后已经没有水汽了,实质已变成干空气团,继续上升时会进行干绝热过程

为什么干绝热线和湿绝热线会是这样的结果?因为干绝热递减率 γ d = g C p , m \gamma_d = \frac{g}{C_{p, m}} γd=Cp,mg 是一个常数,故呈一直线;而湿绝热递减率 γ m = γ d + L C p , m d q s d z \gamma_m = \gamma_d + \frac{L}{C_{p, m}} \frac{\mathrm{d} q_s}{\mathrm{d}z} γm=γd+Cp,mLdzdqs 总是比 γ d \gamma_d γd 小,所以在干绝热线的右方。而且,下部因为温度高, γ m \gamma_m γm 小,上部温度低, γ m \gamma_m γm 大,这样形成上陡下缓的一条曲线。到高层水汽凝结愈来愈多,空气中水汽含量便愈来愈少, γ m \gamma_m γm 愈来愈和 γ d \gamma_d γd 相接近,使干、湿绝热线近于平行。

至于这个图是如何画出来的?后面介绍“位温”和“假相当位温”这个概念的时候会仔细说明。本节内容先到这里。

  • 14
    点赞
  • 10
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值