Computational fluid dynamics simulation of the melting process in the fused filament fabrication additive manufacturing technique

被引:75
作者
Phan D.D. [1 ]
Horner J.S. [1 ]
Swain Z.R. [2 ]
Beris A.N. [1 ]
Mackay M.E. [1 ,2 ]
机构
[1] Department of Chemical and Biomolecular Engineering, University of Delaware, Newark, 19716, DE
[2] Department of Materials Science and Engineering, University of Delaware, Newark, 19716, DE
来源
Additive Manufacturing | 2020年 / 33卷
关键词
3D printing; Fused deposition modeling; Fused filament fabrication; Material extrusion; Numerical simulation;
D O I
10.1016/j.addma.2020.101161
中图分类号
学科分类号
摘要
Numerical simulation is used to understand the melting and pressurization mechanism in fused filament fabrication (FFF). The results show the incoming fiber melts axisymmetrically, forming a cone of unmelted material in the center surrounded by melted polymer. Details of the simulation reveal that a recirculating vortex of melted polymer is formed at the fiber entrance to the hot end. The large viscosity within this vortex acts to effectively seal the system against back-pressures of order 1000 psi (10 MPa), which are typical under standard printing conditions. The Generalized Newtonian Fluid (GNF) model was appropriate for simulation within the region that melts the fiber, however, a viscoelastic model, the Phan-Thien-Tanner (PTT) model, was required to capture flow within the nozzle. This is due to the presence of an elongational flow as molten material transitions from the melting region (diameter of 3 mm) to the nozzle at the exit (diameter of 0.5 mm). Remarkably, almost half the pressure drop occurs over the short capillary (0.5 mm in length) attached to the end of the converging flow region. Increased manufacturing rates are limited by high pressures, necessitating more consideration in the nozzle design of future FFF printers. © 2020 Elsevier B.V.
引用
收藏
相关论文
共 32 条
[21]  
Pyda M., Bopp R.C., Wunderlich B., Heat capacity of poly(lactic acid), J. Chem. Thermodyn., 36, pp. 731-742, (2004)
[22]  
Cox W.P., Merz E.H., Correlation of dynamic and steady flow viscosities, J. Polym. Sci., 28, pp. 619-622, (1958)
[23]  
Morrison F.A., Understanding Rheology, (2001)
[24]  
Boyd J., Buick J.M., Green S., Analysis of the Casson and Carreau-Yasuda non-Newtonian blood models in steady and oscillatory flows using the lattice Boltzmann method, Phys. Fluids, 19, (2007)
[25]  
Phan-Thien N., Tanner R.I., A new constitutive equation derived from network theory, J. Non-Newton. Fluid Mech., 2, pp. 353-365, (1977)
[26]  
Cogswell F.N., Measuring the extensional rheology of polymer melts, Trans. Soc. Rheol., 16, pp. 383-403, (1972)
[27]  
Mackay M.E., Astarita G., Analysis of entry flow to determine elongation flow properties revisited, J. Non-Newton. Fluid Mech., 70, pp. 219-235, (1997)
[28]  
Binding D.M., Walters K., On the use of flow through a contraction in estimating the extensional viscosity of mobile polymer solutions, J. Non-Newton. Fluid Mech., 30, pp. 233-250, (1988)
[29]  
Boger D.V., Binnington R.J., Circular entry flows of fluid M1, J. Non-Newton. Fluid Mech., 35, pp. 339-360, (1990)
[30]  
Debbaut B., Marchal J.M., Crochet M.J., Numerical simulation of highly viscoelastic flows through an abrupt contraction, J. Non-Newton. Fluid Mech., 29, pp. 119-146, (1988)