Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RSM”

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.

17 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

DBD plasma-assisted ethanol steam reforming for green H 2 production: Process optimization through response surface methodology (RSM)

Herein this work investigates ethanol steam reforming (ESR) to produce hydrogen (H 2 ) in a dielectric barrier discharge (DBD) plasma reactor. A five-level, three-factor experiment design was performed using a response surface methodology (RSM) to evaluate the combined effects of the three process parameters, including discharge power, total flow rate, and ethanol-to-water (EtOH/H 2 O) molar ratio on the plasma-assisted ESR reaction. Quadratic regression models were employed in RSM to fit the experimental results and present the correlation between process parameters and targeted responses (EtOH conversion, H 2 yield, H 2 selectivity, and specific energy requirement (SER) for H 2 production). The results suggested that the EtOH/H 2 O molar ratio is considered to have the most significant effect on the EtOH conversion and H 2 , H 2 selectivity, while the total flow rate is the most significant parameter determining SER for H 2 production. Process optimization demonstrated the optimal process conditions, including a discharge power of 55.9 W, a total flow rate of 26.7 ml/min, and an EtOH/H 2 O molar ratio equal to 0.34. A validation test was performed and confirmed the feasibility of the optimization process.

08 HYDROGEN↗
DOE OSTI2020

Vitrification of High-Cr Glass in Research-Scale Melter

This test report describes the experimental results from a small-scale test using the research-scale melter (RSM) at Pacific Northwest National Laboratory to demonstrate processing of a high-Cr simulated feed stream, designated HLW-HCr-16. The RSM is a small, joule-heated melter capable of processing melter feed continuously. The melter is equipped with Inconel ® 693 electrodes, Monofrax ® K-3 refractory, and an Inconel 690 pour spout. An electric kiln surrounds the melter body and minimizes heat loss from the melter body during operation. The RSM is equipped with an offgas treatment system that employs quenching, wet scrubbing, and high-efficiency mist elimination. The glass-discharge section is heated to facilitate pouring of the glass. The melter is fitted with a melt cavity that is ~25 cm (10 in.) in diameter with a nominal glass depth of 8.9 cm (3.5 in.). The melter was operated with a target glass temperature of 1150°C and target plenum temperature between 550°C and 700°C for this test. The air bubbling rate was 4.2 L/min. Overall, during the continuous operation of the melter for ~ 103 hours, ~ 141 kg of glass was produced. At the conclusion of the test, the melter and exhaust lines were visually inspected for particulate deposition and corrosion. Entrained material had adhered to the underside of the melter lid and to the exhaust piping. Enrichments in elements such as Cl, F, B, K, Li, P, Na, and S were measured in these deposits through inductively coupled plasma–optical emission spectroscopy and X-ray fluorescence analysis. When the melter electrodes and air bubbler tube were removed from the glass in the RSM, the electrodes appeared discolored, but no significant loss of metal was observed. The processing of a high-Cr simulant, HLW-HCr-16, in the RSM produced glass at an average rate of 1.36 kg/h, equaling a melter-surface-area normalized glass generation rate of 654 kg/day/m 2 . The resulting glass met the toxicity characteristic leaching procedure requirement. Test results of crystallinity, electrical conductivity, and viscosity showed good processing properties of this high-Cr high-level waste glass. RSM offgas was also sampled and analyzed at periodic intervals during steady-state operating conditions. The total decontamination factor averaged by four sampling periods was 134. The concentrations of CO and NO x in emissions were 237 to 422 and 69 to 94 parts per million by volume, respectively.

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

Full-potential KKR within the removed-sphere method: A practical and accurate solution to the Poisson equation

Here an efficient and accurate generalization of the removed-sphere method (RSM) to solve the Poisson equation for total charge density in a solid with space-filling convex Voronoi polyhedra (VPs) and any symmetry is presented. The generalized RSM avoids the use of multipoles and VP shape functions for cellular integrals, which have associated ill-convergent large, double-internal L sums in spherical-harmonic expansions, so that fast convergence in single-L sums is reached. Our RSM adopts full Ewald formulation to work for all configurations or when symmetry breaking occurs, such as for atomic displacements or elastic constant calculations. The structure-dependent coefficients A L that define RSM can be calculated once for a fixed structure and speed up the whole self-consistent-field procedure. The accuracy and rapid convergence properties are confirmed using two analytic models, including the Coulomb potential and energy. We then implement the full-potential RSM using the Green's function Korringa-Kohn-Rostoker (KKR) method for real applications and compare the results with other first-principle methods and experimental data, showing that they are equally as accurate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Testing and Design of Discriminants for Local Seismic Events Recorded during the Redmond Salt Mine Monitoring Experiment

The Redmond Salt Mine (RSM) Monitoring Experiment in Utah was designed to record seismoacoustic data at distances less than 50 km for algorithm testing and development. During the experiment from October 2017 to July 2019, six broadband seismic stations were operating at a time, with three of them having fixed locations for the duration, whereas the three other stations were moved to different locations every one-and-half to two-and-half months. RSM operations consist of nighttime underground blasting several times per week. The RSM is located in proximity to a belt of active seismicity, allowing direct comparison of natural and anthropogenic sources. Using the recorded data set, we built 1373 events with local magnitude (M L ) of -2.4 and lower to 3.3. For 75 blasts (RMEs) from the Redmond Salt Mine and 206 tectonic earthquakes (EQs), both M L and the coda duration magnitude (MC) are well constrained. We used these events to test and design discriminants that separate the RMEs from the EQs and are effective at local distances. The discriminants consist of M L -M C , low-frequency Sg to high-frequency Sg, Pg/Sg phase-amplitude ratios, and Rg/Sg spectral amplitude ratios, as well as different combinations of two or more of these classifiers. The areas under the receiver operating characteristic curves (AUCs) of 0.92–1.0 for M L -M C , low-frequency Sg to high-frequency Sg, and Rg/Sg indicate that these discriminants are very effective. Conversely, the AUC of only 0.57 for Pg/Sg suggests that this discriminant is only slightly better than a random classifier. Among the effective classifiers, Rg/Sg, shows the lowest likelihood of misclassification (4.3%) for the populations. In conclusion, results of joint discriminant analyses suggest that even the arguably ineffective single classifier, like Pg/Sg in this case, can provide some value when used in combination with others.

58 GEOSCIENCES↗
DOE OSTI2024

Analysis of strain in ion implanted 4H-SiC by fringes observed in synchrotron X-ray topography

A novel high energy implantation system has been successfully developed to fabricate 4H-SiC superjunction devices for medium and high voltage via implantation of dopant atoms with multi-energies ranging from 13 to 66 MeV. The significantly higher levels of energy used compared to conventional implantation processes, necessitates detailed characterization of the lattice damage caused by implantation. To achieve this by employing the novel high energy system, 4H-SiC wafer with 12 μm epilayers were blanket implanted by 13.8–65.7 MeV Al atoms. The lattice damages induced by the implantation were primarily characterized by Synchrotron X-ray Plane Wave Topography (SXPWT) and Reciprocal Space Mapping (RSM). Topographs reveal fringe contrast akin to multiple asymmetric diffraction peaks with an angular separation of only 2″ (arcseconds) observed on rocking curves, indicating inhomogeneous strain distribution across the implanted layer. The strain profile of the implanted layer was extracted from the fringe contrast by applying Rocking-curve Analysis by Dynamical Simulation (RADS). In conclusion, the maximum strain value is similar to that measured on the RSM.

A1. Characterization↗
DOE OSTI2023

Sustainable bioleaching of lithium-ion batteries for critical metal recovery: Process optimization through design of experiments and thermodynamic modeling

Recycling spent lithium-ion batteries (LIBs) could alleviate supply risks for critical metals and be less harmful to the environment compared to new production of metals from mining. Developing a cost-effective LIB bioleaching process could be a promising alternative to traditional energy-intensive recycling technologies. Here, this study aimed to optimize bioleaching conditions for maximum economic competitiveness through design of experiments using iterative response surface methodology (RSM), assisted by thermodynamic modeling. The optimal condition was identified as 2.5% pulp density in 75 mM gluconic acid biolixiviant at 55°C for 30 h which could recover 57%–84% of nickel, 71%–86% of cobalt, and 100% of lithium and manganese, yielding a 17%–26% net profit margin. The recommended pulp density and acid concentrations, together with the observed metal solubilization, were supported by thermodynamic modeling predictions. Our study demonstrated that combining RSM with thermodynamic simulations could be a powerful tool for optimizing bioleaching conditions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Strain Relaxation and Relative Defect Density with Thickness in MBE-Grown Ge 0.85 Sn 0.15 on Ge(001)

Germanium–tin (GeSn) alloys are emerging as promising materials for mid-infrared optoelectronics and silicon-compatible photonic devices, owing to their tunable direct bandgap. However, the growth of high-quality GeSn films with high Sn content remains challenging due to strain-induced defect formation. In this study, we investigate the role of film thickness on strain-induced relaxation, defect density, and Sn segregation. A series of five samples with varying thicknesses and ∼15% Sn-containing GeSn layers were grown, ranging from the critical thickness for strain relaxation to the onset of Sn segregation. All GeSn samples were analyzed using X-ray diffraction reciprocal space mapping (XRD-RSM) to explore the evolution of strain-induced relaxation as a function of thickness. Photoluminescence measurements reveal that increasing the GeSn thickness enhances strain relaxation while reducing defect-related emission, indicating a decrease in effective defect density prior to reaching the threshold thickness of GeSn layer. At a thickness of ∼150 nm, the GeSn layer shows the onset of Sn segregation, evident in the XRD-RSM spectrum, marking the threshold thickness for Sn segregation. This work defines an effective growth window in terms of thickness (35 to 150 nm) for fabricating relaxed, defect-suppressed GeSn layers with 15% Sn content. These findings emphasize the crucial role of thickness control in balancing strain relaxation and defect suppression, advancing the fabrication of high-quality, high Sn-content relaxed GeSn using molecular beam epitaxy.

Defects↗
DOE OSTI2022

Study of carbon nanotube embedded honey as a resistive switching material

In this paper, natural organic honey embedded with carbon nanotubes (CNTs) was studied as a resistive switching material for biodegradable nonvolatile memory in emerging neuromorphic systems. CNTs were dispersed in a honey-water solution with the concentration of 0.2 wt% CNT and 30 wt% honey. The final honey-CNT-water mixture was spin-coated and dried into a thin film sandwiched in between Cu bottom electrode and Al top electrode to form a honey-CNT based resistive switching memory (RSM). Surface morphology, electrical characteristics and current conduction mechanism were investigated. The results show that although CNTs formed agglomerations in the dried honey-CNT film, both switching speed and the stability in SET and RESET process of honey-CNT RSM were improved. The mechanism of current conduction in CNT is governed by Ohm’s law in low-resistance state and the low-voltage range in high-resistance state, but transits to the space charge limited conduction at high voltages approaching the SET voltage.

36 MATERIALS SCIENCE↗
DOE OSTI2020

CONTROLLED END GAS AUTO IGNITION WITH EXHAUST GAS RECIRCULATION ON A STOICHIOMETRIC, SPARK IGNITED, NATURAL GAS ENGINE

The goal of this study is to address fundamental limitations to achieving diesel-like efficiencies in heavy duty on-highway natural gas (NG) engines. Engine knock and misfire are barriers to pathways leading to higher efficiency engines. This study explores enabling technologies for development of high efficiency stoichiometric, spark ignited, natural gas engines. These include design strategies for fast and stable combustion and higher dilution tolerance. Additionally, advanced control methodologies are implemented to maintain stable operation between knock and misfire limits. To implement controlled end- gas autoignition (C-EGAI) strategies a Combustion Intensity Metric (CIM) is used for ignition control with the use of a Woodward large engine control module (LECM). Tests were conducted using a single cylinder, variable compression ratio, cooperative fuel research (CFR) engine with baseline conditions of 900 RPM, engine load of 800 kPa indicated mean effective pressure (IMEP), and stoichiometric air/fuel ratio. Exhaust gas recirculation (EGR) tests were performed using a custom EGR system that simulates a high pressure EGR loop and can provide a range of EGR rates from 0 to 40%. The experimental measurements included the variance of EGR rate, compression ratio, engine speed, IMEP, and CIM. These five variables were optimized through a Modified Box- Benken design Surface Response Method (RSM), with brake efficiency as the merit function. A positive linear correlation between CIM and f-EGAI was identified. Consequently, CIM was used as the feedback control parameter for C-EGAI. As such, implementation of C-EGAI effectively allowed for the utilization of high EGR rates and CRs, controlling combustion between a narrower gap between knock and lean limits. The change from fixed to parametric ignition timing with CIM targeted select values of f-EGAI with an average coefficient of variance (COV) of peak pressure of 5.4. The RSM efficiency optimization concluded with operational conditions of 1080 RPM, 1150 kPa IMEP, 10.55:1 compression ratio, and 17.8% EGR rate with a brake efficiency of 21.3%. At this optimized point of peak performance, a f-EGAI for C-EGAI was observed at 34.1% heat release due to auto ignition, a knock onset crank angle value of 10.3° aTDC and ignition timing of - 24.7° aTDC. This work has demonstrated that combustion at a fixed f-EGAI can be maintained through advanced ignition control of CIM without experiencing heavy knocking events.

Cooperative Fuel Research Engine, Auto Ignition, E↗
DOE OSTI2022

Effect of Annealing Conditions on Recovery of Lattice Damage in a High-Energy-Implanted 4H-SiC Superjunction PIN Diode

A high energy ion implantation system has been recently developed at the Tandem Van de Graaff facility at Brookhaven National Laboratory with tunable energy to 150 MeV capable of multi-step, deep implantation in 4H-SiC wafers with dopant atoms, such as B, P, Al, and N. Medium and high voltage devices with deep junctions can be fabricated using this system. Lattice strain introduced by the implantation process needs to be recovered and dopant atoms activated by appropriate annealing process as the device performance is strongly associated with the extent of recovery of the lattice. Using Synchrotron X-ray Rocking Curve Topography (SXRCT) and Reciprocal Space Mapping (RSM), the strain induced by high energy implantation of Al and N in 4H-SiC in different patterns are measured and mapped. It is observed that the strain levels correlate with the total fluence levels. PiN diodes fabricated on these implanted wafers were then annealed at temperatures ranging from 1700 °C to 2000 °C for 60 min. The SXRCT and RSM analysis of the annealed samples suggests that for the same annealing duration, higher temperature leads to better lattice recovery.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Numerical Optimization of Spray-Guided Spark Assistance for Cold Idle Operation in a Heavy-Duty Gasoline Compression Ignition Engine

This article describes the results of a response surface model (RSM)-based numerical optimization campaign for spray-guided spark assistance at cold operations in a heavy-duty gasoline compression ignition (GCI) engine. On the basis of an earlier work on spark-assisted GCI cold combustion, a space-filling design of experiments (DoE) method was first undertaken to investigate a multitude of hardware design variables and engine operating parameters. The main design variables included the number of injector nozzles, fuel split quantities and injection timings, and spark timing. The objective variables were engine combustion efficiency (ŋ c ), maximum pressure rise rate (MPRR), and engine-out nitrogen oxide (NO x ) emissions. A total of 150 design candidates were automatically generated using the Sobol sequence method provided by the commercial software package, CAESES. Then, closed-cycle computational fluid dynamic (CFD) spark-assisted GCI simulations under cold idling operations were performed. The outcomes from the CFD-DoE design campaign were utilized to construct high-fidelity RSMs that allowed for further design optimization of the spark plug- and fuel injector-related design variables, along with fuel injection strategy parameters. A merit function with respect to objective variables was formulated with an appropriate weight assignment on each objective variable. Finally, the best design candidate was identified from the RSM-based optimization process and further validated in the CFD analysis. The best design candidate showed the potential to significantly improve combustion efficiency (ŋ c > 90%) over the baseline at cold idle while satisfying MPRR and NOx emissions constraints (MPRR < 5 bar/CAD and NO x < 4.5 g/kWh).

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Multi-objective optimization of peel and shear strengths in ultrasonic metal welding using machine learning-based response surface methodology

Ultrasonic metal welding (UMW) is a solid-state joining technique with varied industrial applications. Despite of its numerous advantages, UMW has a relative narrow operating window and is sensitive to variations in process conditions. As such, it is imperative to quantitatively characterize the influence of welding parameters on the resulting joint quality. The quantification model can be subsequently used to optimize the parameters. Conventional response surface methodology (RSM) usually employs linear or polynomial models, which may not be able to capture the intricate, nonlinear input-output relationships in UMW. Furthermore, some UMW applications call for simultaneous optimization of multiple quality indices such as peel strength, shear strength, electrical conductivity, and thermal conductivity. To address these challenges, this paper develops a machine learning (ML)-based RSM to model the input-output relationships in UMW and jointly optimize two quality indices, namely, peel and shear strengths. The performance of various ML methods including spline regression, Gaussian process regression (GPR), support vector regression (SVR), and conventional polynomial regression models with different orders is compared. A case study using experimental data shows that GPR with radial basis function (RBF) kernel and SVR with RBF kernel achieve the best prediction accuracy. The obtained response surface models are then used to optimize a compound joint strength indicator that is defined as the average of normalized shear and peel strengths. In addition, the case study reveals different patterns in the response surfaces of shear and peel strengths, which has not been systematically studied in the literature. While developed for the UMW application, the method can be extended to other manufacturing processes.

42 ENGINEERING↗
DOE OSTI2025

High-fidelity forced convection simulations of the University of Wisconsin–Madison air-cooled reactor cavity cooling system

Among the various types of advanced nuclear technologies that exist, high-temperature gas-cooled reactors (HTGRs) stand out for their co-generation capabilities and exceptional passive safety systems. HTGRs rely on a reactor cavity cooling system (RCCS) to remove decay heat from the reactor pressure vessel (RPV) during the course of long-term transients. Several studies have focused on experimental and numerical modeling of the RCCS; however, the literature reveals a notable gap when it comes to high-fidelity simulations such as direct numerical simulations (DNS) or large eddy simulation (LES) models of the RCCS. An LES model was developed for forced convection inside the RCCS so as to obtain a high-fidelity representation of the RCCS at a scale never before attempted. Furthermore, the LES model was compared against three different Reynolds-averaged Navier–Stokes (RANS) models (Realizable k − ϵ, k − ω SST, and Reynolds stress model [RSM]). This comparison verified the robustness and accuracy of the turbulence modeling approach of the RANS models for the RCCS. Finally, based on the comparison between the LES and the RANS models, the k − ω SST RANS model of the experimental facility was developed and compared against the available experimental data.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Engineering acid-base active sites on an N-doped biochar catalyst for microwave-assisted biodiesel production from waste cooking oil

Acid-base bifunctional catalysts play a crucial role in converting high-free fatty acid (FFA) feedstocks into biodiesel by enabling the simultaneous esterification and transesterification of triglycerides in a single step, thereby minimizing saponification and streamlining the process. In this study, we report the design and development of a novel N-doped acid-base bifunctional biochar catalyst synthesized from renewable biomass precursors- glucose and chitosan. The catalyst’s acidic sites were introduced via surface functionalization with o-sulfobenzoic acid cyclic anhydride, grafting -SO₃H groups onto the carbon framework, while basic sites were generated from amine functionalities (-NH) derived from chitosan, yielding a robust and truly bifunctional surface. The catalyst was characterized and used for the microwave-assisted transesterification of waste cooking oil, and process optimization using response surface methodology (RSM-CCD) yielding a maximum biodiesel yield of 98.8 ± 0.4% under optimal conditions (methanol-to-oil molar ratio 20:1, catalyst loading 5.7 wt%, temperature 99 °C, and reaction time 46.6 min). The catalyst demonstrated good stability, maintaining a yield of 92.2 ± 0.4% after eight reuse cycles. A Life Cycle Assessment (LCA) conducted for 1 kg of biodiesel production revealed a global warming potential of 0.85 kg CO 2 eq and human toxicity potential of 1.59 kg 1,4-DB eq, with the largest environmental contributions arising from catalyst preparation and reuse. In conclusion, this work highlights a green and circular approach to biodiesel production, combining renewable catalyst design, waste feedstock utilization, and quantitative sustainability assessment to guide future developments in sustainable catalysis and biofuel engineering.

Bifunctional catalyst↗
DOE OSTI2023

Mixing in Low Reynolds Number Reacting Impinging Jets in Crossflow

Previous efforts to model uranyl fluoride formation in an impinging jet gas reactor underpredicted spatial mixing and overpredicted chemical conversion into particulates. The previous fluid dynamics model was based on the solution of the Reynolds Averaged Navier Stokes equations. After simulating fluid dynamics, aerosol dynamics were superimposed onto CFD-simulated gas reactant species concentrations. The current work explores the influence of complex unsteady flow features on the overall flow physics and chemistry for a low Reynolds number, opposed flow, impinging jet gas reactor where there is a low Reynolds number cross flow. The objective of this study was to assess the impact of model formulation on scalar mixing and transport. Here, transient flow simulations were performed using Scale Resolving Simulations. Large-Eddy Simulations with the dynamic Smagorinsky turbulence model were performed along with simulations which directly resolved the flow. Average and root-mean-square (RMS) velocities and species concentrations were computed along with modeled and resolved turbulence kinetic energy (TKE), modeled turbulence dissipation, and modeled turbulent viscosity. Lagrangian flow tracers were also used to quantify species concentrations along path lines emanating from the jet tips. Transient simulation data were compared to results from RANS simulations using the k-ω shear stress transport (SST) model and Reynolds Stress Model (RSM). Transient simulations showed spatial mixing patterns which were more consistent with experimental data and helped elucidate the process of particle formation observed in experiments.

42 ENGINEERING↗
DOE OSTI2020

Chemical Impacts of Potential CO2 and Brine Leakage on Groundwater Quality with Quantitative Risk Assessment: A Case Study of the Farnsworth Unit

Potential leakage of reservoir fluids is considered a key risk factor for geologic CO2 sequestration (GCS), with concerns of their chemical impacts on the quality of overlying underground sources of drinking water (USDWs). Effective risk assessment provides useful information to guide GCS activities for protecting USDWs. In this study, we present a quantified risk assessment case study of an active commercial-scale CO2-enhanced oil recovery (CO2-EOR) and sequestration field, the Farnsworth Unit (FWU). Specific objectives of this study include: (1) to quantify potential risks of CO2 and brine leakage to the overlying USDW quality with response surface methodology (RSM); and (2) to identify water chemistry indicators for early detection criteria. Results suggest that trace metals (e.g., arsenic and selenium) are less likely to become a risk due to their adsorption onto clay minerals; no-impact thresholds based on site monitoring data could be a preferable reference for early groundwater quality evaluation; and pH is suggested as an indicator for early detection of a leakage. This study may provide quantitative insight for monitoring strategies on GCS sites to enhance the safety of long-term CO2 sequestration.

58 GEOSCIENCES↗
DOE OSTI2021

Risk Assessment and Management Workflow—An Example of the Southwest Regional Partnership

This paper summarizes the risk assessment and management workflow developed and applied to the Southwest Regional Partnership on Carbon Sequestration (SWP) Phase III Demonstration Project. The risk assessment and management workflow consists of six primary tasks, including management planning, identification, qualitative analysis, quantitative analysis, response planning, and monitoring. Within the workflow, the SWP assembled and iteratively updated a risk registry that identifies risks for all major activities of the project. Risk elements were ranked with respect to the potential impact to the project and the likelihood of occurrence. Both qualitative and quantitative risk analyses were performed. To graphically depict the interactions among risk elements and help building risk scenarios, process influence diagrams were used to represent the interactions. The SWP employed quantitative methods of risk analysis including Response Surface Method (RSM), Polynomial Chaos Expansion (PCE), and the National Risk Assessment Partnership (NRAP) toolset. The SWP also developed risk response planning and performed risk control and monitoring to prevent the risks from affecting the project and ensure the effectiveness of risk management. As part of risk control and monitoring, existing and new risks have been tracked and the response plan was subsequently evaluated. Findings and lessons learned from the SWP’s risk assessment and management efforts will provide valuable information for other commercial geological CO2 storage projects.

Lee, Si-Yong↗