SEARCH · Engineering Papers

Results for “Performance Modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,171 records · Page 65

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Model reduction for hypersonic aerodynamics via conservative LSPG projection and hyper-reduction

High-speed aerospace engineering applications rely heavily on computational fluid dynamics (CFD) models for design and analysis due to the expense and difficulty of flight tests and experiments. This reliance on CFD models necessitates performing accurate and reliable uncertainty quantification (UQ) of the CFD models. However, it is very computationally expensive to run CFD for hypersonic flows due to the fine grid resolution required to capture the strong shocks and large gradients that are typically present. Furthermore, UQ approaches are “many-query” problems requiring many runs with a wide range of input parameters. One way to enable computationally expensive models to be used in such many-query problems is to employ projection-based reduced-order models (ROMs) in lieu of the (high-fidelity) full-order model. In particular, the least-squares Petrov–Galerkin (LSPG) ROM (equipped with hyper-reduction) has demonstrated the ability to significantly reduce simulation costs while retaining high levels of accuracy on a range of problems including subsonic CFD applications. This allows computationally inexpensive LSPG ROM simulations to replace the full-order model simulations in UQ studies, which makes this many-query task tractable, even for large-scale CFD models. This work presents the first application of LSPG to a hypersonic CFD application. In particular, we present results for LSPG ROMs of the HIFiRE-1 in a three-dimensional, turbulent Mach 7.1 flow, showcasing the ability of the ROM to significantly reduce computational costs while maintaining high levels of accuracy in computed quantities of interest.

Fluid Dynamics↗

DOE OSTI2022

High-temperature compressive creep tests of U 3 Si 2 with spark plasma sintering: Experiments and finite element modeling

This paper reports a systematic high-temperature creep behavior of dense U 3 Si 2 pellets using a spark plasma sintering (SPS) apparatus at elevated temperatures and pressures under vacuum conditions. The stress exponent was subsequently derived to be 3.21 at 1173 K and 2.17 at 1223 K, respectively, indicating a grain boundary sliding creep mechanism. The creep activation energy was determined to be 203.6 ± 19.0 kJ/mol, which agrees well with the literature. Finite element modeling was performed using the creep parameters fitted from the strain-time plot. Additionally, the results suggest an excellent match with the experimental data, confirming the validity of the experiments. Microstructure characterizations indicate that the main phase of the specimens after creep tests remains to be U 3 Si 2 , with a 4 μm thick layer of nano-sized particles induced from the diffusion between U 3 Si 2 and alumina disc used to avoid electric current passing through the sample. The successful conduct of creep experiments demonstrates the great potential of SPS to perform high-temperature mechanical testing of nuclear fuels under vacuum conditions. The subsequent finite element modeling exhibits excellent capabilities for accurately predicting material performance in the creep tests and provides a practical tool in evaluating nuclear fuels’ performance for a much-extended time scale.

36 MATERIALS SCIENCE↗

DOE OSTI2022

An Advanced Meso-Scale Peridynamic Modeling Technology using High-Performance Computing for Cost-Effective Product Design and Testing of Carbon Fiber Reinforced Polymer Composites in Light-weight Vehicles

We study a peridynamic composite modeling technology based on the discontinuous Galerkin finite element method, implemented in the commercial LS-DYNA software, for modeling and prediction of failure in carbon fiber reinforced polymer composites. The proposed technology is developed for the material failure analysis at the meso-scale, which provides the prevailing fiber-matrix interaction mechanism, without adoption of the representative volume element method and thus avoiding complicated numerical calibration procedures. Three types of experimental tests—in-plane coupon test, out-of-plane coupon test, and a component crash test—are simulated in a high-performance computing environment to assess the performance of the proposed peridynamic composite modeling technology.

36 MATERIALS SCIENCE↗

DOE OSTI2024

“Understanding Robustness Lottery”: A Geometric Visual Comparative Analysis of Neural Network Pruning Approaches

Deep learning approaches have provided state-of-the-art performance in many applications by relying on large and overparameterized neural networks. However, such networks are very brittle and are difficult to deploy on resource-limited platforms. Model pruning, i.e., reducing the size of the network, is a widely adopted strategy that can lead to a more robust and compact model. Many heuristics exist for model pruning, but our understanding of the pruning process remains limited due to the black-box nature of a neural network model. Empirical studies show that some heuristics improve performance whereas others can make models more brittle. Here, this work aims to shed light on how different pruning methods alter the network’s internal feature representation and the corresponding impact on model performance. To facilitate a comprehensive comparison and characterization of the high-dimensional model feature space, we introduce a visual geometric analysis of feature representations. We evaluated a set of critical geometric concepts decomposed from the commonly adopted classification loss and used them to design a visualization system to compare and highlight the impact of pruning on model performance and feature representation. The proposed tool provides an environment for an in-depth comparison of pruning methods and a comprehensive understanding of how the model responds to common data corruption. By leveraging the proposed visualization, machine learning researchers can reveal the similarities between pruning methods and redundancy in robustness evaluation benchmarks, obtain geometric insights about the differences between pruned models that achieve superior robustness performance, and identify samples that are robust or fragile to model pruning and common data corruption.

Li, Zhimin [Univ. of Utah, Salt Lake City, UT (Uni↗

DOE OSTI2024

Predicting Drug Effects from High-dimensional Asymmetric Drug Data Sets using Graph Neural Networks: A Comprehensive Analysis of Multi-target Drug Effect Prediction

Graph neural networks (GNNs) have emerged as one of the most effective Machine learning (ML) techniques for drug effect prediction from drug molecular graphs. Despite having immense potential, GNN models lack performance when using data sets that contain high dimensional asymmetrically co-occurrent drug effects as targets with complex correlations between them. Training individual learning models for each drug effect and incorporating every prediction result for a wide spectrum of drug effects is beyond practicality. Such an implication provides a testbed to address this challenge as multi-target prediction problems, aiming to predict all drug effects at a time. We develop standard and hybrid graph neural networks (GNNs)to perform two separate tasks that are multi-regression for continuous values and multi-label classification for categorical values contained in our data sets. Since this step makes the target data even more sparse and introduces asymmetric label co-occurrence, the learning of multi-label classification models becomes difficult and heavily impacts the GNN's performance. To address these challenges, we propose a new data oversampling technique to improve multi-label classification performances on all the given imbalanced molecular graph data sets. Using the technique, we improve the data imbalance ratio of the drug effects better than before while protecting the data set's integrity. Finally, we evaluate multi-label classification performance using the best-performant hybrid GNN model on all the oversampled data sets obtained from the proposed oversampling technique. These results outperform those of other ML models including GNN models when they are trained on the original data sets or oversampled data sets using MLSMOTE (a well-known oversampling technique) in all evaluation metrics precision, recall, and F1 score by a significant margin.

Bose, Avishek [ORNL]↗

DOE OSTI2020

Pre-Ionization Controlled Laser Plasma Formation for Ignition Applications (Final Scientific/Technical Report)

This report consists of three grand funded papers. Abstracts for each follow. 1) In this work, we present experimental and modeling studies of air pre-ionization using ultraviolet (UV) laser pulses and its effect on laser breakdown of an overlapped near-infrared (NIR) pulse. Experimental studies are conducted with a 266 nm beam (fourth harmonic of Nd:YAG) for UV pre-ionization and an overlapped 1064 nm NIR beam (fundamental of Nd:YAG), both having pulse duration of 10 ns. Results show that the UV beam produces a pre-ionized volume which assists in breakdown of the NIR beam, leading to reduction in NIR breakdown threshold by factor of >2. Numerical modeling is performed to examine the ionization and breakdown of both beams. The modeled breakdown threshold of the NIR, including assist by pre-ionization, is in reasonable agreement with the experimental results. 2) The present contribution compares the energy absorption, optical emission, temperature, and fluid dynamics of ultraviolet (UV) λ = 266 nm and near infrared (NIR) λ = 1064 nm nanosecond laser induced plasmas in ambient air. For UV pulses at the conditions studied, energy absorption by the plasmas increases relatively gradually with laser pulse energy starting at delivered energy of E~8 mJ. Corresponding measurements of plasma luminosity show that the absorption of UV radiation does not necessarily result in visible plasma emission. For the NIR induced plasmas, the energy absorption profile is far more abrupt and begins at ~55 mJ. In contrast with UV, the absorption of NIR radiation is always accompanied by intense optical emission. The temperatures of both types of plasma have been measured with Rayleigh scattering thermometry (at times after the Thomson signal sufficiently diminishes). The UV plasmas can attain a wider range of temperatures, including lower temperatures, depending on the pulse energy (e.g., T ~400–2000K for E ~7–35 mJ at Δt = 10 μs after the pulse) while the NIR plasmas show only hotter temperatures (e.g., T ~ 12 000 K for E = 75 mJ at Δt = 10 μs after the pulse) as is consistent with the literature. Differences in the fluid dynamics for UV versus NIR pulses are shown with Schlieren imaging. In this work, the contrast in the UV and NIR plasma threshold behavior is attributed to differing roles of avalanche ionization and multiphoton ionization as is also illustrated by a simple numerical model. 3) The present contribution examines the impact of plasma dynamics and plasma-driven fluid dynamics on the flame growth of laser ignited mixtures and shows that a new dual-pulse scheme can be used to control the kernel formation process in ways that extend the lean ignition limit. We perform a comparative study between (conventional) single-pulse laser ignition (λ = 1064 nm) and a novel dual-pulse method based on combining an ultraviolet (UV) pre-ionization pulse (λ = 266 nm) with an overlapped near-infrared (NIR) energy addition pulse (λ = 1064 nm). We employ OH* chemiluminescence to visualize the evolution of the early flame kernel. For single-pulse laser ignition at lean conditions, the flame kernel separates through third lobe detachment, corresponding to high strain rates that extinguish the flame. In this work, we investigate the capabilities of the dual-pulse to control the plasma-driven fluid dynamics by adjusting the axial offset of the two focal points. In particular, we find there exists a beam waist offset whereby the resulting vorticity suppresses formation of the third lobe, consequently reducing flame stretch. With this approach, we demonstrate that the dual-pulse method enables reduced flame speeds (at early times), an extended lean limit, increased combustion efficiency, and decreased laser energy requirements.

Mechanical engineering↗

DOE OSTI2023

High-Fidelity Energy Deposition Ignition Model Coupled with Flame Propagation Models at Engine-like Flow Conditions

With the heightened pressure on car manufacturers to increase the efficiency and reduce the carbon emissions of their fleets, more challenging engine operation has become a viable option. Highly dilute, boosted, and stratified charge, among others, promise engine efficiency gains and emissions reductions. At such demanding engine conditions, the spark-ignition process is a key factor for the flame initiation propagation and the combustion event. From a computational standpoint, there exist multiple spark-ignition models that perform well under conventional conditions but are not truly predictive under strenuous engine operation modes, where the underlying physics needs to be expanded. In this paper, a hybrid Lagrangian-Eulerian spark-ignition (LESI) model is coupled with different turbulence models, grid sizes, and combustion models. The ignition model, previously developed, relies on coupling Eulerian energy deposition with a Lagrangian particle evolution of the spark channel, at every time-step. The spark channel is attached to the electrodes and allowed to elongate at a speed derived from the flow velocity. The LESI model is used to simulate spark ignition in a nonquiescent crossflow environment at engine-like conditions, using converge commercial computational fluid dynamics (CFD) solver. The results highlight the consistency, robustness, and versatility of the model in a range of engine-like setups, from typical with Reynolds-averaged Navier-Stokes (RANS) and a larger grid size to high fidelity with large-eddy simulation (LES) and a finer grid size. The flame kernel growth is then evaluated against Schlieren images from an optical constant volume ignition chamber with a focus on the performance of flame propagation models, such as G-equation and thickened flame model, versus the baseline well-stirred reactor model. Finally, future development details are discussed.

Advanced ignition modeling↗

DOE OSTI2024

Performance Prediction of High‐Entropy Perovskites La 0.8 Sr 0.2 Mn x Co y Fe z O 3 with Automated High‐Throughput Characterization of Combinatorial Libraries and Machine Learning

Perovskite oxides form a large family of materials with applications across various fields, owing to their structural and chemical flexibility. Efficient exploration of this extensive compositional space is now achievable through automated high-throughput experimentation combined with machine learning. In this study, we investigate the composition–structure–performance relationships of high-entropy La 0.8 Sr 0.2 Mn x Co y Fe z O 3±𝞭 perovskite oxides (0 < x, y, z <1; x+y+z≈1) for application as oxygen electrodes in Solid Oxide Cells. Following the deposition of a continuous compositional map using thin-film combinatorial pulsed laser deposition, compositional, structural, and performance properties are characterized using six different techniques with mapping capabilities. Random forests effectively model electrochemical performance, consistently identifying Fe-rich oxides as optimal compounds with the lowest area-specific resistance values for oxygen electrodes at 700 °C. Additionally, the models identify a statistical correlation between oxygen sublattice distortion—derived from spectral analysis of Raman-active modes—and enhanced performance.

high entropy oxides↗

DOE OSTI2023

Mitigation of spatial nonstationarity with vision transformers

Spatial nonstationarity, the location variance of features’ statistical distributions, is ubiquitous in many natural settings. For example, in geological reservoirs rock matrix porosity varies vertically due to geomechanical compaction trends, in mineral deposits grades vary due to sedimentation and concentration processes, in hydrology rainfall varies due to the atmosphere and topography interactions, and in metallurgy crystalline structures vary due to differential cooling. Conventional geostatistical modeling workflows rely on the assumption of stationarity to be able to model spatial features for geostatistical inference. Nevertheless, this is often not a realistic assumption when dealing with nonstationary spatial data and this has motivated a variety of nonstationary spatial modeling workflows such as trend and residual decomposition, cosimulation with secondary features, and spatial segmentation and independent modeling over stationary subdomains. The advent of deep learning technologies has enabled new workflows for modeling spatial relationships. However, there is a paucity of demonstrated best practice and general guidance on mitigation of spatial nonstationarity with deep learning in the geospatial context. We demonstrate the impact of two common types of geostatistical spatial nonstationarity on deep learning model prediction performance and propose the mitigation of such impacts using self-attention (vision transformer) models. We demonstrate the utility of vision transformers for the mitigation of nonstationarity with relative errors as low as 10%, exceeding the performance of alternative deep learning methods such as convolutional neural networks. We establish best practice by demonstrating the ability of self-attention networks for modeling large-scale spatial relationships in the presence of commonly observed geospatial nonstationarity.

58 GEOSCIENCES↗

DOE OSTI2022

Fast and accurate influenza forecasting in the United States with Inferno

Infectious disease forecasting is an emerging field and has the potential to improve public health through anticipatory resource allocation, situational awareness, and mitigation planning. By way of exploring and operationalizing disease forecasting, the U.S. Centers for Disease Control and Prevention (CDC) has hosted FluSight since the 2013/14 flu season, an annual flu forecasting challenge. Since FluSight’s onset, forecasters have developed and improved forecasting models in an effort to provide more timely, reliable, and accurate information about the likely progression of the outbreak. While improving the predictive performance of these forecasting models is often the primary objective, it is also important for a forecasting model to run quickly, facilitating further model development and improvement while providing flexibility when deployed in a real-time setting. In this vein I introduce Inferno, a fast and accurate flu forecasting model inspired by Dante, the top performing model in the 2018/19 FluSight challenge. When pseudoprospectively compared to all models that participated in FluSight 2018/19, Inferno would have placed 2nd in the national and regional challenge as well as the state challenge, behind only Dante. Inferno, however, runs in minutes and is trivially parallelizable, while Dante takes hours to run, representing a significant operational improvement with minimal impact to performance. Forecasting challenges like FluSight should continue to monitor and evaluate how they can be modified and expanded to incentivize the development of forecasting models that benefit public health.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2022

Model Package Report: Composite Analysis Solid Waste Release Model (CASWR Model)

This document describes the implementation of a solid waste form release model in GoldSim for the Hanford Site Composite Analysis (CA) Update. This Composite Analysis Solid Waste Release model (CASWR model) was designed to generate deterministic radionuclide release rates for Hanford’s Central Plateau solid waste disposal sites using single realizations of release model coefficients. Five generalized waste form types are used for the conceptual model of waste release: surplus reactor block, cement, soil-debris, grouted residual waste, and ancillary equipment. The surplus reactor block waste form consists of radionuclide leaching from surplus production graphite reactor core blocks via unspecified processes (White et al., 1984 as cited in PNNL-15965). The cement waste form represents solidified wastes whose permeability is much lower than that of the surrounding soil. The soil-debris waste form type is defined as unconsolidated waste mixed with soil material. Tanks and canyon complexes comprise the grouted residual waste form such that their solid waste will be grouted and capped with a surface barrier at the completion of their cleanup. Finally, the ancillary equipment waste form constitutes contaminant releases from ancillary and auxiliary waste form residues associated with tank farms at closure. Individual sub-models are implemented to numerically represent a respective waste form within the CASWR Model: Surplus Reactor Block Sub-model, Cement Sub-model, SoilDebris Sub-model, Grouted Residual Waste Sub-model, and Ancillary Equipment Submodel. Advection is assumed to be the primary transport process governing the release of radionuclides in the Ancillary Equipment and Soil-Debris Sub-models. Diffusion is assumed to be the primary release process in the Grouted Residual Waste and Cement Sub-models. An unspecified zero-order release process is considered in the Surplus Reactor Block Sub-model due to the lack of information regarding actual processes involved in irradiated graphite leaching. The Surplus Reactor Block, Cement, and SoilDebris Sub-models were compared against analytical solutions (PNNL-11800, Composite Analysis for Low-Level Waste Disposal in the 200 Area Plateau of the Hanford Site). The agreement between results of these analytical solutions and the corresponding waste form models verified their correct implementation in GoldSim. A 1-D modeling abstraction approach for the Grouted Residual Waste and Ancillary Equipment Sub-models was adopted from existing Performance Assessment (PA) models. Despite the simplifications made in these sub-models, they were found to be appropriate representations of the waste forms, similar to what was used in the Waste Management Area C PA model (RPP-ENV-58782, Performance Assessment of Waste Management Area C, Hanford Site, Washington, Rev. 0). A sensitivity analysis was conducted to identify the most influential parameters in each waste form sub-model. Suggestions for considering pH-dependent and redox-dependent release mechanisms are formulated through the development of a conditional constant approach in GoldSim.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Effect of Accurate Modeling of Converter Interfaced Generation (CIG) on a Practical Bulk Power System

With increased penetration of converter interfaced generation (CIG) in power systems, it becomes imperative to adequately model them in both steady-state and dynamic simulations. There are several modeling approaches and simulation platforms for CIG. In addition, it is noted that suitability of modeling approaches and models is also dependent on the system conditions. However, their suitability has yet to be systematically analyzed, and this becomes a challenge in the industry without a systemic method. To tackle these challenges, in this paper, we first develop a methodology for identifying the maximum penetration of aggregated rooftop sources that can be represented as negative loads in both the steady state and dynamic analysis. Then we verify the suitability of positive-sequence performance-based model and detailed electromagnetic transient (EMT) model for utility-scale PV plants under different fault conditions. Based on the results, we perform dynamic simulations on a large practical system with both positive-sequence transient stability(TS) and EMT-TS hybrid simulation. Our study shows that positive-sequence performance based PV plant model is more suitable for study cases with faults relatively far away from the point of connection, and that positive-sequence models and TS simulation approach can produce comparable results with reference to EMT-TS hybrid simulation for normally-cleared single-line-to-ground faults, but the inadequacy of positive-sequence based model and simulation approach becomes obvious for prolonged-tripping of unbalanced faults in proximity of the PV plants.

Huang, Qiuhua↗

DOE OSTI2025

Conflict Detection in Open RAN with Recurrent Neural Networks Using Geometric Manifolds

Allowing third-party applications on Radio Access Network (RAN) Intelligent Controllers (RICs) within the OpenRAN (O-RAN) framework introduces conflicting interactions that are often difficult to detect in advance. These conflicts, occurring between third-party applications in the Near RealTime RIC (Near-RT RIC), known as xApps, can lead to performance degradation and instability in O-RAN if not identified early. Existing conflict detection and mitigation solutions in the literature assume that the conflicts are known beforehand, which is not always accurate due to the complex and often hidden relationships between control parameters and Key Performance Indicators (KPIs). In this paper, we propose a novel Recurrent Neural Network (RNN) to detect both known and unknown conflicts in O-RAN xApps as specified in the O-RAN standards. We model the xApps, control parameters, and KPIs with nodes and edges to create graph structures and use the hidden nonEuclidean geometric properties of the Riemannian manifold to train the RNN model. The performance of this proposed model is validated using evaluation metrics and compared with benchmarks. Results demonstrate that the proposed RNN model, leveraging Riemannian geometric properties, can achieve 100% of the F1-score provided by an optimal solution in just 20 iterations.

5G↗

DOE OSTI2025

Conflict Detection in Open RAN with Recurrent Neural Networks Using Geometric Manifolds

5G↗

DOE OSTI2022

Model development for flash boiling spray and validations with isooctane, hexane, ethanol and their binary mixtures

The objective of this study is to develop a comprehensive model for spray atomization under the consideration of flash boiling, which includes nucleation, bubble growth, and droplet breakup. Although flash boiling has been previously studied, there is still no well-established analytical and numerical model of their mechanisms, or general guidelines on the implementation of this atomization strategy. In this study, detailed formulations are elaborated for the unified spray atomization model under both non-flash and flash boiling conditions. Next, the model was implemented into CFD software and 3-D simulations were performed for model validations. Two fuel blends are tested: ethanol-isooctane blends and hexane-isooctane blends. Simulation results were compared with the experimental measurements conducted using two methodologies, DBI (Diffuse Backlight Imaging) system for spray images captured by a high-speed camera, and a PDA (Phase Doppler Anemometry) system for droplet size measurement during the spray injection process. It was found that the current flash boiling model can accurately capture the spray characteristics under different environment conditions, such as the spray plume interactions and liquid penetration. The averaged Sauter mean diameter (SMD) calculations match the measurement sufficiently well. The spray structures for both single- and multi-component fuels during the injection process are well represented by current model.

42 ENGINEERING↗

DOE OSTI2024

Why Does Atmospheric Radiative Heating Weaken Midlatitude Cyclones?

Abstract Recent work has indicated that atmospheric radiative heating reduces the kinetic energy of large‐scale eddies in the midlatitudes. However, a physical mechanism that connects radiation to the midlatitude eddy kinetic energy is still uncertain. Using a high‐resolution general circulation model we perform an experiment in which the radiative cooling profile at each model time step is overwritten with the climatological mean, computed from a control simulation. This approach separates the mean and transient effects of radiative heating on the extratropical circulation. We find that, when radiative heating is fixed, the globally‐averaged eddy kinetic energy is enhanced by ∼6%. We show that thermal radiation dampens temperature anomalies near the surface and tropopause in low‐pressure systems, destroying eddy available potential energy and eddy kinetic energy. We identify this as a possible mechanism by which atmospheric radiative heating weakens midlatitude cyclones.

Mischell, Eric↗

DOE OSTI2025

Prediction of hydration energies of adsorbates at Pt(111) and liquid water interfaces using machine learning

Aqueous phase heterogeneous catalysis is important to various industrial processes, including biomass conversion, Fischer–Tropsch synthesis, and electrocatalysis. Accurate calculation of solvation thermodynamic properties is essential for modeling the performance of catalysts for these processes. Explicit solvation methods employing multiscale modeling, e.g., involving density functional theory and molecular dynamics have emerged for this purpose. Although accurate, these methods are computationally intensive. This study introduces machine learning (ML) models to predict solvation thermodynamics for adsorbates on a Pt(111) surface, aiming to enhance computational efficiency without compromising accuracy. In particular, ML models are developed using a combination of molecular descriptors and fingerprints and trained on previously published water–adsorbate interaction energies, energies of solvation, and free energies of solvation of adsorbates bound to Pt(111). These models achieve root mean square error values of 0.09 eV for interaction energies, 0.04 eV for energies of solvation, and 0.06 eV for free energies of solvation, demonstrating accuracy within the standard error of multiscale modeling. Feature importance analysis reveals that hydrogen bonding, van der Waals interactions, and solvent density, together with the properties of the adsorbate, are critical factors influencing solvation thermodynamics. Furthermore, these findings suggest that ML models can provide rapid and reliable predictions of solvation properties. This approach not only reduces computational costs but also offers insights into the solvation characteristics of adsorbates at Pt(111)–water interfaces.

Adsorption↗

DOE OSTI2023

Assessment of Carbonated Brine Injection as Low-risk Strategy for Geologic Carbon Storage

The main objective of this early-phase research was to evaluate the techno-economic feasibility and risk associated with combined brine and CO 2 storage in SWD wells using brine dissolution in the North Dakota portion of the Williston Basin. Three simulation studies were conducted to investigate: (1) CO 2 phases at different conditions, (2) wellbore compatibility, and (3) long-term storage fate in reservoir. (1) A simple reservoir model and injection simulations were created using data to represent the BEST (brine extraction and storage test) site, an operational SWD facility located near Watford City, North Dakota. The pressure evolution caused by CO 2 comingled in produced water injectate in a layer cake reservoir was then modeled while tracking aqueous CO 2 throughout the project. The salinity of the injection water, the salinity of the reservoir brine, and the amount of dissolved CO 2 comingled in the injection water were varied. (2) A wellbore corrosion model was performed using the CO 2 concentrations selected based on the reservoir modeling to examine the carbonated produced water impact on wellbore. (3) Reactive transport modeling was conducted with the optimal CO 2 concentration for this injection site to study the rock-fluid interactions and CO 2 fate in the reservoir. Results suggest that CO 2 dissolved in produced water can be injected without appreciably increasing subsurface pressure or leakage risks. Pressure buildup was found to vary with salinity but not with CO 2 mass fraction. Simulation results show that lower CO 2 percent mass fraction leads to a higher amount of CO 2 that can be dissolved at a higher injection salinity. Furthermore, the long-term goal of dissolution trapping in a traditional carbon storage project is accomplished from the start, mitigating risks associated with potential migration of buoyant CO 2 , so long as the reservoir pressure and temperature are used to determine the maximum mass fraction of the dissolved CO 2 .

54 ENVIRONMENTAL SCIENCES↗