Search NASA⌕ Search

SEARCH · Search NASA

Results for “CHARGE DISTRIBUTION”

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.

317 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Multiscale and multidimensional modeling of particle acceleration and transport in solar flares

Multi-messenger, multi-viewpoint, and time-resolved observations of solar flares are now providing unprecedented constraints on particle acceleration sites, energy conversion, and energy transport. The interpretation of current observations, including microwave imaging spectroscopy from EOVSA, hard x-ray (HXR) imaging from Solar Orbiter/STIX, gamma-ray diagnostics from Fermi, and in situ measurements from Parker Solar Probe and Solar Orbiter, collectively demands modeling frameworks that go beyond traditional spatially unresolved, one-zone models or single-mechanism descriptions. This review surveys multiscale and multidimensional modeling approaches, including kinetic, magnetohydrodynamic (MHD), and macroscopic particle models, that are being developed to meet the need. Kinetic simulations reveal that three-dimensional (3D) effects, including field-line chaos and self-generated turbulence, are essential for sustained power-law particle acceleration. MHD simulations now capture flux-rope eruptions, plasmoid-unstable current sheets, and turbulent flare regions in realistic magnetic topologies. Macroscopic models coupling MHD with energetic-particle models produce spatially resolved electron distributions and synthetic HXR and microwave emissions for direct comparison with observations. Despite these advances, outstanding challenges remain in bridging kinetic and global scales, improving MHD simulations and macroscopic particle models, and achieving quantitative model-observation closure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Simulation to a Newborn Supernova Remnant from a Low-mass Iron Core Star

Supernova remnant observations show a high degree of asymmetry, mixing, and inhomogeneity. These asymmetries are seeded during the early seconds of the explosion and are further enhanced and modified as the shock and ejecta move through the stellar progenitor and into the circumstellar medium. We present simulations of a 9.6 M⊙ zero-metallicity progenitor initialized after shock revival and evolved for several years when the ejecta is in the circumstellar medium. A suite of 1D and 2D simulations examines the effects of neutron-star wind and radioactive decay heating. In 1D, decay heating forms a low-density bubble that suppresses the reverse shock. While in 2D, the heating is localized to metal-rich pockets, inflating them and compressing the surrounding material into dense shells. In 3D, the neutron-star wind and decay heating modify the plume morphology, producing more large-scale structures. The extended plume morphology leads to an asymmetrical shock breakout. After breakout, the leading plumes cannot keep up with the shock front, resulting in deceleration and fragmentation by the reverse shock while retaining the large-scale asymmetry. The projected ejecta morphology and velocities are strongly viewing angle dependent. The relatively uniform metal-rich distribution does not resemble the strongly inhomogeneous ejecta structure of Cas A. The 160-isotope decay network shows that 24.4% of the radioactive heating comes from decay chains other than the canonical 56Ni chain. The low explosion energy, low 56Ni yield, and Ni/Fe ratio greater than unity suggest an observational signature similar to an electron capture supernova.

Neopane, Sudarshan [University of Tennessee (UT)]↗
DOE OSTI2024

Case study of UAS ignition of prescribed fire in a mixedwood on the William B. Bankhead National Forest, Alabama

Abstract Background For at least four decades, practitioners have recognized advantages of aerial versus ground ignition for maximizing the effectiveness of prescribed fires. For example, larger areas can be ignited in less time, or ignition energy may be variously targeted over an area in accordance with the uneven distribution of fuels. The maturation of wireless communication, geopositioning systems, and unmanned aerial systems (UAS) has enhanced those advantages, and UAS approaches also provide further advantages relative to helicopter ignitions, such as reduced risk to human safety, lower operating costs, and higher operational flexibility. In a long running study at the Bankhead National Forest in northcentral Alabama, prescribed fire has been used for nearly 20 years. Most of the burns have been hand-ignited via drip torches, while some have been aerially ignited via helicopter. In March 2022, for the first time, a UAS was used to ignite prescribed fires across a landscape that included a long-term research stand. This field note relates comparisons of both fire behavior and fuel consumption metrics for the UAS-ignited burn versus previous burns on the same stand, and versus burns of other research stands in the same year. Results The UAS-ignited prescribed fire experienced burn effects similar to those from ground-ignited prescribed fires on the same stand in previous years, as well as those from ground-ignited prescribed fires on other stands in the same year. Conclusion This post hoc analysis suggests that UAS ignition approaches may be sufficient for achieving prescribed burn goals, thereby enabling practitioners to realize the advantages offered by that ignition mode.

Environmental Sciences & Ecology↗
DOE OSTI2026

High-resolution modeling of indoor radon exposure with uncertainty quantification in Utah

Indoor radon accounts for 37% of population-level exposure to ionizing radiation in the United States. However, radon metrics are typically reported at coarse spatial scales, potentially obscuring meaningful local variation. We developed a high-resolution modeling framework to estimate indoor radon concentrations across Utah while explicitly quantifying predictive uncertainty. A total of 19,497 residential radon measurements collected between 2006 and 2017 were combined with environmental and housing characteristics and analyzed using a geospatial neural network that accommodates spatial dependence and nonlinear associations. Predictions were generated on a uniform hexagonal grid at 0.73 km2 resolution (H3 level 8). Out-of-sample predictions aggregated to the H3 level 8 grid showed good agreement with observed concentrations (Pearson r=0.64), while household-level predictions exhibited more moderate agreement (r=0.45). The model produced well-calibrated uncertainty estimates, with 24.1% of held-out observations exceeding the predicted 75th-percentile threshold. Maps of predicted radon concentrations and the probability of exceeding the U.S. EPA action level of 148 Bq/m3 (4 pCi/L) revealed substantial fine-scale spatial heterogeneity that was not apparent in conventional coarse-resolution summaries, with greater local variability observed in densely monitored urban counties than in sparsely sampled regions. High-resolution radon models that explicitly quantify uncertainty provide a useful framework for characterizing the spatial distribution of indoor radon and identifying areas of elevated exceedance risk. These findings highlight the value of fine-scale monitoring data and uncertainty-aware modeling approaches for radon exposure assessment, environmental risk characterization, and radon-related health research.

Wu, Yunhan [ORNL] (ORCID:0000000178842994)↗
DOE OSTI2026

Prioritizing Uncertainties in Hydrogen Contribution to Risk in Post-Crash Outcomes for Rail

This report presents analysis from Sandia National Laboratories predicting contributions to risk associated with the use of hydrogen technology for rail. Event sequence diagrams are used to describe possible accident scenarios and progressions. Initiating event frequencies and branch event probabilities for each scenario are quantified with uncertainty using distributions fit to Federal Railroad Administration and U.S. Department of Transportation Pipeline and Hazardous Materials Safety Administration data on applicable accidents from 2000 to 2020. Uncertainty is propagated through the event sequence diagram to estimate the frequency and conditional probability of accident end states. The analysis identifies four scenarios with significant contributions to risk from hydrogen that are predicted to occur relatively frequently, which may inform priorities for reducing uncertainty. These scenarios are 1) overpressure events resulting from collisions with hydrogen release due to mechanical damage and delayed ignition, 2) jet fire events resulting from collisions with hydrogen release due to mechanical damage and immediate ignition, 3) jet fires resulting from fire or explosion initiating events involving the hydrogen tank and correct operation of the thermally-activated pressure relief device (TPRD) subsequent to the thermal insult, and 4) pressure burst resulting from fire or explosion initiating events involving the hydrogen tank and failure of the TPRD. Delayed and immediate hydrogen ignition probabilities are identified as being highly uncertain and potential candidates for reducing conservatism in the predicted frequencies for these two scenarios.

08 HYDROGEN↗
DOE OSTI2024

Sustainable Production of Biomass‐Derived Graphite and Graphene Conductive Inks from Biochar

Abstract Graphite is a commonly used raw material across many industries and the demand for high‐quality graphite has been increasing in recent years, especially as a primary component for lithium‐ion batteries. However, graphite production is currently limited by production shortages, uneven geographical distribution, and significant environmental impacts incurred from conventional processing. Here, an efficient method of synthesizing biomass‐derived graphite from biochar is presented as a sustainable alternative to natural and synthetic graphite. The resulting bio‐graphite equals or exceeds quantitative quality metrics of spheroidized natural graphite, achieving a RamanI D /I G ratio of 0.051 and crystallite size parallel to the graphene layers (L a ) of 2.08 µm. This bio‐graphite is directly applied as a raw input to liquid‐phase exfoliation of graphene for the scalable production of conductive inks. The spin‐coated films from the bio‐graphene ink exhibit the highest conductivity among all biomass‐derived graphene or carbon materials, reaching 3.58 ± 0.16 × 10 4 S m −1 . Life cycle assessment demonstrates that this bio‐graphite requires less fossil fuel and produces reduced greenhouse gas emissions compared to incumbent methods for natural, synthesized, and other bio‐derived graphitic materials. This work thus offers a sustainable, locally adaptable solution for producing state‐of‐the‐art graphite that is suitable for bio‐graphene and other high‐value products.

Chemistry↗
DOE OSTI2024

Circumventing data imbalance in magnetic ground state data for magnetic moment predictions

Abstract Magnetic materials play a crucial role in the transition to more sustainable forms of energy and electric vehicles. There is an anticipated shortage in magnetic materials in the future, and as a result there is an urgent need to discover and design new magnetic materials. Computational magnetic material design using density functional theory is daunting because of the challenge in identifying magnetic ground states from a combinatorially large set of possibilities. Machine learning offers a path forward by enabling efficient surrogate models that can more readily enumerate these states, but there is a dearth of training data available, and what is available tends to be imbalanced with too much non-magnetic data. In this work we show that the discrete and previously tackled data imbalance that exists at the level of the magnetic ordering leads to an imbalanced continuous distribution with many zeros when the data is unraveled at the atomic magnetic moment level, which subsequently leads to models with low accuracy for magnetic properties. We mitigate this by using a two-part model framework. Our scheme is able to classify atoms into magnetic and non-magnetic with an F1 score and Matthew’s correlation coefficient (MCC) of ~91% and then to provide an implicit embedding representation that maps directly onto the magnitude of the magnetic moment with a mean absolute error of 0.1 μ B . Beyond screening for new magnetic materials, we demonstrate an additional practical use case of our scheme: the provision of good initial guesses for magnetic moments in first-principles electronic relaxations. Such initialization is shown to lead to faster convergence to configurations that lie closer to the ground state.

Computer Science↗
DOE OSTI2026

Selectivity mechanisms of ion intercalation in Prussian blue analogs

Prussian blue analogs (PBAs) are a family of materials with facile, reversible, and selective ion transport capability for various ions via electrochemical intercalation, owing to their vacancy structure. The large tunable compositional space of PBAs allows for manipulation of intercalation behavior and selectivity by controlling structural vacancy level through choice of transition metal centers and modifications to the synthesis process. However, a lack of understanding of the mechanisms of ion selectivity hinders the material’s design process. Here, for this work, we investigated the origins of ion selectivity using a model PBA, copper hexacyanoferrate, and focused on eight technologically and biologically prominent ions, for which we determined a sequence of selectivity: Rb + > K + > Na + > Ba 2+ > Sr 2+ ≈ Ca 2+ > Mg 2+ > Li + . We provide electrochemical, structural, and redox evidence of strong correlation between the ion identity, the dominant charge-compensating redox, and preferred occupancy site. Specifically, using synchrotron anomalous X-ray diffraction (AXRD), we reveal that monovalent ions exhibit significant association with the corner sites of the unit cell and iron redox, whereas divalent ions display affinity toward the center site with higher ratios of copper redox. Informed by selectivity results, we applied CuHCFe to Li purification and achieved 99.9% purity. Our findings demonstrate an approach to elucidating ion intercalation behavior in order to distinguish and manipulate material properties to optimize separation performance.

Prussian blue analog↗
DOE OSTI2026

Thermally Stable Co@C3N4 Single-Atom Catalysts for CO Oxidation: Atomic-Level Insights into Structure and Activity

Single Co atoms supported on C3N4 (Co@C3N4) have demonstrated high activity and selectivity in photocatalysis. However, the investigation of structure–function relationships and reaction mechanisms under photocatalytic conditions is very challenging due to the complex conditions of light absorption, charge transfer, and catalysis. In this study, we employed thermal CO oxidation as a prototypical probe reaction to benchmark the intrinsic catalytic performance and track the active-site evolution of Co@C3N4. Single Co atoms were identified and shown to be the catalytically active sites for CO oxidation based on control experiments and isotope-labeling experiments. The Co sites remained atomically dispersed before, during, and after the reaction with temperatures up to 400 °C, as established by in situ X-ray absorption fine structure (XAFS) combined with density functional theory (DFT), FDMNES simulations, and dynamic-time-warping (DTW)-assisted X-ray absorption near edge structure (XANES) matching. Together with theoretical calculations, the integrated analysis reveals a stable coordination environment under reaction conditions, which correlates with sustained activity, establishing Co@C3N4 single-atom catalysts as thermally stable CO oxidation catalysts. Beyond these findings, the current study provides a workflow for unambiguously assigning active sites in Co@C3N4 for thermal CO oxidation. This workflow will aid the understanding of their behavior in photocatalysis in the future, where light-driven dynamics obscure direct structure–function links. Notably, this study provides fundamental insights for the rational design of robust single-atom catalysts and a foundation for the broader application of Co@C3N4 catalysts in oxidation reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Physical Modeling and Design of a Nonvolatile Optically Gated High‐Power Diamond Transistor

In this work, we present the theory and modeling framework of a diamond optically gated junction field‐effect transistor (DOGFET). The device utilizes nitrogen substitutional centers in type‐1b diamond to optically modulate a p‐ boron doped diamond channel. Using sub‐gap lasers with intensities as low as 100 W/cm 2 , electrons are optically excited from substitutional nitrogen sites to the conduction band of the diamond substrate, thus enabling the optical gate to exercise control on modulating the space‐charge region at the junction and therefore the channel conductivity. We show that the device can deliver a current of 7 μA/μm, or equivalently 1750 A/cm 2 , while switching at a frequency greater than 100 kHz, in a form factor of 5 μm 2 . The breakdown voltage is found to be greater than 1850 V, with a breakdown field strength of ~13 MV/cm. Moreover, the device supports nonvolatile operation with a “memory effect” enabling single transistor state retention. The presented simulation framework provides a physically grounded insight into the limits and opportunities of optoelectronic diamond systems.

Engineering - Electronic and electrical engineerin↗
DOE OSTI2023

Composite Lithium Metal Structure to Mitigate Pulverization and Enable Long‐Life Batteries

In lithium metal batteries, non‐uniform stripping of lithium results in pit formation, which promotes subsequent non‐uniform, dendritic deposition. This viscous cycle leads to pulverization of lithium which promotes cell shorting or capacity degradation, symptoms further exaggerated by high electrode areal loading and lean electrolytes. Here, to address this challenge, a composite lithium metal anode is engineered that contains uniformly distributed, nanometer‐sized carbon particles. This composite lithium is shown to strip more uniformly since the growth of non‐uniform pits is intercepted by the carbon particles. This mechanism is corroborated by a continuum electrochemical model. Subsequent lithium deposition on carbon particles is also found to be more uniform than on the surface with irregular pits. Notably, the pulverization rate of composite lithium is 26 times slower than that of commercial lithium. Moreover, in a Li‐S battery with sulfurized polyacrylonitrile cathode, the use of the composite anode extends the cycle life by three times when the areal capacity is 8 mAh cm −2 . The approach of using an engineered lithium composite structure to address challenges during both stripping and plating can inform future designs of lithium metal anodes for high areal capacity operations.

high areal capacity↗
DOE OSTI2024

Thermally unstable roosts influence winter torpor patterns in a threatened bat species

Abstract Many hibernating bats in thermally stable, subterranean roosts have experienced precipitous declines from white-nose syndrome (WNS). However, some WNS-affected species also use thermally unstable roosts during winter that may impact their torpor patterns and WNS susceptibility. From November to March 2017–19, we used temperature-sensitive transmitters to document winter torpor patterns of tricolored bats (Perimyotis subflavus) using thermally unstable roosts in the upper Coastal Plain of South Carolina. Daily mean roost temperature was 12.9 ± 4.9°C SD in bridges and 11.0 ± 4.6°C in accessible cavities with daily fluctuations of 4.8 ± 2°C in bridges and 4.0 ± 1.9°C in accessible cavities and maximum fluctuations of 13.8 and 10.5°C, respectively. Mean torpor bout duration was 2.7 ± 2.8 days and was negatively related to ambient temperature and positively related to precipitation. Bats maintained non-random arousal patterns focused near dusk and were active on 33.6% of tracked days. Fifty-one percent of arousals contained passive rewarming. Normothermic bout duration, general activity and activity away from the roost were positively related to ambient temperature, and activity away from the roost was negatively related to barometric pressure. Our results suggest ambient weather conditions influence winter torpor patterns of tricolored bats using thermally unstable roosts. Short torpor bout durations and potential nighttime foraging during winter by tricolored bats in thermally unstable roosts contrasts with behaviors of tricolored bats in thermally stable roosts. Therefore, tricolored bat using thermally unstable roosts may be less susceptible to WNS. More broadly, these results highlight the importance of understanding the effect of roost thermal stability on winter torpor patterns and the physiological flexibility of broadly distributed hibernating species.

Biodiversity & Conservation↗
DOE OSTI2023

CFD Simulation of the Dosing Behavior within the Atomic Layer Deposition Feeding System

The effective operation of atomic layer deposition (ALD) feeding system is the premise of realizing specific ALD processes. In the present work, a detailed computational fluid dynamics (CFD) model of the feeding system has been developed and validated, which accounts for the roles of ALD valves and manifolds. A numerical simulation of the compressible fluid flow and heat/mass transfer within the feeding system was conducted. The dosing amounts and the spatiotemporal distributions of the precursors can be accurately predicted using the CFD model, as validated by experimental results. Different precursors, operating conditions, and structures of the feeding system were simulated and analyzed to examine the operating flexibility of the feeding system. The simulation results can be adopted as the upstream boundary conditions for simulations of the ALD process in the reaction chamber. The substrate-scale simulation indicates that the effect of the feeding system on the film deposition is highly related to the surface kinetics of ALD. The present work can serve as a guide for the development and optimization of different ALD-based processes via proper operation and even the design of the feeding system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

First Wall Design of a Tokamak Pilot Plant Using a Monte Carlo Model for 3-D Heat Flux Deposition

We present a method for calculating the heat fluxes deposited on nonaxisymmetric tokamak first wall components, allowing for a first-of-its-kind model for power handling in the tokamak far scrape-off layer (SOL). The DIV3D Monte Carlo model features strict global power conservation and can calculate the finite cross-field plasma transport into magnetically-shadowed regions, which is significant when dealing with meter-scale shadows introduced by components such as poloidal limiters or antennas. As a case study, we apply the DIV3D model to inform the distribution of first wall poloidal limiters in an ARC-class reactor device. We demonstrate that discrete protection limiters can efficiently reduce peak heat fluxes on recessed breeder wall components in the presence of significant far-SOL plasma fluxes. By varying the toroidal periodicity and radial standoff depth of the limiters, we demonstrate one of the tradeoffs that must be considered in first wall design: more limiters provide greater protection, but at the cost of reduced breeding performance. We also present the impact that radial misalignments between limiters would have on first wall power loading.

Monte Carlo methods↗
DOE OSTI2026

Refractory-based thermal energy storage for industrial process heat: one-dimensional modeling, control, and optimization

The variable and weather-dependent output of wind and solar power plants present a substantial challenge for planning and operating electricity-systems, particularly in the absence of cost-effective and dispatchable energy storage technologies. This study investigates a high-temperature, electrically heated, refractory-based thermal energy storage (RTES) system that stores electrical energy as sensible heat in dense ceramic bricks over the 950–1800 °C range. The stored heat can be discharged as a controlled hot-gas stream for industrial heating, fuel substitution in high-temperature processes, or electricity generation. The main novelty is a comprehensive modelling, control, mapping, and optimization framework that integrates one-dimensional transient gas–solid heat transfer, fan-assisted discharge, bypass-flow regulation, reheating logic, fan-power evaluation, insulation-loss assessment, and genetic-algorithm-based design optimization. The model uses feedback from outlet temperature and delivered power to regulate discharge, while a two-stage genetic algorithm optimizes brick-channel geometry, gas-flow operation, and multilayer insulation thicknesses. Storage capacities below 50 MWh and discharge powers of 5–30 MW are analyzed to evaluate hold time, thermal delivery, fan-power penalty, heat loss, state-of-charge evolution, and indicative capital cost. Results demonstrate that optimized and well-insulated refractory-based thermal energy storage units can provide stable, efficient, and repeatable heat delivery over multiple discharge cycles. The generated performance and cost maps support modular refractory thermal energy storage as a practical option for large-scale integration of wind and solar generation and for high-temperature industrial process heat.

25 ENERGY STORAGE↗
DOE OSTI2025

Decoupling the capacity fade contributions in polymer electrolyte-based high-voltage solid-state batteries

Polymer electrolyte (PE)-based solid-state batteries (PE-SSBs) made with high-voltage cathodes are known to suffer from severe capacity fade, stemming primarily from the poor oxidative stability of most PEs under high-voltage cycling conditions. PEs also suffer from greater ion-transport limitations compared to liquid or solid electrolytes. However, often, these limitations are collectively stated to be responsible for the observed capacity fade, and it is challenging to decouple the contributions of different factors. Herein, a tunable cell fabrication platform was developed to systematically investigate and decouple the two primary capacity fade drivers (cell impedance growth and kinetic limitations), while keeping the other cell parameters constant. Three PE types with distinct transport characteristics were compared. By utilizing a voltage profile analysis method, the contribution of the cell's internal impedance growth was quantitatively decoupled from the kinetic limitations stemming from the high concentration gradient in the polymer catholyte and slow charge transfer reactions. We demonstrate that the high interfacial impedance did not necessarily correlate with the high capacity fade rate. Kinetic limitations that are not reflected by impedance measurements can play a dominant role in causing cumulative capacity decay.

Ock, Ji-young [Oak Ridge National Laboratory (ORNL↗
DOE OSTI2025

Direct observation of distinct bulk and edge nonequilibrium spin accumulation in ultrathin MoTe 2

Low-symmetry two-dimensional (2D) topological materials such as MoTe 2 host efficient charge-to-spin conversion (CSC) mechanisms that can be harnessed for novel electronic and spintronic devices. However, the nature of the various CSC mechanisms and their correlation with underlying crystal symmetries remain unsettled. In this work, we use local spin-sensitive electrochemical potential measurements to directly probe the spatially dependent nonequilibrium spin accumulation in MoTe 2 flakes down to four atomic layers. We are able to clearly disentangle contributions originating from the spin Hall and Rashba-Edelstein effects and uncover an abundance of unconventional spin polarizations that develop uniquely in the sample bulk and edges with decreasing thickness. Using ab-initio calculations, we construct a unified understanding of all the observed CSC components in relation to the material dimensionality and stacking arrangement. Our findings not only illuminate previous CSC results on MoTe 2 but also have important ramifications for future devices that can exploit the local and layer-dependent spin properties of this 2D topological material.

Science & Technology - Other Topics↗
DOE OSTI2026

Multitask graph neural networks for elastoplastic response prediction in dual-phase polycrystals

Microstructure-sensitive prediction of elastoplastic response remains a recurring bottleneck in multiscale damage and fatigue modeling, where large ensembles of statistically distinct polycrystals are required to quantify variability and extreme-value behavior. In this work, we develop a multitask graph neural network (GNN) surrogate that maps dual-phase ferrite–martensite polycrystal microstructures to Statistical Volume Element (SVE)-level elastoplastic Quantities of Interest (QoIs). Each SVE is represented as a grain-adjacency graph, with node features encoding phase, geometry, and crystallographic orientation, and edge features encoding relative misorientation. A message-passing graph convolution generates node embeddings, which are pooled into a graph representation and passed to a multitask regression head that jointly predicts 10 scalar QoIs and vector-valued stress–strain responses in orthogonal loading directions across multiple martensite volume fractions and SVE sizes. Results show high accuracy for scalar QoIs and strong agreement for full stress–strain trajectories, with population envelopes reproducing both median behavior and finite-SVE variability across compositions and partition scales. A unified model trained on pooled volume-fraction data preserves most within-regime accuracy relative to regime-specific models while also capturing the broader cross-regime variation reflected in the pooled test set. Distributional comparisons further demonstrate that the surrogate preserves heterogeneity under SVE partitioning, enabling statistically consistent block-wise random-field construction for mesoscale analyses. Overall, the proposed grain-graph surrogate provides a practical pathway to accelerate ensemble-based studies of SVE-level constitutive variability in dual-phase polycrystals.

Crystal plasticity↗