Search NASA⌕ Search

SEARCH · Search NASA

Results for “predictive 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Peak-Seeking Optimization of Trim for Reduced Fuel Consumption: Architecture and Performance Predictions

A peak-seeking control approach for real-time trim configuration optimization for reduced fuel consumption has been developed by researchers at the National Aeronautics and Space Administration (NASA) Dryden Flight Research Center to address the goals of the NASA Environmentally Responsible Aviation project to reduce fuel burn and emissions. The peak-seeking control approach is based on a steepest-descent algorithm using a time-varying Kalman filter to estimate the gradient of a performance function of fuel flow versus control surface positions. In real-time operation, deflections of symmetric ailerons, trailing-edge flaps, and leading-edge flaps of an F/A-18 airplane (McDonnell Douglas, now The Boeing Company, Chicago, Illinois) are controlled for optimization of fuel flow. This paper presents the design and integration of this peak-seeking controller on a modified NASA F/A-18 airplane with research flight control computers. A research flight was performed to collect data to build a realistic model of the performance function and characterize measurement noise. This model was then implemented into a nonlinear six-degree-of-freedom F/A-18 simulation along with the peak-seeking control algorithm. With the goal of eventual flight tests, the algorithm was first evaluated in the improved simulation environment. Results from the simulation predict good convergence on minimum fuel flow with a 2.5-percent reduction in fuel flow relative to the baseline trim of the aircraft.

performance optimization↗
NASA NTRS2013

Peak-Seeking Optimization of Trim for Reduced Fuel Consumption: Architecture and Performance Predictions

A peak-seeking control approach for real-time trim configuration optimization for reduced fuel consumption has been developed by researchers at the National Aeronautics and Space Administration (NASA) Dryden Flight Research Center to address the goals of the NASA Environmentally Responsible Aviation project to reduce fuel burn and emissions. The peak-seeking control approach is based on a steepest-descent algorithm using a time-varying Kalman filter to estimate the gradient of a performance function of fuel flow versus control surface positions. In real-time operation, deflections of symmetric ailerons, trailing-edge flaps, and leading-edge flaps of an FA-18 airplane (McDonnell Douglas, now The Boeing Company, Chicago, Illinois) are controlled for optimization of fuel flow. This presentation presents the design and integration of this peak-seeking controller on a modified NASA FA-18 airplane with research flight control computers. A research flight was performed to collect data to build a realistic model of the performance function and characterize measurement noise. This model was then implemented into a nonlinear six-degree-of-freedom FA-18 simulation along with the peak-seeking control algorithm. With the goal of eventual flight tests, the algorithm was first evaluated in the improved simulation environment. Results from the simulation predict good convergence on minimum fuel flow with a 2.5-percent reduction in fuel flow relative to the baseline trim of the aircraft.

peak-seeking↗
DOE OSTI2026

Maximum Energy of Particles Accelerated in Gamma-Ray Burst Afterglow Shocks

Particle acceleration in relativistic collisionless shocks remains an open problem in high-energy astrophysics. Particle-in-cell (PIC) simulations predict that electron acceleration in weakly magnetized shocks proceeds via small-angle scattering, leading to a maximum electron energy significantly below the Bohm limit. This upper bound on electron energy manifests observationally as a characteristic synchrotron cutoff, providing a direct probe of the underlying acceleration physics. Gamma-ray burst (GRB) afterglows offer an exceptional laboratory for testing these predictions. Here, we model the spectral evolution of GRB afterglows during the relativistic deceleration phase, incorporating PIC-motivated acceleration prescriptions and self-consistently computing synchrotron and synchrotron self-Compton emission. We find that low-energy bursts in low-density environments, typical of short GRBs, exhibit a pronounced synchrotron cutoff in the GeV band within minutes to hours after the trigger. Applying our framework to GRB 190114C and GRB 130427A, we find that current observations are insufficient to discriminate between PIC-motivated acceleration and the Bohm limit, primarily due to poor photon statistics in the Fermi-Large Area Telescope band. Nevertheless, future MeV–TeV afterglow observations can break model degeneracies and place substantially tighter constraints on the mechanisms responsible for particle acceleration in relativistic shocks. To this end, we simulate a fiducial nearby short GRB as a promising probe of the cutoff location, for which the two acceleration scenarios are cleanly distinguishable and the detection of such an event in the near future remains feasible.

Wu, Zhao-Feng [Purdue University, West Lafayette, ↗
DOE OSTI2025

Effects of high-pressure hydrogen exposure on filler-elastomer adhesion

Elastomers are known to gain enhanced mechanical properties through compounding with nanosized filler particles such as silica or carbon black. Filler dispersion and filler-polymer interfacial strength are key contributing factors to this improvement. The interfacial strength is critical to part lifetime in pressurized gas sealing applications such as O-rings, where weak binding between the filler particle and polymer matrix can lead to internal void structures. With the aim to build a fundamental understanding of precursors to pressurized hydrogen-induced failure in elastomers, we use all-atom molecular dynamics simulations to study the impact of hydrogen oversaturation on filler-polymer interaction strength. We systematically study the interface between a commonly used elastomer, ethylene-propylene-diene monomer (EPDM) and silica by varying gas concentration, crosslink density, and surface chemistry. Our simulations predict that decompression leads to a localization of excess gas near the interface. In conclusion, we demonstrate that this localized gas can weaken interfacial adhesion and quantify the interaction using thermodynamic approaches.

EPDM↗
DOE OSTI2026

Global energy balance simulations of shattered pellet injection

The dynamics of shattered pellet injection (SPI) shutdowns are simulated using a time-dependent global energy balance model, based on a modification of the KPRAD framework. The new SPI particle source in the model calculates the ablation of individual pellet fragments that enter the plasma as a temporally resolved plume, thus capturing the effects of earlier fragments on the ablation of those that follow, which has a significant impact on the overall assimilation. Despite the reduced physics and the global averaging of all quantities, results from a large number of DIII-D, KSTAR, and JET experiments are well reproduced, including the plasma cooling timescales, particle assimilation, and current quench (CQ) rates. Cooling timescales and CQ rates are in good agreement for pellets containing as little as ∼1% neon, while particle assimilations are most accurate for neon fractions above ∼15% by number of atoms. Below this, the assimilation tends to be overestimated due to the lack of radial transport in the particle balance, which becomes important in the low-Z limit. Predictive simulations of mixed-composition dual-SPI shutdowns in ITER are compared against those with the 3D non-linear magnetohydrodynamic code JOREK and are found to reproduce overall trends observed in the higher-fidelity modeling across a range of injection scenarios. The general success of the model points to the critical role of energy balance in determining SPI particle assimilation and the subsequent disruption dynamics and highlights the value of these simulations for experimental interpretation and for optimizing the deployment of computationally expensive, higher-fidelity models.

energy balance model↗
DOE OSTI2026

Cryogenic-Refined MOSFET Modeling for Oscillator, Frequency Divider, and Amplifier Designs Below 4 K

Capturing device characteristic changes at cryogenic temperatures is crucial for cryo-CMOS circuit designs. In this work, we present an isothermal cryogenic-refined modeling approach for CMOS transistors that is simple, low overhead, and easy to implement while offering the required accuracy for predicting circuit performance at the designated temperatures. Guided by die-level measurement data and circuit design principles, the model introduces corrections to only five critical parameters: threshold voltage, carrier mobility, elevated low-frequency flicker noise, dominant high-frequency shot noise, and subthreshold swing (SS). These refinements are implemented around the foundry-provided SPICE model, which is typically validated only down to about 200 K. With these adjustments, the proposed cryogenic-refined model achieves less than 5% error in both large-signal metrics (I–V characteristics) and small-signal parameters (e.g., transconductance) when compared with device measurements at deep-cryogenic temperatures. The methodology is validated in two advanced technologies: TSMC 40-nm CMOS and GlobalFoundries (GF) 45-nm RF-SOI. We further demonstrate its applicability in three representative RF circuits: a 30-GHz LC oscillator, a high-speed current-mode-logic (CML) frequency divider (FD), and a subthreshold Gb/s amplifier, all showing close agreement between simulated predictions and measurements performed at 4 and 2.5 K. Finally, we believe that the proposed approach is implementation-friendly and can significantly accelerate the development of cryo-CMOS integrated circuits.

circuit modeling↗
DOE OSTI2024

Pipeline Hydrogen Decarbonization and Repurposing Analyzers (P-HyDRAs)

The Pipeline Hydrogen Decarbonization and Repurposing Analyzers (P-HyDRAs) are a set of prototype computational tools for simulating and optimizing midstream natural gas pipeline system operations subject to location and time-dependent hydrogen blending. The models can accurately resolve dynamic gas flows through large-scale pipeline networks using non-ideal gas equations of state. The codes can be used as decision support for planning and design decisions involving intra-day energy flow schedules as well as spatiotemporal economic values of natural gas, hydrogen, and net energy delivered to consumers while ensuring that pipeline hydraulic limitations, gas compressor station constraints, operational factors, and pre-existing shipping contracts are satisfied. The inputs to the codes are a model of the pipeline system as well as time-series data that specify boundary conditions on the network. For optimization, the code module requires price and quantity offers for natural gas and hydrogen and price and quantity bids for energy, which are used as time-dependent constraints in an optimal control problem. The outputs are time-series data that provide a predictive simulation of gas flows, mass fractions, and pressures, or with additional degrees of freedom give an approximately optimal solution for gas injections/withdrawals, compressor settings, and sensitivities to the objective function that provide locational values of energy.

Zlotnik, Anatoly↗
DOE OSTI2025

Revealing the Phonon Bottleneck Limit in Negatively Charged CdS Quantum Dots

The capture of photoexcited hot electrons in semiconductors before they lose their excess energy to cooling is a long-standing goal in photon energy conversion. Semiconductor nanocrystals have large electron energy spacings that are expected to slow down electron relaxation by phonon emission, but hot electrons in photoexcited nanocrystals, nevertheless, cool rapidly by energy transfer to holes. This makes the intrinsic phonon-bottleneck-limited hot electron lifetime in nanocrystals elusive. For this work, we used a combination of theory and experiments to probe the hot-electron dynamics of negatively charged cadmium sulfide (CdS) colloidal quantum dots (QDs) in the absence of holes. Experiments found that these hot electrons cooled on a 100 ps time scale. Theoretical simulations predicted that pure phonon-bottleneck-limited hot electron cooling occurs on a similar time scale. This similarity suggests that the experimental measurements reflect the upper limit on the hot-electron lifetimes in these CdS QDs and the lower limit on the rates of processes that can harvest those hot electrons.

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

Substrate-Directed Dimensional and Phase Control of Peptide Assemblies on Two-Dimensional van der Waals Materials

Understanding and controlling biomolecular self-assembly on van der Waals (vdW) materials has the potential to advance hybrid bioelectronic devices by enabling precise tuning of the interface and modulation of the resulting electronic properties of the biomolecule-vdW heterostructure. However, how surface properties of vdW materials direct biomolecule assembly remains poorly understood. To fill this knowledge gap, we investigated the assembly of a peptide known to assemble into two-dimensional (2D) crystalline films on MoS 2 on three representative vdW surfaces: WS 2 , MoS 2 , and highly oriented pyrolytic graphite (HOPG). Using in situ atomic force microscopy (AFM), we find that assembly is substrate-dependent, resulting in multilayers on WS 2 , monolayers on MoS 2 , and multiple coexisting phases on HOPG. WS 2 exhibits a higher negative charge, strong long-range electrostatic interactions, and extensive hydration layering that may promote multilayer stacking. In contrast, MoS 2 has stronger short-range interactions with the peptides but much weaker long-range interactions and hydration structure, which may favor monolayer formation. Molecular dynamics simulations predict a corresponding switch from monolayer to multilayer aggregates of the adsorbed monomers, reflected in their relative mobilities. On hydrophobic HOPG, the peptides bind most strongly and remain as monomers with high surface mobility. The peptide dimers comprising the basic unit of the crystals are more compact on HOPG, which has a smaller lattice constant than WS 2 or MoS 2 , suggesting strain contributes to stabilizing multiple phases. Our results provide mechanistic insights into how surface charge and hydration structure, and the lattice structure of the substrates governs peptide assembly on vdW materials, offering a framework to rationally control the 2D peptide-vdW heterostructures.

Molecular dynamics simulations↗
DOE OSTI2024

Lateral Confinement in 2D Nanoplatelets: A Strategy to Expand the Colloidal Quantum Engineering Toolbox

Among colloidal nanocrystals, 2D nanoplatelets offer a unique set of properties with exceptionally narrow luminescence and low lasing thresholds. Furthermore, their anisotropic shape expands the playground for the complex design of heterostructures where spectra but also scattering rates can be engineered. A challenge that still remains is to combine shell growth which makes NPLs stable, with spectral tunability. Indeed, most reported shelled nanoplatelets end up being red emitters due to a loss of quantum confinement. Here, the combination of both lateral and in-plane confinements within a single heterostructure is explored. A CdS/CdSe/CdS/CdZnS core–crown–crown shell structure that enables yellow emission is grown and that is responsive to a large range of excitation including visible photons, X-ray photons, electron beams, and electrical excitations. k.p simulations predict that emission tunability of up to several 100 s of meV can be obtained in ideal structures. This material also displays stimulated emission resulting from bi-exciton emission with a low threshold. Once integrated into an LED stack, this material is compatible with sub-bandgap excitation and exhibits high luminance. Scaling of the electroluminescence properties by downsizing the pixel size is also investigated.

2D materials↗
DOE OSTI2026

Process-driven roadmap for depositing super duplex stainless steel via wire Arc additive manufacturing

Here, this study systematically investigates the effects of shielding gas, bead spacing, weld mode, and travel speed on the phase balance, porosity, and hardness of wire arc additively manufactured (WAAM) super duplex stainless steel (ER2594). Robotic WAAM was employed to fabricate multilayer walls under systematically varied process conditions, followed by phase transformation simulations, X-ray computed tomography (XCT), electron backscatter diffraction (EBSD), and microhardness evaluation. Thermodynamic simulations predicted rapid cooling of AM process can suppress the potential formation of deleterious precipitates which was later validated via cross-sectional microstructure analyses of printed samples. XCT revealed porosity levels below 0.2% for all deposits, with reduced travel speed significantly lowering defect volume. Microstructural analyses revealed the evolution of various austenite precipitates, including grain boundary austenite (GBA), Widmanstätten austenite (WA), and intergranular austenite (IGA), sequentially upon cooling of the ferrite phase. Among all process parameters, weld transfer mode exhibited the strongest influence on phase balance; pulsed mode promoted higher ferrite retention (~ 36%) compared to RapidX mode. No consistent relationship between stabilized phase fraction and captured microhardness was observed. This work provides critical insights for optimizing WAAM parameters to control phase balance and mechanical performance, which is essential for producing high-integrity super duplex stainless-steel components for nuclear and marine applications.

Grain orientation↗
DOE OSTI2024

Monte Carlo Simulations of 347H Stainless Steel Aging for the Synthetic Generation of Microstructures Under Creep Conditions

Here, a Monte Carlo simulation method capable of replicating the kinetics of M 23 C 6 precipitation in 347H stainless steels was developed for the purpose of producing synthetic microstructures that approximate its microstructural evolution under aging periods of up to 10,000 hours at temperatures between 600 °C and 750 °C. To accomplish this, experimental data from the literature was used to parameterize simulations and replicate the nucleation and growth kinetics of M 23 C 6 particles within 347H and similar austenitic stainless steel alloys. These simulations were found to have considerable fidelity to previous efforts to study the precipitation of M 23 C 6 in other 300 series stainless steel alloys. Synthetic 347H microstructures were then generated that accounted the effects of aging temperature, duration, dislocation density, and the presence of boron within the microstructure. These simulations predict several key trends, those being that (1) the size of M 23 C 6 precipitates decreased with aging temperature and (2) the growth rate of M 23 C 6 particles decreased with aging temperature. Further, while (3) the addition of dislocation density due to creep conditions resulted in increasing intragranular nucleation of M 23 C 6 precipitates with increasing dislocation density and (4) B additions within the microstructure led to modest increases in precipitate size above 700 °C, which indicates that more complex physics are necessary to account for the presence of B.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Genetic algorithm optimization of nuclear criticality experiment for reduction of intermediate-energy 239 Pu nuclear data uncertainties

Nuclear criticality experiments are conducted to investigate specific nuclear data important for safe handling and storage of fissile materials, reactor design and operation, and the validation of radiation transport codes. Incorrect or uncertain nuclear data can prohibitively impact operational safety limits, reactor licensing, and predictive simulation capability; therefore, integral measurements from criticality experiments are necessary and should be performed frequently. To maximize the impact of the integral measurements, it is important to consider experiment geometry, material selection, and component dimensions. When taking these considerations into account, the experiment design process becomes iterative and very time intensive. This work utilizes a genetic algorithm to efficiently explore potential nuclear criticality experiment designs for the Laboratory Directed Research & Development project PARADIGM (PARallel Approach of Differential and InteGral Measurements) at Los Alamos National Laboratory. In this paper, the building blocks of the genetic algorithm are discussed in detail, the genetic algorithm methodology is verified, and the genetic algorithm is used to produce three candidate experiment models for the final PARADIGM design. The three candidate models produced by the genetic algorithm consist of copper-reflected assemblies containing 14 repeating units of alumina, graphite, boron, and plutonium plates. Furthermore, in addition to the optimization results, final design considerations are also discussed for designs with a height and/or weight very close to or slightly above assembly machine operational limits.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Experimental and theoretical investigation into the high pressure deflagration products of 2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

Diamond anvil cell (DAC) laser ignition experiments and reactive ab initio molecular dynamics (AIMD) simulations were performed on the high explosive (HE) LLM-105 to investigate its high pressure (HP) deflagration chemistry. Raman and optical spectroscopy measurements reveal LLM-105 reacts into an opaque carbonaceous product at 4–25 GPa. At pressures >~ 27 GPa, the reaction product consists of an amorphous optically transparent solid and nitrogen (N 2 ) in the solid phase. While not a one-to-one comparison due to the small time and length scales, the HP AIMD simulations show that some of the product is molecular N 2 , in qualitative agreement with experiment, while above 20 GPa most of the product consists of large amorphous C x H y N z O k clusters. Clustering is enhanced with pressure and reduces with temperature. In the experiments with initial sample pressure >~ 25 GPa, the pressure within the DAC decreases with minimal change in DAC cavity area. At initial sample pressures of 43.9 GPa, when quenched to 0 K, simulations predict a product experiencing a lower pressure consistent with the experimental measurement at lower load pressures. In conclusion, the results are important for understanding the HP deflagration chemistry of LLM-105.

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

Solid-state batteries enabled by ultra-high-frequency self-heating

Solid-state batteries (SSBs) are promising next-generation batteries due to their high energy density and enhanced thermal stability and safety. However, their sluggish kinetics and transport at room temperature result in high internal impedance and critically reduce the attainable discharge energy density. Taking advantage of their strong temperature-dependent ionic conductivity, here we introduce ultra-high-frequency (greater than 105 Hz) self-heating (UHFSH) of SSBs, which can rapidly warm up the batteries from room temperature to operating temperature (∼65°C) in less than a minute. As proof of concept, UHFSH experiments were conducted on symmetric solid-state cells with lithium aluminum germanium phosphate electrolyte in different configurations. Using an experimentally validated model, pack-level simulations predict fast heating (50 K/min) and minimized heating energy consumption (less than 4%). Without any modification of the materials or structure of the batteries, our non-intrusive self-heating strategy potentially enables SSBs to discharge more than 2-fold energy in 25°C ambient.

pulse heating↗
DOE OSTI2026

Enhancing microsegregation during rapid directional solidification through ternary microalloying

The nano-cellular dendritic microstructure formed during rapid directional solidification in powder bed fusion additive manufacturing creates unique properties such as simultaneous improvement in strength and ductility. However, process control of microsegregation features remains challenging due to low sensitivity of critical solidification mechanisms to process parameters. This study leverages microalloying to achieve large changes in dendrite composition, microstructure, and interdendritic zone width during laser powder bed fusion without modifying process parameters. CALPHAD simulations predict that the addition of Zr significantly steepens the solidus line of the dilute Cu-Cr alloy system, leading to enhanced Cr rejection into the melt and greater than 95% reduction in solubility of Cr in the solidified Cu matrix. Experimental validation using time-of-flight secondary ion mass spectrometry and Kelvin probe force microscopy reveals that the ternary alloy containing 0.01 wt% Zr exhibited wider interdendritic regions compared to the binary, a significantly higher number of Cr-rich particles within interdendritic regions, near-complete ejection of oxygen impurities from the matrix, and greater nanoscale work function contrast. These features indicate more aggressive Cr segregation in the presence of Zr and a purer Cu matrix and provide a potentially robust method for engineering the nano-cellular dendritic solidification microstructure.

CALPHAD↗
DOE OSTI2024

High strain-rate strength response of single crystal tantalum through in-situ hole closure imaging experiments

The properties of crystalline materials often depend on directionality and operating conditions. Specifically, the strength of materials can depend anisotropically on crystal direction and the loading condition. To probe these effects, a preliminary series of high strain-rate (> 105/s) strength plate-impact hole closure experiments were performed on high purity single crystal Tantalum cubes. The orientation of the single crystals with respect to impact/loading were varied to provide data to inform crystal plasticity modeling efforts. The experiments consist of in-situ high-resolution X-ray radiographic imaging of the hole collapse under dynamic compression conditions to infer the material strength via its resistance to closure at increasing levels of plastic strain. The experiments are compared against hydrocode simulation predictions. Here, a comparison with simple elastic perfectly plastic strength model predictions is presented to elucidate the response of the different crystal orientations at high strain-rate and large plastic strains.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Microstructural evolution and phase stability in Nb-containing interstitial Fe-Mn-Co-Cr-C high-entropy alloys: An in-situ synchrotron X-ray diffraction study during laser melting

The influence of Nb on phase stability and microstructural evolution in an interstitial Fe-Mn-Co-Cr-C high-entropy alloy was investigated using in-situ synchrotron X-ray diffraction (SXRD) during laser melting. Scheil-Gulliver simulations predict the formation of σ and γ-f.c.c. phases in all three alloys, along with NbC in Nb-containing compositions. SXRD confirmed the presence of most predicted phases, but the σ phase was absent. Nb promotes crystallite refinement and increases dislocation density, though excessive additions reduce refinement efficiency due to solubility limits and secondary phase formation. Furthermore, Nb addition also enhances ε-h.c.p. phase formation by reducing stacking fault energy through NbC-induced carbon depletion. Analysis of intensity peak evolution reveals that Nb alters preferred grain orientations, reducing {111} γ intensity while enhancing {220} γ , leading to a more isotropic grain distribution. Overall, Nb plays a key role in phase selection, microstructure refinement, and preferred orientation evolution, allowing the tailored microstructure of high-entropy alloys via rapid solidification.

Interstitial high entropy alloys↗