Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 253 records · Page 14

Lower length scale model for palladium attack of silicon carbide in TRISO fuel

TRistructural ISOtropic (TRISO) particle fuels rely on silicon carbide (SiC) as the primary barrier for metallic fission product (FP) release. Palladium (Pd) generated by fission degrades the SiC layer, resulting in the formation of lamellar layers of palladium silicides (PdxSi) and carbon (C) perpendicular to the direction of attack. The Pd attack has been hypothesized to be responsible for failure of the SiC layer and enhance FP release. To better understand and quantify Pd attack of SiC in TRISO particles, a multiscale, mechanistic model of Pd transport is being developed by the NEAMS program. Previous work provided an initial hypothesis for modeling lamellar microstructure formation in SiC due to Pd attack using a phase-field model. The work described in this report builds on the previous model by using molecular dynamics (MD) simulations to parameterize the phase-field model kinetics, and build a reduced order model in BISON using the improved mesoscale Pd penetration model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Effects of Multiscale Phase-Mixing and Interior Conductance in the Lunar-like Pickup Ion Plasma Wake. First Results from 3-D Hybrid Kinetic Modeling

The study of multiscale pickup ion phase-mixing in the lunar plasma wake with a hybrid model is the main subject of our investigation in this paper. Photoionization and charge exchange of protons with the lunar exosphere are the ionization processes included in our model. The computational model includes the self-consistent dynamics of the light (hydrogen ions or hydrogen molecule ion or helium ion), and heavy (sodium ion ) pickup ions. The electrons are considered as a fluid. The lunar interior is considered as a weakly conducting body. In this paper we considered for the first time the cumulative effect of heavy neutrals in the lunar exosphere (e.g., Aluminum, Argon), an effect which was simulated with one species of sodium ion but with a tenfold increase in total production rates. We find that various species produce various types of plasma tail in the lunar plasma wake. Specifically, sodium ion and helium ion pickup ions form a cycloid-like tail, whereas the hydrogen ion or hydrogen molecule ion pickup ions form a tail with a high density core and saw-like periodic structures in the flank region. The length of these structures varies from 1:5 R(sub M) to 3:3 R(sub M) depending on the value of gyro radius for hydrogen ion or hydrogen molecule ion pickup ions. The light pickup ions produce more symmetrical jump in the density and magnetic field at the Mach cone which is mainly controlled by the conductivity of the interior, an effect previously unappreciated. Although other pickup ion species had little effect on the nature of the interaction of the Moon with the solar wind, the global structure of the lunar tail in these simulations appeared quite different when the hydrogen molecule ion production rate was high.

Lipatov, A.S.↗

Strong Reconnection Electric Fields in Shock-Driven Turbulence

Turbulent magnetic reconnection in a quasi-parallel shock under parameters relevant to the Earth's bow shock is investigated by means of a two-dimensional particle-in-cell simulation. The addressed aspects include the reconnection electric field, the reconnection rate, and the electron and the ion outflow speeds. In the shock transition region, many current sheets are generated in shock-driven turbulence, and electron-only reconnection and reconnection where both ions and electrons are involved can occur in those current sheets. The electron outflow speed in electron-only reconnection shows a positive correlation with the theoretical speed, which is close to the local electron Alfvén speed, and a strong convection electric field is generated by the large electron outflow. As a result, the reconnection electric field becomes much larger than those in the standard magnetopause or magnetotail reconnection. In shock-driven reconnection that involves ion dynamics, both electron outflows and ion outflows can reach of the order of 10 times the Alfvén speed in the X-line rest frame, leading to a reconnection electric field the same order as that in electron-only reconnection. An electron-only reconnection event observed by the magnetospheric multiscale mission downstream of a quasi-parallel shock is qualitatively similar to those in the simulation and shows that the outflow speed reaches approximately half the local electron Alfvén speed, supporting the simulation prediction.

magnetic reconnection↗

The Development of Directional Decohesion Finite Elements for Multiscale Failure Analysis of Metallic Polycrystals

Atomistic simulations of intergranular fracture have indicated that grain-scale crack growth in polycrystalline metals can be direction dependent. At these material length scales, the atomic environment greatly influences the nature of intergranular crack propagation, through either brittle or ductile mechanisms, that are a function of adjacent grain orientation and direction of crack propagation. Methods have been developed to obtain cohesive zone models (CZM) directly from molecular dynamics simulations. These CZMs may be incorporated into decohesion finite element formulations to simulate fracture at larger length scales. A new directional decohesion element is presented that calculates the direction of Mode I opening and incorporates a material criterion for dislocation emission based on the local crystallographic environment to automatically select the CZM that best represents crack growth. The simulation of fracture in 2-D and 3-D aluminum polycrystals is used to illustrate the effect of parameterized CZMs and the effectiveness of directional decohesion finite elements.

Saether, Erik↗

The Magnetospheric Multiscale Constellation

The Magnetospheric Multiscale (MMS) mission is the fourth mission of the Solar Terrestrial Probe (STP) program of the National Aeronautics and Space Administration (NASA). The MMS mission was launched on March 12, 2015. The MMS mission consists of four identically instrumented spin-stabilized observatories which are flown in formation to perform the first definitive study of magnetic reconnection in space. The MMS mission was presented with numerous technical challenges, including the simultaneous construction and launch of four identical large spacecraft with 100 instruments total, stringent electromagnetic cleanliness requirements, closed-loop precision maneuvering and pointing of spinning flexible spacecraft, on-board GPS based orbit determination far above the GPS constellation, and a flight dynamics design that enables formation flying with separation distances as small as 10 km. This paper describes the overall mission design and presents an overview of the design, testing, and early on-orbit operation of the spacecraft systems and instrument suite.

Space mission design↗

A MPET 2 -mPBPK model for subcutaneous injection of biotherapeutics with different molecular weights: From local scale to whole-body scale

Subcutaneous injection of biotherapeutics has attracted considerable attention in the pharmaceutical industry. However, there is limited understanding of the mechanisms underlying the absorption of drugs with different molecular weights and the delivery of drugs from the injection site to the targeted tissue. Here, we propose the MPET 2 -mPBPK model to address this issue. This multiscale model couples the MPET 2 model, which describes subcutaneous injection at the local tissue scale from a biomechanical view, with a post-injection absorption model at injection site and a minimal physiologically-based pharmacokinetic (mPBPK) model at whole-body scale. Utilizing the principles of tissue biomechanics and fluid dynamics, the local MPET 2 model provides solutions that account for tissue deformation and drug absorption in local blood vessels and initial lymphatic vessels during injection. Additionally, we introduce a model accounting for the molecular weight effect on the absorption by blood vessels, and a nonlinear model accounting for the absorption in lymphatic vessels. The post-injection model predicts drug absorption in local blood vessels and initial lymphatic vessels, which are integrated into the whole-body mPBPK model to describe the pharmacokinetic behaviors of the absorbed drug in the circulatory and lymphatic system. We establish a numerical model which links the biomechanical process of subcutaneous injection at local tissue scale and the pharmacokinetic behaviors of injected biotherapeutics at whole-body scale. With the help of the model, we propose an explicit relationship between the reflection coefficient and the molecular weight and predict the bioavalibility of biotherapeutics with varying molecular weights via subcutaneous injection. The considered drug absorption mechanisms enable us to study the differences in local drug absorption and whole-body drug distribution with varying molecular weights. This model enhances the understanding of drug absorption mechanisms and transport routes in the circulatory system for drugs of different molecular weights, and holds the potential to facilitate the application of computational modeling to drug formulation.

59 BASIC BIOLOGICAL SCIENCES↗

Conjunction Assessment Techniques and Operational Results from the Magnetospheric Multiscale Mission

This paper will describe the results that have been obtained to date during the MMS mission concerning conjunction assessment. MMS navigation makes use of a weak-signal GPS-based system: this allows signals to be received even when MMS is flying above the GPS orbits, producing a highly accurate determination of the four MMS orbits. This data is downlinked to the MMS Mission Operations Center (MOC) and used by the Flight Dynamics Operations Area (FDOA) for both maneuver design and conjunction assessment. The MMS fly in tetrahedron formations around apogee, in order to collect simultaneous particles and fields science data. The original plan was to fly tetrahedra between 10 and 160 km in size; however, after Phase 1a of the mission, the science team requested that smaller sizes be flown if feasible. After analysis (to be detailed in a companion paper), a new minimum size of 7 km was decided upon. Flying at this reduced scale size makes conjunction assessment between the MMS spacecraft even more important: the methods that are used by the MMS FDOA to address this problem will be described in the paper, and a summary given of the previous analyses that went into the development of these techniques. Details will also be given of operational experiences to date. Finally, two CA mitigation maneuver types that have been designed (but never yet required to actually be performed) will also be outlined.

conjunction Assessment↗

A Multiscale Cohesive Law for Carbon Fiber Networks

Better predictive models of mechanical failure in low-weight heat shield composites would aid material certification for missions with aggressive atmospheric entry conditions. Here, we develop such a model for the rapid engineering analysis of the failure limits of phenolic impregnated carbon ablator (PICA) - a leading heat shield material whose structural component is a carbon fiber network. We hypothesize inelastic deformation failure mechanisms and model their behavior using molecular dynamics simulations to calculate the binding energy. We then upscale this binding energy to the macroscale using a renormalization argument. The approach delivers insightful and reasonably accurate macroscale predictions that compare favorably to experiments. In application, the model is validated for a particular variety of PICA by comparison to experiment and would then be used to study design scenarios in different entry conditions.

W. SCHILL↗

Influences of shock imprinting on mix in a 3D-printed porous media

Mixing of materials in porous media can cause a significant impact on fusion yield as previously demonstrated by the National Ignition Facility (NIF) MARBLE Campaign. Initially, the reactants are separated, with deuterium in the lattice struts and a tritium gas fill in the voids. Lattice parameters such as the strut thickness and relative pitch, provide a control for the mix parameters in the experiment. Los Alamos National Laboratory’s (LANL) BOSQUE project looks to better understand how the mix of the reactants and shell materials impact the fusion burn and resultant yield on various laser platforms. xRAGE’s Eulerian hydrodynamics and adaptive mesh refinement (AMR) provide the unique ability to study the impacts of multiscale features of complex lattice structures. This modeling provides the ability to measure shock front variations as the wave progress’ through a given media. Initial conditions of the lattice are essential to accurately model mix and burn measured by experiment. By varying the initial orientation and densities of these lattice regions the early time dynamics of how the shock is launched into the system is changed and advocates for the study of resulting effects. In this work, we will study the sensitivities of shock effects in varying 3D printed geometric systems and how these shocks alter the structure and mix in the lattice. We will discuss both preliminary experimental results and simulations to help plan and constrain future experiments where we will study the impact of different lattice geometries and lattice bulk densities. This work concludes with the relative impacts of lattice geometries on shock speeds at different bulk densities and the resultant mix widths due to those shock interactions. We see agreement with theory at the higher end of our bulk density study, however, as we approach lower bulk density systems the dynamics of these interactions begin to change.

2 Photon Polymerization↗

Direct Imaging of Hydrogen‐Driven Dislocation and Strain Field Evolution in a Stainless Steel Grain

Hydrogen embrittlement (HE) poses a significant challenge to the durability of materials used in hydrogen production and utilization. Disentangling the competing nanoscale mechanisms driving HE often relies on simulations and electron-transparent sample techniques, limiting experimental insights into hydrogen-induced dislocation behavior in bulk materials. This study employs in situ Bragg coherent X-ray diffraction imaging to track three-dimensional (3D) dislocation and strain field evolution during hydrogen charging in a bulk grain of austenitic 316 stainless steel. Tracking a single dislocation reveals hydrogen-enhanced mobility and relaxation, consistent with dislocation dynamics simulations. Subsequent observations reveal dislocation unpinning and climb processes, likely driven by osmotic forces. Additionally, nanoscale strain analysis around the dislocation core directly measures hydrogen-induced elastic shielding. These findings experimentally validate theoretical predictions and offer mechanistic insights into hydrogen-driven dislocation behavior. The quantified nanoscale phenomena serve as critical inputs for multiscale modeling frameworks to predict bulk material responses and accelerate the development of HE-resistant alloys.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Millimeter-wave dielectric tunability driven by topological polar structure switching in PbTiO 3 /SrTiO 3 superlattices

Dielectric tunability induced by an external electric field in materials underpins radio frequency signal modulation devices such as phase shifters, which are critical components in wireless communication and sensing systems. However, the tunability and integrability of current devices have yet to be enhanced for emerging applications, particularly at millimeter-wave frequencies. Here, we demonstrate that topological polar structures formed in PbTiO 3 /SrTiO 3 superlattices exhibit large tunable in-plane dielectric properties, as determined by their multiscale structural configurations and polarization switching behaviors. Under a moderate field of 30 kV cm −1 , the dipole wave structure maintains a tunability exceeding 15% at 70 GHz and above 8% over the measured range up to 110 GHz, contrasting with the weakly tunable flux closure structure. Based on in situ structural characterizations and molecular dynamics simulations, we delineate the polarization switching processes and elucidate the mechanisms underlying the observed tunable millimeter-wave dielectric responses. Our results provide new insights into the high-frequency dielectric properties of topological polar phases, potentially broadening the versatility of these materials in next-generation integrated electronic applications.

36 MATERIALS SCIENCE↗

Temperature dependence of vacancy/self-interstitial recombination volumes in copper

We use molecular dynamics to calculate rate coefficients for recombination of vacancies with self-interstitial atoms (SIAs) in Cu at temperatures from 300 to 700 K. From these results, we calculate vacancy/SIA recombination volumes and find that they decrease from around 290 (where is one atomic volume) at 300 K to 160 at 500 K and above. By counting the number of distinct pathways by which a stable SIA may migrate to a site of spontaneous recombination with a nearby vacancy, we find a lower bound estimate of 168 for the recombination volume. We furthermore rationalize its temperature dependence based on differences between the activation energies for recombination and SIA migration. Furthermore, our work sheds light on the fundamental nature of the recombination volume and provides information that may be incorporated into multiscale models of radiation response in solids.

36 MATERIALS SCIENCE↗

Multiscale Modeling of Grain Boundaries in ZrB2: Structure, Energetics, and Thermal Resistance

A combination of ab initio, atomistic and finite element methods (FEM) were used to investigate the structures, energetics and lattice thermal conductance of grain boundaries for the ultra high temperature ceramic ZrB2. Atomic models of idealized boundaries were relaxed using density functional theory. Information about bonding across the interfaces was determined from the electron localization function. The Kapitza conductance of larger scale versions of the boundary models were computed using non-equilibrium molecular dynamics. The interfacial thermal parameters together with single crystal thermal conductivities were used as parameters in microstructural computations. FEM meshes were constructed on top of microstructural images. From these computations, the effective thermal conductivity of the polycrystalline structure was determined.

Lawson, John W.↗

Wave-Particle Energy Exchange Directly Observed in a Kinetic Alfven-Branch Wave

Alfven waves are fundamental plasma wave modes that permeate the universe. At small kinetic scales they provide a critical mechanism for the transfer of energy between electromagnetic fields and charged particles. These waves are important not only in planetary magnetospheres, heliospheres, and astrophysical systems, but also in laboratory plasma experiments and fusion reactors. Through measurement of charged particles and electromagnetic fields with NASAs Magnetospheric Multiscale (MMS) mission, we utilize Earths magnetosphere as a plasma physics laboratory. Here we confirm the conservative energy exchange between the electromagnetic field fluctuations and the charged particles that comprise an undamped kinetic Alfven wave. Electrons confined between adjacent wave peaks may have contributed to saturation of damping effects via non-linear particle trapping. The investigation of these detailed wave dynamics has been unexplored territory in experimental plasma physics and is only recently enabled by high-resolution MMS observations.

magnetohydrodynamic↗

Tandem bulk oxygen diffusion and surface reactions in reducible metal oxides control redox cycle dynamics

The interplay between bulk oxygen diffusion and surface reactions in reducible metal oxides is key in heterogeneous catalysts, but direct measurements of oxygen mobility, transient kinetics, and in situ spectroscopies have been lacking. Here, we reveal complex dynamic behavior of ceria-zirconia by H 2 using transient kinetics via mass spectrometry and in situ Raman and near-ambient pressure x-ray photoelectron spectroscopies. Molecular dynamics simulations with a machine learning potential delineate competitive oxygen diffusion mechanisms, with an optimal mobility at intermediate reductions. We expose a compensation between vacancy availability and lattice distortion at intermediate to high reductions and Frenkel defects at low reductions, underscoring a potential deficiency of 16 O/ 18 O exchange experiments in deducing oxygen mobility. Vacancies in proximity require electron localization on Ce atoms further away. The continuous replenishment of surface oxygen results in a varying reduction rate, with H 2 dissociation being the rate-limiting step. Multiscale transient simulations, consistent with experiments, indicate catalysts of potentially spatially varying oxidation states. The approach is broadly applicable to reducible oxide materials.

36 MATERIALS SCIENCE↗

Crystallization Mechanisms of Poly(vinylidene Fluoride- co -chlorotrifluoroethylene): Nucleation Transitions, Growth Kinetics, and Microstructure Evolution

The crystallization kinetics of FK-800, a commercial semicrystalline copolymer of chlorotrifluoroethylene and vinylidene fluoride (poly(CTFE-co-VDF)), dictate its microstructure and functional performance. Crystallization occurs between the glass transition (T g ≈ 31 °C) and melting point (T m ≈ 110 °C), with grain boundaries playing a key role in applications such as memristors, where they regulate filament growth. This study integrates in situ atomic force microscopy (AFM), grazing-incidence wide-angle X-ray scattering (GIWAXS), and differential scanning calorimetry (DSC) to characterize nucleation, growth, and morphology across the full crystallization window. Hot-stage AFM reveals a transition from homogeneous to heterogeneous nucleation near 45 °C and a shift from reaction-limited to diffusion-limited growth above 60 °C. Kinetics are analyzed using Lauritzen–Hoffman, Turnbull–Fisher, and Avrami models to extract surface free energies and activation barriers. A multiscale modeling framework is developed in which nucleation densities and growth rates measured by AFM are used to reconstruct Avrami kinetics, which are then validated against coverage data, and extrapolated to predict bulk crystallization behavior observed by DSC. GIWAXS confirms a primarily flat-on chain orientation at all temperatures, linking orientation to domain morphology. This integrative approach quantitatively connects nanoscale crystallization dynamics with macroscopic phase evolution, establishing a quantitative framework for modeling crystallization across time, temperature, and length scales.

Crystallization↗

Collaborative Research: Vlasov-Maxwell Simulations to Resolve Electron Heating and Dissipation, in Quasi-Perpendicular Shocks (Final Report)

Collisionless shocks are a grand challenge problem in plasma physics and have been the subject of study for more than six decades. Shocks are ubiquitous phenomena in the universe and are responsible for transforming high energy and momentum flows into thermal energy and energetic particles. Understanding how shocks operate is of primary importance to understand the sun-earth coupling, protecting manned missions and spacecrafts from high energy particles, inertial confinement fusion, and radiation observed from astrophysical plasmas, such as supernova remnants and astrophysical jets. A major unanswered question on this frontier is, how does a collisionless plasma transform flow energy and momentum into electron thermal energy? Many potential mechanisms have been proposed to perform the conversion between flow and thermal energy in a collisionless plasma, but the answer has been elusive. To solve the mystery, non-relativistic perpendicular and quasi-perpendicular shocks are studied via a coordinated approach employing cutting-edge numerical simulations with the continuum kinetic code Gkeyll and spacecraft observations supplied by NASA's Magnetospheric Multiscale (MMS) mission. In both sets of data, we employ novel techniques to identify how and where energy is exchanged between electromagnetic fields and plasma particles. By combining cutting-edge spacecraft data and numerical simulations with these diagnostics, we have been able to identify some of the primary mechanisms responsible electron heating and dissipation in heliospheric shocks. In addition to electron dynamics, we have complete access to proton simulation and in situ data, enabling us to characterize proton heating and acceleration in quasi-perpendicular shocks. Many astrophysical quasi-perpendicular shocks have similar parameters to those in this study, including supernova remnants and shocks in galaxy clusters. Additionally, the simulation code, Gkeyll, and associated analysis scripts are openly available for community use. This award has supported a collaborative effort between Princeton University, the University of Arizona - Tucson, and the University of Maryland - College Park. The award funded three early career scientists, including the funding and training of one female postdoctoral researcher, who was involved in all aspects of the MMS data analysis and some portions of the simulation analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Kalman Filter for Mass Property and Thrust Identification of the Spin-Stabilized Magnetospheric Multiscale Formation

The Magnetospheric Multiscale (MMS) mission consists of four identically instrumented, spin-stabilized observatories, elliptically orbiting the Earth in a tetrahedron formation. For the operational success of the mission, on-board systems must be able to deliver high-precision orbital adjustment maneuvers. On MMS, this is accomplished using feedback from on-board star sensors in tandem with accelerometers whose measurements are dynamically corrected for errors associated with a spinning platform. In order to determine the required corrections to the measured acceleration, precise estimates of attitude, rate, and mass-properties are necessary. To this end, both an on-board and ground-based Multiplicative Extended Kalman Filter (MEKF) were formulated and implemented in order to estimate the dynamic and quasi-static properties of the spacecraft.

Formation Flying↗