Simulation of fiber-induced melt pressure fluctuations within large scale polymer composite deposition beads

Abstract

The process-structure-property relationship in Large Area Additive Manufacturing (LAAM) technology is an ongoing area of research as the inherent microstructural properties (chiefly fibers and voids) affect the performance of printed parts. Unfortunately, we currently lack adequate understanding of micro void nucleation and evolution during the LAAM and fused deposition modelling (FDM) processes. Modeling of the polymer melt flow during the extrusion process is important in understanding the underlying microstructural formation and associated properties of the print, that determines the part performance in service conditions. In this paper we compute fiber-induced local pressure fluctuations which may promote void formation in the bead’s microstructure. On a macro-scale, we determine flow fields of a purely viscous, Newtonian planar polymer deposition flow through a LAAM nozzle which are utilized on a micro-scale model where we simulate the evolution of a single ellipsoidal fiber along streamlines obtained from the macro-model. On the micro-scale, we determine instantaneous values of the translational and rotational velocities of the rigid ellipsoidal fiber that satisfies a balance of hydrodynamic forces and couples on the fiber’s surface based on a Newton Raphson algorithm and we track the fiber’s motion along the flow path via an explicit numerical iterative algorithm. Model verification is achieved by benchmarking results with solutions from well-known Jeffery’s equation of motion of a particle in homogeneous simple shear flow. We account for rotary diffusivity due to short-range fiber-fiber interaction in the FEA simulation by determining an effective fluid domain size representative of the interaction coefficient of the melt flow through a correlation analysis that yields an equivalent steady state orientation based on the Advani-Tucker equation. We also consider different possible motions of the fiber along individual LAAM flow paths from a given set of random initial fiber conditions to determine pressure bounds on the fiber surface along each streamline. For improved computational efficiency, calculations are carried out with respect to the fiber’s local coordinate axes to overcome the rigor of adaptive remeshing during the quasi-transient analysis. Results show low pressure extremes near the fiber’s surface which varies across the printed bead as well as through its thickness. Discussion is provided to gain insight into the effect of low-pressure extremes on micro void formation, particularly at the nozzle exit and during die swell/expansion.