化工学报 ›› 2016, Vol. 67 ›› Issue (7): 3047-3054.DOI: 10.11949/j.issn.0438-1157.20151900

• 材料化学工程与纳米技术 • 上一篇    下一篇

注塑成型Hele-Shaw流动模拟中热对流的异步长求解

王超房, 黄明, 石宪章, 申长雨, 赵振峰   

  1. 郑州大学橡塑模具国家工程研究中心, 河南 郑州 450002
  • 收稿日期:2015-12-14 修回日期:2016-03-19 出版日期:2016-07-05 发布日期:2016-07-05
  • 通讯作者: 黄明
  • 基金资助:

    国家自然科学基金重点项目(11432003);河南省高等学校重点科研项目(15A430009,15A430045)。

Thermal convection calculation with variable time step in Hele-Shaw flow simulation of injection molding

WANG Chaofang, HUANG Ming, SHI Xianzhang, SHEN Changyu, ZHAO Zhenfeng   

  1. National Engineering Research Center for Advanced Polymer Processing Technology, Zhengzhou University, Zhengzhou 450002, Henan, China
  • Received:2015-12-14 Revised:2016-03-19 Online:2016-07-05 Published:2016-07-05
  • Supported by:

    supported by the Key Program of National Natural Science Foundation of China (11432003) and the Key Research Project for Henan Universities (15A430009, 15A430045).

摘要:

在Hele-Shaw流动数值模拟中,速度是压力的后处理结果。如果是点浇口,则浇口附近速度会随单元尺寸缩小而趋于无穷大,导致能量方程作为一个整体求解时,时间步长必须非常小,否则会产生很大误差;而根据热对流物理意义分步求解,则需追踪当前物质在上一时刻位置,当单元速度很高、逆向搜索需穿透多个单元时,搜索可能会失败。鉴于此,基于分步求解法,研究提出一种变长度子时间步长方法处理对流项,确保搜索路径局限在当前单元内,并采用二分法确定子时间步数量,使算法简洁有效。算例表明,该方法在保证计算精度和求解稳定性的同时,可以明显减少计算时间,提高计算效率。

关键词: 热对流, 传热, Hele-Shaw, 注塑成型, 数值模拟, 聚合物加工

Abstract:

In Hele-Shaw flow simulation, the directly solved variable is pressure and the velocity is only the post-treating result of pressure. Around the injection gate, the velocity may be very high along with reducing elemental size. This means that when the energy conservation equation is solved as a whole, the time step must be very short, otherwise, the error in the heat convection is unavoidable. The above problem can be overcome by using the operator-splitting method, in which the material at the current computing points needs to be tracked back to its position in the last time-step. However, this may lead to a new difficulty. If the elemental velocity is very high, the tracking needs to pass through a few elements and the reversetracking may fail. To solve the problem, a new algorithm named sub time step with variable size was suggested to deal with the thermalconvection in this paper, in which the sub time step using dichotomy was introduced that bounded the tracking path in a certain element. The new algorithm made the computation more simple and effective. The numerical examples showed that the new method had same accuracy as one using uniform small time step and high solving stability, but calculating time was dramatically reduced.

Key words: thermal convection, heat transfer, Hele-Shaw, injection molding, numerical simulation, polymer processing

中图分类号: