Structural integrity and quality of short fiber composite parts produced by Big Area Additive Manufacturing (BAAM) are largely affected by inherent bead microstructural features such as voids. Unfortunately, our understanding of void nucleation and evolution during polymer deposition process is lacking. Flow modeling focused on the associated microstructural formation provides a means for better understanding the process-structure-properties relations in large area extrusion deposition additive manufacturing of fiber reinforced composites. Our prior computational effort that investigated mechanisms that may promote micro-void formation was based on 2-dimensional planar models of a single ellipsoidal fiber motion in purely viscous polymer extrusion/deposition flow through a BAAM nozzle. Here we present a 3D finite element modelling approach to simulate single fiber out-of-plane rotations utilizing velocity and velocity gradient values computed along streamlines obtained from a 3D extrusion/deposition simulation of the BAAM polymer deposition process. The pressure distribution on the fiber’s surface along the flow path provides new insight into potential micro-void nucleation mechanism. Results show low pressure regions occur near the fiber’s surface which varies across the printed bead and through its thickness. 
                        more » 
                        « less   
                    This content will become publicly available on July 29, 2026
                            
                            Computational Fluid Dynamics Modeling of Material Deposition in Pneumatic Micro-Extrusion Additive Manufacturing
                        
                    
    
            Pneumatic micro-extrusion (PME) presents inherent challenges in achieving precise control over material deposition due to the small scale and high viscosity of the material flow. Traditional experimental methods often fall short in capturing the complex mechanisms, interactions, and dynamics governing material transport and deposition in PME, highlighting the necessity for advanced computational approaches to unveil these intricate physical phenomena. The primary objective of this study is to develop a robust computational framework for simulating material transport and deposition in PME, offering detailed insights into the fluid dynamics, flow complexities, and deposition characteristics intrinsic to the PME process. This involves systematically investigating the influence of key process parameters, such as material properties, print speed, and flow pressure, on deposition dynamics. Three-dimensional (3D) computational fluid dynamics (CFD) modeling was employed using ANSYS Fluent, with boundary conditions defined to replicate pneumatic extrusion on a moving substrate. The CFD simulations captured three distinct deposition regimes: (i) under-extrusion, (ii) normal extrusion, and (iii) over-extrusion. Results show that decreasing print speed at constant inlet velocity produced over-extrusion, with increased bead height and width due to material overflow. At normal extrusion, bead geometry remained stable, while under-extrusion reduced bead width and ultimately led to discontinuities when flow stresses exceeded cohesive strength. Notably, bead width decreased significantly between 10 mm/s and 15 mm/s, but showed little difference between 15 mm/s and 20 mm/s. Instead, the higher speed produced discontinuities consistent with experimental observations. Centerline velocity profiles revealed that the flow inside the cartridge was very slow at approximately 0.058 mm/s, whereas the velocity increased sharply within the nozzle throat, reaching nearly 12.88 mm/s. These predictions were further supported by validation experiments, where the measured nozzle velocity of 12.95 mm/s closely matched the CFD-simulated value of 12.88 mm/s, demonstrating strong agreement between simulation and experiment. Additionally, pressure decreased slightly with increasing print speed due to reduced backflow, while nozzle velocity and wall shear stress remained unchanged under fixed inlet velocity conditions. Overall, the outcomes of this study are expected to inform parameter optimization strategies and enhance PME process efficiency, providing critical insights into how variations in PME process dynamics influence print morphology and quality, advancing the path toward optimized fabrication of scaffolds for tissue engineering applications. 
        more » 
        « less   
        
    
                            - Award ID(s):
- 2327460
- PAR ID:
- 10638363
- Publisher / Repository:
- American Society of Mechanical Engineers (ASME) - International Mechanical Engineering Congress and Exposition (IMECE 2025)
- Date Published:
- Subject(s) / Keyword(s):
- Advanced Manufacturing Bone Tissue Engineering Regenerative Engineering Bone Pathology Natural Origin Materials.
- Format(s):
- Medium: X
- Location:
- Memphis, TN, USA
- Sponsoring Org:
- National Science Foundation
More Like this
- 
            
- 
            In this study we characterized the process–structure interactions in melt extrusion-based 3D bioplotting of polycaprolactone (PCL) and developed predictive models to enable the efficient design and processing of scaffolds for tissue engineering applications. First, the effects of pneumatic extrusion pressure (0.3, 0.4, 0.5, 0.6 N/mm2), nozzle speed (0.1, 0.4, 1.0, 1.4 mm/s), strand lay orientation (0°, 45°, 90°, 135°), and strand length (10, 20, 30 mm) on the strand width were investigated and a regression model was developed to map strand width to the two significant parameters (extrusion pressure and nozzle speed; p < 0.05). Then, proliferation of NIH/3T3 fibroblast cells in scaffolds with two different stand widths fabricated with different combinations of the two significant parameters was assessed over 7 days, which showed that the strand width had a significant effect on proliferation (p < 0.05). The effect of strand lay orientation (0° and 90°) on tensile properties of non-porous PCL specimens was determined and was found to be significantly higher for specimens with 0° lay orientation (p < 0.05). Finally, these data were used to develop and experimentally validate a finite element model for a porous PCL specimen with 1:1 ratio of inter-strand spacing to strand width.more » « less
- 
            Babski-Reeves, K; Eksioglu, B; Hampton, D. (Ed.)Extrusion-based 3D bioprinting is a promising method for repairing patient-specific tissues and organs due to its inherent capacity to release biocompatible materials containing living cells in a preset area. The filament geometry and width mostly determine the scaffold architecture. Extrusion pressure, print speed, print distance, nozzle diameter, and material viscosity are just a few of the process variables that can be carefully chosen to affect the filament shape and width, ultimately verifying the user-defined scaffold porosity. To maintain defined filament width variation for various hydrogels within an acceptable range and to confirm the overall geometric fidelity of the scaffold, in this paper, filament width for a set of biomaterial compositions was determined using an image processing technique and an analytical relationship, including various process parameters, was developed.more » « less
- 
            The formation and development of micro-voids within the bead microstructure of a polymer composite during the extrusion/deposition additive manufacturing process continues to be of interest given the adverse effect these features have on part quality. A computational method is employed here to investigate potential volatile-induced micro-void nucleation mechanism which simulates the evolution of a single rigid ellipsoidal fiber in purely viscous polymer extrusion/deposition flow through a Large Area Additive Manufacturing (LAAM) nozzle. Our previous studies on potential micro-void nucleation mechanisms have assumed a Newtonian fluid property definition for the polymer melt flow, the current study assesses the effect of assuming a generalized Newtonian fluid (GNF) model on the fiber’s response. Preliminary findings based on Jeffery’s flow assumption reveal the fiber’s orientation kinetics are unaffected by the shear thinning fluid behavior, however there is a reduction in the pressure distribution on the fiber’s surface as the power law index is decreased which is expected to reduce the likelihood for microvoid nucleation.more » « less
- 
            The addition of short carbon fibers to the feedstock of large-scale polymer extrusion/deposition additive manufacturing results in significant increases in mechanical properties dependent on the fiber distribution and orientation in the beads. In order to analyze those factors, a coupled computational fluid dynamics (CFD) and discrete element modeling (DEM) approach is developed to simulate the behavior of fibers in an extrusion/deposition nozzle flow after calibrations in simple shear flows. The DEM model uses bonded discrete particles to make up flexible and breakable fibers that are first calibrated to match Jeffery’s orbit and to produce interactions that are consistent with Advani-Tucker orientation tensor predictions. The DEM/CFD model is then used to simulate the processing of fiber suspensions in the variable flow and geometries present in extrusion/deposition nozzles. The computed results provide enhanced insight into the evolution of fiber orientation and distribution during extrusion/deposition as compared to existing models through individual fiber tracking over time and space on multiple parameters of interest such as orientation, flexure, and contact forces.more » « less
 An official website of the United States government
An official website of the United States government 
				
			 
					 
					
