Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective diffusivity coefficient”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

Study of quasi-collisional effects in laboratory and astrophysical plasmas

High-amplitude turbulence excited in plasmas at small-scales can leave imprint on the radiation produced by (accelerated) plasma particles -- mostly electrons in electron-ion plasmas and both electrons and positrons in lepton pair plasmas. Furthermore, turbulence is known to introduce "effective collisional" effects -- anomalous resistivity, dissipation, diffusion, etc. -- in the otherwise collisionless plasmas. These effective collisions affect radiation transfer and modify optical and magneto-optic effects -- the transmittance and reflectivity coefficients, Faraday rotation, etc. Note that the effective collisions considered in this project are not resonant wave-particle interactions, but much less conventional randomization of particles’ paths. Therefore, we colloquially refer them to as "quasi-collisions". The study of quasi-collisions on transport, radiative and optical properties of plasmas is the main focus of this project.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hydration and transport properties of cesium hydroxide and mixed cesium hydroxide–sodium nitrite aqueous solutions

Here, this study explores the hydration and transport properties of aqueous cesium hydroxide (CsOH) solution, with or without 1 molar (M) sodium nitrite (NaNO 2 ). Historic studies of electrolyte solutions indicate that Cs + ions decrease viscosity and increase diffusion rates, whereas OH − ions have the opposite effect. Here, the influence of OH − was dominant in CsOH solutions, leading to increased viscosity and reduced diffusion rates. There was a linear relationship between diffusion coefficients and water activity, emphasizing the significant role of ion–water interactions in determining transport properties. This may be because the interaction between Cs + and the anions is weak even when they are in direct contact with each other. The weakness of the ion-pairing was established through thermodynamic analysis. The findings suggest that ion-pairing is not the only important interaction controlling transport properties when ion-pairing is weak. Nonetheless, ion-pairing or obstructions did result in more sluggish transport properties as electrolyte concentrations increased. Overall, the research enhances the understanding of the complexities underlying ion interactions in multicomponent solutions.

Reynolds, Jacob G. [Hanford Site (HNF), Richland, ↗

Electrochemical measurement and analysis of YCl 3 , ScCl 3 , GdCl 3 and MgCl 2 in molten eutectic LiCl-KCl

For this work, electrochemical measurements of YCl 3 , ScCl 3 , GdCl 3 and MgCl 2 were obtained in molten LiCl-KCl eutectic at 773 K using a three-electrode configuration. Tungsten rods were used as the working and counter electrodes. Two quasi- and one reference electrode(s) (glassy carbon, nickel wire, and Ni/NiO) were used to collect electrochemical measurements of the analyte salts. Cyclic voltammetry data were analyzed to calculate values of diffusion coefficient, exchange current density, and charge transfer resistance and were determined to be on the order of 10 –5 cm 2 s –1 , 10 –2 A cm –2 , and 10 0 Ω; respectively. In general, the values of diffusion coefficients were found to be consistent with those reported of high-temperature molten salts. Relatively large values of exchange current density corresponded with smaller values of charge transfer resistance. These values were found to be reasonable in comparison to results available in the literature. The measured cyclic voltammograms were normalized with respect to both lithium reduction and chloride ion oxidation potentials. Such a normalization technique is effective for comparing experimentally obtained cyclic voltammetry data to those that have been published. An analysis of experimentally obtained results indicates the relation of electrode reactions (e.g., differences of reduction potentials) are independent of the choice of reference electrode. Additionally, the choice of reference electrode did not affect the electrochemical window, exchange current density, or diffusion coefficient values.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Physical and chemical effects of isopropanol exchange in cement-based materials

Drying cement-based materials is a necessary step to characterize the microstructure by microscopy, adsorption, or mercury intrusion. The isopropanol (IPA) exchange method is commonly used to replace the pore solution prior to drying, as it reduces the capillary pressure and thereby helps to preserve the microstructure of the hydrated cementitious materials. However, some physical and chemical effects of IPA on hydration products have been reported. These effects cannot be completely avoided, but can be reduced by shortening the exchange duration, in particular for effects on ettringite. This study carried out experiments with different exchange durations and IPA removal techniques to address research gaps in the literature. For 1-mm cubes of cement paste with IPA diffusion coefficient of 1 × 10{sup −11} m{sup 2}/s, the compared exchange durations are 5 min, 15 min (value from the literature), 40 min (when the mixture at the sample center achieves 95% of the drop in surface tension from water to IPA), 80 min (when the IPA/water mixture at the sample center reaches an azeotrope), 5 h (when over 99% IPA is replaced at the sample center), and 24 h (prolonged exchange). After IPA exchange, samples were either directly dried at 40 °C with flushing N{sub 2} or washed by diethyl ether prior to N{sub 2} drying. Even though the exchange duration of 40 and 80 min, in theory, can reduce the damage from capillary pressure, the nitrogen adsorption results do not show any advantage in preventing microstructural alteration. The 5-h exchange provided results similar to the 24-h exchange, but the latter induced slightly more serious chemical effects resulting from the longer contact of IPA with hydration products. In addition, the compared IPA removal techniques do not show significant differences. Therefore, this study suggests using small samples (~1 mm) for microstructural studies following IPA exchange for a period corresponding to the time when 99% of the water (mole fraction) at the sample center is replaced.

36 MATERIALS SCIENCE↗

High‐Throughput Evaluation of Hardening Coefficients of Eight Alloying Elements in Magnesium

Liquid–solid diffusion couples (LSDCs) are employed to generate a composition gradient in the single‐phase hexagonal closed‐packed (hcp) solid solution with compositions up to the solubility limit of various solutes in Mg. Nanoindentation scanning across the composition gradient in LSDCs allows effective evaluation of composition‐dependent hardness of eight alloying elements (Al, Ca, Ce, Gd, Li, Sn, Y, and Zn) in the hcp Mg phase. The hardening coefficients, an indicator of the potency of solid‐solution hardening, are evaluated from the measured composition‐hardness data and correlated with various materials properties such as atomic radius, shear modulus, and elastic modulus of the solutes. The rank of hardening potency of Al, Gd, Sn, Y, and Zn measured by nanoindentation is in good agreement with that measured by microindentation reported in the literature. The hardening coefficient (potency) from the strongest to the weakest is Ce > Ca > Y ≈ Gd > Zn > Al ≈ Sn > Li in Mg‐based hcp binary solid solutions. The hardening coefficient is found to be closely correlated with the strengthening potency.

36 MATERIALS SCIENCE↗

On the minimum transport required to passively suppress runaway electrons in SPARC disruptions

In Izzo et al (2022 Nucl. Fusion 62 096029), state-of-the-art modeling of thermal and current quench (CQ) magnetohydrodynamics (MHD) coupled with a self-consistent evolution of runaway electron (RE) generation and transport showed that a non-axisymmetric (n = 1) in-vessel coil could passively prevent RE beam formation during disruptions in SPARC, a compact high-field tokamak projected to achieve a fusion gain Q > 2 in DT plasmas. However, such suppression requires finite transport of REs within magnetic islands and re-healed flux surfaces; conservatively assuming zero transport in these regions leads to an upper bound of RE current ~1MA compared to ~8.7MA of pre-disruption plasma current. Further investigation finds that core-localized electrons, within r/a < 0.4 and with kinetic energies ~0.2–15 MeV, contribute most to the RE plateau formation. Yet only a relatively small amount of transport, i.e. a diffusion coefficient ~18 m 2 s -1 , is needed in the core to fully mitigate these REs. Properly accounting for (a) the CQ electric field's effect on RE transport in islands and (b) the contribution of significant RE currents to disruption MHD may help achieve this.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigating Morphology and Diffusion in Simulations of Precise Anion-Conducting Polymers

Using atomistic molecular dynamics simulations, we investigate the morphology and transport properties of a new class of polymers which are functionalized with quaternary ammonium groups for use as anion exchange membranes. The polymers are precision polyolefins with either a trimethylammonium (p5CNMe3) or a dimethyl-hexyl ammonium (p5CNMe2Hx) pendant group at every fifth carbon along a polyethylene backbone. Simulations are performed at hydration levels of 5, 10, 15, and 20 water molecules per ammonium group. The hydrated polymers form nanoscale, percolated hydrophilic domains (water channels) in the hydrophobic polymer matrix that become wider with increasing water content. Water and hydroxide anion diffusion coefficients also increase with increasing water content. The morphology of the water domains is similar in both polymers, while the diffusion coefficients are somewhat lower in p5CNMe2Hx at fixed water content. Furthermore, the diffusion coefficients in both polymers fall on the same curve as a function of the fractal dimension of the percolated water channels, which appears to be a useful scalar measure of the effects of the nanoscale morphology on water and hydroxide anion transport.

Anion exchange membrane↗

Xe Recovery from Nuclear Power Plants Off-Gas Streams: Molecular Simulations of Gas Permeation through DD3R Zeolite Membrane

Recent experimental work has shown zeolite membrane-based separation as a promising potential technology for Kr/Xe gas mixtures due to its much lower energy requirements in comparison to cryogenic distillation, the conventional separation method for such mixtures. Such a separation is also economically rewarding because Xe is in high demand, as a valuable product for many applications/processes. In this work, we have used Molecular Dynamics (MD) simulations to study the effects of different conditions, i.e., temperature, pressure, and gas feed composition, on Kr/Xe separation performance via DD3R zeolite membranes. We provide a comprehensive study of the permeation of the different gas species, density profiles, and diffusion coefficients. Molecular simulations show that if the feed is changed from pure Kr/Xe to an equimolar mixture, the Kr/Xe separation factor increases, which agrees with experiments. In addition, when Ar is introduced as a sweep gas, the adsorption of both Kr and Xe increases, while the permeation of pure Kr increases. A similar behavior is observed with equimolar mixtures of Kr/Xe with Ar as the sweep gas. High-separation Kr/Xe selectivity is observed at 50 atm and 425 K but with low total permeation rates. Changing pressure and temperature are found to have profound effects on optimizing the separation selectivity and the permeation throughput.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Convective and Turbulent Motions in Nonprecipitating Cu. Part III: Characteristics of Turbulence Motions

Velocity field in a nonprecipitating Cu under BOMEX conditions, simulated by SAM with 10-m resolution and spectral bin microphysics is separated into the convective part and the turbulent part, using a wavelet filtering. In Part II of the study properties of convective motions of this Cu were investigated. Here in Part III of the study, the parameters of cloud turbulence are calculated in the cloud updraft zone at different stages of cloud development. The main points of this study are (i) application of a fine-scale LES model of a single convective cloud allowed a direct estimation of turbulence parameters using the resolved flow in the cloud and (ii) the separation of the resolved flow into the turbulence flow and the nonturbulence flow allowed us to estimate different turbulent parameters with sufficient statistical accuracy. We calculated height and time dependences of the main turbulent parameters such as turbulence kinetic energy (TKE), spectra of TKE, dissipation rate, and the turbulent coefficient. It was found that the main source of turbulence in the cloud is buoyancy whose contribution is described by the buoyancy production term (BPT). The shear production term (SPT) increases with height and reaches its maximum near cloud top, and so does BPT. In agreement with the behavior of BPT and SPT, turbulence in the lower cloud part (below the inversion level) is weak and hardly affects the processes of mixing and entrainment. The fact that BPT is larger than SPT determines many properties of cloud turbulence. For instance, the turbulence is nonisotropic, so the vertical component of TKE is substantially larger than the horizontal components. Another consequence of the fact that BPT is larger than STP manifests itself in the finding that the turbulence spectrum largely obeys the -11/5 Bolgiano–Obukhov scaling. The classical Kolmogorov -5/3 scaling dominates for the low part of a cloud largely at the dissolving stage of cloud evolution. Using the spectra obtained we evaluated an “effective” dissipation rate which increases with height from nearly zero at cloud base up to 20 cm 2 s -3 near cloud top. The coefficient of turbulent diffusion was found to increase with height and ranged from 5 m 2 s -1 near cloud base to 25 m 2 s -1 near cloud top. In conclusion, the possible role of turbulence in the process of lateral entrainment and mixing is discussed.

54 ENVIRONMENTAL SCIENCES↗

Effect of CO 2 -brine-rock reactions on pore architecture and permeability in dolostone: Implications for CO 2 storage and EOR

Geologic carbon sequestration (GCS) is considered a feasible technology for storing substantive volumes of greenhouse gases in subsurface geological formations. In the reservoir, far from carbon dioxide (CO 2 ) injection wells or in post-injection scenarios, diffusion dominates over advection. This condition conjoins with spatially distributed geochemical reactions to induce heterogeneous changes in pore architecture, i.e. pore body and throat sizes or surface roughness. These changes can affect CO 2 transport properties and storage capacity. In this work, we investigated mineral dissolution and precipitation in dolomite samples saturated with a CO 2 -saturated brine at 93 °C and 34.5 MPa, aged without flow. Two rock types samples, i.e. intergranular- and vuggy-dominant, were selected to investigate changes in pore size, porosity and permeability under reactive conditions. Mineral dissolution and precipitation were characterized using scanning electron microscopy. Changes in pore size were quantified via time-domain nuclear magnetic resonance (TD-NMR) transverse relaxation time (T 2 ) and diffusion coefficient (D) distributions. We show that mineral dissolution likely occurs in highly permeable pathways. These observations are confirmed through analysis of (T 2 ) and diffusion coefficient (D) distributions. In contrast to results during CO 2 -enriched brine continuous injection, mineral precipitation was observed in micropores. The leftward shift of the T 2 peaks, corresponding to micropores, also evidenced mineral precipitation in lowpermeability zones. However, microscale alterations resulted only in a subtle increase in porosity and permeability. Results in this study shed light on effects of geochemical reactions on alteration of rock properties in diffusion-dominated regions during CO 2 storage.

58 GEOSCIENCES↗

HIGH-FIDELITY SIMULATION OF SOOT FORMATION AND THERMAL RADIATION IN A LABORATORY-SCALE RICH-QUENCH-LEAN BURNER

High-fidelity simulations of a swirl-stabilized turbulent spray flame in a laboratory-scale aero-combustor have been performed to evaluate the predictability of state-of-the-art models in capturing soot formation. The simulations employ a complex chemical mechanism developed for Jet-A with PAH chemistry, coupled with the Hybrid Method of Moments (HMOM) soot model, and a Lagrangian dilute spray model for the fuel injection. Two simulations are performed to compare the results when thermal radiation is neglected or included in the solution with a mean spectral model. Modeling closures for the soot differential diffusion effects in mixture fraction space, as well as turbulence-radiation interaction are also evaluated using the data generated by the simulations. Given the degree of complexity of the simulation, the results showed good agreement with experimental measurements of the spatial distribution of the soot volume fraction ensemble average. A closer agreement with the experiment is observed when thermal radiation is included in the solution. Thermal radiation is observed to reduce the flame temperature and increase the flame intermittency, denoted by the increase in the temperature standard deviation in mixture fraction space. The reduction in temperature also leads to a reduction in PAH production and soot volume fraction. Turbulence is observed to have different effects on radiative emission depending on the mixture fraction. Turbulent scalar fluctuations significantly enhance radiative emission in fuel lean mixtures and can also play a role for fuel rich conditions. The statistical description of the turbulence-radiation interaction, previously proposed in the literature, was observed to correctly reproduce the high-fidelity results. Model coefficients were provided for swirl-stabilized flames. The soot differential diffusion model, previously proposed in the literature, based on the residual between the exact term and its model approximation, was also evaluated. The residual correction term further improved the agreement with exact differential diffusion term evaluated with the high-fidelity simulation data in mixture fraction space. The results suggest that the effective turbulent Lewis number can be equal to unity in simulations of turbulent non-premixed recirculating flames.

Soriano, Bruno [Sandia National Laboratories (SNL)↗

Pore connectivity influences mass transport in natural rocks: Pore structure, gas diffusion and batch sorption studies

For this work, six rocks (one granodiorite, one limestone, two chalks, one mudstone, and one dolostone) with different extents of heterogeneity at six different particle sizes (from 75 to 8000 μm) were studied to describe the effects of pore connectivity on mass transport. The methods applied were (i) porosity measurement of granular rocks, (ii) analyses of gas-phase diffusive transport in a bed of packed particles, along with a solid quartz method at these six particle sizes being developed to identify the contribution of intraparticle diffusion, and (iii) batch sorption tests of multiple ions (anions and cations) with subsequent analyses of inductively coupled plasma-mass spectrometry. Granular porosity measurement results reveal that with decreasing particle sizes, the effective porosities for the “heterogenous” group of rocks (Grimsel granodiorite and Edwards limestone) increase, whereas the porosities of another “homogeneous” group (two Israel chalk samples, Japan mudstone, and Wyoming dolostone) remain constant. Gas diffusion results show that the intraparticle gas diffusion coefficient among these two sample groups, varying in the magnitude of 10 -8 to 10 -6 m 2 /s, are not directly correlated to the porosity differences. Moreover, the batch sorption work displays a different affinity of rocks for various tracers. For Grimsel granodiorite, Japan mudstone, and Wyoming dolostone, the adsorption capacity of Sm 3+ and Eu 3+ increases as the particle size decreases. In general, this integrated research of grain size distribution, granular rock porosity, intraparticle diffusivity, and ionic sorption capacity gives insights into the pore connectivity effect on both physical and chemical transport behaviors for different lithologies and/or different particle sizes.

58 GEOSCIENCES↗

Molecular simulation of lignin-related aromatic compound permeation through gram-negative bacterial outer membranes

Lignin, an abundant aromatic heteropolymer in secondary plant cell walls, is the single largest source of renewable aromatics in the biosphere. Leveraging this resource for renewable bioproducts through targeted microbial action depends on lignin fragment uptake by microbial hosts and subsequent enzymatic action to obtain the desired product. Recent computational work has emphasized that bacterial inner membranes are permeable to many aromatic compounds expected from lignin depolymerization processes. In this study, we expand on these findings through simulations for 42 lignin-related compounds across a gram-negative bacterial outer membrane model. Unbiased simulation trajectories indicate that spontaneous crossing for the full outer membrane is relatively rare at molecular simulation timescales, primarily due to preferential membrane partitioning and slow diffusion within the lipopolysaccharide layer within the outer membrane. Membrane partitioning and permeability coefficients were determined through replica exchange umbrella sampling simulations to overcome sampling limitations. We find that the glycosylated lipopolysaccharides found in the outer membrane increase the permeation barrier to many lignin-related compounds, particularly the most hydrophobic compounds. However, the effect is relatively modest; at industrially relevant concentrations, uncharged lignin-related compounds will readily diffuse across the outer membrane without the need for specific porins. Together, our results provide insight into the permeability of the bacterial outer membrane for assessing lignin fragment uptake and the future production of renewable bioproducts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling the Effects of Drift Orbit Bifurcation on Radiation Belt Electrons

Abstract Drift orbit bifurcation (DOB) has been suggested to play a major role in the loss and transport of radiation belt electrons since it violates the second adiabatic invariants of particles and makes the third invariant undefined. Results from our guiding‐center test particle simulations using the Tsyganenko‐1989c magnetic field model show that the DOB could affect a broad region of the outer radiation belt, which can penetrate inside the geosynchronous orbit at Kp ≥ 3, and its effects are more significant further away from the Earth, at a higher Kp and for higher electron energies. Specifically, the short‐term simulation results after one electron drift show both traditional and nontraditional DOB transport of electrons, with the nontraditional DOB, caused by a third minimum of the magnetic field strength near the equator, reported here for the first time. Moreover, our results show large ballistic jumps in the second invariant and radial distance for electrons at high equatorial pitch angles after one drift. In addition, long‐term DOB transport coefficients of electrons over many drifts are calculated based on our simulation results. We find that the pitch angle and radial diffusion coefficients of electrons due to DOB could be comparable to or even larger than those caused by electron interactions with chorus and Ultra‐Low‐Frequency waves, respectively. In sum, our results demonstrate that DOB could cause effective loss and transport of radiation belt electrons even in the absence of waves.

Astronomy & Astrophysics↗

Inner Belt Electron Decay Timescales: A Comparison of Van Allen Probes and DREAM3D Losses Following the June 2015 Storm

Abstract NASA's Van Allen Probes observed significant, long‐lived fluxes of inner belt electrons up to ∼1 MeV after geomagnetic storms in March and June 2015. Reanalysis of Magnetic Electron Ion Spectrometer (MagEIS) data with improved background correction showed a clearer picture of the relativistic electron population that persisted through 2016 and into 2017 above the Fennell et al. (2015, https://doi.org/10.1002/2014gl062874 ) limit. The intensity and duration of these enhancements allow estimation of decay timescales for comparison with simulated decay rates and theoretical lifetimes. We compare decay timescales from these data and DREAM3D simulations based on them using geomagnetic activity‐dependent pitch angle diffusion coefficients derived from plasmapause‐indexed wave data (Malaspina et al., 2016, https://doi.org/10.1002/2016gl069982 , 2018, https://doi.org/10.1029/2018gl078564 ) and phase space densities derived from MagEIS observations. Simulated decay rates match observed decay rates more closely than the theoretical lifetime due to significantly nonequilibrium pitch angle distributions in simulation and data. We conclude that nonequilibrium effects, rather than a missing diffusion or loss process, account for observed short decay rates.

79 ASTRONOMY AND ASTROPHYSICS↗

Phase-field model of freeze casting

Directional solidification of water-based solutions has emerged as a versatile technique for templating hierarchical porous materials. However, the underlying mechanisms of pattern formation remain incompletely understood. In this work, we present a detailed derivation and analysis of a quantitative phase-field model for simulating this nonequilibrium process. The phase-field model extends the thin-interface formulation of dilute binary alloy solidification with antitrapping to incorporate the highly anisotropic energetic and kinetic properties of the partially faceted ice-water interface. This interface is faceted in the basal plane normal to the ⟨0001⟩ directions and atomically rough in other directions within the basal plane. On the basal plane, the model reproduces a linear or nonlinear relationship between the interface growth rate and the kinetic undercooling that can be linked to experimental measurements. In both cases, spontaneous parity breaking of the solidification front is observed when the preferred growth direction is aligned with the temperature gradient. This phenomenon leads to the formation of partially faceted ice lamellae that drift laterally in one of the ⟨0001⟩ directions. Here, we demonstrate that the drifting velocity of the ice lamellae is controlled by the kinetics on the basal plane and converges as the thickness of the diffuse solid-liquid interface decreases. Furthermore, we examine the effect of the form of the kinetic anisotropy, which is chosen here such that the inverse of the kinetic coefficient varies linearly from a finite value in the ⟨0001⟩ directions to zero in all other directions within the basal plane, consistent with the assumption that the interface grows in local thermodynamic equilibrium in this plane. Our results indicate that the drifting velocity of ice lamellae is not affected by the slope of this linear relation, and the radius and undercooling at the tip of an ice lamella converge at relatively small slope values. Consequently, the phase-field simulations remain quantitative with computationally tractable choices of both the interface thickness and the slope assumed in the form of the kinetic anisotropy.

Materials science↗

Chemical Mixing Induced by Internal Gravity Waves in Intermediate-mass Stars

Internal gravity waves can cause mixing in the radiative interiors of stars. We study this mixing by introducing tracer particles into 2D hydrodynamic simulations. Following the work of Rogers & McElwaine, we extend our study to different masses (3, 7, and 20 M ⊙ ) and ages (ZAMS, midMS, and TAMS). The diffusion profiles of these models are influenced by various parameters such as the Brunt–Väisälä frequency, density, thermal damping, the geometric effect, and the frequencies of waves contributing to these mixing profiles. We find that the mixing profile changes dramatically across age. In younger stars, we noted that the diffusion coefficient increases toward the surface, whereas in older stars the initial increase in the diffusion profile is followed by a decreasing trend. We also find that mixing is stronger in more massive stars. Hence, future stellar evolution models should include this variation. In order to aid the inclusion of this mixing in 1D stellar evolution models, we determine the dominant waves contributing to these mixing profiles and present a prescription that can be included in 1D models.

79 ASTRONOMY AND ASTROPHYSICS↗

MODELLING FIBER ORIENTATION DURING ADDITIVE MANUFACTURING COMPRESSION MOLDING PROCESSES

The production of high-performance thermoplastic composites reinforced with short carbon fibers can be achieved by a novel “additive manufacturing-compression molding” technique. An advantage of such combination is two-fold: controlled fiber orientation in additive manufacturing and less void content by compression molding. In this study, a computational fluid dynamics model has been developed to predict the behavior of printed layers during fiber-reinforced thermoplastic extrusion and subsequent compression molding. The fiber orientation was modelled with a simple quadratic closure model. The interaction between the fibers is included using a rotary diffusion coefficient which becomes significant in concentrated regimes. Finally, the second order orientation tensor is coupled with the momentum equation as an anisotropic part of the stress term. The effect of processing parameters on the behavior of printed layers was investigated to determine the favorable printing scenarios. The developed numerical model enables design of high-performance composites with tunable mechanical properties.

Seta, Berin↗