SEARCH · Search NASA

Results for “parallel state estimation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2023

tih_vqe [SWR-23-32]

This software supports the paper, "Exploring the scaling limitations of the variational quantum eigensolver with the bond dissociation of hydride diatomic molecules," published in the International Journal of Quantum Chemistry, whose abstract is as follows: Materials simulations involving strongly correlated electrons pose fundamental challenges to state-of-the-art electronic structure methods but are hypothesized to be the ideal use case for quantum computing. To date, no quantum computer has simulated a molecule of a size and complexity relevant to real-world applications, despite the fact that the variational quantum eigensolver (VQE) algorithm can predict chemically accurate total energies. Nevertheless, because of the many applications of moderately-sized, strongly correlated systems, such as molecular catalysts, the successful use of the VQE stands as an important waypoint in the advancement toward useful chemical modeling on near-term quantum processors. In this paper, we take a significant step in this direction. We lay out the steps, write, and run parallel code for an (emulated) quantum computer to compute the bond dissociation curves of the TiH, LiH, NaH, and KH diatomic hydride molecules using VQE. TiH was chosen as a relatively simple chemical system that incorporates d orbitals and strong electron correlation. Because current VQE implementations on existing quantum hardware are limited by qubit error rates, the number of qubits available, and the allowable gate depth, recent studies have focused on chemical systems involving s and p block elements. Through VQE + UCCSD calculations of TiH, we evaluate the near-term feasibility of modeling a molecule with d-orbitals on real quantum hardware. We demonstrate that the inclusion of d-orbitals and the use of the UCCSD ansatz, which are both necessary to capture the correct TiH physics, dramatically increase the cost of this problem. We estimate the approximate error rates necessary to model TiH on current quantum computing hardware using VQE+UCCSD and show them to likely be prohibitive until significant improvements in hardware and error correction algorithms are available.

Graf, Peter↗

DOE OSTI2013

Impact of the MLC on the MRI field distortion of a prototype MRI-linac

Purpose: To cope with intrafraction tumor motion, integrated MRI-linac systems for real-time image guidance are currently under development. The multileaf collimator (MLC) is a key component in every state-of-the-art radiotherapy treatment system, allowing for accurate field shaping and tumor tracking. This work quantifies the magnetic impact of a widely used MLC on the MRI field homogeneity for such a modality.Methods: The finite element method was employed to model a MRI-linac assembly comprised of a 1.0 T split-bore MRI magnet and the key ferromagnetic components of a Varian Millennium 120 MLC, namely, the leaves and motors. Full 3D magnetic field maps of the system were generated. From these field maps, the peak-to-peak distortion within the MRI imaging volume was evaluated over a 30 cm diameter sphere volume (DSV) around the isocenter and compared to a maximum preshim inhomogeneity of 300 μT. Five parametric studies were performed: (1) The source-to-isocenter distance (SID) was varied from 100 to 200 cm, to span the range of a compact system to that with lower magnetic coupling. (2) The MLC model was changed from leaves only to leaves with motors, to determine the contribution to the total distortion caused by MLC leaves and motors separately. (3) The system was configured in the inline or perpendicular orientation, i.e., the linac treatment beam was oriented parallel or perpendicular to the magnetic field direction. (4) The treatment field size was varied from 0 × 0 to 20×20 cm{sup 2}, to span the range of clinical treatment fields. (5) The coil currents were scaled linearly to produce magnetic field strengths B{sub 0} of 0.5, 1.0, and 1.5 T, to estimate how the MLC impact changes with B{sub 0}.Results: (1) The MLC-induced MRI field distortion fell continuously with increasing SID. (2) MLC leaves and motors were found to contribute to the distortion in approximately equal measure. (3) Due to faster falloff of the fringe field, the field distortion was generally smaller in the perpendicular beam orientation. The peak-to-peak DSV distortion was below 300 μT at SID≥130 cm (perpendicular) and SID≥140 cm (inline) for the 1.0 T design. (4) The simulation of different treatment fields was identified to cause dynamic changes in the field distribution. However, the estimated residual distortion was below 1.2 mm geometric distortion at SID≥120 cm (perpendicular) and SID≥130 cm (inline) for a 10 mT/m frequency-encoding gradient. (5) Due to magnetic saturation of the MLC materials, the field distortion remained constant at B{sub 0}>1.0 T.Conclusions: This work shows that the MRI field distortions caused by the MLC cannot be ignored and must be thoroughly investigated for any MRI-linac system. The numeric distortion values obtained for our 1.0 T magnet may vary for other magnet designs with substantially different fringe fields, however the concept of modest increases in the SID to reduce the distortion to a shimmable level is generally applicable.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2025

Particle-based modelling of axisymmetric tandem mirror devices

In this work, we describe the use of a 1D-2V quasi-neutral hybrid electrostatic PIC with Monte-Carlo Coulomb collisions and non-uniform magnetic field to model the parallel transport and confinement in an axisymmetric tandem mirror device. End-plugs, based on simple-mirrors, are positioned at each end of the device and fueled with neutral beams (25 and 100 keV) to produce a sloshing ion population and increase the density of the end-plugs relative to the central cell. Results show the formation of a potential difference barrier between the central cell and the end-plugs. This potential confines a large fraction of the low energy thermal ions in the central cell which would otherwise be lost in a simple mirror, demonstrating the advantage of the beam-driven tandem mirror configuration relative to simple mirrors. In addition, we explore the effect of end-plug electron temperature on the confinement time of the device and compare it with theoretical estimates. Finally, we discuss the limitations of the code in its present form and describe the next logical steps to improve its predictive capability such as a fully nonlinear Fokker–Planck collision operator, multiply nested flux surface solutions and modeling the exhaust region up to the wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

Surrogate Models for TREAT Transient Calculations

The Transient Reactor Test Facility (TREAT) at Idaho National Laboratory (INL) is currently preparing for resumption of operation after more than 20 years in a long-term standby state. In parallel with preparation for operations and a new transient testing program, the DOENE Advanced Modeling and Simulation program (NEAMS) has been supporting the development of three-dimensional transient multi-physics modeling capabilities within the reactor physics application MAMMOTH. Recently, the MAMMOTH reactor physics team has begun development of uncertainty and sensitivity analysis for the TREAT core. This research will be used to support the analysis of reactor operation and experiment planning. Traditional approaches, i.e., Monte Carlo simulations, for such tasks require a large number of model evaluations under both steady state and transient conditions. However, exhaustive sets of sampling runs are simply not practical for TREAT since evaluations of each model can take a substantial amount of time. A solution to this problem is to build a surrogate model employing a smaller number of model evaluations. The surrogate model is designed to provide an accurate and fast approximation to the detailed TREAT model. In this manner, the surrogate model can be used in place of the TREAT transient model for uncertainty and sensitivity analysis at a fraction of the computational cost. The high dimensional model representation (HDMR) method is employed here for approximating the responses of interest in terms of functions of lower dimensions. This method has been implemented in RAVEN and all TREAT transient calculations are performed using the multiphysics code MAMMOTH. In addition, a code interface has also been developed in order to couple RAVEN with MAMMOTH. The surrogate model is then tested and compared for TREAT startup testing program Transient Test no. 15 core, and its accuracy is measured by its ability to estimate the mean, coefficient of variance, skewness, excess kurtosis and probability density function of the outputs of Transient Test no. 15 model, as well as root mean square error metric.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2013

Environmental Impacts of Wind Power Development on the Population Biology of Greater Prairie-Chickens

1. We investigated the impacts of wind power development on the demography, movements, and population genetics of Greater Prairie-Chickens (Tympanuchus cupido) at three sites in northcentral and eastern Kansas for a 7-year period. Only 1 of 3 sites was developed for wind power, the 201MW Meridan Way Wind Power Facility at the Smoky Hills site in northcentral Kansas. Our project report is based on population data for prairie chickens collected during a 2-year preconstruction period (2007-2008), a 3-year postconstruction period (2009-2011) and one final year of lek surveys (2012). Where relevant, we present preconstruction data from our field studies at reference sites in the northern Flint Hills (2007-2009) and southern Flint Hills (2006-2008). 2. We addressed seven potential impacts of wind power development on prairie chickens: lek attendance, mating behavior, use of breeding habitat, fecundity rates, natal dispersal, survival rates, and population numbers. Our analyses of pre- and postconstruction impacts are based on an analysis of covariance design where we modeled population performance as a function of treatment period, distance to eventual or actual site of the nearest wind turbine, and the interaction of these factors. Our demographic and movement data from the 6-year study period at the Smoky Hills site included 23 lek sites, 251 radio-marked females monitored for 287 bird-years, and 264 nesting attempts. Our genetic data were based on genotypes of 1,760 females, males and chicks that were screened with a set of 27 microsatellite markers that were optimized in the lab. 3. In our analyses of lek attendance, the annual probability of lek persistence during the preconstruction period was ~0.9. During the postconstruction period, distance to nearest turbine did not have a significant effect on the probability of lek persistence. However, the probability of lek persistence increased from 0.69 at 0 m to 0.89 at 30 km from turbines, and most abandoned lek sites were located <5 km from turbines. Probability of lek persistence was significantly related to habitat and number of males. Leks had a higher probability of persistence in grasslands than agricultural fields, and increased from ~0.2 for leks of 5 males, to >0.9 for leks of 10 or more males. Large leks in grasslands should be a higher priority for conservation. Overall, wind power development had a weak effect on the annual probability of lek persistence. 3. We used molecular methods to investigate the mating behavior of prairie chickens. The prevailing view for lek-mating grouse is that females mate once to fertilize the clutch and that conspecific nest parasitism is rare. We found evidence that females mate multiple times to fertilize the clutch (8-18% of broods, 4-38% of chicks) and will parasitize nests of other females during egg-laying (~17% of nests). Variable rates of parentage were highest in the fragmented landscapes at the Smoky Hills field site, and were lower at the Flint Hills field site. Comparisons of the pre- and postconstruction periods showed that wind energy development did not affect the mating behaviors of prairie chickens. 4. We examined use of breeding habitats by radio-marked females and conducted separate analyses for nest site selection, and movements of females not attending nests or broods. The landscape was a mix of native prairie and agricultural habitats, and nest site selection was not random because females preferred to nest in grasslands. Nests tended to be closer to turbines during the postconstruction period and there was no evidence of behavioral avoidance of turbines by females during nest site selection. Movements of females not attending nests or broods showed that females crossed the site of the wind power development at higher rates during the preconstruction period (20%) than the postconstruction period (11%), and that movements away from turbines were more frequent during the postconstruction period. Thus, wind power development appears to affect movements in breeding habitats but not nest site selection of female prairie chickens during the breeding season. 5. We tested the effects of wind power development on five components of female fecundity: timing of clutch initiation, clutch size of first nests and renests, nest survival, and hatchability of eggs. Average date of clutch initiation was 26 April, clutch size was 12.7 and 10.6 eggs for first nests and renests, probability of nest survival was low at 0.18, but egg hatchability was high at 0.79. Wind power development had no impact on reproductive effort or nesting success, and all five components of fecundity were not related to treatment period or distance to turbine. Nest survival was the main factor limiting reproductive output of female prairie chickens and most losses were due to predation. Daily nest survival was strongly related to vegetative cover at the nest. Changes to rangeland management practices that would double nesting cover from 2.5 to 5 dm would triple the probability of nest survival from 0.17 to 0.52. Grass and forb cover had weak positive effects on daily nest survival whereas shrub cover, proximity to woodlands, and recent rainfall had negative effects. Reproductive performance of prairie chickens is low in managed rangelands in northcentral Kansas and efforts to improve range conditions and reduce predator activity would aid recovery of prairie chicken populations. 6. We used molecular methods to investigate patterns of natal dispersal in prairie chickens. High rates of nest failure limited the number of young that we could sample. Direct detections of natal dispersal were limited because survival of newly hatched chicks to become adults were low and because we were unable to detect dispersal distances outside of our study area. Direct observations of natal movements were limited and were inadequate to make conclusions about the potential impacts of wind energy development on natal dispersal. Spatial correlograms of genetic distance among males at leks were a more sensitive measure of population structure, and indicated a weak effect of wind energy development on the spatial genetic structure of prairie chickens. 7. We tested the effects of wind power development on female survival with time-to-event models, and on residual body mass of males with analyses of covariance. Distance to turbine and the interaction of distance and treatment period had no effect on female survival. Contrary to predictions of negative impacts of wind power development, the probability of female was lowest during the preconstruction period (0.274) and increased significantly during the postconstruction period (0.543). Inspection of hazard functions indicated that the difference in annual survival could be attributed to a higher risk of mortality during the lekking season in the preconstruction period. We suggest that wind power development may have improved ecological conditions for prairie chickens by disrupting the foraging behavior of diurnal raptors that kill prairie chickens at lek sites. In support of this idea, raptor kills tended to be farther from turbines during the postconstruction period whereas mammalian kills were closer. Analyses of the major causes of mortality did not support our hypothesis because the odds of raptor predation were greater after development was completed. Most mortality losses of radio-marked females were due to predation, and losses to collision mortality or harvest were rare events. Low rates of natural mortality during fall and winter imply that harvest is likely to be additive mortality in prairie chickens. Wind energy development reduced the residual body mass of male Greater Prairie-Chickens at lek sites near turbines. Low values of residual body mass could have a negative impact on individual survival or fecundity rates, or may be related to predation risk and flight performance of males displaying at open lek sites. 8. We tested for impacts of wind power development on population numbers of prairie chickens with monitoring of male numbers at leks, and with genetic measures of population structure. Lek counts indicated that wind power development did not affect the population size of prairie chickens. Peak counts of males at leks were recorded the first year after construction was completed and the highest rates of population change were observed during the interval when the wind power facility was constructed. Population numbers of prairie chickens near and distant from turbines appeared to covary in parallel, probably because bird numbers were linked by dispersal movements. Estimates of population viability based on genetic diversity, effective population size and rates of population exchange did not show annual changes and were unaffected by wind development during our study. Estimates of relatedness among males at the same and different leks suggested that wind power development has either reduced dispersal rates or changed settlement patterns, leading to higher rates of relatedness among males displaying at the same lek site. 9. Greater Prairie-Chickens were not strongly affected by wind power development in Kansas. Negative impacts of wind power development included a trend for reductions in lek persistence near turbines, behavioral avoidance of turbines by females during their breeding season movements, and changes in the genetic structure of males at leks that were consistent with reduced dispersal or recruitment rates. We found no impacts of wind power development on nest site selection, female reproductive effort or nesting success, or population numbers. Positive impacts of wind power development included an increase in female survival rates. We hypothesized that the unexpected increase in female survival was related to changes in trophic interactions and disruption of the foraging behavior of raptors that kill prairie chickens at lek site. 10. Research funding for this project included a grant from the 20% Wind by 2030 Program of the Department of Energy (this final report), and grants from the Kansas Department of Wildlife, Parks and Tourism, the National Fish and Wildlife Federation, and initial funding from the National Wind Coordinating Collaborative. Research products from data collected during the preconstruction period have included five peer-reviewed research articles and two PhD dissertations at Kansas State University. Additional manuscripts are in review for possible publication in 2013-2014.

17 WIND ENERGY↗

DOE OSTI2017

Surrogate Models for TREAT Transient Calculations

The Transient Reactor Test Facility (TREAT) at Idaho National Laboratory (INL) is currently preparing for resumption of operation after more than 20 years in a long-term standby state. In parallel with preparation for operations and a new transient testing program, the DOE-NE Advanced Modeling and Simulation program (NEAMS) has been supporting the development of three-dimensional transient multi-physics modeling capabilities within the reactor physics application MAMMOTH [1]. Recently, the MAMMOTH reactor physics team has begun development of uncertainty and sensitivity analysis for the TREAT core. This research will be used to support the analysis of reactor operation and experiment planning. Traditional approaches, i.e., Monte Carlo simulations, for such tasks require a large number of model evaluations under both steady state and transient conditions. However, exhaustive sets of sampling runs are simply not practical here since evaluations of each model can take a substantial amount of time. A solution to this problem is to build a surrogate model employing a smaller number of model evaluations. The surrogate model is designed to provide an accurate and fast approximation to the detailed TREAT model. In this manner, the surrogate model can be used in place of the TREAT transient model for uncertainty and sensitivity analysis at a fraction of the computational cost. The high dimensional model representation (HDMR) [2] method is employed here for approximating the responses of interest in terms of functions of lower dimensions. This method has been implemented in RAVEN [3] and all TREAT transient calculations are performed using the multiphysics code MAMMOTH. In addition, a code interface has also been developed in order to couple RAVEN with MAMMOTH. The surrogate model is then tested and compared for the MAMMOTH model of the TREAT startup testing program Transient Test #15 core [4], and its accuracy is measured by its ability to estimate the mean, coefficient of variance, skewness, excess kurtosis and probability density function of the outputs of Transient Test #15 model, as well as root mean square error metric.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Wellbore Fracture Imaging Using Inflow Detection Measurements

One of the most striking measurements taken during DOE’s EGS Collab project at the 4850-foot depth location was the so-called ‘sewer cam’, which enabled direct visualization of the flow of water into the production well through fractures during the stimulation. The ability to see directly which fractures were flowing and (roughly) how much was a breakthrough in understanding the topology of the created fracture network. Achieving this kind of fracture flow imaging at FORGE would be more challenging because of the 225°C temperature, but equally or even more valuable if it could be achieved. In 2017, a joint project between Sandia and Stanford developed a downhole tool concept to measure the enthalpy of multiphase fluid entering a geothermal well from individual fractures (Gao et al., 2017). For the FORGE project, measuring enthalpy is of less interest because the fluid is expected to be single-phase liquid water. However, the foundation of the device was the measurement of chloride ion concentration, which could form the basis for a direct measurement of inflow from fractures. During the 2017 project, this novel chloride sensing system was implemented into a laboratory test instrument, and we confirmed the capability of the system to measure the ion concentration of fluid entering a model wellbore through a small entry port. The wellbore was a 6-inch diameter model well, and the port was approximately 0.08 inch (2mm) in diameter. The device could measure the chloride concentration accurately even when the well was flowing in a bubbly flow. Given its accuracy, the tool should be able to identify locations of water entering the wellbore even if the ion concentration differs only slightly from that of the water in the well. It is likely that different fractures may flow slightly different chloride concentrations, which would make it feasible to detect individual fractures as well as to estimate the volume of their flow. Ultimately, we could also recognize different fractures flowing back significantly different ion concentrations after fracturing in the FORGE wells. This could be realized by adding different ions in the fracturing fluids in different fractures created at different stages of stimulation (and modifying the tool to include different ion specificity). Sandia’s tool was shown during the study to have the capability to withstand the 225°C temperature, and the electrochemical sensing elements were tested in the laboratory to 225°C at 1500 psia for 24 hours. An early implementation of the fully integrated downhole electrochemical tool, including high-temperature electronics, robust housing, and wireline truck interface, had previously been constructed and tested successfully at Sandia; thus, hardware development tasks focused on advancing the technology readiness level (TRL) of this promising technology for FORGE deployment, rather than on developing a new scientific basis for its operation. The data collection electronics in this tool allowed for several other sensors (pressure, temperature, flow spinner) to be implemented in parallel as well. The research was a new collaboration between Stanford and Sandia to modify and refine the tool for FORGE deployment, to make the downhole measurements, and to characterize the evolving fractures.

15 GEOTHERMAL ENERGY↗

DOE OSTI2020

10 MWE CDCL Large Pilot Plant – pre-FEED Study

Coal Direct Chemical Looping (CDCL) is an advanced oxy-combustion technology that has the potential to substantially reduce the energy penalty and the cost of electricity (COE) for coal-fired power generation with CO₂ capture. The Babcock & Wilcox Company (B&W) and The Ohio State University (OSU) have been collaborating on the development of an iron oxide oxygen-carrier based chemical looping technology for clean power generation with inherent carbon capture. In this process, coal is dried and pulverized prior to being transported into a moving-bed reducer. In the reducer, coal reacts with the oxygen-carrier particles, forming combustion byproducts, predominantly CO₂ and H₂O, while reducing the iron oxide oxidation state from Fe₂O₃ to a mixture of FeO and Fe. The reduced state particles are then transported to a combustor reactor and re-oxidized with air. Following the oxidation, the oxygen-carrier particles are regenerated, and a large amount of heat is released for steam production. The produced steam is sent to a turbine for electricity generation. Meanwhile, the CO₂–rich stream leaving the reducer is cooled, cleaned, and compressed for subsequent pipeline transportation and sequestration. By combining air separation and fuel conversion into a single system, the CDCL technology enables the intensification of oxy-combustion processes by eliminating the energy and cost intensive cryogenic air separation unit and thereby results in higher overall plant efficiencies and lower COE’s. The use of a moving-bed reducer results in high conversions of volatile hydrocarbons and high CO₂ purity, which reduces the cost of downstream CO₂ purification for sequestration or utilization. The Babcock & Wilcox Company in collaboration with The Ohio State University, Johnson & Matthey, The Electric Power Research Institute, and Dover Light & Power completed a Preliminary Front-End Engineering and Design (Pre-FEED) study of a 10 MWe coal-direct chemical looping (CDCL) pilot plant. The planned system is a modular 10 MWe CDCL large pilot facility consisting in 4 modules of 2.5 MWe each, working in parallel and to be hosted within the current structure at the City of Dover’s Municipal Power Plant. The CDCL system can achieve auto-thermal operation and includes a sub-critical steam cycle for power generation. The pilot system was designed to demonstrate full commercial operation at a reduced scale. The coal distribution per plan area is a representative slice of larger commercial arrangements. The system includes a CO₂ recycle system, but it does not include a compression system. The large pilot includes all the environmental control equipment and oxygen carrier and ash handling systems. As part of the project, the Team performed laboratory testing and carried out multiple pilot test campaigns to obtain design and performance information at the 250 kWth CDCL pilot facility at the Babcock & Wilcox Company’s Research Center. The Team demonstrated sustained operations at designed coal inputs, high coal conversion, high CO₂ purity, heat generation on the combustor, low carbon carryover between reactors and low particle attrition. Emissions generated in the reducer reactor were identified as SO₂ and NO x . The commercial manufacturing cost of oxygen carrier particle was evaluated by JM. A particle manufacturing report was generated and submitted to the DOE. Based on the results from the pilot tests and the pre-FEED design efforts, a techno-economic analysis was performed. The study shows that the CDCL process is a promising carbon-friendly technology capable of producing electricity with high efficiency. The estimated COE of the supercritical CDCL plant is $83.3 / MWh, which meets DOE’s target of less than 30% increase in COE when compared to a supercritical PC plant without CO₂ capture. This is the lowest among the existing carbon capture technologies (post-combustion and oxy-combustion) for fossil fuel power plant. CDCL is evaluated to be the most promising technology for carbon capture from the economic aspect.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2018

Final Report: Optimal Model Complexity in Geological Carbon Sequestration: A Response Surface Uncertainty Analysis

The critical component of a risk assessment study in evaluating GCS is an analysis of uncertainty in CO2 modeling. In such analyses, direct numerical simulation of CO2 flow and leakage requires many time-consuming model runs. Alternatively, analytical methods have been developed which allow fast and efficient estimation of CO2 storage and leakage, although restrictive assumptions on formation rock and fluid properties are employed. In this study, an intermediate approach is proposed based on the Design of Experiment and Response Surface methodology, which consists of using a limited number of numerical simulations to estimate a prediction outcome as a combination of the most influential uncertain site properties. The methodology can be implemented within a Monte Carlo framework to efficiently assess parameter and prediction uncertainty while honoring the accuracy of numerical simulations. The choice of the uncertain properties is flexible and can include geologic parameters that influence reservoir heterogeneity, engineering parameters that influence gas trapping and migration, and reactive parameters that influence the extent of fluid/rock reactions. The method was tested and verified on modeling long-term CO2 flow, non-isothermal heat transport, and CO2 dissolution storage by coupling two-phase flow with explicit miscibility calculation using an accurate equation of state that gives rise to convective mixing of formation brine variably saturated with CO2. All simulations were performed using three-dimensional high-resolution models including a target deep saline aquifer, overlying caprock, and a shallow aquifer. To evaluate the uncertainty in representing reservoir permeability, sediment hierarchy of a heterogeneous digital stratigraphy was mapped to create multiple irregularly shape stratigraphic models of decreasing geologic resolutions: heterogeneous (reference), lithofacies, depositional environment, and a (homogeneous) geologic formation. To ensure model equivalency, all the stratigraphic models were successfully upscaled from the reference heterogeneous model for bulk flow and transport predictions (Zhang & Zhang, 2015). GCS simulation was then simulated with all models, yielding insights into the level of parameterization complexity that is needed for the accurate simulation of reservoir pore pressure, CO2 storage, leakage, footprint, and dissolution over both short (i.e., injection) and longer (monitoring) time scales. Important uncertainty parameters that impact these key performance metrics were identified for the stratigraphic models as well as for the heterogeneous model, leading to the development of reduced/simplified models at lower characterization cost that can be used for the reservoir uncertainty analysis. All the CO2 modeling was conducted using PFLOTRAN – a massively parallel, multiphase, multi-component, and reactive transport simulator developed by a multi-laboratory DOE/SciDAC (Scientific Discovery through Advanced Computing) project (Zhang et al., 2017, in review). Within the uncertainty analysis framework, increasing reservoir depth were investigated to explore its effect on the uncertainty outcomes and the potential for developing gravity-stable injection with increased storage security (Dai et al., 20126; Dai et al., 2017, in review). Finally, to accurately model CO2 fluid-rock reactions and resulting long-term storage as secondary carbonate minerals, a modified kinetic rate law for general mineral dissolution and precipitation was proposed and verified that is invariant to a scale transformation of the mineral formula weight. This new formulation will lead to more accurate assessment of mineral storage over geologic time scales (Lichtner, 2016).

58 GEOSCIENCES↗

DOE OSTI2025

Parallel-in-time quantum simulation via Page and Wootters quantum time

In the past few decades, researchers have created a veritable zoo of quantum algorithms by drawing inspiration from classical computing, information theory, and even from physical phenomena. Here, we present quantum algorithms for parallel-in-time simulations that are inspired by the Page and Wootters formalism. In this framework, and thus in our algorithms, the classical time variable of quantum mechanics is promoted to the quantum realm by introducing a Hilbert space of “clock” qubits that are then entangled with the “system” qubits. We show that our algorithms can compute temporal properties over 𝑁 different times of many-body systems by only using log⁡(𝑁) clock qubits. As such, we achieve an exponential trade-off between time and spatial complexities. In addition, we rigorously prove that the entanglement created between the system qubits and the clock qubits has operational meaning, as it encodes valuable information about the system’s dynamics. We also provide a circuit depth estimation of all the protocols, showing a running time advantage in computation times over traditional sequential-in-time algorithms. In particular, for the case when the dynamics are determined by the Aubry-Andre model, we present a hybrid method for which our algorithms have a depth that only scales as 𝒪⁡(log⁡(𝑁)⁢𝑛). As a by-product, we can relate the previous schemes to the problem of equilibration of an isolated quantum system, thus indicating that our framework enables a new dimension for studying dynamical properties of many-body systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

The value of integrating a geothermal district heating system into a microgrid

As electrical grids increasingly rely on variable renewable energy, maintaining reliability and cost efficiency becomes more complex. To address these challenges, this study analyzed the integration of geothermal district heating as a grid-responsive thermal resource within a microgrid in Tuttle, Oklahoma. Building energy modeling using EnergyPlus estimated annual district heating demand at 2.9 GWh, with a peak load of 2.8 MW th . Techno-economic analyses were conducted to meet the heating demand under three geothermal scenarios, varying by production depth, flow rate, and thermal output, each supplemented by natural gas peaking boilers. In parallel, equivalent electrical load profiles were developed using typical coefficients of performance (COPs) for air-source heat pumps and electric boilers to establish an electrified baseline scenario. A complete end-use electrical load profile was also developed for the microgrid using Cambium dataset. The modeling results demonstrated reliable and economic operation of the geothermal systems over 30 years, with COPs ranging from 2.6 to 8.9 and the lowest levelized heating cost at $\$$54.6/MWh. Geothermal integration reduced electricity consumption by up to 94.7 % compared to the non-geothermal base case, yielding annual energy savings of up to $\$$803 k. Avoided grid costs ranged from $\$$65 k–$\$$147 k per year, with individual events avoiding up to $\$$4,863 per hour. Grid-responsive operation further reduced wholesale energy costs by 53–56 %. These findings demonstrate geothermal heating, traditionally treated as a non-grid-responsive thermal resource, can be reconfigured to support dynamic grid services, offering a scalable pathway to enhance reliability and reduce costs in renewable-rich microgrids and district heating networks.

15 GEOTHERMAL ENERGY↗

DOE OSTI2013

COMMERCIALIZATION OF AN ATMOSPHERIC IRON-BASED CDCL PROCESS FOR POWER PRODUCTION. PHASE I: TECHNOECONOMIC ANALYSIS

Coal Direct Chemical Looping (CDCL) is an advanced oxy-combustion technology that has potential to enable substantial reductions in the cost and energy penalty associated with carbon dioxide (CO2) capture from coal-fired power plants. Through collaborative efforts, the Babcock & Wilcox Power Generation Group (B&W) and The Ohio State University (OSU) developed a conceptual design for a 550 MWe (net) supercritical CDCL power plant with greater than 90% CO2 capture and compression. Process simulations were completed to enable an initial assessment of its technical performance. A cost estimate was developed following DOE’s guidelines as outlined in NETL’s report “Quality Guidelines for Energy System Studies: Cost Estimation Methodology for NETL Assessments of Power Plant Performance”, (2011/1455). The cost of electricity for the CDCL plant without CO2 Transportation and Storage cost resulted in $ $102.67 per MWh, which corresponds to a 26.8 % increase in cost of electricity (COE) when compared to an air-fired pulverized-coal supercritical power plant. The cost of electricity is strongly depending on the total plant cost and cost of the oxygen carrier particles. The CDCL process could capture further potential savings by increasing the performance of the particles and reducing the plant size. During the techno-economic analysis, the team identified technology and engineering gaps that need to be closed to bring the technology to commercialization. The technology gaps were focused in five critical areas: (i) moving bed reducer reactor, (ii) fluidized bed combustor, (iii) particle riser, (iv) oxygen-carrier particle properties, and (v) process operation. The key technology gaps are related to particle performance, particle manufacturing cost, and the operation of the reducer reactor. These technology gaps are to be addressed during Phase II of project. The project team is proposing additional lab testing to be completed on the particle and a 3MWth pilot facility be built to evaluate the reducer reactor performance among other aspects of the technology. A Phase II proposal was prepared and submitted to DOE. The project team proposed a three year program in Phase II. Year 1 includes lab testing and particle development work aimed at improving the chemical and mechanical properties of the oxygen carrier particle. In parallel, B&W will design the 3MWt pilot plant. Any improvements to the particle performance discovered in year 1 that would impact the design of the pilot will be incorporated into the final design. Year 2 will focus on procurement of materials and equipment, and construction of the pilot plant. Year 3 will include, commissioning, start-up, and testing in the pilot. Phase I work was successfully completed and a design and operating philosophy for a 550 MWe commercial scale coal-direct chemical looping power plant was developed. Based on the results of the techno-economic evaluation, B&W projects that the CDCL process can achieve 96.5% CO2 capture with a

Coal Combustion↗

DOE OSTI2025

Introduction: Neuromorphic Materials

The explosive growth in data collection and the need to process it efficiently, as well as the desire to automate increasingly complex tasks in transportation, medical care, manufacturing, security and many other fields have motivated a growing interest in neuromorphic computing. Unlike the binary, transistorbased ON/OFF logic gates and separate logic and memory functionalities employed in digital computing, neuromorphic computing is inspired by animal brains that use interconnected synapses and neurons to perform processing, storage and transmission of information at the same location, while only consuming ~20 W or less of power. Motivated by the brain’s efficiency, adaptability, self-learning and resiliency qualities, neuromorphic computing can be broadly defined as an approach to processing and storing information using hardware and algorithms inspired by models of biological neural systems. Present research in neuromorphic computing encompasses approaches that vary significantly in their degree of neuro-inspiration, from systems that only incorporate features such as asynchronous, event-driven operation or use crossbar arrays of non-volatile memory (NVM) elements to accelerate deep neural networks (DNNs), to designs that embrace the extreme parallelism, sparsity, reconfigurability, adaptability, complexity and stochasticity observed in nervous systems. The term ‘neuromorphic’ computing is often credited to Carver Mead, who in the 1980s investigated Si-based analog electronics to replicate functions of the animal retina. Earlier important advances in this field include the work of Frank Rosenblatt, who proposed the concept of the perceptron, Bernard Widrow, who used this concept to build one of the first analog neural networks, the Adaline and many other researchers (see ref. 6 for an historical perspective on neuromorphic computing). With the recent increase in the use of artificial intelligence and large language models, and rising concerns over the associated energy costs, interest in neuromorphic hardware has expanded rapidly. According to some estimates, driven largely by the drastic growth in the training use of artificial intelligence (AI) models using the current computing architectures, the energy cost of computing is projected to reach the energy supply worldwide by 2045. Furthermore, while this is not a realistic outcome, it means that, if more efficient computing technologies are not developed -- soon -- the world will soon become one where demand for energy and market constraints limit the continued increase of societal access to AI and cloud services from data centers. Data centers used for training and use of these models consume hundreds of terawatt hours of electricity, already past 4% of the US electricity demand.

Circuits↗

DOE OSTI2025

Status Report on Design of In-situ Thermomechanical Testing at LANSCE

Nuclear fuel encounters severe thermomechanical environments in which its mechanical response is determined by its microstructure, temperature and stress level histories. Simulating the response of such microstructures is crucial for predicting both performance and transient fuel mechanical responses and experimental verification of such predictions is therefore of great interest. While most of the deformation in a nuclear fuel rod occurs in the cladding, deformation of the fuel itself is still of interest with deformation mechanisms at operating temperature and above including creep, swelling, cracking as well as pellet-clad interaction. Characterization of these properties and understanding of the underlying deformation phenomena at operating or excursion temperatures is therefore of great importance for development and ultimately licensing of improved and novel nuclear fuel forms. Diffraction techniques offer unique insight on the atomistic (e.g. crystal structure) and microstructure (e.g. phase transformations, texture, defects) length scales and have a long history of providing unique data to inform relevant deformation models that enable the required predictive capabilities. For example, dislocations lead to diffraction peak broadening that can be characterized to estimate the dislocation density and study the role of dislocations on the deformation while measuring lattice strains allows to studie load sharing in two phase materials. In this report the requirements for a sample environment for high temperature deformation of nuclear fuels are defined. The HIPPO neutron time-of-flight diffractometer at LANSCE will host this sample environment and is also described. This instrument covers diffraction angles from 140° to 40° and is also equipped with an event-mode neutron imaging detector system, enabling energy-resolved neutron imaging in parallel with the diffraction that could measure sample temperature from Doppler broadening of neutron absorption resonances or measure pore densities from changes in the attenuation. Designs of devices to characterize thermomechanical properties of nuclear fuel without diffraction are also considered to guide the design. While this report is focused on applications for nuclear fuels, the device can also characterize cladding, moderator or structural materials and therefore contribute to other fields of research and development for advanced reactors. The temperatures planned to be reached are above 2000℃, thus enabling characterization of LWR reactor fuels under accident scenarios but also reaching temperatures of fuels developed for nuclear thermal propulsion and providing opportunities to characterize those. In conjunction with the energy-resolved neutron imaging detector, this setup would allow to measure neutron cross-sections at high temperatures, filling a gap towards development of reactors operating at high temperatures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Hybrid PDES Simulation of HPC Networks Using Zombie Packets

Although high-fidelity network simulations have proven to be reliable and cost-effective tools to peer into architectural questions for high-performance computing (HPC) networks, they incur a high resource cost. The time spent in simulating a single millisecond of network traffic in the highest detail can take hours, even for static, well-behaved traffic patterns such as uniform random. Surrogate models offer a significant reduction in runtime, yet they cannot serve as complete replacements and should only be used when appropriate. Thus, there is a need for hybrid modeling, where high-fidelity simulation and surrogates run side-by-side. Here, we present a surrogate model for HPC networks in which: packets bypass the network, while the network state is left untouched, i.e., suspended. To bypass the network, we use historical data to estimate the arrival time at which every packet should be scheduled at; to suspend the network, all in-flight packets are scheduled to arrive at their destinations, and are kept in the system to awaken as zombies when switching back to high-fidelity. Speedup for a hybrid model is relative to the proportion of surrogate to high-fidelity. This light-weight surrogate obtained up to 76× speedup. Keeping the zombies in the network showed an increase in the accuracy of the high-fidelity simulation on restart when compared to restarting the network from an empty state.

HPC networks↗

DOE OSTI2017

FINAL REPORT - Development of High Pressure Hydrogen Storage Tank for Storage and Gaseous Truck Delivery

The “Development of High Pressure Hydrogen Storage Tanks for Storage and Gaseous Truck Delivery” project [DE-FG36-08GO18062] was initiated on 01 July 2008. Hexagon Lincoln (then Lincoln Composites) received grant funding from the U.S. Department of Energy to support the design and development of an improved bulk hauling and storage solution for hydrogen in terms of cost, safety, weight and volumetric efficiency. The development of this capability required parallel development and qualification of large all-composites pressure vessels, a custom ISO container to transport and store said tanks, and performance of trade studies to identify optimal operating pressure for the system. Qualification of the 250 bar TITAN® module was completed in 2009 with supervision from the American Bureau of Shipping [ABS], and the equipment has been used internationally for bulk transportation of fuel gases since 2010. Phase 1 of the project was successfully completed in 2012 with the issuance of USDOT SP 14951, the special permit authorizing the manufacture, marking, sale and use of TITAN® Mobile Pipeline® equipment in the United States. The introduction of tube trailers with light weight composite tankage has meant that 2 to 3 times as much gaseous fuel can be transported with each trip. This increased hauling efficiency offers dramatically reduced operating costs and has enabled a profitable business model for over-the-road compressed natural gas delivery. The economic drivers of this business opportunity vary from country to country and region to region, but in many places gas distribution companies have realized profitable operations. Additional testing was performed in 2015 to characterize hydrogen-specific operating protocols for use of TITAN® systems in CHG service at 250 bar. This program demonstrated that existing compression and decompression methodologies can efficiently and safely fill and unload lightweight bulk hauling systems. Hexagon Lincoln and U.S. DOE agreed to continue into Phase 2 of the project without pursuing the development of higher pressure capabilities as originally planned. At 250 bar, development of equipment for hydrogen transport is supported by strong activity in the adjacent natural gas transportation sector. Trade studies performed since 2011 indicate optimization of hauling efficiency and system cost for hydrogen transport at about 350 bar (5076 psi). However, due to reduced efficiency of compression of natural gas above 250 bar, 350 bar operation is not an attractive option for natural gas transportation. The CHG market is not developed at this time, and it is difficult to forecast the arrival of significant revenues. On the investment side, the cost to fully qualify a large tank module at 350 bar is estimated at $3MM to $5MM. There is insufficient CHG market definition to support a stand-alone business case for this investment without near term revenue in the adjacent CNG transportation market. Therefore development of a 350 bar TITAN® system was deferred and not pursued under this project. Hexagon Lincoln continues to support the development of tankage and equipment for operation at 350 bar and above; with 700 bar vehicle tanks and 950 bar tanks for ground storage applications. Phase 2 activities were focused on reducing system cost, increasing system capacity, increasing system safety and characterization of polymer material performance specific to hydrogen pressure vessel usage. With the successful launch of TITAN® modules and trailers in natural gas transportation, over 600 units have been produced through the end of 2016, resulting in improved purchasing power for raw materials and manufactured components. This has allowed Hexagon Lincoln to approach the current project goals for system cost. At $590/kg of compressed hydrogen delivered, the system cost of the baseline TITAN® module is below the project’s 2015 target of $730/kg H2 delivered, and very close to the project’s 2020 target of $575/kg H2 delivered. [Based on product pricing in 1Q2017.] Emphasis was placed on configuration of larger capacity systems within the vehicle weights and dimensions allowed on federal and state highways in the United States and other countries. These activities resulted in the design and development of integrated tube trailer systems that have increased delivery capacities by 45%. The hydrogen delivery capacity of our largest system is 845 kg, exceeding the project’s 2015 target of 700 kg H2 delivered. Emerging technologies offering improvement of the safety systems used on the equipment were investigated, with particular focus on improving the reliability and cost of the emergency venting system for fire protection. Finally, investment in our materials laboratory improved detection and characterization of hydrogen-induced damage in polymer materials, supporting the development of operational protocols to avoid damage to pressure vessel liners and valve components.

08 HYDROGEN↗

DOE OSTI2026

A biosynthetic gene cluster for three post-chorismate pathways in Arabidopsis

Chorismate is a branch-point metabolite in the biosynthesis of aromatic amino acids, vitamins, antibiotics and various other aromatic products in bacteria, fungi and plants. Although 13 chorismate-utilizing enzymes have been identified in bacteria, only 6 have been described in plants, where an estimated 30% of all photosynthetically fixed carbon passes through chorismate. Here, in this study, we describe a biosynthetic gene cluster (BGC) consisting of five core genes, including two reductases, two methyltransferases and one glucosyltransferase. Genetic and biochemical evidence shows that these five enzymes collectively give rise to three biosynthetic pathways, each originating from chorismate: two parallel pathways produce a class of non-aromatic, isomeric compounds abundant in the roots of Arabidopsis thaliana, whereas the third pathway produces methylated and glucosylated chorismate derivatives that subsequently react non-enzymatically with glutathione. Genome analysis revealed that variants of this BGC are present in some but not all species in the Brassicaceae family. Taken together, our study uncovered a BGC, containing three chorismate-utilizing enzymes, that controls three distinct post-chorismate pathways in A. thaliana. This work not only advances our understanding of carbon flow in this model plant but also highlights that the biochemical complexity encoded by plant BGCs is greater than previously appreciated.

Peng, Meng [Ghent Univ. (Belgium); Flemish Institu↗

DOE OSTI2013

Analysis & Tools to Spur Increased Deployment of “Waste Heat” Rejection/Recycling Hybrid Ground-source Heat Pump Systems in Hot, Arid or Semiarid Climates Like Texas

This project team analyzed supplemental heat rejection/recovery (SHR) devices or systems that could be used in hybrid ground source heat pump (HGHP) systems located in arid or semi-arid regions in southwestern U.S. Identification of effective SHR solutions would enhance the deployment of ground source heat pumps (GHP) in these regions. In a parallel effort, the team developed integrated GHP models that coupled the building load, heat pump, and ground loop subsystems and which could be applied to residential and commercial office buildings. Then GHP and HGHP performances could be compared in terms of operational performance and life-cycle costs. Several potential SHR devices were analyzed by applying two strategies: 1) to remove heat directly from the water in the ground loop before it enters the ground and 2) to remove heat in the refrigerant loop of the vapor compression cycle (VCC) of the heat pump so less heat is transferred to the water loop at the condenser of the VCC. Cooling towers, adsorption coolers, and thermoelectric liquid coolers were included in strategy 1, and expanded desuperheaters, thermosyphons, and an optimized VCC were included in strategy 2. Of all SHR devices analyzed, only the cooling tower provided a cost-effective performance enhancement. For the integrated GHP model, the project team selected the building load model HAMBASE and its powerful computational Simulink/MatLab platform, empirical performance map models of the heat pumps based upon manufacturers’ performance data, and a ground loop model developed by Oklahoma State University and rewritten for this project in Simulink/MatLab. The design process used GLHEPRO, also from Oklahoma State University, to size the borehole fields. The building load and ground loop models were compared with simulations from eQuest, ASHRAE 140-2008 standards, EnergyPlus, and GLHEPRO and were found to predict those subsystems’ performance well. The integrated GHP model was applied to a 195m 2 (2100ft 2 ) residential building and a 4,982m 2 (53,628ft 2 ) three-story commercial office building, and it ran 10-15 year simulations. The integrated GHP model and its Simulink platform provided residential data, ranging from seconds to years, and commercial office building data, ranging from minutes to years. A cooling tower model was coupled to the base case integrated GHP model for the residential building and the resulting HGHP system provided a cost-effective solution for the Austin, TX location. Simulations for both the residential and commercial building models were run with varying degrees of SHR (device/system not identified) and the results were found to significantly decrease installation costs, increase heat pump efficiency (lower entering water temperature), and prolong the lifetime of the borehole field. Lifetime cycle costs were estimated from the simulation results. Sensitivity studies on system operating performance and lifetime costs were performed on design parameters, such as construction materials, borehole length, borehole configuration and spacing, grout conductivity, and effects of SHR. While some of the results are intuitive, these studies provided quantitative estimates of improved performance and cost. One of the most important results of this sensitivity study is that overall system performance is very sensitive to these design parameters and that modeling and simulation are essential tools to design cost-effective systems.

15 GEOTHERMAL ENERGY↗