Search NASA⌕ Search

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.

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

At least 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

An Improved Plastically Dilatant Unified Viscoplastic Constitutive Formulation for Multiscale Analysis of Polymer Matrix Composites Under High Strain Rate Loading

Polymer matrix composites are commonly used to fabricate energy-absorbing structures expected to experience impact loading. As such, a detailed understanding of the dynamic response of the constituent materials is necessary. Since the rate, temperature, and pressure dependence of carbon fiber reinforced polymer matrix composites are primarily manifestations of the rate, temperature, and pressure dependence of the polymer matrix, it is crucial that the constitutive behavior of the matrix be accurately characterized. In this work, an existing unified viscoplastic constitutive formulation is extended to ensure thermodynamic consistency and to more accurately account for the tension-compression asymmetry observed in the response of polymeric materials. A new plastic potential function is proposed, and elementary loading conditions are utilized to determine relations between model constants to ensure nonnegative plastic dissipation, a necessary thermodynamic requirement. Expressions for plastic Poisson’s ratios are derived and are bounded by enforcing nonnegative plastic dissipation. The model is calibrated against available experimental data from tests conducted over a range of strain rates, temperatures, and loading cases on a representative thermoset epoxy; good correlation between simulations and experimental data is obtained. Temperature rises due to the conversion of plastic work to heat are computed via the adiabatic heat energy equation. The viscoplastic polymer model is then used as a constitutive model in the generalized method of cells micromechanics theory to investigate the effects of matrix adiabatic heating on the high strain rate response of a unidirectional composite. The thermodynamic consistency of the model ensures plastic dissipation can only cause an increase in temperature. Simulation results indicate that significant thermal softening due to the conversion of plastic work to heat is observed in the composite for matrix dominated deformation modes.

Impact loading↗
DOE OSTI2025

Atomic-scale understanding of oxide growth and dissolution kinetics of Ni-Cr alloys

Aqueous corrosion of metals is governed by formation and dissolution of a passivating, multi-component surface oxide. Unfortunately, a detailed atomistic description is challenging due to the compositional complexity and the need to consider multiple kinetic factors simultaneously. To this end, we combine experiments with a first-principles-derived, multiscale computational framework that transcends thermodynamic descriptions to explicitly simulate the kinetic evolution of surface oxides of Ni-Cr alloys as a function of composition, temperature, pH, and applied voltage. In the absence of pitting, we identify three distinct voltage regimes, which are kinetically dominated by oxide growth, dissolution, and competitive dissolution and reprecipitation. Evolving compositional gradients and oxide thickness are revealed, including a transition between a metastable Ni-Cr mixed oxide and a thick, porous Ni-dominated oxide. Beyond elucidating the underlying physics, we highlight the need for competing kinetics in models to properly predict the transition from passivation to corrosion. Our results provide a key step towards co-design of alloy composition alongside environmental conditions for sustainable use across a variety of critical energy and infrastructure applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

pH-Dependent Vibrational Dynamics Drives Excited-State Quenching in the Phycobiliprotein Complex PC645

Phycocyanin 645 (PC645) is a closed-form lightharvesting complex found in the lumen of the photosynthetic membrane of cryptophyte algae. These peripheral antenna complexes contain bilin chromophores that absorb sunlight and transfer excitation energy to the core antenna complexes embedded in the thylakoid membrane. The location of cryptophyte antenna complex on the luminal side of the membrane is unusual. During photosynthetic activity, the pH of the lumen drops, by up to two pH units. There is little known about how this pH-change affects the light-harvesting complexes. In this study, we report multiscale simulations using a computationally efficient density functional tight-binding framework to investigate the spectroscopy and excitation energy transfer in the PC645 complex. Complementary experiments were conducted using both steady-state and time-resolved spectroscopic measurements at low, neutral, and high pH values. Our study shows that (de)protonation of specific bilin pigments, namely, the mesobiliverdins (MBVs), modulates the excitation energies, excitonic couplings, and spectral densities. These changes cause excitation transfer rates to increase by up to a factor of two to three, leading to pH-dependent energy transfer pathways in the complex. Using this model, we calculated the pH-dependent fluorescence quantum yield of the system, obtaining quantitative agreement with the experimental results. These computational simulations, supported by experiments, identify MBVs as a more prominent excitation sink than previously realized, and that this role is tuned by pH.

Maity, Sayan [Constructor Univ., Bremen (Germany);↗
DOE OSTI2025

Spin correlations in the nematic quantum disordered state of FeSe

The quantum-disordered state in FeSe, intertwined with superconductivity and nematicity, has been a research focus in iron-based superconductors. However, the intrinsic spin excitations across the entire Brillouin zone in detwinned FeSe, crucial for understanding its magnetism and superconductivity, have remained unresolved. Using inelastic neutron scattering, we reveal that stripe spin excitations (Q = (1, 0)/(0, 1)) exhibit the C 2 symmetry, while Néel spin excitations (Q = (1, 1)) retain C 4 symmetry within the nematic state. Temperature-dependent differences between Q = (1, 0) and (0, 1) spin excitations above the structural transition unambiguously reveals the nematic quantum disordered state. Comparison with NaFeAs suggests the Néel excitations originate from enhanced 3d xy orbital correlations. Modeling the stripe dispersions using a J 1 -K-J 2 Heisenberg Hamiltonian, we establish a spin-interaction phase diagram, positioning FeSe near a crossover regime between the antiferroquadrupolar, Néel, and stripe orders. Our results provide key insights into the microscopic spin interactions and their role in the intertwined orders in iron-based superconductors.

Electronic properties and materials↗
NASA NTRS2021

Micromechanics Modeling of Textiles for Re-Entry Parachute Applications

Recent flight test projects and NASA missions have highlighted the challenges associated with accurately and efficiently modeling the behavior of parachute deployment systems needed for parachute design. Moreover, parachute deployment has been identified as one of the higher risk components for such missions. The analysis of textile fabrics used for atmospheric entry is inherently complex due to the multiple scales present in the fabric structure, including individual fiber filaments at the microscale, yarn bundles of fibers at the mesoscale, and the overall woven fabric at the macroscale. Computational tools for simulating fabric behavior must be able to account for the different mechanisms present at each scale without sacrificing computational efficiency. This work examines the generalized multiscale method of cells micromechanics theory, which has previously been used for the analysis of reinforced composite structures, to unreinforced textile fabrics. Modifications to the existing composite multiscale framework, implemented in NASA’s Multiscale Analysis Tool (NASMAT), include the specific mechanics unique to unreinforced textile fabrics, and overcoming the assumptions of a fixed fiber angle. It looks to assess the feasibility of using the NASMAT tool for efficient prediction of the response of unreinforced fabrics to loading such that it can ultimately be applied to fluid structure interaction tools for the prediction of parachute deployment systems. In this work, fabric behavior is simulated in NASMAT through homogenization of a triply periodic repeating unit cell, where the geometry of the subcells can change as a function of loading to represent the relative rotation and uncrimping that can occur in fabric tows. Predictions from the amended NASMAT code are compared to experimental data for uniaxial and off-axis tension to verify the ability of the code to incorporate lower-scale mechanics in prediction of unreinforced fabrics under loading.

Micromechanics↗
NASA NTRS2012

Expected Navigation Flight Performance for the Magnetospheric Multiscale (MMS) Mission

The Magnetospheric Multiscale (MMS) mission consists of four formation-flying spacecraft placed in highly eccentric elliptical orbits about the Earth. The primary scientific mission objective is to study magnetic reconnection within the Earth s magnetosphere. The baseline navigation concept is the independent estimation of each spacecraft state using GPS pseudorange measurements (referenced to an onboard Ultra Stable Oscillator) and accelerometer measurements during maneuvers. State estimation for the MMS spacecraft is performed onboard each vehicle using the Goddard Enhanced Onboard Navigation System, which is embedded in the Navigator GPS receiver. This paper describes the latest efforts to characterize expected navigation flight performance using upgraded simulation models derived from recent analyses.

Olson, Corwin↗
DOE OSTI2025

Integrated Hydro-terrestrial Modeling 2.0: Progress and Path Forward on Building a National Capability

It is the role of the U.S. federal government and its supporting agencies, including academia and future scientists, to ensure that its people have sustained and equitable freshwater services, as well as the critical knowledge necessary to make decisions about the future as it relates to freshwater services. Clear and consistent information and guidance from federal agencies is critical. Integrated Hydro-Terrestrial Modeling (IHTM), as a United States (U.S.) national capability, focuses on understanding, quantifying, and managing the replenishment of water supply through hydrologic cycle processes and their governing forces. To provide that information, we need enhanced IHTM capabilities that capitalize on the strengths of each U.S. governmental agency and its core mission. The first IHTM workshop was held in 2019, and its subsequent report was published in 2020. The U.S. Global Change Research Program (USGCRP) and member agencies held a second IHTM workshop (IHTM 2.0) from October 31 to November 2, 2023 in Reston, Virginia. The IHTM 2.0 workshop focused on the need to support a multiscale framework to accelerate research insights, better integrate operational and planning perspectives, and bridge national-to-regional capabilities to address major interdependent societal water challenges. The workshop was organized according to a “WHAT” and “HOW” framework, with the common underlying “WHY” being the integrated water resource challenges and the “WHO” defined through interagency and cooperating academic partners. The report provides a summary of plenary presentations and breakout discussions, and a road map that focuses on near-term activities.

99 GENERAL AND MISCELLANEOUS↗
NASA NTRS2009

Multiscale Failure Analysis of Laminated Composite Panels Subjected to Blast Loading Using FEAMAC/Explicit

This preliminary report demonstrates the capabilities of the recently developed software implementation that links the Generalized Method of Cells to explicit finite element analysis by extending a previous development which tied the generalized method of cells to implicit finite elements. The multiscale framework, which uses explicit finite elements at the global-scale and the generalized method of cells at the microscale is detailed. This implementation is suitable for both dynamic mechanics problems and static problems exhibiting drastic and sudden changes in material properties, which often encounter convergence issues with commercial implicit solvers. Progressive failure analysis of stiffened and un-stiffened fiber-reinforced laminates subjected to normal blast pressure loads was performed and is used to demonstrate the capabilities of this framework. The focus of this report is to document the development of the software implementation; thus, no comparison between the results of the models and experimental data is drawn. However, the validity of the results are assessed qualitatively through the observation of failure paths, stress contours, and the distribution of system energies.

Pineda, Evan J.↗
DOE OSTI2024

Anomaly Detection in Materials Digital Twins with Multiscale ICME for Additive Manufacturing

Detecting anomaly in fatigue and fracture experimental materials science is an interesting yet challenging topic. The reasons are threefold. First, the anomalous microstructure feature that gives rise to structural failure is small, sometimes in the order of 10 -7 of the interrogated volume. This, in turn, results in a highly imbalanced classification problem in machine learning (ML). Second, the consequence is high, in the sense that the test specimen is destructed in such case. Third, the convolution between microstructure stochasticity and the small probability of void nucleation, growth, and coalescence makes failure and fracture a hard-to-predict and challenging problem in materials science due to its irreproducibility, even experimentally. In this paper, we developed a materials digital twin and applied anomaly detection methods to detect voids and anomaly in additive manufacturing (AM). The materials digital twin is driven by two integrated computational materials engineering (ICME) models, which are kinetic Monte Carlo (kMC) and crystal plasticity finite element method (CPFEM). In conclusion, we demonstrated that by using anomaly detection, it is possible to detect voids and other defects in materials digital twin, which paves way for future research in integrating materials digital twin with its physical counterpart.

ICME↗
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.↗