Search NASA⌕ Search

SEARCH · Search NASA

Results for “proxy”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Implementing a neural network interatomic model with performance portability for emerging exascale architectures

The two main thrusts of computational science are increasingly accurate predictions and faster calculations; to this end, the zeitgeist in molecular dynamics (MD) simulations is pursuing machine learned and data driven interatomic models, e.g. neural network potentials, and novel hardware architectures, e.g. GPUs. Current implementations of neural network potentials are orders of magnitude slower than traditional interatomic models and while looming exascale computing offers the ability to run large, accurate simulations with these models, achieving portable performance for MD with new and varied exascale hardware requires rethinking traditional algorithms, using novel data structures, and library solutions. We re-implement a neural network interatomic model in CabanaMD, an MD proxy application, built on libraries developed for performance portability. Our implementation shows significantly improved thread scaling in this complex kernel as compared to a current LAMMPS implementation, across both strong and weak scaling. Our single-source solution enables simulations up to 20 million atoms on a single CPU node and 4 million atoms with improved performance on a single GPU. Furthermore, we also explore parallelism and data layout choices (using flexible data structures called AoSoAs) and their effect on performance, seeing up to ~50% and ~5% improvements in performance on a GPU by choosing the right level of parallelism and data layout respectively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Monitoring and modeling hydrologic conditions in Ukraine for hydropower generation

Study region: The Dnieper and Dniester Rivers of Ukraine. Study focus: The ongoing conflict in Ukraine has caused disruptions to electricity generation, of which hydroelectric sources contribute approximately 9 % to the country’s needs. With the takeover of the Zaporizhzhia nuclear power plant by enemy forces, the loss of the Kakhovka hydroelectric dam, and the future impacts of the conflict on electricity generation unclear, it may be valuable for the Ukrainian government to better understand how it could leverage hydroelectric power sources in the near future. Unfortunately, measurements of river discharge throughout Ukraine ceased data collection in the late 1980’s to early 1990’s. To address this data gap, we developed a protocol that combined satellite-based time-series measurements of river width at seven locations throughout Ukraine from 2013 to 2023 with reanalysis data, climate-model predictions, and hydrologic models to both provide a means of monitoring a proxy for near-real-time discharge and also predict near-term (i.e., 2023–2030) hydrologic patterns for the region. New hydrological insights for the region: We ran new algorithms on 144 WorldView-2 and WorldView-3 satellite images to map rivers and extract width, one of which was validated against river gauge data located along the same river but in a neighboring country. Hydrologic models using two climate scenarios found minimal change in annual discharge at all sites, but magnitude and timing of peak discharge showed a moderate trend. The results suggest that hydropower is underutilized in Ukraine.

13 HYDRO ENERGY↗
DOE OSTI2020

Spatio-temporal impacts of a utility’s efficiency portfolio on the distribution grid

Energy Efficiency has historically focused on delivering savings to offset growth in energy supply. Today’s growing emphasis on decarbonization of the energy supply is driving renewables adoption and increased interest in electrification. As a result, energy efficiency is being assessed not just in its ability to offset load growth, but also for its ability to alleviate location-specific constraints on transmission and distribution infrastructure. This work demonstrates that advanced measurement and verification modeling techniques can be used to estimate the spatio-temporal grid impact of a portfolio of energy efficiency programs. It extends measurement-based methods to an entire Demand Side Management portfolio and uses a single model to predict annual as well as seasonal building energy use with near-zero bias. In addition, new metrics are introduced to assess grid level impacts of energy efficiency. The results show that the efficiency program portfolio delivers savings of over 12% at the territory-wide proxy level, with substation and feeder level savings ranging from 0.4% to 26%, and ~5%-42% respectively. Overall, these savings impacted 1.0%–1.4% of the energy used at these locations in the grid. This work provides a methodology with potential to connect efficiency with distribution planning, carrying implications for non-wires alternatives and targeted delivery of efficiency programs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Impact of operating reserve rules on electricity prices with high penetrations of renewable energy

In competitive wholesale electricity markets, significant effort is devoted to designing markets that set efficient prices for maintaining supply-demand balance. One factor that can impact prices is administratively-set scarcity pricing, which sets prices to a preset level when the market is not able to meet operating reserve or energy requirements. When energy and operating reserves are co-optimized, assumptions surrounding operating reserve requirements and scarcity pricing can impact system-wide price outcomes for both operating reserves and energy. This study uses production cost modeling of an ERCOT-like system to evaluate the impact of operating reserve eligibility, scarcity pricing, and quantity rules on electricity prices, and therefore also on generator revenues. Results reveal economic and operational benefits with allowing open participation in reserve markets, as this enables greater access to the full set of capable resources at lowest cost. Furthermore, both energy and reserve prices are strongly impacted by reserve scarcity pricing events, which reveals that reserve scarcity pricing assumptions can impact price outcomes even for units not providing reserves. This study highlights the importance of operating reserve scarcity pricing rules because of the strong coupling between energy and reserve prices and because these rules serve as proxies for true price responsive demand.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

No detectable redox exchange between sulfur and iron during rapid cooling of basalts

The valence states of iron and sulfur in mantle-derived melts influence and respond to the composition of the mantle as well as the subsequent evolution of those melts in Earth's crust. Characterizing the valence states of iron and sulfur in silicate melts that have quenched to form glasses is therefore critical to understanding the composition of the mantle and the magmatic processes that take place before and during volcanic eruptions. Glasses quenched from silicate melts offer the opportunity to precisely measure sulfur and iron valence states, but whether the transition from melt to glass affects the valence state of sulfur and iron is uncertain. Here, we use the glassy margins of two pillow basalts, one from a mid-ocean ridge and one from a back-arc basin, as a natural experiment to test the effect of quench rate on sulfur and iron valence states. We report micro-X-ray Absorption Near Edge Structure (XANES) measurements of S 6+ /ΣS (S 6+ /[S 6+ +S 2- ]) and Fe 3+ /ΣFe (Fe 3+ /[Fe 2+ +Fe 3+ ]) along transects from the rapidly quenched rims of the pillows down into the slowly cooled crystal-rich interiors. The range of average quench rates estimated for our selected analysis areas range from 1.2 x 10 5 °C/sec at the pillow margin to 1.1 °C/sec in the interior and span quench rates experienced by natural and experimental glasses. On average, S6+/ΣS = 0.09 ± 0.01 and Fe 3+ /ΣFe = 0.147 ± 0.002 in the mid-ocean ridge sample and S 6+ /ΣS = 0.21 ± 0.01and Fe 3+ /ΣFe = 0.153 ± 0.004 in the back-arc sample, where the quoted uncertainty is the 1σ standard deviation of n= 9-16 analyses. We find that the S 6+ /ΣS and Fe 3+ /ΣFe values along transects in each sample show no systematic changes with depth, and by proxy, no systematic changes with quench rate. In each pillow, all measured S 6+ /ΣS and Fe 3+ /ΣFe ratios are within the absolute uncertainty of each technique. We propose that these observations are consistent with either sluggish reaction kinetics between sulfur and iron, or an insignificant temperature effect on sulfur-iron redox exchange equilibria. These observations demonstrate that sulfur and iron valence states in glasses with varying quench rates can be directly compared, and that changes in sulfur and iron valence during quench are below the limits of detection with XANES.

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

On the feasibility of using physics-informed machine learning for underground reservoir pressure management

In this work, we evaluate the feasibility of using physics-informed machine learning (PIML) for underground energy-related pressure management. To this end, we develop a PIML framework to manage underground reservoir pressures by training neural networks to determine fluid extraction rates for dedicated extraction wells during fluid injection operations given a range of reservoir conditions (e.g., transmissivity and storativity). We implement an automatically-differentiable analytical physics model of fluid flow in porous media within the PIML framework as a proxy for more complicated models. This allows us to execute a sufficient number of training scenarios to fully evaluate the feasibility of using PIML to support pressure management activities. We quantify the number of physics-model parameters required for automatic differentiation to become more efficient than finite-difference gradient calculations. We use a simple scenario with a single injector, extractor, and critical location for our feasibility analysis. We evaluate the effect of the size of the training dataset (i.e., the number of reservoir condition samples) on the accuracy and efficiency of the PIML framework. For an equivalent number of model evaluations, the larger training dataset took less time to train and produced a neural network that was able to more accurately manage reservoir pressures. We also evaluate the effect of the training dataset batch size (i.e., number of reservoir condition samples used to update the neural network coefficients during training; i.e., how the training dataset is partitioned). While training ran faster with larger batch sizes, they produced neural networks that managed pressures less accurately. We demonstrate the approach on a more complex scenario involving 10 injectors, 10 extractors, and 4 critical locations (a relatively high well density of 20 wells/km2). We provide the number of forward and adjoint model evaluations required in each case as an indication of the feasibility of using PIML for pressure management when more complicated physics models with longer execution times are used.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Optical diagnostic design for measuring the radiation front with a mid-leg pumped divertor on DIII-D

The design of an optical diagnostic system to localize the radiation front in a mid-leg pumped divertor configuration on DIII-D is presented. Divertor detachment is a key mechanism for handling power exhaust in tokamaks, and mid-leg pumping offers a promising approach to radiatively dissipate power while maintaining acceptable core performance. To predict the location of the radiation front and design a spectroscopic diagnostic to measure its position, a database of SOLPS-ITER simulations across a range of input powers and gas puffing rates representative of DIII-D operation was employed. These simulations provide self-consistent plasma backgrounds for a Cherab-Raysect synthetic diagnostic framework, which incorporates detailed tokamak geometry and physically accurate ray-tracing. Synthetic simulations of optical sightlines and viewing cone geometry were used to calculate line-integrated emission from the plasma, particularly of the C III 465 nm line, to serve as a proxy for the T e ≈ 7 − 10 eV temperature region associated with the onset and evolution of divertor detachment. Neutral deuterium emission is also evaluated for comparison. The synthetic diagnostic is used to assess the performance of existing DIII-D optical systems, including filterscopes and the Multichord Divertor Spectrometer, and to optimize line-of-sight placement within mechanical and installation constraints. The results provide quantitative guidance for diagnostic implementation on DIII-D and demonstrate the advantage of integrated synthetic diagnostics for divertor design studies and future advanced divertor concepts.

Cherab↗
DOE OSTI2020

Deposition of sulfate aerosols with positive Δ 33 S in the Neoarchean

Anomalous sulfur isotope compositions present in Archean rocks have been intensely scrutinized over the last 20 years because they record key aspects of Earth's atmospheric composition prior to the appearance of free molecular oxygen ca. 2.3 billion years ago. These isotopic compositions can be described as mass anomalous fractionations (MAF) and are produced in the atmosphere as UV light interacts with SO 2 molecules. Most interpretations suggest that atmospheric processes generate a reduced S-phase with a positive ( 33 S-enriched) MAF signature, as measured in pyrites, and an oxidized S-phase with a negative anomaly, as measured in bedded barite deposits. However, recent data for carbonate-associated sulfate (CAS) — a direct proxy for the isotopic composition of sulfur from seawater sulfate — in Neoarchean rocks showed no such negative values, but rather the opposite. To understand if the positive MAF anomalies we measured in Neoarchean CAS reflect secondary processes (diagenetic, metamorphic, handling) instead of original signals of Archean seawater sulfate, we collected additional sample suites with various degrees of preservation and metamorphic alteration across the Campbellrand-Malmani platform in South Africa. Results illustrate that within this comprehensive suite, less-altered samples all contain positive MAF values while secondary processes tend to either remove CAS from the sample and/or decrease the 33 S-enrichment. This positive MAF signal in sulfate is therefore reasonably interpreted as a primary depositional origin, and implies that the assumption that sulfate always carries a negative MAF anomaly throughout the Archean rock record needs to be reconsidered. Our CAS observations suggest that future experiments and calculations should also consider atmospheric and/or sulfur cycling processes that can produce oxidized sulfur with a positive MAF signature.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Past fires and post-fire impacts reconstructed from a southwest Australian stalagmite

Speleothem records of past environmental change provide an important opportunity to explore fire frequency and intensity in the past, and the antecedent climatic conditions leading to fire events. Here, fire sensitive geochemical signals in a stalagmite from Yonderup Cave, a shallow cave in Western Australia, are compared to well-documented wildfire events that occurred in recent decades. The results are extended to identify wildfires during the growth interval of the speleothem (1760 CE–2005 CE). Principal component analysis of the stalagmite time-series revealed distinct peaks in a combination of phosphorus and metal (aluminium, zinc, copper and lead) concentrations in response to known fire events, which are interpreted to have come from ash. Varying responses in the geochemical signal in the pre- and post-European colonisation period are likely linked to changes in land management, fire frequency and fire intensities. Of note is what we infer to be a particularly intense fire event concluded to have occurred in 1897 ± 5 CE, which climate sensitive proxies (δ18O and bedrock-derived elements) indicate was preceded by a multi-decadal dry period that began in the late 1860s. The intensity of the deduced fire event is supported by a peak in P that is 6.3 times higher than elsewhere in the speleothem, accompanied by a peak in Zn. It was also inferred that changes to the surface-cave hydrology occurred as a result of this event, indicated by increased input of colloidal organic material onto the stalagmite that was interpreted to be caused by shallow karst bedrock fracturing from the heat-induced deformation that is observed after intense fires. These findings show the potential for speleothems to provide dated records of fire intensity and recurrence intervals. Further development could lead to a better understanding of the climate-fire relationship and the effects of land-management practices on wildfire frequency and intensity.

58 GEOSCIENCES↗
DOE OSTI2022

Mineralogical, magnetic and geochemical data constrain the pathways and extent of weathering of mineralized sedimentary rocks

The oxidative weathering of sulfidic rock can profoundly impact watersheds through the resulting export of acidity and metals. Weathering leaves a record of mineral transformation, particularly involving minor redox-sensitive phases, that can inform the development of conceptual and quantitative models. In sulfidic sedimentary rocks, however, variations in depositional history, diagenesis and mineralization can change or overprint the distributions of these trace minerals, complicating the interpretation of weathering signatures. Here we show that a combination of bulk mineralogical and geochemical techniques, micrometer-resolution X-ray fluorescence microprobe analysis and rock magnetic measurements, applied to drill core samples and single weathered fractures, can provide data that enable the development of a geochemically consistent weathering model. This work focused on one watershed in the Upper Colorado River Basin sitting within the Mesaverde Formation, a sedimentary sandstone bedrock with disseminated sulfide minerals, including pyrite and sphalerite, that were introduced during diagenesis and subsequent magmatic-hydrothermal mineralization. Combined analytical methods revealed the pathways of iron (Fe), carbonate and silicate mineral weathering and showed how pH controls element retention or release from the actively weathering fractured sandstone. Drill core logging, whole rock X-ray diffraction, and geochemical measurements document the progression from unweathered rock at depth to weathered rock at the surface. X-ray microprobe analyses of a 1-cm size weathering profile along a fracture surface are consistent with the mobilization of Fe(II) and Fe(III) into acidic pore water from the dissolution of primary pyrite, Fe-sphalerite, chlorite, and minor siderite and pyrrhotite. These reactions are followed by the precipitation of secondary minerals such as of goethite and jarosite, a Fe-(oxyhydr)oxide and hydrous Fe(III) sulfate, respectively. Microscale analyses also helped explain the weathering reactions responsible for the mineralogical transformations observed in the top and most weathered section of the drill core. For example, dissolution of feldspar and chlorite neutralizes the acidity generated by Fe and sulfide mineral oxidation, oversaturating the solution in both Fe-oxides. The combination of X-ray spectromicroscopy and magnetic measurements show that the Fe(III) product is goethite, mainly present either as a coating on fracture surfaces in the actively weathering region of the core or more homogeneously contained within the unconsolidated regolith at the top of the core. Low-temperature magnetic data reveal the presence of ferromagnetic Fe-sulfide pyrrhotite that, although it occurs at trace concentrations, could provide a qualitative proxy for unweathered sulfide minerals because the loss of pyrrhotite is associated with the onset of oxidative weathering. Pyrrhotite loss and goethite formation are detectable through room-temperature magnetic coercivity changes, suggesting that rock magnetic measurements can determine weathering intensity in rock samples at many scales. In conclusion, this work contributes evidence that the weathering of sulfidic sedimentary rocks follows a geochemical pattern in which the abundance of sulfide minerals controls the generation of acidity and dissolved elements, and the pH-dependent mobility of these elements controls their export to the ground- and surface-water.

58 GEOSCIENCES↗
DOE OSTI2022

Dissolved pyrogenic carbon leaching in soil: Effects of soil depth and pyrolysis temperature

We report Pyrogenic organic carbon (PyC) is found in soils as a heterogeneous mixture of thermally altered plant residues that range in their susceptibility to losses via mineralization, leaching, and erosion. Leaching of PyC as DOC (DPyC) within the soil profile is likely influenced by the chemical composition of solid PyC and the dominant soil processes and properties at a particular depth. Here we report the results of a 2-year laboratory decomposition and leaching study designed to investigate the interactive effects of pyrolysis temperature levels (no pyrolysis, 300 °C, and 450 °C) and soil depth (surface soil, 0–10 cm depth; subsurface soil, 50–70 cm depth) on losses of 13 C-labeled jack pine (Pinus banksiana) wood, wood pyrolyzed at 300 °C (PyC300) and wood pyrolyzed at 450 °C (PyC450) as dissolved organic carbon (DOC). Losses of wood and PyC in the 13 C- DOC pool were measured in leachates drained from soils once a month over the course of 1 year, and at one single leaching event after 2 years. We found that pyrolysis temperature levels interacted with time and soil depth to affect losses of wood and PyC as DOC throughout the 1-year incubation and leaching study, with greater DOC losses from wood in surface soil (0.73 ± 0.076% of added C) than in subsurface soil (0.40 ± 0.063% of added C) averaged across sampling time. Monthly DOC losses from PyC were not affected by soil depth. Cumulative data indicated a small contribution of wood (2.86 ± 0.07%), PyC300 (0.40 ± 0.04%), and PyC450 (0.16 ± 0.01%) to total DOC leached from soils. DOC losses from wood and PyC300 (as a proportion of added C) were greater in surface soil than in subsurface soil, whereas DOC losses from PyC450 were unaffected by differences between surface and subsurface soil. Losses of DOC were greater from PyC300 than from PyC450 in surface soil, with no significant differences in DOC losses between PyC300 and PyC450 in subsurface soil. One single leaching event at the end of the 2-year decomposition study resulted in higher DOC losses from wood in subsurface soil than in surface soil likely due to desorption, and no differences in DOC losses between PyC300 and PyC450 regardless of soil depth. Our results suggest strong interactions between the initial physicochemical composition of organic C inputs and soil properties (soil depth as a proxy) that control the mobility and transport of PyC and should be better represented in Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

"Hidden" hydrothermal technical potential & technoeconomics: Revealing permeability & fluids with more data

Historical hydrothermal estimates have largely relied on temperature or heat flow estimates ignoring the need for natural flowing fluids. More accurate hydrothermal estimates require some indication of permeability and fluids that naturally exist in the subsurface. This paper describes a novel approach that includes proxies of permeability and fluids in hydrothermal estimates by leveraging the relatively data-rich Great Basin. Specifically, nameplate capacities (megawatts) of operating geothermal plants, negative (0 megawatt) locations and 48 geophysical and geologic features are used to used in eXtreme Gradient Boosting (XGBoost) regression to make hydrothermal capacity predictions. Additionally, this work inputs the XGBoost-based hydrothermal predictions into the Renewable Energy Potential (reV) model to quantify technical capacity, its uncertainty and techno-economics. Compared to historical hydrothermal estimates, these predictions adhere to the 37 operating geothermal plants and negative locations. We present a method for subsampling the negative sites to bring the labels into balance that uses the geologic domain knowledge to proportionally represent negatives. Overall, the distributions of the hydrothermal technical capacity and the site levelized cost of energy are respectively much tighter, lower and more accurate than the previous estimates for the Great Basin, as they include geological and geophysical surrogates for permeability and fluids. Percentile (50th and 90th, median and high estimate, respectively) models provide bookends for these metrics.

13 HYDRO ENERGY↗
DOE OSTI2022

The effect of original and initial saturation on residual nonwetting phase capillary trapping efficiency

Injection of supercritical carbon dioxide (CO 2 ) into geological formations is a strategy for both atmospheric greenhouse gas reduction (climate change mitigation) and enhanced oil recovery. To understand CO 2 trapping efficiency, the capillary trapping behaviors that immobilize subsurface fluids were analyzed at the pore-scale using pairs of proxy fluids representing the range of in situ nonwetting and wetting fluid properties encountered in geologic storage reservoirs. The pairs of fluids were cycled through imbibition and drainage processes using a flow cell apparatus containing a sintered glass bead column. Computed x-ray microtomography (microCT) was used to identify immobilized nonwetting fluid after imbibition and drainage events. From microCT images, the trapped residual (post-secondary imbibition) nonwetting phase was spatially correlated to both the original (post-primary imbibition) and the initial (post-primary drainage) nonwetting phase; this relationship is referred to here as the original saturation dependence (SO-dependence) and initial saturation dependence (SI-dependence), respectively. Significant trends of decreasing SO- and SI-dependence with increasing wetting and nonwetting fluid phase viscosities were observed. This finding implies that the amount of CO 2 injected and ultimately trapped is dependent on the nonwetting phase (e.g. oil or gas) already present in the formation, as well as on the manner in which supercritical CO 2 is initially injected, which are factors not considered in current trapping models. To explore the potential effect that varying viscosity and interfacial tension (IFT) might have on this process, we report on a variety of fluid pairs with different viscosities and IFT.

42 ENGINEERING↗
DOE OSTI2026

A trait syndrome ties cell morphology to glycolysis across the yeast subphylum

Co-variation of traits provides fundamental insights into constraints governing their evolution. An inverse correlation between glucose uptake rates (GURs) and cell surface area-to-volume (SA:V) ratios across 11 yeast species was recently reported. Here, we expand the analysis to 282 species to assess the generalizability of this correlation across Saccharomycotina yeasts and the contribution of shared evolutionary history to the co-variation of these traits. Using phylogenetic regression models, we found extracellular acidification rates (ECARs, used as a proxy for GURs) were weakly, but significantly, correlated with SA:V across Saccharomycotina. ECARs were also correlated with genome sizes and growth rates. Our findings support the reported correlation between GURs and SA:V ratios, but suggest other associated traits, including genome sizes. Specifically, yeasts that consume glucose faster tend to have lower SA:V, faster growth rates, and larger genomes, suggesting a trait syndrome governing several metabolic, genomic, and morphological traits across the yeast subphylum.

biological sciences↗
DOE OSTI2023

Diurnal variation of the fair-weather atmospheric electric field in Binchuan, China

With the purpose of understanding the variation of fair-weather atmospheric electric field, the vertical component of atmospheric electric field (recorded as Ez) has been measured in the town Binchuan (25°50'N, 100°34'E, 1438.7 m asl) of Dali city, Yunnan, China from June 14, 2017 to December 1, 2020. Automatic measured visual range is used as a proxy of aerosol concentration. The mean diurnal variation curve of the fair-weather Ez in Binchuan, having two peaks occurred at around 9LT and 20LT, shows little correlation with the well-known Carnegie Curve and other two high-latitude sites' Curves. While strong negative correlation between fair-weather Ez and visual range is found, indicating the variation of fair-weather is dominated by aerosol concentration. The main source of aerosol concentration variation probably sunrise and sunset effect. We support our findings by comparing Ez, visual range with sunrise and sunset times. The peaks (troughs) of Ez (visual range) are proved to have weak correlations with anthropogenic pollution through comparing the Ez, visual range under workdays and holidays.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Overset meshes for incompressible flows: On preserving accuracy of underlying discretizations

This study on overset meshes for incompressible-flow simulations is motivated by accurate prediction of wind farm aerodynamics involving large motions and deformations of components with complex geometry. Specifically, using first-order hyperbolic and elliptic equation proxies for the incompressible Navier-Stokes (NS) equations, we investigate the influence of information exchange between overset meshes on numerical performance where the underlying discretization is second-order accurate. The first aspect of information exchange surrounds interpolation of solution where we examine Lagrange and point-cloud-based interpolation for creating constraint equations between overset meshes. To maintain overall second-order accuracy, higher-order interpolation is required for elliptic problems, but linear interpolation is sufficient for hyperbolic problems in first-order form. Higher-order point-cloud-based interpolation provides a pathway to maintaining accuracy in unstructured meshes, but at higher complexity. The second aspect of information exchange focuses on comparing the approaches of overset single system (OSS) and overset Additive Schwarz (OAS) for coupling the linear systems of the overlapping meshes. While the former involves a single linear system, in the latter the discrete linear systems are solved separately, and solving the global system is accomplished through outer iterations and sequential information exchange in a Jacobi fashion. For the test cases studied, accuracy for hyperbolic problems is maintained by performing two outer iterations, whereas many outer iterations are required for elliptic systems. The order-of-accuracy studies explored here are critical for verifying the overset-mesh coupling algorithms used in engineering simulations. Accuracy of these simulations themselves is, however, quantified using engineering quantities of interest such as drag, power, etc. Consequently, we conclude with numerical experiments using NS equations for incompressible flows where we show that linear interpolation and few outer iterations are sufficient for achieving asymptotic convergence of engineering quantities of interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A family of independent Variable Eddington Factor methods with efficient preconditioned iterative solvers

We present a family of discretizations for the Variable Eddington Factor (VEF) equations that have high-order accuracy on curved meshes and efficient preconditioned iterative solvers. The VEF discretizations are combined with the Discontinuous Galerkin transport discretization from to form effective high-order, linear transport methods. The VEF discretizations are derived by extending the unified analysis of Discontinuous Galerkin methods for elliptic problems presented by Arnold et al. to the VEF equations. This framework is used to define analogs of the interior penalty, second method of Bassi and Rebay, minimal dissipation local Discontinuous Galerkin, and continuous finite element methods. The analysis of subspace correction preconditioners, which use a continuous operator to iteratively precondition the discontinuous discretization, is extended to the case of the non-symmetric VEF system. Numerical results demonstrate that the VEF discretizations have arbitrary-order accuracy on curved meshes, preserve the thick diffusion limit, and are effective on a proxy problem from thermal radiative transfer in both outer transport iterations and inner preconditioned linear solver iterations. We demonstrate that the VEF solution converges to the S N transport solution as the mesh is refined on both problems with smooth and non-smooth behavior in angle. Parallel performance studies show that the interior penalty VEF discretization's linear solve weak scales out to 1024 processors and strong scales well on a single node. Particular attention is paid to the parallel performance of the VEF algorithm when used in combination with a parallel block Jacobi transport sweep.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A hybrid Monte Carlo-deterministic second moment method with efficient variance reduction

In this work, we present a hybrid method that combines Monte Carlo with deterministic finite element methods to solve a linear Boltzmann transport equation. Our hybrid method runs orders of magnitude faster than Monte Carlo, without sacrificing accuracy, for a proxy problem from radiative transfer that contains both optically-thick and optically-thin material. We believe that this is the first demonstration of a hybrid Second Moment Method in more than one spatial dimension, the first to consider more than one material, and the first to use variance reduction. Our variance reduction approach arises from an asymptotic analysis in which we show that the magnitude of the scattering source grows without bound. We transform the problem to compute the deviation of the radiation intensity from isotropy. The magnitude of the source in the transformed problem is bounded, and the quality of the hybrid method solution is dramatically improved by a substantial reduction in the variance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗