Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatial interpolation”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Co-Firing Switchgrass and Waste Coal in A Power Plant: A Techno-Economic and Life Cycle Evaluation for The Ohio River Valley (SWITCH) (Final Technical Report for Ohio State/FE0032204)

Abandoned coal mine lands (AMLs) represent one of the most persistent environmental challenges in the United States. Prior to the enactment of the Surface Mining Control and Reclamation Act (SMCRA) in 1977, coal mining operations were not legally required to reclaim disturbed lands, leaving behind approximately 500,000 AML sites nationwide. These sites pose severe environmental and health risks, including acid mine drainage, soil and water contamination, and spontaneous combustion of waste coal piles. Millions of Americans live within one mile of these AMLs, underscoring the urgency of remediation. Traditional reclamation practices, such as planting cool-season grasses, often fail to fully restore ecological function or leverage the economic potential of these lands. This project addressed these challenges by developing integrated strategies for resource recovery, land reclamation, and sustainable energy production. This project evaluated an integrated strategy to convert this liability into an opportunity by recovering waste coal and co-firing it with switchgrass (Panicum virgatum L.) cultivated on reclaimed or marginal AML areas in existing coal-fired power plants. Switchgrass not only provides a renewable feedstock but also aids in land reclamation and carbon sequestration. 1) Remote Sensing and Machine Learning for Waste Coal Identification Using Sentinel-2 satellite imagery and supervised classification, we applied four machine learning models to detect historical waste coal piles. Random Forest achieved the highest accuracy (precision: 86%, recall: 77%). Time-series analysis revealed gradual vegetation recovery since 1986, indicating natural reclamation processes in historical sites, while active mining areas showed ongoing disturbance. This workflow enables scalable monitoring and prioritization of reclamation efforts. 2) UAS-Based Stockpile Volume Estimation To quantify recoverable waste coal, we evaluated Unmanned Aerial Systems (UAS) equipped with Light Detection and Ranging (LiDAR) and multispectral sensors. Structure-from-Motion (SfM) photogrammetry combined with interpolated Digital Terrain Models (DTMs) achieved strong agreement with LiDAR reference volumes (Root Mean Square Error (RMSE) ≈147 m 3 , Mean Absolute Percentage Error (MAPE) ≈2%). Sensitivity analysis confirmed that spatial resolution significantly influences accuracy, emphasizing the need for high-resolution data for precise volume estimation. This approach offers a scalable, cost-effective, and accurate alternative to conventional ground-based surveys. 3) Switchgrass Cultivation for Bioenergy and Water Quality Improvement We assessed the hydrological and water quality impacts of converting AMLs to switchgrass production areas using the Soil and Water Assessment Tool (SWAT). Results showed that converting 10% of the watershed area into the switchgrass production zone reduced streamflow by 3.1%, total suspended solids by 18.1%, total nitrogen by 7.6%, and total phosphorus by 6.2%, while achieving biomass yields of 8.6–9.2 metric tons per hectare. These findings highlight switchgrass as a dual-benefit strategy for land reclamation and bioenergy feedstock production. 4) Integrated Co-Firing and CCS for Carbon-Negative Power Generation We modeled co-firing scenarios using the Power Plant Flexible Model (PPFM) to evaluate plant efficiency, greenhouse gas (GHG) emissions, and levelized cost of electricity (LCOE). Without carbon capture and storage (CCS), increasing switchgrass co-firing ratios reduced LCOE from $\$$150/MWh at 0% biomass to $\$$110/MWh at full substitution. Under CCS, costs remained higher (~$\$$250/MWh at 0% biomass) but decreased to $\$$200/MWh at 100% biomass, while enabling net-zero or carbon-negative electricity due to switchgrass sequestration benefits. Although CCS introduced efficiency penalties, pairing it with biomass co-firing offset these impacts and maximized climate benefits. Overall, optimizing co-firing ratios between 60-100%, supported by reliable logistics and storage strategies, emerged as a practical pathway to balance affordability, sustainability, and net-zero or negative GHG emissions while promoting productive reuse of AMLs.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Beyond Taylor’s hypothesis: a novel volumetric reconstruction of velocity and density fields for variable‑density and shear flows

Abstract This work presents a novel numerical procedure for reconstructing volumetric density and velocity fields from planar laser-induced fluorescence (PLIF) and stereoscopic particle image velocimetry (SPIV) data. This new method is theoretically and practically demonstrated to provide more accurate 3D vortical structures and density fields in high shear flows than reconstruction methods based on the mean convective velocity. While Taylor’s hypothesis of frozen turbulence is commonly applied by using the local mean streamwise velocity, the proposed algorithm uses the measured local instantaneous velocity for data convection. It consists of a step-by-step reconstruction based on a mixed Lagrangian–Eulerian solver that includes the 3D interpolation of scattered flow data and that relaxes the Taylor’s hypothesis by iterative enforcement of the incompressibility constraint on the velocity field. This methodology provides 3D fields with temporal resolution, spatial resolution, and accuracy comparable to that of real 3D snapshots, thus providing a practical alternative to tomographic measurements. The procedure is validated using numerical data of the constant-density channel flow available on the Johns Hopkins University Turbulence Database (JHTDB), showing the accurate reconstruction of the 3D velocity field. The algorithm is applied to an experimental dataset of PLIF and SPIV measurements of a variable-density jet flow, demonstrating its capability to provide 3D velocity and density fields that are more consistent with the Navier–Stokes equations compared to the mean flow convective method. Graphic abstract

42 ENGINEERING↗
DOE OSTI2016

The concentration dependence of the galaxy–halo connection: Modeling assembly bias with abundance matching

Empirical methods for connecting galaxies to their dark matter halos have become essential for interpreting measurements of the spatial statistics of galaxies. In this work, we present a novel approach for parameterizing the degree of concentration dependence in the abundance matching method. Furthermore, this new parameterization provides a smooth interpolation between two commonly used matching proxies: the peak halo mass and the peak halo maximal circular velocity. This parameterization controls the amount of dependence of galaxy luminosity on halo concentration at a fixed halo mass. Effectively this interpolation scheme enables abundance matching models to have adjustable assembly bias in the resulting galaxy catalogs. With the new $400\,\mathrm{Mpc}\,{h}^{-1}$ DarkSky Simulation, whose larger volume provides lower sample variance, we further show that low-redshift two-point clustering and satellite fraction measurements from SDSS can already provide a joint constraint on this concentration dependence and the scatter within the abundance matching framework.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Finite spatial-grid effects in energy-conserving particle-in-cell algorithms

Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. We also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length. Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. Finally, we also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Real-Space Quantification of Exciton Localization in Acene Crystals Using Wannier Function Decomposition

Here, we introduce the Wannier function decomposition of excitons (WFDX) method to quantify exciton localization in solids within the ab initio Bethe-Salpeter equation framework. By decomposing each Bloch exciton wave function into products of single-particle electron and hole maximally localized Wannier functions, this real-space approach provides well-defined orbital- and spatial-resolved measures of both Frenkel and charge-transfer excitons at low computational cost. We apply WFDX to excitons in acene crystals, quantifying how the number of rings, the exciton spin state, and the center-of-mass momentum affect spatial localization. Additionally, we show how this real-space representation reflects structural nonsymmorphic symmetries that are hidden in standard reciprocal-space descriptions. We demonstrate how the WFDX framework can be used to efficiently interpolate exciton expansion coefficients in reciprocal-space and outline how it may facilitate evaluation of observables involving position operators, highlighting its potential as a general tool for both analyzing and computing excitonic properties in solids.

Tao, Zui [University of California, Berkeley, CA (↗
DOE OSTI2025

High-resolution national mapping of natural gas composition substantially updates methane leakage impacts

Methane is emitted from oil and gas operations alongside heavier hydrocarbons and non-hydrocarbon gases, shaping emissions management decision-making, including air quality impacts. Yet, most assessments assume fixed gas composition, overlooking significant spatial and temporal variations. Here, we generate a high-resolution, data-driven map of natural gas composition across the United States, reconstructing methane, heavier hydrocarbons, and non-hydrocarbon species using spatio-temporal interpolation and oil-and-gas production patterns. Our approach is able to reduce composition prediction errors by 39% in terms of Mean Absolute Error (MAE) compared to standard techniques and reveals that methane loss rates have been underestimated by more than 50% in some regions. Beyond methane, we uncover substantial variability in co-emitted gases, exposing blind spots in current emissions inventories and emissions management frameworks. Our work enables more accurate emissions assessments, guides targeted measurement strategies, and informs emissions management decision-making. It also provides a general framework for prediction in environmental applications that integrate sparse measurements with auxiliary variables.

03 NATURAL GAS↗
DOE OSTI2019

BonVoision: Leveraging Spatial Data Smoothness For Recovery From Memory Soft Errors

The increasing soft error rates in memory systems raise an emerging concern for modern computing systems. As a result, detectable but uncorrectable errors (DUEs) become potentially more frequent and affect HPC applications. Today, upon encountering a DUE, HPC applications crash, incurring significant performance, storage, and energy overheads. In this paper, we propose a technique to continue application execution past a DUE through the repair of the corrupted memory data by leveraging spatial data smoothness. We present BonVoision, a run-time system that intercepts DUE events, analyzes the binary to identify data elements in the structural neighborhood of the event, and fixes the corrupted data elements by interpolating from the values in their neighborhood. Our evaluation demonstrates that BonVoision incurs negligible overhead and outperforms other recovery strategies by a factor of 2×, on average. We demonstrate that BonVoision also improves the efficiency of existing checkpointing/restart schemes by approximately increasing the optimal checkpoint interval by 23%.

Fang, Bo↗
DOE OSTI2017

The Curvature-Augmented Closest Point method with vesicle inextensibility application

Here, the Closest Point method, initially developed by Ruuth and Merriman, allows for the numerical solution of surface partial differential equations without the need for a parameterization of the surface itself. Surface quantities are embedded into the surrounding domain by assigning each value at a given spatial location to the corresponding value at the closest point on the surface. This embedding allows for surface derivatives to be replaced by their Cartesian counterparts (e.g. ∇ s =∇). This equivalence is only valid on the surface, and thus, interpolation is used to enforce what is known as the side condition away from the surface. To improve upon the method, this work derives an operator embedding that incorporates curvature information, making it valid in a neighborhood of the surface. With this, direct enforcement of the side condition is no longer needed. Comparisons in R 2 and R 3 show that the resulting Curvature-Augmented Closest Point method has better accuracy and requires less memory, through increased matrix sparsity, than the Closest Point method, while maintaining similar matrix condition numbers. To demonstrate the utility of the method in a physical application, simulations of inextensible, bi-lipid vesicles evolving toward equilibrium shapes are also included.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Finite grid instability and spectral fidelity of the electrostatic Particle-In-Cell algorithm

The origin of the Finite Grid Instability (FGI) is studied by resolving the dynamics in the 1D electrostatic Particle-In-Cell (PIC) model in the spectral domain at the single particle level and at the collective motion level. The spectral fidelity of the PIC model is contrasted with the underlying physical system or the gridless model. The systematic spectral phase and amplitude errors from the charge deposition and field interpolation are quantified for common particle shapes used in the PIC models. Lastly, it is shown through such analysis and in simulations that the lack of spectral fidelity relative to the physical system due to the existence of aliased spatial modes is the major cause of the FGI in the PIC model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

87 Sr/ 86 Sr Concentrations in the Appalachian Basin: A Review

This document reviews 87 Sr/ 86 Sr isotope data across the Appalachian Basin from existing literature to show spatial and temporal variation. Isotope geochemistry presents a means of understanding the geochemical effects hydraulic fracturing may have on shallow ground substrates. Isotope fractionation is a naturally occurring phenomenon brought about by physical, chemical, and biological processes that partition isotopes between substances; therefore, stable isotope geochemistry allows geoscientists to understand several processes that shape the natural world. Strontium isotopes can be used as a tool to answer an array of geological and environmental inquiries. In some cases, strontium isotopes are sensitive to the introduction of a non-native fluid into a system. This ability allows strontium isotopes to serve as tracers in certain systems. Recently, it has been demonstrated that strontium isotopes can serve as a monitoring tool for groundwater and surface water systems that may be affected by hydraulic fracturing fluids (Chapman et al., 2013; Kolesar Kohl et al., 2014). These studies demonstrated that 87 Sr/ 86 Sr values have the potential to monitor subsurface fluid migration in regions where extraction of Marcellus Shale gas is occurring. This document reviews publicly available strontium isotope data from 39 sample locations in the Appalachian Basin (Hamel et al., 2010; Chapman et al., 2012; Osborn et al., 2012; Chapman et al., 2013; Capo et al., 2014; Kolesar Kohl et al., 2014). The data is divided into two sets: stratigraphic (Upper Devonian/Lower Mississippi, Middle Devonian, and Silurian) and groundwater. ArcMap™ (ESRI, Inc.) was used to complete inverse distance weighting (IDW) analyses for each dataset to create interpolated surfaces in an attempt to find regional trends or variations in strontium isotopic values across the Appalachian Basin. 87 Sr/ 86 Sr varies up to ~ 0.011 across the Appalachian Basin, but the current publicly available data is limited in frequency and regional extent, causing artifacts and high uncertainty when interpolating data for locations far from sampling sites. These factors highlight the need for additional strontium isotope sampling across the region. Identifying potential contamination from hydraulic fracturing fluid in Appalachian Basin groundwater using strontium isotopes would require additional sampling. For a more comprehensive strontium isotope database, samples would need to be collected during prefracturing, syn-fracturing, and post-fracturing stages. This would add a temporal component to the spatial data and make tracing of fluid migration with strontium isotopes more accurate. Future research and modeling that incorporates subsurface geology and watershed data would also serve to increase the accuracy and certainty of the interpolations of these analyses. Prospective geospatial Appalachian Basin isotope studies would also benefit from the integration of geologic mapping because surface and subsurface geology influences observed strontium isotope values.

58 GEOSCIENCES↗
DOE OSTI2023

Spatio-Temporal Surrogates for Interaction of a Jet with High Explosives: Part I - Analysis with a Small Sample Size

Computer simulations, especially of complex phenomena, can be expensive, requiring high-performance computing resources. Often, to understand a phenomenon, multiple simulations are run, each with a different set of simulation input parameters. These data are then used to create an interpolant, or surrogate, relating the simulation outputs to the corresponding inputs. When the inputs and outputs are scalars, a simple machine learning model can suffice. However, when the simulation outputs are vector valued, available at locations in two or three spatial dimensions, often with a temporal component, creating a surrogate is more challenging. In this report, we use a two-dimensional problem of a jet interacting with high explosives to understand how we can build high-quality surrogates. The characteristics of our data set are unique - the vector-valued outputs from each simulation are available at over two million spatial locations; each simulation is run for a relatively small number of time steps; the size of the computational domain varies with each simulation; and resource constraints limit the number of simulations we can run. We show how we analyze these extremely large data-sets, set the parameters for the algorithms used in the analysis, and use simple ways to improve the accuracy of the spatio-temporal surrogates without substantially increasing the number of simulations required.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Systematic analysis of melt pool dynamics in laser processing of mixed powder feedstocks

Functionally graded materials (FGMs) fabricated via additive manufacturing of blended powders offer the potential to spatially tailor properties for new technologies, such as fusion first-wall systems, turbine blades, and spacecraft. However, processing these materials is difficult due to the multiplicity of processing parameters to optimize, all of which must be changed as substrate material, powder feedstock compositions, and melt pool dynamics evolve. Here, this work systematically evaluates the qualitative and quantitative effects of these variables on the melt pool size, shape, composition, and particle distribution in an exemplar Ti-Ta system, and connects the experimental results to Marangoni flow behavior and phenomena observed in other systems. Increasing laser power linearly increases melt pool size and layer thickness, driving engineering considerations such as part/geometrical tolerances. Decreasing laser velocity changes the melt pool shape from lenticular to convex and reduces chemical homogeneity due to extreme thermal and compositional gradients between the melt pool center and boundaries. Thermophysical property differences between the powder feedstock and substrate material, as well as the directionality of the gradient, affect dilution and melt pool dynamics, which in turn affect the melt pool boundary characteristics, shape, and uniformity. Mixed powder feedstocks of intermediate compositions do not behave according to linear interpolations between single-material endpoints, instead building taller and wider melt pools. As such, it is recommended to quantify process maps for at least one intermediate composition in the FGM or multi-material system of interest to ensure optimized processing parameters, predictable melt pool sizes and shapes, and compositional and spatial precision.

Dissimilar↗
DOE OSTI2023

Fitting Matérn smoothness parameters using automatic differentiation

The Mat$\acute{e}$rn covariance function is ubiquitous in the application of Gaussian processes to spatial statistics and beyond. Perhaps the most important reason for this is that the smoothness parameter $\nu$ gives complete control over the mean-square differentiability of the process, which has significant implications for the behavior of estimated quantities such as interpolants and forecasts. Unfortunately, derivatives of the Mat$\acute{e}$rn covariance function with respect to $\nu$ require derivatives of the modified second-kind Bessel function $K$ $\nu$ with respect to $\nu$. While closed form expressions of these derivatives do exist, they are prohibitively difficult and expensive to compute. For this reason, many software packages require fixing $\nu$ as opposed to estimating it, and all existing software packages that attempt to offer the functionality of estimating $\nu$ use finite difference estimates for $\partial$ $\nu$ $K$ $\nu$ . In this work, we introduce a new implementation of $K$$\nu$ that has been designed to provide derivatives via automatic differentiation (AD), and whose resulting derivatives are significantly faster and more accurate than those computed using finite differences. Here, we provide comprehensive testing for both speed and accuracy and show that our AD solution can be used to build accurate Hessian matrices for second-order maximum likelihood estimation in settings where Hessians built with finite difference approximations completely fail.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

An a priori evaluation of a principal component and artificial neural network based combustion model in diesel engine conditions

A principal component analysis (PCA) and artificial neural network (ANN) based chemistry tabulation approach is presented. ANNs are used to map the thermochemical state onto a low-dimensional manifold consisting of five control variables that have been identified using PCA. Three canonical configurations are considered to train the PCA-ANN model: a series of homogeneous reactors, a nonpremixed flamelet, and a two-dimensional lifted flame. The performance of the model in predicting the thermochemical manifold of a spatially-developing turbulent jet flame in diesel engine thermochemical conditions is a priori evaluated using direct numerical simulation (DNS) data. The PCA-ANN approach is compared with a conventional tabulation approach (tabulation using ad hoc defined control variables and linear interpolation). The PCA-ANN model provides higher accuracy and requires several orders of magnitude less memory. Here, these observations indicate that the PCA-ANN model is superior for chemistry tabulation, especially for modelling complex chemistries that present multiple combustion modes as observed in diesel combustion. The performance of the PCA-ANN model is then compared to the optimal estimator, i.e. the conditional mean from the DNS. The results indicate that the PCA-ANN model gives high prediction accuracy, comparable to the optimal estimator, especially for major species and the thermophysical properties. Higher errors are observed for the minor species and reaction rate predictions when compared to the optimal estimator. It is shown that the prediction of minor species and reaction rates can be improved by using training data that exhibits a variation of parameters as observed in the turbulent flame. The output of the ANN is analysed to assess mass conservation. It is observed that the ANN incurs a mean absolute error of 0.05% in mass conservation. Furthermore, it is demonstrated that this error can be reduced by modifying the cost function of the ANN to penalise for deviation from mass conservation.

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

MCNP6 Study of Spatial Fidelity Required for ROC Curve Convergence

The detection of special nuclear material at ports and borders is one of the primary missions of the Department of Homeland Security's (DHS) Domestic Nuclear Detection Office (DNDO). In support of this mission a large-scale modeling and simulation effort has been undertaken with the goal of providing a comprehensive tool-kit that addresses the complexity of problems involved with developing and maintaining technologies deployed to global nuclear detection architecture (GNDA) nodes. An important objective of DNDO's mission is to characterize detector system performance before deployment. In the effort to support this mission, the latest version of MCNP has been used to perform detectability studies of highly enriched uranium (HEU) to populate a database that will be used for the tool-kit. The characterization method used in the implementation of the database is the receiver operating characteristic (ROC). A ROC curve shows the limits of detectability of a given system as well as an overall metric of the efficacy of a detection system. This is accomplished by comparing total detector count distributions to the count distribution produced by background alone. For an in-depth explanation of the theory behind the MCNP production of ROC curves please refer to the paper, 'How To Use The Receiver Operating Characteristic Tally Option in MCNP6'. Building a useful database requires an expansive parameter space. The current tool-kit use nine parameters, which if varied with just 5 options each, results in nearly two million simulations to run. Depending on the set of parameters these simulations can take a few hours or a few weeks to complete. This computational cost necessitates the use of scaling algorithms to interpolate between the data contained in the database. The tool-kit employs a discrete event method to simulate the movement and detector response of special nuclear material (SNM) through a GNDA node. Given the nature of a discrete event simulation, constant motion of SNM must be split up in to discrete spatial detection points. The fidelity required of the spatial splitting to achieve a properly converged answer is investigated here for a neutron detection system. The detector used was a SNAP-3 neutron detector in front of a sphere of SNM in the presence of background radiation. A parameter study has been performed to determine the spatial fidelity required for ROC curve convergence for moving objects in a discrete event simulator. It appears that 5 points is the minimum number of detection points with 7 points needed for one of the parameter sets. The convergence worsens with increases to detection length but appears to improve with increasing detector distance. These results will help future work involving quantifying the accuracy based on given tolerances and parameters. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Distributed Fiber Optics Measurements of Rock Deformation and Failure in Triaxial Tests

The implementation of distributed strain measurement methods in triaxial and uniaxial tests have demonstrated the development of strain localization, even at early stages. This implies that single point measurement methods are location dependent. The use of distributed methods is required not only to improve the interpreted constitutive parameters obtained from triaxial tests, but also to understand the implications of strain localization in the failure process. Here, this work uses optical fibers in triaxial tests. The developed distributed measurement method was implemented on granite, gabbro, and sandstone samples and tested under different confining pressures, reaching 200 MPa on the granite sample. Using a temporal resolution of 0.25 Hz and 5 mm of spatial resolution, the strain evolution at over 300 locations at the sample surface was measured during testing. When compared to point sensing methods, the use of optical fiber greatly increases the number of measurements at the surface of the sample, and their interpolation provides the entire deformation of the sample surface. In all the tests performed, strain localizations were revealed before failure. A three-dimensional interpretation of a test, combining an optical scan of the sample and the distributed measurements, show good correlation between the fractures and the strain localization.

58 GEOSCIENCES↗
DOE OSTI2026

Integrated Neutronics Modeling for Inertial Fusion Energy Systems: Development and Application to LD-FIRST

Lawrence Livermore National Laboratory (LLNL) is proposing a new Laser Driven Fusion Integration Research and Science Test Facility (LD-FIRST) with the goal of providing an experimental testbed for future Inertial Fusion Energy (IFE) systems. However, IFE systems require detailed and accurate multiphysics modeling to quantify material damage, thermal loading, and tritium breeding within complex chamber environments. This article presents the first step in an integrated multiphysics framework that couples meshed CAD-based geometry within Monte Carlo neutronic simulations to enable high-fidelity analysis of IFE chamber concepts, with future coupling to external codes. The neutronics workflow utilizes OpenMC and its third-party capability to use CAD-based geometries through DAGMC and tally on unstructured meshes with Libmesh to evaluate neutron transport behavior, geometric fidelity, and material performance under reactor-relevant conditions. The use of tailored tallies on unstructured meshes in this framework allows direct transfer without interpolating to CFD simulation tools. Two IFE chambers were evaluated, both conceived by LLNL: HYLIFE-II and Laser IFE (LIFE). This work produced high-fidelity conformal surface and volumetric meshes of the HYLIFE-II and LIFE chambers with mapped spatial insight into material damage, thermal loading, and tritium breeding. The HYLIFE-II model was built utilizing available resources and used as a test case to verify that the neutronics framework can handle complex geometries. The LIFE chamber CAD was provided by LLNL and was the main focus of this work. This work analyzes multiple ternary alloy breeding materials for the LIFE chamber, across different 6 Li enrichments to produce data relevant to the LD-FIRST project. This work also investigates the level of model fidelity for the LIFE chamber, and results show that inclusion of detailed first wall and coolant structures increased the predicted tritium breeding ratio (TBR) by ~30%, highlighting the sensitivity of tritium breeding and the need for a high-fidelity simulation framework for IFE chambers. These developments provide a scalable toolset for the design and optimization of next-generation IFE chambers, forming a solid foundation for future coupled multiphysics analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗