Search NASA⌕ Search

SEARCH · Search NASA

Results for “Space simulators”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Electron Influence on the Parallel Proton Firehose Instability in 10-moment, Multifluid Simulations

Instabilities driven by pressure anisotropy play a critical role in modulating the energy transfer in space and astrophysical plasmas. For the first time, we simulate the evolution and saturation of the parallel proton firehose instability using a multifluid model without adding artificial viscosity. These simulations are performed using a 10-moment, multifluid model with local and gradient relaxation heat-flux closures in high-β proton–electron plasmas. When these higher-order moments are included and pressure anisotropy is permitted to develop in all species, we find that the electrons have a significant impact on the saturation of the parallel proton firehose instability, modulating the proton pressure anisotropy as the instability saturates. Even for lower β's more relevant to heliospheric plasmas, we observe a pronounced electron energization in simulations using the gradient relaxation closure. Our results indicate that resolving the electron pressure anisotropy is important to correctly describe the behavior of multispecies plasma systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

TCAD-Machine Learning Enabled TID Compact Model Development for Commercial SiC MOSFET

We propose a TCAD (Technology Computer Aided Design)-machine learning coupled approach that combines a TCAD tool (Charon), optimization/uncertainty quantification tool (Dakota), surrogate models, and Bayesian learning capabilities. The coupling approach is used for accurate modeling and calibration of total ionizing dose (TID) induced threshold voltage (V th ) shifts in Commercial-Off-The-Shelf (COTS) semiconductor devices and to develop physics-informed TID compact models. This versatile approach is applied to model the TID effect in an exemplar COTS 3.3 kV SiC power MOSFET (Metal-Oxide-Semiconductor Field-Effect Transistor). With the Charon-Dakota coupling, we can determine key device geometry and doping values based on device physics, which are difficult to obtain or not available for COTS devices but important for TCAD simulation; additionally, we can efficiently generate thousands of simulation results in a large parameter space, which makes it possible to develop data-driven surrogate models and perform Bayesian calibration. Utilizing the full tool-coupling approach, we achieve calibrated TCAD simulation models that accurately capture the average TID-induced V th shifts behavior with total doses and V th shifts saturation at high doses as observed in experimental data. More importantly, the calibrated TCAD simulations are obtained with determined TID model parameters (e.g., hole trap density and capture cross section) values that contain well quantified uncertainties. Furthermore, we can isolate and quantify the noises that are not captured by the TCAD models but exist in the measured data due to measurements and devices variabilities. Lastly, the calibrated surrogate models are used to develop physics-informed TID compact models. The method is generalizable to other devices and/or radiation conditions with few modifications and can provide well-determined uncertainties.

COTS↗
DOE OSTI2021

PSTN-051: Survey Strategy and Cadence Choices for the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST)

A summary of survey strategy and cadence choices, simulated and evaluated by the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) Scheduler Team, prepared for the Survey Cadence and Optimization Committee (SCOC). A large telescope survey, covering the entire visible sky repeatedly every few days in multiple bandpasses over the course of ten years, is the core idea of the LSST. An area of about 20,000 square degrees observed under a wide range of conditions to deep coadded limiting magnitudes in bandpasses ugrizy enables cosmological studies and studies of the Milky Way structure with unprecedented precision; the same survey, when cadenced well, can serve to open new windows into our understanding of transient and variable stars, and extend our knowledge of small bodies throughout the Solar System by orders of magnitude. The outlines of, and some basic necessary requirements for these goals are outlined in the LSST Science Requirements Document (SRD). Finding options for the observing strategy to meet more detailed needs of an even wider range of science goals, as well as building the LSST Scheduler and Metrics Analysis Framework, has been the work of the LSST Scheduler Team with support and input from the astronomical community, including the Community Observing Strategy Evaluation Paper (COSEP), the Call for White Papers, numerous performance metrics, and guidance from the LSST Science Advisory Committee in their Recommendations for Operations Simulator Experiments. This report will help enable the SCOC to make its recommendation for the preferred survey strategy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A velocity space hybridization-based Boltzmann equation solver

In the present research, a new method for simulation of rarefied gas flows is proposed, a velocity-space hybrid of both a DSMC representation of particles and a discrete velocity quasi-particle representation of the distribution function. The hybridization scheme is discussed in detail, and is numerically verified for two test-cases: the BKW relaxation problem and a stationary Maxwellian distribution. It is demonstrated that such a velocity-space hybridization can provide computational benefits when compared to a pure discrete velocity method or pure DSMC approach, while retaining some of the more attractive properties of discrete velocity methods. Further possible improvements to the velocity-space hybrid approach are discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Film Thickness and Composition Effects in Symmetric Ternary Block Copolymer/Homopolymer Blend Films: Domain Spacing and Orientation

We report thickness and homopolymer molar mass dependencies in thin film blends of a compositionally symmetric 75 kg/mol lamellar polystyrene-block-poly(methyl methacrylate) diblock copolymer (PS-b-PMMA) with near-equal volume fractions of PS and PMMA homopolymers on “neutral” surfaces that promote a vertical domain orientation. PS and PMMA homopolymers with equal molar masses of 1, 3, and 22 kg/mol are blended in a 1:1 mass ratio. To efficiently explore this parameter space, we prepared combinatorial samples with gradient homopolymer concentrations using electrospray deposition (ESD) and gradient film thicknesses using flow coating; these samples were subsequently characterized by grazing-incidence small-angle X-ray scattering and scanning electron microscopy, respectively. Our results demonstrate that the added homopolymers increase or reduce domain spacings with respect to the neat (unblended) block copolymer depending on their molar masses, in line with previous reports; a simple heuristic to estimate the domain scaling in ternary blends based on the homopolymer molar mass agrees well with data for the thinnest blend films studied here. More surprisingly, we observe that the measured domain spacing in blends also depends on film thickness, with thicker films exhibiting larger domain spacings than thinner films at the same homopolymer concentration. Furthermore, thickness gradients at fixed homopolymer concentrations reveal a change in lamellae orientation from vertical to horizontal as film thickness is increased beyond the nominal lamellar domain spacing. Coarse-grained molecular dynamics simulations indicate that this change in orientation is induced by homopolymer segregation to film surfaces. Finally, tying these effects together, we hypothesize that the observed dependence of domain spacing on film thickness is a consequence of changes in the vertical homopolymer composition profile in thicker films during lamellae reorientation and the conformational asymmetry between PMMA and PS segments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Provably accurate simulation of gauge theories and bosonic systems

Quantum many-body systems involving bosonic modes or gauge fields have infinite-dimensional local Hilbert spaces which must be truncated to perform simulations of real-time dynamics on classical or quantum computers. To analyze the truncation error, we develop methods for bounding the rate of growth of local quantum numbers such as the occupation number of a mode at a lattice site, or the electric field at a lattice link. Our approach applies to various models of bosons interacting with spins or fermions, and also to both abelian and non-abelian gauge theories. We show that if states in these models are truncated by imposing an upper limit Λ on each local quantum number, and if the initial state has low local quantum numbers, then an error at most ϵ can be achieved by choosing Λ to scale polylogarithmically with ϵ − 1 , an exponential improvement over previous bounds based on energy conservation. For the Hubbard-Holstein model, we numerically compute a bound on Λ that achieves accuracy ϵ , obtaining significantly improved estimates in various parameter regimes. We also establish a criterion for truncating the Hamiltonian with a provable guarantee on the accuracy of time evolution. Building on that result, we formulate quantum algorithms for dynamical simulation of lattice gauge theories and of models with bosonic modes; the gate complexity depends almost linearly on spacetime volume in the former case, and almost quadratically on time in the latter case. We establish a lower bound showing that there are systems involving bosons for which this quadratic scaling with time cannot be improved. By applying our result on the truncation error in time evolution, we also prove that spectrally isolated energy eigenstates can be approximated with accuracy ϵ by truncating local quantum numbers at Λ = polylog ( ϵ − 1 ) .

Tong, Yu↗
DOE OSTI2021

Manipulation and statistical analysis of the fluid flow of polymer semiconductor solutions during meniscus-guided coating

Recent work in structure–processing relationships of polymer semiconductors have demonstrated the versatility and control of thin-film microstructure offered by meniscusguided coating (MGC) techniques. Here, we analyze the qualitative and quantitative aspects of solution shearing, a model MGC method, using coating blades augmented with arrays of pillars. The pillars induce local regions of high strain rates—both shear and extensional—not otherwise possible with unmodified blades, and we use fluid mechanical simulations to model and study a variety of pillar spacings and densities. We then perform a statistical analysis of 130 simulation variables to find correlations with three dependent variables of interest: thin-film degree of crystallinity and transistor field-effect mobilities for charge-transport parallel (Μ para ) and perpendicular (Μ perp ) to the coating direction. Our study suggests that simple fluid mechanical models can reproduce substantive correlations between the induced fluid flow and important performance metrics, providing a methodology for optimizing blade design.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optimized Quantum Phase Estimation for Simulating Electronic States in Various Energy Regimes

While quantum algorithms for simulation exhibit better asymptotic scaling than their classical counterparts, they currently cannot be implemented on real-world devices. Instead, chemists and computer scientists rely on costly classical simulations of these quantum algorithms. In particular, the quantum phase estimation (QPE) algorithm is among several approaches that have attracted much attention in recent years for its genuine quantum character. However, it is memory-intensive to simulate and intractable for moderate system sizes. Herein, this paper discusses the performance and applicability of QPESIM, a new simulation of the QPE algorithm designed to take advantage of modest computational resources. In particular, we demonstrate the versatility of QPESIM in simulating various electronic states by examining the ground and core-level states of H 2 O. For these states, we also discuss the effect of the active-space size on the quality of the calculated energies. For the high-energy core-level states, we demonstrate that new QPE simulations for active spaces defined by 15 active orbitals significantly reduce the errors in core-level excitation energies compared to earlier QPE simulations using smaller active spaces.

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

Ab Initio Molecular Dynamics Simulations of Amorphous Calcium Carbonate: Interpretation of Pair Distribution Function and X-Ray Absorption Spectroscopy Data

The structure and transformation of hydrous amorphous calcium carbonate (ACC) are key to understanding biomineralization pathways and their relationship with the properties of the resulting material. Quantitative interpretation of scattering experiments aimed at elucidating the structure of ACC is challenging, due to the amorphous nature of this material, and, therefore, requires models for the structure and the scattering physics. Here, we generate physically realistic ensembles of hydrated ACC structures and their vibrational disorder from ab initio molecular dynamics (AIMD) simulations with an emphasis on enabling the consistent interpretation of the finer details of three complementary structural probes: neutron and x-ray pair distribution experiments and x-ray absorption spectroscopy (XAS). In each case, we simulate the signal directly in reciprocal space and then manipulate it into the real-space pair distribution function (PDF) or spectrum using the same procedures for the experimental and theoretical data. Good agreement with experiment was obtained across the three techniques with the simulations accounting well for all features in the spectra. Remaining small discrepancies pointed to differences between real samples and the idealized simulated systems such as deviations from the nominal CaCO3·nH2O stoichiometry. Additionally, the simulations offered a more accurate description of the local coordination environment of calcium than previous shell-by-shell fits to spectra of synthetic ACC and classical molecular dynamics simulations. This work demonstrates that AIMD is a powerful approach for extracting detailed structural information from neutron PDF, x-ray PDF, and XAS of amorphous carbonate phases.

Prange, Micah P.↗
DOE OSTI2026

Optimizing fluvial flood mitigation strategies: A multi-objective approach for cost-effective and socially-aware infrastructure feasibility analysis

Effective levee planning must balance capital cost, risk reduction, and community priorities. These objectives are rarely optimized together. This study presents a feasibility phase, simulationin-the-loop framework that couples terrain-based flood modeling with a socially aware multiobjective optimizer. Flood risk is measured as Expected Annual Exposed Population (EAEP), obtained by integrating exposure over Annual Exceedance Probability (AEP) nodes, mirroring the Hydrologic Engineering Center's Flood Damage Reduction Analysis (HEC-FDA) expected-annual formulation but with people rather than dollars. Exposure per scenario is computed by overlaying binary inundation masks with a population surface at the tract level. Distributional fairness is encoded through a Group Benefit Share (GBS) constraint that requires high-SVI tracts to receive at least a baseline share of annualized benefits. Capital cost is represented by a height-dependent unit-cost model suitable for screening. This study addresses the two-objective problem, minimize cost and expected annual exposure subject to the GBS constraint, using Non-Dominated Sorting Genetic Algorithm II (NSGA-II) and leveraging Pareto front for feasibility phase decision making. Implemented with terrain-based flood modeling, GeoFlood, for rapid scenario evaluation, the framework is demonstrated in Southeast Texas. The results reveal clear trade-offs among cost, risk, and social benefits and identify non-dominated levee height configurations that satisfy the benefit-share floor. The contributions are a scalable decision support method that operationalizes expected annual population-based risk, embeds enforceable benefit-sharing guarantees, and uses lightweight simulation to explore large design spaces before higher fidelity design stages.

Flood mitigation↗
DOE OSTI2024

Lithiophilic CoF 2 @C hollow spheres towards spatial lithium deposition for stable lithium metal batteries

Lithium metal (LM) is a promising anode for next-generation batteries due to its high theoretical capacity and low electrode potential. Nonetheless, side reactions, volume change, and unwanted lithium dendrite growth seriously limit the practical application of LM. Herein, with the aid of a hard template approach, a novel lithiophilic CoF 2 -carbon hollow sphere (CoF 2 @C-HS) composite material is successfully prepared via a facile in-situ fluorination and etching strategy. The lithiophilic CoF 2 acts as nucleation sites to reduce nucleation overpotential as well as induces the spatial Li deposition and the formation of LiF-rich solid electrolyte interphase (SEI), and the hollow carbon matrix can enhance the electrical conductivity and offer free space for LM deposition. Theoretical simulations reveal that the synergistic effect of lithiophilic CoF 2 and hollow carbon matrix homogenizes the electric field distribution and Li + flux. Benefiting from these advantages, the CoF 2 @C-HS-modified copper substrate electrode delivers an enhanced Coulombic efficiency (CE) of 93.7% for 280 cycles at 1 mA cm –2 and 1 mA h cm –2 . The symmetrical cell using CoF 2 @C-HS can stably cycle more than 1800 h with a low voltage hysteresis of 11 mV at a current density of 0.5 mA cm –2 and an areal capacity of 0.5 mA h cm –2 . Moreover, the Li@CoF 2 @C-HS composite anode enables more than 300 stable cycles at 1 C with a capacity retention of 95% in LiFePO 4 -based full cell and 110 stable cycles at 1 C in LiNi 0.8 Co 0.1 Mn 0.1 O 2 (NCM811)-based high-voltage full cell. Finally, this work might shed a new light on designing lithiophilic hosts to spatially confine LM deposition, realizing dendrite-free LM anodes and the practical applications of LM batteries.

25 ENERGY STORAGE↗
DOE OSTI2026

A robust alloy design (RAD) strategy for next-generation (IV) nuclear fission reactors

Next-generation nuclear reactors demand structural materials capable of withstanding extreme conditions, including high temperatures, intense neutron flux, and corrosive environments. Multi-Principal Element Alloys (MPEAs) have emerged as promising candidates due to their exceptional radiation tolerance, thermal stability, and compositional flexibility. This study introduces a versatile and customizable Robust Alloy Design (RAD) strategy for systematically designing MPEAs for GEN-IV reactor fuel cladding. The RAD framework integrates nuclear-relevant selection criteria, empirical parameter assessments, and high-throughput CALPHAD simulations to efficiently narrow compositional space and identify stable alloys. A unified RAD score developed for the first time, combines key performance metrics, including fuel-clad chemical interaction (FCCI), neutron absorption cross-section (NAC), valence electron configuration (VEC), and melting point factor (MPF), into a flexible ranking system adaptable to reactor-specific priorities. Among 724 candidates, V555(5Al–5Cr–5Fe–85V) emerged as the top alloy, validated experimentally with a homogeneous single-phase BCC microstructure and superior mechanical properties (nano-indentation: 3.389 ± 0.258 GPa; Vickers hardness: 240 ± 6.7 HV), significantly outperforming Zircaloy-4 and V-4Cr-4Ti. Importantly, the RAD strategy is not limited to nuclear applications; its customizable weighting system enables scalability to other extreme environments. This adaptability positions RAD strategy as a versatile tool for advanced materials design across multiple industries.

Alloy design↗
DOE OSTI2024

Transition-Metal-Related Quantum Emitters in Wurtzite AlN and GaN

Transition-metal centers exhibit a paramagnetic ground state in wide-bandgap semiconductors and are promising for nanophotonics and quantum information processing. Specifically, there is a growing interest in discovering prominent paramagnetic spin defects that can be manipulated using optical methods. Here, we investigate the electronic structure and magneto-optical properties of Cr and Mn substitutional centers in wurtzite AlN and GaN. We use state-of-the-art hybrid density functional theory calculations to determine level structure, stability, optical signatures, and magnetic properties of these centers. The excitation energies are calculated using the constrained occupation approach and rigorously verified with the complete active space configuration interaction approach. Our simulations of the photoluminescence spectra indicate that $Cr$$^{1+}_{Al}$ in AlN and $Cr$$^{1+}_{Ga}$ in GaN are responsible for the observed narrow quantum emission near 1.2 eV. We compute the zero-field splitting (ZFS) parameters and outline an optical spin polarization protocol for $Cr$$^{1+}_{Al}$ and $Cr$$^{1+}_{Ga}$. Our results demonstrate that these centers are promising candidates for spin qubits.

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

Modeling the Dynamics of Energetic Protons in Earth's Inner Magnetosphere

Recently, Van Allen Probes data have revealed the interesting energy-dependent dynamics of energetic protons in Earth's inner magnetosphere. To quantify the role of radial diffusion and charge exchange to the observed energy-dependent dynamics of ring current protons, a 1D radial diffusion model with charge exchange loss is implemented. The observed proton flux over the long-term period of November 2012 to September 2013 is first converted to phase space densities, which are then simulated by our model driven by an outer boundary condition at L* = 5.5 derived from the data. The simulation results show that our model generally captures the transport and acceleration of energetic protons at μ = 30, 50, and 80 MeV/G and K = 0.11 G 1/2 R E , suggesting that radial diffusion is the dominant source mechanism for >75 keV protons at L* = 3.5–5.5. Additionally, the observed fast decay of protons at lower μ and slow decay at higher μ are well captured by the model, demonstrating the dominant role of charge exchange in explaining the observed energy-dependent proton decay. For higher μ protons, prompt losses of protons on the time scale of hours are observed over a wide range of L*, which is too fast to be explained by charge exchange. Some of the prompt losses at high L* regions are reproduced by the model with outward radial diffusion to the outer boundary. However, many prompt losses observed at lower L* regions are not captured by the model, which could be due to other loss mechanisms including electromagnetic ion cyclotron wave scattering and field line curvature scattering.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Contrasting Responses of Hailstorms to Anthropogenic Climate Change in Different Synoptic Weather Systems

Hailstones and extreme precipitation generate substantial economic losses across the United States (US) and the globe. Their strong association with short-lived, intense convective storms poses a great challenge in predicting their future changes. Here, we conducted model simulations at 1.2 km grid spacing for severe convective storms with large hail and heavy precipitation that occurred in two typical types of synoptic-scale environments in spring seasons over the central US under both current and future climate conditions. We find that the responses of large hail (diameters >2.5 cm) to anthropogenic climate change (ACC) are markedly different between the hailstorms developed in the two types of synoptic-scale environments, with over 110% increase in large hail occurrences for the frontal systems, whereas less than 30% increase for the Great Plains low-level jet (GPLLJ) systems. This is explained by the larger increase in convective intensity and updraft width and a smaller increase in warm cloud depth in the frontal storms compared with the GPLLJ storms. Interestingly, the occurrences and intensity of heavy precipitation (rain rate >20 mm hr-1) in both types of systems are similarly sensitive to ACC (e.g., 40% and 33% increases in the occurrences for the frontal and GPLLJ systems, respectively). These results advance our knowledge of hail projection and have important implications for managing risks for future hail.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Uncertainty quantification for Bayesian active learning in rupture life prediction of ferritic steels

Abstract Three probabilistic methodologies are developed for predicting the long-term creep rupture life of 9–12 wt%Cr ferritic-martensitic steels using their chemical and processing parameters. The framework developed in this research strives to simultaneously make efficient inference along with associated risk, i.e., the uncertainty of estimation. The study highlights the limitations of applying probabilistic machine learning to model creep life and provides suggestions as to how this might be alleviated to make an efficient and accurate model with the evaluation of epistemic uncertainty of each prediction. Based on extensive experimentation, Gaussian Process Regression yielded more accurate inference ( $$Pearson\;correlation\;coefficent> 0.95$$ P e a r s o n c o r r e l a t i o n c o e f f i c e n t > 0.95 for the holdout test set) in addition to meaningful uncertainty estimate (i.e., coverage ranges from 94 to 98% for the test set) as compared to quantile regression and natural gradient boosting algorithm. Furthermore, the possibility of an active learning framework to iteratively explore the material space intelligently was demonstrated by simulating the experimental data collection process. This framework can be subsequently deployed to improve model performance or to explore new alloy domains with minimal experimental effort.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Symmetry breaking as predicted by a phase space Hamiltonian with a spin Coriolis potential

Here, we perform electronic structure calculations for a set of molecules with degenerate spin-dependent ground states ( 3 CH 2 , 2 CH$^{•}_{3}$, 3 O 2 ) going beyond the Born–Oppenheimer approximation and accounting for nuclear motion. According to a phase space approach that parameterizes electronic states (|Φ⟩) and electronic energies (E) by nuclear position and momentum [i.e., |Φ(R, P)⟩ and E(R, P)], we find that the presence of degenerate spin degrees of freedom leads to broken symmetry ground states. More precisely, rather than a single degenerate minimum at (R, P) = (R min , 0), the ground state energy has two minima at (R,P)=(R' min ,±P min ) (where R' min is close to R min ), dramatically contradicting the notion that the total energy of the system can be written in separable form as E = $\frac{P^2}{2M}$ + V el . Although we find that the broken symmetry solutions have small barriers between them for the small molecules, we hypothesize that the barriers should be macroscopically large for metallic solids, thus offering up a new phase-space potential energy surface for simulating the Einstein–de Haas effect.

Berry connection↗
DOE OSTI2025

Spatiotemporal Energy Cascade in Three-Dimensional Magnetohydrodynamic Turbulence

We present a new scale decomposition method to investigate turbulence in wavenumber-frequency space. Using 3D magnetohydrodynamic turbulence simulations, we show that magnetic fluctuations with time scales longer than the nonlinear time exhibit an inverse cascade toward even smaller frequencies. Low frequency magnetic fluctuations support turbulence, acting as an energy reservoir that is converted into plasma kinetic energy, the latter cascading toward large wavenumbers and frequencies, where it is dissipated. Our results shed new light on the spatiotemporal properties of turbulence, potentially explaining the origin and role of low frequency turbulent fluctuations in the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗