Numerical simulation of reactive flow in liquid composite molding using flux-corrected transport (FCT) based finite element/control volume (FE/CV) method
Numerical simulation of reactive flow in liquid composite molding using flux-corrected transport (FCT) based finite element/control volume (FE/CV) method
复制标题
DOI：
10.1016/j.ijheatmasstransfer.2009.12.003
复制
发表时间：
2010-04
影响因子：
5.2
通讯作者：
H. Tan;K. Pillai
中科院分区：
文献类型：
--
作者：
H. Tan;K. Pillai
复合材料液态成型的充模过程模拟对于优化这种高性价比的复合材料制造工艺具有重要意义。LCM中的流动是对流占主导地位的反应，非等温流动，涉及移动边界，因此Galerkin有限元法（FEM）必须与稳定技术相适应，以解决此类问题。流线迎风Petrov-Galerkin（SUPG）方法是最常用的稳定化方法之一。然而，SUPG的使用仍然会导致局部数值摆动附近的尖锐的解决方案梯度，这是经常遇到的三维模具填充模拟LCM厚的零件。在这项研究中，我们提出使用通量校正输运（FCT）的有限元求解一组高度对流输运方程。本文提出的数值算例表明，FCT为基础的FEM在抑制寄生振荡的区域中的陡峭的解决方案梯度，而不是在这样的区域中的数值不稳定的SUPG的优异性能。首次将基于有限体积法的有限元法与控制体积法相结合，对LCM的非等温充模过程进行了数值模拟。我们已经开发了一个模拟程序PORE-FLOW©的基础上提出的方案的研究。数值研究证明了基于FCT的有限元法在模拟LCM充型过程中的稳定性。
The mold-filling simulation of liquid composite molding (LCM) is of great importance in optimizing this cost-effective polymer-composites manufacturing process. The flow in LCM is a convection-dominated reactive, non-isothermal flow involving a moving boundary, so the Galerkin finite element method (FEM) has to be adapted with stabilization techniques to solve such problems. The streamline-upwind Petrov–Galerkin (SUPG) method is one of the most popular stabilized methods. However, the use of SUPG still leads to localized numerical wiggles in the vicinity of sharp solution gradients, which is often encountered in the 3D mold-filling simulation of LCM for thick parts. In this study, we propose to use the flux-corrected transport (FCT) based FEM to solve a set of highly convective transport equations. The numerical examples presented in this paper demonstrate the excellent performance of FCT based FEM in suppressing spurious oscillations in the regions of steep solution-gradients as opposed to the numerical instability of SUPG in such regions. For the first time, the FCT based FEM combined with the control-volume method is employed to simulate the non-isothermal mold-filling process in LCM. We have developed a simulation code PORE-FLOW© based on the scheme proposed in the study. Numerical studies have proven the stability of the FCT based FEM while modeling the mold-filling process of LCM.