skip to main content


Title: Non-isothermal non-Newtonian three-dimensional flow simulation of fused filament fabrication
This work investigates three-dimensional simulation of fused filament fabrication using the Cross-WLF model for the non-isothermal and shear thinning behavior of the melt. To realistically simulate the deposition flow, the acceleration, viscosity evolution, and flow front tracking models have been included with the pressure gradient in the deposited road and boundary modeling of the melt and air interface. The results indicate that the non-isothermal and shear thinning behaviors greatly affect the geometry of the deposited roads including the flow front and trailing cross-section shapes. The thermal footprint of the interface between the deposited melt and the substrate is also predicted as a function of the thermal contact conductance. The pressure distribution within the deposited road is also modeled and is found to be not symmetric with respect to the nozzle center-line. Rather, the pressure peak shifts slightly downstream due to redirection of the melt around a stagnation point opposite the nozzle exit. Furthermore, a negative stress is observed downstream the exterior nozzle face associated with the free expansion of the melt as the extruded material climbs and releases from the exterior nozzle face. The developed simulation is verified by comparison with experimental results providing contact pressures ranging from 5 to 132 kPa.  more » « less
Award ID(s):
1914651
NSF-PAR ID:
10341311
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Additive manufacturing
Volume:
55
Issue:
102833
ISSN:
2214-7810
Page Range / eLocation ID:
1-15
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    The nozzle pressure was monitored in a fused filament fabrication process for the printing of high impact polystyrene. The contact pressure, defined as the pressure applied by the newly deposited layer onto the previous layer, is experimentally calculated as the difference between the pressure during printing and open discharge at the same volumetric flow rates. An analytical method for estimating the contact pressure, assuming one-dimensional steady isothermal flow, is derived for the Newtonian, power-law, and Cross model dependence of shear rates. A design of experiments was performed to characterize the contact pressure as a function of the road width, road height, and print speed. Statistical analysis of the results suggests that the contribution of the pressure driven flow is about twice that of the drag flow in determining contact pressure, which together describe about 60% of the variation in the observed contact pressure behavior. Modeling of the elastic and normal stresses at the nozzle orifice explains an additional 30% of the observed behavior, indicating that careful rheological modeling is required to successfully predict contact pressure. 
    more » « less
  2. Part design and process parameters directly influence the spatiotemporal distribution of temperature and associated heat transfer in parts made using additive manufacturing (AM) processes. The temporal evolution of temperature in AM parts is termed herein as thermal profile or thermal history. The thermal profile of the part, in turn, governs the formation of defects, such as porosity and shape distortion. Accordingly, the goal of this work is to understand the effect of the process parameters and the geometry on the thermal profile in AM parts. As a step towards this goal, the objectives of this work are two-fold: (1) to develop and apply a finite element-based framework that captures the transient thermal phenomena in the fused filament fabrication (FFF) additive manufacturing of acrylonitrile butadiene styrene (ABS) parts, and (2) validate the model-derived thermal profiles with experimental in-process measurements of the temperature trends obtained under different feed rate settings (viz., the translation velocity, also called scan speed or deposition speed, of the extruder on the FFF machine). In the specific context of FFF, this foray is the critical first-step towards understanding how and why the thermal profile directly affects the degree of bonding between adjacent roads (linear track of deposited material), which in turn determines the strength of the part, as well as, propensity to form defects, such as delamination. From the experimental validation perspective, we instrumented a Hyrel Hydra FFF machine with three non-contact infrared temperature sensors (thermocouples) located near the nozzle (extruder) of the machine. These sensors measure the surface temperature of a road as it is deposited. Test parts are printed under three different settings of feed rate, and subsequently, the temperature profiles acquired from the infrared thermocouples are juxtaposed against the model-derived temperature profiles. Comparison of the experimental and model-derived thermal profiles confirms a high-degree of correlation therein, with maximum absolute error less than 10%. This work thus presents one of the first efforts in validation of thermal profiles in FFF via in-process sensing. In our future work, we will focus on predicting defects, such as delamination and inter-road porosity based on the thermal profile.

     
    more » « less
  3. Abstract Part design and process parameters directly influence the instantaneous spatiotemporal distribution of temperature in parts made using additive manufacturing (AM) processes. The temporal evolution of temperature in AM parts is termed herein as the thermal profile or thermal history. The thermal profile of the part, in turn, governs the formation of defects, such as porosity and shape distortion. Accordingly, the goal of this work is to understand the effect of the process parameters and the geometry on the thermal profile in AM parts. As a step toward this goal, the objectives of this work are two-fold. First, to develop and apply a finite element-based framework that captures the transient thermal phenomena in the fused filament fabrication (FFF) additive manufacturing of acrylonitrile butadiene styrene (ABS) parts. Second, validate the model-derived thermal profiles with experimental in-process measurements of the temperature trends obtained under different material deposition speeds. In the specific context of FFF, this foray is the critical first-step toward understanding how and why the thermal profile directly affects the degree of bonding between adjacent roads (linear track of deposited material), which in turn determines the strength of the part, as well as, propensity to form defects, such as delamination. From the experimental validation perspective, we instrumented a Hyrel Hydra FFF machine with three non-contact infrared temperature sensors (thermocouples) located near the nozzle (extruder) of the machine. These sensors measure the surface temperature of a road as it is deposited. Test parts are printed under three different settings of feed rate, and subsequently, the temperature profiles acquired from the infrared thermocouples are juxtaposed against the model-derived temperature profiles. Comparison of the experimental and model-derived thermal profiles confirms a high degree of correlation therein, with a mean absolute percentage error less than 6% (root mean squared error <6 °C). This work thus presents one of the first efforts in validating thermal profiles in FFF via direct in situ measurement of the temperature. In our future work, we will focus on predicting defects, such as delamination and inter-road porosity based on the thermal profile. 
    more » « less
  4. Abstract The energy and momentum balance of an abyssal overflow across a major sill in the Samoan Passage is estimated from two highly resolved towed sections, set 16 months apart, and results from a two-dimensional numerical simulation. Driven by the density anomaly across the sill, the flow is relatively steady. The system gains energy from divergence of horizontal pressure work and flux of available potential energy . Approximately half of these gains are transferred into kinetic energy while the other half is lost to turbulent dissipation, bottom drag, and divergence in vertical pressure work. Small-scale internal waves emanating downstream of the sill within the overflow layer radiate upward but dissipate most of their energy within the dense overflow layer and at its upper interface. The strongly sheared and highly stratified upper interface acts as a critical layer inhibiting any appreciable upward radiation of energy via topographically generated lee waves. Form drag of , estimated from the pressure drop across the sill, is consistent with energy lost to dissipation and internal wave fluxes. The topographic drag removes momentum from the mean flow, slowing it down and feeding a countercurrent aloft. The processes discussed in this study combine to convert about one-third of the energy released from the cross-sill density difference into turbulent mixing within the overflow and at its upper interface. The observed and modeled vertical momentum flux divergence sustains gradients in shear and stratification, thereby maintaining an efficient route for abyssal water mass transformation downstream of this Samoan Passage sill. 
    more » « less
  5. ABSTRACT

    We analyse the internal structure and dynamics of cosmic-web filaments connecting massive high-z haloes. Our analysis is based on a high-resolution arepo cosmological simulation zooming-in on three Mpc-scale filaments feeding three massive haloes of $\sim 10^{12}\, \text{M}_\odot$ at z ∼ 4, embedded in a large-scale sheet. Each filament is surrounded by a cylindrical accretion shock of radius $r_{\rm shock} \sim 50 \, {\rm kpc}$. The post-shock gas is in virial equilibrium within the potential well set by an isothermal dark-matter filament. The filament line-mass is $\sim 9\times 10^8\, \text{M}_\odot \, {\rm kpc}^{-1}$, the gas fraction within rshock is the universal baryon fraction, and the virial temperature is ∼7 × 105 K. These all match expectations from analytical models for filament properties as a function of halo mass and redshift. The filament cross-section has three radial zones. In the outer ‘thermal’ (T) zone, $r \ge 0.65 \, r_{\rm shock}$, inward gravity, and ram-pressure forces are overbalanced by outward thermal pressure forces, decelerating the inflowing gas and expanding the shock outwards. In the intermediate ‘vortex’ (V) zone, 0.25 ≤ r/rshock ≤ 0.65, the velocity field is dominated by a quadrupolar vortex structure due to offset inflow along the sheet through the post-shock gas. The outward force is dominated by centrifugal forces associated with these vortices, with additional contributions from global rotation and thermal pressure. Shear and turbulent forces associated with the vortices act inwards. The inner ‘stream’ (S) zone, $r \lt 0.25 \, r_{\rm shock}$, is a dense isothermal core, $T\sim 3 \times 10^4 \, {\rm K}$ and $n_{\rm H}\sim 0.01 \, {\rm cm^{-3}}$, defining the cold streams that feed galaxies. The core is formed by an isobaric cooling flow and is associated with a decrease in outward forces, though exhibiting both inflows and outflows.

     
    more » « less