PLIC方法在浮法成型模拟中的应用.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：093421 PLIC 方法在浮法成型模拟中的应用 张勤 陈泽敬 李志信 （清华大学航天航空学院，热科学与动力工程教育部重点实验室，北京 100084） （电话：010 Email: ZQ410@163.com） 摘 要：通过分析浮法玻璃生产工艺，建立了玻璃液流动的数学模型，引入CSF 模型计算玻璃液运动界 面上的表面张力，并使用CTT 模型对物性间断进行连续性处理，提高了数值计算的鲁棒性。在此基础 上，采用PLIC 方法重构玻璃液、锡液、保护气体之间的自由界面，追踪玻璃液的运动。对匹兹堡公司 的LB 法锡槽入口的模拟表明，PLIC 方法能获得具有较高精细性的自由界面；LB 锡槽坎砖与槽底交接 处存在流动停滞区，容易造成玻璃液的积聚，降低锡槽使用寿命。 关键词：PLIC， 浮法成型，模拟 1 引言 我国浮法玻璃产量占世界第一，但由于产品结构不合理、质量低，优质浮法玻璃仍 然依赖进口。目前浮法玻璃工业正在从中低档的普通浮法向高附加值、低能耗的优质浮 法发展。成型过程是浮法玻璃质量的决定性因素，而锡槽是浮法玻璃成型的主要设备。 由于锡液的密度大于玻璃液，且玻璃液与锡液互不浸润、互无化学反应，玻璃液在锡槽 中漂浮于熔融锡液表面，完成摊平、展薄、抛光、冷却等过程。为避免锡在高温下氧化 而污染玻璃表面，通常采用惰性或还原性气体作为保护气体[1]。浮法玻璃成型数值模拟 技术由于成本低，产品开发周期快，能对产品质量进行预测，在工艺控制方面也起到重 要作用，因此日益受到重视。 整个成型过程中，玻璃液下部与锡液直接接触、上部与保护气体直接接触，涉及到 移动边界的流动问题。移动边界的准确定位是预测浮法玻璃厚薄差、成品率、渗锡量等 参数的前提。尤其是玻璃液刚进入锡槽时，由于本身温度高(约1100℃)，流动性较好， 玻璃液、锡液、保护气体间的界面模拟是整个成型过程精确模拟的难点[2]。VOF(Volume of Fluid)方法是目前常用的自由界面跟踪技术，该方法只需将运动界面在空间网格内定义 成流体体积函数，在此基础上构造流体体积函数的发展过程，从而确定该运动界面的位 置、形状和变形方向。PLIC(Piecewise Linear Interface Construction，分段线性) 是VOF 方法的一种，在单个网格内，采用直线段模拟近似界面，具有较高的精度且计算 [3] 量小 。 鉴于此，基于N-S 方程建立了玻璃液流动的数学模型，针对玻璃液与锡液、保护气 体之间的自由界面，引入Youngs 具有二阶精度的PLIC 自由界面重构技术，编制成型过 程的流动与界面重构程序，模拟玻璃液的流动，并追踪玻璃液的运动界面。为了提高模 拟精度，程序中使用CSF(continuous surface force)模型处理表面张力。为了克服界 面上不同流体工质物性的不连续性，采用CTT(连续变化处理)模型。运用该方法对处于 LB 法锡槽入口的玻璃液自由面进行了追踪。1 1资助项目：国家科技支撑计划（2006BAF02A27 ）和973 项目（2007CB206901 ） 2 数学模型 2.1 PLIC 方法 VOF 方法的基本原理是通过研究网格单元内流体体积函数 来确定自由面，从而 F 追踪流体的变化，VOF 方法可以处理自由面重入等强非线性现象，计算时间短、存储量 小。VOF 方法根据体积函数来构造和追踪自由面，体积函数F 表示网格内流体的填充 状态，若F 0 ，则该网格为无指定相流体单元，若F 1，则说明该网格内充满了指 定相流体，当0 F 1，则该网格为交界面单元。两维体积函数以守恒形式的传输方 程表示为， （ ） （ ） ∂F ∂ uF ∂ vF + + 0