Search NASA⌕ Search

SEARCH · Search NASA

Results for “supercritical CO 2”

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 73 records · Page 4

The Impact of Impurities on Alloy Behavior in Supercritical CO 2 at 700 °C

As part of round robin testing, specimens of commercial alloys 316, 120, 625, and 740 were exposed to 20 MPa research grade CO 2 at 700 °C for up to 1500 h. The first set of specimens had higher mass gain likely due to impurities not flushed from the autoclave at startup. After this issue was corrected, an identical set of specimens exhibited lower mass gains for both the Fe- and Ni-based alloys. The differences in reaction products were characterized to understand the effect of impurities under these conditions. As suggested by the mass change, thicker oxides were formed in each case, primarily for the Fe-based alloys. Furthermore, for the Ni-based alloys, the difference was primarily in the transient stage of oxidation with no change in rate constant. Alloy 120 exhibited increased internal attack, and differences in the scale phases formed were characterized to better understand the effect of impurities. Alloy 740 primarily exhibited a greater depth of internal attack in the first run.

36 MATERIALS SCIENCE↗

Measurement of Convective Heat Transfer Coefficients With Supercritical CO 2 Using the Wilson-Plot Technique

This paper describes the measurement of convective heat transfer coefficients and friction factors for sCO 2 flowing in a smooth tube and compares the results with published correlations for validation. The paper also describes the Heat Exchange and Experimental Testing (HEET) rig recently designed and built at the U.S. Department of Energy’s (DoE’s) National Energy Technology Laboratory (NETL) in Morgantown, WV. Here, the Wilson-plot technique used for measuring the heat transfer coefficients is described along with the data reduction process. The Wilson-plot technique was chosen as the basis for the design of NETL’s HEET rig. Advantages of the Wilson-plot technique include the (1) ability to measure high convective heat transfer coefficients accurately, (2) ability to measure average heat transfer coefficient for complicated heat exchange geometries like those produced using additive manufacturing, (3) ability to measure heat transfer coefficients on both sides of a heat exchanger independently, and (4) simplicity of experimental setup. Capabilities of the HEET rig include pressure to 24 MPa (3500 psig), temperature to 538 °C (1000 °F), mass flow rate to 1.5 kg/s (3 lb/s), and Re to 500,000. The rig is designed to operate with pure CO 2 or a mixture of CO 2 and up to 10% N2 by volume to study the impact of gas mixtures typical of direct-fired sCO 2 power cycles on the convective heat transfer and pressure drop. Preliminary tests in the HEET rig were performed with smooth stainless-steel tube and pure CO 2 , and the results were compared with published correlations for Nusselt number (Nu) and friction factor. Over a Reynolds number (Re) range from 58,000 to 228,000, measured Nu was compared to predictions using the Dittus and Boelter equation (Kreith and Bohn, 1993, “Principles of Heat Transfer, West Publishing Company”) within 5% and measured friction factors were compared to predictions using the McAdams correlation (“McAdams, 1954, “Heat Transmission,” 3rd ed., McGraw Hill, New York)” for smooth tube to be within 5%.

42 ENGINEERING↗

Investigation of Gas Turbine Internal Cooling Using Supercritical CO 2 —Effect of Surface Roughness and Channel Aspect Ratio

In this paper, an experimental and numerical investigation of internal cooling channels with rib turbulators is presented with sCO 2 as the working fluid at process conditions (pressure-20.7 MPa and temperature up to 150 °C). The effect of channel aspect ratio up to 2:1 on thermal-hydraulic performance is explored in additively manufactured rectangular channels and square channels, both with and without 60 deg ribs on the top and bottom sides. The Wilson-plot method is employed to experimentally measure channel-averaged Nusselt number over a Reynolds number range up to 370,000. The friction factor is calculated from pressure drop and mass flow rate and additionally, the overall thermal performance factor (TPF) is reported. A companion computational fluid dynamics (CFD) simulation is performed for the rib turbulated cooling configurations reported in the experiments using the Reynolds average Navier–Stokes-based turbulence model. The objective of the numerical study is to gain insight into the local heat transfer augmentation in the ribbed channels as a result of varying the aspect ratio, channel configuration (square versus rectangular), operating conditions (Reynolds number) and the surface roughness, an inherent outcome of the additive manufacturing process. Surface roughness is simulated using sand grain roughness height (K S ) calculated from the experimental data, and a comparison is presented with the corresponding channel configuration with varying surface roughness heights starting from smooth surfaces (K S = 0). Experimental results indicate that the heat transfer augmentation is negligible in the rectangular channels with ribs on the long side compared to the square channel. However, it is enhanced by 60% in comparison to placing ribs on the shorter side. The TPF remains constant at around 1 for the entire range of Reynolds numbers consistent with prior work at the National Energy Technology Laboratory (NETL). Here, the simulation results highlight that increased surface roughness can have a favorable considerable influence on Nusselt number and overall thermal performance enhancement.

42 ENGINEERING↗

Technology Development for Modular, Low-Cost, High-Temperature Recuperators for Supercritical CO 2 Power Cycles (Final Report)

The program’s focus was to evaluate, advance, and demonstrate recuperator design concepts, materials, and fabrication methods that facilitate the commercial availability of compact and low-cost recuperators. Emphasis was placed on recompression closed Brayton cycles having a turbine inlet temperature exceeding 700°C and differential pressures on the order of 200 bar and were scalable to accommodate plant sizes from 10 to 1,000 MWe. The primary objectives of the first budget period, which included conceptual design, preliminary development, and prototype fabrication and testing, were accomplished. A new recuperator design concept, stacked-sheet heat exchanger (SSHX), was realized during the project. The design was taken from paper to hardware and successfully tested. The SSHX prototype recuperator met or exceeded program thermal, hydraulic and performance/cost metrics goals. The first budget period was extended to retire risk in scaling the recuperator from the kWt to MWt thermal capacity scale. The fabrication and design requirements for the budget period one extension were met. Delays associated with the Covid-19 pandemic prevented the fabrication and testing of the industrial scale prototype.

14 SOLAR ENERGY↗

Modeling supercritical CO 2 flow and mineralization in reactive host rocks with PFLOTRAN v7.0

Understanding the flow and reactivity of CO 2 injected into geological reservoirs is important for many subsurface applications including secure geologic carbon storage (GCS), critical mineral extraction, enhanced geothermal systems (EGS), and enhanced oil recovery (EOR). Traditionally, subsurface CO 2 injection for GCS applications has focused on geologic formations with favorable subsurface configurations for CO 2 migration and trapping through non-reactive mechanisms such as structural, solubility, and petrophysical trapping. Recently, CO 2 -reactive rocks such as mafic and ultramafic basalts have been investigated for their potential to react with injected CO 2 in situ to simultaneously dissolve host rock minerals and mineralize CO 2 as carbonates. Engineering rapid CO 2 mineralization in the subsurface is attractive because of the increased density of stored CO 2 , the additional safety factors associated with solidification, and the potential to extract valuable critical minerals. Here we present recent developments in the parallel flow and reactive transport simulator PFLOTRAN to model coupled CO 2 -brine flow and reactive transport for a wide range of injection and production applications involving reactive CO 2 -brine systems. These developments are based on the well established and trusted CO 2 flow capabilities in the STOMP-CO 2 simulator. New capabilities added to PFLOTRAN include new CO 2 -brine equations of state with optional thermal coupling, several new constitutive relationships like capillary pressure smoothing and scanning path hysteresis, a fully implicit well model, and native linkage with PFLOTRAN's well-established reactive transport libraries. A series of benchmarks between PFLOTRAN and STOMP-CO 2 verify the newly developed CO 2 -brine flow capabilities. Demonstrations of coupled CO 2 -brine flow modeling and reactive transport show how CO 2 mineralization can be engineered in reactive host rocks. Finally, an example use case involving copper leaching by CO 2 and critical mineral extraction is presented to showcase the strengths of this new implementation. Several limitations still remain, including limited availability of field data to parameterize models. Future work should constrain the evolution of mineral surface area during mineralization and the temperature and/or pH dependence of geochemical reactions for specific systems of interest.

Critical Minerals↗

High temperature thermochemical energy storage system

A thermochemical energy storage system and method of storing thermal energy are described. The energy storing system described herein comprises a reactor comprising: a) a reactor with a CO 2 sorbent including MgO; and b) a supercritical CO 2 source with supercritical CO 2 and H 2 O, wherein the supercritical CO 2 source is in fluid communication with the reactor and the CO 2 sorbent including MgO to allow flow of the supercritical CO 2 and H 2 O between the supercritical CO 2 source and the reactor, thereby allowing contact of CO 2 with the CO 2 sorbent comprising MgO.

Muto, Andrew Jerome↗

User's Manual for the FE/NETL Onshore CO 2 EOR Cost Model, Version 1

This user's manual describes the conceptual and mathematical basis for the FE/NETL Onshore CO 2 EOR Cost Model (a Fortran program). The model performs a cash flow analysis to estimate the cost of implementing CO 2 EOR using supercritical CO 2 by incorporating oil field performance outputs for a pattern from the FE/NETL CO 2 Prophet Model (available on NETL's website along with its associated user's manuals under the Collection Name: FE/NETL CO 2 Prophet Model) and implementing patterns to develop an oil field for CO 2 EOR. The model calculates capital costs, operation and maintenance costs, and financing costs. The user’s manual also describes how to run the FE/NETL Onshore CO 2 EOR Cost Model, along with the model’s file structure, inputs and outputs. The FE/NETL Onshore CO 2 EOR Cost Model is available on NETL's website under the Collection Name: FE/NETL Onshore CO 2 EOR Cost Model.

54 ENVIRONMENTAL SCIENCES↗

Pressurized in situ X-ray diffraction insights into super/subcritical carbonation reaction pathways of steelmaking slags and constituent silicate minerals

This study explores mineral carbonation of industrial stainless steelmaking slags and relevant synthetic constituent minerals via in situ pressurized X-ray diffraction, to clarify carbonation reaction pathways and efficiency for carbon storage and waste valorization. The primary mineral phases of Argon Oxygen Decarburization (AOD) and Continuous Casting (CC) slags, namely, åkermanite, bredigite, cuspidine, merwinite, and β- and γ-C 2 S, were reacted in a custom-built beryllium-capped XRD reactor filled with either water-saturated (wet) subcritical CO 2 (g), or wet supercritical CO 2(SC) , in a series of carbonation experiments. Formation of calcite and aragonite was observed for most Ca-bearing minerals, transient precipitation of metastable vaterite was observed, while carbonation of AOD and CC slags in CO 2(SC) resulted in hydrated crystalline calcium carbonates and indirect evidence of amorphous carbonate. Overall, results obtained suggest that low-pressure subcritical carbonation routes are attractive for standalone slag carbonation processes, while high-pressure supercritical carbonation routes are amiable to symbiotic integration with CO 2(SC) -generating green technologies.

42 ENGINEERING↗

Thermohydraulic experiments on a supercritical carbon dioxide–air microtube heat exchanger

Heat exchangers are critical components in supercritical CO 2 Brayton cycles and provide necessary heat input, recovery, and dissipation. Tubular heat exchangers with unconventionally small tube sizes (tube diameters less than 5 mm) are promising components for supercritical CO 2 cycles and potentially provide excellent structural stability with wide scope of application. This paper provides essential design and fabrication guidelines for a compact shell-and-tube heat exchanger with microtubes (with an inner diameter equal to 1.75 mm). A heat exchanger test rig is used to evaluate the thermohydraulic performance of this heat exchanger with supercritical CO 2 and air as working fluids. Thermohydraulic data are reported for more than forty sets of experiments with varying Reynolds numbers for shell and tube flows. Critical performance metrics are calculated from the data and compared with predictions from a previously described numerical model. The average deviations between the experimental and model results fall within 10% for all critical metrics. This excellent agreement validates the numerical model for supercritical CO 2 heat exchanger optimization and scale-up.

42 ENGINEERING↗

Hydrodynamic instabilities and heat transfer characteristics in the duct flow of a fluid in the supercritical thermodynamic regime

The behavior of fluids at supercritical thermodynamic conditions is inherently complex due to large variations in thermodynamic and transport properties. Recent numerical and experimental investigations illustrate ongoing interest for these fluids, especially supercritical CO 2 and supercritical water, for a variety of applications. For example, supercritical water reactors (SCWR) operate in this extreme condition of high-pressure and temperature, resulting in highly dynamic flow fields and unexpected heat transfer regimes. The potential heat transfer benefits in this regime are directly associated with the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. This work characterizes the hydrodynamic instabilities that arise for fluids at supercritical thermodynamic conditions when buoyancy forces are significant. Two specific configurations are considered, a natural convection cavity flow, and a mixed convection, heated, horizontal channel flow. Natural convection flow in a cavity is a classical configuration with expected behavior below the critical point. This configuration aids in characterizing the effect of the variable properties in the supercritical thermodynamic regime. Further, limited studies in the existing literature have been conducted for low-Reynolds and intermediate-Rayleigh numbers, mixed-convection channel flows for supercritical water, which is the focus of the channel flow configuration. To investigate the thermally driven hydrodynamic instabilities in this regime, a high-order fully-implicit numerical method is used. Such strong variations in thermophysical properties (in particular, density) are difficult to simulate and an altogether compressible framework is needed. Therefore, the compressible Navier-Stokes equations are solved without any additional assumptions. The fully implicit, high-order in space and time, reconstructed discontinuous Galerkin method as implemented within the multi-physics code called ALE3D (Arbitrary Lagrangian and Eulerian in 2D and 3D), developed at Lawrence Livermore National Laboratory (LLNL), is used. This fully implicit, L-stable method accurately captures the compressible nature of the ow in the limit of very low Mach number. It has been widely accepted that above the critical point, only one phase is observed. However, recent research has indicated the existence of the distinct gas-like and liquid-like regions separated by the Widom line, the locus of the maxima of the specific heat. Along the Widom line, density decreases 6-fold, viscosity drops by a factor of 2, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film boiling. A pseudo-film at the heated wall of the cavity and the horizontal channel is observed where buoyancy forces induce mixing through the specific configurations. Further the local Rayleigh and Richardson numbers provide maps of the flow field and the buoyancy forces driving the microscopic mixing. In the first chapter, I describe a background of supercritical fluid and the various applications. The second chapter focuses on the mathematical model and numerical method used for simulations, where a description of the equation of state for supercritical water is described. The third chapter focuses on the natural convection cavity with a heated bottom wall. In this cavity a gas-like and a liquid-like flow within the supercritical thermodynamic regime are observed. The fourth chapter focuses on a forced convection, horizontal channel, distinguishing between the gas-like, liquid-like, and mixed flow regimes. Mixed convection flow, with the addition of gravitational forces in the horizontal channel show the influence of variable properties on the hydrodynamic development, heat transfer, and rising instabilities. The last chapter of this research focuses on characterizing the unstable hydrodynamics through time-averaging processes and analysis of the movement of energy through the developing plumes.

42 ENGINEERING↗

Two-component dynamics in supercritical $\text {CO}_2$ from inelastic X-ray scattering

Supercritical fluids are characterized by unique thermodynamic properties. One of these properties is the existence of two-component dynamics that is associated with distinct low-frequency and high-frequency vibrational responses of the fluid. However, the origin of this behavior remains unknown. By combining inelastic X-ray scattering and molecular dynamics simulations, we show that this behavior can be connected to density heterogeneities arising from molecular clusters. Analyses of measurements and molecular trajectories suggest that the two-component dynamics emerges due to distinct momentum fluctuations of clustered and unbound molecules. This connection between clusters and two-component dynamics highlights the importance of molecular-structural heterogeneities in supercritical fluids, colloids, and condensed-matter systems.

biophysics↗

Gigaton commercial-scale carbon storage and mineralization potential in stacked Columbia River basalt reservoirs

This work presents a detailed supercritical CO 2 storage resource estimation for the stacked basalt reservoirs in the Grande Ronde Basalt of the Columbia River Basalt Group in eastern Washington and Oregon. The assessment aims to derisk the commercialization potential of geologic carbon storage in basalt by leveraging both structural and mineralization trapping of CO 2 in basalt. The structural closures formed by anticlinal ridges and synclinal valleys in Yakima Fold Belt are excellent physical traps to accommodate injected supercritical CO 2 . Rigorous hydraulic testing, well logs and simulation results from the Wallula Basalt Pilot #1 well showed the occurrence of 17 suitable permeable injection zones (up to 2,496 mD) intercalated with dense seals (~2.6E-10 mD) in the Grand Ronde Basalt. In addition, geochemical studies showed fast reactions between supercritical CO 2 and dissolved basalt minerals to form stable carbonates. In conclusion, our calculation indicates up to 40 gigatons (P90) of mineralization storage resources exist in the Grande Ronde Basalt reservoirs.

58 GEOSCIENCES↗

Compact Counter Flow Fluidized Bed Heat Exchanger (Final Report)

Generation 3 concentrating solar power systems which utilize solid granular media as the heat transfer and storage media hold promise in reducing the levelized cost of electricity relative to generation 2 systems which utilize molten salt. The particles can be heated to temperatures exceeding 800 °C which enables the use of high efficiency supercritical CO 2 Brayton power cycles. A key challenge associate with the use of solid media is its low bulk thermal conductivity which limits the ability to transfer heat from the particles to a working fluid, such as supercritical CO 2 . Tradition shell and tube and shell and plate heat exchangers feature gravity driven particle flow which limits the maximum particle side heat transfer coefficient in these systems to approximately 400 W/m 2 K. A fluidized bed heat exchanger does not rely solely on gravity for particle flow. Instead, enhanced mixing and horizontal conveyance is achieved through fluidization induced through air injection at the bottom of a bed of particles. A tube bundle within the bed conveys the supercritical CO 2 which is heated via the turbulent particle motion at the tube wall.

14 SOLAR ENERGY↗

Creep of MARBN-type 9Cr martensitic steel in gaseous CO 2 environment

Here, with an interest in moving from steam to supercritical CO 2 as the working fluid in advanced energy systems, it is important to study the mechanical response of structural alloys to CO 2 -containing environments. MARBN-type 9Cr martensitic steel was originally developed for application in boilers for supercritical and ultra-supercritical power plants where the steam temperature reaches 650 °C. In this research, a MARBN-type 9Cr martensitic steel designed and manufactured at NETL was creep tested in a gaseous CO 2 environment (i.e., 0.1 MPa at 650 °C), and the results are directly compared to creep tests conducted in air. It was found that environmentally-assisted cracking facilitated by carbide formation beneath the growing oxide accelerated failure in gaseous CO 2 . This work indicates that creep-oxidation interactions in advanced martensitic steels represent an important consideration for materials selection in supercritical CO 2 power cycles. Furthermore, this work confirms that more oxidation-resistant alloys will be required for the highest temperature portions of these systems.

36 MATERIALS SCIENCE↗

Simulated CO 2 storage efficiency factors for saline formations of various lithologies and depositional environments using new experimental relative permeability data

Saline formations are attractive geologic reservoirs for permanent carbon dioxide (CO 2 ) storage. Here, the U.S. Department of Energy's National Energy Technology Laboratory (DOE-NETL) has worked to develop and refine methods and tools for the calculation of CO 2 storage potential in subsurface reservoirs. DOE-NETL's CO 2 -SCREEN provides an online tool for executing these storage methods. CO 2 storage efficiency terms are input parameters in DOE-NETL's methods and equations embedded in the CO 2 -SCREEN, which assesses pore space available for CO 2 storage. In this work, a modeling workflow was initiated to refine two CO 2 storage efficiency terms - volumetric displacement (E V ) and microscopic displacement (E d ). The models are based on new experimental relative permeability data that are specific to homogenous lithology and depositional environments of key subsurface saline formations targeted for CO 2 storage. In future work, heterogenous features will be added to this initial modeling effort to update efficiency factors as described in DOE-NETL's methods and CO 2 -SCREEN tool. E V accounts for the volume utilized in the reservoir under the areal plume, while E d accounts for saturation values in the plume to assess efficiency of CO 2 storage at the pore scale. The results of this work are significant in that prior values were based on a limited geologically non-specific relative permeability data set that were collected prior to 2009. Specifically, we applied numerical simulations using TOUGH3 models to update CO 2 storage efficiency values for supercritical CO 2 injection into brine-saturated reservoirs for three lithologies (clastics, limestone, dolomite) and six depositional environments (Marginal Marine, Strand Plain, Deltaic Complex Fluvial, Aeolian, Shallow Marine, and Reef) that have a high potential for geologic CO 2 storage. Experimental relative permeability data in cores from these environments were utilized in the models with corresponding rock type/sedimentary environment. Results of this study showed that dolomite followed by limestone generated higher ranges of storage efficiency compared to clastics. The updated values provided a tighter efficiency range for clastics, lower P 10 but higher P 90 range for limestone, and higher P 10 and P 90 for dolomite. In general, tighter reservoirs with relatively low permeability and porosity were associated with higher E V and E d , showing efficient reservoir and pore utilization in these scenarios. High reservoir pressure and temperature associated with increasing depth increased the E V , and high CO 2 injection rates resulted in increases in E V and E d , while the impact of permeability anisotropy was minimal after the 30-year injection period.

03 NATURAL GAS↗