来源公众号:COMSOL仿真交流
COMSOL
建模仿真
马铃薯在冷冻干燥过程中主要发生热传递和湿传递,采用 COMSOL Multiphysics 对马铃薯切片进行模拟,综合利用传热学、传质学和有关学科知识,建立马铃薯切片在升华干燥过程中的物理模型和数学模型,模拟马铃薯片的升华干燥过程。
+ ✦ +
物理模型
假设马铃薯切片在冻干机的环境下进行升华干燥模拟。市场上选择的新鲜马铃薯的初始含水率 81.85%左右,马铃薯的共晶点温度为-22.5℃,一般预冷冻温度低于共晶点 5-10℃,所以选取马铃薯切片预冷冻温度为-30℃。为了加快马铃薯的干燥,在不影响马铃薯冻干质量的情况下把冻干箱隔板温度提高到-10℃,环境温度提高到 30℃。以减少大量升华阶段的时间。
考虑到模拟是连续性模拟,马铃薯内部具有各向同性,因此选取马铃薯圆柱一个过圆心的半截面来模拟马铃薯升华干燥过程中传热传质 。
图1采用 COMSOL建立马铃薯切片物理模型,在模型开发器下面的几何中建立直径 50mm高 10mm 的二维轴对称瞬态模型。同时加入了达西定律、多孔介质的传热、变形几何接口等。并在上述接口中建立了模型的初始状态和边界条件。用变形几何跟踪升华界面,计算耦合的质量平衡和热平衡。升华干燥模拟采用了干燥层和冷冻层方法,干燥初先假设切片有一个微小的干燥层 ,厚度为0.05mm,并假设升华产生的蒸汽全部进入干燥层并通过上表面流出圆柱。
如图1所示,上层为干燥层,下层为冷冻层,中间那条线代表升华界面。
+ ✦ +
数学模型
为了简化模型,做出如下假设:
(1)将马铃薯切片假设为各向同性连续材料,假设马铃薯切片内部的初始温度和冰含量是均匀分布的。
(2)升华界面上水蒸气的浓度与冰平衡。
(3)冰升华后,形成的固体基质是可渗透的,允许蒸汽通过。综合利用传热学、传质学和有关学科知识,建立马铃薯切片在升华干燥过程中传热传质的数学模型。
1、传热方程
在冷冻层没有对流,只有热传导,传热方程如下:
为冷冻层马铃薯切片的有效密度(kg/m3),𝐶1,𝑝为冷冻层马铃薯切片有效比热(J/(kg·K));λ1为冷冻层马铃薯切片的有效导热系数(W/(m•K));T为马铃薯切片的温度(K),t 为干燥时间。
干燥层因蒸汽的扩散和渗透有传热和对流。传热方程如下:
为马铃薯骨架和水蒸气混合的有效密度(kg/m3),𝐶2,𝑝为马铃薯骨架和水蒸气混合的有效比热(J/(kg·K));λ2为马铃薯骨架和水蒸气混合的有效导热系数(W/(m•K));𝑝𝑣为蒸汽密度 (kg/m3) , 𝐶𝑝,𝑣 为蒸汽比热(J/(kg·K));u 是根据达西定律得到的速度。
在蒸汽域中,密度是用蒸汽的摩尔质量𝑀𝑣和理想气体定律计算的:
其中 R 是通用气体常数。假设相变边界上的热力学平衡,则升华前沿的温度𝑇𝑠如下:
其中𝐿𝑠是升华潜热,𝑝𝑣是相变边界上的蒸汽压。升华界面处的干燥层和冷冻层温度相等,两边热流通过以下条件连接:
其中𝑉𝑠是升华界面速度,(m/s);𝑁𝑣,𝑛为垂直升华方向的蒸汽通量,g/ (m2·h);𝑄𝑠是界面处法向热通量的突变,J。
2、传质方程
仅在干燥层有蒸汽的扩散和渗透,因此马铃薯切片传质仅发生在干燥层。传质方程为:
式中ε是产品的孔隙率;𝑝2,𝑣为干燥区水蒸气分压,Pa;ρ2,𝑝为干燥区多孔介质的密度,kg/𝑚3;
根据质量传递和达西定律,能得到通过孔隙的蒸汽流率𝑄𝑚和质量通量𝑁𝑣:
其中ρ𝑣是蒸汽密度,kg/𝑚3; μ𝑣是蒸汽动力黏度, Pa∙s; 𝑝𝑣是蒸汽分压,Pa;k 是干燥产品渗透率, 𝑚3;ρ𝑖𝑐𝑒是冰的密度,kg/𝑚3。
3、初始条件和边界条件
初始条件:
t=0,T=𝑇𝑖= 243.15𝐾,p=𝑝𝑐= 24𝑃𝑎
边界条件:
t> 0,0<z≤Z0,𝑇α=303.15K,𝑇𝑠=263.15K
λ2*𝜕𝑇/𝜕𝑡 + 𝐿𝑠𝑄𝑚 = 𝑞𝑟 + 𝑞𝑐 + 𝑞𝑑
式中:p为马铃薯表面压力,Pa;z为升华界面位移,cm;𝑞𝑟、 𝑞𝑐、𝑞𝑑分别为辐射、对流、导热换热量,J。
4、热物性参数
马铃薯的导热系数计算公式如下:
式中:𝑀𝑖𝑣为马铃薯各成分的体积分数,%;𝑀𝑖马铃薯各成分的质量,g; ρ𝑖为马铃薯各成分的密度,kg/𝑚3; λ𝑖为马铃薯各成分的导热系数,W/(m∙K)。连续性模拟所需的主要物性参数如表 1所示。
+ ✦ +
模拟分析
升华干燥总时长为 4.8383h。冰质量与初始时刻冰质量之比随着干燥时间增加逐渐减少,脱水速率随着干燥时间增加逐渐减小。考虑到模拟忽略了蒸汽从切片四周流出的实际情况,实际上升华干燥能除去马铃薯切片85%以上甚至 90%以上的水分,干燥效果较好。
投稿专栏
“COMSOL仿真交流”微信公众号自创办以来得到了广大科研工作者和研究生的广泛关注和支持。为更好地服务仿真学习与研究,本公众号现因业务需要长期招聘供稿作者。您可发送简历到likeapoem@163.com,邮件主题请注明“姓名+供稿作者”。
欢迎“光声力热”类相关专业的科研爱好者加入“COMSOL仿真交流”团队。