Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “integral experiment”

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

SATS Transient Fission Gas Release Test with Irradiated Fuel under Loss of Coolant Accident Conditions

The transient fission gas release (tFGR) during the temperature ramp associated with a loss-of-coolant accident (LOCA) in light-water reactors (LWRs) is likely a significant contribution to the total pressure in a fuel rod and may cause an unexpected rod burst. The lack of data related to tFGR continues to be a key gap in understanding LWR cladding burst behavior under LOCA conditions. To fully characterize this behavior, tFGR data must be collected from several different systems that can capture all relevant testing conditions. Oak Ridge National Laboratory (ORNL) has developed a system to measure the integral tFGR from irradiated fuel segments. This system was designed to integrate with the existing Severe Accident Test Station (SATS) and to build upon decades of experience capturing fission gas to characterize fuel behavior. The tFGR system consists of a sweep gas system to transport gases from the in-cell SATS apparatus to an out-of-cell fission gas detection system composed of a series of cold traps to capture the off-gas from the heating tests and a gamma spectrometry system to detect and measure 85 Kr. Initial system testing operations were completed during which 85 Kr collection and measurement were verified along with the ability to detect stable inert gases. A tFGR test with a high-burnup fuel specimen was successfully conducted by the in-cell SATS-tFGR system at ORNL. The posttest examination is under way, the result of which will be reported in FY24.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Grid Operator Analytics and Assessment Tools for Inverter- Based Resources Dominated Grid (GOAAT-IBR) Project Update

This presentation provides an update on the OPTIMA GOAAT project, with emphasis on the cloud-native data platform developed in-house to ingest, manage, and operationalize high-resolution power system data. Since our last NASPI presentation, accessible via OSTI ID #2671437, the project team advanced the design and deployment of a scalable architecture capable of handling both synchronized and non-synchronized streams, including PMU, point-on-wave (POW), COMTRADE, and SCADA data. These materials review the project status, recent progress, and key lessons learned. The core of the presentation examines the architecture and engineering of our cloud-native ingestion and data management platform. We then explain how pipelines were designed to collect, normalize, time-align, store, and serve heterogeneous data at scale. We will discuss design choices such as data models, streaming versus batch ingestion, storage tiers, and interoperability with analytics applications. Practical experiences with cloud-native technologies were shared during the event, including benefits, limitations, and integration challenges in a utility environment, along with methods used to improve performance, reduce latency, and optimize resource usage. The presentation also showcases user interface designs and visualization tools that convert raw measurements and analytics results into intuitive, actionable insights for operators and engineers. During the presentation examples were provided demonstrating how visualization, event views, and summarized analytics enhance situational awareness and support operational decision-making. These use cases illustrate how a well-designed data infrastructure can bridge the gap between high-volume measurements and practical grid operations.

Aminifar, Farrokh↗

Sensitivity Studies of Premixed Flame Response to Transverse, High Frequency Disturbances

This study analyzes the response of premixed flames to high frequency transverse acoustic disturbances. A significant body of experiment, modeling and simulation efforts has been performed to understand the response of the spatially integrated heat release of premixed flames to low-frequency disturbances. These studies quantified this sensitivity through the Flame Transfer Function (FTF), as the global heat release is the key quantity influencing thermoacoustic stability for acoustically compact flames. Recent efforts have expanded these models to consider high-frequency acoustic disturbances, and the sensitivity of premixed flames through the Rayleigh integral (RI). The full Rayleigh Integral, as opposed to simply the spatially integrated heat release, must be considered in these cases where the flame is not acoustically compact. Prior modeling studies considered several effects - stability of axisymmetric flames centered in a cylindrical combustion chamber and the effect of flame offset from the combustor axis and the effect of mean flame asymmetries in the flame centered in the combustor. These studies showed that the flame's stability was sensitive to all these parameters, as well as to the transverse mode in question. In this paper, we quantify these different parametric sensitivities of premixed flame stability in the presence of high-frequency transverse acoustic modes - specifically sensitivity to flame length and nozzle location within the acoustic mode shape for 2 different transverse modes. The sensitivity of the RI to these control parameters provides an improved understanding of which flame parameters are critical to improve the overall thermoacoustic stability of the system.

Acharya, Vishal S.↗

Identifying Urban Pluvial Frequency Flooding Hotspots Using the Topographic Control Index and Remote Sensing Radar Images for Early Warning Systems

Identifying areas that frequently experience post-rainfall ponding is essential for effective flood mitigation and planning. This study integrates Sentinel-1 radar imagery and the Topographic Control Index (TCI) to identify 378 flood-prone urban depressions in Beaumont, Texas. Out of 159 major rainfall events, only six had Sentinel-1 radar imagery acquired within six hours of peak rainfall, and these were used to generate the flood frequency map; the ground-based flood sensor data were used to verify that these selected events corresponded to actual peak rainfall and to validate radar-detected water pixels. Validation results showed 100% precision, 70.87% recall, an F1-score of 82.95%, and 71.32% overall accuracy. Approximately 84% of medium-to-high TCI depressions overlapped with Beaumont’s two-year inundation map, confirming a strong relationship between TCI and observed flooding. A total of 124 depressions retained significant water, and after excluding 25 engineered detention ponds, 99 natural depressions remained flood vulnerable. Among these, 74 depressions with medium or high TCI were identified as the highest-priority nuisance flooding hotspots. The results demonstrate that combining TCI with radar imagery provides a reliable and cost-effective approach for identifying areas prone to frequent urban ponding. This framework supports practical decision-making for drainage improvements, hotspot identification, and early-warning system development in urban flood-prone regions.

Sentinel-1 radar imagery↗

Long-range transported continental aerosol in the eastern North Atlantic: three multiday event regimes influence cloud condensation nuclei

Abstract. The eastern North Atlantic (ENA) is a region dominated by pristine marine environment and subtropical marine boundary layer clouds. Under unperturbed atmospheric conditions, the regional aerosol regime in the ENA varies seasonally due to different seasonal surface-ocean biogenic emissions, removal processes, and meteorological regimes. However, during periods when the marine boundary layer aerosol in the ENA is impacted by particles transported from continental sources, aerosol properties within the marine boundary layer change significantly, affecting the concentration of cloud condensation nuclei (CCN). Here, we investigate the impact of long-range transported continental aerosol on the regional aerosol regime in the ENA using data collected at the U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) user facility on Graciosa Island in 2017 during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) campaign. We develop an algorithm that integrates number concentrations of particles with optical particle dry diameter (Dp) between 100 and 1000 nm, single scattering albedo, and black carbon concentration to identify multiday events (with duration >24 consecutive hours) of long-range continental aerosol transport in the ENA. In 2017, we detected nine multiday events of long-range transported particles that correspond to ∼ 7.5 % of the year. For each event, we perform HYSPLIT 10 d backward trajectories analysis, and we evaluate CALIPSO aerosol products to assess, respectively, the origins and compositions of aerosol particles arriving at the ENA site. Subsequently, we group the events into three categories, (1) mixture of dust and marine aerosols, (2) mixture of marine and polluted continental aerosols from industrialized areas, and (3) biomass burning aerosol from North America and Canada, and we evaluate their influence on aerosol population and cloud condensation nuclei in terms of potential activation fraction and concentrations at supersaturation of 0.1 % and 0.2 %. The arrival of plumes dominated by the mixture of dust and marine aerosol in the ENA in the winter caused significant increases in baseline Ntot. Simultaneously, the baseline particle size modes and CCN potential activation fraction remained almost unvaried, while cloud condensation nuclei concentrations increased proportionally to Ntot. Events dominated by a mixture of marine and polluted continental aerosols in spring, fall, and winter led to a statistically significant increase in baseline Ntot, a shift towards larger particular sizes, a higher CCN potential activation fractions, and cloud condensation nuclei concentrations of >170 % and up to 240 % higher than during baseline regime. Finally, the transported aerosol plumes characterized by elevated concentration of biomass burning aerosol from continental wildfires detected in the summertime did not statistically contribute to increase baseline aerosol particle concentrations in the ENA. However, particle diameters were larger than under baseline conditions, and CCN potential activation fractions were >75 % higher. Consequentially, cloud concentration nuclei concentrations increased by ∼ 115 % during the period affected by the biomass burning events. Our results suggest that, through the year, multiday events of long-range continental aerosol transport periodically affect the ENA and represent a significant source of CCN in the marine boundary layer. Based on our analysis, in 2017, the multiday aerosol plume transport dominated by a mixture of dust and marine aerosol, a mixture of marine and polluted continental aerosols, and biomass burning aerosols caused increases in the NCCN baseline regime of, respectively, 6.6 %, 8 %, and 7.4 % at SS 0.1 % (and, respectively, 6.5 %, 8.2 %, and 7.3 % at SS 0.2 %) in the ENA.

54 ENVIRONMENTAL SCIENCES↗

Refactoring the elastic–viscous–plastic solver from the sea ice model CICE v6.5.1 for improved performance

This study focuses on the performance of the elastic–viscous–plastic (EVP) dynamical solver within the sea ice model, CICE v6.5.1. The study has been conducted in two steps. First, the standard EVP solver was extracted from CICE for experiments with refactored versions, which are used for performance testing. Second, one refactored version was integrated and tested in the full CICE model to demonstrate that the new algorithms do not significantly impact the physical results. The study reveals two dominant bottlenecks, namely (1) the number of Message Parsing Interface (MPI) and Open Multi-Processing (OpenMP) synchronization points required for halo exchanges during each time step combined with the irregular domain of active sea ice points and (2) the lack of single-instruction, multiple-data (SIMD) code generation. The standard EVP solver has been refactored based on two generic patterns. The first pattern exposes how general finite differences on masked multi-dimensional arrays can be expressed in order to produce significantly better code generation by changing the memory access pattern from random access to direct access. The second pattern takes an alternative approach to handle static grid properties. The measured single-core performance improvement is more than a factor of 5 compared to the standard implementation. The refactored implementation of strong scales on the Intel® Xeon® Scalable Processors series node until the available bandwidth of the node is used. For the Intel® Xeon® CPU Max series, there is sufficient bandwidth to allow the strong scaling to continue for all the cores on the node, resulting in a single-node improvement factor of 35 over the standard implementation. This study also demonstrates improved performance on GPU processors.

58 GEOSCIENCES↗

Integrated GW Farm ABM

This Data Repository includes data used for the integrated groundwater- farm ABM model, raw model output from scenario ensemble, and processed outputs that isolate the groundwater storage depletion outcomes for the 35,000 farm cells. Model Inputs: Farm ABM Inputs: This folder contains the input data used by the integrated groundwater - farm ABM modelling script (Python file) used for the high performance computing (HPC) experiments. The sub-folder "data inputs" contains all of the farm attribute data, while the three files in the folder have the hydrogeological data lookup table (NLDAS Cost Curve Attributes.csv), a lookup table (Theis well function table.csv) for the groundwater cost curve function, and the farm indexes and corresponding NLDAS ids for all of the cells run in this experiment (nldas farms subset final.csv). NLDAS Cost curve hydrogeological data: Hydrogeological data aggregated to 1/8 degree resolution and aligned with the NLDAS grid. Parameters include: water depth below ground surface [meters], subsurface porosity [unitless], aquifer depth from ground surface to aquifer bottom [meters], annual average recharge (USGS: mm, Doll: meters), and three different hydraulic conductivity (K) values (meters/day). The three K values represent the mean value from Gleeson et al. (2018), one standard deviation above the mean from Gleeson et al. (2018), and the de Graaf et al. 2020 modifications to certain lithologies. Additional information about these datasets and their processing are documented in the supplement to Yoon et al. 2025 (in review). Output: Raw outputs: This folder contains a .zip file that has model outputs for the entire scenario ensemble. There is one csv for each farm id, using the format "farm farmid cases.csv". The relationship between the farm id and NLDAS id is defined by the "nldas farms subset final.csv" located in the Farm ABM Inputs folder. Each csv has 625 rows, corresponding to 625 combinations of different scenario parameter values. Each row (scenario) represents the outcome of a 100 year simulation. Columns define scenario settings and summary statistics for each scenario. The first four columns define the scenario settings: "hydro ratio," "econ ratio," "K scenario," and "gamma scenario." The hydro and econ ratios are values passed to the modeling script that influence multipliers for other model parameters, as documented in the supplement to Yoon et al. 2025 (in review). The gamma multiplier is a coefficient multiplier applied to the baseline gamma values (values below 1 represent lower unobserved costs compared to baseline, values above 1 represent higher costs). The K scenario names represent K values of: "low": 0.5 m/d, "int 1": 2.5 m/d, "int 2": 10 m/d, "high": 50 m/d, and "gleeson": mean Gleeson K value. "Perc vol depleted" is the fraction of groundwater depleted at the end of the 100 simulation. Processed Output: Derived depletion outcomes from raw outputs: All of the individual csv files from the Raw outputs were aggregated into a single file that has the scenario settings and fraction depletion "Perc vol depleted" for every farm cell, for every scenario. The other two files define relationships between the farm id, NLDAS id, and local and major aquifer units, used for aquifer-level depletion analysis.

Agent based modeling↗

Introduction of a Variable Inductance Transformer for the Design of Resonant Power Converters

Magnetic integration is a hot topic in power electronics that concerns the use of a transformer’s leakage and magnetizing inductances purposefully in isolated power electronic converters, thereby giving the opportunity to save the cost and footprint of any additional inductor. This is of prime interest, especially in CLLLC resonant converters which require up to three inductors. For a complete integration of these inductances, the concept of a variable inductance transformer (VIT) is introduced in this thesis. A VIT is an adaptive magnetic structure that facilitates an easy adjustment of both magnetizing and leakage inductances to meet their desired values. However, for a more promising design, an accurate estimation of these inductances is necessary. While the evaluation of magnetizing inductance is quite straightforward, the calculation of leakage inductance is rather convoluted, because the leakage inductance is influenced by both the winding layout and the operating frequency. In this thesis, three new semi-analytical methods for calculating the frequency-independent leakage inductance, and a novel semi-analytical method for evaluating the frequency-dependent leakage inductance are proposed. These methods can calculate the respective leakage inductances of a VIT within an outstanding ±5% uncertainty. Finally, a bidirectional CLLLC resonant dc-dc converter is investigated for the constant current constant voltage (CCCV) charging of the next-generation 900 V traction battery of an electric vehicle. A new voltage gain equation is derived for designing the CLLLC resonant tank, and a small-signal model is presented for designing the variable-frequency feedback controller. Furthermore, a new methodology to design a VIT is developed to overcome the challenges associated with small coupling coefficients and guarantee a complete magnetic integration of the tank inductances. All theoretical results presented herein are verified through simulations and experiments performed on hardware prototypes designed in the lab.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Extended magnetohydrodynamics simulations of thin-foil Z-pinch implosions with comparison to experiments

Cylindrical foil liners, with foil thicknesses on the order of 400 nm, are often used in university-scale Z-pinch experiments (~1 MA in 100 ns) to study physics relevant to inertial confinement fusion efforts on larger-scale facilities (e.g., the MagLIF effort on the 25-MA Z facility at Sandia National Laboratories). The use of ultrathin foil liners typically requires a central support rod to maintain the structural integrity of the liner target assembly prior to implosion. The radius of this support rod sets a limit on the maximum convergence ratio achievable for the implosion. In recent experiments with a support rod and a preimposed axial magnetic field, helical instability structures in the imploding foil plasma were found to persist as the foil plasma stagnated on the rod and subsequently expanded away from the rod. We have now used the 3D extended magnetohydrodynamics simulation code PERSEUS (which includes Hall physics) to study these experiments. The results suggest that it is the support rod which is responsible for the helical structures persisting beyond stagnation. Furthermore, we find that as the radius of the support rod decreases (i.e., as the convergence ratio increases), the integrity and persistence of the helical modes diminish. In the limit with no support rod, we find that the structure of the final stagnation column is governed by the structure of the central precursor plasma column. Furthermore, these simulation results and their comparisons to experiment are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Lead tungstate calorimeter of the Jefferson Lab Eta Factory experiment

A new electromagnetic calorimeter (ECAL) consisting of 1596 lead tungstate PbWO 4 scintillating crystals has been fabricated and installed in the experimental Hall D at Jefferson Lab (JLab). The highgranularity, high-resolution calorimeter is required by the JLab Eta Factory experiment, whose main physics goal is to study rare decays of eta mesons. The ECAL replaced the inner part of the forward lead glass calorimeter of the GlueX detector. Signals from the detector will be digitized using twelve-bit flash analog-to-digital converters operated at a sampling rate of 250 MHz. The ECAL is integrated into the trigger system of the GlueX detector using electronics modules designed at JLab. The ECAL is currently at the commissioning stage and should be ready for the physics run in January 2025. We will give an overview of the JEF experiment, the design and construction of the ECAL, and the integration of the detector and its infrastructure into the Hall D experimental setup.

Somov, Alexander [Thomas Jefferson National Accele↗

Multivariable degradation modeling and life prediction using multivariate fractional Brownian motion

In system prognostics and health management, multivariable degradation models have been widely developed to predict the life of complex systems using degradation data of multiple Performance Characteristics (PCs). Recent studies have detected a Long-Term Memory (LTM) effect among the degradation process of various PCs, implying a strong coupling phenomenon between the future degradation behavior and historical degradation trajectory. Although the LTM has been widely integrated into single-PC-based degradation modeling, it has not been considered in multi-PC-based scenarios. To capture LTM among multiple PCs, this article proposes a novel LTM-integrated Multivariate Degradation Model (MDM) for system life prediction based on multivariate fractional Brownian motion, which simultaneously incorporates the cross-correlation among different PCs. To estimate parameters of the LTM-integrated MDM, a maximum likelihood method is developed. Here, two likelihood-ratio hypothesis tests are developed to test the existence of the overall and individual LTM effect among multiple PCs. Both simulation studies and physical experiments on the performance degradation of solar energy conversion and storage devices are conducted to validate the proposed model. Results reveal that the proposed LTM-integrated MDM significantly outperforms existing MDMs in life prediction, while the lifetime uncertainty is heavily underestimated by those traditional approaches that neglect the LTM.

42 ENGINEERING↗

Automatic Generation of High-Performance Convolution Kernels on ARM CPUs for Deep Learning

In this work, we present FastConv, a template-based code auto-generation open source library that can automatically generate high-performance deep learning convolution kernels of arbitrary matrices/tensors shapes. FastConv is based on the Winograd algorithm, which is reportedly the highest performing algorithm for the time-consuming convolution layers of convolutional neural networks. ARM CPUs cover a wide range designs and specifications, from embedded devices to HPC-grade CPUs. The leads to the dilemma of how to consistently optimize Winograd-based convolution solvers for convolution layers of different shapes. FastConv addresses this problem by using templates to auto-generate multiple shapes of tuned kernels variants suitable for skinny tall matrices. As a performance portable library, FastConv transparently searches for the best combination of kernel shapes, cache tiles, scheduling of loop orders, packing strategies, access patterns, and online/offline computations. Auto-tuning is used to search the parameter configuration space for the best performance for a given target architecture and problem size. The experiments with layer-wise evaluation on the VGG--16 model confirms a 1.25x performance gains is got by tuning the Winograd library. Integrated comparison results shows 1.02x to 1.40x, 1.14x to 2.17x, and 1.22x and 2.48x speedup is achieved over NNPACK, Arm NN, and FeatherCNN on the Kunpeng 920 beside few cases. Furthermore, problem size performance portability experiments with various convolution shapes shows that FastConv achieves 1.2x to 1.7x speedup and 2x to 22x speedup over NNPACK and ARM NN inference engine using Winograd on Kunpeng 920 . CPU performance portability evaluation on the VGG--16 show an average speedup over NNPACK of 1.42x, 1.21x, 1.26x, 1.37x, 2.26x, and 11.02x is observed on Kunpeng 920, Snapdragon 835, 855, 888, Apple M1, and AWS Graviton2, respectively.

97 MATHEMATICS AND COMPUTING↗

Partitioned exponential methods for coupled multiphysics systems

Multiphysics problems involving two or more coupled physical phenomena are ubiquitous in science and engineering. This work develops a new partitioned exponential approach for the time integration of multiphysics problems. After a possible semi-discretization in space, the class of problems under consideration is modeled by a system of ordinary differential equations where the right-hand side is a summation of two component functions, each corresponding to a given set of physical processes. The partitioned-exponential methods proposed herein evolve each component of the system via an exponential integrator, and information between partitions is exchanged via coupling terms. Here, the traditional approach to constructing exponential methods, based on the variation-of-constants formula, is not directly applicable to partitioned systems. Rather, our approach to developing new partitioned-exponential families is based on a general-structure additive formulation of the schemes. Two method formulations are considered, one based on a linear-nonlinear splitting of the right hand component functions, and another based on approximate Jacobians. The paper develops classical (non-stiff) order conditions theory for partitioned exponential schemes based on particular families of T-trees and B-series theory. Several practical methods of third order are constructed that extend the Rosenbrock-type and EPIRK families of exponential integrators. Several implementation optimizations specific to the application of these methods to reaction-diffusion systems are also discussed. Numerical experiments reveal that the new partitioned-exponential methods can perform better than traditional unpartitioned exponential methods on some problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

BISON Validation to In-Situ Cladding Burst Test and High Burnup LOCA Experiments

The process to develop and qualify nuclear fuels for commercial nuclear application requires fundamental material development, characterization, and design; out-of-pile testing on unirradiated materials; integral fuel rod irradiations, testing, and post irradiation examinations (PIEs); and transient analyses. The historical approach depends on the generation of large empirical datasets and series of integral fuel rod irradiations, and this approach ultimately takes approximately 20 years—or sometimes longer—to acquire data through extensive sequential testing. Thus, the qualification and eventual deployment of new fuel systems constitute a long and drawn-out process. However, recent technological advancements have provided researchers the opportunity to perform out-of-cell, in-situ measurements to assess material performance for the duration of the experiment. One such example of this capability is the use of digital image coordination and thermal imaging to assess Zircaloy cladding performance under a simulated loss-of-coolant accident (LOCA) transient condition. In general, the in-situ measurements provide high-fidelity strain, strain rates, and temperature surface maps. This is critical for the United States nuclear industry, which is actively developing a technical basis to support extending the peak rod average burnup from 62 to ~75 GWd/tU and the deployment of accident-tolerant fuel. However, the Nuclear Regulatory Commission, through their research information letter, outlined a number of technical issues for the industry to address before extending burnup. One topic of specific interest is understanding the cladding balloon and rupture geometry during the LOCA heatup phase. Leveraging these advanced in-situ capabilities, this work intends to use insitu data generated from a simulated LOCA in the Severe Accident Test Station at Oak Ridge National Laboratory to better understand high-temperature creep and its impact on Zircaloy balloon and rupture performance. This work used the BISON fuel performance code to compare BISON’s high-temperature creep model predictions to in-situ data and identify limitations and gaps within the model. Additionally, BISON will subsequently be used to simulate relevant high-burnup LOCA tests. The results will be analyzed and compared to the available post-test data in a manner consistent with the approach outlined in the Nuclear Regulatory Commission research information letter.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Integrating DOE ASCR Computing into HEPCloud through GlideinWMS

Fermilab's HEPCloud facility expands the laboratory's computing capacity by provisioning resources beyond the local grid, using GlideinWMS to deliver pilots to where experiments such as CMS and DUNE run. The High-Performance Computing (HPC) facilities of the DOE Office of Advanced Scientific Computing Research (ASCR) are a growing part of that pool. HEPCloud currently provisions NERSC over SSH, but NERSC is moving away from that path as it adopts multi-factor authentication and directs automated access to its Superfacility API and the DOE Integrated Research Infrastructure (IRI) APIs. Maintaining and extending access across the ASCR ecosystem now requires provisioning through these interfaces. This work adds new pilot submission paths to GlideinWMS for the NERSC Superfacility API, IRI, and Globus Compute. Each uses the provisioning model GlideinWMS already applies to batch resources, so experiments can run on ASCR computing resources without changes to their existing workflows. This work finally presents a comparison of the paths to guide which interfaces are best suited for different workflows.

Majumder, Meghanto [U. Houston (main)]↗

Investigation of the Effect of Injected CO 2 on the Morrow B Sandstone through Laboratory Batch Reaction Experiments: Implications for CO 2 Sequestration in the Farnsworth Unit, Northern Texas, USA

About one million tons of CO 2 have been injected into the Farnsworth unit to date. The target reservoir for CO 2 injection is the Morrow B Sandstone, which is primarily made of quartz with lesser amounts of albite, calcite, chlorite, and clay minerals. The impact of CO 2 injection on the mineralogy, porosity, and pore water composition of the Morrow B Sandstone is a major concern. Although numerical modeling studies suggest that porosity changes will be minimal, significant alterations to mineralogy and pore water composition are expected. Given the implications for CO 2 storage effectiveness and risk assessment, it is crucial to verify the accuracy of theoretical model predictions through laboratory experiments. To this end, batch reaction experiments were conducted to model conditions near an injection well in the Morrow B Sandstone and at locations further away, where the CO 2 has been diluted by formation water. The laboratory experiments involved submerging thin sections of both coarse- and fine-grained facies of the Morrow B Sandstone in formation water samples with varying levels of CO 2 . The experiments were conducted at the reservoir temperature of 75 °C. Two experimental runs were conducted, one lasting for 61 days and the other for 72 days. The initial fluid composition used in the second run was the same as in the first. The mineralogy changes in the thin sections were analyzed using SEM and the Tescan Integrated Mineral Analyzer (TIMA), while changes in the composition of the formation water were determined using ICP-AES. During each experiment, a thin layer of white fine-grained particles consisting mainly of dolomite and silica formed on the surface of the thin sections, leading to significant reductions in Ca, Mg, and Sr in the formation water. This outcome is consistent with numerical model predictions that dolomite would be the primary mineral that would react with injected CO 2 and that silica would be oversaturated in the formation water. Changes in mineral abundance in the thin sections themselves were much less systematic than in the theoretical modeling experiments, perhaps reflecting heterogeneities in the mineral grain size surface area to volume ratios and mineral distributions in the thin sections not considered in the numerical models.

58 GEOSCIENCES↗

Report on FY 2023 Research and Development on Specially Designed Creep-fatigue Experiments on Alloy 617 in Support of Improving Creep-fatigue Evaluation Approaches

Experimental and numerical studies in support of developing the integrated Elastic–Perfectly Plastic (EPP) plus Simplified Model Test (SMT) design methodology, referred to as the EPP+SMT method, continued in FY 2023. This report focuses on the methods for extrapolating the EPP+SMT creep-fatigue life curves at long hold times and low strain ranges at elevated temperatures. In this work, the available uniaxial creep-fatigue failure data on Alloy 617 at temperatures of 950°C and 850°C were analyzed to provide a guidance on the development of the extrapolation method and the creep-fatigue failure criteria. A viscoplastic constitutive model for Alloy 617 was adopted to extrapolate the mechanical responses to low strain ranges and long hold times. A set of design curves of Alloy 617 at temperatures of 950°C, 850°C, and 800°C with tensile hold times of 1 hr, 100 hr, and 1,000 hr are developed. Furthermore, creep-fatigue testing on two notch specimen geometries, shallow-notch and sharp V-notch, on Alloy 617 was performed 950°C to understand the multiaxial stress relaxation behavior. The experimental and numerical results on the notch specimens were compared with those on the standard uniaxial smooth bar specimens. The effect of multi-axial stress state combined with elastic follow-up on the stress relaxation behavior was investigated in this report.

36 MATERIALS SCIENCE↗

Enabling CO 2 Isothermal Compression Using Liquid Piston and Integrated Gas Cooler

This report outlines the development of a near-isothermal liquid piston compressor. Initially, a proof-of-concept experiment was conducted to validate the feasibility and potential of near-isothermal compression. A first-of-its-kind near-isothermal liquid piston compressor was then successfully built and tested. Additionally, a Matlab-based transient model was developed to capture performance trends accurately. These efforts laid the groundwork for developing a second-generation prototype with enhanced efficiency and design. Subsequent field tests at the Helix Innovation Center in Ohio demonstrated a cooling capacity of 2 kW with a COP of 1.64 at an ambient temperature of 25°C. An improved model was also created to suggest different heat exchanger geometries for better compression efficiency.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗