SEARCH · Engineering Papers

Results for “dependence”

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

SourceResearch connectionFrom yearThrough year

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

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Band structure dependent electronic localization in macroscopic films of single-chirality single-wall carbon nanotubes

Significant understanding has been achieved over the last few decades regarding chirality-dependent properties of single-wall carbon nanotubes (SWCNTs), primarily through single-tube studies. However, macroscopic manifestations of chirality dependence have been limited, especially in electronic transport, despite the fact that such distinct behaviors are needed for many applications of SWCNT-based devices. In addition, developing reliable transport theory is challenging since a description of localization phenomena in an assembly of nanoobjects requires precise knowledge of disorder on multiple spatial scales, particularly if the ensemble is heterogeneous. Here, we report an observation of pronounced chirality-dependent electronic localization in temperature and magnetic field dependent conductivity measurements on macroscopic films of single-chirality SWCNTs. The samples included large-gap semiconducting (6,5) and (10,3) films, narrow-gap semiconducting (7,4) and (8,5) films, and armchair metallic (6,6) films. Experimental data and theoretical calculations revealed Mott variable-range-hopping dominated transport in all samples, while localization lengths fall into three distinct categories depending on their band gaps. Armchair films have the largest localization length. Here our detailed analyses on electronic transport properties of single-chirality SWCNT films provide significant new insight into electronic transport in ensembles of nanoobjects, offering foundations for designing and deploying macroscopic SWCNT solid-state devices.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Revealing the unusual rate-dependent mechanical behaviors of nematic liquid crystal elastomers

Liquid crystal elastomers (LCEs) exhibit unique mechanical properties of soft elasticity and enhanced energy dissipation with rate dependency. They are potentially transformative materials for applications in mechanical impact mitigation and vibration isolation. However, previous studies have primarily focused on the mechanics of LCEs under equilibrium and quasistatic loading conditions. Critical knowledge gaps exist in understanding their rate-dependent behaviors, which are a complex mixture of traditional network viscoelasticity and the soft elastic behaviors with changes in the mesogen orientation and order parameter. Together, these inelastic mechanisms lead to unusual rate-dependent energy absorption responses of LCEs. In this work, we developed a viscoelastic constitutive theory for monodomain nematic LCEs to investigate how multiple underlying sources of inelasticity manifest in the rate-dependent and dissipative behaviors of monodomain LCEs. Here, the theoretical modeling framework combines the neo-classical network theory with evolution rules for the mesogen orientation and order parameter with conventional viscoelasticity. The model is calibrated with uniaxial tension and compression data spanning six decades of strain rates. The established 3D constitutive model enables general loading predictions taking the initial mesogen orientation and order parameter as inputs. Additionally, parametric studies were performed to further understand the rate dependence of monodomain LCEs in relation to their energy absorption characteristics. Based on the parametric studies, particularly loading scenarios are identified as conditions where LCEs outperform conventional elastomers regarding energy absorption.

42 ENGINEERING↗

DOE OSTI2024

A non-isothermal breakage-damage model for plastic-bonded granular materials incorporating temperature, pressure, and rate dependencies

Plastic-bonded granular materials (PBM) are widely used in industrial sectors, including building construction, abrasive applications, and defense applications such as plastic-bonded explosives. The mechanical behavior of PBM is highly nonlinear, irreversible, rate dependent, and temperature sensitive governed by various micromechanical attributions such as grain crushing and binder damage. This paper presents a thermodynamically consistent, microstructure-informed constitutive model to capture these characteristic behaviors of PBM. Key features of the model include a breakage internal variable to upscale the grain-scale information to the continuum level and to predict grain size evolution under mechanical loading. In addition, a damage internal state variable is introduced to account for the damage, deterioration, and debonding of the binder matrix upon loading. Temperature is taken as a fundamental external state variable to handle non-isothermal loading paths. The proposed model is able to capture with good accuracy several important aspects of the mechanical properties of PBM, such as pressure-dependent elasticity, pressure-dependent yield strength, brittle-to-ductile transition, temperature dependency, and rate dependency in the post-yielding regime. Furthermore, the model is validated against multiple published datasets obtained from confined and unconfined compression tests, covering various PBM compositions, confining pressures, temperatures, and strain rates.

Breakage↗

DOE OSTI2022

Observation of a multiplicity dependence in the $\mathcal{p}_T$-differential charm baryon-to-meson ratios in proton–proton collisions at $\sqrt{s}$ = 13 TeV

The production of prompt D 0 , D$^+_s$, and Λ$^+_c$ hadrons, and their ratios, D$^+_s$/D 0 and Λ$^+_c$/D 0 , are measured in proton–proton collisions at $\sqrt{s}$ = 13 TeV at midrapidity (|y| < 0.5) with the ALICE detector at the LHC. The measurements are performed as a function of the charm-hadron transverse momentum (p T ) in intervals of charged-particle multiplicity, measured with two multiplicity estimators covering different pseudorapidity regions. While the strange to non-strange D$^+_s$/D 0 ratio indicates no significant multiplicity dependence, the baryon-to-meson p T -differential Λ$^+_c$/D 0 ratio shows a multiplicity-dependent enhancement, with a significance of 5.3σ for 1 < p T < 12 GeV/c, comparing the highest multiplicity interval with respect to the lowest one. The measurements are compared with a theoretical model that explains the multiplicity dependence by a canonical treatment of quantum charges in the statistical hadronisation approach, and with predictions from event generators that implement colour reconnection mechanisms beyond the leading colour approximation to model the hadronisation process. The Λ$^+_c$/D 0 ratios as a function of p T present a similar shape and magnitude as the Λ/K$^0_s$ ratios in comparable multiplicity intervals, suggesting a potential common mechanism for light- and charm-hadron formation, with analogous multiplicity dependence. The p T -integrated ratios, extrapolated down to p T = 0, do not show a significant dependence on multiplicity within the uncertainties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Probing O 2 dependence of hydroperoxy-butyl reactions via isomer-resolved speciation

Degenerate chain-branching mechanisms of n-alkanes are centered on the formation of hydroperoxy-alkyl radicals ($\dot{Q}$OOH), formed via $\dot{R}$ + O 2 reactions, and the ensuing competition between unimolecular decomposition and second-O 2 -addition. Quantitative measurements of partially oxidized intermediates formed via reactions of $\dot{Q}$OOH provide critical constraints that are required for accurate modeling of combustion chemistry. To examine the influence of temperature and oxygen concentration on intermediates from unimolecular decomposition of $\dot{Q}$OOH, isomer-resolved speciation measurements were conducted on n-butane oxidation at 835 Torr in a jet-stirred reactor (JSR) from 500 – 900 K. Resulting from negative-temperature coefficient behavior, cyclic ether formation peaked at two temperatures, 650 K and 800 K, which were selected for separate experiments to quantify the O 2 -dependence of species profiles using O 2 concentrations of 4.2 · 10 17 – 1.1 · 10 19 molecules cm –3 . Utilizing vacuum-ultraviolet absorption spectroscopy and electron-impact mass spectrometry, cyclic ether isomers were quantified separately, including explicit resolution of cis– and trans– isomers of 2,3-dimethyloxirane. Stereoisomers of 2-butene were also quantified explicitly. For all cyclic ethers, a common trend in O 2 -dependence emerged: species concentrations reach a maximum near 3.0 · 10 18 molecules cm –3 (equivalence ratio of 0.5). Although quantitative disparities are evident, chemical kinetics modeling qualitatively reproduces the O 2 dependence of species at 650 K. However, at 800 K, weak dependence on O 2 is predicted, which is in contrast with the measurements. Two carbonyls, diacetyl and methyl vinyl ketone, were also quantified and follow similar dependence on [O 2 ] and temperature as the cyclic ethers, which indicates some fraction forms via $\dot{Q}$OOH-mediated reactions. The discrepancies between the measured and model-predicted species profiles indicate that sub-mechanisms for important intermediates may require additional elementary reactions, including stereochemical-specific reactions, to improve the fidelity of n-alkane combustion modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Size Dependence of Lattice Parameter and Electronic Structure in CeO2 Nanoparticles

Intrinsic properties of a compound (e.g. electronic structure, crystallographic structure, optical and magnetic properties) define notably its chemical and physical behavior. In the case of nanomaterials, these fundamental properties depend on the occurrence of quantum mechanical size effects and on the considerable increase of the surface to bulk ratio. However, the literature on this size-dependence and on the involved mechanisms is quite elusive and scarce. Here, we explore the size-dependence of both crystal and electronic properties of CeO2 nanoparticles (NPs) with different sizes by state-of-the art spectroscopic techniques. XRD, XPS and HERFD-XANES demonstrate that the as-synthesized NPs crystallize in the fluorite structure and they are predominantly composed of CeIV ions. The strong dependence of the lattice parameter with the NPs size was attributed to the presence of adsorbed species at the NPs surface thanks to FTIR and TGA measurements. In addition, the size-dependence of the eg level in the Ce LIII XANES spectra was experimentally observed by HERFD-XANES and confirmed by theoretical calculations.

Lanthanide, CeO2, HEFRD-XANES, Electronic structur↗

DOE OSTI2021

Influence of Ions on the Size Dependent Morphology of Aerosol Particles

The study of aerosol particles composed of mixtures of organic and inorganic compounds provides insights for understanding chemistry in the atmosphere as well as information about phase transitions of systems under confinement. In the submicron size regime, we have previously found that liquid-liquid phase separation can be inhibited at sufficiently small particle diameters leading to phase separated particles at larger sizes and homogeneous particles below a threshold diameter. In this paper, we have investigated the influence of cations and anions in the inorganic compound (NH 4 + , Na + , SO 4 2- , HSO 4 - , and Cl-) on the phase separation of submicron aerosol particles. Each of five salts were studied with two different organic compounds. Surprisingly, a strong dependence on the identity of the cation is evident in the size dependence of the particle morphology and no dependence on the anion was found. Sodium containing samples exhibit phase separation in particles below 20 nm in diameter whereas ammonium containing samples cease to undergo phase separation in particles between 45 and 65 nm in diameter. The separation relative humidity for supermicron droplets also depends on the identity of the cation in the inorganic component if the anion is the same across the salts compared. In addition, the correlation between the separation relative humidity and size dependence was investigated and displays a trend, but is not a dominant feature of the data. We explain the results in terms of hard and soft ions, where the harder cation (Na + ) leads to phase separation down to smaller sizes while the softer cation (NH 4 + ) prevents phase separation and causes particles up to a larger diameter to remain homogenous. This study has implications for atmospheric chemistry in regions dominated by sea spray aerosol, which contains sodium as the primary cation, when compared to continental aerosol, where ammonium is an abundant cation. Additionally, these findings can be used to understand the influence of cations on the phase transitions of confined materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Crystallographic Evidence of Size-Dependent Bond Flexibility in Metal–Organic Framework Nanocrystals

Size-dependent electronic, magnetic, and optical behavior suggests that metal–organic frameworks become softer materials as their particle sizes decrease, but direct evidence is lacking. Here, we report variable-temperature powder X-ray diffraction data of Fe(1,2,3-triazolate)2 particles that offer crystallographic insight into size-dependent bond flexibility. Rietveld refinement reveals size-dependent positive thermal expansion upon downsizing the crystalline domains from 178 to 9 nm, with a 6-fold increase from 16 MK –1 to 96 MK –1 . Here, this behavior occurs in tandem with size-dependent elongation of metal–ligand bonds and increasing thermal displacement parameters, consistent with pronounced metal-linker bond lability. We propose that these effects, as well as size-dependent annealing of crystallite sizes, originate from the high charge density and surface stress of smaller particles. Taken together, these results provide structural evidence that size reduction serves as a synthetic route to controlling the dynamic response of materials to external stimuli.

Crystals↗

DOE OSTI2023

Voltage-dependent excitation dynamics in UV-absorbing organic photovoltaics with efficient charge transfer exciton emission

Intermolecular charge-transfer excitons play a central role in determining the performance of organic solar cells as their voltage-dependent formation, dissociation, and recombination dynamics contribute to photocurrent generation, radiative/nonradiative voltage losses, and photovoltaic fill factor. Here, we explore the properties of brightly-emitting wide energy gap (>2 eV) charge transfer excitons by measuring the voltage-dependent photoluminescence, photocurrent, and ultrafast pump–probe transient absorption spectra of organic solar cells employing five UV-absorbing donor molecules that differ only by the length of the oligophenylene or acene group at their core. We find that organic solar cells with a strong correlation between their voltage-dependent photocurrent and charge-transfer exciton photoluminescence have low photovoltaic fill factors as they require voltage to facilitate efficient charge-transfer exciton dissociation. In contrast, solar cells that are efficient can readily generate charges without an applied field and have a separate population of tightly-bound charge-transfer excitons that are responsible for emission. Furthermore, considering that the sum of all excitation loss rates (i.e., recombination and charge extraction) must be equal to the excitation generation rate in the steady state, these voltage-dependent data allow us to solve for the voltage-dependent fate of all excitations in the solar cells and estimate upper and lower bounds for geminate and non-geminate recombination, respectively.

14 SOLAR ENERGY↗

DOE OSTI2022

Dissociation-induced depletion of high-energy reactant molecules as a mechanism for pressure-dependent rate constants for bimolecular reactions

In 1922, Lindemann proposed the now-well-known mechanism for pressure-dependent rate constants for unimolecular reactions: reactant molecules with sufficiently high energies dissociate more quickly than collisions can reestablish the Boltzmann distribution of the internal energies of the molecule during its dissociation at low pressures – yielding pressure-dependent rate constants for unimolecular reactions due to the preferential depletion of the high energy states capable of dissociation. In the last century, incredible progress has been made in achieving a far greater understanding of and quantitative predictions for unimolecular and association reactions. In the modern era, pressure-dependent phenomenological rate constants are now nearly universally used to describe the rates of unimolecular and associative reactions in phenomenological kinetic modeling. However, there is a second, more indirect, implication of Lindemann’s mechanism that relates to how these dissociation-induced non-equilibrium distributions impact bimolecular reactions, including non-associative bimolecular reactions – which are generally not considered to have pressure-dependent rate constants. Yet, as we show herein, the same high energy states depleted due to dissociation would otherwise react most rapidly in high-activation-energy bimolecular reactions – yielding a mechanism for pressure-dependent rate constants for bimolecular reactions (including non-associative reactions). Here, we present results from a case study for CH 2 O dissociation, isomerization, and bimolecular reaction with O2 to explore this question. Results from our master equation calculations indicate that the effect of dissociation-induced non-equilibrium distributions on bimolecular reactions can be substantial – even when chemical timescales are well separated from internal energy relaxational timescales (i.e. when the traditional rate constant description would be thought to apply). This effect is found to be more pronounced – and more complex – for bimolecular reactions involving molecular entities whose chemical timescales are merged with the internal energy relaxational timescales. Lastly, we present some ideas for discussion regarding what should be considered as “chemical species” in phenomenological kinetic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Assessing models of force-dependent unbinding rates via infrequent metadynamics

Protein–ligand interactions are crucial for a wide range of physiological processes. Many cellular functions result in these non-covalent “bonds” being mechanically strained, and this can be integral to proper cellular function. Broadly, two classes of force dependence have been observed—slip bonds, where the unbinding rate increases, and catch bonds, where the unbinding rate decreases. Despite much theoretical work, we cannot predict for which protein–ligand pairs, pulling coordinates, and forces a particular rate dependence will appear. Here, we assess the ability of MD simulations combined with enhanced sampling techniques to probe the force dependence of unbinding rates. We show that the infrequent metadynamics technique correctly produces both catch and slip bonding kinetics for model potentials. We then apply it to the well-studied case of a buckyball in a hydrophobic cavity, which appears to exhibit an ideal slip bond. Finally, we compute the force-dependent unbinding rate of biotin–streptavidin. Here, the complex nature of the unbinding process causes the infrequent metadynamics method to begin to break down due to the presence of unbinding intermediates, despite the use of a previously optimized sampling coordinate. Allowing for this limitation, a combination of kinetic and free energy computations predicts an overall slip bond for larger forces consistent with prior experimental results although there are substantial deviations at small forces that require further investigation. Furthermore, this work demonstrates the promise of predicting force-dependent unbinding rates using enhanced sampling MD techniques while also revealing the methodological barriers that must be overcome to tackle more complex targets in the future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Origin of the injection-dependent emission blueshift and linewidth broadening of III-nitride light-emitting diodes

III-nitride light-emitting diodes (LEDs) exhibit an injection-dependent emission blueshift and linewidth broadening that is severely detrimental to their color purity. By using first-principles multi-scale modeling that accurately captures the competition between polarization-charge screening, phase-space filling, and many-body plasma renormalization, we explain the current-dependent spectral characteristics of polar III-nitride LEDs fabricated with state-of-the-art quantum wells. Our analysis uncovers a fundamental connection between carrier dynamics and the injection-dependent spectral characteristics of light-emitting materials. For example, polar III-nitride LEDs offer poor control over their injection-dependent color purity due to their poor hole transport and slow carrier-recombination dynamics, which forces them to operate at or near degenerate carrier densities. Designs that accelerate carrier recombination and transport and reduce the carrier density required to operate LEDs at a given current density lessen their injection-dependent wavelength shift and linewidth broadening.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Space-time dependent thermal conductivity in nonlocal thermal transport

Nonlocal thermal transport is generally described by the Peierls-Boltzmann transport equation (PBE). However, solving the PBE for a general space-time dependent problem remains a challenging task due to the high dimensionality of the integro-differential equation. In this work, we present a direct solution to the space-time dependent PBE with a linearized collision matrix using an eigendecomposition method. We show that there exists a generalized Fourier-type relation that links heat flux to the local temperature, and this constitutive relation defines a thermal conductivity that depends on both time and space. Combining this approach with ab initio calculations of phonon properties, we demonstrate that the space-time dependent thermal conductivity gives rise to an oscillatory response in temperature in a transient grating geometry in high thermal conductivity materials. The present solution method allows us to extend the reach of our computational capability for heat conduction to space-time dependent nondiffusive transport regimes. Here, this capability will not only enable a more accurate interpretation of thermal measurements that observe nonlocal thermal transport, but also enhance our physical understanding of nonlocal thermal transport in high thermal conductivity materials that are promising candidates for nanoscale thermal management applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Unraveling the structural dependency of Weyl nodes in Co 2 MnGa

Conventionally, the modulation of the intrinsic Weyl nodes in Weyl semimetals is challenging, due to topological protection. Here we report the structural dependence of the Weyl nodes in a Co2⁢MnGa Heusler thin film via a temperature-dependent tetragonal distortion. The ability to manipulate these Weyl nodes allows for the control of the intrinsic electromagnetic properties. Temperature-dependent x-ray diffraction (XRD) measurements identify a compressive tetragonal distortion with decreasing temperature from 300 to 20 K. The calculated Weyl properties can be directly compared with experimental parameters through the temperature-dependent XRD measurements which show the intrinsic correlation between Weyl properties and important magnetic parameters. Further, the microscopic momentum space properties of Weyl nodes such as the distance (d W ), solid angle (Ω W ), tilt (φ W ), and nodal point energy (E W ) directly affect the macroscopic observable properties such as exchange stiffness (A), magnetization (M), and effective anisotropy field $H^{eff}_{K}$, as shown via structure-dependent density functional theory calculations. These predictions are experimentally observed as large variations in the bulk magnetization and effective anisotropy field as a function of temperature. These results highlight a unique degree of freedom in the control of macroscopic magnetic properties via the modulation of the intrinsic properties of Weyl nodes through structural distortions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Simple model to investigate jet quenching and correlated errors for centrality-dependent nuclear modification factors in relativistic heavy-ion collisions

Here, we apply Bayesian techniques to compare a simple, empirical model for jet quenching in heavy-ion collisions to centrality-dependent jet R AA measured by ATLAS for Pb + Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. We find that the R AA values for central collisions are adequately described with a model for the mean p T -dependent jet energy loss using only two parameters. This model is extended by incorporating two-dimensional initial geometry information from TRENTo and compared to centrality-dependent R AA values. We find that the results are sensitive to the value of the jet-quenching formation time, τ ƒ , and that the optimal value of τ ƒ varies with the assumed path-length dependence of the energy loss. We construct a covariance error matrix for the data from the p T -dependent contributions to the ATLAS systematic errors and perform Bayesian calibrations for several different assumptions for the systematic error correlations. We show that the most-probable functions and $χ^2_d$ values are sensitive to assumptions made when fitting to correlated errors. This work demonstrates the utility of a simple model that can quickly demonstrate the constraining power of jet-quenching observables with corresponding uncertainties and guide future studies using more sophisticated models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Computing nuclear response functions with time-dependent coupled-cluster theory

We compute nuclear response functions by solving the time-dependent 𝐴-body Schrödinger equation, recording the time-dependent transition moment and extracting spectral information via Fourier transforms. The solution of the time-dependent many-body problem accounts for correlations on top of the mean field by taking advantage of a time-dependent formulation of coupled-cluster theory. As a validation, we focus on electric dipole transitions in 4 He and 16 O and compare moments of the response function distribution to the results of an equivalent static framework, finding negligible discrepancies. We investigate how proton and neutron densities evolve in time, and we see the traditional picture of soft and giant dipole resonances as collective oscillations of protons and neutrons emerging from our calculations in 16 O and 24 O. Furthermore, this method also allows us to investigate the behavior of the nucleus in the presence of a strong electric field. In that regime, the behavior of the system becomes chaotic. Qualitatively, the spectral information obtained in this limit is in line with previous time-dependent mean-field results.

Ab initio calculations↗

DOE OSTI2020

Altered climate leads to positive density‐dependent feedbacks in a tropical wet forest

Abstract Climate change is predicted to result in warmer and drier Neotropical forests relative to current conditions. Negative density‐dependent feedbacks, mediated by natural enemies, are key to maintaining the high diversity of tree species found in the tropics, yet we have little understanding of how projected changes in climate are likely to affect these critical controls. Over 3 years, we evaluated the effects of a natural drought and in situ experimental warming on density‐dependent feedbacks on seedling demography in a wet tropical forest in Puerto Rico. In the +4°C warming treatment, we found that seedling survival increased with increasing density of the same species (conspecific). These positive density‐dependent feedbacks were not associated with a decrease in aboveground natural enemy pressure. If positive density‐dependent feedbacks are not transient, the diversity of tropical wet forests, which may rely on negative density dependence to drive diversity, could decline in a future warmer, drier world.

Bachelot, Benedicte↗

DOE OSTI2024

Impact of Pressure-Dependent Interfacial Tension and Contact Angle on Capillary Heterogeneity Trapping of CO2 in Storage Aquifers

Summary Carbon dioxide (CO2) capillary trapping increases the total amount of CO2 that can be effectively immobilized in storage aquifers. This trapping, manifesting itself as accumulated CO2 columns at a continuum scale, is because of capillary threshold effects that occur below low-permeability barriers. Considering that capillary pressure is dictated by heterogeneous pore throat size, the trapped CO2 column height and associated CO2 saturation will vary spatially within a storage aquifer. This variation will be influenced by two pressure-dependent interfacial parameters—CO2/brine interfacial tension (IFT) and CO2/brine/rock contact angle. Our objective is to understand how the pressure dependence of these two parameters affects the heterogeneity of capillary trapped CO2 at a continuum scale. Our conceptual model is a 1D two-zone system with the upper zone being a flow barrier (low permeability) and the lower zone being a flow path (high permeability). The inputs to this model include microfacies-dependent capillary pressure vs. saturation curves and permeability values. The input capillary pressure curves were collected in the literature that represents carbonate microfacies (e.g., dolograinstone) in a prevalent formation in the Permian Basin. We then used the Leverett j-function to scale the capillary pressure curve for the two zones that are assigned with the same or different microfacies. During scaling, we considered the influence of pressure on both the IFT and contact angle of CO2/brine/dolomite systems. We varied the zone permeability contrast ratio from 2 to 50. We then assumed capillary gravity equilibriums and calculated the CO2 saturation buildup corresponding to various trapped CO2 column heights. The CO2 saturation buildup is defined as the CO2 saturation in the lower layer minus that in the upper one. We found that the saturation buildup can be doubled when varying pressure in a storage aquifer, after considering pressure-dependent IFT and contact angles. Thus, assuming these two parameters to be constant across such aquifers would cause large errors in the quantification of capillary trapping of CO2. The whole study demonstrates the importance of considering pressure-dependent interfacial properties in predicting the vertical distribution of capillary trapped CO2. It has important implications in developing a better understanding of leakage risks and consequent storage safety.

Engineering↗