SEARCH · Engineering Papers
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.
With the growing importance of developing new and improved methodologies for plastic recycling, conducting reproducible research and ensuring that results are transferable across labs are increasingly important. This Voices article reflects on how academia and industry view the path forward for strengthening reproducibility to advance science and enable a circular plastics economy.
Crystal plasticity (CP) is a powerful meso-scale technique for deformation modeling in polycrystalline materials. Crystal plasticity models typically do not include an explicit description of grain boundary (GB) sliding, which could lead to inaccurate predictions of strain distributions, especially in the vicinity of GBs. In the present study, a CP model is developed that includes GB sliding as an additional deformation mechanism and models the interaction between slip and GB sliding. Atomistic simulations are used to formulate the constitutive model for GB sliding and its interaction with incident slip. The deformation fields and structure of the GBs are obtained directly from experimental characterization and are faithfully reproduced as bicrystal systems for molecular dynamics simulations. The GBs are modeled as random-type (non-coincidence, asymmetrical) boundaries as observed from experimental data. The underlying atomistic-scale mechanism for pure sliding of random GBs was found to be analogous to fluid-flow. The interaction between incident slip and a sliding GB caused local increases in stress concentration, which further led to a momentary increase in the local sliding rate. The displacement profiles at the sliding GBs computed from the CP model with sliding accommodation are 38% more accurate than the baseline CP model as quantified by the mean squared error between the simulations and experiment. Here these results help improve the sophistication and accuracy of deformation modeling by including physics-based descriptions for GB sliding and its interaction with slip, eventually leading to more reliable predictions of micromechanical quantities.
Non-conservative processes play a fundamental role in plasticity and are behind important macroscopic phenomena such as creep, dynamic strain aging, loop raft formation, etc. In the most general case, vacancy-induced dislocation climb is the operating unit mechanism. While dislocation/vacancy interactions have been modeled in the literature using a variety of methods, the approaches developed rely on continuum descriptions of both the vacancy population and its fluxes. However, there are numerous situations in physics where point defect populations display heterogeneous concentrations and/or non-smooth kinetics. Here, a kinetic Monte Carlo (kMC) approach for modeling vacancy transport in response to arbitrary stress fields is used. Vacancies are treated as point particles and are coupled to the dislocation substructure representing a deformed material via an advection term defined by the local stress gradients. The stress fields and the dislocation substructure are evolved using a discrete dislocation dynamics (DDD) module. To extend the coupled model to the treatment of large systems, we have implemented it in the massively-parallel DDD code ParaDiS. To avoid numerical incompatibilities associated with merging deterministic (DDD) and stochastic (kMC) integration algorithms, we cast the entire elasto-plastic-diffusive problem within a single stochastic framework, taking advantage of a parallel kMC algorithm to evolve the system as a single event-driven process. The large-scale implementation enables the study of the evolution of a variety of dislocation-defect scenarios governed by non-conservative transport kinetics. After carrying out an exhaustive numerical and computational analysis of our parallel algorithm, we show results that emphasize situations where inhomogeneous vacancy dynamics are of relevance, and compare discrete kinetics to continuum solutions for several cases.
Lithium-ion batteries (LIBs) often suffer from capacity fading and poor cyclic performance due to mechanical degradation of the solid-electrolyte interphase (SEI). Here we perform numerical simulations and theoretical analysis to elucidate the role of plasticity in wrinkling and ratcheting behaviors of an SEI/electrode system. A coupled diffusion and finite deformation framework is formulated and numerically implemented as a user-element subroutine (UEL) to describe transient lithium diffusion and accompanying elastic–viscoplastic deformation of the electrode. It is found that concentration dependent properties and plastic deformation facilitate wrinkling in such a system. A wrinkled morphology may further lead to ratcheting and related failure under cycling. A phase diagram of four types of cyclic behaviorsis identified in terms of the charging rate and time. Our analysis suggests several potential strategies to avoid wrinkling and ratcheting instabilities, such as charging/discharging the electrode at a sufficiently slow rate, and/or introducing a thick artificial SEI with a pre-tension.
Classically, crystal plasticity modeling has used a range of constitutive equations, in which the incorporation of additional physics-based relationships typically results in additional model parameters. These additional parameters need to be reliably calibrated, which often necessitates the use of a range of experimental data acquired at multiple length scales. In this work, a crystal plasticity based finite element (CPFE) model for a dual-phase Titanium alloy, Ti–6Al–4V, is developed. The α and β phases of the microstructure are explicitly modeled. The model is calibrated using a systematic optimization routine and experimental data that consist of macroscopic stress-strain curves coupled with lattice strains on different crystallographic planes for the two phases. These experimental data were obtained from in situ high energy X-ray diffraction experiments for multiple material pedigrees, with varying crystallographic orientation distribution and β volume fractions. Depending on the thermomechanical-processing route and the heat treatment used to manufacture the alloy, Ti–6Al–4V can exist in a wide number of microstructural forms, which often results in the α and β phases either having well aligned slip systems (following the Burgers orientation relationship (BOR)) or possessing no alignment of the slip systems across the interphase boundary (not following the BOR). In this study, the fully-calibrated CPFE model is used to gain a comprehensive understanding of the deformation behavior of Ti–6Al–4V, specifically, the effect of microstructures that follow the BOR (or not) on time-dependent cyclic loading (including the effects of dwell hold times).
A new finite strain thermomechanical model for the high-rate deformation of the β-polymorph of cyclotetramethylene tetranitramine (β-HMX) has been developed and applied to simulations of plate impact experiments. The crystal plasticity model is based on a model developed previously for RDX (Luscher et al. 2017), which is extended to incorporate deformation twinning. Twinning during normal plate impacts is simulated with a phase-field twin model. First, material parameters governing the kinetics of dislocation slip are calibrated on the subset of simulations which had negative Schmid factors for the twin system. Second, a parametric study of the twin material parameters was performed to find suitable values. Overall, the results of the simulations with the phase-field twinning model are reported for impacts on several crystal orientations. We find that the twin growth decreases with increasing distance from the impact surface because of dissipation of the shock front via dislocation-mediated plasticity, and that the simulated interface velocity with and without phase-field twinning do not show appreciable differences. These modeling results suggest that the significance of twinning in β-HMX cannot be determined with traditional loading configurations and diagnostics.
Various methods for calibrating crystal plasticity finite element (CPFE) models lead to non-unique input parameter values, which subsequently introduce uncertainty in the predicted mechanical response. Sensitivity analysis (SA) conducted on crystal plasticity models is used to identify how variability in these parameters contribute to output uncertainty. Traditional SA on CPFE parameters uses simplified surrogate models to save computational time. However, the accuracy of the surrogate models depends on the quantity of training data used, and any modeling error can propagate into the SA results, potentially affecting their reliability. In this work, the elementary effects test (EET) method, a global SA technique using direct CPFE simulations was employed, and the results obtained were compared with the First Order Second Moment (FOSM) method. ExaConstit, an open-source GPU-enabled CPFE code, was used to perform the simulations and direct SA. The EET method was accurately able to capture the non-linear effects of all the input parameters on the output and is a valuable approach for reliably attributing parameter sensitivities in CPFE models. Based on the results, efficient strategies to perform future parameter calibration and SA are discussed. Additionally, the SA trends observed in different single crystal orientations closely mirrored the activity of the slip systems.
In this work, the influence of stress relief on the plasticity and fracture behavior of Inconel 625 fabricated through laser powder bed fusion additive manufacturing (AM) was investigated. The as-built versus stress relieved microstructures were compared, showing similar grain structures but the presence of ~10 vol % δ phase in the stress relieved condition, and no δ phase in the as-built condition. Mechanical tests under plane strain tension were performed on the stress relieved samples, and an anisotropic plasticity model was calibrated and validated using finite element simulations. Uniaxial and notched tension tests were performed on both as-built and stress relieved samples to probe the effect of stress relief on stress state- and direction-dependent fracture behavior. It was found that on average, the fracture strain of the stress relieved samples along the build direction was 30% higher than that along the perpendicular build direction in the stress state range studied, and the stress relief heat treatment resulted in a 45% decrease in fracture strain. The fracture strain in stress relieved samples was more strongly dependent on stress state than in as-built samples.
We studied the Vickers microhardness HV0.025 of CdZnTe and CdZnTeSe samples upon illumination with quasi-monochromatic light with wavelengths ranging between 750 and 1540 nm. We observed that light with a wavelength close to or above the bandgap of the samples results in a positive photo-plastic effect (hardening), and light in the spectral range well below the band gap can lead both to positive or negative effects depending on the parameters of deep levels in the bandgap. We performed measurements of spectrally dependent photoconductivity and compared them with HV0.025. Based on the analysis, we propose a complex qualitative model explaining the spectral dependence of HV0.025 by the charging and neutralization of deep levels by optical transitions between the valence and conduction bands and between the bands and deep levels in the band gap. Tuning of parameters of deep levels could be used to control the sign and magnitude of the photo-plastic effect.
Passive daytime radiative cooling, requiring zero external energy consumption, is a promising cooling strategy achieved by simultaneously reflecting solar irradiance and thermally radiating heat into the cold outer space (~3 K) through the atmospheric transparency window. However, current materials for passive radiative cooling face huge challenges, such as complicated fabrication approaches, expensive raw materials, and environmental requirements for practical applications. In line with the urgent need for plastic recycling to curb global environmental pollution, the recycled plastics are used to fabricate a passive radiative cooling material. Herein, the foam-paper composite (FPC) with excellent self-cooling capability is fabricated by a simple crushing-and-mixing procedure using recycled polystyrene (PS) foam and printer paper. The superhydrophobic PS foam particles not only protect the FPC from water damage for field applications but also reinforce its solar reflectivity via their porous structure. The cellulose fibers in printer paper can efficiently emit infrared thermal radiation into the cold outer space and bond dispersed PS foam particles together, further increasing its mechanical strength. The combination of highly diffusely reflective PS foam particles and fiber-based printer paper results in a reflectivity of 96% in the solar spectrum, a sub-ambient cooling performance of 8.4 degrees C, and a maximum radiative cooling power of 90 W/m2 during a 24-h cycle. Meanwhile, the FPC with high humidity can maintain its high solar reflectivity, which promotes its application in humid subtropical areas. Further, the low material cost and ease of fabrication will provide a path for effective daytime radiative cooling, especially in less developed areas.
The scintillation process in plastic scintillators has been studied with mixed systems containing high concentrations of multiple fluorescence dyes. It has been shown that the triplet–triplet interaction phenomena leading to the production of delayed light and pulse shape discrimination (PSD) can vary depending on the triplet energies of the interacting dyes. At small differences between the first excited triplet state energies, found to be below 0.1 eV, a pair of interacting dyes of different molecular species can be involved in heterogeneous excitation migration and triplet–triplet annihilation that result in the enhancement of delayed light and PSD, due to the engagement of both dyes in the process. At differences above 0.27 eV, excitation trapping on the lower energy triplet molecules leads to conditions when triplet–triplet migration and annihilation can proceed only homogeneously between molecules of the same species. The results explain the different effects produced by the selection of fluorescent dyes in previously reported PSD systems containing one primary and one secondary dyes. They also demonstrate a non-traditional approach to the design of new plastic scintillators with multiple dyes that, at improved scintillation and PSD performance, may simplify the production process.
Microwave (MW)-assisted catalytic pyrolysis offers a promising pathway for efficient plastic upcycling. This work develops an integrated modeling framework combining dynamic data reconciliation, a temperature-dependent rate model, and a yield model to represent the time-varying production rate of components in MW-assisted LDPE pyrolysis conducted in a batch reactor. An Arrhenius-type rate model with a temperature-dependent reaction order is developed. A biexponential correlation is proposed for the yield of gaseous products that enables to capture the evolving product formation behavior during conversion. In the yield correlation, one term is used to represent the initial increase in yield, reflecting the rapid formation of intermediate or primary products at the early stages of the reaction when a larger fraction of the reactant remains available. As conversion progresses, the influence of this term gradually diminishes. The other term accounts for the subsequent decrease in the predicted yield, representing secondary reactions such as further cracking or coke formation that reduce the concentration of certain products at higher conversion. The model is found to accurately represent reconciled experimental flow rate profiles from an in-house MW-assisted catalytic batch reactor for major products, including ethylene, ethane, 1-butene, and benzene, across 250−350 °C. Ethylene remains the dominant product but decreases from about 41.95% at 250 °C to 30.14% at 350 °C, while heavier products increase significantly, with 1-butene rising to nearly 8.37% and benzene reaching 2.17% at intermediate temperatures. The model shows that the ethylene production rate can be maximized at around 270 °C. The models developed in this work can be utilized for process optimization, reactor design and scale-up of microwave-assisted plastic conversion technologies, and economic analysis.
Many thermoplastic polymers are ductile by combining strength and large deformations. These deformations are irreversible - known as plastic deformation. Elastomers can deform reversibly but have low strength. To this end, we developed glassy and ductile polyamide networks capable of large plastic deformation (>200% strain) and high strength (~50 MPa tensile strength and ~1500 MPa Young’s modulus), similar to those of polyolefins and Nylon-66. We discovered that hydrogen bonding between meta-phthalamide groups was essential to the ductility. Since these polyamide networks are covalently bonded, we demonstrated their unique durability by repeatable elastic recovery at elevated temperatures, exhibiting indifferent tensile properties in each cycle. Furthermore, when we fixed the strain during the elastic recovery, these polyamide networks actuated stresses of 9–18 MPa, among the highest reported in shape-memory polymer actuators. Here, we envision these ductile, glassy polymer networks as promising alternatives to ductile thermoplastics, given the combined benefit of ductility and durability.
Earth System Models generally predict increasing upper ocean stratification from 21st century warming, which will cause a decrease in the vertical nutrient flux forcing declines in marine net primary productivity (NPP) and carbon export. Recent advances in quantifying marine ecosystem carbon to nutrient stoichiometry have identified large latitudinal and biome variability, with low-latitude oligotrophic systems harboring pico-sized phytoplankton exhibiting large phosphorus to carbon cellular plasticity. The climate forced changes in nutrient flux stoichiometry and phytoplankton community composition are thus likely to alter the ocean's biogeochemical response and feedback with the carbon-climate system. We have added three pico-phytoplankton functional types within the Biogeochemical Elemental Cycling component of the Community Earth System Model while incorporating variable cellular phosphorus to carbon stoichiometry for all represented phytoplankton types. The model simulates Prochlorococcus and Synechococcus populations that dominate the productivity and sinking carbon export of the tropical and subtropical ocean, and pico-eukaryote populations that contribute significantly to productivity and export within the subtropical to mid-latitude transition zone, with the western subtropical regions of each basin supporting the most P-poor stoichiometries. Subtropical gyre recirculation regions along the poleward flanks of surface western boundary currents are identified as regional hotspots of enhanced carbon export exhibiting C-rich/P-poor stoichiometry, preferentially inhabited by pico-eukaryotes and diatoms. Collectively, pico-phytoplankton contribute ∼58% of global NPP and ∼46% of global particulate organic carbon export below 100 m through direct and ecosystem processing pathways. Biodiversity and cellular nutrient plasticity in marine pico-phytoplankton combine to increase their contributions to ocean productivity and the biological carbon pump.
Abstract The chemical deconstruction of polyolefins to fuels, lubricants, and waxes offers a promising strategy for mitigating their accumulation in landfills and the environment. Yet, achieving true recyclability of polyolefins into C 2 -C 4 monomers with high yields, low energy demand, and low carbon dioxide emissions under realistic polymer-to-catalyst ratios remains elusive. Here, we demonstrate a single-step electrified approach utilizing Rapid Joule Heating over an H-ZSM-5 catalyst to efficiently deconstruct polyolefin plastic waste into light olefins (C 2 -C 4 ) in milliseconds, with high productivity at much higher polymer-to-catalyst ratio than prior work. The catalyst is essential in producing a narrow distribution of light olefins. Pulsed operation and steam co-feeding enable highly selective deconstruction (product fraction of >90% towards C 2 -C 4 hydrocarbons) with minimal catalyst deactivation compared to Continuous Joule Heating. This laboratory-scale approach demonstrates effective deconstruction of real-life waste materials, resilience to additives and impurities, and versatility for circular polyolefin plastic waste management.
Cultivation of microorganisms in space has enormous potential to enable in-situ resource utilization (ISRU) Here, we develop an autonomous payload with fully programmable serial passaging and sample preservation, termed the Modular Open Biological Platform (MOBP), and flight-test the MOBP aboard the International Space Station (ISS) by conducting enzymatic and microbial plastics upcycling experiments. The MOBP is a compact, modular bioreactor system that allows for sustained microbial growth via automated media transfers, such as those for sample collection and storage for terrestrial analyses, and precise data monitoring from integrated sensors. The MOBP was flight-tested with two experiments designed to evaluate biological upcycling of the plastic poly(ethylene terephthalate) (PET). The bioproduct βKA can be polymerized into a nylon-6,6 analog with improved properties for use in the production of a variety of materials. We posit the MOBP will aid in democratizing the execution of synthetic biology in spaceflight towards enabling ISRU.
The corrosion dealloying behavior of cold-worked (CW) Ni20Cr alloy (wt%) was studied in molten LiF-NaF-KF (or FLiNaK) salts at 600 °C, equal to a homologous temperature (TH) of 0.52. Alloys were cold-rolled to achieve reductions of thickness of 10%, 30%, and 50% introducing plastic deformation and a high density of dislocations. Potentiostatic holds (Eapplied) were applied in two different electrode potential regimes. At 1.75VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$, Cr dealloying to Cr(II) and Cr(III) is predominant, while at 1.90 VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$, both Ni and Cr are oxidized in molten FLiNaK at 600 °C. In these potential regimes, dealloyed Ni20Cr displayed bicontinuous porosity within the grain interior and at grain boundaries, driven by the high driving force for Cr dissolution and sustained by defect mediated outward solid state diffusion of Cr in parallel with surface diffusion of Ni. The bicontinuous porous structure developed was observed to undergo further coarsening and densification of the Ni-rich ligaments at a higher electrode potential. The main effect of CW observed is the introduction of plastic deformation and dislocation substructures that serve as short-circuit paths for Cr solid state diffusion to surfaces exposed to FLiNaK. This modified the evolution of the bicontinuous porous structure which increased with CW substantially. Kinetic analysis reveals that the Cr dealloying at +1.75 VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$ and 1.90 VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$ is initially charge transfer controlled, except for the 50% CW condition at 1.90 VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$, where the process becomes limited by slow Cr defect mediated bulk diffusion. The rate determining factors are explored and compared to experimental results.
Abstract While there exists in nature abundant examples of materials with site-specific gradients in microstructures and properties, engineers and designers have traditionally used monolithic materials with discrete properties. Now, however, additive manufacturing (AM) offers the possibility of creating structures that mimic some aspects of nature. One example that has attracted attention in the recent years is the hierarchical structure in bamboo. The hierarchical architecture in bamboo is characterized by spatial gradients in properties and microstructures and is well suited to accommodate and survive complex stress states, severe mechanical forces, and large deformations. While AM has been used routinely to fabricate functionally graded materials, this study distinguishes itself by leveraging AM and physical metallurgy concepts to trigger cascading deformation in a single sample. Specifically, we have been successful in using AM to fabricate steel with unique spatial hierarchies in structure and property to emulate the structure and deformation mechanisms in natural materials. This study shows an improvement in the strength and ductility of the nature-inspired “hierarchical steel” compared with conventional cast stainless steels. In situ characterization proves that this improvement is due to the sequential activation of multiple deformation mechanisms namely twinning, transformation-induced plasticity, and dislocation-based plasticity. While significantly higher strengths can be achieved by refining the chemical and processing technique, this study sets the stage to achieve the paradigm of using AM to fabricate structures which emulate the flexibility in mechanical properties of natural materials and are able to adapt to in-service conditions.