Skip to content

Latest commit

 

History

History
54 lines (35 loc) · 1.24 KB

README.md

File metadata and controls

54 lines (35 loc) · 1.24 KB
Raw

Psychrometric

Mathematica语言实现的湿空气状态点计算程序

Example

绘制焓湿图

焓湿图

根据任意两参数求状态点

变量 含义
t 干球温度(°C)
d 含湿量(g/kg)
φ 相对湿度(%)
h 焓(kJ/kg)
ts 湿球温度(°C)
tl 露点温度(°C) 
pt1=psy[t->25,\[CurlyPhi]->100]
pt2=psy["h"->77,"\[CurlyPhi]"->80]


(*<|t->25,d->20.0829,\[CurlyPhi]->100,h->76.2561,ts->25.,tl->25.|>
<|t->28.0022,d->19.1485,\[CurlyPhi]->80,h->77,ts->25.1769,tl->24.2275|>*)

求两状态点混合

ratio=75;(*混合比例*)

psy["h"->(pt1["h"]*ratio+pt2["h"]*(1-ratio))/100,"d"->(pt1["d"]*ratio+pt2["d"]*(1-ratio))/100]


(*<|t->-2.00542,d->0.892314,\[CurlyPhi]->28.0496,h->0.212045,ts->-5.6097,tl->-16.4011|>*)

求送风状态点

\[CurlyEpsilon]=10000;(*热湿比(J/g)*)
\[CurlyPhi]l=95;(*机器露点*)

psy@@FindRoot[{dp[\[CurlyPhi]l/100Pqb[t]]==df,d==df}/.
NSolve[{(hf-pt2["h"])/(df-pt2["d"])==\[CurlyEpsilon]/1000,htd[t,df]==hf},{hf,df}][[1]],
	{{t,pt2["t"]},{d,pt2["d"]}}]


(*<|t->24.688,d->18.686,\[CurlyPhi]->95.,h->72.3752,ts->24.0545,tl->23.8324|>*)