Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “neutral model”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Activation of CO 2 by Actinide Cations (Th + , U + , Pu + , and Am + ) as Studied by Guided Ion Beam and Triple Quadrupole Mass Spectrometry

Reactions of CO 2 with Th+ have been studied using guided ion beam tandem mass spectrometry (GIBMS) and with An + (An + = Th + , U + , Pu + , and Am + ) using triple quadrupole inductively coupled plasma mass spectrometry (QQQ-ICP-MS). Additionally, the reactions ThO + + CO and ThO + + CO 2 were examined using GIBMS. Modeling the kinetic energy dependent GIBMS data allowed determination of bond dissociation energies (BDEs) for D o (Th + -O) and D o (OTh + -O) that are in reasonable agreement with previous GIBMS measurements. The QQQ-ICP-MS reactions were studied at higher pressures where multiple collisions between An + and the neutral CO 2 occur. As a consequence, both AnO + and AnO 2 + products were observed for all An + except Am + , where only AmO + was observed. Here, the relative abundances of the observed monoxides compared to the dioxides are consistent with previous reports of the AnO n + (n = 1, 2) BDEs. Comparison of the periodic trends of the group 4 transition metal, lanthanide (Ln), and actinide atomic cations in reactions with CO 2 (a formally spin-forbidden reaction for most M + ground states), and O 2 (a spin unrestricted reaction) indicate that spin conservation plays a minor role, if any, for the heavier An + metals. Further correlation of Ln + and An + + CO 2 reaction efficiencies with the promotion energy (E p ) to the first electronic state with two valence d-electrons (E p (5d 2 ) for Ln + and E p (6d 2 ) for An + ) indicates that the primary limitation in the activation of CO 2 is the energetic cost to promote from the electronic ground state of the atomic metal ion to a reactive state.

bond activation↗
DOE OSTI2023

Two-fluid theory of composite bosons and fermions and the quantum Hall proximity effect

Here we propose a two-fluid description of fractional quantum Hall systems, in which one component is a condensate of composite bosons and the other a Fermi liquid formed by composite fermions (or simply electrons). We employ the theory to model the interface between a fractional quantum Hall liquid and a (composite) Fermi liquid metal, where we find a penetration of quantum Hall condensate into the metallic region reminiscent of the proximity effect in superconductor-metal interfaces. We also find a novel and physically reasonable set of gapped quasielectron and neutral modes in fractional quantum Hall liquids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Feasibility of using nuclear microreactor process heat for bioconversion and agricultural processes

Introduction There is a global goal to reduce greenhouse gas emissions by 43% by 2023. Nuclear microreactors, a subset of small modular reactors, offer a potential solution due to their compact size, transportability, and carbon-neutral power generation capabilities. Methods This study explores the feasibility of using heat from nuclear microreactors for bioconversion and agricultural processes, including transforming biomass into energy carriers and products such as syngas, bio-oil, and pasteurized milk. Operating requirements for gasification, pyrolysis, hydrothermal carbonization, hydrothermal liquefaction, hydrothermal gasification, ethanol production, anaerobic digestion, and pasteurization were obtained through a literature review. A Brayton cycle model based on the eVinci TM microreactor was developed to assess the feasibility of powering these processes using nuclear microreactor heat. Results and Discussion Exergetic efficiency values for high-temperature processes ranged from 72% to 100%, whereas lower-temperature processes ranged from 2% to 53%. These efficiencies depend on the available source temperature for each microreactor design. There were trade-offs between producing net power and using process heat, particularly for high-temperature processes. Three heat exchanger locations were considered: before the turbine (600 ℃ ), between the turbine and regenerator (370 ℃ ), and after the regenerator (192 ℃ ). High-temperature processes like gasification require temperatures too high for feasibility. Middle temperature processes are better suited to a heat exchanger between the turbine and regenerator, while also operable before the turbine. Lower-temperature processes like pasteurization and anaerobic digestion can use waste heat after the regenerator and do not impact power production. These findings are valuable for optimizing nuclear microreactor heat use and aligning with global climate initiatives.

09 BIOMASS FUELS↗
DOE OSTI2023

Heterogeneity or illusion? Track the carbon Kuznets curve of global residential building operations

Residential buildings, the “last mile” sector in the global decarbonization, have become the most significant uncertain factor hindering carbon neutrality with increasing household energy demand. To track the operational carbon in buildings, this study investigates the carbon Kuznets curve (CKC) and the corresponding decoupling status of residential building operations at four emission scales by using the data of 30 countries from 2000 to 2019. The results show that (1) the CKC model can fit more than half of the samples. Most curves have an inverted U-shape, with 76% of emission per household and 82% of total emissions. (2) In the presence of the CKC, over four-fifths of global residential buildings peak regardless of any emission scale. The analysis denotes that the carbon emissions of developed countries reach their peaks earlier. In the total emissions, the samples’ peaking proportion is 20% and 25% with income per capita < 20,000 United States dollars (USD) and 20,000–40,000 USD, respectively. (3) The Tapio decoupling analysis and the threshold regression effectively verify the robustness and the heterogeneity of CKCs, respectively. Strong decoupling effects of CKCs in most countries are demonstrated at the scales of emission per floor space and the total emissions, and the heterogeneity proves the classic inverted U-shaped relationship between economy and emissions doesn’t exist in all emitters. Overall, this study tracks the historical carbon emission trajectories of residential building operations at a global scale, providing reference for different economies to simulate the dynamic of building carbon emissions along with the economic booming.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Kinetics of low temperature plasma assisted NH 3 /H 2 oxidation in a nanosecond-pulsed discharge

Ammonia (NH 3 ) has been widely recognized as one of the carbon-neutral fuels. However, ammonia combustion suffers low reactivity and high N 2 O/NO x emissions. Here, to overcome these issues, this work reports plasma assisted NH 3 /H 2 oxidation and unveils the kinetics of fuel oxidation and N 2 O/NO x formation by combining time-resolved laser diagnostics with plasma modeling. Firstly, we found that the NH 3 consumption is promoted with a H 2 blending ratio of 0.3, due to enhancements of H and OH formation by plasma assisted H 2 dissociation. Secondly, at a high reduced electric field, when the H 2 blending ratio increases, the NH 3 oxidation is promoted due to both the HO 2 formation and strong NO kinetic enhancement via NO-HO 2 and NO 2 -H pathways. In the meantime, it is shown that the NO mole fraction also increases with H 2 blending ratio, because the NO formation is enhanced via N( 2 D)-O 2 pathways, and the DeNO x chemistry is weakened with less NH 2 production. By contrast, at a lower reduced electric field, when the H 2 blending ratio increases, the decreased N( 2 D) formation does not produce enough NO to replenish the NO formation drop caused by lower NH 3 concentration. Thirdly, the reduced electric field non-monotonically affects fuel consumption and N 2 O/NO x formation by manipulating electron energy deposition pathways. The NH 3 consumption is maximized with an optimal reduced electric field where N 2 * excitation and O 2 dissociation are most efficient. When the reduced electric field deviates from its optimum, the NH 3 consumption decreases due to the discharge energy deposition to either vibrational excitation or dissociation of N 2 . The N 2 O/NO x emissions governed by the NH 3 oxidation follow the above NH 3 consumption trend.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A Post-translational Histidine–Histidine Cross-Link Enhances Enzymatic Oxygen Reduction Activity with Greater pH Adaptability

Cross-linked protein residues exist as enzyme cofactors to enable or enhance catalytic activities. Despite their importance in nature, the chemical identity of the cross-links is limited to certain amino acid combinations, whose function and the formation mechanism remain insufficiently understood due to the difficulty in isolating native enzymes without the cross-links. Herein, we report the formation and characterization of both His-Tyr and His-His cross-links under oxidative enzymatic turnover conditions in L29H/F33Y/F43H Mb, a structural and functional model of heme-copper oxidase (HCO). The connectivity of the cross-link was characterized as Ne ε2 (His29)-Cd δ2 (His43) by mass spectrometry (LC-MS/MS) and nuclear magnetic resonance (NMR). Interestingly, formation of the cross-link significantly enhances the oxygen reduction activity of the enzyme at neutral or basic pH with higher product specificity. X-ray crystallography has identified a novel Tyr-His cross-link through a Tyr-O-His linkage. Our mechanistic studies indicate the involvement of high-valent heme-iron and the neighboring tyrosine in an oxidative self-processing pathway to generate the cross-link. In conclusion, this work serves as a new example while providing insights into the enzyme cross-link formation, allowing the design of artificial biocatalysts containing these novel cross-links with higher activity and pH adaptability.

Liu, Yiwei [Univ. of Texas, Austin, TX (United Sta↗
DOE OSTI2021

Critical adjustment of land mitigation pathways for assessing countries' climate progress

Collective progress towards the Paris Agreement’s temperature targets will be assessed by the Agreement’s Global Stocktake process. Global emission pathways by Integrated Assessment Models (IAMs) are a likely benchmark against which aggregated countries’ greenhouse gas mitigation pledges will be assessed. However, this approach requires comparability, which is currently hampered by a 5 GtCO2yr-1 mismatch between the estimated global land anthropogenic flux of IAMs and countries’ GHG inventories (GHGIs) for 2005-2015. We show that this gap is mostly the result of differences in how the anthropogenic forest sink is defined. Countries define managed forest more broadly than IAMs, and on this larger area typically consider anthropogenic also fluxes due to human-induced environmental change, which IAMs do not include their pathways. Using five IAMs and a dynamic global vegetation model (DGVM), we present and apply a novel method that adjusts IAM results to the countries’ definition, closing the gap with GHGIs at global and regional scale. When expressed in a way compatible with GHGIs, adjusted 1.5°C and well-below 2°C pathways have a cumulative CO2 balance until carbon neutrality that is 110–176 GtCO2 lower than the original IAMs pathways. The use of these country-compatible emission pathways under the Global Stocktake would therefore be essential for an accurate assessment of the collective progress towards the Paris Agreement’s climate goals.

Grassi, Giacomo↗
DOE OSTI2022

Experimental and theoretical study of hole scattering in RF sputtered p-type Cu 2 O thin films

One of the key materials of interest for p-type oxide semiconductor thin film electronics is cuprous oxide (Cu 2 O), due to its relatively high hole mobility. Here, in this work, we use experiments, analytical models, and density functional theory calculations to study the scattering mechanisms that determine Hall mobility in two Cu 2 O samples. First, we examine a polycrystalline Cu 2 O thin film deposited by RF magnetron sputtering, and second, a single-crystalline Cu2O bulk substrate. Temperature-dependent Hall measurements indicate that neutral impurity and grain boundary scattering are dominant for the polycrystalline Cu 2 O thin film, while phonon scattering is dominant for single-crystalline Cu 2 O. Our first-principles calculations show that the room-temperature intrinsic hole mobility of Cu 2 O is 106 cm 2 V –1 s –1 , indicating the great promise of the material for p-type electronic devices. This intrinsic mobility is limited by phonon scattering, with the most dominant scattering modes having phonon energies of 88.4 and 17.1 meV. These results indicate that the key pathways to increase the hole mobility in Cu 2 O thin films are by reducing the impurity concentration and by increasing grain size. Our work thus sets the stage for the future development of high performance Cu 2 O-based p-type thin film transistors.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Neoclassical transport of tungsten ion bundles in total-f neoclassical gyrokinetic simulations of a whole-volume JET-like plasma

Here, the application of a bundling technique to model the diverse charge states of tungsten impurity species in total-f gyrokinetic simulations is demonstrated. The gyrokinetic bundling method strategically groups tungsten ions of similar charge, optimizing computational efficiency. The initial radial configuration of these bundles and their respective charges are derived from a coronal approximation and the quasi-neutrality of the plasma. A low-density JET H-mode like plasma is simulated using the neoclassical version of XGC across the entire plasma volume, spanning from the magnetic axis to the divertor. An accumulation of tungsten is observed at the pedestal top, as a result of low-Z tungsten ions moving inward from the scrape-off-layer into the core region and high-Z tungsten ions moving outward from the core into the pedestal. This organization of the fluxes cannot be captured by a single tungsten-ion simulation. Large up-down poloidal asymmetries of tungsten form in the pedestal and strongly influence the direction of neoclassical fluxes. The temperature screening effect and its correlation with asymmetries are analyzed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Neutrally-buoyant balloons experiment – EMRTC 2025

This memo announces the successful execution of the neutrally-buoyant balloon experiment (EMRTC 2025), conducted at New Mexico Institute of Mining and Technology’s Energetic Materials Research and Testing Center (EMRTC) in Socorro, NM between April 3 and May 9, 2025. The experiment aimed to build on EMRTC 2022 balloon experiment through testing the effects of lateral and vertical displacement of release locations in a different terrain than EMRTC 2022. EMRTC 2025 employed dual GPS trackers on each balloon to provide both detailed near-source motion tracking in addition to the low-frequency long-distance trackers used in EMRTC 2022. A total of 156 neutrally-buoyant balloons were released over the week of May 5-9, which was approximately 2.5 times the number released in the EMRTC 2022 campaign. The dual tracking system will allow characterizing complex three-dimensional wind structures close to terrain features and diurnal wind patterns to validate wind and plume dispersion models for treaty verification.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Probabilistic feasibility assessment of sequestration reliance for climate targets

Countries worldwide are transforming their energy systems to achieve Carbon-Neutrality. Investing in renewable resources-based technologies and implementing Carbon Capture and Storage (CCS) are common strategies to achieve higher sequestration levels. Negative emissions through Bioenergy with CCS are expected to play an essential role in the transition to full decarbonization. On top of that, biomass is a limited resource that depends on environmental factors, which create uncertainties related to the amount that can be sustainably provided to the energy system. Here this study emphasizes the relevance of variability in carbon sequestration for achieving climate targets by 2050. Here this paper proposes a probabilistic approach that integrates the Global Change Analysis Model for Latin America (GCAM-LA) with a chance constraint approach. GCAM-LA is used to assess the impact on the energy sector of different limits of sustainable biomass and carbon budget scenarios. The risk associated with exceeding the sequestration capacity of a given region is modeled via Chance Constraint. Results show that electrification is an appropriate long-term decarbonization strategy. It smoothes the effects of uncertainty in sequestration capacity and responds to end-user demands. For this case study, higher levels of electrification are obtained at likelihood levels >66% for end-use sectors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Power and particle exhaust for the ARC fusion power plant

To successfully show that fusion is an attractive energy source, the ARCTM fusion power plant will need to operate with a robust, integrated power and particle exhaust solution. To maximise ARC’s fusion power output while avoiding excessive erosion of the plasma-facing components, we will need to radiatively dissipate most of the power crossing the last-closed flux surface, injecting radiating impurities such as argon or neon to access divertor detachment. Divertor detachment will need to be integrated with a high-performance core plasma, and with efficient impurity pumping to prevent the accumulation of helium ash in the core. To access and control detachment in high-performance plasmas, we have designed ARC with up–down-symmetric divertors, with secondary X-points in long, tightly baffled outer legs. Using a core-edge modelling workflow, we predict that with this divertor design, ARC will access detachment with modest argon seeding in the divertor (c Ar,div ∼0.9%), which should have minimal impact on the core ( Z eff,core <0.5) for reasonable argon enrichment (c Ar,div /c Ar,core =6). Due to the high predicted divertor neutral pressure (p div ∼20 Pa), we predict that ARC will sufficiently pump helium to limit ash accumulation in the core (c He,core <2%) for a helium enrichment of c He,div /c He,core =0.4. ARC’s divertor design is expected to increase the stability of a detachment front in the outer divertor leg, helping to prevent divertor reattachment during smaller heat-flux transients such as scrape-off-layer filaments associated with the quasi-continuous exhaust regime. However, this buffering will not be sufficient to prevent divertor reattachment during large type-I edge-localised modes (ELMs), and as such these will need to be avoided on ARC. Experiments on SPARC will be used to select an integrated scenario which avoids or mitigates type-I-ELMs while maintaining access to detachment, good core fusion performance and sufficient impurity exhaust. SPARC experiments will also be used to finalise ARC’s divertor design, by studying the impact of magnetic and first-wall geometry on detachment stability, impurity enrichment and neutral baffling under conditions similar to those expected for ARC. In conclusion, our analysis finds that ARC will have a viable power and particle exhaust solution which is compatible with high-power operations, and this solution will be validated in experiments on SPARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Switching Modes for Reduction of Peak Voltage Transients in GaN-Based Three Level ANPC Inverter

The switching pole voltage can transition between P,0, and N states in a three-level active neutral point clamped (ANPC) inverter. State transitions between 0-P and 0-N can be realized with different switching modes. Three switching modes- Short 1 , Short 2 , and Full are studied considering the effect of capacitive current paths. The role of clamping and inner switch on/off conditions is determined using simplified equivalent circuit models. Short 2 with clamping switch off is beneficial for turn-on overvoltage suppression but results in high turn-off overvoltage. Full mode with a parallel current sharing path is effective for turn-off overvoltage suppression but leads to high turn-on overvoltage. Hence, a Modified Full mode is proposed to achieve the simultaneous objective of overvoltage suppression at both turn-on and turn-off transient during high load currents. This benefits 3L ANPC operation at a low power factor with a high load current at fundamental voltage zero crossing points. Here, a 650V GaN-based three-level ANPC inverter prototype is designed and used to evaluate the switching modes and overvoltage suppression strategies through double pulse tests and continuous operation.

42 ENGINEERING↗
DOE OSTI2025

Charged Higgs search in 2HDM

In this paper, we present a comprehensive study of the collider search limits on the charged Higgses in the four types of Two Higgs Double Models (2HDM). In addition to constraints from flavor physics measurements, we include both the LEP charged Higgs search channels, as well as the LHC search results on the light and heavy charged Higgses. We consider both the conventional charged Higgs search channels of H ± → τν, cs, cb, tb, and the latest search results on the exotic decay channels H ± → AW ± /HW ± . We find that H ± → AW ± /HW ± are complementary to the conventional fermionic channels for m H ± < m t . For heavy H ± , H ± → AW ± /HW ± extend the reach of tan β beyond that of H ± → tb in the Type-L 2HDM. We also present the combined reach of all the neutral and charged Higgs searches.

2HDM↗
DOE OSTI2023

Recent DIII-D progress toward validating models of tungsten erosion, re-deposition, and migration for application to next-step fusion devices

Abstract Fundamental mechanisms governing the erosion and prompt re-deposition of tungsten impurities in tokamak divertors are identified and analyzed to inform the lifetime of tungsten plasma-facing components in ITER and other future devices. Various experiments conducted at DIII-D to benchmark predictive models are presented, leveraging the DiMES removable sample exposure probe capability and the Metal Rings Campaign, in which toroidally symmetric rows of tungsten-coated tiles were installed in the DIII-D divertor. In tokamak divertors, the width of the electric sheath is of the order of the main ion Larmor radius, and a vast majority of sputtered tungsten impurities are typically ionized within the sheath. Therefore, W prompt redeposition is mainly governed by the ratio of the characteristic ionization mean-free path of neutral tungsten to the width of the sheath. In-situ monitoring of the prompt redeposition of tungsten impurities in divertors is demonstrated via the use of WII/WI line ratios and the ionizations/photon (S/XB) method in L-mode discharges. Even with this relatively limited set of emission measurements, net erosion measurements were found to be a consistent upper bound to an analytic scaling based on the ratio of the W ionization length, λ iz , and the width of the magnetic sheath rather than the ratio of λ iz and the W + gyro-radius. In the far-scrape-off layer (SOL) of the ITER divertor, however, it is calculated that the measurement of photon emissions associated with the ionization of tungsten impurities up to W 5 + may be required. Finally, W deposition patterns on DiMES collector probes, interpreted via DIVIMP-WallDYN modelling, reveal the key roles of progressive W erosion/re-deposition staps and E × B drifts in regulating long-range high-Z material migration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Vortex splitting in two-dimensional fluids and non-neutral electron plasmas with smooth vorticity profiles

Initially, elliptical, quasi-two-dimensional (2D) fluid vortices can split into multiple pieces if the aspect ratio is sufficiently large due to the growth and saturation of perturbations known as Love modes on the vortex edge. Presented here are experiments and numerical simulations, showing that the aspect ratio threshold for vortex splitting is significantly higher for vortices with realistic, smooth edges than that predicted by a simple “vortex patch” model, where the vorticity is treated as piecewise constant inside a deformable boundary. The experiments are conducted by exploiting the E × B drift dynamics of collisionless, pure electron plasmas in a Penning–Malmberg trap, which closely model 2D vortex dynamics due to an isomorphism between the Drift–Poisson equations describing the plasmas and the Euler equations describing ideal fluids. The simulations use a particle-in-cell method to model the evolution of a set of point vortices. The aspect ratio splitting threshold ranges up to about twice as large as the vortex patch prediction and depends on the edge vorticity gradient. This is thought to be due to spatial Landau damping, which decreases the vortex aspect ratio over time and, thus, stabilizes the Love modes. Near the threshold, asymmetric splitting events are observed in which one of the split products contains much less circulation than the other. These results are relevant to a wide range of quasi-2D fluid systems, including geophysical fluids, astrophysical disks, and drift-wave eddies in tokamak plasmas.

Physics↗
DOE OSTI2025

Modeling fast ion losses due to tearing and internal kink perturbations in MAST-U

Fast ion (FI) loss properties in the presence of tearing mode and internal kink perturbations are numerically investigated for discharges in the MAST-U spherical tokamak, utilizing the MARS-F magnetohydrodynamic stability code and the REORBIT test particle guiding-center orbit-following module. Here, tracing about 100 000 particle markers sampled from the equilibrium distribution of the neutral-beam injection induced FIs, it is found that about 10% out of the total strike the limiting surface (including the divertor surface) in MAST-U discharge 46943, assuming a maximum perturbation of 100 G inside the plasma (corresponding to ~6 G at the Mirnov probe location at the outboard mid-plane). Detailed particle tracing, assuming a uniform initial distribution in the 2D phase space (at given radial locations), reveals that initially counter-current FIs launched near the plasma edge are subject to significant prompt losses, while almost all initially co-current ions remain well confined at the assumed perturbation level. Most lost FIs strike the lower-half of the limiting surface. Finite gyro-radius effects prevent lost ions from striking the top-outer corner of the super-X divertor chamber. A scan of the perturbation level (based on discharge 45163) reveals, not surprisingly, an approximately linear scaling of the particle loss fraction (for counter-current FIs) with respect to the perturbation amplitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Investigation of enhanced CO 2 storage in deep saline aquifers by WAG and brine extraction in the Minnelusa sandstone, Wyoming

Geological CO 2 sequestration in deep saline aquifers has been extendedly investigated to reach the goal of carbon neutral as large amount of CO 2 can be reduced in a short time. Here, this study investigates the feasibility of water-alternating-gas (WAG) injection and brine extraction on enhancing CO 2 storage in a target deep saline aquifer, Minnelusa Sandstone in the Powder River Basin (PRB) of Wyoming. An integrated numerical model is developed to account for four scenarios: (1) CO 2 continuous injection (CI) through one injection well, (2) WAG injection through one injection well, (3) CO 2 CI through one injection well and brine extraction through one producer, (4) WAG injection through one injection well and brine extraction through one producer. Simulation results suggest that WAG injection and brine extraction corporately make it possible to enhance CO 2 injectivity while securing CO 2 storage safety. For instance, WAG injection considerably reduces structural trapping contribution while enhances dissolved and residual trapping contributions. As for brine extraction, it can decrease the maximum averaged reservoir pressure by 37% and 22% for CI and WAG injection, respectively. Besides, sensitive analyses of the operational parameters for the fourth scenarios are performed. Results reveal that, when the amount of total CO 2 injection has been predetermined, CO 2 injection time and rate within each period of WAG should be small whereas it is preferable to have large water injection time and rate. This is to secure desirable dissolved and residual trapping contributions (store CO 2 with safety). Sensitive degree results confirm that CO 2 injection rate, water injection time, and CO 2 injection time of WAG injection exhibit significant effects on CO 2 storage whereas the impacts of producer bottom-hole pressure (BHP), water injection rate, and CO 2 -water injection time are relatively weaker. This study not only sheds light on achieving double-win goal: enhance CO 2 injectivity and store CO 2 with safety, but also can be a critical reference for other CO 2 storage in deep saline aquifers.

58 GEOSCIENCES↗