Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Effective diffusivity”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Transient fuel performance analysis for the preliminary fuel concept of general atomics fast modular reactor

This study investigates the transient fuel performance of General Atomics Fast Modular Reactor (GA-FMR) during accident scenarios, focusing on the behavior of its innovative fuel system that combines high-assay low enriched uranium dioxide (HALEUO2) fuel with SiGA® ceramic matrix composite silicon carbide cladding. The preliminary fuel design’s response was analyzed during reactivity-initiated accidents (RIA) and loss of coolant accidents (LOCA) using BISON fuel performance analysis code, which included both the diffusion enhanced and BISON-FASTGRASS coupled UO 2 models. The RIA analysis demonstrated that effective reactivity control reduced fuel temperature, though with transient fission gas release resulting in additional tensile stress state on the cladding. LOCA simulations revealed differing predictions between the two models: the BISON UO 2 model showed more transient fission gas release but minimal pellet expansion, while the BISON-FASTGRASS UO 2 model predicted less pronounced fission gas release but more fuel swelling and thermal expansion, potentially leading to pellet-cladding mechanical interaction. Here, these findings highlight critical areas for fuel design optimization and identify knowledge gaps requiring further experimental and computational investigation to advance GA-FMR fuel development.

Lee, Soon K. [Argonne National Laboratory (ANL), A↗
DOE OSTI2024

Readout of Oriented Triplet Excitons in Linear Acenes via Room-Temperature Electrically Detected Magnetic Resonance

In this study, optically generated molecular spin centers offer an attractive platform for room-temperature spintronic and quantum applications. The linear acene family of molecules are especially good candidates due to their efficient generation of highly polarized triplet excitons via singlet fission. However, the direct detection and manipulation of these spin centers in thin films via the electrical means desirable for ultimate microelectronic devices has proven challenging. In particular, highly oriented triplet features have previously been detected in crystalline anthracene but longer acenes reveal only doublet features in Electrically-Detected Magnetic Resonance (EDMR). In this work we present EDMR spectra of highly oriented triplet excitons in pentacene for the first time, using a host-guest style device made of tetracene and pentacene. The guest acts as an energetic trap site, permitting the isolation and detection of molecular triplets at room temperature. Modeling of these results shows that the observed resonance features correspond to triplet sublevel transitions on isolated pentacene guest molecules. Rotation of the applied field confirms the tendency of the linear acenes to self-orient with the longest molecular axis perpendicular to the device substrate. Lastly, we find the disappearance of resonant triplet features in the neat acenes is not primarily due to the effects of exciton delocalization, but a broader mechanism of spin relaxation primarily influenced by exciton diffusivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Contrasting SnTe–NaSbTe 2 and SnTe–NaBiTe 2 Thermoelectric Alloys: High Performance Facilitated by Increased Cation Vacancies and Lattice Softening

Defect chemistry is critical to designing high performance thermoelectric materials. In SnTe, the naturally large density of cation vacancies results in excessive hole doping and frustrates the ability to control the thermoelectric properties. Yet, recent work also associates the vacancies with suppressed sound velocities and low lattice thermal conductivity, underscoring the need to understand the interplay between alloying, vacancies, and the transport properties of SnTe. Here, we report solid solutions of SnTe with NaSbTe 2 and NaBiTe 2 (NaSn m SbTe m +2 and NaSn m BiTe m +2 , respectively) and focus on the impact of the ternary alloys on the cation vacancies and thermoelectric properties. We find introduction of NaSbTe 2 , but not NaBiTe 2 , into SnTe nearly doubles the natural concentration of Sn vacancies. Furthermore, DFT calculations suggest that both NaSbTe 2 and NaBiTe 2 facilitate valence band convergence and simultaneously narrow the band gap. These effects improve the power factors but also make the alloys more prone to detrimental bipolar diffusion. Indeed, the performance of NaSn m BiTe m +2 is limited by strong bipolar transport and only exhibits modest maximum ZTs ≈ 0.85 at 900 K. In NaSn m SbTe m +2 however, the doubled vacancy concentration raises the charge carrier density and suppresses bipolar diffusion, resulting in superior power factors than those of the Bi-containing analogues. Lastly, NaSbTe 2 incorporation lowers the sound velocity of SnTe to give glasslike lattice thermal conductivities. Facilitated by the favorable impacts of band convergence, vacancy-augmented hole concentration, and lattice softening, NaSn m SbTe m +2 reaches high ZT ≈ 1.2 at 800–900 K and a competitive average ZT avg of 0.7 over 300–873 K. The difference in ZT between two chemically similar compounds underscores the importance of intrinsic defects in engineering high-performance thermoelectrics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

MONACO: accurate biological network alignment through optimal neighborhood matching between focal nodes

Motivation: Alignment of protein–protein interaction networks can be used for the unsupervised prediction of functional modules, such as protein complexes and signaling pathways, that are conserved across different species. To date, various algorithms have been proposed for biological network alignment, many of which attempt to incorporate topological similarity between the networks into the alignment process with the goal of constructing accurate and biologically meaningful alignments. Especially, random walk models have been shown to be effective for quantifying the global topological relatedness between nodes that belong to different networks by diffusing node-level similarity along the interaction edges. However, these schemes are not ideal for capturing the local topological similarity between nodes. Results: Here, we propose MONACO, a novel and versatile network alignment algorithm that finds highly accurate pairwise and multiple network alignments through the iterative optimal matching of ‘local’ neighborhoods around focal nodes. Extensive performance assessment based on real networks as well as synthetic networks, for which the ground truth is known, demonstrates that MONACO clearly and consistently outperforms all other state-of-the-art network alignment algorithms that we have tested, in terms of accuracy, coherence and topological quality of the aligned network regions. Furthermore, despite the sharply enhanced alignment accuracy, MONACO remains computationally efficient and it scales well with increasing size and number of networks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

CO2 refixation is higher in leaves of woody species with high mesophyll and stomatal resistances to CO2 diffusion

Abstract The percentage of respiratory and photorespiratory CO2 refixed in leaves (Pr) represents part of the CO2 used in photosynthesis. The importance of Pr as well as differences between species and functional types are still not well investigated. In this study, we examine how Pr differs between six temperate and boreal woody species: Betula pendula, Quercus robur, Larix decidua, Pinus sylvestris, Picea abies and Vaccinium vitis-idaea. The study covers early and late successional species, deciduous broadleaves, deciduous conifers, evergreen conifers and evergreen broadleaves. We investigated whether some species or functional types had higher refixation percentages than others, whether leaf traits could predict higher Pr and whether these traits and their impact on Pr changed during growing seasons. Photosynthesis CO2 response (A/Ci)-curves, measured early, mid and late season, were used to estimate and compare Pr, mesophyll resistance (rm) and stomatal resistance (rs) to CO2 diffusion. Additionally, light images and transmission electron microscope images were used to approximate the fraction of intercellular airspace and cell wall thickness. We found that evergreens, especially late successional species, refixed a significantly higher amount of CO2 than the other species throughout the entire growing season. In addition, rm, rs and leaf mass per area, traits that typically are higher in evergreen species, were also significantly, positively correlated with Pr. We suggest that this is due to higher rm decreasing diffusion of (photo) respiratory CO2 out of the leaf. Cell wall thickness had a positive effect on Pr and rm, while the fraction of intercellular airspace had no effect. Both were significantly different between evergreen conifers and other types. Our findings suggest that species with a higher rm use a greater fraction of mitochondria-derived CO2, especially when stomatal conductance is low. This should be taken into account when modeling the overall CO2 fertilization effect for terrestrial ecosystems dominated by high rm species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Prediction and observation of intermodulation sidebands from anharmonic phonons in NaBr

A quantum Langevin model, like models from optomechanics, was developed for phonons. It predicts intermodulation phonon sidebands (IPSs) in anharmonic crystals. Ab initio calculations of anharmonic phonons in rock-salt NaBr showed these spectral features as many-body effects. Modern inelastic neutron scattering measurements on a crystal of NaBr at 300 K revealed diffuse intensity at high phonon energy from a predicted upper IPS. The transverse optical (TO) part of the new features originates from phonon intermodulation between the transverse acoustic (TA) and TO phonons. The longitudinal optical spectral features originate from three-phonon coupling between the TA modes and the TO lattice modes. The partner lower IPS proves to be an intrinsic localized mode. Interactions with the thermal bath broaden and redistribute the spectral weight of the IPS pair. These sidebands are a probe of the anharmonicity and quantum noise of phonons in NaBr and suggest novel interactions between photons and phonons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Incentivizing Adoption of Software Quality Practices

Although many software teams across the laboratories comply with yearly software quality engineering (SQE) assessments, the practice of introducing quality into each phase of the software lifecycle, or the team processes, may vary substantially. Even with the support of a quality engineer, many teams struggle to adapt and right-size software engineering best practices in quality to fit their context, and these activities aren’t framed in a way that motivates teams to take action. In short, software quality is often a “check the box for compliance” activity instead of a cultural practice that both values software quality and knows how to achieve it. In this report, we present the results of our 6600 VISTA Innovation Tournament project, "Incentivizing and Motivating High Confidence and Research Software Teams to Adopt the Practice of Quality." We present our findings and roadmap for future work based on 1) a rapid review of relevant literature, 2) lessons learned from an internal design thinking workshop, and 3) an external Collegeville 2021 workshop. These activities provided an opportunity for team ideation and community engagement/feedback. Based on our findings, we believe a coordinated effort (e.g. strategic communication campaign) aimed at diffusing the innovation of the practice of quality across Sandia National Laboratories could over time effect meaningful organizational change. As such, our roadmap addresses strategies for motivating and incentivizing individuals ranging from early career to seasoned software developers/scientists.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2022

Thermoacoustic Soret separation

In a fluid mixture in a channel with an axial time-averaged temperature gradient, high-amplitude oscillating flow can greatly increase the axial flux of thermal diffusion (Soret) separation of the components of the mixture. The enhancement occurs when the oscillating lateral temperature gradient greatly exceeds the axial gradient, causing a large oscillating concentration that can be favorably time-phased with the oscillating flow. This process can occur even with a negligible pressure oscillation or with a negligible temperature response to pressure, as is the case in most liquid solutions. The thermal boundary condition imposed by realistic solids on thermoacoustic liquids is imperfect, adding mathematical complications that are absent for typical gases, for which the solid surface is temporally isothermal. Compared with gas mixtures, the high Lewis number in typical liquid solutions reduces the separation flux associated with the time-averaged temperature gradient, but it also reduces the remixing associated with the time-averaged mole-fraction gradient. For large enough channels, the second-law separation efficiency is only slightly reduced from that of steady liquid Soret separation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Computational diagnostics for flame acceleration and transition to detonation in a hydrogen/air mixture

A new computational diagnostic method for pressure-induced compressibility is proposed by projecting its local contribution to the chemical explosive mode (CEM) in the chemical explosive mode analysis (CEMA) framework. The new method is validated for the study of detonation development during the deflagration-to-detonation transition (DDT) process. The flame characteristics are identified through the quantification of individual CEM contributions of chemical reaction, diffusion, and pressure-induced compressibility. Numerical simulations are performed to investigate the DDT processes in a stoichiometric hydrogen-air mixture. A Godunov algorithm, fifth-order in space, and third-order in time are used to solve the fully compressible Navier-Stokes equations on a dynamically adapting mesh. A single-step, calibrated chemical diffusive model (CDM) described by Arrhenius kinetics is used for energy release and conservation between the fuel and the product. The new diagnostic method is first applied to onedimensional (1D) canonical flame configurations followed by two-dimensional (2D) simulations of DDT in an obstructed channel where different detonation initiation scenarios are examined using the new CEMA projection formulation. Detailed examinations of the idealized configuration of detonation initiation through shock focusing mechanism at a flame front are also studied using the new formulation. A comparison of the currently proposed CEMA projection and the original formulation by the authors suggests that including the pressure-induced compressibility is essential for the use of CEMA in DDT process. The results also show that the new formulation of CEMA projection can successively capture the detonation initiation through either a gradient mechanism or a direct initiation mechanism, and therefore can be used as an effective local analytical tool for the computational diagnostics of detonation initiation in a DDT process. It was found that detonation development is characterized by a strong contribution of chemistry role to the CEM which is pivotal to the initiation of detonation. The role of compressibility is found enhanced at the edge of the detonation front where diffusion was found to have minimal effects on detonation development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Comparison of mercury (Hg) bioaccumulation with mono- and mixed Lemna minor and Spirodela polyrhiza cultures

Mercury (Hg) is a prevalent and harmful contaminant that persists in the environment. For phytoremediation, it is important to discover which plants can bioaccumulate meaningful amounts of Hg while also tolerating its toxicity. Additionally, increasing biodiversity could create a more resilient and self-sustaining system for remediation. This study explores whether mixed populations of Lemna minor and Spirodela polyrhiza can better bioaccumulate and tolerate Hg than monocultures. Mono- and mixed cultures of L. minor and S. polyrhiza were grown in mesocosms of 0.5 μg/L or 100 μg/L Hg ( HgCl2) spiked water for 96 h. Change in weight of duckweed was used to assess Hg tolerance. Diffusive gradients in thin-films (DGTs) were used as surrogate monitoring devices for bioavailable levels of Hg. For biomass growth, the mixed culture of the L. minor was greater than the monoculture at the high dose. The L. minor accumulated more Hg in the mixed culture at the low dose while the S. polyrhiza was higher in the mixed at the high dose. Hg speciation in water was modeled using Windermere Humic Aqueous Model 7 (WHAM7) to compare the bioavailable species indicated by the DGTs. Potentially due to the controlled conditions, the WHAM7 output of bioavailable Hg was almost 1:1 to that estimated by the DGTs, indicating good predictive capability of geochemical modeling and passive sampler DGT on metal bioavailability. Altogether, the mixed cultures statistically performed as well as or better than the monocultures when tolerating and bioaccumulating Hg. However, there needs to be further work to see if the significant differences translate into practical differences worth the extra resources to maintain multiple species.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2021

Cs retention and diffusion in C-S-H at different Ca/Si ratio

Cement and concrete are commonly used in the construction of repository sites for radioactive wastes. The correct isolation of those contaminants requires good adsorption and low diffusion rates. Both parameters are highly affected by many factors, such as pH, temperature or composition. The large variability of experimental conditions and formulations makes it extraordinarily difficult to tackle the influence of each of them independently in experimental samples. To this effect, molecular dynamics simulations have been employed in this study to investigate the role of the composition in the capacity to retain Cs and diffusivity of these ions in calcium silicate hydrate (C-S-H) gel pores. The results indicate that the adsorption of Cs ions is worsened at high Ca/Si ratios due to a lower interaction of the cations with the C-S-H surface, while hydrophilicity of the C-S-H nanopore rises, resulting in higher long-range ordering and lower diffusion coefficients of Cs ions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Runaway electron deconfinement in SPARC and DIII-D by a passive 3D coil

The operation of a 3D coil—passively driven by the current quench (CQ) loop voltage—for the deconfinement of runaway electrons (REs) is modeled for disruption scenarios in the SPARC and DIII-D tokamaks. Nonlinear magnetohydrodynamic (MHD) modeling is carried out with the NIMROD code including time-dependent magnetic field boundary conditions to simulate the effect of the coil. Further modeling in some cases uses the ASCOT5 code to calculate advection and diffusion coefficients for REs based on the NIMROD-calculated fields, and the DREAM code to compute the runaway evolution in the presence of these transport coefficients. Compared with similar modeling in Tinguely et al (2021 Nucl. Fusion 61 124003), considerably more conservative assumptions are made with the ASCOT5 results, zeroing low levels of transport, particularly in regions in which closed flux surfaces have reformed. Of three coil geometries considered in SPARC, only the n = 1 coil is found to have sufficient resonant components to suppress the runaway current growth. Without the new conservative transport assumptions, full suppression of the RE current is maintained when the thermal quench MHD is included in the simulation or when the RE current is limited to 250kA, but when transport in closed flux regions is fully suppressed, these scenarios allow RE beams on the order of 1–2 MA to appear. Additional modeling is performed to consider the effects of the close ideal wall. In DIII-D, the CQ is modeled for both limited and diverted equilibrium shapes. In the limited shape, the onset of stochasticity is found to be insensitive to the coil current amplitude and governed largely by the evolution of the safety-factor profile. In both devices, prediction of the q-profile evolution is seen to be critical to predicting the later time effects of the coil.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

An Evaluation of the Global Effects of Tritium Emissions from Nuclear Fusion Power

We report that tritium, like all hydrogen isotopes, is difficult to confine and easily diffuses through most materials. As currently planned, fusion power plants will process and handle large quantities of deuterium and tritium as fuel, and therefore, will become sources of tritium input into the environment. Tritium releases from a worldwide distribution of tritium sources (fusion or fission) will lead to higher tritium levels overall and have global impact. This report investigates the hydrologic partitioning of yearly tritium releases assuming 1 g/yr loss per 500 MW of nuclear power generation for varying scales of power generation. On large scales, it is found that the worldwide average levels of tritium could exceed some public health goals and regulatory guidelines. Tritium concentrations at such levels may also have other impacts on the environment through increased ion-pair generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Effect of iron and silica nanoparticles’ size on in vitro human skin binding and penetration

This in vitro skin study determined absorption, diffusion, and binding rates of four [14C]-labeled nanoparticles (NPs): 12 nm Fe 3 O 4 , 32 nm Fe 3 O 4 @SiO 2 , 33 nm SiO 2 , and 78 nm SiO 2 in each layer of human cadaver skin. In vitro microdialysis device and flow-through skin diffusion system were used to measure the binding affinity to the stratum corneum (SC) and permeability into/through skin layer of the four NPs with different physical–chemical properties, respectively, in short (30 min) and/or long (24 hours) exposures. Results show that NP size is an important factor affecting NP percutaneous absorption. The 12 nm Fe 3 O 4 NPs reached the SC and viable epidermis; 32 nm Fe 3 O 4 @SiO 2 core/shell NPs only reached SC. However, 33 nm and 78 nm silica NPs did not permeate SC. Similar patterns were observed for NP binding affinity to SC and dermatopharmacokinetic analysis using the tape stripping method. The binding affinity determination may be a useful method to efficiently screen skin penetration of NPs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Gas diffusion through variably-water-saturated zeolitic tuff: Implications for transport following a subsurface nuclear event

Noble gas transport through geologic media has important applications in the characterization of underground nuclear explosions (UNEs). Without accurate transport models, it is nearly impossible to distinguish between xenon signatures originating from civilian nuclear facilities and UNEs. Understanding xenon transport time through the earth is a key parameter for interpreting measured xenon isotopic ratios. One of the most challenging aspects of modeling gas transport time is accounting for the effect of variable water saturation of geological media. In this study, we utilize bench-scale laboratory experiments to characterize the diffusion of krypton, xenon, and sulfur hexafluoride (SF6) through intact zeolitic tuff under different saturations. Here, we demonstrate that the water in rock cores with low partial saturation dramatically affects xenon transport time compared to that of krypton and SF6 by blocking sites in zeolitic tuff that preferentially adsorb xenon. This leads to breakthrough trends that are strongly influenced by the degree of the rock saturation. Xenon is especially susceptible to this phenomenon, a finding that is crucial to incorporate in subsurface gas transport models used for nuclear event identification. We also find that the breakthrough of SF6 diverges significantly from that of noble gases within our system. When developing field scale models, it is important to understand how the behavior of xenon deviates from chemical tracers used in the field, such as SF6 (Carrigan et al., 1996). These new insights demonstrate the critical need to consider the interplay between rock saturation and fission product sorption during transport modeling, and the importance of evaluating specific interactions between geomedia and gases of interest, which may differ from geomedia interactions with chemical tracers.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2025

Single and Dual-Solute Transport in Quaternary Ammonium-Functionalized Anion Exchange Membranes: Isolating the Impact of Water Volume Fraction

Membranes are essential components of photoelectrochemical CO 2 reduction cells (PEC-CRC), as they regulate solute transport and impact device efficiency. These cells are a promising approach for converting carbon dioxide into valuable chemicals. However, electrochemical reactions involve not just a single electrolyte; membrane transport typically includes multiple electrolytes. When a mixture of solutes is present, the presence of one solute in the membrane can impact the diffusion and sorption of other solutes. In complex systems where multiple solutes are diffusing simultaneously through a hydrated, dense polymer membrane, the resulting transport behavior is poorly understood. This study highlights the individual and two solute transport behavior of potassium formate, acetate, and bicarbonate ions in neutral membranes and quaternary ammonium-functionalized anion exchange membranes (AEMs) based on phenyl acrylate polymer backbones. By isolating the effect of water content, we emphasize how fixed charge density and solute-membrane interactions influence transport behavior. In neutral membranes, the diffusivity follows the order of potassium formate > potassium acetate > potassium bicarbonate. However, in AEMs, the order shifts due to changes in ion-membrane interactions, with potassium bicarbonate diffusivity surpassing that of potassium acetate. These results underscore the importance of ion exchange membranes and ion-membrane interaction in modulating transport properties, providing valuable insights for optimizing membranes in electrochemical applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Hierarchical Membrane Dynamics in Phase-Separated Model Membranes

Lipid membranes, the primary matrix of cell membranes, exhibit a hierarchy of dynamics from molecular motions to collective undulations, which control spatiotemporal membrane phenomena and regulate various cellular functions. Thus, understanding how various dynamic modes influence membrane behavior is crucial to uncovering the design rules of cell membranes and implementing them in future technologies. Previous membrane dynamics studies using neutron spin-echo spectroscopy (NSE) have shed light on the dependence of collective membrane fluctuations on viscoelastic membrane properties, demonstrating that thickness fluctuation modes (100ns) are dictated by membrane viscosity, which directly governs fast (ps) molecular diffusive motions. Here, we focus on phase-separating lipid membranes, composed of DMPC:DSPC mixtures, to understand the effect of domain formation/growth on hierarchical matrix dynamics using selective lipid deuteration and different neutron spectroscopy techniques. NSE measurements on DMPC:DSPC-d83 (70:30 mol%) show that bending rigidity of the DMPC-rich matrix graduallyincreases with decreasing temperature - as the fully-fluid membrane approaches the upper phase transition (i.e. with the formation of transient DSPC clusters) - and increases further in the phase coexistence region with the growth of DSPC domains. Comparatively, NSE measurements on DMPC-d54:DSPC-d83 vesicles show pronounced matrix thickness fluctuations in the fully-fluid phase and drastic suppression of this fluctuation mode in the gel-fluid coexistence phase, even when the DMPC-rich matrix is still in its fluid phase. This suppression in fluctuations suggests an increase in the matrix viscosity with the growth of gel DSPC domains, which is further investigated by selectively measuring lipid diffusion in the DMPC-rich matrix using neutron backscattering spectroscopy. Comparisons of the experimental results with coarse-grained MD simulationsillustrate the underlying coupling mechanisms. Overall, these observations indicate a clear mechanical and viscoelastic coupling betweenlipid domains and their host matrix, suggesting that domain formation influences membrane properties beyond local effects.

Carrillo, Jan Michael↗
DOE OSTI2021

Effect of CO 2 -brine-rock reactions on pore architecture and permeability in dolostone: Implications for CO 2 storage and EOR

Geologic carbon sequestration (GCS) is considered a feasible technology for storing substantive volumes of greenhouse gases in subsurface geological formations. In the reservoir, far from carbon dioxide (CO 2 ) injection wells or in post-injection scenarios, diffusion dominates over advection. This condition conjoins with spatially distributed geochemical reactions to induce heterogeneous changes in pore architecture, i.e. pore body and throat sizes or surface roughness. These changes can affect CO 2 transport properties and storage capacity. In this work, we investigated mineral dissolution and precipitation in dolomite samples saturated with a CO 2 -saturated brine at 93 °C and 34.5 MPa, aged without flow. Two rock types samples, i.e. intergranular- and vuggy-dominant, were selected to investigate changes in pore size, porosity and permeability under reactive conditions. Mineral dissolution and precipitation were characterized using scanning electron microscopy. Changes in pore size were quantified via time-domain nuclear magnetic resonance (TD-NMR) transverse relaxation time (T 2 ) and diffusion coefficient (D) distributions. We show that mineral dissolution likely occurs in highly permeable pathways. These observations are confirmed through analysis of (T 2 ) and diffusion coefficient (D) distributions. In contrast to results during CO 2 -enriched brine continuous injection, mineral precipitation was observed in micropores. The leftward shift of the T 2 peaks, corresponding to micropores, also evidenced mineral precipitation in lowpermeability zones. However, microscale alterations resulted only in a subtle increase in porosity and permeability. Results in this study shed light on effects of geochemical reactions on alteration of rock properties in diffusion-dominated regions during CO 2 storage.

58 GEOSCIENCES↗