浮法玻璃摊平过程温度场的模拟研究.pdf

中国工程热物理学会 传热传质学 学术会议论文 编号：083414 浮法玻璃摊平过程温度场的模拟 张勤 陈泽敬 李志信 （清华大学 航天航空学院，热科学与动力工程教育部重点实验室，北京 100084） (Tel: 010806 Email: ZQ410@163.com) 摘 要：在分析浮法玻璃生产工艺及建立玻璃液摊平过程数学模型的基础上，综合考虑玻璃液物性与流 动特性，采用经验公式计算了玻璃液自由表面的高度与宽度，并预设玻璃液的实际活动空间，建立了 自由表面标高－SIMPLE计算方法，并采用该方法计算了摊平过程中的温度场。 关键词：浮法玻璃， 摊平，温度场 0 前言 玻璃液漂浮在熔融锡液表面，通过自身重力、表面张力、粘滞力及拉力综合作用完 成抛光与成型过程，称为浮法工艺。为了保证锡液不被氧化而污染玻璃，锡槽空间充满 保护气体，并保持微正压。玻璃液刚流入锡槽时，受到自身重力和液－液－气三相系统 表面张力的作用。随着玻璃液的流入，在自身重力影响下，玻璃液沿锡液表面摊开，并 在锡液面形成玻璃带成型源流[1]。摊平过程中，温度场对玻璃液的平整度和抛光效果有 直接影响，正确的温度场描述亦能对成品玻璃的缺陷进行预测。 目前，摊平及成型加工参数的调节依赖于工程人员的直觉经验，无法实现工艺参数 重复的一致性。浮法玻璃生产仅凭经验设计成为我国平板玻璃生产的软肋，已无法适应 浮法玻璃市场的发展与竞争。浮法玻璃成型工艺中玻璃液存在变形率大、材料物性高度 非线性、边界确定困难等特点，通过控制方程也无法获得解析解。因而数值模拟方法是 实现浮法工艺性能预测及产品升级的最好技术手段。 唐卫龙[2]使用二维有限元法，提出了回推模型，将玻璃熔体假设为弹性体，模拟结 果与实际情况有一定误差。A.milutinvic[3]在浮法工艺一维数学模拟基础上，考虑玻璃 物性的变化，模拟了玻璃带减薄过程。邵宏根[4]在忽略锡液摩擦力和玻璃液表面张力， 并认为玻璃宽度及厚度方向等温的假设下，用刚粘塑性理论建立了浮法成型过程的数学 模型，并采用 FIDAP 有限元软件模拟了玻璃带在拉边机下的减薄过程。但这些方法忽视 了摊平过程中玻璃液的流动性，并在物性、表面处理方面做了多种假设，无法对摊平过 程中的温度场进行准确计算。 鉴于此，本文采用 SIMPLE 方法对玻璃液摊平过程中的温度场和速度场进行模拟计 算。玻璃液在摊平过程中流动平缓，并且不与锡槽表面型壁接触，但与锡液、保护气体 存在运动界面，因此采用标高法来描述玻璃液的自由表面，即采用经验公式直接计算玻 璃液的轮廓，降低玻璃液摊平模拟中的计算量，并定义滑动对流边界条件，最后经 1 SIMPLE 算法计算，获得了玻璃液的最终温度场。 1数学模型与模拟方法 1.1 数学模型 资资助项目：国家科技支撑计划（2006BAF02A27 ） 锡槽的截面如图 1 所示，玻璃液下方与锡液直接接触，上方通常采用惰性或还原性 气体，以避免锡在高温下氧化而污染玻璃表面，造成质量缺陷。国内外一般使用 N +H ， 2 2 国内浮法生产线 N 一般占 92%～96%，H 占 4%～8%。保护气体直接与玻璃上表面接触， 2 2 其成份会影响玻璃液表面张力和表面粘附性。玻璃液与锡液、保护气体之间存在传热过 程，因此计算中需采用第三类边界条件确定玻璃液表面温度。 图 1 锡槽中玻璃液摊平简图 图2 玻璃液浸润情况 摊平过程中，流体粘度高，运动速度慢，Re 1 ，因此忽略惯性力，则玻璃液的数 学模型如式(1)-(3)所示； 连续方程 ∂ρg +∇(ρ v) 0