看起来,您的程序正在使用基于$q \ approapprow = w_s*RH$的近似值,并使用clausus - clapeyron近似值来找到$w_s$。看看RH,T和P的一些值,你的近似值与解析答案非常接近(+/- 5%)。根据您引用的输出,看起来您提供了不正确的RH值。注意在你的例行程序的评论中它说:
相对湿度(比例,不是%)
这意味着你需要提供RH比例,而不是百分比。例如,除以100——RH=1代表100%,RH=0.5代表50%,等等。
如果您调整了输入数据,您应该能够按原样使用代码。如果您希望将其与其他东西进行比较,您可以参考下面的解决方案。
如果你得到$RH$(范围[0,1]),$T$ (K)和$p$ (Pa),你可以按以下步骤进行。
知道$ $ RH = \ dfrac {e} {e_s}, $ $ $ $ w = \ dfrac {e \ R_d} {R_v(次)},$ $ $ $ q = \ dfrac {w} {w + 1} $ $
然后我们可以解出比湿度q。
与其将其组合成一个单独的公式并求解,不如将其增量地表示为更直接的方法。
首先,找到e_s美元(T),美元美元e_s (T) = e_ {s0} \ exp \离开\[左(\ dfrac {L_v (T)} {R_v} \) \离开(\ dfrac {1} {T_0} - \ dfrac {1} {T} \) \右]$ $,然后从第一个公式发现e美元(e = RH * e_s美元)。然后将$e$代入$w$的公式,然后将结果代入$q$的公式。
变量的使用:
比湿度或水蒸气与总空气的质量混合比(无量纲)
水蒸气与干空气的质量混合比(无量纲)
$e_s(T)$饱和蒸汽压(Pa)
$T_0$ (Pa)时的饱和蒸汽压
$R_d$干空气比气体常数(J kg$^{-1}$ K$^{-1}$)
$R_v$水蒸气比气体常数(J kg$^{-1}$ K$^{-1}$)
压力(Pa)
$L_v(T)$比蒸发焓(J kg$^{-1}$)
温度(K)
参考温度(一般为273.16 K) (K)