5 超疏水表面冰核凝固特性数值研究.docVIP

超疏水表面冰核凝固特性数值研究 研究目的 模拟液滴在水平壁面上的凝固成冰核的过程，并分析壁面的浸润性对冰核形状、冰核与壁面接触角、凝固速率等因素的影响。 研究方案 现存的LBM相变模型主要有两种：基于多组分流动的相变模型（Mulity-componment method with phase change）[128-130]和基于焓法的相变模型（Enthalpy method with phase change）[131-132]。基于多相流动的相变模型不能模拟非等温的相变过程，不适合用于本文的研究。下面将主要介绍焓法的相变模型（Enthalpy method with phase change）[131-132] 基于焓法的相变模型 本文以液体凝固成固体的相变过程为例，给出基于焓法的相变模型的数学方程。 引入固相体积分数定义来描述固体的分布情况，定义如下： （1） 其中，为流体系统在温度下所对应的焓值，和所对应的温度和分别为液体开始凝固和凝固结束所对应的温度。考虑到凝固过程中的相变潜热，焓值的计算式如下： （2） 考虑相变潜热对温度场的影响，凝固过程中温度场的宏观控制方程可以写成如下形式： （3） 通过BGK近似以及泰勒展开，使用D2Q4格子，将宏观温度方程离散为LB温度演化方程。 （4） 式（1）、（2）和（4）构成了焓法相变的LB数学模型，显然这是一套封闭的方程组，包括三个方程和三个未知变量（固相体积分数，焓值，温度分布函数），求解这三个变量时，每一个时间步长下都需要进行数值迭代。假设时间为时刻的所有的固体体积分数、温度以及温度分布函数均已知，现需要求解时间为时刻时的变量。下面给出具体的数值求解方法。 （1）时刻第步的温度分布函数按下式计算， （5） （2）温度分布函数已知，计算时刻第步的温度； （3）时刻第步的焓值可以按下式计算： （6） （4）时刻第步的固体体积分数可按下式划分： （7） （5）重复上述（1）-（4）步，直至满足下面的收敛条件： （8） 其中，为判断收敛的小量，一般取。 使用上述基于焓法的LB相变模型模拟矩形空间内液体的凝固过程。物理模型如图1所示。初始状态下，矩形空间内充满静止的液体，温度为（高于或等于凝固温度）。四周为固体壁面，左壁面和下壁面保持为低温（低于凝固温度）。 图1 矩形空间内液体凝固过程数值模拟物理模型 Fig.1 Physical model for solidification of liquid in a rectangular space 计算区域长宽比为，采用D2Q4格子进行温度离散，网格划分为。设置液体初始温度，初始固相体积分数，下壁面和左壁面温度，凝固开始温度，凝固结束温度，比热为，相变潜热。四周壁面均采用非平衡外推格式处理。 模拟得到了不同时刻下温度分布和固相体积分数情况。温度沿着右上角至左下角的方向逐渐减低，凝固过程首先发生在流体区域左下角，然后随着时间的推移，固相界面逐渐向右上角方向推进，直至整个流体区域全部凝固成固体。 建立研究模型 格子Boltzmann凝固模型 本文所研究问题的物理模型如图2所示。假设在水平的壁面上方存在一球形饱和液滴，仅有底部直径最大的一点与壁面相接触。壁面保持恒定低温，温度低于冰点。液滴外部空间充满气体，初始温度为和壁面一样的低温。液滴温度不断降低，低于冰点时会凝固成冰。流体系统初始状态内两相（气、液）共存，当液滴逐渐凝固成冰核时，系统内三相（气、液、固）共存。模拟液滴在具有不同浸润性壁面上的凝固过程，分析壁面浸润性对凝固过程的影响是本文的研究重点。因此，模型在考虑相变过程对流场与温度场的影响的同时，还需要考虑气液两相之间的作用力以及固体壁面对气体、液滴和冰核的作用力。 图2饱和液滴在低温壁面上凝固过程物理模型示意图 Fig.2 Physical model of solidition of saturated droplet onto a cold solid surface 本文在基于焓法的相变模型，模型并没有考虑到相变过程对流动的影响，对其进行修正改进，使用固体体积分数修正流场演化过程中迁移步的粒子分布函数，并使用变粘度的处理方法，改进相变模型，使其更适合用来模拟高温液滴在低温壁面上的凝固过程。下面给出具体的数学模型。 使用伪势模型处理气液两相之间的作用力以及壁面对气体和液体之间的作用力。至于壁面对凝固而成的冰核的作用力，由于其过于复杂，本文暂不考虑。忽略浮升力的作用，使用伪势模型来求解流场，使用基于焓法的相变模型来求解温度场以及相变过程，需要说明的是，这里固相体积分数的定义仅适用于液相区域，对于气相流体，。 液滴