相变界面流动技术的精度.pdfVIP

相变界面流动技术的精度 柯道友 清华大学热能工程系，北京。l∞O“ ● 有多种方式甩于楣变问患的分析中，其中界面藏动方法将系坑分捌成两个部分，在两 部分中分别末—羹量方程．在高热流情况下使用一阶有限差分模型模掇运动坐标jl簪会 ● 在燕蕾计算中引入凝差．一瞬中心差分将带来发敢问题，而一阶曩风差分也台在解中 引入误差． 本文分折了界面流动技术中模拟运动界面的方法，比较了一阶和二阶有限差分模拟运 动坐檬的慵提。圉■说明了高热流情况下一阶模型的局限性和使用二阶模型的必要性． 2．理论 在界嗣舞硪|技术中，目撼棱分割为固相和液相两千域．‘在分领j过程中，假定嗣相初始 置度为翱受■废．田爵仅需考察}蝴城，在边界条件为界面■度嗣定札的情况下求 表面量度固定在T_+AT．椴再醐运动由界面扮始量平 单元格尺寸和位置的变化；因此，需要籽麓量方程中的时 问差分项替换成全囊分。井减去代表移动网格的项，从而将单元格变换到同定网格坐 标中． 们’Dr田’搬 百2百一面i (1) 3．数值方莹 用二阶中心有限差分代替传导项，用cmk-N．啪tIs∞或是臆式格式代表瞬态项。 式(1冲用一阶和二输质风有限差分代替了温度对x的一阶导数项．这里使用顺风差分是 因为在此导数项甫有负号．由于系数为负。因而在某些迭代求解的情况下中心差分项 将不会收敛。一阶曩风差分为： ∞ 融 、出 7 一 罢=≥≠+等窘c鲁∞缸 2苏2 一阶晨风差分的误差项与能量方程同样包含二阶导散，因此使用一阶差分将在解中引 入发散项．在口i体力学文献中，广泛使用二阶顺风差分格式求解高精度问题，顺风方 向二阶差分为： 缸 2缸 、击 7 聊… 罢=墨鼍＆+竽窘(}o)3苏’ 负方向速度的公式相似． 假设在初始条件下，在极短时阿冉有极小区域的同相融解为液相，将这一区域分翻为N 个单元并进行计算。使用二阶界面流动技术．在l酽s的时间步长下，对璧面沮度超 堕融解量度l∞oc和lo∞oc的情况，分别给出融解全长的误差。以璧面韧始长度为自 本文由中国自然科学基金f，9辨5550-3)资助 172 变量．在初始长度大于箍长1％的慵醌下。最终误差随初始长度的I叫、而茂小，在韧始 长度小于总长l％的情况下，初始长度的制、不再对最鳝误差产生膨响．在麓解长度变 化的情况下融解时间相应变化，相应总的融解长度不对最终的误差百分值产生影响。 c棚止．N．啪Ison和隐式格式都是无条件收敛的，而显式格式仅在傅利叶散 Fo电6“6x，l汗仉5的情况下才会收敛．但是cmkjNj∞I湖和纛式咯式的精度并不符台 要隶。因为=阶导散井不连续．在初始情况下，大量的热量将从热蠼传人，面没有 热量从冷端传出，热靖最前靖的单元格将被加热封—个极高的假瀑度。如果时闻步长 的最终温度低于壁面温度，cr∞k-N耙ols∞方法收敛的条件为F昨l。求解结果表明当傅 利叶数大于0．5时，使用cr曲k二Nicol∞啾求解，在每—迭代时闻步中结果均在真值 周围振荡。当单元格和壁面间存在不同的温度差时均会发生这种振荡，如图2．但是当 温度的