环形双层液体内热毛细对流过程的数值模拟

为了了解环形双层液体内热毛细对流的基本特征,利用有限容积法进行了数值模拟,液池外壁被加热,内壁被冷却,液封流体为氧化硼,熔体为砷化镓,其宽深比为4.模拟结果表明,液封能极大地抑制熔体内流动,使其速度降低约两个数量级;浮力会削弱液封层内的流动,强化熔体内流动.     

Level Set方法及其在两相流数值模拟研究中的应用

本文介绍了一种崭新的捕捉两相流相界面的 Ｌｅｖｅｌ Ｓｅｔ方法,并应用该方法对两种互不相容流体间的平面界面波动现象进行了数值模拟,捕捉到了两相流界面波动过程中的一些复杂现象。计算结果与实验观测现象符合良好,表明这种方法是一种很有前途的两相流数值模拟方法,很值得进一步研究和推广。     

环形池内双层流体热毛细对流不稳定现象实验研究

为了了解水平温度梯度作用下双层流体内热毛细对流的特征,实验观察了环形池内由1cSt硅油和氟化液FC-40组成的双层流体内热毛细对流失稳后的不稳定现象。实验表明,随着上层或下层流体厚度的增加,流动稳定性减弱。当上层厚度为1.5 mm时,内外环间温差超过临界温差后在外环区域产生三维稳态的直条纹、并在内环区域产生热流体波,波纹产生于上层流体。当上层流体厚度为2mm时,液层内三维稳态直条纹占主导。环形池旋转使稳态条纹变为弯曲波纹,并且液层内只有一组传播方向与环形池的旋转方向相反的热流体波。     

大尺度环形液池内双层热毛细对流不稳定性

假设双层流体的交界面不发生变形,热毛细力作用于此交界面,三维数值研究了大尺度环形液池中双层流体系统在内外壁面温差加热下的热毛细对流不稳定性,其中外壁面维持高温,内壁面维持低温。计算结果显示,上下层流体的流动特性受Marangoni效应和浮力效应的影响;热毛细对流的振荡产生于内壁面附近,并沿着温度梯度的方向传播;随着温差的增大,热毛细对流的振荡逐渐增强,温度振荡波数增大。     

界面捕捉Level Set方法的(AMR)数值模拟

在流体力学方程的计算中采用高精度WENO格式,用AMR(adaptive mesh refinement)方法提高流场局部分辨率,在采用Level Set函数标定物质界面的计算中用GFM(ghost fluid method)方法进行界面处理,尝试将AMR技术与界面追踪技术相互融合并应用于数值模拟,对不同的模拟结果进行了比较.     

深径比对双层热毛细-浮力对流的影响

三维数值研究了环形液池内双层流体热毛细-浮力对流的流动特性,考察了液池深径比对流动的影响,其中外壁面维持高温,内壁面维持低温。结果显示,上、下层流体热毛细-浮力对流的流动结构依赖于液池深径比的大小,其内在机制是由浮力和热毛细对流相对强度所决定;随着深径比的减小,热毛细-浮力对流不稳定性的周向温度波数先增大后减小,存在发生流动不稳定性的最佳深径比;热毛细对流的振荡周期随着深径比的减小逐渐减小,而周向速度平均扰动幅度逐渐增大.     

绕汽柱热毛细对流的数值模拟

本文建立了一个在常重力及微重力情况下绕汽柱热毛细对流的数学模型,并采用有限差分法来积分非稳态动量方程和能量方程,得到了不同参数下热毛细对流的数值解．通过数值模拟,得到了不同参数下热毛细对流的速度矢量场,了解了有关无量纲参数对热毛细对流的影响规律．得出结论：热毛细对流不仅在微重力条件下存在,即使在重力场中,在特定情况下也是不能忽略的。     

水平磁场对双层流体热毛细对流的影响

三维数值研究了零重力时水平温度梯度作用下,B2O3封闭液与InP熔液组成的双层流体系统在水平磁场作用下的热毛细对流.结果显示,当磁场强度较小时,上层流体中对流涡的结构变化其微,而下层流体中的逆对流涡得到了抑制;随着磁场的增强,温度场分布逐渐趋于均匀,涡量强度逐渐减小,表明热毛细对流强度逐步地得到削弱.当磁场增强列Bx=0.4 T时,上层流体中的对流涡较均匀地充满上层流体区,而下层流体中对流涡紧贴着交界面.     

液封液桥内振荡热毛细对流的三维数值模拟

为了了解微重力下液封液桥内热毛细对流的基本特性,利用有限差分法进行了非稳态三维数值模拟,液桥高为(1-3)mm,直径为2mm和3 mm,液封外直径为(4-7)mm。模拟结果表明,当Marangoni数较小时,液封液桥内的热毛细对流为稳定的轴对称运动,当Marangoni数超过某一临界值后,流动将转化为三维振荡流动;为此,确定了发生振荡的临界Marangoni数,分析了各种条件下热毛细对流的振荡特性,计算了相应的振荡频率。     

环形浅液池内中等Pr数流体的热毛细对流

为了了解微重力下水平温度梯度作用时环形浅液池内的热毛细对流特性,利用有限差分法进行了非稳态二维数值模拟,环形液池外壁被加热,内壁被冷却,流体为 0.65 cSt的硅油,其Pr 数为 6.7。结果表明,当温度梯度较小时,流动为稳态流动,随着温度梯度的增加,流动将会失去其稳定性,转化成各种振荡流动,模拟结果与实验结果基本吻合。     

微重力下两同轴液柱内的热毛细对流的渐近解及数值模拟

1引言获得高质量、大体积单晶体对理论研究和实际工业应用都极其重要。由温度梯度引起的表面张力梯度所诱发的热毛细对流对晶体生长不利。zer等人山的空间实验证实，时相关的热毛细对流将引起熔质分凝，带来杂质条纹。人们急需有效抑制热毛细对流的方法。近年发展起来的液封技术在抑制热毛细对流方面表现出诱人的前景［‘－’］。本文着重研究两不相溶混的轴对称同轴液柱内的热毛细对流，通过分析求解及数值模拟，找出了外层流体对内层流体运动的抑制规律，为浮区晶体生长提供一定的理论和实践指导。2物理模型和数学模型如图1所示，距离为L…     

环形浅液池内热毛细对流的热力学特性

为了了解径向温度梯度作用下环形浅液池内硅熔体热毛细对流的热力学特性,利用有限差分法进行了非稳态三维数值模拟。液池外半径r0=50 mm,内半径ri=15 mm,深度为d=3 mm。结果表明,当温度梯度较小时,流动为稳定轴对称流动,系统总熵产较小;随着温度梯度的增加,流动将失去其稳定性,首先转化为径向脉动波,此时系统总熵产呈周期性变化;温度梯度再增加时,流动转化为热流体波,系统总熵产较大,但不再随时间变化。     

硅熔体CZ结构浅池内热毛细对流转变滞后特性

为了了解水平温度梯度作用时Czochralski(CZ)结构浅池内硅熔体热毛细对流的转变滞后特性,利用有限差分法进行了非稳态三维数值模拟,坩埚外壁被加热,半径为50mm,晶体半径为15mm,液池深度为3mm,坩埚外壁与晶体生长界面温差变化范围为6-27K。模拟结果表明,当逐渐增加温差时,在△T=9K处,二维轴对称流动转变为三维稳态流动,在△T=20.6K处,三维稳态流动转变为三维振荡流动;当逐渐减小温差时,在△T=19.5K处,三维振荡流动才转变为三维稳态流动,因此,二次流动转变存在滞后,滞后温差约为1.1K。     

浮力对环形液池内热毛细对流的影响

为了了解浮力的影响,对水平温度梯度作用时环形液池内的热毛细对流进行了非稳态三维数值模拟,环形液池外壁被加热,半径为40 mm,内壁被冷却,半径为20 mm,液池深度为(1-17)mm,流体为0．65cSt的硅油,其Pτ数为6．7．模拟结果表明,当水平温度梯度较小时,流动为轴对称稳态流动,随着温度梯度的增加,流动将会失去其稳定性,在浅液池内,转化为热流体波,浮力对失稳后的流型无影响,但会使热流体波的振幅下降;在深液池内,在常重力条件下,转化成三维稳定流动,在微重力和小重力条件下,转化为三维振荡流动．     

具有隔板的平行通道内空气混合对流换热数值模拟

本文以锅炉干排渣装置为背景,对抽象的理论模型具有隔板的平行通道内空气混合对流换热进行了数值模拟.数值计算表明,在Re＞1000时应采用非稳态数学模型进行数值模拟;在Re＞500时,自然对流机制对流动和换热的影响基本可以忽略.数值计算给出了不同Re时的进出口无量纲压差、局部的Nux和平均Nu以及流线图.这些结果可为深入研究干排渣装置中流动和换热特性提供参考.     

分离结晶过程中熔体热毛细对流的转变

为了了解微重力条件下新型分离结晶生长过程中熔体热毛细对流的基本特征,利用有限差分法进行了数值模拟,熔体深径比A取1和2,自由界面无因次宽度B分别取0.05、0.075和0.1.当熔体上表面为自由表面时,得到了分离结晶Bridgman生长过程中熔体热毛细对流的流函数和温度分布.计算结果表明:当Ma数较小时,在上下两个自由表面的表面张力的驱动下,熔体内部产生了两个流动方向相反的流胞,流动为稳态流动,随着Ma数的增加,上下自由表面速度增大,温度分布的非线性增加;当Ma数超过某一临界值后,流动将转化为非稳态流动;与熔体上表面为固壁时相比,A=1时的临界Ma数减小,而A=2时的临界Ma数增大;流动失稳的物理机制是流速的变化和阻力的变化之间存在滞后.     

方形空间内混合对流换热的数值研究

以建筑物内人工环境控制为应用背景,对有对称空气射流的方形空间内混合对流换热进行了数值模拟,探讨了这种具有对称结构的混合对流换热解的分岔问题。数值结果表明,Ｒｅｙｎｏｌｄｓ数、强制通风气流的射流角度、以及方形空间的宽高比都会影响解的分岔。当Ｒｅ数超过某一临界值时,会出现非对称数值解。宽高比减小,出现非对称解的临界Ｒｅ数也随之减小。Ｒｅ数、宽高比一定,仅当通风气流的射流角度在某个范围内时,能够得到非对称的数值解。