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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Hybrid data-driven and model-informed online tool wear detection in milling machines

Precision machining tool wear is responsible for low product throughput and quality. Monitoring the tool wear online is vital to prevent degradation in machining quality. However, direct real-time tool wear measurement is not practical. This paper presents residual-based anomaly detection models, combining a hybrid model comprised of a physics-based model and a data-driven model (a decision tree or a neural network) to predict signals of interest (e.g., power or forces) under nominal conditions, followed by Page’s cumulative sum test for detecting tool wear on-line using the computer numerical control machine measurements. The most informative features are ranked using dynamic programming and its approximation variants from real-time measurements and machine settings, such as the width of cut, depth of cut, feed rate and spindle speed, that serve as inputs to the predictive models. The baseline nominal model is incrementally updated with experimental data via a gradient boosted adaptation model to generate the residuals that account for discrepancies between the actual machine data under normal conditions and the baseline nominal model predictions. The hybrid model is validated against 20 Mazak milling machine experimental tests and one Haas run-to-failure experiment. The proposed anomaly detector is applied to synthetic data from simulations of the physics-based model at different operating conditions, measurement noise levels, and tool wear levels, and the methods were able to achieve an overall 92% accuracy in data with 1% noise. The anomaly detection methods based on hybrid model reduced the false alarms of either the data-driven or physical-based models alone, and are found to be capable of good online detection of tool wear.

42 ENGINEERING↗
DOE OSTI2020

Max-independent set and the quantum alternating operator ansatz

he maximum-independent set (MIS) problem of graph theory using the quantum alternating operator ansatz is studied. We perform simulations on the Rigetti Forest simulator for the square ring, K 2,3 , and K3,3 graphs and analyze the dependence of the algorithm on the depth of the circuit and initial states. The probability distribution of observation of the feasible states representing maximum-independent sets is observed to be asymmetric for the MIS problem, which is unlike the Max-Cut problem where the probability distribution of feasible states is symmetric. For asymmetric graphs, it is shown that the algorithm clearly favors the independent set with the larger number of elements even for finite circuit depth. Finally, we also compare the approximation ratios for the algorithm when we choose different initial states for the square ring graph and show that it is dependent on the choice of the initial state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Code and Solution Verification Assessment of the CTF Thermal Hydraulic Subchannel Code

CTF is a thermal-hydraulics subchannel code jointly developed by Oak Ridge National Laboratory and North Carolina State University. Over the past seven years, the Consortium for Advanced Simulation of Light Water Reactors (CASL) has made a significant investment in developing CTF so it can be used to model light water reactors, including nominal operating conditions, departure from nucleate boiling analysis, and transients ranging from loss of flow to reactivity insertion accidents. In addition to implementing new modeling capabilities and developing the user input and output interface, extensive work has been performed to improve the code’s quality assurance program, resulting in a development process that conforms with NQA-1 requirements. The CASL program follows the Predictive Capability Maturity Model (PCMM) approach for assessing code quality, which emphasizes performing code verification(ensuring the code converges to the correct answer) and solution verification (ensuring the code converges for the intended application). Code and solution verification are used to identify uncertainty errors introduced by numerical approximations in the code and are important for demonstrating that the model is coded without error, which is an important aspect of the Best Estimate plus Uncertainty method. This paper presents a comprehensive overview of the code and solution verification testing that has been performed on CTF. A top-down approach is taken in which the intended CTF applications are presented, followed by the code features required for their modeling. These features are then linked to the applicable code and solution verification tests that demonstrate proper functioning. Past testing efforts are summarized, and new tests are added to help close gaps in the presented test matrix. Rather than performing “one-off” exercises, these tests are added to the automated CTF regression test suite to ensure continual code quality.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

An ASME-Compliant Helium-4 Evaporation Refrigerator for the SpinQuest Experiment

This paper presents the design, safety basis, and commissioning results of a 1 K liquid helium-4 (4He) evaporation refrigerator developed for the Fermilab SpinQuest Experiment (E1039). The system represents the first high power helium evaporation refrigerator operated in a fixed target scattering experiment at Fermilab and was engineered to comply with the Fermilab ES&H Manual (FESHM) requirements governing pressure vessels, piping, cryogenic systems, and vacuum vessels. The design is mapped to ASME B31.3 (Process Piping) and the ASME Boiler and Pressure Vessel Code (BPVC) for pressure boundary integrity and overpressure protection, with documented compliance to FESHM Chapters 5031 (Pressure Vessels), 5031.1 (Piping Systems), and 5033 (Vacuum Vessels). This work documents the methodology used to reach compliance and approval for the 4He evaporation refrigerator at Fermilab which the field lacks. Design considerations specific to the high radiation target-cave environment including remotely located instrumentation approximately 20 m from the cryostat are summarized, together with the relief-system sizing methodology used to accommodate transient heat loads from dynamic nuclear polarization microwaves and the high-intensity proton beam. Commissioning data from July 2024 confirms that the system satisfies all thermal performance and safety objectives.

Roberts, Jordan D. [Virginia U.]↗
DOE OSTI2021

The diffractive contribution to deep inelastic lepton-proton scattering: Implications for QCD momentum sum rules and parton distributions

The cross section for deep inelastic lepton-proton scattering (DIS) ℓp→ℓ'p'X includes a diffractive deep inelastic (DDIS) contribution ℓp→ℓ'p'X, in which the proton remains intact with a large longitudinal momentum fraction x F greater than 0.9 and small transverse momentum. The DDIS events, which can be identified with Pomeron exchange in the t-channel, account for approximately 10% of all of the DIS events. Thus, when one measures DIS, one automatically includes the leading-twist Bjorken-scaling DDIS events as a contribution to the DIS cross section, whether or not the final-state proton p' is detected. In such events, the missing momentum fraction x p ' ~0.9 carried by the final-state proton p' in the DDIS events could be misidentified with the light-front momentum fraction carried by sea quarks or gluons in the protons' Fock structure. As we shall show in this article, the underlying QCD Pomeron-exchange amplitude which produces the DDIS events does not obey the operator product expansion nor satisfy momentum sum rules. Thus we conclude that the quark and gluon distributions measured in DIS experiments will be misidentified, unless the measurements explicitly exclude the DDIS events and that a correct determination of the parton distribution functions (PDFs) derived from the DIS data requires the explicit subtraction of the DDIS contribution from the full DIS cross section.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Initial-State Dependent Optimization of Controlled Gate Operations with Quantum Computer

There is no unique way to encode a quantum algorithm into a quantum circuit. With limited qubit counts, connectivity, and coherence times, a quantum circuit optimization is essential to make the best use of near-term quantum devices. We introduce a new circuit optimizer called AQCEL, which aims to remove redundant controlled operations from controlled gates, depending on initial states of the circuit. Especially, the AQCEL can remove unnecessary qubit controls from multi-controlled gates in polynomial computational resources, even when all the relevant qubits are entangled, by identifying zero-amplitude computational basis states using a quantum computer. As a benchmark, the AQCEL is deployed on a quantum algorithm designed to model final state radiation in high energy physics. For this benchmark, we have demonstrated that the AQCEL-optimized circuit can produce equivalent final states with much smaller number of gates. Moreover, when deploying AQCEL with a noisy intermediate scale quantum computer, it efficiently produces a quantum circuit that approximates the original circuit with high fidelity by truncating low-amplitude computational basis states below certain thresholds. Our technique is useful for a wide variety of quantum algorithms, opening up new possibilities to further simplify quantum circuits to be more effective for real devices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Non-Fouling, Low Cost Electrolytic Coagulation & Disinfection for Treating Flowback and Produced Water for Reuse

Executive Summary This Final Report is composed of two major sections. The first section presents experimental results obtained in the laboratory during Budget Period 1. The second section presents results from the field test conducted during Budget Period 2. Laboratory Results This research investigated a novel electrochemical process for producing a ferric iron coagulant for use in treating flowback and produced water from hydraulic fracturing and oil production operations. The treatment system improves the effectiveness and lowers the cost of coagulation processes using Fe3+ as the coagulant. The electrolytic coagulant generation (ECG) system uses an electrochemical cell to produce acid and base from oilfield brine solutions. The acid is used to dissolve scrap iron to provide a Fe3+ coagulating agent. The base is used to neutralize the treated water. Compared to conventional electrocoagulation (EC), the main advantage of the ECG system is an order of magnitude lower cost for the source of iron. The second advantage over conventional EC is that it can deliver Fe3+ doses greater than 1 mM, since it is not limited by the amount of dissolved oxygen in the water required to oxidize ferrous to ferric iron. The capital costs for conventional EC and the ECG system are similar, but the operational costs for the ECG system are an order of magnitude lower than conventional EC. The combined costs for iron and electrical energy for treating 1 m3 of FPW with 1 mM Fe3+ is estimated to be $0.87 for conventional EC, and $0.087 for the ECG system. The estimated all-in cost for treating FPW with a 2 mM Fe3+ dose is $0.73/m3 ($0.12/bbl). Field Test Results The field test was conducted at the Paul Foster Central Tank Battery (CTB) in Lea County, New Mexico from November 10, 2022 through December 15, 2022. The feed water to the system was produced water from the Tatanka 1H formation. After approximately two weeks of testing, the initial batch of produced water had been treated and no untreated produced water was available. Thus, after this time, the feed water to the system consisted of previously treated water (i.e., recycled water). The recycled water had nearly all colloidal particles removed, and had a much lower alkalinity due to precipitation of carbonate minerals during the first pass through the system. Although the recycled water was not an ideal test solution due to its low particulate concentrations, its lower alkalinity did allow us to identify the main problem with the treatment system. The main problem with the treatment system was caused by the high alkalinity of the initial feed water (5.4 meq/L) that consumed a significant fraction of the electrochemically generated acid. This resulted in pH values exiting the iron contact tank that were too high to dissolve enough iron to effectively treat the produced water. Tests performed with recycled water with lower alkalinity did not have this issue, and dissolved iron concentrations greater than 20 mM could be achieved. One consistent observation was that effluent water from the iron contact tank was always free of particulates, even when fed with circumneutral solutions. This suggests that there is no need to dissolve high concentrations of iron if all the water to be treated is passed through the scrap iron canister. In this case, dissolved O2 and hypochlorous acid can promote sufficient iron corrosion to provide an effective coagulating agent – even in neutral pH water. This solves the problem resulting from highly alkaline produced water. Modifications to the design of the treatment system were made based on the field test results. These modifications will add minimal additional cost, and were tested in bench-scale laboratory experiments. In short-term testing, the modified treatment process was able to remove colloidal FeS particulates to levels below detection. Long-term, steady state testing will be required to determine whether the modified process is suitable for commercial treatment systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Challenging the Curse of Dimensionality in Multidimensional Numerical Integration by Using a Low-Rank Tensor-Train Format

Numerical integration is a basic step in the implementation of more complex numerical algorithms suitable, for example, to solve ordinary and partial differential equations. The straightforward extension of a one-dimensional integration rule to a multidimensional grid by the tensor product of the spatial directions is deemed to be practically infeasible beyond a relatively small number of dimensions, e.g., three or four. In fact, the computational burden in terms of storage and floating point operations scales exponentially with the number of dimensions. This phenomenon is known as the curse of dimensionality and motivated the development of alternative methods such as the Monte Carlo method. The tensor product approach can be very effective for high-dimensional numerical integration if we can resort to an accurate low-rank tensor-train representation of the integrand function. In this work, we discuss this approach and present numerical evidence showing that it is very competitive with the Monte Carlo method in terms of accuracy and computational costs up to several hundredths of dimensions if the integrand function is regular enough and a sufficiently accurate low-rank approximation is available.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Effects of grain size and porosity on cladding failure in high-burnup UO 2 : A sensitivity and uncertainty study

Isotopic taggants are being studied to aid in the provenance assessment of nuclear materials. However, these taggants must be selected such that they do not adversely affect fuel performance during normal operation or accident scenarios. Taggants are known to affect the fuel’s grain size and porosity. In the work described in this paper, the BISON fuel performance code was used to assess the potential effects of taggants (i.e., grain size and porosity) on fuel rod behavior and cladding failure during a high-burnup, large-break loss-of-coolant accident. Here, 281 individual fuel rods from the same reactor core were modeled for a sensitivity study, a parametric study, and uncertainty quantification. The cladding failure predictions often exhibited stochastic behavior. After additional study, it was found that the cladding failure model is highly sensitive to residual error inherent to numerical approximation solvers. Some strategies to mitigate this sensitivity are discussed. The study found no relationship between known taggant effects and cladding failure status. However, taggants were found to affect the time and location of failure in certain rods. In conclusion, future work to continue investigating and validating these findings is briefly discussed.

Doped UO 2↗
DOE OSTI2021

Collective Analysis of Alpha Particle Losses Due to Self-Absorption by Mass Loading on Radioactive Particulate Glass Fiber Filters

In this study, we derived a relationship between filter mass loading and the percent loss during analysis using the mass loading data collected from six previous studies of self-absorption. Components of mass loading include particulate dust, radioactive particulates, and filter material. In a research report published in 1984, Higby calculated a minimum burial depth for an alpha particle to be lost due to absorption (100% loss) of about 3.7 mg/cm 2 based on calculations for the range of 239-Pu alpha particles in glass fiber filters. From there, Higby concluded that a correction factor of 0.85 assumes approximately 15% losses in the count rate of both alpha and beta particles. In 2000, Luetzelschwab et al. recommended assuming a 40% loss at a loading of 3.3 mg/cm 2 and a 28% loss for a loading of 2.3 mg/cm 2 which included the frontal face mass of the filter. More recently, the 100% losses due to absorption were reported to be in the 10 mg/cm 2 range. Presented here is a trinomial relationship method of relating percent loss due to self-absorption to filter mass loading, based on data reported by Higby, Luetzelschwab et al., Huang et al., Barnett et al., Smith et al., and Hogue et al. Under normal operating conditions at the stacks monitored by Effluent Management, the mass loading of sample filters averages 0.09 ± 0.12 (2s) mg/cm 2 (excluding negative values and outliers) and ranges from 0 mg/cm 2 to 0.24 mg/cm 2 . Based on current mass loading results for Effluent Management stack sample filters, the forced-zero trinomial relationship method estimated self-absorption losses of less than 5%. Because American National Standards Institute/Health Physics Society N13.1-2011 guidelines indicate a correction factor should be used when the penetration of radioactive material into the collection media or self-absorption of radiation by the material collected would reduce the count rate by more than 5%, it is possible continued application of a correction factor to the Effluent Management stack samples is no longer necessary. Nevertheless, continuing to assign a correction factor at the 5% threshold (i.e., 0.95) would be a conservative approach.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Enabling control co-design of the next generation of wind power plants

Abstract. Layout design and wake steering through wind plant control are important and complex components in the design and operation of modern wind power plants. They are currently optimized separately, but with more and more computational and experimental studies demonstrating the gains possible through wake steering, there is a growing need from industry and regulating bodies to combine the layout and control optimization in a co-design process. However, combining these two optimization problems is currently infeasible due to the excessive number of design variables and large solution space. In this article, we present a method that enables the coupled optimization of wind power plant layout and wake steering with no additional computational expense than a traditional layout optimization. We developed a geometric relationship between wind turbines to find an approximate optimal yaw angle, bypassing the need for either a nested or coupled wind plant control optimization. It also provides a significant and immediate improvement to wind power plant design by enabling the co-design of turbine layout and yaw control for wake steering. A small co-designed plant shown in this article produces 0.8 % more energy than its sequentially designed counterpart. This additional energy production comes with no additional infrastructure, turbine hardware, or control software; it is simply the outcome of optimizing the turbine layout and yaw control together, resulting in millions of dollars of additional revenue for the wind power plants of the future.

17 WIND ENERGY↗
DOE OSTI2023

Short-Depth QAOA circuits and Quantum Annealing on Higher-Order Ising Models (Rev.2)

The Quantum Alternating Operator Ansatz (QAOA) and Quantum Annealing (QA) are quantum algorithms that are both based on the adiabatic theorem and both have the goal of sampling the optimal solution(s) of combinatorial optimization problems. Quantum annealing has been physically instantiated on D-Wave devices using superconducting flux qubits, and QAOA can be programmed on digital gate-model quantum computers such as the programmable superconducting transmon qubits devices of the IBMQ series, for instance ibm washington. QAOA and QA address the same types of problems, but it is unclear how they will scale to large problem sizes and to larger and higher-fidelity quantum computers. In this article, we present a direct comparison between QAOA, one and two rounds, run on all 127 qubits of ibm washington and QA run on D-Wave Advantage system4.1 and Advantage system6.1. The problems which allow for this comparison are random Ising model problems whose connectivity matches the heavy hexagonal lattice topology of ibm washington and the Pegasus graph connectivity of the two D-Wave devices. We create two classes of problem instances for this comparison: one with higher order terms (ZZZ variable interactions), linear terms, and quadratic terms, and a separate problem type with only linear and quadratic terms. Our QAOA circuits are novel and extremely short depth, with a CNOT depth of 6 per round, which allows whole chip usage of ibm washington’s heavy hexagonal lattice and can be applied to future heavy-hex chips. We also test the effectiveness of the error suppression technique digital dynamical decoupling on the QAOA circuits. The QAOA circuits compiled to ibm washington are composed of several thousand circuit instructions, approximately 3, 000 depending on the details of the circuit, making these some the largest quantum circuits ever executed on a digital quantum processor. QAOA and QA are compared against the classical heuristic algorithm of simulated annealing and all problem instances are exactly solved using CPLEX in order to evaluate which samplers, if any, correctly found the ground state solution(s) of the problem instances. We find that (i) QA outperforms QAOA on all problem instances, (ii) QAOA samples the problems better than random sampling, and (iii) QAOA angle computation exhibits clear parameter concentration across the ensemble of Ising models.

127 Qubits↗
DOE OSTI2025

First-ever deployment of a SiPM-on-tile calorimeter in a collider: a parasitic test with 200 GeV p p collisions at RHIC.

We describe the testing of a prototype SiPM-on-tile iron-scintillator calorimeter at the Relativistic Heavy Ion Collider (RHIC) during its 200 GeV p p run in 2024. The prototype, measuring 20 × 20 cm 2 and 24 radiation lengths in depth, was positioned in the STAR experimental hall, approximately 8 m from the interaction point and 65 cm from the beam line, covering a pseudorapidity range of about 3.1 < η < 3.4. By using the dark current of a reference SiPM as a radiation monitor, we estimate that the prototype was exposed to a fluence of about 10 10 1-MeV n eq /cm 2 . Channel-by-channel calibration was performed in a data-driven way with the signature from minimum-ionizing particles during beam-on conditions. AGeant4 detector simulation, with inputs from thePythia8 event generator, describes measurements of energy spectra and hit multiplicities reasonably well. These results mark the first deployment, commissioning, calibration, and long-term operation of a SiPM-on-tile calorimeter in a collider environment. This experimental campaign will guide detector designs and operational strategies for the ePIC detector at the future EIC, as well as other applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Fe-Coated Optical Fiber for Distributed Corrosion Monitoring in Soil and Aqueous Environment

Natural gas pipeline corrosion represents a substantial cost and safety concern during normal operations. The effective and real-time monitoring of corrosion is important to detect and mitigate pipeline risks before corrosion-related catastrophic events happen. Here, we describe the use of iron (Fe) coated optical fiber sensors (OFS) for distributed corrosion monitoring where Fe acts as a corrosion proxy. By using an optical backscattering reflectometer (OBR), corrosion was monitored based on the increase in the backscattered intensity amplitude of the light being passed as Fe underwent corrosion. The Fe-coated OFSs were prepared with a film thickness between 25–225 nanometers (nm) by an electroless plating approach and corroded by a carbon dioxide (CO2)-saturated acidic electrolyte. The corrosion rate (CR) was approximately 2.5 millimeters (mm)/year for Fe with a film thickness of 30 nm and increased with film thickness. Backscattering-based CRs were supported by visible light transmission measurements, which have been previously demonstrated. Additionally, a correlation between the intermittent transmission and residual Fe film thickness during the corrosion of Fe was established. The corrosion sensor was also tested in soil at > 1 foot depth. Fe-coated OFS was loaded into a Draka cable, which provides mechanical support during fiber deployment and installation into either acidified soil, top-soil, or sandy soil (50/50).

corrosion↗
DOE OSTI2025

Fe-Coated Optical Fiber for Distributed Corrosion Monitoring in Soil and Aqueous Environments

Natural gas pipeline corrosion represents a substantial cost and safety concern during normal operations. The effective and real-time monitoring of corrosion is important to detect and mitigate pipeline risks before corrosion-related catastrophic events happen. Here, we describe the use of iron (Fe) coated optical fiber sensors (OFS) for distributed corrosion monitoring where Fe acts as a corrosion proxy. By using an optical backscattering reflectometer (OBR), corrosion was monitored based on the increase in the backscattered intensity amplitude of the light being passed as Fe underwent corrosion. The Fe-coated OFSs were prepared with a film thickness between 25–225 nanometers (nm) by an electroless plating approach and corroded by a carbon dioxide (CO2)-saturated acidic electrolyte. The corrosion rate (CR) was approximately 2.5 millimeters (mm)/year for Fe with a film thickness of 30 nm and increased with film thickness. Backscattering-based CRs were supported by visible light transmission measurements, which have been previously demonstrated. Additionally, a correlation between the intermittent transmission and residual Fe film thickness during the corrosion of Fe was established. The corrosion was measurable out to > 100 meters (m), which is a significant improvement over our previous work, which showed corrosion sensing at < 10 m. The corrosion sensor was also tested in soil at > 1 foot depth.

corrosion↗
DOE OSTI2025

Fe-Coated Optical Fiber for Distributed Corrosion Monitoring in Soil and Aqueous Environments

Natural gas pipeline corrosion represents a substantial cost and safety concern during normal operations. The effective and real-time monitoring of corrosion is important to detect and mitigate pipeline risks before corrosion-related catastrophic events happen. Here, we describe the use of iron (Fe) coated optical fiber sensors (OFS) for distributed corrosion monitoring where Fe acts as a corrosion proxy. By using an optical backscattering reflectometer (OBR), corrosion was monitored based on the increase in the backscattered intensity amplitude of the light being passed as Fe underwent corrosion. The Fe-coated OFSs were prepared with a film thickness between 25–225 nanometers (nm) by an electroless plating approach and corroded by a carbon dioxide (CO2)-saturated acidic electrolyte. The corrosion rate (CR) was approximately 2.5 millimeters (mm)/year for Fe with a film thickness of 30 nm and increased with film thickness. Backscattering-based CRs were supported by visible light transmission measurements, which have been previously demonstrated. Additionally, a correlation between the intermittent transmission and residual Fe film thickness during the corrosion of Fe was established. The corrosion was measurable out to > 100 meters (m), which is a significant improvement over our previous work, which showed corrosion sensing at < 10 m. The corrosion sensor was also tested in soil at > 1 foot depth.

aqueous environment↗
DOE OSTI2026

Basin-scale analysis of Mokelumne River Formation for multi-well CO 2 injection

Large-scale carbon sequestration will likely require multiple projects injecting CO 2 into the same subsurface formation, raising concerns about safe operation and efficient use of storage capacity. This study evaluates the long-term response of the Mokelumne River Formation in California’s Sacramento Basin to multi-megaton CO 2 injection using three geologic models of the formation and the open-source simulator GEOS. The analysis focuses on three aspects of reservoir performance: (1) average pressure increase and dissolved CO 2 mass after 30 years for varying well counts and injection rates, (2) pressure interference in a multi-well configuration, and (3) dynamic storage capacity with identification of overpressure-prone regions. The results show that average formation pressure increases linearly with injected mass, while CO 2 dissolution exhibits mixed scaling: approximately linear with the number of wells but sublinear with injection rate, indicating that distributing injection across more wells enhances dissolution more effectively than increasing per-well rates. Pressure-interference effects are significant, with lower-permeability conditions delaying their onset but amplifying their magnitude at later times. Dynamic capacity, defined by the first occurrence of pressure exceeding the local overburden-based limit anywhere in the formation, varies across geologic models and assumed overburden pressure gradients. A lower fidelity geologic model predicts nearly twice the storage capacity of the two higher fidelity models, which consistently estimate approximately 1 Gt under the upper-bound overburden pressure gradient considered for the Sacramento Basin. In all model scenarios, overpressure develops away from injection wells, particularly in higher-elevation regions, highlighting the importance of basin-scale modelling for identifying risks beyond the immediate well vicinity.

Basin-scale↗
DOE OSTI2024

Investigation of precooling unit design options in hydrogen refueling station for heavy-duty fuel-cell electric vehicles

Precooling gaseous hydrogen fuel to a cold temperature before refueling a heavy-duty (HD) hydrogen fuel cell electric vehicle (HFCEV) is essential to avoid overheating the vehicle tank, as well as achieving a high state of charge (SOC). Because a large volume of hydrogen is dispensed during each fill, the need for a shorter fill time amplifies the need to precool each load for refueling a HFCEV. Thus, the design and operation of a precooling unit (PCU), as well as the associated capital and operating costs, plays a pivotal role in any plans for heavy-duty hydrogen refueling stations. Here, in this paper, we present a thermodynamic and technoeconomic analysis of a PCU in a gaseous hydrogen refueling station (HRS) for HD HFCEVs. By employing Argonne National Laboratory's hydrogen station cost optimization and performance evaluation (H2SCOPE) model, the refueling of 50 kg of hydrogen on-board type IV tank at ambient temperatures of 15–45 °C and varying fill rates is simulated. The required degree of precooling temperature to obtain either 100% or the maximum possible SOC is obtained from the simulation. Additionally, the simulation results demonstrate that the average flow rate of hydrogen is approximately 40% lower than the maximum flow rate during a typical fill; which motivates further evaluation of the instantaneous hydrogen mass flow rate profile and suggests the scope of improving precooling unit design. Accordingly, a hybrid strategy of precooling hydrogen has been proposed to address the cooling load by sizing the refrigeration unit for the average flow rate of hydrogen, while supplementing the above average peak hydrogen flow cooling load through thermal buffering. The combined technique enables the downsizing of the original PCU capacity by 25–40% and demonstrates a potential cost reduction of the PCU by approximately 30%, which translates to an installed cost reduction of ∼$125,000 per dispenser.

25 ENERGY STORAGE↗