液体真空蒸发冻结过程的动力学研究(2)

2010-01-12 张世伟东北大学机械工程与自动化学院

2.4、冷凝室的动力学模型

　　由真空室流入冷凝室内的气体,其中永久气体和少量的蒸汽被真空泵抽走,其余绝大部分液体蒸汽则被冷阱盘管所凝附。冷凝室内永久气体和液体蒸汽的压力控制方程分别为

　　由式(10)和(11)还可得出混合气体的总压力控制方程为

　　液体蒸汽在冷阱盘管上的凝结速率为

　　冷凝室内混合气体向真空泵流动的总质量迁移速率为

　　其中液体蒸汽和永久气体的各自质量迁移速率分别为

　　其中P _2a 、P_2W和P₂ —分别为冷凝室的永久气体压力、液体蒸汽压力和总压力,Pa ; P_i —固态溶液的饱和蒸汽压, Pa ; q_L2 —冷凝室的漏气量, Pa ·m³ ·s ^{- 1} ;q _2a 、q_2W和q₂ —分别为冷凝室流向真空泵的永久气体流量、液体蒸汽流量和总流量,Pa·m³·s ^{- 1} ; V₂ —冷凝室的有效体积,m³ ; qi—为冷凝管的液体蒸汽凝结速率,Pa·m³·s -¹ ; A_i —冷凝盘管的凝结表面积,m² ;T_i —凝结盘管的表面温度,K; S_e —真空泵在冷凝室出口的有效抽速,m³·s ^{- 1} ;α₂—液体在冷阱盘管上的凝结系数。

3、模拟计算与实验测试

3.1、模拟计算

　　依据上述模型,以纯水的真空蒸发冻结过程为示例,利用Matlab 软件进行了模拟计算。选择纯水为示例,主要是考虑其物性参数数据易于获得。模拟计算中,取水的汽化潜热λ2 = 2.499 ×10⁶J·kg - 1 ,比热CL = 4181J·(kg·K) -¹ ;冰的升华潜热λi和比热Ci 随温度的变化则利用热工数据表中的数据按照两温度间线性插值处理;水的饱和蒸汽压PW利用拟合公式计算。同时近似取q _L1 = 0 , q _L2 = 0 , QS= 0 ; Se 、CC MC 和u1 取为常数,它们的取值是在模拟计算过程中通过试算确定的。模拟计算中,所有动态参数均从所给定的初始条件出发,采用追赶法依次计算下一时刻的取值。计算时间步长取为1ms ,每隔10ms 保存一次数据。

3.2、实验测试系统

　　本文专门组建了实验系统进行实验测试,并比较了计算和测试的结果。实验装置是以GLZ-0.4型真空冷冻干燥机为主体,通过增设高精度、快速响应的温度、质量和真空度测量系统,专门改造而成。能够同步获取液体质量、温度以及真空室内真空度随时间变化的数据。

3.3、模拟计算与实验测试的结果比较

　　图2～图5 给出一组模拟计算与实验测试结果的对照图。模拟计算的结果与实验测试数据在变化趋势和量值上均很接近,表明所建立的数学模型和采用的模拟方法是正确的。质量、压力实验测试曲线上出现的尖峰值对应实验中液体发生暴沸的现象,造成电子称和真空计出现冲击读数。温度曲线在冰点处观测到过冷现象,这是数学模型中没有考虑因此无法模拟出的。