它看起来像你的程序是一个近似的基础上使用q \大约w = w_s * RH美元克劳修斯——克拉珀龙方程的近似w_s美元。我还没有测试看着几个值的RH, T和P,你近似是很接近(我将当我找到时间+ / - 5%),但这是来一个解析解决方案的答案。根据你报的输出看起来你是RH的提供不正确的值。注意在评论你的日常生活:
@param rh相对湿度(比例,而不是%)
这意味着您需要提供RH比例,而不是百分比。例如除以100,RH = 1 100%, 50% RH = 0.5等。
如果你调整你的输入数据应该能够使用你原有的代码。如果你想把它比作什么,你可以参考下面的解决方案。
如果你给出了RH美元(在[0范围,1)),元新台币(K)和$ p $(Pa)你可以进行如下。
知道$ $ RH = \ dfrac {e} {e_s}, $ $$ $ w = \ dfrac {e \ R_d} {R_v(次)},$ $和$ $ q = \ dfrac {w} {w + 1} $ $
然后我们可以解出特定的湿度问美元。
而不是结合成一个单一的公式和解决,这是更直接的呈现递增。
首先,找到美元美元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美元然后,结果到的公式问美元。
变量的使用:
问美元特定的湿度或水蒸气的大规模混合比总空气(无量纲)
w美元大量的水汽混合比干燥的空气(无量纲)
美元美元e_s (T)饱和蒸汽压(Pa)
美元e_ {s0} $饱和蒸汽压在T_0美元(Pa)
R_d美元干燥的空气(J公斤的特定气体常数$ ^ {1}$K$ ^ {1}$)
R_v美元具体为水蒸气气体常数(J公斤$ ^ {1}$K$ ^ {1}$)
$ p $压力(Pa)
美元美元L_v (T)特定的蒸发焓(J公斤$ ^ {1}$)
元新台币温度(K)
T_0美元参考温度(通常为273.16 K) (K)