Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “experimental evolution”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

GSAS Tools

SAND2023-06684O GSAS Tools is a web application that manages user access to modeling and simulation tools and promotional material. This software, which is a spiking neural network (SNN) simulator, represents an SNN as a graph of stochastic differential equations and simulates the time-evolution of these equations. It has the capability of reading inputs from file and writing outputs to file, and generally supports experimentation, analysis, and algorithm development using SNNs. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Noel, Todd↗
DOE OSTI2020

Intense Pulse Light Annealing of Perovskite Photovoltaics Using Gradient Flashes

Perovskite solar cells (PSCs) have been fabricated through high-speed and low-cost depositions but often require long annealing. Intense pulse light (IPL) can anneal thin films in seconds after deposition by inducing very high temperatures lasting for milliseconds, and multiple flashes can be used to tune the temperature profile. In this study, a gradient flash annealing (GFA) approach is introduced and compared to uniform flash annealing (UFA) by investigating the crystallinity, morphology, and phase evolution of the CH 3 NH 3 PbI 3 perovskite films and their impact on PSC performance. Unlike UFA, low-intensity pulsed irradiation during the pre-annealing stage of GFA played a significant role on enhancing the PSC performance by forming pure-phase CH 3 NH 3 PbI 3 perovskite thin films with superior morphology and high crystallinity. To understand the kinetics, a transient thermal simulation using ANSYS was developed and confirmed with an experimental setup. The results combined with microscopy and spectrophotometry were used to visualize how duration, delay time, photon flux, and flash count parameters participated in crystallization, phase, and morphology evolution. The IPL annealing induced rapid surface temperature increase, reaching as high as 800 °C, while produced well-developed and bound perovskite grains without any surface defects in an uncontrolled ambient environment with high humidity (>60%). The study resulted in PSCs with maximum efficiency and fill factor of 9.27 and 69.92% for the UFA and 11.75 and 68% when annealed through the GFA approach, respectively. Lastly, this work utilized IPL as the sole thermal source for post-deposition annealing to rapidly fabricate efficient PSCs in only 10 s, which opens the pathway for high-speed, low-cost, and large-scale automated fabrication of perovskite photovoltaics and semiconductors.

14 SOLAR ENERGY↗
DOE OSTI2022

Sample Preparation and Experimental Design for In Situ Multi-Beam Transmission Electron Microscopy Irradiation Experiments

We discuss the need to understand materials exposed to overlapping extreme environments such as high temperature, radiation, or mechanical stress. When these stressors are combined there may be synergistic effects that enable unique microstructural evolution mechanisms to activate. Understanding of these mechanisms is necessary for the input and refinement of predictive models and critical for engineering of next generation materials. The basic physics and underlying mechanisms require advanced tools to be investigated. The in situ ion irradiation transmission electron microscope (I 3 TEM) is designed to explore these principles. To quantitatively probe the complex dynamic interactions in materials, careful preparation of samples and consideration of experimental design is required. Particular handling or preparation of samples can easily introduce damage or features that obfuscate the measurements. There is no one correct way to prepare a sample; however, many mistakes can be made. The most common errors and things to consider are highlighted within. The I 3 TEM has many adjustable variables and a large potential experimental space, therefore it is best to design experiments with a specific scientific question or questions in mind. Experiments have been performed on large number of sample geometries, material classes, and with many irradiation conditions. The following are a subset of examples that demonstrate unique in situ capabilities utilizing the I 3 TEM. Au nanoparticles prepared by drop casting have been used to investigate the effects of single ion strikes. Au thin films have been used in studies on the effects of multibeam irradiation on microstructure evolution. Zr films have been exposed to irradiation and mechanical tension to examine creep. Ag nanopillars were subjected to simultaneous high temperature, mechanical compression, and ion irradiation to study irradiation induced creep as well. These results impact fields including: structural materials, nuclear energy, energy storage, catalysis, and microelectronics in space environments.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Engineering grain boundary anisotropy to elucidate grain growth behavior in alumina

Current grain growth models have evolved to account for the relationship between grain boundary energy/mobility anisotropy and the five degrees of grain boundary character. However, the role of grain boundary networks on overall growth kinetics remains poorly understood. To experimentally investigate this problem, a highly textured Al 2 O 3 was fabricated by colloidal casting in a strong magnetic field to engineer a unique spatial distribution of grain boundary character. Microstructural evolution was quantified and compared to an untextured sample. From this comparison, a prevalence of (0001)/(0001) terminated grain boundaries with anisotropic networks were identified in the textured sample. These boundaries and their networks were found to be driving grain growth at a faster rate than predicted by models. These findings will allow better modelling of grain growth in real systems by experimentally exploring the impact thereon of grain boundary plane anisotropy and relative energy/mobility differences between neighboring boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Time correlations from steady-state expectation values

Recovering properties of correlation functions is typically challenging. On one hand, experimentally, it requires measurements with a temporal resolution finer than the system's dynamics. On the other hand, analytical or numerical analysis requires solving the system evolution. Here, we use recent results of quantum metrology with continuous measurements to derive general lower bounds on the relaxation and second-order correlation times that are both easy to calculate and measure. These bounds are based solely on steady-state expectation values and their derivatives with respect to a control parameter, and can be readily extended to the autocorrelation of arbitrary observables. We validate our method on two examples of critical quantum systems: a critical driven-dissipative resonator, where the bound matches analytical results for the dynamics, and the infinite-range Ising model, where only the steady state is solvable and thus the bound provides information beyond the reach of existing analytical approaches. Our results can be applied to experimentally characterize ultrafast systems, and to theoretically analyze many-body models with dynamics that are analytically or numerically hard.

Górecki, Wojciech [INFN, Pavia]↗
DOE OSTI2025

Time correlations from steady-state expectation values

Recovering properties of correlation functions is typically challenging. On one hand, experimentally, it requires measurements with a temporal resolution finer than the system's dynamics. On the other hand, analytical or numerical analysis requires solving the system evolution. Here, we use recent results of quantum metrology with continuous measurements to derive general lower bounds on the relaxation and second-order correlation times that are both easy to calculate and measure. These bounds are based solely on steady-state expectation values and their derivatives with respect to a control parameter, and can be readily extended to the autocorrelation of arbitrary observables. We validate our method on two examples of critical quantum systems: a critical driven-dissipative resonator, where the bound matches analytical results for the dynamics, and the infinite-range Ising model, where only the steady state is solvable and thus the bound provides information beyond the reach of existing analytical approaches. Our results can be applied to experimentally characterize ultrafast systems, and to theoretically analyze many-body models with dynamics that are analytically or numerically hard.

Górecki, Wojciech [INFN, Pavia]↗
DOE OSTI2025

A Route to Design Novel Functional Peptides by Applying a Denoising Diffusional Model to mRNA Display Libraries

In vitro directed evolution techniques, such as mRNA display, enable peptide ligand discovery and optimization. However, physical libraries that rely on a genetic code can only search a small fraction of sequence space due to inherent biases in the genetic code and experimental limitations. To address this challenge, denoising diffusion implicit models (DDIMs) are applied to generate novel peptide ligands against B‐cell lymphoma extra‐large (Bcl‐x L ), a key cancer target. Starting with high‐throughput sequencing data from previous selections, a DDIM is trained to produce novel sequences with high affinity binding. Experimental validation confirms that most generated sequences are functionally equivalent to the original library members for Bcl‐x L binding and demonstrated comparable binding kinetics and affinity relative to the wildtype and nearest original neighbors. Importantly, this approach generated rare sequences not easily accessible via mutation and directed evolution. These results indicate that DDIMs can complement and expand directed evolution data, efficiently exploring underrepresented regions of sequence space. This approach provides a broadly applicable framework for accelerating ligand discovery and optimizing molecular properties across diverse targets.

Qi, Pearl [Mork Family Department of Chemical Engi↗
DOE OSTI2026

Time Correlations from Steady-State Expectation Values

Recovering properties of correlation functions is typically challenging. On the one hand, experimentally, it requires measurements with a temporal resolution finer than the system’s dynamics. On the other hand, analytical or numerical analysis requires solving the system evolution. Here, we use recent results of quantum metrology with continuous measurements to derive general lower bounds on the relaxation and second-order correlation times that are both easy to calculate and measure. These bounds are based solely on steady-state expectation values and their derivatives with respect to a system parameter, and can be readily extended to the autocorrelation of arbitrary observables. We validate our method on two examples of critical quantum systems: a critical driven-dissipative resonator, where the bound matches analytical results for the dynamics, and the infinite-range Ising model, where only the steady state is solvable, and thus the bound provides information beyond the reach of existing analytical approaches. Our results can be applied to the experimental characterization of ultrafast systems and to the theoretical analysis of many-body models whose dynamics are hard to compute.

Górecki, Wojciech [INFN, Pavia] (ORCID:00000001991↗
DOE OSTI2023

Modeling shock-induced void collapse in single-crystal Ta systems at the mesoscales

Understanding the role of microstructural heterogeneities on the shock wave propagation and defect evolution behavior is essential to predicting the dynamic response of metals. Heterogeneities, such as voids, provide challenges to understanding the wave propagation behavior as the shock-void interaction can collapse the void and result in large plastic strains and significant heating (hot spot formation) in the metal. Accurate modeling of this phenomenon requires predicting the void collapse mechanisms and the related heat generation and dissipation mechanisms (hotspot formation) that determine the microstructure evolution. While molecular dynamics (MD) simulations can model the void collapse behavior, the time/length scale capabilities pose a challenge to connect with continuum models or the experimental scales. Here, this study presents the capability of the newly developed quasi-coarse-grained dynamics (QCGD) method that extends the MD simulations to larger system sizes and longer times to model this phenomenon. This study uses QCGD simulations to investigate the mechanisms of shock wave interactions with pre-existing voids in single-crystal Ta microstructures for variations in shock pressures, void size, and loading orientations. For a given orientation, the plasticity contributions and rates of void collapse are observed to vary with shock pressures and void size. The larger void sizes and higher pressures result in increased temperatures (hot spots) and faster void collapse rates and unravel the variations in the plasticity contributions. In addition, QCGD simulations investigate the post-collapse microstructure evolution as a release wave travels through the hot spot region. The simulations reveal that the reduction in temperatures due to heat dissipation initiates the dynamic recrystallization behavior in the hotspot regions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Model Parameter Development for Complex Materials: Species-Specific Diffusion Barriers in 316 Stainless Steel from Systematic DFT Calculations

Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computational evolution of high-performing unfused non-fullerene acceptors for organic solar cells

Materials optimization for organic solar cells (OSCs) is a highly active field, with many approaches using empirical experimental synthesis, computational brute force to screen a subset of chemical space, or generative machine learning methods that often require significant training sets. While these methods may find high-performing materials, they can be inefficient and time-consuming. Genetic algorithms (GAs) are an alternative approach, allowing for the “virtual synthesis” of molecules and a prediction of their “fitness” for some property, with new candidates suggested based on good characteristics of previously generated molecules. In this work, a GA is used to discover high-performing unfused non-fullerene acceptors (NFAs) based on an empirical prediction of power conversion efficiency (PCE) and provides design rules for future work. The electron-withdrawing/donating strength, as well as the sequence and symmetry, of those units are examined. The utilization of a GA over a brute-force approach resulted in speedups up to 1.8 × 10 12 . New types of units, not frequently seen in OSCs, are suggested, and in total 5426 NFAs are discovered with the GA. Of these, 1087 NFAs are predicted to have a PCE greater than 18%, which is roughly the current record efficiency. While the symmetry of the sequence showed no correlation with PCE, analysis of the sequence arrangement revealed that higher performance can be achieved with a donor core and acceptor end groups. Future NFA designs should consider this strategy as an alternative to the current A-D-A'-D-A architecture.

14 SOLAR ENERGY↗
DOE OSTI2021

In-Situ Calibrated Modeling of Residual Stresses Induced in Machining under Various Cooling and Lubricating Environments

Although many functional characteristics, such as fatigue life and damage resistance depend on residual stresses, there are currently no industrially viable ‘Digital Process Twin’ models (DPTs) capable of efficiently and quickly predicting machining-induced stresses. By leveraging advances in ultra-high-speed in-situ experimental characterization of machining and finishing processes under plane strain (orthogonal/2D) conditions, we have developed a set of physics-based semi-analytical models to predict residual stress evolution in light of the extreme gradients of stress, strain and temperature, which are unique to these thermo-mechanical processes. Initial validation trials of this novel paradigm were carried out in Ti-6Al4V and AISI 4340 alloy steel. A variety dry, cryogenically cooled and oil lubricated conditions were evaluated to determine the model’s ability to capture the tribological changes induced due to lubrication and cooling. The preliminarily calibrated and validated model exhibited an average correlation of better than 20% between the predicted stresses and experimental data, with calculation times of less than a second. Based on such fast-acting DPTs, the authors envision future capabilities in pro-active surface engineering of advanced structural components (e.g., turbine blades).

36 MATERIALS SCIENCE↗
DOE OSTI2022

Double echo symmetry-based REDOR and RESPDOR pulse sequences for proton detected measurements of heteronuclear dipolar coupling constants

1 H{X} symmetry-based rotational echo double resonance pulse sequences (S-REDOR) and symmetry-based rotational echo saturation pulse double resonance (S-RESPDOR) solid-state NMR experiments have found widespread application for 1 H detected measurements of difference NMR spectra, dipolar coupling constants, and internuclear distances under conditions of fast magic angle spinning (MAS). In these experiments the supercycled $\text{R4}^2_1 (\text{SR4}^2_1)$ symmetry-based recoupling pulse sequence is typically applied to the 1 H spins to reintroduce heteronuclear dipolar couplings. However, the timing of $\text{SR4}^2_1$ and other symmetry-based pulse sequences must be precisely synchronized with the rotation of the sample, otherwise, the evolution of 1 H CSA and other interactions will not be properly refocused. For this reason, significant distortions are often observed in experimental dipolar dephasing difference curves obtained with S-REDOR or S-RESPDOR pulse sequences. Here we introduce a family of double echo (DE) S-REDOR/S-RESPDOR pulse sequences that function in an analogous manner to the recently introduced $t_1$-noise eliminated (TONE) family of dipolar heteronuclear multiple quantum coherence (D-HMQC) pulse sequences. Through numerical simulations and experiments the DE S-REDOR/S-RESPDOR sequences are shown to provide dephasing difference curves similar to those obtained with S-REDOR/S-RESPDOR. However, the DE sequences are more robust to the deviations of the MAS frequency from the ideal value that occurs during typical solid-state NMR experiments. In this work, the DE sequences are shown to provide more reliable 1 H detected dipolar dephasing difference curves for nuclei such as 15 N (with isotopic labelling), 183 W and 35 Cl. The double echo sequences are therefore recommended to be used in place of conventional S-REDOR/S-RESPDOR sequences for measurement of weak dipolar coupling constants and long-range distances.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Back to functional hydrides: Effects of neutron-irradiated microstructure on hydrogen retention in yttrium hydride

Functional hydrides are promising candidates for advanced nuclear reactors, particularly in portable or transportable applications, due to their high hydrogen-retention capabilities, enabling efficient neutron moderation, and compact reactor design. However, hydrogen mobility in hydrides at elevated irradiation temperatures poses significant technological challenges, necessitating a comprehensive understanding of their irradiation behavior. Furthermore, this study investigated the microstructural and chemical stability of neutron-irradiated yttrium hydrides to assess their hydrogen-retention capacity. A targeted literature review was also conducted to contextualize neutron-irradiation effects on functional hydrides with regards to structural stability and hydrogen retention. Experimental characterizations revealed that, at high temperatures, irradiated hydrides retained their phase stability, which was likely enhanced by irradiation-induced microstructure evolution. Notably, an amorphous yttrium and oxygen -rich surface layer was present at the free surface of the hydride. Its thickness decreased while a continuous crystalline Y-O-rich layer was formed with increasing neutron damage. Additionally, the number density of dislocation loops and cavities generally increased as a function of neutron dose. First-principles calculations of hydrogen behavior within yttrium vacancy clusters in yttrium hydrides demonstrated vacancy-size-dependent hydrogen stability and configuration, highlighting the role of vacancy geometry in regulating hydrogen retention. Thus, the presence of irradiation-induced dislocation loops and cavities were hypothesized to improve hydrogen retention. Collectively, these findings advance the understanding of hydride behavior under neutron irradiation as well as their technological readiness for portable or transportable nuclear reactors.

Defect clusters↗
DOE OSTI2022

Detailed modeling of hydrogen release and particle shrinkage during pyrolysis of inhomogeneous wood

Hydrogen release during pyrolysis of woody biomass is studied considering anisotropicity and inhomogeneity of wood structure. A new anisotropic shrinkage model is proposed based on the decomposition of main wood constituents, i.e., cellulose, hemicellulose, and lignin. The new shrinkage model can predict the temporal evolution of the wood structure, and the differences between axial and radial shrinkage during pyrolysis. The model agrees very well with several experimental data from the literature. Based on particle temperature during conversion, the pyrolysis is partitioned into four stages, and the hydrogen release and H 2 formation from each stage are investigated. Stage (IV) of pyrolysis, from 1000 to 1273 K, is found to be efficient for H 2 production owing to the production of considerable mass of H with a minimal amount of tar species. Furthermore, the char quality is found to be different at the end of stages (II), (III), and (IV), where around 67.7, 80.5, and 93.4% wt. of solid residue is made of carbon, respectively. The model is also used to explain how the heating rate affects the temperature distribution inside the particle and how it shifts the peak of hydrogen release. Finally, the pyrolysis of two inhomogeneous wood samples — a beech twig with bark and a beech dowel with growth rings — are investigated. The bark can affect the pyrolysis rate, products, and flow pattern inside the particle. The growth rings do not have a considerable effect on the pyrolysis rate and products, but they have a significant impact on the flow pattern. This has an important implication for char conversion studies where the internal surface area and porosity field distribution have a significant effect on the gasification and oxidation rates.

09 BIOMASS FUELS↗
DOE OSTI2021

Electron–Nuclear Dynamics Accompanying Proton-Coupled Electron Transfer

Although photoinduced proton-coupled electron transfer (PCET) plays an essential role in photosynthesis, a full understanding of the mechanism is still lacking due to the complex nonequilibrium dynamics arising from the strongly coupled electronic and nuclear degrees of freedom. Here we report the photoinduced PCET dynamics of a biomimetic model system investigated by means of transient IR and two-dimensional electronic–vibrational (2DEV) spectroscopies, IR spectroelectrochemistry (IRSEC), and calculations utilizing long-range-corrected hybrid density functionals. This collective experimental and theoretical effort provides a nuanced picture of the complicated dynamics and synergistic motions involved in photoinduced PCET. In particular, the evolution of the 2DEV line shape, which is highly sensitive to the mixing of vibronic states, is interpreted by accurate computational modeling of the charge separated state and is shown to represent a gradual change in electron density distribution associated with a dihedral twist that occurs on a 120 fs time scale.

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

3D simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK

In recent years, the nonlinear 3D magnetohydrodynamic codes JOREK, M3D-C1, and NIMROD developed the capability of modeling realistic 3D vertical displacement events (VDEs) including resistive walls. In this paper, a comprehensive 3D VDE benchmark is presented between these state-of-the-art codes. The simulated case is based on an experimental NSTX plasma but with a simplified rectangular wall. There are differences between the physics models and numerical methods, and the VDE evolution leads to sensitivities on the initial conditions that cannot be avoided as can be done in edge localized modes (ELM) and sawtooth simulations (due to the non-cyclical nature of VDEs). Nonetheless, the comparison serves to quantify the level of agreement in the relevant quantities used to characterize disruptions, such as the 3D wall forces and energy decay. The results bring confidence regarding the use of the mentioned codes for disruption studies, and they distinguish aspects that are specific to the models used (e.g., reduced vs full MHD models). Here, the simulations show important 3D fea- tures for a NSTX plasma, such as the self-consistent evolution of the halo current and the origin of the wall forces. In contrast to other reduced MHD models based on an ordering in the aspect ratio, the ansatz-based JOREK reduced MHD model allows capturing many aspects of the 3D dynamics even in the spherical tokamak limit considered here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Determination of Thermal Expansion, Defect Formation Energy, and Defect-Induced Strain of α-U Via ab Initio Molecular Dynamics

Uranium (U) is often alloyed with molybdenum (Mo) or zirconium (Zr) in order to stabilize its high-temperature body-centered cubic phase for use in nuclear reactors. However, in all metallic fuel forms, the α phase of U remains in some fraction. This phase decomposition due to temperature or compositional variance can play an outsized role on fuel performance and microstructural evolution. Relatively little is known about fundamental point defect properties in α -U at non-zero temperatures, from either computational or experimental studies. This work performs the first thorough evaluation of the α phase of U via ab initio molecular dynamics (AIMD). A number of thermophysical properties are calculated as a function of temperature, including equilibrium lattice parameters, thermal expansion, and heat capacity. These results indicate a two-region behavior, with the transition at 400 K. The thermal expansion/contraction in the a / b direction occurs rapidly from 100 up to 400 K, after which a more linear and gradual change in the lattice constant takes place. The volumetric expansion matches experiments quantitatively, but the individual lattice constant expansion only matches experiments qualitatively. Point defect formation energies and induced lattice strains are also determined as a function of temperature, providing insight on defect populations and the fundamentals of irradiation growth in α -U. Interstitials induce significantly more strain than vacancies, and the nature of that strain is highly dependent on the individual lattice directions. The direction of point defect-induced lattice strain is contrary to the irradiation growth behavior of α -U. This work shows that isolated point defects cannot be the primary driving force responsible for the significant irradiation-induced growth of α -U observed experimentally.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗