聚合物全三维非稳态非等温多相分层流动成型过程的理论模型和数值.PDF

第 40 卷第 9 期 机 械 工 程 学 报 Vol.40 No.9 2004 年 9 月 CHINESE JOURNAL OF MECHANICAL ENGINEERING Sep. 2 0 0 4 聚合物全三维非稳态非等温多相分层流动 成型过程的理论模型和数值模拟* 周国发 孙 懋 (南昌大学环境科学与工程学院 南昌 330029) 摘要：基于聚合物多组分成型技术的工程背景，建立了全三维非稳态非等温多相分层充模流动的理论模型。在综 合分析了该理论模型产生数值解的不稳定及发散的主要原因基础上，提出了求解理论模型的稳定快速收敛的数值 算法。采用罚函数法，以及速度场分析、温度场分析和流体体积分数分析相分离等方法来降低对计算机的 CPU 和 存储能力的需求，而通过 SUPG 法、罚函数法和 + 三维单元实现有限元数值分析的稳定性。并基于罚函数法 P1 / P1 和 SUPG 法，推导出求解 N-S 方程、能量方程和全三维多相分层流动成型移动前沿界面和分层界面追踪方程的有 限元数值模型，并探讨了全三维非稳态非等温多相分层流动移动前沿界面和分层界面重构技术。最后揭示了粘性 包围和界面不稳定的产生机理。 关键词：多相分层流动 多组分成型技术 数值模拟 粘性包围 界面不稳定 中图分类号：TQ320 合理的多相分层流动成型过程的理论模型必须综 0 前言 1 合考虑粘性包围和界面不稳定相互作用，而要研 究并模拟粘性包围和界面不稳定等关键问题，必 多组分成型技术广泛应用于多层复合基材等高 须是基于全三维非稳态非等温多相分层流动成型 性能复合塑料制品的成型加工，其典型成型工艺如 过程的理论模型和数值算法，为此提出了聚合物 共注成型、气辅共注成型等。其成型过程均属全三 全三维非稳态非等温多相分层充模流动成型的理 维多相分层流动成型过程，而成型机理极为复杂， 论模型和数值算法，并依此研究了聚合物多相分 涉及粘性包围、界面不稳定和不平衡流动模式以及 层充模流动成型过程产生粘性包围和界面不稳定 粘弹性等多因素作用耦合。迄今，有关多相分层充 的产生原因。 模流动成型的机理研究尚不完善。聚合物熔体多相 分层流动成型技术问题主要来源于不同流变性能和 1 数学模型 物理性能的聚合物熔体在成型过程中所引起的粘性 包围和界面不稳定两种特有现象，这两种现象必然 1.1 控制方程 引起多相分层流动层厚再分布，并产生三维波动分 本数学模型是基于流体体积法，并认为各变量 层界面，这种不理想的层厚分布不仅严重影响最终 场(压力场、速度场和温度场等)均由多相分层流动 成品的性能，而且影响多相分层流动成型的生产效 成型的各相通过体积加权平均来共享[1～3] 。为此， 率[1]。由此可见，为了建立聚合物熔体多相分层流 确定整个计算区域各相的体积分数ϕ 极为重要。ϕ α α 动成型技术的理论模型和成型工艺的设计原则，必