Search NASA⌕ Search

SEARCH · Search NASA

Results for “multiscale interactions”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2026

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗
DOE OSTI2025

Large-eddy simulation of an atmospheric bore and associated gravity wave effects on wind farm performance in the southern Great Plains

Gravity waves are a common occurrence in the atmosphere, with a variety of generation mechanisms. Their impact on wind farms has only recently gained attention, with most studies focused on wind farm-induced gravity waves. In this study, the interaction between a wind farm and gravity waves generated by an atmospheric bore event is assessed using multiscale large-eddy simulations. The atmospheric bore is created by a thunderstorm downdraft from a nocturnal mesoscale convective system (MCS). The associated gravity waves impact the wind resource and power production at a nearby wind farm during the American Wake Experiment (AWAKEN) in the US southern Great Plains. A two-domain nested setup (Δx=300 and 20 m) is used in the Weather Research and Forecasting (WRF) model, forced with data from the High-Resolution Rapid Refresh model, to capture both the formation of the bore and its interaction with individual wind turbines. The MCS is resolved on the large outer domain, where the structure of the bore and the associated gravity waves are found to be especially sensitive to parameterized microphysics processes. On the finer inner domain, gravity wave interactions with individual wind turbines are resolved; wake dynamics are captured using a generalized actuator disk parameterization in WRF. The gravity waves are found to have a strong effect on the atmosphere above the wind farm; however, the effect of the waves is more nuanced closer to the surface where there is additional turbulence, both ambient and wake-generated. Notably, the gravity waves modulate the mesoscale environment by weakening and dissipating the preexisting low-level jet, which reduces hub-height wind speed and hence the simulated power output, which is confirmed by the observed supervisory control and data acquisition (SCADA) power data. Additionally, the gravity waves induce local wind direction variations correlated with fluctuations in pressure, which lead to fluctuations in the simulated power output as various turbines within the farm are subjected to waking from nearby turbines.

17 WIND ENERGY↗
DOE OSTI2022

hPIC2: A hardware-accelerated, hybrid particle-in-cell code for dynamic plasma-material interactions

The exascale era of high performance computing promises to bring the field of computational plasma physics ever closer to the goal of accurate multiscale modeling. Such computers will rely on hardware acceleration to offload work to dedicated components, notably general-purpose graphics processing units (GPUs). However, devices from different manufacturers require software to be written with different parallel programming models, greatly increasing the code maintenance burden of applications designed to perform on more than one such device. hPIC2 is a hybrid plasma simulation code developed with the Kokkos performance portability framework to target the architectures that will drive exascale computing for the foreseeable future. As a hybrid simulation code, hPIC2 investigates the simultaneous use of various plasma models on the same domain, at the same time. hPIC2 also optionally couples to RustBCA, which accurately models ion-material interactions using the binary collision approximation (BCA) method. In conclusion, hPIC2 therefore achieves scalable performance on a variety of computing architectures when simulating complex and diverse plasmas, particularly near plasma-material interfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2004

The Role of Auroral Imaging in Understanding Ionosphere-Inner Magnetosphere Interactions

The more ways we probe the ionosphere and inner magnetosphere, the better we can understand their interaction. For example, the multifaceted imaging of geospace with the IMAGE mission complements the more traditional in situ measurements made with many previous missions. Together they have enabled new knowledge of the ionosphere-magnetosphere (IM) coupling. The role of imaging the aurora in understanding this interaction has received renewed attention recently. Based on in situ data, such as FAST or DMSP, and our recent theories, we believe that imaging multiscale features of the aurora is a key component to gaining insight into the processes and mechanisms at work. This talk will explore how auroral imaging can be used to provide improved insight of the dynamics of IM interaction on micro and meso scales, with an emphasis on the current limitations and future possibilities of quantitative analyses.

Spann, Jim↗
DOE OSTI2022

Multiscale spatial analysis of fracture arrangement and pattern reconstruction using Ripley's K-function

This work presents novel multiscale spatial data analytics using Ripley's K-function, as a measure of spatial interaction, to study one-dimensional arrangement of fractures. Fracture spatial arrangements are classified into clustered, anticlustered, or indistinguishable from random by testing statistical significance of the calculated Ripley's K-function. Characterizations of fracture arrangements are performed as a function of length scale and position. Analysis of the K-function along the study interval identifies where fracture clustering and anticlustering occur. Here, a simulation technique is also introduced here to statistically reconstruct spatial arrangements and to generate fracture realizations that are spatially similar to the fractures observed in the field. With this simulation technique, one can also fill spatial gaps in fracture measurements where data are absent, unreliable, or unused. Synthetic as well as field-measured 1D fracture datasets are used for testing and demonstration. Methods introduced in this work can be readily applied to fracture datasets observed in outcrops, borehole image logs, and cores.

58 GEOSCIENCES↗
DOE OSTI2022

The Critical Stack Pressure to Alter Void Generation at Li/Solid-Electrolyte Interfaces during Stripping

The lithium stripping process generates vacancies, which may accumulate as voids and lead to uneven current distribution and dendrite growth in the following plating cycles. A stack pressure is typically required during stripping, but how to optimize the stack pressure is not clear. In this work, extremely lithiophilic Li/Li 2 O and lithiophobic Li/LiF interfaces were used to reveal the combining effect of interface interaction and stack pressure induced lithium creep on the stripping critical current density (CCD). A multiscale simulation scheme with Density Functional Theory (DFT), kinetic Monte Carlo (KMC) simulations, and an analytical model was developed. The analytical model predicted lithiophobic interfaces require a higher stack pressure than lithiophilic interfaces to reach the same CCD. The KMC simulations also showed higher stack pressure is needed at lithiophobic interfaces to accelerate Li vacancy diffusion into the bulk and maintain a flat surface. This stack pressure needs to be high enough to alter the Li forward-and-backward hopping barriers at the interface. This multiscale simulation scheme illustrates the importance to include the chemical-mechanical effects during Li stripping morphology evolution. It can be used to design ideal interlayer coating materials to maintain a flat Li surface during cycling.

25 ENERGY STORAGE↗
DOE OSTI2025

Influence of simple terrain on the spatial variability of a low-level jet and wind farm performance in the AWAKEN field campaign

In wind energy research, scientific challenges are often associated with complex terrain sites, where orography, vegetation, and buildings disrupt flow uniformity. However, even sites characterized as simple terrain can exhibit significant spatial variability in wind speed, particularly during stable boundary layers (SBLs) and low-level jets (LLJs). This study investigates these terrain interactions using both simulations and observations from the American WAKe ExperimeNt (AWAKEN). We employ a multiscale Weather Research and Forecasting (WRF) model simulation, integrating mesoscale forcing in the coarse domains and representing three rows of turbines from the King Plains wind farm as generalized actuator disks (GAD) in the large-eddy simulation (LES) domains. During a nocturnal LLJ event on 3 April 2023, the downstream, wake-affected turbine rows outperformed the upstream, unwaked row by 25 %–51 %. This counterintuitive result arises from terrain-induced streamwise variations in hub-height wind speed of approximately 4 m s −1 over 5 km – equivalent to ∼50 % of the upstream reference speed. This enhancement outweighs the wake-induced reduction in mean wind speed (∼12 %) and global blockage effects reported in the literature (∼1 %–3.4 %). The multiscale simulations capture the intra-farm spatial variability in power performance observed in SCADA data. Terrain-induced vertical displacement of the LLJ, coupled with large wind shear below the jet maximum, drives the substantial streamwise acceleration within the wind farm. These findings underscore the importance of accounting for spatial variability related to terrain, even in simple landscapes, particularly during LLJ conditions. Incorporating such effects into reduced-order modeling frameworks for wind farm design and control could significantly enhance their effectiveness.

17 WIND ENERGY↗
DOE OSTI2023

Understanding grain boundary segregation under irradiation: A mesoscale modeling study

Grain boundaries (GBs) are an important component in the design of irradiation resistant materials as they play a crucial role as point defect sinks for non-equilibrium vacancies and self-interstitials. However, in the presence of sustained point defect fluxes, the coupled transport of atomic components with inherently different diffusivities often leads to radiation-induced segregation (RIS), with undesired effects on properties such as intergranular corrosion and stress-assisted cracking. In this talk, we investigate the fundamental mechanisms leading to GB segregation in BCC FeCr-based alloys using multiscale modeling and experiments. A phase-field model is employed to account for both thermodynamic and kinetic segregation mechanisms. The full Onsager transport matrix is derived using kinetic Monte Carlo, molecular dynamics, and kinetic cluster expansion (KineCluE) calculations. The thermodynamics of GB segregation is described via an atomic density-based modification to the CALPHAD free energy. The modeling predictions are then compared against segregation data from atom probe characterization for various temperatures, compositions, and neutron irradiation conditions. We will conclude by discussing the utility of the multiscale modeling approach to shed light on the fundamental mechanisms leading to solute-GB interaction.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multiscale Microtubule Dynamics in Active Nematics

In microtubule-based active nematics, motor-driven extensile motion of microtubule bundles powers chaotic large-scale dynamics. We quantify the interfilament sliding motion both in isolated bundles and in a dense active nematic. The extension speed of an isolated microtubule pair is comparable to the molecular motor stepping speed. In contrast, the net extension in dense 2D active nematics is significantly slower; the interfilament sliding speeds are widely distributed about the average and the filaments exhibit both contractile and extensile relative motion. These measurements highlight the challenge of connecting the extension rate of isolated bundles to the multimotor and multifilament interactions present in a dense 2D active nematic. Furthermore, they also provide quantitative data that is essential for building multiscale models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Automated Symbolic Upscaling: 2. Model Generation for Extended Applicability Regimes

Abstract In this second part of the two paper series, we detail an algorithmic procedure for systematically implementing the generalized closure form strategy presented in Part 1. This strategy extends the applicability of homogenized models with respect to classical homogenization theory, as demonstrated in Part 1 where upscaled models are rigorously derived in moderately reactive physical regimes. After encoding the algorithm into Symbolica, an automated upscaling framework, we upscale two reactive mass transport problems and numerically validate the resulting nonlinear homogenized models by showing the absolute error estimates predicted by homogenization theory are satisfied. In both problems, nontrivial closure forms and closure problems are automatically formulated using the encoded strategy with no human interaction, nor prior knowledge regarding the closure required for the systems. We hope these demonstrations spark further interest in automated analytical frameworks for multiscale modeling, as such capabilities are invaluable for generating rigorous multiscale models of complex phenomena in porous media.

Pietrzyk, Kyle↗
DOE OSTI2025

Solute clustering in polycrystals: Unveiling the interplay of grain boundary junction and long-range solute attraction effects

Spectral analysis of local atomic environments has become a powerful tool for studying solute atom segregation and interactions at grain boundaries in nanocrystalline alloys. When applied to individual grain boundaries, the spectral analysis has shown that solute-solute interaction can be either attractive or repulsive, with long-range relative attraction enhancing the likelihood for solute atoms to begin clustering. In this article, we combine this analysis with a new grain-boundary structure descriptor based on grain-boundary atom coordination, to investigate the impact of grain-boundary junctions on solute atom segregation in polycrystals. Specifically, we systematically characterize the tendency of solute clusters to begin forming at various types of ordinary grain boundaries, triple junctions, and high-order junctions in Ag polycrystals containing either Ni or Cu solute atoms. Our findings demonstrate that the formation of solute clusters at grain boundaries is primarily driven by long-range relative solute attraction, rather than short-range solute-solute interactions. Furthermore, this effect is most pronounced near grain-boundary junctions. Our study highlights the multiscale nature of solute segregation at crystalline interfaces and provides new insights into the complex phenomena governing heterogeneous solute segregation in grain-boundary networks.

Grain-boundary junction↗
DOE OSTI2020

Tracking ion intercalation into layered Ti 3 C 2 MXene films across length scales

Enhancing the energy stored and power delivered by layered materials relies strongly on improved understanding of the intricate interplay of electrolyte ions, solvents, and electrode interactions as well as the role of confinement. Here we report a highly integrated study with multiscale theory/modelling and experiments to track the intercalation of aqueous Li + , Na + , K + , Cs + , and Mg 2+ ions into Ti 3 C 2 MXene. The integrated analysis of experiments assisted by theory/modelling allows for a deep understanding of energy storage processes highlighting the importance of the dynamics of cations, their positionings between MXene sheets, their effects on mechanical properties and capacitive energy storage. Computational simulations and operando calorimetry measurements prove the processes involving cation dehydration and H+ rehydration, showing a good correlation for heat variations between experiments and theory. Operando liquid AFM mapped energy dissipation of ions appears non-uniformly across the MXene surface, indicating heterogeneities of ions inside the MXene and confirming partially the ion behaviour obtained in theory. We directly demonstrate that the average distance between the cation and MXene surface follows a modified two-sided Helmholtz model when plotted versus the open circuit potential capacitance, revealing a different electrical double layer mechanism in confinement. This new fundamental understanding lays the foundation for improved functional devices utilizing electrodes and membranes made of two-dimensional materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

SciDAC ISEP: Integrated Simulation of Energetic Particles in Burning Plasmas

The objective of the SciDAC Center for Integrated Simulation of Energetic Particles in Burning Plasmas (ISEP) is to improve physics understanding of energetic particle (EP) confinement and EP interactions with burning thermal plasmas through large-scale simulations. The ISEP center will develop a multiscale and multiphysics ISEP framework for a predictive capability of EP physics and deliver an EP module incorporating both first-principles simulations and high fidelity reduced transport models to the fusion whole device modeling (WDM) project. The ISEP framework will enable us to perform long time, global kinetic simulations of EP physics in burning plasmas, by utilizing the full power of the next generation supercomputers. Our research and development activities will build on fruitful collaborations with computer scientists and applied mathematicians to offer enabling technologies for performance scalability, portability, solvers, coupling for integration with the fusion WDM project, and long-term preservation of data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Assimilation of Multiscale Data into Multifidelity Biogeochemical Models (Final Report)

Quantitative predictions of subsurface processes rely on computational models that capture, with different degrees of fidelity, complex interactions between hydrologic and biogeochemical processes. Molecular- and pore-scale models provide a high-fidelity representation of these processes but are impractical at the field scale. Reduced complexity (e.g., field-scale or data-driven) models sacrifice some degree of fidelity in favor of computational efficiency. Uncertainty and assimilation of data into model predictions, pose a question of model selection: Given a significant difference in computational cost, can a lower-fidelity model be preferable to its higher-fidelity counterpart? Since a predictive model must be computable in reasonable time it has to operate at the field scale, with molecular- and pore-scale data (obtained either from simulations or measurements) determining both the model's structure and parameters. Availability of such multi-resolution data raises a question hitherto undressed in hydrogeology: Can a coarse model dynamically ``learn'' its own structure (e.g., adjust the reaction pathways in its transport module) as more fine-scale data become available during simulations? Our results in multifidelity simulations, data assimilation, and machine learning led us to hypothesize that the answer to these questions is ``yes''. The overarching goal of this project was to confirm this hypothesis by developing scalable, computationally efficient tools for assimilation of data into multiscale models, in which fine-scale data and simulations dynamically inform and autonomously modify coarse-scale models. Our research led to nine manuscripts, three of which have been published and the other six are currently under review.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2025

Integrating Ultra-Coarse-Grained Protein Models into Accessible Workflows for Multiscale Molecular Dynamics

To capture protein conformational transitions using molecular dynamics (MD), several simulation resolutions covering different spatial and temporal scales are typically needed. All-atom (AA) simulations provide fine resolution, but are computationally infeasible for large systems over longer durations. Coarse-grained (CG) and ultra-coarse-grained (UCG) models have a lower resolution and computational cost while still being able to conserve essential protein features. Prior work on a Multiscale Machinelearned Modeling Infrastructure (MuMMI) combined both AA and CG simulations to study RAS-RAF protein interactions, leveraging CG models for longer time scales and using AA to investigate unusual conformations in greater detail. However, MuMMI is still resource-intensive, and this study aims to maximize exploration of the protein conformational space while reducing computational cost. In this paper, we build on prior work that integrates UCG models based on heterogeneous elastic network modeling (hENM) into the MuMMI workflow. We demonstrate that UCG models enable accurate sampling of protein conformations, focusing on simulating RAS-RAF protein interactions. Using higher-resolution CG Martini simulation data, we can automatically refine intramolecular interactions in UCG models. We present a scalable Python package that uses fluctuations observed in higher-resolution CG Martini simulations to estimate bond coefficients of the UCG model. We built novel machine learning-based backmapping methods to recover more detailed CG Martini structures from UCG structures, using diffusion models to learn the mapping between scales. Finally, we present UCG-mini-MuMMI, an accessible and less compute-intensive version of MuMMI as a resource for the scientific community. Incorporating UCG models into MD studies is applicable to a broad range of systems and proteins, and our study offers insights into the advantages and limitations of these methods.

Chemical structure↗
DOE OSTI2020

Calculation of the detonation state of HN 3 with quantum accuracy

HN 3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN 3 with the accuracy of Kohn–Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN 3 . We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn–Sham density-functional theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

The Use of CASES-97 Observations to Assess and Parameterize the Impact of Land-Surface Heterogeneity on Area-Average Surface Heat Fluxes for Large-Scale Coupled Atmosphere-Hydrology Models

To understand the effects of land-surface heterogeneity and the interactions between the land-surface and the planetary boundary layer at different scales, we develop a multiscale data set. This data set, based on the Cooperative Atmosphere-Surface Exchange Study (CASES97) observations, includes atmospheric, surface, and sub-surface observations obtained from a dense observation network covering a large region on the order of 100 km. We use this data set to drive three land-surface models (LSMs) to generate multi-scale (with three resolutions of 1, 5, and 10 kilometers) gridded surface heat flux maps for the CASES area. Upon validating these flux maps with measurements from surface station and aircraft, we utilize them to investigate several approaches for estimating the area-integrated surface heat flux for the CASES97 domain of 71x74 square kilometers, which is crucial for land surface model development/validation and area water and energy budget studies. This research is aimed at understanding the relative contribution of random turbulence versus organized mesoscale circulations to the area-integrated surface flux at the scale of 100 kilometers, and identifying the most important effective parameters for characterizing the subgrid-scale variability for large-scale atmosphere-hydrology models.

Chen, Fei↗
NASA NTRS2023

Foreshock Ion Motion Across Discontinuities: Formation of Foreshock Transients

In the ion foreshock, hot flow anomalies (HFAs) and foreshock bubbles (FBs) are two types of foreshock transients that have the strongest fluctuations, which can disturb the magnetosphere-ionosphere system and increase shock acceleration efficiency. They form due to interaction between the foreshock ions and solar wind discontinuities: the direction of the foreshock ion-driven current and whether it decreases or increases the magnetic field strength behind the discontinuity determine whether the transient's formation can be promoted or suppressed. Thus, to predict the HFA and FB formation and forecast their space weather effects, it is necessary to predict the foreshock ion-driven current direction. In this study, we derive analytical equations of foreshock ion velocities within discontinuities to estimate foreshock ion-driven current direction, which provides a quantitative criterion of HFA and FB formation. To validate the criterion, we use Acceleration Reconnection Turbulence & Electrodynamics of Moon's Interaction with the Sun to observe pristine solar wind discontinuities and calculate discontinuity parameters. We use Magnetospheric Multiscale to observe the foreshock ion motion around the discontinuities and show that the data support our model. This study is another step toward a predictive model of HFA and FB formation so that we can forecast their space weather effects at Earth using solar wind observations at lunar orbit or L1.

Terry Z. Liu↗