第六章果汁冷冻浓缩过程冰晶生长的相场法模拟2.docVIP

第六章 果汁冷冻浓缩过程冰晶生长的相场法模拟 0. 引言 上一章对过冷纯水的冰晶生长过程进行了相场模拟，并证实了模拟的可行性。本章利用相场法对果汁冷冻浓缩过程进行模拟。众所周知，不同的果汁其组成成分不同，且任何果汁中都含有多种组成物质。为了简化处理，把果汁简单视为只含有糖和水两种成分，即探索糖和水的二元物质冷冻浓缩过程的冰晶生长的数值模拟。 二元糖水的冰晶生长的相场模型 1.1 相场控制方程 利用第五章相场模型建立的基本原理，通过采用kim等提出的模型并应用稀溶液近似时的基础方程，提出了关于糖水的相场法的相场模型。相场法是利用自由能密度的形式导出相场方程和扩散方程的。设、分别为固相和液相的自由能密度，且其固、液相分数分别为和，为界面的过剩自由能，则固液混合物的自由能密度可以由（6-1）式表达： （6-1） 其中是溶质的摩尔分数，；。由此可得到以下相场方程： （6-2） 是（6-1）式对的偏微分，当利用稀溶液近似时，其表达式为： （6-3） 式中：与界面动力学有关的相场参数；为相场参数，为自由能密度，为式（6-1）对的偏导数；为气体常数；为温度；为摩尔体积；为平衡浓度；上、下标和分别表示固相和液相。 1.2溶质扩散方程 对于扩散方程，仍利用自用能密度的形式可改写成下面的形式： （6-4） 式中为溶质扩散系数，的下脚标、分别为浓度的一阶、二阶偏导数。另外，在数值计算中，不是将（6-4）式差分化，而是先计算括号内的溶质通量成分，然后利用其和求出各单元内的浓度变化。这时应用稀溶液近似处理，（6-4）式的右边括号内的各项分别如（6-5）、（6-6）所示。 (6-5) (6-6) 式中，为气体常数，T为温度，为摩尔体积。 在界面区域的溶质浓度是固相和液相的质量分数的和，即： （6-7） （6-8） 相场参数的确定 相场参数和与界面能及界面厚度有关，则与界面动力学系数有关，它们可以表示为：， （6-9） （6-10） 式中，为界面能，为界面厚度，为动力系数，为气体常数，为温度，为摩尔体积，为液相线斜率，为平衡常数。 本章考虑的是等温状态。等温模拟忽略了冰晶生长过程潜热的释放，而将过冷液的温度场视为在空间和时间上恒定不变的体系，从而使其相场模型中只存在相场控制方程和溶质传输方程的藕合。这种简化处理称为等温近似。利用这种近似，来初步研究过冷糖水溶液中冰晶的生长行为。 1.4 数值计算 1）初始条件和边界条件 设初始晶核半径为，以下列条件来确定初始条件和边界条件，即 （6-11） 式中：和分别分别是模拟中的横坐标和纵坐标；为初始温度。在计算区域边界，相场和溶质场边界条件采用Zero-Neumann边界条件，即 （6-12） 2）数值计算方法 采用有限差分法将方程（6-2）和（6-4）求解，均匀剖分网格，空间步长，时间步长为，和满足以下关系： （6-13） 冰晶长大过程的形貌分析 图6-1显示了迭代时间增加情况下糖水冰晶二维形貌演化图。图形的内部区域是固相，外部区域是液相，两者之间的过渡区域是固、液界面区域。图（a）为无量纲时间，冰晶小，四周比较光滑，晶臂界面稳定；图（b）为无量纲时间，半径增大，冰晶尖端有些凹陷，呈枝状；图(c)为无量纲时间，影响程度加大，冰晶枝状明显，枝端开始变长、变细，枝晶尖端半径减小，生长速率增大；图(d)为无量纲时间，冰晶枝状明显增大，枝端继续变长、变细。 如果模拟时间继续增大，将会出现二次枝晶臂,而且还可能出现三次枝晶臂。这个模拟结果与第五章纯水的模拟结果一致，很好地解释了在冷冻浓缩过程中，冷冻时间对冰晶生长的影响，控制适当的冷冻时间可以有效地降低冰晶的夹带率。 （a） （b） （c） （d） 图6-1不同时间下的冰晶形貌演化图 Fig.6-1 Ice crystal pattern evolvement corresponding to different time 图6-2 是冰晶生长过程中温度对冰晶尖端生长速率和曲率半径的影响曲线。 从图中可以看出，随着温度的升高，冰晶尖端的生长速度减小，而冰晶尖端的曲率半径增大，基本上是近似的线性变化关系。还可以看到，在不同温度下的冰晶主干的生长长度是不同的，随着初始温度越低，其长度越长，即枝晶的生长速率越快。这表明过冷度大有利于枝晶的生长，这和凝固理论中