SEARCH · Engineering Papers

Results for “multiscale dynamics”

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 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS2006

Modeling the Self-organized Critical Behavior of the Plasma Sheet Reconnection Dynamics

Analyses of Polar UVI auroral image data reviewed in our other presentation at this meeting (V. Uritsky, A. Klimas) show that bright night-side high-latitude UV emissions exhibit so many of the key properties of systems in self-organized criticality (SOC) that an alternate interpretation has become virtually impossible. It is now necessary to find and model the source of this behavior. We note that the most common models of self-organized criticality are numerical sandpiles. These are, at root, models that govern the transport of some quantity from a region where it is loaded to another where it is unloaded. Transport is enabled by the excitation of a local threshold instability; it is intermittent and bursty, and it exhibits a number of scale-free statistical properties. Searching for a system in the magnetosphere that is analogous and that, in addition, is known to produce auroral signatures, we focus on the reconnection dynamics of the plasma sheet. In our previous work, a driven reconnection model has been constructed and has been under study. The transport of electromagnetic (primarily magnetic) energy carried by the Poynting flux into the reconnection region of the model has been examined. All of the analysis techniques, and more, that have been applied to the auroral image data have also been applied to this Poynting flux. Here, we report new results showing that this model also exhibits so many of the key properties of systems in self-organized criticality that an alternate interpretation is implausible. Further, we find a strong correlation between these key properties of the model and those of the auroral UV emissions. We suggest that, in general, the driven reconnection model is an important step toward a realistic plasma physical model of self-organized criticality and we conclude, more specifically, that it is also a step in the right direction toward modeling the multiscale reconnection dynamics of the magnetotail.

Klimas, Alex↗

NASA NTRS2006

Modeling the Self-organized Critical Behavior of Earth's Plasma Sheet Reconnection Dynamics

Analyses of Polar UVI auroral image data show that bright night-side high-latitude W emissions exhibit so many of the key properties of systems in self-organized criticality that an alternate interpretation has become virtually impossible. These analyses will be reviewed. It is now necessary to find and model the source of this behavior. We note that the most common models of self-organized criticality are numerical sandpiles. These are, at root, models that govern the transport of some quantity from a region where it is loaded to another where it is unloaded. Transport is enabled by the excitation of a local threshold instability; it is intermittent and bursty, and it exhibits a number of scale-free statistical properties. Searching for a system in the magnetosphere that is analogous and that, in addition, is known to produce auroral signatures, we focus on the reconnection dynamics of the magnetotail plasma sheet. In our previous work, a driven reconnection model has been constructed and has been under study. The transport of electromagnetic (primarily magnetic) energy carried by the Poynting flux into the reconnection region of the model has been examined. All of the analysis techniques (and more) that have been applied to the auroral image data have also been applied to this Poynting flux. New results will be presented showing that this model also exhibits so many of the key properties of systems in self-organized criticality that an alternate interpretation is implausible. A strong correlation between these key properties of the model and those of the auroral UV emissions will be demonstrated. We suggest that, in general, the driven reconnection model is an important step toward a realistic plasma physical model of self-organized criticality and we conclude, more specifically, that it is also a step in the right direction toward modeling the multiscale reconnection dynamics of the magnetotail.

Klimas, Alexander J.↗

DOE OSTI2025

Dynamics and activation of membrane-bound B cell receptor assembly

B-cell receptor (BCR) complexes are expressed on the surface of a B-cell and are critical in antigen recognition and modulating the adaptive immune response. Even though the relevance of antibodies has been known for almost a hundred years, the antigen-dependent activation mechanism of B-cells has remained elusive. Several models have been proposed for BCR activation, including cross-linking, conformation-induced oligomerization, and dissociation activation models. Recently, the first cryo-EM structures of the human B-cell antigen receptor of the IgM and IgG isotypes have been published that validates the asymmetric organization of the BCR complex. Here, we carry out extensive molecular dynamics simulations to probe the conformational changes upon antigen binding and the influence of the membrane lipids. We identify two critical dynamical events that could be associated with antigen-dependent activation of BCR. First, antigen binding causes increased flexibility in regions distal to the antigen binding site. Second, antigen binding alters the rearrangement of IgM transmembrane helices, including the relative interaction of Igα/Igβ that mediates intracellular signaling. Furthermore, these transmembrane rearrangements lead to changes in localized lipid composition. Our work indirectly supports the conformational-change induced models of BCR activation and contributes to the understanding of the antigen-dependent activation mechanism of BCRs.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Binding and Translocation of Substrate Allosterically Promotes Functional Interactions Within the AlkB–AlkG Electron Transfer Complex

The alkane monooxygenase AlkB and rubredoxin AlkG form an electron transfer complex that hydroxylates terminal alkanes to produce alcohols. The recent cryoEM study of Fontimonas thermophila AlkB-AlkG complex revealed its architecture, including a dodecane (D12) substrate at the active site. However, FtAlkBG molecular mechanism of action of remains unknown. Here, in this study, we examined its dynamics and interactions by multiscale computations, including molecular dynamics simulations, elastic network models, and QM/MM of the oxygen activation mechanism at the AlkB catalytic site. D12 maintained stable interactions within the catalytic site during two MD runs, coordinated by hydrophobic residues L263-L264, I267, I133. A third extended run revealed that D12 could translocate to a membrane-exposed site near S49/F46 along a hydrophobic channel gated by I54. During this translocation, D12 was temporarily stabilized at intermediate sites IS1 (lined by I27/L30-G31/G50/L53-I54/P59/S124/A127-V128) and IS2 (I33-G34/L37/L45-F46/S49) before nearly exiting the protein, and diffused back to the active site, assisted by L30. Substrate binding and translocation across those intermediate sites affects the coupling between the iron centers in AlkBG, and interfacial interactions between AlkB-AlkG. The channel was further connected to the cytosol, near two surface-exposed arginines, potentially allowing for O 2 passage. The allosteric effects between D12 putative entry site, catalytic site and AlkB-AlkG interface were analyzed by ENM-based methods which confirmed the cooperative perturbation-responses and strongly correlated movements of residues belonging to those distal regions. Our study provides new mechanistic insights into key sites and their interactions that could be targeted for developing AlkB-variants with desirable alkane conversion functions.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Seamlessly joining length scales: From atomistic thermal graphs to anisotropic continuum conductivity

Thermal transport in complex solids is governed by local structure, defects, and anisotropy, yet most continuum models still rely on oversimplified and homogenized conductivities. Here, we bridge atomistic and continuum descriptions by building finite element (FE) models directly from the site-projected thermal conductivity (SPTC), an atomic-level decomposition of the Green–Kubo thermal conductivity. We introduce a toolkit, the “Simulator Collection for Atomic-to-Continuum Scales (SCACS)”, which uses a graph neural network to predict SPTC on large atomic structures, coarse-grains these fields into anisotropic conductivity tensors, and embeds them into the heat-flow FE equation with a customized, anisotropy-aware adaptive mesh refinement scheme. Applied to silicon nanostructures, the resulting FE models act as representative volume elements, reproduce bulk conductivities, and capture interfacial and defect-driven anisotropy while maintaining thermodynamic consistency. Additionally, SCACS predicts experimental conductance trends and fields. This work demonstrates a general route for transferring atomistic transport information into device-scale thermal simulations with physics-based approximations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

iMMD

Iterative Multiscale Molecular Dynamics in OpenMM: Accelerating Conformational Sampling of Biomolecular Systems by Iterating All-Atom and Coarse-Grained Molecular Dynamics Simulations

Do, Hung Nguyen [Los Alamos National Laboratory]↗

NASA NTRS1997

Parallel Multiscale Algorithms for Astrophysical Fluid Dynamics Simulations

Our goal is to develop software libraries and applications for astrophysical fluid dynamics simulations in multidimensions that will enable us to resolve the large spatial and temporal variations that inevitably arise due to gravity, fronts and microphysical phenomena. The software must run efficiently on parallel computers and be general enough to allow the incorporation of a wide variety of physics. Cosmological structure formation with realistic gas physics is the primary application driver in this work. Accurate simulations of e.g. galaxy formation require a spatial dynamic range (i.e., ratio of system scale to smallest resolved feature) of 104 or more in three dimensions in arbitrary topologies. We take this as our technical requirement. We have achieved, and in fact, surpassed these goals.

Norman, Michael L.↗

NASA NTRS2008

A New Concurrent Multiscale Methodology for Coupling Molecular Dynamics and Finite Element Analyses

The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.

Yamakov, Vesselin↗

NASA NTRS2018

Multiscale Modeling of PEEK Using Reactive Molecular Dynamics Modeling and Micromechanics

Polyether ether ketone (PEEK) is a high-performance, semi-crystalline thermoplastic that is used in a wide range of engineering applications, including some structural components of aircraft. The design of new PEEK-based materials requires a precise understanding of the multiscale structure and behavior of semi-crystalline PEEK. Molecular Dynamics (MD) modeling can efficiently predict bulk-level properties of single phase polymers, and micromechanics can be used to homogenize those phases based on the overall polymer microstructure. In this study, MD modeling was used to predict the mechanical properties of the amorphous and crystalline phases of PEEK. The hierarchical microstructure of PEEK, which combines the aforementioned phases, was modeled using a multiscale modeling approach facilitated by NASA's MSGMC. The bulk mechanical properties of semi-crystalline PEEK predicted using MD modeling and MSGMC agree well with vendor data, thus validating the multiscale modeling approach.

Pisani, William A.↗

DOE OSTI2026

Hydration-Controlled Proton Transport in Respiratory Complex I

Proton pumping by respiratory complex I is one essential element for generating the proton motive force that drives ATP synthesis in mitochondria. Although it is understood that electrons from NADH reduce ubiquinone at the peripheral arm and that four protons are transferred in the membrane domain, the mechanism by which this redox reaction initiates proton translocation remains unclear. A lateral pathway linking the quinone binding site to the membrane domain via ND1, ND3, and ND4L subunits has been proposed as a possible initial path of an excess proton. However, experimental structures indicate that the hydration connectivity between D66 ND3 and E34 ND4L is comparatively weaker than in neighboring segments, suggesting a potential regulatory point for proton transfer. Using multiscale reactive molecular dynamics (MS-RMD) and a water wire connectivity metric, we directly simulate proton transport through this region as coupled to the hydration by water molecules. Our results reveal that proton transfer is thermodynamically feasible when transient hydration aligns with the presence of an excess proton, revealing the strong coupling between hydration and proton transfer (PT) in this region of Complex I. These findings support a model where proton injection enhances local hydration, dynamically opening the pathway for proton transfer and regulating the onset of proton pumping in Complex I.

Chemical structure↗

DOE OSTI2025

Multi-scale Interaction Mechanism for Edge-Localized-Mode Suppression in the Tokamak Edge

A central challenge in fusion energy is reconciling the high-confinement mode required for reactor performance with the intense intermittent relaxation events it produces, known as edge-localized modes. These instabilities arise in the steep pressure pedestal at the plasma edge when magnetohydrodynamic thresholds are crossed, inflicting damaging heat loads on reactor components. Here, we show that multiscale interactions between microscopic turbulence and macroscopic magnetohydrodynamic modes provide encouraging prospects for self-organized edge-localized modes regulation. Using direct quantitative measurements of multiscale modes, eddy dynamics, and turbulent flux, we show that small-scale electron drift wave turbulence actively scatters the large-scale peeling-ballooning modes. This scattering decorrelates the pressure and velocity fields of the instability, so arresting its growth. Our modeling and theoretical analysis confirm this suppression mechanism is effective even when conventional linear stability thresholds are exceeded. This work establishes a nonlinear principle for edge-localized modes stability, revealing how ambient micro-turbulence can be leveraged to maintain a macro-stable, high-performance pedestal for future fusion reactors.

experimental nuclear physics↗

DOE OSTI2024

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗

NASA NTRS2017

Realistic 3D Radiative Modeling of Turbulent Structure of Moderate-Mass Stars and Sun

Understanding the turbulent dynamics of the Sun and stars is a critical element for interpreting observed processes and phenomena on different scales and for predicting extreme events such as flares and super flares. High-‐resolution observations of the Sun and high-‐fidelity radiative MHD numerical simulations have substantially advanced our understanding of solar and stellar local dynamics and magnetism from the upper convection zone to the atmosphere and corona. However, global modeling of the Sun with such a high degree of realism is currently not affordable due to the extremely high computational cost of resolving the scales in the convection zone. The physics of the deep solar dynamics can be effectively addressed through modeling more massive solar-‐type stars where the convection zone is shallower and the convective overturning time is much shorter than those on the Sun. We present recent 3D realistic simulation results of moderate-‐mass stars and discuss links between solar and stellar dynamics, such as the multiscale structure of granulation, convective overshoot, and others. In particular, these simulations have provided better understanding of the dynamics of the tachocline (the overshoot layer at the bottom of the convection zone) and have explained long-‐ standing results from helioseismology.

Realistic↗

NASA NTRSDate not supplied

Dissecting the Roles of Stability and Dynamic Regimes in Tropical Cloud Feedbacks Simulated from an Upgraded Multiscale Modeling Framework (MMF)

In this study, we analyze the tropical cloud feedback from a 2 K sea surface temperature perturbation experiment performed with a multiscale modeling framework (MMF). We sort the monthly-mean cloud and radiative properties according to circulation and stability regimes. It is found that the dynamic component of cloud feedback for a given dynamic regime dominates the thermodynamic component in terms of the absolute magnitude.

Kuan-Man Xu↗

DOE OSTI2024

FY24 Advanced Computing HPC Annual Report for allocation "vtocei"

This project aims to elucidate the mechanisms underlying the formation of the cathodeelectrolyte interphase (CEI) and its impact on the performance of Li-ion batteries during electrochemical cycling. Our approach integrates advanced multiscale modeling with multimodal characterization to offer a comprehensive understanding of CEI dynamics. We employ a validated multiscale modeling framework to analyze microstructure-dependent transport properties, alongside joint theory-experiment protocols for detailed resolution of interfacial chemistry via spectroscopy.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Application of a temporal multiscale method for efficient simulation of degradation in PEM Water Electrolysis under dynamic operating conditions

Hydrogen is emerging as a vital energy carrier, driven by the need to reduce carbon emissions. Proton Electrolyte Membrane Water Electrolysis (PEMWE) enables hydrogen production under fluctuating renewable power conditions but requires improved understanding and stability of the anode catalyst layer under dynamic operating conditions, especially with low noble metal loadings. Long-term degradation experiments are both time-consuming and costly; therefore, a systematic, model-aided approach is essential. In the present work, a temporal multiscale method is applied to reduce the computational effort of simulating long-term degradation processes in PEMWE, with an exemplary focus on catalyst dissolution. A mechanistic model incorporating the oxygen evolution reaction, catalyst dissolution, and hydrogen permeation from the cathode to the anode was hypothesized and implemented. In this way, the local periodicity of transport and reaction processes in dynamic PEMWE operation, which influence the gradual degradation of the catalyst layer, is captured. The temporal multiscale method significantly reduces the computational effort of simulation, decreasing processing time from hours to mere minutes. This efficiency gain is attributed to the limited evolution of Slow-Scale variables during each period of time P of the Fast-Scale variables. Consequently, simulation is required only until local periodicity is achieved within each Slow-Scale time step. Hence, the fully resolved dynamic problem is decoupled into these two scales, employing a heterogeneous multiscale technique. The developed approach effectively accelerates parameter estimation and predictive simulations, supporting systematic modeling of PEMWE degradation under dynamic conditions.

08 HYDROGEN↗

NASA NTRS2005

Model-Based, Multiscale Self-Tuning Controller Developed for Active Combustion Control

New challenges concerning system health-monitoring and life-extending robust controls for the Ultra-Efficient Engine Technology Project, as well as other advanced engine and power system concepts at NASA and elsewhere, have renewed the control community s interest in smart, model-based methods. In particular, these challenges have further motivated efforts at the NASA Glenn Research Center to exploit the versatility and superiority of the dynamic features extraction of multiscale analysis for controls--such as with "wavelets" and "wavelet filter-banks.' The accomplishments reported herein pertain to the active suppression of combustion instabilities in liquid-fuel combustors via fuel modulation. The fundamentals and initial success of this innovation were reported for a unique demonstration of active combustion control (a research collaboration of NASA Glenn with Pratt & Whitney and the United Technologies Research Center, UTRC). This demonstration, conducted in 2002 at UTRC on the NASA single nozzle rig (SNR) combustor, was the first known suppression of high-frequency instability with a liquid-fueled combustor. The SNR is based on a high-powered military engine combustor that exhibited well-known instabilities.

Le, Dzu K.↗

NASA NTRS2016

Initial Satellite Formation Flight Results from the Magnetospheric Multiscale Mission

This paper describes the underlying dynamics of formation flying in a high-eccentricity orbit such as that of the Magnetospheric Multiscale mission. The GPS-based results used for MMS navigation is summarized, as well as the procedures that are used to design the maneuvers used to place the spacecraft into a tetrahedron formation and then maintain it. The details of how to carry out these maneuvers are then discussed. Finally, the numerical results that have been obtained concerning formation flying for the MMS mission to date (e.g. tetrahedron sizes flown, maneuver execution error, fuel usage, etc.) are presented in detail.

Conjunction↗