SEARCH · Search NASA

Results for “Multiscale Modeling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41

All sources NASA NTRS DOE OSTI

DOE OSTI2025

From Electronic Structure to Ion Transport: Photoelectron Spectroscopy and Molecular Dynamics Simulations Reveal the Role of Anions in Lithium Battery Electrolytes

Electrolyte anions are pivotal for lithium battery performance, yet their fundamental electronic structural properties are not well understood. In this work, we employ a combination of negative-ion photoelectron spectroscopy (NIPES), ab initio calculations, and molecular dynamics (MD) simulations to investigate the electronic structures of three representative electrolyte anions. This multiscale approach enables us to elucidate how their intrinsic electronic properties govern anion–solvent interactions in gas-phase clusters, as well as lithium-ion (Li + ) solvation structures and ion transport behavior in the condensed phase. NIPES reveals that difluoro(oxalato)borate (DFOB – ), bis(fluorosulfonyl)imide (FSI – ), and bis(oxalato)borate (BOB – ) all exhibit high electron binding energies, with vertical/adiabatic detachment energies increasing from DFOB – (6.09/5.70 eV) to FSI – (6.80/6.10 eV) to BOB – (6.82/6.40 eV), correlating with enhanced oxidation stability. Ab initio calculations reveal that DFOB – /FSI – –solvent complexes bind Li + ∼ 10 kcal/mol stronger than BOB – series, aligning with the strength of a Li + –anion model. DFOB – exhibits pronounced charge localization on both oxygen and fluorine atoms, enabling their involvement in Li + coordination. In contrast, fluorine atoms in FSI – are largely electron-depleted and remain excluded from direct Li + binding. MD simulations further demonstrate that LiDFOB and LiFSI systems exhibit Li + diffusion coefficients three and five times higher than those of LiBOB across four common solvents. Notably, LiFSI salt in acetonitrile (AN) exhibits the fastest Li + diffusion among 12 electrolyte systems, highlighting the synergistic effect of FSI – and AN in promoting ion mobility. In conclusion, these findings provide a molecular-level understanding of the critical roles of anion and its microsolvation in optimizing Li + diffusion dynamics, once again emphasizing the positioning of FSI – and DFOB – as prime candidates for next-generation electrolytes.

25 ENERGY STORAGE↗

DOE OSTI2024

The hierarchical structure of organic mixed ionic–electronic conductors and its evolution in water

Polymeric organic mixed ionic–electronic conductors underpin several technologies in which their electrochemical properties are desirable. These properties, however, depend on the microstructure that develops in their aqueous operational environment. We investigated the structure of a model organic mixed ionic–electronic conductor across multiple length scales using cryogenic four-dimensional scanning transmission electron microscopy in both its dry and hydrated states. Four-dimensional scanning transmission electron microscopy allows us to identify the prevalent defects in the polymer crystalline regions and to analyse the liquid crystalline nature of the polymer. The orientation maps of the dry and hydrated polymers show that swelling-induced disorder is mostly localized in discrete regions, thereby largely preserving the liquid crystalline order. Therefore, the liquid crystalline mesostructure makes electronic transport robust to electrolyte ingress. Furthermore, this study demonstrates that cryogenic four-dimensional scanning transmission electron microscopy provides multiscale structural insights into complex, hierarchical structures such as polymeric organic mixed ionic–electronic conductors, even in their hydrated operating state.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Trilinos: Enabling Scientific Computing across Diverse Hardware Architectures at Scale

Trilinos is a community-developed, open-source software framework that facilitates building large-scale, complex, multiscale, multiphysics simulation code bases for scientific and engineering problems. Since the Trilinos framework has undergone substantial changes to support new applications and new hardware architectures, this document is an update to “An Overview of the Trilinos project” by Heroux et al. (ACM Transactions on Mathematical Software, 31(3):397–423, 2005). It describes the design of Trilinos, introduces its new organization in product areas, and highlights established and new features available in Trilinos. Particular focus is put on the modernized software stack based on the Kokkos ecosystem to deliver performance portability across heterogeneous hardware architectures. This article also outlines the organization of the Trilinos community and the contribution model to help onboard interested users and contributors.

Heterogeneous Hardware Architectures↗

NASA NTRS2015

The Parameterization of Top-Hat Particle Sensors with Microchannel-Plate-Based Detection Systems and its Application to the Fast Plasma Investigation on NASA's Magnetospheric MultiScale Mission

The most common instrument for low energy plasmas consists of a top-hat electrostatic analyzer geometry coupled with a microchannel-plate (MCP)-based detection system. While the electrostatic optics for such sensors are readily simulated and parameterized during the laboratory calibration process, the detection system is often less well characterized. Furthermore, due to finite resources, for large sensor suites such as the Fast Plasma Investigation (FPI) on NASA's Magnetospheric Multiscale (MMS) mission, calibration data are increasingly sparse. Measurements must be interpolated and extrapolated to understand instrument behavior for untestable operating modes and yet sensor inter-calibration is critical to mission success. To characterize instruments from a minimal set of parameters we have developed the first comprehensive mathematical description of both sensor electrostatic optics and particle detection systems. We include effects of MCP efficiency, gain, scattering, capacitive crosstalk, and charge cloud spreading at the detector output. Our parameterization enables the interpolation and extrapolation of instrument response to all relevant particle energies, detector high voltage settings, and polar angles from a small set of calibration data. We apply this model to the 32 sensor heads in the Dual Electron Sensor (DES) and 32 sensor heads in the Dual Ion Sensor (DIS) instruments on the 4 MMS observatories and use least squares fitting of calibration data to extract all key instrument parameters. Parameters that will evolve in flight, namely MCP gain, will be determined daily through application of this model to specifically tailored in-flight calibration activities, providing a robust characterization of sensor suite performance throughout mission lifetime. Beyond FPI, our model provides a valuable framework for the simulation and evaluation of future detection system designs and can be used to maximize instrument understanding with minimal calibration resources.

Instrument-model↗

DOE OSTI2024

Stabilization of Charge-Transfer Excited States in Biological Systems: A Computational Focus on the Special Pair in Photosystem II Reaction Centers

Charge-transfer (CT) excited states play an important role in many biological processes. However, many computational approaches often inadequately address the equilibration effects of nuclear and environmental degrees of freedom on these states. One prominent example of systems in which CT states are of utmost importance is reaction centers (RC) in photosystems. Here we use a multiscale approach combined with time-dependent density functional theory to explore the lowest CT excited state of the special pair P D1 –P D2 in the Photosystem II-RC of a cyanobacterium. We find that the nonequilibrium CT excited state resides near the Soret band, making an exciton the lowest-energy excited state. However, accounting for nuclear and state-specific dielectric equilibration along the CT potential energy surface (PES), the CT state P D1 – –P D2 + stabilizes energetically below the excitonic state. Importantly, this underscores the crucial role of state-specific solvation in mapping the PES of CT states, as demonstrated in a simplified dimer model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2009

The Magnetospheric Multiscale Missions Fast Plasma Investigations Dual Electron Spectrometer Development

The Magnetospheric Multiscale mission (MMS) is designed to examine magnetic reconnection that occurs on both the Earths dayside magnetopause and in the magnetotail region on Earths night side. In order to resolve fine structures of the three dimensional electron distributions in both regions, the Fast Plasma Investigation's (FPI) Dual Electron Spectrometer (DES) is designed to measure electron distributions with a time resolution of 30 ms. In order to achieve this unprecedented sampling rate, the DES will have eight individual spectrometers each sampling 180 x 45 degree sections of the sky. Because of the field of view limitations of top hat analyzers, each spectrometer will use electro-static deflectors to change its look direction. The engineering model of the DES has been fabricated and tested. We will present the results of measurements for fields of view, angular FVVHM responses, dE/E, analyzer constant, and geometric factors for all deflection states. We will compare these results to simulation results and discuss causes of the response variations.

Shappirio, M.↗

NASA NTRS2016

Stable Reconnection at the Dusk Flank Magnetopause

The dusk flank magnetopause was surveyed with instruments on board the Magnetospheric Multiscale (MMS) spacecraft on 28 August 2015 between 13:55 UT and 14:15 UT during a period of persistent southward interplanetary magnetic field (IMF) with varying dawn-dusk component. Plasma measurements (500 eV electrons are greater than 2 keV ions) revealed the existence of at least one active reconnection region that persisted throughout the interval. The reconnection region convected equatorward despite the poleward and tailward magnetosheath flow, which ranged from slightly sub-Alfvenic to slightly super-Alfvenic throughout the interval. These results suggest that magnetic reconnection moved in response to changes in the IMF clock angle rather than the magnetosheath flow, which is corroborated using predictions of the maximum magnetic shear model.

Gomez, R. G.↗

NASA NTRS2006

Multiscale Analysis of Delamination of Carbon Fiber-Epoxy Laminates with Carbon Nanotubes

A multi-scale analysis is presented to parametrically describe the Mode I delamination of a carbon fiber/epoxy laminate. In the midplane of the laminate, carbon nanotubes are included for the purposes of selectively enhancing the fracture toughness of the laminate. To analyze carbon fiber epoxy carbon nanotube laminate, the multi-scale methodology presented here links a series of parameterizations taken at various length scales ranging from the atomistic through the micromechanical to the structural level. At the atomistic scale molecular dynamics simulations are performed in conjunction with an equivalent continuum approach to develop constitutive properties for representative volume elements of the molecular structure of components of the laminate. The molecular-level constitutive results are then used in the Mori-Tanaka micromechanics to develop bulk properties for the epoxy-carbon nanotube matrix system. In order to demonstrate a possible application of this multi-scale methodology, a double cantilever beam specimen is modeled. An existing analysis is employed which uses discrete springs to model the fiber bridging affect during delamination propagation. In the absence of empirical data or a damage mechanics model describing the effect of CNTs on fracture toughness, several tractions laws are postulated, linking CNT volume fraction to fiber bridging in a DCB specimen. Results from this demonstration are presented in terms of DCB specimen load-displacement responses.

Riddick, Jaret C.↗

NASA NTRS2016

Energy Limits of Electron Acceleration in the Plasma Sheet During Substorms: A Case Study with the Magnetospheric Multiscale (MMS) Mission

We present multipoint observations of earthward moving dipolarization fronts and energetic particle injections from NASAs Magnetospheric Multiscale mission with a focus on electron acceleration. From a case study during a substorm on 02 August 2015, we find that electrons are only accelerated over a finite energy range, from a lower energy threshold at approx. 7-9 keV up to an upper energy cutoff in the hundreds of keV range. At energies lower than the threshold energy, electron fluxes decrease, potentially due to precipitation by strong parallel electrostatic wavefields or initial sources in the lobes. Electrons at energies higher than the threshold are accelerated cumulatively by a series of impulsive magnetic dipolarization events. This case demonstrates how the upper energy cutoff increases, in this case from approx. 130 keV to >500 keV, with each depolarization/injection during sustained activity. We also present a simple model accounting for these energy limits that reveals that electron energization is dominated by betatron acceleration.

Turner, D. L.↗

NASA NTRS2021

The Dynamics of a High Mach Number Quasi-Perpendicular Shock: MMS Observations

Shock parameters at Earth’s bow shock, in rare instances, can approach the Mach numbers predicted at astrophysical shocks and supernova remnants. We present our analysis of a high Alfv ́en Mach number (MA= 27) shock, by utilizing multipoint measurements from the Magnetospheric Multiscale (MMS) spacecraft during a crossing of Earth’s quasi-perpendicular bow shock. We find that the shock dynamics are mostly driven by reflected ions, perturbations that they generate, and nonlinear amplification of the perturbations. Our analyses indicate that reflected ions create modest magnetic enhancements upstream of the shock front which evolve in a nonlinear manner as they traverse the shock foot. They can transform into proto-shocks that propagate at small angles to the magnetic field and towards the bow shock. The nonstationary bow shock shows signatures of both reformation and surface ripples. Our observations indicate that although shock reformation occurs, the main shock layer never disappears. These observations are at high plasmaβ, a parameter regime which has not been well explored by numerical models.

H Madanian↗

NASA NTRSDate not supplied

The Lunar GNSS Receiver Experiment (LuGRE)

The Lunar GNSS Receiver Experiment (LuGRE) is a joint NASA-Italian Space Agency (ASI) payload on the Firefly Blue Ghost Mission 1 (BGM1) with the goal to demonstrate GNSS-based positioning, navigation, and timing at the Moon. LuGRE was chosen by the NASA Commercial Lunar Payload Services (CLPS) program as one of ten payloads on its “19D” task order for delivery to the lunar surface in 2023. The LuGRE payload consists of a weak-signal GNSS receiver, a high-gain L-band patch antenna, a low-noise amplifier, and an RF filter. The receiver will track GPS L1 C/A and L5, and Galileo E1 and E5a signals and will return pseudorange, carrier phase, and Doppler measurements to the ground. It will also calculate least-squares point solutions and Kalman-filter based navigation solutions onboard. In addition, the receiver features the capability to record raw I/Q baseband samples for downlink and ground processing. LuGRE will build on the legacy of prior missions in the Space Service Volume (SSV) including the initial experiments by AMSAT-OSCAR 40 and others, the GOES-R series of geostationary weather satellites, and the NASA Magnetospheric Multiscale (MMS) mission currently operating on GPS-based navigation at nearly 50% of lunar distance. Further, LuGRE will be one of the very first demonstrations of GNSS signal reception and navigation in the lunar environment and on the lunar surface, paving the way for operational use by future lunar missions such as Orion, Gateway, robotic and human landers, and surface rovers. Ultimately, all LuGRE science data will be released to a public data archive for the benefit of the GNSS and space communities. This paper provides a detailed overview of the LuGRE payload, including its design, concept of operations, and its predicted ability to meet its core science objectives. The baseline science investigations and priorities are outlined. Simulated performance results are shown based on the latest calibrated models including signal strength, signal availability, onboard navigation performance and convergence properties, and ground-based post-processed navigation performance.

LuGRE↗

DOE OSTI2026

Multiphysics simulation of TRISO fuel compacts and the effects of homogenization on silver release predictions

This work studies the impact of explicit and homogenized modeling approaches on the multiphysics simulation of TRistructural ISOtropic (TRISO) fuel compacts in prismatic High Temperature Gas Reactors (HTGRs) and the silver release predictions. TRISO fuel compacts exhibit complex double heterogeneity that significantly affects heat conduction, neutron transport, and fission products release. In this work, we use Cardinal, a multiphysics tool based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, to couple neutron transport and heat conduction. OpenMC is used for neutron transport and the MOOSE heat transfer module is used for the heat conduction. Then, we use BISON for the silver release predictions. Two different modeling approaches—explicit modeling of individual TRISO particles and homogenized representation using effective thermal properties—are compared at high TRISO packing fractions (20% and 40%) across varying power densities. Results demonstrate that homogenization significantly under-predicts peak temperatures, for the case of high power/TRISO, there is a difference of 83.73 K in the maximum temperature. Additionally, homogenization underestimates the silver release predictions compared to explicit modeling. Specifically, the temperature and power density differences lead to significant differences in silver release predictions. This work demonstrates the importance of accurately modeling heterogeneity of the TRISO particles to reliably predict fission product release and assess reactor safety margins.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2015

High Frequency Design Considerations for the Large Detector Number and Small Form Factor Dual Electron Spectrometer of the Fast Plasma Investigation on NASA's Magnetospheric Multiscale Mission

Each half of the Dual Electron Spectrometer (DES) of the Fast Plasma Investigation (FPI) on NASA's Magnetospheric MultiScale (MMS) mission utilizes a microchannel plate Chevron stack feeding 16 separate detection channels each with a dedicated anode and amplifier/discriminator chip. The desire to detect events on a single channel with a temporal spacing of 100 ns and a fixed dead-time drove our decision to use an amplifier/discriminator with a very fast (GHz class) front end. Since the inherent frequency response of each pulse in the output of the DES microchannel plate system also has frequency components above a GHz, this produced a number of design constraints not normally expected in electronic systems operating at peak speeds of 10 MHz. Additional constraints are imposed by the geometry of the instrument requiring all 16 channels along with each anode and amplifier/discriminator to be packaged in a relatively small space. We developed an electrical model for board level interactions between the detector channels to allow us to design a board topology which gave us the best detection sensitivity and lowest channel to channel crosstalk. The amplifier/discriminator output was designed to prevent the outputs from one channel from producing triggers on the inputs of other channels. A number of Radio Frequency design techniques were then applied to prevent signals from other subsystems (e.g. the high voltage power supply, command and data handling board, and Ultraviolet stimulation for the MCP) from generating false events. These techniques enabled us to operate the board at its highest sensitivity when operated in isolation and at very high sensitivity when placed into the overall system.

Microchannel plates high frequency performance↗

DOE OSTI2024

Scalable Multiphysics Block Preconditioning for Low Mach Number Compressible Resistive MHD with Application to Magnetic Confinement Fusion

This study investigates multiphysics block preconditioners that are critical in devising scalable Newton–Krylov iterative solvers for longer time-scale fully implicit fluid plasma models. The specific model of interest is the visco-resistive, low Mach number, compressible magnetohydrodynamics (MHD) model. This model describes the dynamics of conducting fluids in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics. The specific application of interest that motivates this study is the macroscopic simulation of longer time-scale stability and disruptions of magnetic confinement fusion devices, specifically the ITER Tokamak. The computational solution of the governing balance equations for mass, momentum, heat transfer, and magnetic induction for resistive MHD systems can be extremely challenging. These difficulties arise from both the strong nonlinear, nonsymmetric coupling of fluid and electromagnetic phenomena as well as the significant range of time and length scales that the interactions of these physical mechanisms produce. To handle the range of time and spatial scales of interest, a fully implicit unstructured variational multiscale finite element formulation is employed. For the scalable solution of the Newton linearized systems, fully coupled block preconditioners are designed to leverage algebraic multigrid subsolves. In conclusion, results are presented for the strong and weak scaling of the method as well as the robustness of these techniques for a large range of Lundquist numbers.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Improving High Resolution Offshore Wind Resource Assessments and Forecasts using Observations in the MA/RI Lease Areas

The third Wind Forecast Improvement Project (WFIP3) sought to improve understanding of the physical phenomena in the atmosphere and ocean that dictate the structure and variability of wind and thermodynamic fields within the Marine Atmospheric Boundary Layer (MABL). WFIP3 conducted a comprehensive 18-month observational study over the Northeast U.S. outer continental shelf, a high use coastal zone, using a 3-D multiscale sensor array to highly resolve the temporal, vertical, and horizontal structure of the coupled atmospheric and oceanic boundary layers. Multiple land-based study sites adjacent to the coastal ocean observed surface meteorology and vertical profiles of atmospheric properties via passive infrared and microwave radiometers, active lidars and radars, and radiosondes. At sea, an array of surface flux buoys and two vertical profiling lidar buoys observed both atmospheric and oceanic properties, augmented by land-based oceanographic radar systems and routine ship-based surveys. Intensive observations of the MABL over the ocean was done from an air-sea interaction flux tower and extended deployments of a large barge platform. WFIP3 focused on mesoscale and sub-mesoscale flows -- including sea breezes, low-level jets, low-level clouds, and coastal storms -- and the ability of advanced numerical model parameterizations to represent them within fully coupled oceanic and atmospheric modeling systems and foundational weather forecast models. Numerous critical forecasting phenomena were observed, however the project was terminated prior to the completion of the field observational period and the analysis period.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Accelerating multiscale electronic stopping power predictions with time-dependent density functional theory and machine learning

Knowing the rate at which particle radiation releases energy in a material, the “stopping power,” is key to designing nuclear reactors, medical treatments, semiconductor and quantum materials, and many other technologies. While the nuclear contribution to stopping power, i.e., elastic scattering between atoms, is well understood in the literature, the route for gathering data on the electronic contribution has for decades remained costly and reliant on many simplifying assumptions, including that materials are isotropic. We establish a method that combines time-dependent density functional theory (TDDFT) and machine learning to reduce the time to assess new materials to hours on a supercomputer and provide valuable data on how atomic details influence electronic stopping. Our approach uses TDDFT to compute the electronic stopping from first principles in several directions and then machine learning to interpolate to other directions at a cost of 10 million times fewer core-hours. We demonstrate the combined approach in a study of proton irradiation in aluminum and employ it to predict how the depth of maximum energy deposition, the “Bragg Peak,” varies depending on the incident angle—a quantity otherwise inaccessible to modelers and far outside the scales of quantum mechanical simulations. The lack of any experimental information requirement makes our method applicable to most materials, and its speed makes it a prime candidate for enabling quantum-to-continuum models of radiation damage. The prospect of reusing valuable TDDFT data for training the model makes our approach appealing for applications in the age of materials data science.

36 MATERIALS SCIENCE↗

NASA NTRS2021

Temporal Interpolation of Geostationary Satellite Imagery With Optical Flow

Applications of satellite data in areas such as weather tracking and modeling, ecosystem monitoring, wildfire detection, and land-cover change are heavily dependent on the tradeoffs to spatial, spectral, and temporal resolutions of observations. In weather tracking, high-frequency temporal observations are critical and used to improve forecasts, study severe events, and extract atmospheric motion, among others. However, while the current generation of geostationary (GEO) satellites has hemispheric coverage at 10-15-min intervals, higher temporal frequency observations are ideal for studying mesoscale severe weather events. In this work, we present a novel application of deep learning-based optical flow to temporal upsampling of GEO satellite imagery. We apply this technique to 16 bands of the GOES-R/Advanced Baseline Imager mesoscale dataset to temporally enhance full-disk hemispheric snapshots of different spatial resolutions from 10 to 1 min. Experiments show the effectiveness of task-specific optical flow and multiscale blocks for interpolating high-frequency severe weather events relative to bilinear and global optical flow baselines. Finally, we demonstrate strong performance in capturing variability during convective precipitation events.

Image processing↗

NASA NTRSDate not supplied

Detecting Risk and Anomalies in Airplane Dynamics Through Entropic Analysis of Time Series Data

Despite recent efforts to move away from traditional threshold exceedance detection methods for aircraft state monitoring, modern aircraft still rely on safety thresholds to communicate to pilots the identification of an anomaly in the aircraft when a threshold is surpassed. Current anomaly detection methods mainly depend on uninterpretable machine learning models to learn complex patterns and relationships contained in the time series data of aircraft. Although these methods are capable of identifying known anomalies, their deficiency in interpretability presents a challenge when translating them to different aircraft. To overcome this deficiency, entropic analysis of aircraft dynamics seeks to characterize the complexity, or lack thereof, of the aircraft dynamics prior to the development of a risk scenario. This complexity characterization provides a more straightforward summary of state changes in the dynamics of flight variables. To build a foundation for entropic analysis, we analyzed the complexity of unstable approaches, an anomalous event present in many of today’s aviation accidents. The analysis revealed a statistically significant difference in the complexity distribution of flight variables under a stable approach versus an unstable approach. These differences in complexity were especially notable minutes before an approach was identified as unstable. Moreover, the multiscale entropic analysis revealed the presence of signal complexity at multiple time scales across multiple time windows before landing. By capturing state changes and corrections in the aircraft dynamics using entropy, advanced, yet still interpretable, sensor systems based on entropic frameworks from this study can be constructed in the future using classical machine learning approaches.

Risk detection↗