Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatial refinement”

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 199 records · Page 11

Seeing the Disturbed Forest for the Trees: Remote Sensing Is Underutilized to Quantify Critical Zone Response to Unprecedented Disturbance

Understanding the severity and extent of near surface critical zone (CZ) disturbances and their ecosystem response is a pressing concern in the face of increasing human and natural disturbances. Predicting disturbance severity and recovery in a changing climate requires comprehensive understanding of ecosystem feedbacks among vegetation and the surrounding environment, including climate, hydrology, geomorphology, and biogeochemistry. Field surveys and satellite remote sensing have limited ability to effectively capture the spatial and temporal variability of disturbance and CZ properties. Technological advances in remote sensing using new sensors and new platforms have improved observations of changes in vegetation canopy structure and productivity; however, integrating measures of forest disturbance from various sensing platforms is complex. By connecting the potential for remote sensing technologies to observe different CZ disturbance vectors, we show that lower severity disturbance and slower vegetation recovery are more difficult to quantify. Case studies in montane forests from the western United States highlight new opportunities, including evaluating post-disturbance forest recovery at multiple scales, shedding light on understory vegetation regrowth, detecting specific physiological responses, and refining ecohydrological modeling. Learning from regional CZ disturbance case studies, we propose future directions to synthesize fragmented findings with (a) new data analysis using new or existing sensors, (b) data fusion across multiple sensors and platforms, (c) increasing the value of ground-based observations, (d) disturbance modeling, and (e) synthesis to improve understanding of disturbance.

54 ENVIRONMENTAL SCIENCES↗

STILGAR End-of-Project Report

The Subsurface Tunnel Imaging LeveraGed by Analysis of Rayleigh wave ellipticity (STILGAR) project demonstrated an integrated geophysical approach for detecting, locating, and characterizing underground structural changes using dense seismic arrays and advanced inversion techniques. Field campaigns were conducted at two operational mines—the Redmond salt mine (Utah) and Graymont Pleasant Gap limestone mine (Pennsylvania)—providing real-world testbeds for monitoring anthropogenic subsurface activity. At the Redmond salt mine, seismic interferometry combined with back-projection inversion successfully identified continuous, low-amplitude signals from mining operations. The approach differentiated stationary from migrating anthropogenic sources, captured daily operational cycles, and validated the potential of passive seismic monitoring for remote detection of underground activity. At the Graymont Pleasant Gap mine, two dense seismic deployments in the spring and fall of 2023 generated over 4 TB of high-resolution data. Key outcomes included the relocation of 199 underground and 8 surface explosions with accuracies within tens of meters and the development of a 3D P-wave velocity model using the triple-difference tomography algorithm (tomoTD) that resolved major structural features such as the mine entrance, low-velocity tunnels, and roof-collapse areas. Ambient noise cross-correlation and back-projection analyses revealed persistent sources linked to ongoing mining activity, whereas horizontal-to-vertical spectral ratio (HVSR) and ellipticity studies confirmed stable site responses across seasons and identified soil thickness trends consistent with regional erosional and depositional processes. Checkerboard and sensitivity tests further validated the robustness of the tomographic results. Overall, the findings emphasize that although significant progress has been made in subsurface imaging, further work is needed to enhance the detection and localization of underground structures. Accurate imaging requires higher frequencies, yet anthropogenic sources tend to dominate the seismic record at those frequencies, and high-frequency surface waves are affected by higher modes that complicate interpretation. The improved detection and localization of human-induced signals enabled detailed temporal and spatial mapping of daily mine operations, demonstrating the feasibility of continuous anthropogenic source monitoring. Sensitivity to signals from nontraditional sources, such as fan operations, highlights the broader applicability of this approach to other industrial environments in which continuous and impulsive signals are present. The field campaigns produced a substantial volume of high-quality seismic data, supporting the development and testing of new methods for seismic source characterization and subsurface imaging. Future deployments should include sensors capable of recording lower frequencies to probe deeper structures, increase bandwidth to enhance resolution and sensitivity to both shallow and deep targets, and collect additional large-scale datasets to refine imaging and source characterization techniques. Moreover, conducting 3D modeling studies of seismic wavefields at higher frequencies will provide a better understanding of wave scattering and cavity–wavefield interactions in complex underground environments. In conclusion, the STILGAR project demonstrated that integrated seismic monitoring can effectively characterize underground operations, capturing both natural and anthropogenic signals. The approaches developed provide a foundation for improved detection, localization, and imaging of subsurface structures and are directly transferable to broader industrial monitoring applications.

58 GEOSCIENCES↗

Refined views of ancient ocean chemistry: Tracking trace element incorporation in pyrite framboids using atom probe tomography

The trace element chemistry of pyrite can be used to determine the origin, timing, and conditions of formation of ore deposits; as a vectoring tool for mineral exploration; and to determine the evolution of the Earth’s oceanic and atmospheric chemistry. However, little is known about whether trace elements are held with the pyrite structure or within nano-inclusions of other phases. This distinction is important for two primary reasons. First, trace element incorporation into the pyrite structure can affect the partitioning of other trace elements. Second, if trace elements are held within nano-inclusions, the partition coefficient of the mineral phases that make up the nano-inclusion, rather than pyrite, are the critical consideration in related interpretations. Previous studies addressing this topic have used laser ablation inductively coupled plasma mass spectrometry; however, the resulting data do not provide sufficient spatial resolution to delineate trace element distributions unless the inclusions are large. Further, they use these flat element profiles in time resolved laser ablation output graphs to argue that pyrite trace element content provides a direct relationship between trace element content of seawater and pyrite trace element content. To improve resolution, we have used atom probe tomography to characterize trace element distributions in pyrite framboids from the Cariaco Basin and Demerara Rise in three-dimensions at sub-nanometer resolution. Manganese was found to be concentrated in the pyrite nanocrystal part of the original framboid structure. In contrast, Ni was mostly found along the grain boundary, though it still appeared to be contained within the pyrite structure. Copper was concentrated in later pyrite overgrowths, and As varied in its location. These observations suggest that some important trace elements are incorporated into pyrite during early diagenesis, even in euxinic settings dominated by water-column pyrite formation. Statistical analysis was used to determine whether trace elements were incorporated in the lattice or within nanoscale inclusions (referred to here are nano-inclusions). We found that As, Ni, Cu, and Mn were commonly held within the pyrite structure, but As, Ni, and Cu can also be held as nano-inclusions or within grain interfaces. Incorporation of As is known to enhance the incorporation of other trace elements and in this case appears to correlated to elevated Ni and Cu concentrations in the Cariaco Basin samples. Furthermore, understanding these relationships strongly impacts our ability to utilize pyrite trace metal concentrations to analyze and quantify early ocean chemistry and its evolution through time.

54 ENVIRONMENTAL SCIENCES↗

Entanglement Hamiltonian of Many-Body Dynamics in Strongly Correlated Systems

A powerful perspective in understanding nonequilibrium quantum dynamics is through the time evolution of its entanglement content. Yet apart from a few guiding principles for the entanglement entropy,to date, much less is known about the refined characteristics of entanglement propagation. In this paper we unveil signatures of the entanglement evolving and information propagating out of equilibrium, from the view of the entanglement Hamiltonian. We investigate quantum quench dynamics of prototypical Bose-Hubbard model using state-of-the-art numerical technique combined with conformal field theory. Before reaching equilibrium, it is found that a current operator emerges in the entanglement Hamiltonian, implying that entanglement spreading is carried by particle flow. In the long-time limit the subsystem enters a steady phase, evidenced by the dynamic convergence of the entanglement Hamiltonian to the expectation of a thermal ensemble. Importantly, the entanglement temperature in steady state is spatially independent,which provides an intuitive trait of equilibrium. These findings not only provide crucial information on how equilibrium statistical mechanics emerges in many-body dynamics, but also add a tool to exploring quantum dynamics from the perspective of the entanglement Hamiltonian.

36 MATERIALS SCIENCE↗

C HIMERA : A Massively Parallel Code for Core-collapse Supernova Simulations

Herein we provide a detailed description of the C HIMERA code, a code developed to model core collapse supernovae (CCSNe) in multiple spatial dimensions. The CCSN explosion mechanism remains the subject of intense research. Progress to date demonstrates that it involves a complex interplay of neutrino production, transport, and interaction in the stellar core, three-dimensional stellar core fluid dynamics and its associated instabilities, nuclear burning, and the fundamental physics of the neutrino–stellar core weak interactions and the equations of state of all stellar core constituents—particularly, the nuclear equation of state associated with core nucleons, both free and bound in nuclei. C HIMERA , by incorporating detailed neutrino transport, realistic neutrino–matter interactions, three-dimensional hydrodynamics, realistic nuclear, leptonic, and photonic equations of state, and a nuclear reaction network, along with other refinements, can be used to study the role of neutrino radiation, hydrodynamic instabilities, and a variety of input physics in the explosion mechanism itself. It can also be used to compute observables such as neutrino signatures, gravitational radiation, and the products of nucleosynthesis associated with CCSNe. The code contains modules for neutrino transport, multidimensional compressible hydrodynamics, nuclear reactions, a variety of neutrino interactions, equations of state, and modules to provide data for post-processing observables such as the products of nucleosynthesis, and gravitational radiation. C HIMERA is an evolving code, being updated periodically with improved input physics and numerical refinements. We detail here the current version of the code, from which future improvements will stem, which can in turn be described as needed in future publications.

79 ASTRONOMY AND ASTROPHYSICS↗

Benchmark microgravity experiments and computations for 3D dendritic-array stability in directional solidification

In this study, we present a comprehensive quantitative analysis of stability bands for dendritic arrays during directional solidification of a transparent succinonitrile-0.46 wt % camphor alloy, spanning a broad range of pulling velocities. Taking advantage of the microgravity environment aboard the International Space Station where most convection effects are suppressed, we obtain unique measurements that quantify the stable primary spacing range of spatially extended three-dimensional dendritic array structures under purely diffusive growth conditions. Through carefully designed velocity jump experiments and detailed examination of sub-grain boundary dynamics, we characterize key instabilities, including elimination and tertiary branching, shedding new light on the mechanisms governing dynamic dendritic spacing selection in extended 3D arrays. Phase field simulations are performed to characterize the stability limits of dendritic array structures for quantitative comparison with the flight experiments. Although the simulations capture general trends, significant deviations are noted at the upper stability boundary, indicating the influence of additional, unexplored factors. These findings contribute to a deeper understanding of dendritic growth dynamics and offer valuable benchmark data that could aid in refining predictive models and improving control of dendritic microstructures in metallurgical applications.

36 MATERIALS SCIENCE↗

Localizing tetrahedral aluminum in nitrate-bearing gibbsite to constrain defect-impurity coupling

The enhanced radiolytic stability of gibbsite (α-Al(OH) 3 ) containing trace nitrate (NO 3 − ) is a phenomenon in nuclear waste management, but its structural origins remain unresolved. Motivated by the detection of minority tetrahedral aluminum (T d ) defects in synthetic gibbsite, we hypothesized that these sites may participate in NO 3 − retention or mediate H 2 suppression. To evaluate this, we combined orthogonal techniques comprised of spatially selective solid-state 27 Al MAS NMR, comparative spectroscopy, and density functional theory (DFT) modeling. Paramagnetic editing and dynamic nuclear polarization (DNP) MAS NMR confirm that T d defects are confined to the particle interior. DFT calculations reveal no energetic stabilization of NO 3 − near T d sites. Comparative NMR analysis shows that T d is also present in chloride-bearing gibbsite, which exhibits high radiolytic hydrogen yields. These three independent disqualifications rule out T d as a structural contributor to nitrate-mediated suppression and narrow the scope of defect-driven explanations. The findings redirect mechanistic attention away from coordination defects and toward redox-active impurity pathways, providing a refined foundation for understanding radiation tolerance in Al(OH) 3 .

Graham, Trent R. [Pacific Northwest National Labor↗

Performance and Stability Characterization of the ICARUS Light Detection System

The ICARUS detector, a key component of the Short Baseline Neutrino (SBN) Program at Fermilab, consists of two identical T300 modules filled with liquid argon. It is equipped with a Light Detection System (LDS) based on 360 8-inch Hamamatsu R5912-MOD photomultiplier tubes (PMTs) arranged behind the wire planes to collect Vacuum Ultraviolet (VUV, $\sim$ 128 nm) scintillation light. Operating under cryogenic conditions ( $\sim$ 87 K), the LDS is essential for determining the event start time (t0) with nanosecond precision for beam spill synchronization, improving longitudinal spatial resolution, and contributing to the event trigger and cosmic-ray mitigation. The performance of the LDS was investigated addressing both hardware and data analysis aspects. Following a progressive degradation in PMT gain observed during operations at FNAL, systematic gain measurements were first carried out from room temperature down to low temperatures. The results show stable performance at room temperature but a significant, irreversible reduction in gain at low temperatures. Based on these findings, a series of mitigation strategies were implemented in the ICARUS detector to preserve PMT performance and ensure reliable cryogenic operation. Currently, ongoing waveform analysis of the PMT signals is being performed to characterize signal shape, charge integration, and timing properties, aiming to refine and improve the agreement between experimental data and Monte Carlo simulations.

Saia, Clara [U. Catania (main); INAF, Catania; Cat↗

Uncertainties in Simulating Flooding During Hurricane Harvey Using 2D Shallow Water Equations

Abstract Flooding is one of the most impactful weather‐related natural hazards. Numerical models that solve the two dimensional (2D) shallow water equations (SWE) represent the first‐principles approach to simulate all types of spatial flooding, such as pluvial, fluvial, and coastal flooding, and their compound dynamics. High spatial resolution (e.g., () m) is needed in 2D SWE simulations to capture flood dynamics accurately, resulting in formidable computational challenges. Thus, relatively coarser spatial resolutions are used for large‐scale simulations of flooding, which introduce uncertainties in the results. It is unclear how the uncertainty associated with the model resolution compares to the uncertainties in precipitation data sets and assumptions regarding boundary conditions when channelized flows interact with other water bodies. In this study, we compare these three sources of uncertainties in 2D SWE simulations for the 2017 Houston flooding event. Our results show that precipitation uncertainty and mesh resolution have more significant impacts on the simulated streamflow and inundation dynamics than the choice of the downstream boundary condition at the watershed outlet. We point out the viability to confine the uncertainty of coarsening mesh resolution by using the variable resolution mesh (VRM) which refines critical topographic features with far fewer grid cells. Specifically, in simulations with VRM, the simulated inundation depths over the refined region are comparable to that use the finest uniform mesh. This study contributes to understanding the challenges and pathways for applying 2D SWE models to improve the realism of flood simulations over large scales.

54 ENVIRONMENTAL SCIENCES↗

Emu v1.1

Emu is a particle-in-cell code for solving the neutrino quantum kinetic equations in 1, 2, or 3 spatial dimensions with arbitrary angular resolution in order to simulate neutrino flavor transformation in neutron star merger and core collapse supernova environments. Emu represents the neutrino distribution function as a set of particles, each of which represent a collection of neutrinos and antineutrinos with unique position and momentum. Each particle carries two density matrices to define the flavor state of the neutrinos and antineutrinos it represents. Emu includes the vacuum, matter, and neutrino self-interaction potentials. Emu calculates the self-interaction potential using PIC deposition and interpolation algorithms that efficiently compute the local number density and flux at particle locations. Emu is implemented in C++ and is based on the AMReX library for high-performance, block-structured adaptive mesh refinement. Emu is parallelized with MPI + OpenMP for CPUs and MPI + CUDA for GPUs.

Willcox, Donald↗

Steady-state mixing state of black carbon aerosols from a particle-resolved model

Abstract. Black carbon (BC) exerts a notable warming effect due to its strong light absorption, largely influenced by its “mixing state”. However, due to computational constraints, the mixing state is challenging to accurately represent in large-scale models. In this study, we employ a particle-resolved model to simulate the evolution of BC mixing state based on field observation. Our result shows that aerosol compositions, coating thickness (CT) distribution, and optical properties of BC aerosols all exhibit a tendency toward a steady state with a characteristic timescale of less than 1 d, considerably shorter than the BC atmospheric lifetime. The rapid attainment of a steady state suggests that it is reasonable to disregard this pre-steady-state period and instead concentrate on the average properties of BC across extensive spatial and temporal scales. The distribution of CT follows an exponential linear distribution and can be characterized by a single slope parameter k. This distribution is independent of the BC core's distribution. In the model simulation, the mean CT, equivalent to the 1/k, is 62 nm, which is consistent with the statistical results indicating a mean CT of 63 nm. Utilizing the slope parameter k, which effectively characterizes the CT distribution under the steady-state simplifying assumption, the BC absorption enhancement closely corresponds to the results obtained via the particle-resolved method. This study simplifies the BC mixing state description and yields a precise evaluation of the BC optical properties, which has the potential utility for modeling efforts in the refinement of the assessment of BC's radiative effects.

Zhang, Zhouyang↗

Deformation mechanisms in an additively manufactured dual-phase eutectic high-entropy alloy

Nanostructured metals and alloys often exhibit high strengths but at the expense of reduced ductility. Through harnessing the far-from-equilibrium processing conditions of laser powder-bed fusion (L-PBF) additive manufacturing, we develop a dual-phase nanolamellar structure comprised of FCC/L1 2 and BCC/B2 phases in a Ni 40 Co 20 Fe 10 Cr 10 Al 18 W 2 eutectic high-entropy alloy (EHEA), which exhibits a combination of ultrahigh yield strength (>1.4 GPa) and large tensile ductility (~17%). The deformation mechanisms of the additively manufactured EHEA are studied via in-situ synchrotron X-ray diffraction and high-resolution transmission electron microscopy. The high yield strength mainly results from effective blockage of dislocation motion by the high density of lamellar interfaces. The refined nanolamellar structures and low stacking fault energy (SFE) promote stacking fault (SF)-mediated deformation in FCC/L1 2 nanolamellae. The accumulation of abundant dislocations and SFs at lamellar interfaces can effectively elevate local stresses to promote dislocation multiplication and martensitic transformation in BCC/B2 nanolamellae. The cooperative deformation of the dual phases, assisted by the semi-coherent lamellar interfaces, gives rise to the large ductility of the as-printed EHEA. In addition, here we also demonstrate that post-printing heat treatment allows us to tune the non-equilibrium microstructures and deformation mechanisms. After annealing, the significantly reduced SFE and thicknesses of the FCC nanolamellae facilitate the formation of massive SFs. The dissolution of nano-precipitates in the BCC/B2 nanolamellae reduces spatial confinement and further promotes martensitic transformation to enhance work hardening. Our work provides fundamental insights into the rich variety of deformation mechanisms underlying the exceptional mechanical properties of the additively manufactured dual-phase nanolamellar EHEAs.

36 MATERIALS SCIENCE↗

Local Resource Assessment: Beetle-Killed Spruce in Alaska

This report synthesizes forest inventory and spatial data from federal and state sources to present a snapshot of beetle-killed spruce across Southcentral and Interior Alaska. This natural resource is under investigation for use in cellulose-based products, particularly building insulation. Research funded by the Department of Energy's ARPA-E program, has successfully developed insulation - Celium - using beetle-killed spruce combined with mycelium, the root network of fungi. The resulting insulation product can be manufactured locally and used across all Alaskan climate zones for residential, commercial, and shipping applications. Removing standing dead wood also contributes to wildfire risk reduction and mitigates associated public health impacts. This Local Resource Assessment, funded by the Denali Commission, contextualizes and reconciles unaligned forestry data sources, establishing a pathway to refined, replicable analysis and, ultimately, presenting estimates of beetle-killed spruce biomass to inform the potential scale and development of a Celium industry.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Modeling wave propagation in elastic solids via high-order accurate implicit-mesh discontinuous Galerkin methods

Here, a high-order accurate implicit-mesh discontinuous Galerkin framework for wave propagation in single-phase and bi-phase solids is presented. The framework belongs to the embedded-boundary techniques and its novelty regards the spatial discretization, which enables boundary and interface conditions to be enforced with high-order accuracy on curved embedded geometries. High-order accuracy is achieved via high-order quadrature rules for implicitly-defined domains and boundaries, whilst a cell-merging strategy addresses the presence of small cut cells. The framework is used to discretize the governing equations of elastodynamics, written using a first-order hyperbolic momentum-strain formulation, and an exact Riemann solver is employed to compute the numerical flux at the interface between dissimilar materials with general anisotropic properties. The space-discretized equations are then advanced in time using explicit high-order Runge–Kutta algorithms. Several two- and three-dimensional numerical tests including dynamic adaptive mesh refinement are presented to demonstrate the high-order accuracy and the capability of the method in the elastodynamic analysis of single- and bi-phases solids containing complex geometries.

42 ENGINEERING↗

2D Magnetohydrodynamic Simulations of the Electrothermal Instability in Metallic Liners

The Virginia Tech (VT) Plasma Dynamics Laboratory Computational (PDCL) and Lawrence Livermore National Laboratory (LLNL) are performing two dimensional (2D) simulations of the electrothermal instability (ETI) using the LLNL multi-physics code Ares. Ares is a multi-physics arbitrary–Lagrangian-Eulerian (ALE) code developed by LLNL and is of particular use in studying magnetohydrodynamic (MHD) instabilities like the ETI due to its resistive MHD, magnetic diffusion, and radiative-hydrodynamics packages. Among its capabilities, it has the ability to model material strength, perform adaptive mesh refinement (AMR), and incorporate a wide variety of equations of state models and conductivity models. The 2D Ares simulation model created by VT-LLNL for studying the development and growth of the electrothermal instability has been configured with initial conditions based on the Mykonos Electrothermal Instability II (METI-II) experiments described by this grant and conducted by team members at the University of Nevada (UNR), the University of New Mexico (UNM), and Sandia National Laboratories. Previously, preliminary 2D Ares simulations of the ETI had been run to approximately 80ns. The rods in these simulations were initiated with sinusoidal perturbations at a similar order of magnitude to those measured on the aluminum rods used for the Mykonos experiment. This model has been improved by increasing the spatial resolution of the simulations and running the simulations further in time. In addition to the simulation run-times extending, the preliminary sinusoidal perturbation has been replaced with a perturbation derived from amplitude measurements by the experimental team, thereby correlating the simulation inputs better to the experimental runs. These new runs are capable of reaching 120ns of simulated time for the uncoated cases and to 200ns the 41 μm coated cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Continuous surface-to-distributed acoustic sensor snapshots explain reactivation of individual natural fractures during an unconventional reservoir stimulation

ABSTRACT Fiber-optic sensing technologies allow petroleum engineering teams to detect hydraulic fracture interaction with boreholes during unconventional reservoir stimulation. In combination with high-repeatability seismic sources, the same distributed acoustic sensors (DASs) enable vertical seismic profiling (VSP) of the fracture evolution away from the boreholes. We discovered clear signatures of seismic scattering on activated fractures during nine days of continuous seismic monitoring of the fracturing stages at the Austin Chalk/Eagle Ford Field Laboratory. The present study applies a novel approach for quantitative analysis of the scattering events in terms of the evolution of the geometry and elastic stiffness of individual fractures. Our characterization strategy sequentially refines the fracture models: from a stack of 1D soft layers to 3D rectangular inclusions. First, we estimate the number of fracture locations and reflectivity using a modified sparse-spike deconvolution of the stacked VSP traces. The fracture set consists of five fractures spaced by 15–30 m with a reflectivity of approximately 1%. Then, we develop a scattering integral method to refine these estimates along with an inversion of the fracture top and bottom for each monitoring vintage. We find that, initially, some of the fractures are located above the monitoring fiber with the height of approximately 100 m. Then we integrate the seismic interpretation with the low-frequency DAS and pressure and microseismic monitoring to reconstruct the activation process of the fractures. Most likely, some of the natural fractures slowly grew downward to the monitoring fiber as a result of fluid injections in the stimulated well. This led to bright strain anomalies but did not trigger seismicity. The top of the fractures remained almost constant and were limited by a lithologic boundary/stress barrier. To our knowledge, this is the first time VSP data enabled tracking of the fracture evolution with such high spatial and temporal resolution, which was previously only available for crosswell surveys and at a much smaller scale.

Glubokovskikh, Stanislav↗

Mapping Competitive Reduction upon Charging in LiNi 0.8 Co 0.15 Al 0.05 O 2 Primary Particles

Side reactions involving surface reduction play a critical role in the failure of LiNi 0.8 Co 0.15 Al 0.05 O 2 to reach its theoretical capacity as a cathode material for Li-ion batteries. While macroscopic consequences are known, the underlying nanoscopic mechanisms are not fully elucidated. By coupling X-ray spectroscopy with several X-ray microscopy modalities, we have spatially resolved the extent of Ni oxidation at several states of charge and uncovered heterogeneity that is hidden when considering ensemble measurements alone. The use of morphologically controlled particles enabled high-resolution imaging of these materials, uncovering gradients of Ni oxidation states within individual primary particles. At high states of charge, these gradients revealed regions of possible oxygen deficiency extending deeper into the particle than previously observed. Surface-sensitive X-ray coupled scanning tunneling microscopy allows oxidation states to be measured at the material's surface, showing predominantly Ni II in the first atomic layer and mixtures of Ni II with Ni III /Ni IV already appearing 1.5 nm into the particle. These results reveal the subtle interplay between irreversible surface transformations and the bulk reactions that ultimately define function, which will refine strategies of surface passivation that are key to overcoming current performance limitations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Approaches for the Simulation of Coupled Processes in Evolving Fractured Porous Media Enabled by Exascale Computing

Models have historically represented fractured porous media with continuum descriptions that characterize the media using bulk parameters. The impact of small-scale features is not captured in these models, although they may be controlling the performance of subsurface applications. Pore-scale models can simulate processes in small-scale features by representing the pore space geometry explicitly but are computationally expensive for large domains. The alternative multiscale approach entails the combination of pore-scale and continuum-scale descriptions in a single framework. We use Chombo-Crunch, a computational capability that discretizes complex geometries with an adaptive, embedded boundary method to contrast these two approaches. Chombo-Crunch takes advantage of recent computational performance and memory bandwidth improvements resulting from the emergence of exascale computing resources. These combined improvements enable the efficient simulation of reactive transport in fractured media with a high degree of fidelity and the ability to capture the control small-scale processes exert on the overall medium evolution.

42 ENGINEERING↗