Search NASA⌕ Search

SEARCH · Search NASA

Results for “Liquid structure simulation”

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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Supramolecular meso -Trick: Ambidextrous Mirror Symmetry Breaking in a Liquid Crystalline Network with Tetragonal Symmetry

Bicontinuous and multicontinuous network phases are among nature's most complex structures in soft matter systems. Here, a chiral bicontinuous tetragonal phase is reported as a new stable liquid crystalline intermediate phase at the transition between two cubic phases, the achiral double gyroid and the chiral triple network cubic phase with an I 23 space group, both formed by dynamic networks of helices. The mirror symmetry of the double gyroid, representing a meso -structure of two enantiomorphic networks, is broken at the transition to this tetragonal phase by retaining uniform helicity only along one network while losing it along the other one. This leads to a conglomerate of enantiomorphic tetragonal space groups, P 4 1 2 1 2 and P 4 3 2 1 2. Phase structures and chirality were analyzed by small-angle X-ray scattering (SAXS), grazing-incidence small-angle X-ray scattering (GISAXS), resonant soft X-ray scattering (RSoXS) at the carbon K-edge, and model-dependent SAXS/RSoXS simulation. In conclusion, our findings not only lead to a new bicontinuous network-type three-dimensional mesophase but also reveal a mechanism of mirror symmetry breaking in soft matter by partial meso -structure racemization at the transition from enantiophilic to enantiophobic interhelical self-assembly.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Influence of Surfactants on Noble Metal Liquid-Gas Interface Mass Transfer in Molten Salt Reactors

Noble metals, a group of insoluble fission products (FPs), can accumulate in molten salt reactors (MSRs) and influence system performance and safety through deposition. Upon generation in fission, noble metals are transported by circulating salt and may either deposit on structural surfaces or enter the circulating gas bubble interfaces. While sparging was primarily designed for xenon and krypton removal, the mass transfer of noble metals into bubble interfaces highlights a potential pathway for their transport and partial extraction. The properties of the liquid-gas interface therefore play an important role in noble metal deposition and removal behaviors. Observations from earlier studies suggest that noble metals accumulated at the bubble interface can act as surfactants, making the bubble interface partially immobile. This behavior reduces the mass transfer capacity of circulating bubbles and affects both the overall distribution of noble metals in the loop and their removal efficiency. To investigate the influence of surfactants on noble metal transport, deposition, and extraction behaviors in MSRs, a previously developed and validated species transport model incorporating noble metal interphase mass transfer was used to simulate these behaviors. Two designed cases were studied to predict the distribution of 132Te along the Molten Salt Reactor Experiment (MSRE) loop under different bubble interface conditions. The results show that the mass transfer coefficients at the liquid-gas interface, determined by bubble interface characteristics, significantly influence noble metal distribution within the reactor loop. These findings emphasize the importance of continuously introducing fresh circulating bubbles into the molten salt to improve the removal efficiency of insoluble FPs from the MSR primary loop.

MSR↗
DOE OSTI2021

Imaging moiré flat bands in three-dimensional reconstructed WSe 2 /WS 2 superlattices

Moiré superlattices in transition metal dichalcogenide (TMD) heterostructures can host novel correlated quantum phenomena due to the interplay of narrow moiré flat bands and strong, long-range Coulomb interactions. However, microscopic knowledge of the atomically reconstructed moiré superlattice and resulting flat bands is still lacking, which is critical for fundamental understanding and control of the correlated moiré phenomena. Here we quantitatively study the moiré flat bands in three-dimensional (3D) reconstructed WSe 2 /WS 2 moiré superlattices by comparing scanning tunnelling spectroscopy (STS) of high-quality exfoliated TMD heterostructure devices with ab initio simulations of TMD moiré superlattices. A strong 3D buckling reconstruction accompanied by large in-plane strain redistribution is identified in our WSe 2 /WS 2 moiré heterostructures. STS imaging demonstrates that this results in a remarkably narrow and highly localized K-point moiré flat band at the valence band edge of the heterostructure. A series of moiré flat bands are observed at different energies that exhibit varying degrees of localization. Our observations contradict previous simplified theoretical models but agree quantitatively with ab initio simulations that fully capture the 3D structural reconstruction. Furthermore, our results reveal that the strain redistribution and 3D buckling in TMD heterostructures dominate the effective moiré potential and the corresponding moiré flat bands at the Brillouin zone K points.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Abiotic phosphorus recycling from adsorbed ribonucleotides on a ferrihydrite-type mineral: Probing solution and surface species

Iron (Fe) (oxyhydr)oxide minerals, which are amongst most reactive minerals in soils and sediments, are known to exhibit strong adsorption of inorganic phosphate (P i ) and organophosphate (P o ) compounds. Beyond synthetic P o compounds, much still remains unknown about the reactivity of these minerals to transform naturally-occurring P o compounds to P i , particularly with respect to solution versus surface speciation of P o hydrolysis. To investigate this reactivity with a ferrihydrite-type mineral and ribonucleotides, we employed high-resolution liquid chromatography-mass spectrometry (LC-MS), X-ray absorption near-edge structure (XANES), Fourier-transform infrared (FTIR) spectroscopy, and molecular modeling. Kinetic experiments were conducted with the mineral (1 g L -1 ) reacted with adenosine monophosphate, diphosphate, or triphosphate (respectively AMP, ADP, ATP; 50 µM). Analysis of solution organic species by LC-MS implied that only adsorption occurred with AMP and ADP but both adsorption and dephosphorylation of ATP were evident. Maximum adsorption capacities per gram of mineral were 40.6 ± 0.8 µmol AMP, 35.7 ± 1.6 µmol ADP, and 10.9 ± 1.0 µmol ATP; solution dephosphorylated by-products accounted for 15% of initial ATP. Subsequent XANES analysis of the surface species revealed that 16% of adsorbed AMP and 30% of adsorbed ATP were subjected to dephosphorylation, which was not fully quantifiable from the solution measurements. Molecular simulations predicted that ADP and ATP were complexed mainly via the phosphate groups whereas AMP binding also involved multiple hydrogen bonds with the adenosine moiety; our FTIR data confirmed these binding confirmations. Lastly, our findings thus imply that specific adsorption mechanisms dictate the recycling and subsequent trapping of Pi from ribonucleotide-like biomolecules reacted with Fe (oxyhydr)oxide minerals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Polarizable Water Potential Derived from a Model Electron Density

A new empirical potential for efficient, large scale molecular dynamics simulation of water is presented. The HIPPO (Hydrogen-like Intermolecular Polarizable POtential) force field is based upon the model electron density of a hydrogen-like atom. This framework is used to derive and parametrize individual terms describing charge penetration damped permanent electrostatics, damped polarization, charge transfer, anisotropic Pauli repulsion, and damped dispersion interactions. Initial parameter values were fit to Symmetry Adapted Perturbation Theory (SAPT) energy components for ten water dimer configurations, as well as the radial and angular dependence of the canonical dimer. The SAPT-based parameters were then systematically refined to extend the treatment to water bulk phases. The final HIPPO water model provides a balanced representation of a wide variety of properties of gas phase clusters, liquid water, and ice polymorphs, across a range of temperatures and pressures. This water potential yields a rationalization of water structure, dynamics, and thermodynamics explicitly correlated with an ab initio energy decomposition, while providing a level of accuracy comparable or superior to previous polarizable atomic multipole force fields. The HIPPO water model serves as a cornerstone around which similarly detailed physics-based models can be developed for additional molecular species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Genetic control of morphological transitions in a coacervating protein template

Nature routinely exploits liquid–liquid phase separation (LLPS) of proteins to control the assembly and mineralization of hybrid materials. Here, we show that fusion of the Car9 silica-binding peptide to an elastin-like polypeptide (ELP) yields temperature- and sequence-programmable soft matter templates for the synthesis of silicified architectures ranging in size from nanometers to micrometers. Specifically, we demonstrate unprecedented control over the diameter of silica nanoparticles (SiNP) in the 30–60 nm range with 4 nm precision, show that a single arginine residue (R4) in the Car9 sequence underpins the transition from micelles to proteinosomes, and find that substitutions in other basic residues modulate electrostatic repulsion and solvation to enable access to kinetically trapped species. These structures, which include interconnected micelles, small (∼200 nm) and large (>5 µm) vesicles, are readily visualized by SEM imaging following silicification. Molecular dynamics (MD) simulations and AlphaFold predictions reveal that mutations in positively charged residues alter interfacial packing, hydration, and conformational freedom of the silica-binding segments. Overall, our results establish sequence and thermal energy as synergistic levers for morphological control across length scales using solid-binding ELPs and establish mineralization as a powerful tool to visualize the structure of dynamic soft matter assemblies.

hierarchy↗
DOE OSTI2018

Decoding the spectroscopic features and time scales of aqueous proton defects

Acid solutions demonstrate a variety of complex structural and dynamical features arising from the presence of multiple interacting reactive proton defects and counterions. Yet, disentangling the transient structural motifs of proton defects in the water hydrogen bond network and the mechanisms for their interconversion remains a formidable challenge. In this work, we use simulations treating the quantum nature of both the electrons and nuclei to show how the experimentally observed spectroscopic features and relaxation time scales can be elucidated using a physically transparent coordinate that encodes the overall asymmetry of the solvation environment of the proton defect. We demonstrate that this coordinate can be used both to discriminate the extremities of the features observed in the linear vibrational spectrum and to explain the molecular motions that give rise to the interconversion time scales observed in recent nonlinear experiments. This analysis provides a unified condensed-phase picture of the proton structure and dynamics that, at its extrema, encompasses proton sharing and spectroscopic features resembling the limiting Eigen [H 3 O(H 2 O) 3 ]+ and Zundel [H(H 2 O) 2 ]+ gas-phase structures, while also describing the rich variety of interconverting environments in the liquid phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Liquid to crystal Si growth simulation using machine learning force field

Machine learning force field (ML-FF) has emerged as a potential promising approach to simulate various material phenomena for large systems with ab initio accuracy. However, most ML-FFs have been used to study the phenomena relatively close to the equilibrium ground states. In this work, we have studied a far from equilibrium system of liquid to crystal Si growth using ML-FF. Here, we found that our ML-FF based on ab initio decomposed atomic energy can reproduce all the aspects of ab initio simulated growth, from local energy fluctuations to transition temperatures, to diffusion constant, and growth rates. We have also compared the growth simulation with the Stillinger-Weber classical force field and found significant differences. A procedure is also provided to correct a systematic fitting bias in the ML-FF training process, which exists in all training models, otherwise critical results like transition temperature will be wrong.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Large-scale physically accurate modelling of real proton exchange membrane fuel cell with deep learning

Proton exchange membrane fuel cells, consuming hydrogen and oxygen to generate clean electricity and water, suffer acute liquid water challenges. Accurate liquid water modelling is inherently challenging due to the multi-phase, multi-component, reactive dynamics within multi-scale, multi-layered porous media. In addition, currently inadequate imaging and modelling capabilities are limiting simulations to small areas (<1 mm 2 ) or simplified architectures. Herein, an advancement in water modelling is achieved using X-ray micro-computed tomography, deep learned super-resolution, multi-label segmentation, and direct multi-phase simulation. The resulting image is the most resolved domain (16 mm 2 with 700 nm voxel resolution) and the largest direct multi-phase flow simulation of a fuel cell. This generalisable approach unveils multi-scale water clustering and transport mechanisms over large dry and flooded areas in the gas diffusion layer and flow fields, paving the way for next generation proton exchange membrane fuel cells with optimised structures and wettabilities.

25 ENERGY STORAGE↗
NASA NTRS1996

Organics In Meteorites

The variety of classes of organic compounds that occur in carbonaceous meteorites suggests a rich pre-planetary chemistry with possible connections to interstellar, solar nebular and parent body processes. Structural diversity prevails within all classes examined in detail. Among amino acids for instance, all possible isomers are found up to species containing 4-6 carbon atoms, with abundances decreasing with increasing molecular weight. Such diversity seems limited to those carbonaceous meteorites which show evidence of having been exposed to liquid water; meteorites lacking such evidence also show much lower abundances and less structural diversity in their organic contents. This apparent dependency on water suggests a role for cometary ices in the chemical evolution of organic compounds on parent bodies. Measurements of the stable isotope compositions of C, H, N and S in classes of compounds and at the individual compound level show strong deviations from average chondritic values. These deviations are difficult to explain by solar system or parent body processes, and precedents for some of these isotopic anomalies exist in interstellar (e.g., high D/H ratios) and circumstellar chemistry. Therefore, presolar origins for much if not all of the meteoritic organic compounds (or their precursors) is a distinct possibility. In contrast, evidence of solar nebular origins is either lacking or suspect. Results from molecular and isotopic analyses of meteoritic organics, from laboratory simulations and from a model of interstellar grain reactions will be used to flesh out the hypothesis that this material originated with interstellar chemistry, was distributed within the early solar system as cometary ices, and was subsequently altered on meteorite parent bodies to yield the observed compounds.

Chang, Sherwood↗
DOE OSTI2023

Localized High-Concentration Electrolytes Get More Localized

Liquid electrolytes in batteries are typically treated as macroscopically homogeneous ionic transport media, despite having complex chemical composition and atomistic solvation structures, thus leaving a knowledge gap of microstructural characteristics. A promising electrolyte system design has recently been developed with the title of a localized high-concentration electrolyte (LHCE). These electrolytes are a summation of three components – a charge-carrying salt, a solvent that solvates the salt, and a diluent that is sparingly interactive with the salt in terms of ion solvation. This creates the perception of a high-concentration electrolyte (HCE) with high anion-cation interactions that supports a salt-derived, stable solid-electrolyte interphase, while an LHCE improves bulk electrolyte properties, such as ionic conductivity and viscosity, when compared to an HCE. Much trial by error design of this LHCE concept has been supported by simulations to confirm ionic interactions, cycling data to confirm cell stability, and post-mortem analysis to provide evidence of mechanisms. However, a systematic design to extend the capabilities of LHCEs is still to be desired. Here, we analyzed a specific LHCE mixture to dictate the best viable mixture that extends the localization of solvation structures while maintaining viable cell-level capabilities. We seek to observe and understand solvent-to-salt interactions that impact electrolyte performance. These interactions range from single-ion solvation, multidentate interactions (i.e., solvation-driven clustering), to more widespread network-level interactions. With the support of both experimental and computational analysis of different electrolyte features, we have begun a methodology of formulating the most ideal LHCE mixtures for improving the cyclability of high-capacity battery chemistries.

25 ENERGY STORAGE↗
DOE OSTI2019

Patterning non-equilibrium morphologies in stimuli-responsive gels through topographical confinement

Stimuli-responsive “smart” polymers have generated significant interest for introducing dynamic control into the properties of antifouling coatings, smart membranes, switchable adhesives and cell manipulation substrates. Switchable surface morphologies formed by confining stimuli-responsive gels to topographically structured substrates have shown potential for a variety of interfacial applications. Beyond patterning the equilibrium swelling behavior of gels, subjecting stimuli-responsive gels to topographical confinement could also introduce spatial gradients in the various timescales associated with gel deformation, giving rise to novel non-equilibrium morphologies. In this work, we show how by curing poly(N-isopropylacrylamide) (pNIPAAm)-based gel under confinement to a rigid, bumpy substrate, we can not only induce the surface curvature to invert with temperature, but also program the transient, non-equilibrium morphologies that emerge during the inversion process through changing the heating path. Finite element simulations show that the emergence of these transient morphologies is correlated with confinement-induced gradients in polymer concentration and position-dependent hydrostatic pressure within the gel. To illustrate the relevance of such morphologies in interfacial applications, we show how they enable us to control the gravity-induced assembly of colloidal particles and microalgae. Finally, we show how more complex arrangements in particle assembly can be created through controlling the thickness of the temperature-responsive gel over the bumps. Patterning stimuli-responsive gels on topographically-structured surfaces not only enables switching between two invertible topographies, but could also create opportunities for stimuli ramp-dependent control over the local curvature of the surface and emergence of unique transient morphologies. Harnessing these features could have potential in the design of multifunctional, actuatable materials for switchable adhesion, antifouling, cell manipulation, and liquid and particle transport surfaces.

42 ENGINEERING↗
DOE OSTI2025

Medium-range order and compositional correlation in metallic glasses

The compositional atomic ordering in metallic glasses was studied by simulation focusing on the medium-range order (MRO). Many metallic alloy liquids and glasses show MRO characterized by the oscillations in the atomic pair-distribution function (PDF) beyond the first peak, which decay exponentially with distance. To study the effects of the local chemical order on MRO, we examine the compositionally resolved PDF and its MRO for models of various binary metallic alloy glasses. We show that compositional ordering is limited mostly to the nearest-neighbor atoms and the MRO is largely independent of the compositional order. For some elements that strongly repel each other in the alloy, a second MRO periodicity is observed owing to the distinct correlations among them. These results are discussed in light of the idea that the MRO oscillations in the PDF describe the correlations in the atomic density fluctuations, rather than the detailed local atomic structure.

Atomic structure↗
DOE OSTI2026

Prospects for detecting UF6 hydrolysis intermediates by mass spectrometry, resonance Raman, and NQR spectroscopy

Simulations are performed to consider several spectrometric and spectroscopic candidates for elucidating the mechanism of the hydrolysis of uranium hexafluoride (UF6). This study is among the first to benchmark the def-mTZVP basis sets for actinide-containing molecules, and it is shown to be a suitable basis set for surveying geometrical structures and vibrational spectra when used in conjunction with density functional theory. An experiment is proposed coupling mass spectrometric ion selection with vibrational spectroscopy, and supporting infrared and Raman spectral simulations demonstrate that ionization blue-shifts bands and can change their qualitative features. Ultraviolet resonance Raman is shown to have good prospects for discriminating U–O–U bridged intermediates, evidence for which was observed recently [L. E. McNamara et al., J. Phys. Chem. A 130, 775–786 (2026)]. As a more speculative approach, we also consider the prospect of addressing 233U or 235U nuclei by quadrupole resonance (NQR) spectroscopy for assigning early stage intermediate complexes. In doing so, we provide a first order-of-magnitude estimate for the collision-induced NQR signal for the UF6 dimer, which is of interest for the interpretation of a fast decoherence time observed in liquid-phase nuclear magnetic resonance. Having provided critical insights into the formation and stability of UF6 hydrolysis intermediates in previous studies, this computational spectroscopy survey is expected to help guide and expedite future laboratory campaigns.

Lutz, Jesse J. [Center for Computing Research, San↗
NASA NTRS1994

A double-moment multiple-phase four-class bulk ice scheme, Part 1: Description

A detailed ice-phase bulk microphysical scheme has been developed for simulating the hydrometeor distributions of convective and stratiform precipitation in different large-scale environmental conditions. The proposed scheme involves 90 distinct microphysical processes, which predict the mixing ratios and the number concentrations of small ice crystals, snow, graupel, and frozen drops/hail, as well as the mixing ratios of liquid water on wet precipitation ice (snow, graupel, frozen drops). The number of adjustable coefficients has been significantly reduced in comparison with other bulk schemes. Additional improvements have been made to the parameterization in the following areas: (1) representing small ice crystals with nonzero terminal fall velocities and dispersive size distributions, (2) accurate and computationally efficient calculations of precipitation collection processes, (3) reformulating the collection equation to prevent unrealistically large accretion rates, (4) more realistic conversion by riming between different classes of precipitation ice, (5) preventing unrealistically large rates of raindrop freezing and freezing of liquid water on ice, (6) detailed treatment of various rime-splintering ice multiplication mechanisms, (7) a simple representation of the Hobbs-Rangno ice enhancement process, (8) aggregation of small ice crystals and snow, and (9) allowing explicit competition between cloud water condensation and ice deposition rates rather than using saturation adjustment techniques. For the purposes of conserving the higher moments of the particle distributions, preserving the spectral widths (or slopes) of the particle spectra is shown to be more important than strict conservation of particle number concentration when parameterizing changes in ice-particle number concentrations due to melting, vapor transfer processes (sublimation of dry ice, evaporation from wet ice), and conversion between different hydrometeor species. The microphysical scheme is incorporated into a nonhydrostatic cloud model in Part 2 of this study. The model performed well in simulating the radar and microphysical structures of a midlatitude-continental squall lines and a tropical-maritime squall system with minimal tuning of the parameterization, even though the vertical profiles of radar reflectivity differed substantially between these storms.

Ferrier, Brad Shoenberg↗
NASA NTRS2008

25th Space Simulation Conference. Environmental Testing: The Earth-Space Connection

Topics covered include: Methods of Helium Injection and Removal for Heat Transfer Augmentation; The ESA Large Space Simulator Mechanical Ground Support Equipment for Spacecraft Testing; Temperature Stability and Control Requirements for Thermal Vacuum/Thermal Balance Testing of the Aquarius Radiometer; The Liquid Nitrogen System for Chamber A: A Change from Original Forced Flow Design to a Natural Flow (Thermo Siphon) System; Return to Mercury: A Comparison of Solar Simulation and Flight Data for the MESSENGER Spacecraft; Floating Pressure Conversion and Equipment Upgrades of Two 3.5kw, 20k, Helium Refrigerators; Affect of Air Leakage into a Thermal-Vacuum Chamber on Helium Refrigeration Heat Load; Special ISO Class 6 Cleanroom for the Lunar Reconnaissance Orbiter (LRO) Project; A State-of-the-Art Contamination Effects Research and Test Facility Martian Dust Simulator; Cleanroom Design Practices and Their Influence on Particle Counts; Extra Terrestrial Environmental Chamber Design; Contamination Sources Effects Analysis (CSEA) - A Tool to Balance Cost/Schedule While Managing Facility Availability; SES and Acoustics at GSFC; HST Super Lightweight Interchangeable Carrier (SLIC) Static Test; Virtual Shaker Testing: Simulation Technology Improves Vibration Test Performance; Estimating Shock Spectra: Extensions beyond GEVS; Structural Dynamic Analysis of a Spacecraft Multi-DOF Shaker Table; Direct Field Acoustic Testing; Manufacture of Cryoshroud Surfaces for Space Simulation Chambers; The New LOTIS Test Facility; Thermal Vacuum Control Systems Options for Test Facilities; Extremely High Vacuum Chamber for Low Outgassing Processing at NASA Goddard; Precision Cleaning - Path to Premier; The New Anechoic Shielded Chambers Designed for Space and Commercial Applications at LIT; Extraction of Thermal Performance Values from Samples in the Lunar Dust Adhesion Bell Jar; Thermal (Silicon Diode) Data Acquisition System; Aquarius's Instrument Science Data System (ISDS) Automated to Acquire, Process, Trend Data and Produce Radiometric System Assessment Reports; Exhaustive Thresholds and Resistance Checkpoints; Reconfigurable HIL Testing of Earth Satellites; FPGA Control System for the Automated Test of MicroShutters; Ongoing Capabilities and Developments of Re-Entry Plasma Ground Tests at EADS-ASTRIUM; Operationally Responsive Space Standard Bus Battery Thermal Balance Testing and Heat Dissipation Analysis; Galileo - The Serial-Production AIT Challenge; The Space Systems Environmental Test Facility Database (SSETFD), Website Development Status; Simulated Reentry Heating by Torching; Micro-Vibration Measurements on Thermally Loaded Multi-Layer Insulation Samples in Vacuum; High Temperature Life Testing of 80Ni-20Cr Wire in a Simulated Mars Atmosphere for the Sample Analysis at Mars (SAM) Instrument Suit Gas Processing System (GPS) Carbon Dioxide Scrubber; The Planning and Implementation of Test Facility Improvements; and Development of a Silicon Carbide Molecular Beam Nozzle for Simulation Planetary Flybys and Low-Earth Orbit.

Packard, Edward↗
DOE OSTI2021

Powder spreading, densification, and part deformation in binder jetting additive manufacturing

Binder jetting additive manufacturing (AM) can print complex structures in economical and scalable manner. Binder jetting AM comprises of deposition and weak binding of particles, known as green part, at room temperature and subsequent binder removal and sintering densification at high temperatures. However, during the densification (i.e., sintering), the part significantly deforms due to volume shrinkage. The deformation during sintering is difficult to predict, which prevents the widespread application of this technology. In this research, powder spreading simulation using discrete element method (DEM) was performed first to capture local variations in powder bed configuration. Second, finite element method (FEM) with a phenomenological continuum constitutive model was used to predict part shrinkage during the sintering process. DEM simulation showed variations in packing density, particle segregation, formation of uneven powder bed surface, and shift in particle size distribution (PSD). The sintering simulation modeled part deformation with a reasonable accuracy of 3% for solid-state sintering and intermediate liquid phase sintering. A demonstration case with non-uniform initial packing density showed that inhomogeneous green part density and PSD should be accounted for prediction of part deformation in binder jetting AM.

36 MATERIALS SCIENCE↗