Search NASA⌕ Search

SEARCH · Search NASA

Results for “forward rate”

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

Long-lived axionlike particles from tau decays

Axionlike particles (ALPs) are well-motivated examples of light, weakly coupled particles in theories beyond the Standard Model. In this work, we study long-lived ALPs coupled exclusively to leptons in the mass range between 2⁢𝑚 𝑒 and 𝑚 𝜏 − 𝑚 𝑒 . For anarchic flavor structure the leptophilic ALP production in tau decays or from ALP-tau bremsstrahlung is enhanced thanks to derivative couplings of the ALP and can surpass production from electron and muon channels, especially for ALPs heavier than 𝑚 𝜇 . Using past data from high-energy fixed-target experiments such as CHARM and BEBC we place new constraints on the ALP decay constant 𝑓 𝑎 , reaching scales as high as 𝒪⁡(10 8 ) GeV in lepton-flavor-violating channels and 𝑓 𝑎 ∼ 𝒪⁡(10 2 ) GeV in lepton-flavor-conserving ones. We also present projections for the event-rate sensitivity of current and future detectors to ALPs produced at the Fermilab Main Injector, the CERN SPS, and in the forward direction of the LHC. We show that SHiP will be sensitive to 𝑓 𝑎 values that are over an order of magnitude above the existing constraints.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Novel features of attractors and transseries in nonconformal Bjorken flows

In this work we investigate the impact of conformal symmetry breaking on hydrodynamization of a far-from-equilibrium fluid. We find a new kind of transseries solutions for the nonconformal hydrodynamic equations of a longitudinal boost invariant expanding plasma. The new transseries solutions unveil a rich physical structure which arises due to the interplay of different physical scales. In the perfect fluid case the nonconformal speed of sound slows down the cooling of the temperature due to the emergence of logarithmic corrections that depends on the mass of the particle. Further, these terms propagate into the perturbative and nonperturbative sectors of the transseries once viscous corrections are included. The logarithmic mass contributions increase the asymptotic value of the Knudsen number while decreasing the damping rate of the transient nonhydrodynamic modes and thus, yielding to an extremely slow hydrodynamization process where flow lines merge to their forward attractor at extremely late times. The early time free streaming expansion is modified and receives logarithmic mass corrections induced by the shear-bulk couplings. The global flow structure and numerical analyses carried out in our work demonstrate the existence of local attraction at early and late times for the shear viscous tensor and bulk viscous pressure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Parametric Study of SANSMIC Input and Resulting Impact on Predicted Cavern Geometry and Leaching Efficiency

The Sandia Solution Mining Code (SANSMIC) has been used for many years to examine the development of salt cavern geometry, both in a confirmatory manner with comparisons made to real-world sonar data and in a predictive manner when updated sonar data are not available. SANSMIC models require some modeling choices in order to incorporate real-world data. Key modeling choices include the vertical resolution of cavern geometry to implement, as well as how to incorporate daily raw water injection data into the SANSMIC model. This report documents five studies that address the impact of the modeling choices on the predicted cavern geometries and calculated leaching efficiencies. In most cases, hypothetical cylindrical initial cavern geometries are used to provide a common baseline against which to test the systematic variation of input variables including cavern radius, oil-brine-interface (OBI) depth, vertical cell size, raw water injection rate, raw water injection duration, workover time, and number of leaching stages. The use of smaller cell sizes is recommended moving forward to provide a better one-to-one relationship between sonar data and the modeled cavern. A new methodology for incorporating raw water injection data is also recommended, in order to more closely model real-world injection and workover times. Overall, the systematic studies performed here have increased our confidence in previous SANSMIC model results, as well future use of the code for predicting leaching effects on cavern geometries. Some minor changes to modeling choices are recommended, which can easily be applied with the version of SANSMIC currently under development.

42 ENGINEERING↗
DOE OSTI2021

Development and validation of a gamma-ray detection response model for neutron-diagnosed subcritical experiments

In a neutron-diagnosed subcritical experiment, a highly sensitive and fast detector is required to accurately determine the true fission γ-ray emission rate from an object composed of special nuclear material under interrogation, in addition to a proper detector response model for use in forward modeling. We present the development of a detector response model that is based in GEANT4 and uses post-processing techniques to include photomultiplier tube time and pulse height response effects. We also define a set of detector performance evaluation metrics that emphasized our experimental objectives. A suite of time and pulse height response measurements for several diverse detector assemblies are described. These assemblies were simulated with this response model and benchmarked against the laboratory measurements. This response model reproduces the pulse height curve shapes and performance hierarchy of the measured detectors well. Additionally, the simulated time response curves exhibit reasonable agreement with counterpart measurements. These measurements and simulation results together suggest that a large volume EJ-299-49 plastic scintillator that is mated to an ultraviolet transmitting plastic light guide and relies on total internal reflection for light transport to the photocathode will provide significantly faster timing, increase pulse height by a factor of 3, and require 25% less γ rays to meet experimental needs, compared to the candidate detector design developed prior to the availability of this response model. Application of this response model and associated performance metrics in the context of large-scale design sensitivity studies is also discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Catalytic Resonance Theory: The Catalytic Mechanics of Programmable Ratchets

Catalytic reaction networks of multiple elementary steps operating under dynamic conditions via a programmed input oscillation are difficult to interpret and optimize due to reaction system complexity. To understand these dynamic systems, individual elementary catalytic reactions oscillating between catalyst states were evaluated to identify their three fundamental characteristics that define their ability to promote reactions away from equilibrium. First, elementary catalytic reactions exhibit directionality to promote reactions forward or backward from equilibrium as determined by a ratchet directionality metric comprised of the input oscillation duty cycle and the reaction rate constants. Second, catalytic ratchets are defined by the catalyst state of strong or weak binding that permit reactants to proceed through the transition state. Third, elementary catalytic ratchets exhibit a cutoff frequency which defines the transition in applied frequency for which the catalytic ratchet functions to promote chemistry away from equilibrium. All three ratchet characteristics are calculated from chemical reaction parameters including rate constants derived from linear scaling parameters, reaction conditions, and catalyst state. The characteristics of the reaction network’s constituent elementary catalytic reactions provided an interpretation of complex reaction networks and a method of predicting the behavior of dynamic surface chemistry on oscillating catalysts.

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

Search for resonance-enhanced C P and angular asymmetries in the Λ c + → p μ + μ − decay at LHCb

The first measurement of the C P asymmetry of the decay rate ( A C P ) and the C P average ( Σ A FB ) and C P asymmetry ( Δ A FB ) of the forward-backward asymmetry in the muon system of Λ c + → p μ + μ − decays is reported. The measurement is performed using a data sample of proton-proton collisions, recorded by the LHCb experiment from 2016 to 2018 at a center-of-mass energy of 13 TeV, which corresponds to an integrated luminosity of 5.4 fb − 1 . The asymmetries are measured in two regions of dimuon mass near the ϕ -meson mass peak. The dimuon-mass integrated results are A C P = ( − 1.1 ± 4.0 ± 0.5 ) % , Σ A FB = ( 3.9 ± 4.0 ± 0.6 ) % , Δ A FB = ( 3.1 ± 4.0 ± 0.4 ) % , where the first uncertainty is statistical and the second systematic. The results are consistent with the conservation of C P symmetry and the Standard Model expectations. © 2025 CERN, for the LHCb Collaboration 2025 CERN

Aaij, R. (ORCID:0000000305331952)↗
DOE OSTI2024

Precision Measurements of the Neutron Magnetic Form Factor to High Momentum Transfer using Durand's Method

Protons and neutrons, collectively known as nucleons, along with electrons, constitute the funda- mental building blocks of the visible universe. Understanding their internal structure is crucial for addressing key scientific questions about our origin and existence. Elastic electron-nucleon scatter- ing provides insights into the spatial distributions of charge and current within nucleons through their electromagnetic form factors. Accurate knowledge of these form factors over a broad range of Q2, the squared four-momentum transfer in the scattering process, reveals details about the nucleon’s internal structure. However, high-Q2 data of the nucleon electromagnetic form factor is scarce due to the challenges associated with such measurements. This thesis reports preliminary results from high-precision measurements of the neutron magnetic form factor (Gn M ) to unprecedented Q2 using Durand’s method, also known as the “ratio” method. Systematic errors are greatly reduced by extracting Gn M from the ratio of neutron-coincident (D(e, e'n)) to proton-coincident (D(e, e'p)) quasi-elastic electron scattering from deuteron. The scattered electrons were detected in the BigBite spectrometer, which features multiple Gas Elec- tron Multiplier (GEM) layers with large active area for high-precision tracking at very high rates. Simultaneous nucleon detection was performed by the Super BigBite spectrometer, which utilizes a dipole magnet with large solid angle acceptance at forward angles and a novel hadron calorimeter with very high and comparable detection efficiencies for both protons and neutrons. This setup could handle very high luminosity, making high-Q2 measurements feasible. Data were collected at five Q2 points: 3, 4.5, 7.4, 9.9, and 13.6 (GeV/c)2. Preliminary results are reported for all, with the lowest two Q2 points in good agreement with existing world data, while the higher points significantly extend the Q2 range in which Gn M is known accurately. The precision of the highest Q2 point is expected to remain unmatched for years to come.

Datta, Provakar↗
DOE OSTI2023

Revealing the closed pore formation of waste wood-derived hard carbon for advanced sodium-ion battery

Abstract Although the closed pore structure plays a key role in contributing low-voltage plateau capacity of hard carbon anode for sodium-ion batteries, the formation mechanism of closed pores is still under debate. Here, we employ waste wood-derived hard carbon as a template to systematically establish the formation mechanisms of closed pores and their effect on sodium storage performance. We find that the high crystallinity cellulose in nature wood decomposes to long-range carbon layers as the wall of closed pore, and the amorphous component can hinder the graphitization of carbon layer and induce the crispation of long-range carbon layers. The optimized sample demonstrates a high reversible capacity of 430 mAh g −1 at 20 mA g −1 (plateau capacity of 293 mAh g −1 for the second cycle), as well as good rate and stable cycling performances (85.4% after 400 cycles at 500 mA g −1 ). Deep insights into the closed pore formation will greatly forward the rational design of hard carbon anode with high capacity.

25 ENERGY STORAGE↗
DOE OSTI2023

Pressure Transient Analyses and Poroelastic Modeling of Hydraulic Fracture Dilation for Multiple Injections at the Devine Fracture Pilot Site

Our team has conducted electromagnetic (EM) surveys for the past six years to monitor hydraulic-fracture behavior at the Devine Fracture Pilot Site (DFPS). The sub-horizontal orientation of a shallow hydraulic fracture at the DFPS provides uniform access to the fracture area for interrogation and data collection. Ahmadian et al. (2023) suggested a possible correlation between spatiotemporal changes in the flow rate, bottomhole pressure (BHP), and the observed surface recorded electric field at the DFPS. In this paper, we present the development of poroelastic forward models and pressure transient analyses (PTAs) to support the development of a multiphysics inverse model for these EM surveys. First, we conducted PTAs of the shut-in periods after six injections out of 10 to determine the fracture closure pressure (FCP) or the overburden pressure used in a poroelastic fracture reopening model. Second, we developed a finite-element poroelastic model throughout five injection cycles to include the effect of the cumulative injected volumes due to the previous injections on current fracture dilation in the presence of highly permeable unpropped and propped zones adjacent to the cohesive layer that models fracture reopening. Fracture reopening in this poroelastic model is based on a calibrated traction-separation response using the bottomhole pressure collected in two injection campaigns in 2020 and 2022. We used the outcomes of a previous simulation study of the primary hydraulic-fracturing stimulation to define the dimension of an unpropped fracture zone ahead of the propped fracture area. The PTAs led to FCPs consistent with those obtained using the injection data collected at the DFPS in 2020. Further, these analyses showed that at later injections, the fracture closure occurred at a later time with respect to the shut-in time, inferring the effect of cumulative injected volumes in previous injections. The simulation results show that considering the propped and unpropped fracture zones improves our poroelastic model in predicting the injection-well BHP. The numerical simulation results demonstrate a significant excess pore pressure near the fracture because of the preceding formation loadings by the previous injections. The obtained fracture dilation area and fluid pressure distribution provide a basis to improve the development of a multiphysics inverse model. Furthermore, in an iteratively coupled scheme, this pressure distribution can be introduced into EM models to render a holistic view of the causative mechanisms for the surface signal anomalies.

02 PETROLEUM↗
DOE OSTI2022

Lossy checkpoint compression in full waveform inversion: a case study with ZFPv0.5.5 and the overthrust model

This paper proposes a new method that combines checkpointing methods with error-controlled lossy compression for large-scale high-performance full-waveform inversion (FWI), an inverse problem commonly used in geophysical exploration. This combination can significantly reduce data movement, allowing a reduction in run time as well as peak memory. In the exascale computing era, frequent data transfer (e.g., memory bandwidth, PCIe bandwidth for GPUs, or network) is the performance bottleneck rather than the peak FLOPS of the processing unit. Like many other adjoint-based optimization problems, FWI is costly in terms of the number of floating-point operations, large memory footprint during backpropagation, and data transfer overheads. Past work for adjoint methods has developed checkpointing methods that reduce the peak memory requirements during backpropagation at the cost of additional floating-point computations. Combining this traditional checkpointing with error-controlled lossy compression, we explore the three-way tradeoff between memory, precision, and time to solution. We investigate how approximation errors introduced by lossy compression of the forward solution impact the objective function gradient and final inverted solution. Empirical results from these numerical experiments indicate that high lossy-compression rates (compression factors ranging up to 100) have a relatively minor impact on convergence rates and the quality of the final solution.

58 GEOSCIENCES↗
DOE OSTI2023

Evidence for Modification of $b$ Quark Hadronization in High-Multiplicity $pp$ Collisions at $\sqrt{s}$ = 13 TeV

The production rate of $B^{0}_{s}$ mesons relative to $B^0$ mesons is measured by the LHCb experiment in $pp$ collisions at a center-of-mass energy $\sqrt{s}$ = 13 TeV over the forward rapidity interval 2 < $y$ < 4.5 as a function of the charged particle multiplicity measured in the event. Evidence at the 3.4$σ$ level is found for an increase of the ratio of $B^{0}_{s}$ to $B^0$ cross sections with multiplicity at transverse momenta below 6 GeV /$c$, with no significant multiplicity dependence at higher transverse momentum. Comparison with data from $e^+$ $e^–$ collisions implies that the density of the hadronic medium may affect the production rates of $B$ mesons. This is qualitatively consistent with the emergence of quark coalescence as an additional hadronization mechanism in high-multiplicity collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Hyphae move matter and microbes to mineral microsites: Integrating the hyphosphere into conceptual models of soil organic matter stabilization

Associations between soil minerals and microbially derived organic matter (often referred to as mineral-associated organic matter or MAOM) form a large pool of slowly cycling carbon (C). Here, the rhizosphere, soil immediately adjacent to roots, is thought to control the spatial extent of MAOM formation because it is the dominant entry point of new C inputs to soil. However, emphasis on the rhizosphere implicitly assumes that microbial redistribution of C into bulk (non-rhizosphere) soils is minimal. We question this assumption, arguing that because of extensive fungal exploration and rapid hyphal turnover, fungal redistribution of soil C from the rhizosphere to bulk soil minerals is common, and encourages MAOM formation. First, we summarize published estimates of fungal hyphal length density and turnover rates and demonstrate that fungal C inputs are high throughout the rhizosphere–bulk soil continuum. Second, because colonization of hyphal surfaces is a common dispersal mechanism for soil bacteria, we argue that hyphal exploration allows for the non-random colonization of mineral surfaces by hyphae-associated taxa. Third, these bacterial communities and their fungal hosts determine the chemical form of organic matter deposited on colonized mineral surfaces. Collectively, our analysis demonstrates that omission of the hyphosphere from conceptual models of soil C flow overlooks key mechanisms for MAOM formation in bulk soils. Moving forward, there is a clear need for spatially explicit, quantitative research characterizing the environmental drivers of hyphal exploration and hyphosphere community composition across systems, as these are important controls over the rate and organic chemistry of C deposited on minerals.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Life Cycle Assessment and Techno-Economic Assessment of Lithium Recovery from Geothermal Brine

Lithium-ion batteries (LIB) play an essential role in the electrification of the transportation sector, and battery demand for lithium compounds will see a significant increase in the coming decades. This has raised concerns on the supply of lithium, and as a result, technologies are being developed to process unconventional lithium sources. One promising technology is to extract lithium from geothermal brine using lithium-aluminum-layered double hydroxide chloride (LDH) sorbent and forward osmosis. A combined life cycle assessment (LCA) and techno-economic assessment (TEA) is conducted to evaluate the environmental and economic performance of this technology. It is assumed that the lithium extraction unit is an add-on to a 50 MW geothermal power plant located in California. The analysis is based on lab-scale experimental data and stoichiometry while considering the economy of scale for an industrial system. LCA results suggest that, compared with conventional LiOH and Li 2 CO 3 production pathways, the new technology achieves 1–95% reduction in environmental impacts. Even higher reduction can be achieved for LiOH produced via electrolysis. This add-on unit for lithium extraction could achieve a payback period of less than 1 year and reach net present values of $454M and $315M and internal rates of return of 792 and 1130% for LiOH and Li 2 CO 3 production pathways, respectively. The favorable environmental and economic performance suggests that LDH sorption coupled with forward osmosis has great potential to enable the domestic production of battery lithium compounds and that further development should be carried out.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Commercial, industrial, and institutional discount rate estimation for efficiency standards analysis: Sector-level data 1998–2022

Underlying each of the U.S. Department of Energy’s (DOE’s) federal appliance and equipment energy conservation standards are a set of complex analyses of the projected costs and benefits of regulation. Any new or amended standard must be designed to achieve significant additional energy conservation, provided that it is technologically feasible and economically justified (42 U.S.C. 6295(o)(2)(A)). DOE determines economic justification based on whether the benefits exceed the burdens, considering a variety of factors, including the economic impact of the standard on consumers of the product and the savings in lifetime operating cost compared to any increase in price or maintenance expenses (42 U.S.C. 6295(o)(2)(B)). As part of this determination, DOE conducts a life-cycle cost (LCC) analysis, which models the combined impact of appliance first cost and operating cost changes on a representative commercial building sample to identify the fraction of customers achieving LCC savings or incurring net cost at the considered efficiency levels. Thus, the commercial discount rate value(s) used to calculate the present value of energy cost savings within the LCC model implicitly plays a role in estimating the economic impact of potential standard levels. This report provides an in-depth discussion of the commercial discount rate estimation process. It is an update to previous reports on estimating commercial discount rates from firm-level and sector-level financial data (e.g., Fujita, 2021, 2016). Major topics covered in this report include the following: -Discount rate estimation methods and rationale -Data sources used and data limitations -Discount rate distributions for use in standards analysis -Discount rate estimation methods and distributions specific to the small business subgroup analysis A version of this analysis was most recently released in 2022. Going forward, this report will be updated as data allow and analyses necessitate.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

High event rate analysis technique for the dual-axis duo-lateral position-sensitive silicon detectors of FAUST

The dual-axis duo-lateral (DADL) position-sensitive silicon detector was developed to obtain precise position and energy information for detected charged particles. The Forward Array Using Silicon Technology (FAUST) is currently equipped with 68 DADL detectors backed by CsI(Tl) scintillators for the study of charged particle correlations in heavy-ion collisions where precise position and energy information is essential. When conventional signal processing electronics were used for the DADL detectors, a position dependence of the measured energy as well as distortions in the calculated particle positions were observed. In previous work, waveforms from the detector after preamplification were studied to better understand the features that give rise to these distortions; therein, a waveform analysis technique was developed to improve the energy resolution and linearity in position reconstruction. However, the reading and writing of waveforms for an entire detector array limits data collection rates and adds significant burden in data storage and analysis speed. In this work, the integrators of a Struck SIS3316 ADC were utilized to process 228 Th source data to develop and optimize a new analysis method that captures the benefits of the waveform analysis technique while circumventing the waveform writing requirement. This integrator method – capable of 59 keV (FWHM) energy resolution – was used in the collection of 35 MeV/nucleon 28 Si + 12 C collision data using FAUST to investigate exotic decays of highly excited highly deformed nuclei. In this data, a position resolution of 0.4 mm (FWHM) was obtained for 25 MeV α-particles; for α-particles near this energy that originate from 8 Be ground state decays, a 8 Be ground state width of 30 keV (FWHM) was obtained. The impact of the energy-dependent DADL position resolution emergent from electronic noise on the quality of excited state measurement was modeled and compared to the experimental data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Enhancing photoionization rate calculations in low-temperature plasmas using spectral methods

Photoionization plays a central role in the development of streamer discharges and other non-equilibrium plasma phenomena. It creates seed electrons, which are essential for positive streamer propagation, allowing the ionization front to move forward. Because of this, accurate modeling of photoionization is very important for predicting streamer behavior and plasma evolution. The photoionization process in air (N 2 – O 2 mixture) is often described by the Zheleznyak model (1982). This model is usually solved through Helmholtz-type equations that approximate the Zheleznyak photoionization model (Zheleznyak et al. 1982) as Partial Differential Equations (PDEs). Conventional numerical methods, such as the Finite Difference Method (FDM) or Finite Volume Method (FVM), are widely used to solve these equations. Although they are prevalent, the computational cost of these methods due to their need for matrix operations and iterative solver is demanding. To address this challenge, this work develops a spectral solver based on the Fast Fourier Transform (FFT) combined with Discrete Cosine Transform (DCT) and Discrete Sine Transform (DST) to calculate the photoionization rate efficiently in an axisymmetric cylindrical domain. This method naturally satisfies the boundary conditions used in the model and converts the PDE into algebraic ones in spectral space. Thus, avoids the need for iterative matrix solvers. When compared with FDM results, it is demonstrated that the new solver not only maintains accuracy, but also reduces the computational cost, showing a performance increase of approximately 100 compared to FDM over a wide range of problem sizes. The method is parallelized using Message Passing Interface (MPI) and has been integrated into a fluid plasma model for streamer simulation. Here, this FFT-based approach provides a fast and reliable alternative for calculating photoionization in fluid models, helping large-scale plasma simulations run faster and efficiently, and allows higher-resolution simulation without extra computational cost.

Axisymmetric system↗
DOE OSTI2024

CCUS 2024, Interpreting the strain tensor Larry Murdoch Interpreting strain tensor data to characterize and monitor reservoirs for CO2 storage and other applications

Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and this has opened the door to new opportunities for characterization and monitoring during CCUS. We have demonstrated this method by deploying strainmeters at shallow depths (30 to 40m) and then conducting injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. We have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS. Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and our objective was to evaluate opportunities for strain monitoring during characterization and monitoring for CCUS. Our approach was to deploy strainmeters at shallow depths (30 to 40m) and then conduct injection well tests in an underlying reservoir at 530m depth. The results indicate that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. In conclusion, we have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS.

Murdoch, Larry↗
DOE OSTI2024

Renewable low-Z wall for fusion reactors with built-in tritium recovery (Final Technical Report)

This project pursued development of a novel renewable plasma-facing wall technology for fusion reactors. The technology is based on a slurry which can be easily delivered by delivery tubes to the reactor wall. The slurry dries at the hot reactor wall into pebble rods which are extruded out into the hot plasma, where the pebbles break off and fall along the reactor wall and can be recovered by gravity and re-used. The falling pebbles carry away heat and tritium and also protect the wall against large scale erosion or redeposition of material. The research focused on carbon-based pebble rods and demonstrated that pebble rods could be produced from slurry with tolerable levels of outgassing on a reactor-relevant timescale (< 5 minutes). Steady-state handling of reactor relevant (up to 50 MW/m 2 ) normal-incidence heat loads was demonstrated. Pebble release velocities were found to be sufficiently small (< 1 m/s) to allow recovery below the vacuum chamber. Tunability of the pebble rod breaking rate was demonstrated by changing the fill fraction of the interpebble matrix which binds the pebbles together. This work could benefit the public by helping move forward the design of commercially viable fusion energy reactors. Designing a first wall for magnetic fusion reactors which can handle the huge heat loads present and also avoid buildup of tritium-containing deposits is extremely challenging and requires novel approaches like the one being investigated here.

36 MATERIALS SCIENCE↗