磁场作用下双扩散液层热毛细对流的数值模拟.docxVIP

磁场作用下双扩散液层热毛细对流的数值模拟

热盐双扩散对流的稳定性

在晶体生长和其他材料生长期间，祖母通常具有晶体和硬化的特征。在这一过程中，除了通过温度和铜质浓度引起的浮力对比外，从表面张力驱动的运动也对晶体的生长产生了重要影响。1981年，schwabe分析了关闭容器和开放容器之间的晶体生长差异，并解释了热纺织热作用对晶体生长的影响。因此，许多研究人员研究了单-双流的单-扩散源对矩阵熔渣的影响。thop等人和chen等人首次报道了在水平高差的作用下，在同一层中，温度和溶解浓度高的液体从温度的高度边界流动，而温度和溶解浓度低的液体从冷壁的表面流出。这种特殊的温度和浓度梯度的分布容易导致盐的不稳定性。之后，景观设计证实了热盐双扩展系统的非线性稳定性。结果表明，当表面的张力损失满足剪切条件时，矩形腔的液体流动也被分离。结果表明，当热态平衡性条件下的扩散热流的表面传递时，它对浮区的流场和浓度场的影响是明显的。王守元通过线性稳定性方法研究了波动损失和热态平衡性界面的基本特征。结果表明，由于表面张力的驱动作用，浮区长的扩散波的非常灵活，在很大程度上影响了浮区的热流和浓度。结果表明，由于主要区域的渐近解和数值模拟，大多数区域的渐近解和数值模拟结果具有良好的一致性。chen等人使用线性稳定方法研究了波动和表面积下两个薄液层的热流的两个扩散热流的表面传递。结果表明，由于主要表面积的逐步解和数值模拟，大多数领域的渐近解和数值模拟结果具有良好的平滑性。结果表明，由于主要表面积，离散扩散和表面积的热流的软格度很好，结果表明，庄内波的nmol-xiushi等人通过线性稳定方法研究了动态和静态傅立叶变换的波的不稳定性。结果表明，在很大的maranoni数下，提出了致动器热容量损失的初始、对称的不稳定性和无序扩散。然而，关于双重扩散热流的详细对称，没有文献充分解释。与此同时,通过引入外部磁场对热对流进行抑制在晶体生长过程中已被广泛采用.因此,本文在现有的双扩散热毛细对流研究的基础上,进一步对外部磁场作用下的双扩散热毛细对流进行研究,这对于了解磁场作用下的双扩散热毛细对流的流动特性及优化晶体生长过程都具有重要意义.

1液层中双扩散热毛细对流的标准显微控制方程

图1显示的是简化的研究模型,液层的厚度为d=0.01m,长度为L=0.04m.液池的左右壁面保持恒定的温度和溶质浓度,且右侧壁面的温度(T2)和浓度(C2)高于左侧壁面的温度(T1)和浓度(C1);自由表面和底面为绝热壁面,所有壁面都满足绝缘条件.施加的外部磁场为沿着y轴方向的均匀磁场.同时,假设流体为不可压缩的牛顿型黏性流体;液层的自由表面为水平且不发生变形;液层在Marangoni效应的作用下流动仍是层流;流体表面张力是温度和溶质浓度的线性函数,且表面张力随温度和溶质浓度的变化关系表示为

其中,γT为表面张力温度系数,γs为表面张力溶质浓度系数,初始表面张力σ0=σ(T0,C0).

基于以上假设,液层中双扩散热毛细对流的无量纲控制方程可以表示为

式中分别以U0=v/L和L作为特征速度和特征长度,V是速度矢量.各无量纲变量表示为:速度,温度,浓度坐标,时间,压力,诱导电流J*=J/(B0σmU0),磁场强度B*=B/B0.无量纲数为:Hartmann数,;Pr数,数,；磁雷诺数,Rem=μmσmU0L.α为热扩散系数,μ为动力黏度,v为运动黏度,D为质量扩散系数,μm为介电常数,σm为电导率.

流场和磁场相互耦合的控制方程可以表示为

其中B*为无量纲磁场强度,包括外加磁场和诱导磁场b*.诱导电流可以表示为:.

边界条件及初始条件:

式中,,.流体系统中的热Marangoni数表示为Ma=ReT·Pr,浓度Marangoni数表示为Mas=Res·Pr.热Marangoni效应和浓度Marangoni效应相对大小表示为.本文计算中所采用物性参数如表1所示.

2程序结果验证

采用非均匀交错网格对控制方程进行空间离散,其中在左右壁面、自由表面、底面附近对网格进行局部加密.采用有限体积法离散方程,对流项采用对流项的二次迎风插值(quadraticupstreaminterpolationforconvectivekinematics,QUICK)格式,扩散项采用中心差分,压力隐式分裂算子(pressure-implicitsplittingofoperators,PISO)算法用于求解N-S方程,时间项采用二阶隐式推进法.诱导磁场方程的离散和求解采用与动量方程相同的格式,动量方程中洛伦兹力的求解是采用上一个时间步的诱导电流和磁场强度.当速度、温度、压力等变量相对变化率小于10-5时认为计算收敛.本文计算区域划分为121(x)×51(y)个网格,网格的独立性已得到验证,计算中采用时间步长