Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Operator approximation”

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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A buffering heat exchanger/thermal energy storage system for desalination applications

A buffering heat exchanger/thermal energy storage (BHXTES) system was designed using a fast-response phase-change material/graphite foam (PCM/GF) medium for desalination and potentially other industrial applications. Unlike a store-now-and-use-later TES system, this study focused on developing a dynamically balanced, continuously used, system that optimized three functions. A lab-scale prototype was designed and fabricated with a PureTemp 151/GF storage medium. The experimental data showed excellent thermal performance repeatability indicating minimal effects of PCM expansion and contraction on the GF as well as negligible PCM redistribution effects. The experimental data also served as a validation of the numerical model of the BHXTES system, after which optimization of the BHXTES system was conducted through numerical simulations for both PureTemp 151/GF and solar salt/GF as the storage medium. The results show that, with a Therminol 55 heat transfer fluid, continuous operation of the system can be realized with a solar energy source through cycles of 8-h changing for desalination and energy storage with the supplied heat followed by 16-h discharging for desalination with the stored energy. As a result, the working fluid outlet temperatures are generally in the ranges of approximately 128–160 °C for a PureTemp 151/GF medium and 198–220 °C for a solar salt/GF medium.

Buffering↗
DOE OSTI2021

Automatic Differentiation in MetaPhysicL and Its Applications in MOOSE

Efficient solution via Newton’s method of nonlinear systems of equations requires an accurate representation of the Jacobian, corresponding to the derivatives of the component residual equations with respect to the degrees of freedom. In practice these systems of equations often arise from spatial discretization of partial differential equations used to model physical phenomena. These equations may involve domain motion or material equations that are complex functions of the systems’ degrees of freedom. Computing the Jacobian by hand in these situations is arduous and prone to error. Finite difference approximations of the Jacobian or its action are prone to truncation error, especially in multiphysics settings. Symbolic differentiation packages may be used, but often result in an excessive number of terms in realistic model scenarios. An alternative to symbolic and numerical differentiation is automatic differentiation (AD), which propagates derivatives with every elementary operation of a computer program, corresponding to continual application of the chain rule. Automatic differentiation offers the guarantee of an exact Jacobian at a relatively small overhead cost. In this work, we outline the adoption of AD in the Multiphysics Object Oriented Simulation Environment (MOOSE) via the MetaPhysicL package. We describe the application of MOOSE’s AD capability to several sets of physics that were previously infeasible to model via hand-coded or Jacobian-free simulation techniques, including arbitrary Lagrangian-Eulerian and level-set simulations of laser melt pools, phase-field simulations with free energies provided through neural networks, and metallic nuclear fuel simulations that require inner Newton loop calculation of nonlinear material properties.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Operator Relaxation and the Optimal Depth of Classical Shadows

Classical shadows are a powerful method for learning many properties of quantum states in a sample-efficient manner, by making use of randomized measurements. Here we study the sample complexity of learning the expectation value of Pauli operators via “shallow shadows,” a recently proposed version of classical shadows in which the randomization step is effected by a local unitary circuit of variable depth t. Here we show that the shadow norm (the quantity controlling the sample complexity) is expressed in terms of properties of the Heisenberg time evolution of operators under the randomizing (“twirling”) circuit—namely the evolution of the weight distribution characterizing the number of sites on which an operator acts nontrivially. For spatially contiguous Pauli operators of weight k, this entails a competition between two processes: operator spreading (whereby the support of an operator grows over time, increasing its weight) and operator relaxation (whereby the bulk of the operator develops an equilibrium density of identity operators, decreasing its weight). From this simple picture we derive (i) an upper bound on the shadow norm which, for depth t~log⁡(k), guarantees an exponential gain in sample complexity over the t=0 protocol in any spatial dimension, and (ii) quantitative results in one dimension within a mean-field approximation, including a universal subleading correction to the optimal depth, found to be in excellent agreement with infinite matrix product state numerical simulations. Our Letter connects fundamental ideas in quantum many-body dynamics to applications in quantum information science, and paves the way to highly optimized protocols for learning different properties of quantum states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Zero-crosstalk silicon photonic refractive index sensor with subwavelength gratings

Abstract Silicon photonic index sensors have received significant attention for label-free bio and gas-sensing applications, offering cost-effective and scalable solutions. Here, we introduce an ultra-compact silicon photonic refractive index sensor that leverages zero-crosstalk singularity responses enabled by subwavelength gratings. The subwavelength gratings are precisely engineered to achieve an anisotropic perturbation-led zero-crosstalk, resulting in a single transmission dip singularity in the spectrum that is independent of device length. The sensor is optimized for the transverse magnetic mode operation, where the subwavelength gratings are arranged perpendicular to the propagation direction to support a leaky-like mode and maximize the evanescent field interaction with the analyte space. Experimental results demonstrate a high wavelength sensitivity of − 410 nm/RIU and an intensity sensitivity of 395 dB/RIU, with a compact device footprint of approximately 82.8 μm 2 . Distinct from other resonant and interferometric sensors, our approach provides an FSR-free single-dip spectral response on a small device footprint, overcoming common challenges faced by traditional sensors, such as signal/phase ambiguity, sensitivity fading, limited detection range, and the necessity for large device footprints. This makes our sensor ideal for simplified intensity interrogation. The proposed sensor holds promise for a range of on-chip refractive index sensing applications, from gas to biochemical detection, representing a significant step towards efficient and miniaturized photonic sensing solutions. Graphical Abstract

Materials Science↗
DOE OSTI2020

Best Practices for Shielding Analyses of Activated Metals and Spent Resins from Reactor Operation

This report investigated best practices for performing shielding evaluations of Type B waste packages, as defined in 10 CFR Part 71 regulations on packaging and transportation of radioactive material, or packages for which the contents are not defined before loading and may include a broad range of nuclides, geometries, and non-fuel materials. The following non-fuel waste streams were analyzed: Activated metals from decommissioned commercial power reactors, including Type 304 stainless steel, reactor vessel steel, and Inconel, Control blades from boiling water reactors, Neutron-activated corrosion products on surfaces of activated metals, and Spent resins from power plant operations. Measured elemental compositions, including major constituents and impurities, for steel and Inconel samples from commercial power reactors were used in activation calculations to determine radionuclide inventories in activated metals. For a simplified cask model, 60 Co contribution to the total external package dose rate at 30 days after shutdown varied from approximately 60% to 95%, depending on the activated metal, initial cobalt impurity concentration in the metal, and the thickness of the overpack gamma shield. Its maximum contribution to the total external package dose rate of approximately 100% was reached within the time interval of 2 to 5 years after shutdown and was maintained for up to 45 to 60 years after shutdown, depending on material, initial cobalt impurity concentration, and shield thickness. Thereafter, the 60 Co contribution to external package dose rate decreased with increasing decay time. Cobalt-60 is primarily produced by neutron reactions with the cobalt impurity in steel and Inconel. Other important radionuclides in activated metals contributing to package external dose rate are radionuclides with relatively short decay times, including 51 Cr, 59 Fe, 58 Co, and 54 Mn. These radionuclides may be represented as an equivalent 60 Co activity/source because 60 Co gamma ray emissions are bounding in terms of source strength and energy to other important radionuclides identified in the analyzed activated metals. Approaches for modeling the neutron-activated corrosion products that may be attached to activated reactor components were analyzed in this report. It was demonstrated that a surface source is more conservative than a uniform volumetric source for the treatment of neutron-activated corrosion products with respect to external package dose rates. An analysis of the maximum radionuclide loadings reported on spent resins identified the radionuclides 137 Cs, 60 Co, 134 Cs, 65 Zn, and 58 Co as the primary contributors to external package dose rate. For a resin cooled for 3.08 years, the external package dose rate was entirely produced by the reported 137 Cs and 60 Co inventory. The neutron sources from actinides found on spent resins or activated metals produced negligible dose rates and may be ignored in dose rate analyses. Effects of idealized waste material, source geometry, and spatial material/source distributions on external package dose rates were determined based on dose rate results for a simplified cask model under normal conditions of transport. Type 304 stainless steel, zirconium, and aluminum with adjusted mass densities based on a maximum content weight were analyzed for material modeling effects on external package dose rate. These materials produced identical external package dose rates within the statistical uncertainties of the dose rate estimates. Among four different source geometry configurations with homogeneous material of different mass densities, uniform volumetric source distribution, and the same source strength, the geometry configuration with lowest mass density (i.e., minimum self-shielding effect among the four cases) was most conservative. Spatial source distributions that better represent localized peak 60 Co activity values were more conservative than a uniform volumetric source distribution, assuming the same weight and total source strength per package. The increase in external dose rate caused by localized source peaks can be as much as the ratio between source peak activity density to average activity density, depending on the location of the activated metal with peak activity density. Therefore, the shielding analysis may be simplified if localized peak activities and the average activity per package can be measured/determined and documented at the time of cask loading. For simplicity, external package dose rates may be determined based on average source activity and a uniform volumetric source distribution. The dose rate results from that calculation model multiplied by the ratio of peak activity density to average activity density will produce maximum dose rate values for conservative estimates.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Energy and economic impacts of U.S. federal energy and water conservation standards adopted from 1987 through 2024

This paper presents estimates of the key impacts of U.S. national energy and water conservation standards adopted from 1987 through 2024. The standards for consumer products, commercial and industrial equipment, lighting products, and plumbing products include those set by legislation as well as standards adopted by the Department of Energy (DOE) through rulemaking. In 2024 alone, these standards are estimated to have saved an estimated 6.0 quads of primary energy, which is equivalent to 6.5% of total U.S. energy consumption, and 1.7 trillion gallons of water, which is equivalent to approximately 12% of the annual water withdrawals for public supply in the U.S in 2015. The estimated reduction in CO 2 emissions associated with the standards in 2024 was 270 million metric tons, which is equivalent to 5.6% of total annual U.S. CO 2 emissions from energy consumption. The annual savings in operating costs for households and businesses totaled $\$$105 billion, and the average household saved $\$$576 in operating costs as a result of standards on residential appliances and plumbing products. The estimated cumulative past and future energy and water savings from these standards amount to 307 quads of energy and almost 53 trillion gallons of water. The estimated cumulative CO 2 emissions reduction from the standards come to 10.3 billion metric tons. Accounting for the increased upfront costs of more-efficient products and the energy and water cost savings over the products’ lifetime, the standards have a cumulative net present value of benefit of $\$$3.2 trillion using 3 percent discount rate and $\$$3.5 trillion using 7 percent discount rate when discounting past and future benefits to 2024.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

FPGA-Based Spill Regulation System for the Muon Delivery Ring at Fermilab

The Muon to Electron Experiment (\acrshort{Mu2e}) requires a uniform beam profile from the Muon Delivery Ring to meet their experimental needs. A specialized Spill Regulation System (\acrshort{SRS}) has been developed to help achieve consistent spill uniformity. The system is based on a custom-designed carrier board featuring an Arria 10 SoC, capable of executing real-time feedback control. The FPGA processes beam pulses of approximately 200 ns every 1.695 microseconds, allowing for continuous monitoring of the extracted spill intensity through fast bunch integration. The system directly controls three quadrupole magnets, which work in conjunction with sextupole magnets to achieve third-order resonant extraction. Furthermore, the board interfaces with Fermilab s Accelerator Control Network (ACNET), enabling operators to modify spill regulation settings in real-time via the control network while providing diagnostic waveforms. These waveforms help operators monitor the process and fine-tune the feedback mechanisms. This paper presents an overview of the board's architecture and its initial progress toward regulating beam extraction. This initial version of the regulation system aims to evaluate baseline performance to inform future system improvements.

Berlioz, Jose Rene [Fermilab]↗
DOE OSTI2021

Developing Practical Models of Complex Salts for Molten Salt Reactors

Molten salt reactors (MSRs) utilize salts as coolant or as the fuel and coolant together with fissile isotopes dissolved in the salt. It is necessary to therefore understand the behavior of the salts to effectively design, operate, and regulate such reactors, and thus there is a need for thermodynamic models for the salt systems. Molten salts, however, are difficult to represent as they exhibit short-range order that is dependent on both composition and temperature. A widely useful approach is the modified quasichemical model in the quadruplet approximation that provides for consideration of first- and second-nearest-neighbor coordination and interactions. Its use in the CALPHAD approach to system modeling requires fitting parameters using standard thermodynamic data such as phase equilibria, heat capacity, and others. A shortcoming of the model is its inability to directly vary coordination numbers with composition or temperature. Another issue is the difficulty in fitting model parameters using regression methods without already having very good initial values. The proposed paper will discuss these issues and note some practical methods for the effective generation of useful models.

Besmann, Theodore M. (ORCID:0000000155980550)↗
DOE OSTI2022

Geophysical Characterization and Monitoring of the 200 Area Treated Effluent Disposal Facility to Support Permit Renewal; TEDF Geophysical Monitoring

The Treated Effluent Disposal Facility (TEDF), located in the 200 East Area of the Hanford Site, is a site where non-hazardous and non-radioactive liquid wastes are disposed into two state-permitted infiltration basins. In 2016, the Washington State Department of Ecology denied a permit renewal request for TEDF due to the inability to adequately assess the impact of TEDF discharge water on the underlying groundwater quality. TEDF overlies the relatively impermeable Ringold Lower Mud (RLM) unit, whose upper contact lies in the vadose zone approximately 30 m below ground surface and approximately 10 m above the water table boundary. The RLM is assumed to isolate TEDF discharge water from the natural groundwater aquifer, which is monitored using wellbores screened below the RLM. Therefore, samples collected from monitoring wells near TEDF are not considered representative of TEDF discharge water. Rather, TEDF discharge water is assumed to mound on top of the RLM to form a perched aquifer. To support permitting of the TEDF, a new monitoring well is required that can be used to sample water from the presumed perched water aquifer above the RLM. Ideally, the screened section of the well would be located at the peak of the perched water mound(s), which presumably occurs at the point of maximum vertical flux from TEDF to the RLM, or equivalently where the dominant infiltration flow paths reach the RLM. This report describes how time-lapse 3D electrical resistivity tomography (ERT) was used in conjunction with nominal TEDF discharge operations to image the dominant flow paths from each pond to the RLM. Results are summarized in Figure ES.1. Figure ES.1A shows a satellite image of the TEDF overlain by an array of surface ERT electrodes. The solid and dashed black circles denote the zones of maximum vertical flux at the RLM within the south and north ponds, respectively, and presumably the regions where perched water peaks during discharge. Figure ES.1B shows time-lapse difference imaging results approximately 15 days after switching discharge from the south pond to the north pond. Blue iso-surfaces beneath the north pond denote regions of increasing bulk electrical conductivity caused by increasing saturation due to infiltrating water. Red iso-surfaces beneath the south pond denote regions that were previously saturated during the south pond discharge and are now de-saturating, causing a decrease in bulk electrical conductivity. In both cases, the zones of maximum change mark the dominant flow paths to the RLM. Figure ES.1C shows time-lapse imaging results approximately 15 days after switching discharge from the north pond to the south pond. In this case, blue iso-surfaces mark the dominant flow paths to the RLM from the south pond. The red iso-surfaces mark the dominant flow paths that existed during discharge to the north pond. The regions of maximum change in bulk conductivity (due to saturation or desaturation) that mark the primary flow paths are equivalent in both cases and denoted by the dashed solid and black circles. If a perched water zone forms on the RLM, it is likely to mound within or near the dashed circle during discharge to the north pond, and in the solid circle during discharge to the south pond. In other words, if perched water mounding occurs, the circles mark the optimum locations suggested by the ERT imaging for monitoring boreholes to be placed, enabling samples to be collected that are representative of TEDF discharge impacts on groundwater quality.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Distinguishing fissile uranium isotopes using an active well neutron coincidence counter

Proposed thorium-based nuclear fuel cycles are likely to require quantification and verification of 233 U within nuclear material. Because of their similar fission cross sections, active neutron nondestructive assay (NDA) systems may respond similarly to 233 U and 235 U. Traditional safeguards equipment has been optimized for 235 U and 238 U quantification associated with conventional uranium/plutonium fuel cycles and may not be directly applicable to 233 U quantification when mixed with other actinides. This work used models of the large volume active well coincidence counter (LV-AWCC) at Oak Ridge National Laboratory to evaluate the performance of this neutron NDA system to differentiate fissile uranium isotopes. The models were developed to simulate NDA system performance in response to a number of triangular radiation signature training device sources within the central cavity or well. This work predicted that the LV-AWCC can effectively differentiate 233 U from 235 U in certain modes of operation. In active mode, the LV-AWCC with the cadmium liner results in different doubles count rates between the fissile isotopes for a given fissile uranium mass. Without the cadmium liner, the uranium isotopes provide a statistically indistinguishable doubles count rate response for the fissile masses considered in this work (up to approximately 150 g). The cadmium liner serves to harden the neutron interrogation spectrum, which better exploits the notable difference in the 233 U and 235 U fission cross sections at approximately 1 eV. In passive mode, the two fissile isotopes exhibit different doubles and singles count rates regardless of liner presence because the passive source strength of 233 U is approximately 2 orders of magnitude stronger than that of 235 U due to the shorter half-life and correspondingly higher (α, n) yield. We conclude that using neutron interrogation in the LV-AWCC, two measurements are needed to quantify 233 U content in mixed uranium items. The first measurement is used to determine the total fissile uranium mass using a mode that cannot distinguish fissile isotopes (i.e., where a similar response is observed for both fissile uranium isotopes such as active doubles without cadmium or using a thermal neutron interrogation source). In conclusion, the second measurement is used to determine the 233 U content by using a differentiating technique (e.g., passive doubles, passive doubles to singles ratio, active doubles with cadmium).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Understanding Inlet Concentration Effects on the Electrocatalytic Conversion of CO 2 to Formic Acid in Gas-Fed Electrolyzers

The electrochemical CO 2 reduction reaction (CO2RR) to produce value-added products remains a developing technology for utilizing waste CO 2 streams. Most device-level CO2RR studies use pure CO 2 gas feeds; however, the effect of dilute CO 2 on the electrolyzer performance is an important consideration for large-scale electrolyzer operation, single-pass conversion, and real-world CO 2 source utilization. This work investigates the effect that the CO 2 concentration has on the performance of formic acid (HCOOH) producing tin oxide (SnO 2 ) and bismuth oxide (Bi 2 O 3 ) catalysts in an electrolyzer device setting. Surprisingly, SnO2 demonstrated an approximately 20% increase in HCOOH selectivity (Faradaic efficiency) when the CO 2 concentration decreased from 100 to 20%. In contrast, Bi 2 O 3 consistently demonstrated high selectivity toward HCOOH across the same CO 2 concentration range. The effects of the CO 2 concentration on selectivity were further investigated with half-cell experiments and in situ Raman spectroscopy, which revealed dynamic changes in the cathodic overpotential and chemical state of the catalyst that depended on the CO 2 concentration. Density functional theory calculations showed how changes in the surface oxidation state of Sn, varying from fully oxidized SnO 2 to metallic Sn(0), affect the thermodynamic barriers of the three main observed products: HCOOH, CO, and H 2 . Our results indicate that dilute CO 2 concentrations required larger cathodic overpotentials to sustain a fixed current density, which, in turn, pushed the Sn-based catalyst toward a more reduced surface that was favorable to HCOOH formation. On the other hand, the Bi-based catalyst remained in a metallic state at CO2RR-relevant potentials and demonstrated a consistent product selectivity regardless of CO 2 concentration. These findings highlight how varying the CO 2 inlet gas concentrations affects the chemical state of catalysts and the resulting performance metrics.

42 ENGINEERING↗
DOE OSTI2020

Excited-state electronic structure of molecules using many-body Green’s functions: Quasiparticles and electron–hole excitations with VOTCA-XTP

Here, we present the open-source VOTCA-XTP software for the calculation of the excited-state electronic structure of molecules using many-body Green’s function theory in the GW approximation with the Bethe–Salpeter equation (BSE). This work provides a summary of the underlying theory and discusses the details of its implementation based on Gaussian orbitals, including resolution-of-identity techniques and different approaches to the frequency integration of the self-energy or acceleration by offloading compute-intensive matrix operations using graphics processing units in a hybrid OpenMP/Cuda scheme. A distinctive feature of VOTCA–XTP is the capability to couple the calculation of electronic excitations to a classical polarizable environment on an atomistic level in a coupled quantum- and molecular-mechanics (QM/MM) scheme, where a complex morphology can be imported from Molecular Dynamics simulations. The capabilities and limitations of the GW–BSE implementation are illustrated with two examples. First, we study the dependence of optically active electron–hole excitations in a series of diketopyrrolopyrrole-based oligomers on molecular-architecture modifications and the number of repeat units. Second, we use the GW–BSE/MM setup to investigate the effect of polarization on localized and intermolecular charge-transfer excited states in morphologies of low-donor content rubrene–fullerene mixtures. These showcases demonstrate that our implementation currently allows us to treat systems with up to 2500 basis functions on regular shared-memory workstations, providing accurate descriptions of quasiparticle and coupled electron–hole excited states of various characters on an equal footing.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Precision Reconstruction of Rational Conformal Field Theory from Exact Fixed-Point Tensor Network

The novel concept of entanglement renormalization and its corresponding tensor network renormalization technique have been highly successful in developing a controlled real-space renormalization group (RG) scheme. Numerically approximate fixed-point (FP) tensors are widely used to extract the conformal data of the underlying conformal field theory (CFT) describing critical phenomena. In this paper, we present an explicit analytical construction of the FP tensor for 2D rational CFT. We define it as a correlation function between the “boundary-changing operators” (BCO) on triangles. Our construction fully captures all the real-space RG conditions. We also provide concrete examples, such as Ising, Yang-Lee, and tricritical Ising models, to compute the scaling dimensions explicitly based on the corresponding FP tensor. The BCO descendants turn out to be an optimal basis such that truncation in bond dimensions naturally produces comparable accuracies with the leading existing FP algorithms. Interestingly, our construction of FP tensors is closely related to a strange correlator, where the holographic picture naturally emerges. Our results also open a new door toward understanding CFT in higher dimensions. Published by the American Physical Society 2025

Cheng, Gong (ORCID:0009000891587404)↗
DOE OSTI2022

Impact of Storage Dispatch Assumptions on Resource Adequacy and Capacity Credit

Representing energy-limited resources in power system probabilistic resource adequacy assessment introduces new considerations that invalidate classical modeling assumptions. In particular, such resources have multi-period operating objectives and constraints that in real systems are addressed via a sequence of rolling intertemporal optimizations. Ideally, adequacy models would develop dispatch decisions by solving a similar sequence of problems, but this approach has historically been too computationally intensive for practical use in Monte Carlo simulations, with studies making use of simplifying approximations instead. These simplifications have the potential to distort the assessed value of energy-limited resources on the system.This work describes three classes of storage dispatch assumptions in current use and discusses their theoretical differences. It then provides an empirical analysis of their differences on test systems with different levels of storage, assessing the potential for a study's modeling assumptions to influence the perceived contribution of energy-limited resources.

capacity credit↗
DOE OSTI2025

Cryogenic Design and Thermal Analysis of EIC Central Detector (ePIC) Solenoid Magnet (MARCO)

The Electron Ion Collider (EIC) physics program utilizes a 2.0 T superconducting magnet at the heart of its ePIC detector system. This approximately 3.5 m long and 2.84 m diameter warm bore magnet has a 20 tons cold mass which is conduction-cooled using liquid helium at 4.5 K. A closed loop active thermosiphon system is chosen to facilitate the cooling and to maintain a minimum of 2 K temperature margin on the peak operating temperature (4.7 K) of the superconductor. Here, this paper presents the cryogenic design of the cooling system and the thermal analysis of the solenoid. A Computational fluid dynamics (CFD) model was developed to quantify the performance of the two–phase closed thermosiphon system and predict the temperature gradient on the cold mass.

Gopinath, Sandesh [Thomas Jefferson National Accel↗
DOE OSTI2025

FPGA-Based Spill Regulation System for the Muon Delivery Ring at Fermilab

The Muon to Electron Experiment (Mu2e) requires a uniform beam profile from the Muon Delivery Ring to meet their experimental needs. A specialized Spill Regulation System (SRS) has been developed to help achieve consistent spill uniformity. The system is based on a custom-designed carrier board featuring an Arria 10 SoC, capable of executing real-time feedback control. The FPGA processes beam pulses of approximately 200 ns every 1.695 $μ$s, allowing for continuous monitoring of the extracted spill intensity through fast bunch integration. The system directly controls three quadrupole magnets, which work in conjunction with sextupole magnets to achieve third-order resonant extraction. Furthermore, the board interfaces with Fermilab's Accelerator Control Network (ACNET), enabling operators to modify spill regulation settings in real-time via the control network while providing diagnostic waveforms. These waveforms help operators monitor the process and fine-tune the feedback mechanisms. This paper presents an overview of the board's architecture and its initial progress toward regulating beam extraction. This initial version of the regulation system aims to evaluate baseline performance to inform future system improvements.

Berlioz, J. R. [Fermilab]↗
DOE OSTI2023

Development of Novel Ferritic-Martensitic Steels with Superior Creep Properties for Power Plant Applications

A creep resistant martensitic steel, CPJ-7, was developed with an operating temperature approaching 650°C. Subsequently, another creep resistant martensitic steel, JMP, was designed with the potential to operate at, or slightly above, 650°C. This report describes the development of these alloys from early iterations of CPJ steels to the final CPJ-7 formulation as well as the more heat resistant JMP steel formulation. The design originated from computational modeling for phase stability and precipitate strengthening using fifteen constituent elements. Approximately forty heats of CPJ and ten of JMP, each weighing ~7 kg, were vacuum induction melted. A computationally optimized heat treatment schedule was developed to homogenize the ingots prior to hot forging and rolling prior to final normalization and tempering. Overall, wrought and cast versions of CPJ-7 present superior creep properties when compared to wrought and cast versions of COST alloys for steam turbine and wrought and cast versions of P91/92 for boiler applications. For instance, the Larson Miller Parameter curve for CPJ-7 at 650°C almost coincides with that of COST E at 620°C. The prolonged creep life was attributed to slowing down the process of the destabilization of the MX and M 23 C 6 precipitates at 650°C. On the other hand, the cast version of CPJ-7 also revealed superior mechanical performance, well above commercially available cast 9% Cr martensitic steel or derivatives, especially for creep life at 650°C. The casting process employed slow cooling to simulate the conditions of a thick wall full-size steam turbine casing but utilized a separate homogenization step prior to final normalization and tempering. To advance the development of CPJ-7 for commercial applications, a process was used to scale up the production of the alloy using vacuum induction melting (VIM) and electroslag remelting (ESR), which underlined the importance of melt processing control of intentionally designed minor and trace elements in these advanced alloys. Following the work on CPJ-7, the JMP steels were designed with higher Co for increased solid solution strengthening, Si for oxidation resistance and increased W (with low Mo content) for matrix strength and stability as well as solid solution strengthening. The JMP steels showed increases in creep life compared to CPJ-7 between 118 to 150% at 650°C for testing at various stresses between 138 MPa and 207 MPa. On a Larson-Miller plot, the performance of the JMP steels surpasses that of state-of-the-art MARBN and other MARBN-type steels. The influence of various elements within the composition of the alloys on the microstructure and mechanical properties are discussed. This report presents approximately 420,000 h of in-house creep testing, the equivalent of almost 50 years of cumulative creep tests.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Time-Lapse Electromagnetic Methods for Monitoring Plume Development in a Carbon Storage Reservoir

Conference presentation at International Meeting for Applied Geoscience & Energy (IMAGE), Houston, Texas, August 25–28, 2025. The Energy & Environmental Research Center (EERC) is leading applied research on electromagnetic (EM) monitoring methods at an active carbon storage site in North Dakota. Injection operations at the site began in February 2024, with a permitted injection rate of up to 2.7 million tonnes of CO 2 annually using six injection wells. CO 2 is captured on-site and injected into the Broom Creek Formation, a predominantly sandstone reservoir and saline aquifer located at a depth of approximately 1800 meters. The EERC led acquisition of multiple active- and passive-source EM techniques between August and October of 2024 to provide a thorough understanding of the resistivity profile at the site.

02 PETROLEUM↗