SEARCH · Engineering Papers
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.
Farnsworth Unit (FWU) CO2EOR Eclipse compositional model: Hydrostratigraphic Region 4 model uses the Morrow2 (UNOCAL) relative permeability curve with capillary pressure assigned heterogeneously by hydrostratigraphic unit
Spontaneous imbibition in an angular tube with corner films is a fundamental problem in many scientific and engineering processes. In this paper, a modified interacting capillary bundle model is developed to describe the liquid imbibition dynamics in a square tube with corner films. The square tube is decomposed into several interacting sub-capillaries and the local capillary pressure in each sub-capillary is derived based on the specific shape of its meniscus. The conductance of each sub-capillary is calculated using single-phase lattice Boltzmann simulation. The modified interacting capillary bundle model and color-gradient lattice Boltzmann method are used to simulate the liquid imbibition dynamics in the square tube with different fluid properties. The predictions by the modified interacting capillary bundle model match well with the lattice Boltzmann simulation results for different conditions, demonstrating the accuracy and robustness of the interacting capillary bundle model to describe the imbibition dynamics with corner films. In addition, the interacting capillary bundle model is helpful to investigate the mechanisms during spontaneous imbibition and the influences of fluid viscosity, surface tension, wetting phase contact angle and gravity on imbibition dynamics. Finally, a universal scaling law of imbibition dynamics for the main meniscus is developed and the scaling law for arc meniscus is also analyzed.
Capillary Pressure to Relative Permeability water-wet algorithm.
Farnsworth Unit (FWU) CO2EOR Eclipse compositional model: Hydrostratigraphic Region 1 model uses the HS1 relative permeability curves assigned heterogeneously by hydrostratigraphic unit. The model with capillary pressure applies the HS1 curves heterogeneously by hydrostratigraphic unit.
These excel spreadsheets were used to develop the two-phase and three-phase methods for converting capillary pressure to relative permeability. Data was collected by the SWP and literature sources cited in the spreadsheets.
Rough surfaces in contact with a flow of fluid exhibit alternating no-slip and free shear boundary conditions at the solid–liquid and air–liquid interfaces, respectively, thereby potentially offering drag reduction benefits. The balance between the dynamic pressure in the flow and the restoring capillary pressure in the interasperity spaces determines the stability of the Cassie state of wettability and is a function of the relative extent of no-slip and free shear regions per unit surface area. In the present study, using a fractal representation of rough surface topography, an analytical model is developed to quantify the stability of the Cassie state of wettability as well as drag reduction and the friction factor for laminar flow in a rectangular channel between nonwetting multiscale rough surfaces. A systematic study is conducted to quantify the effects of fractal parameters of the surfaces and the flow Reynolds number on drag reduction and the friction factor. The studies are used to develop friction factor curves extending the classical Moody diagram to hydrophobic and superhydrophobic surfaces. On the basis of the studies, regime maps are derived for estimating the extent of drag reduction offered by hydrophobic and superhydrophobic surfaces, revealing that superhydrophobic surfaces do not always offer the best drag reduction performance. The application of the fractal model to practical topographies of nonwetting surfaces of copper, aluminum, and zinc oxide fabricated via electrodeposition and etching is also discussed.
Here, this paper investigates the validity of the interpretation of results from testing saturated axisymmetric triaxial compression (ATC) sand specimens utilising three-dimensional (3D) synchrotron micro-computed tomography (SMT) to probe localised events that are completely missed or misinterpreted when analysing ATC measurements based on global standard measurements. Drained and undrained experiments were conducted at low and high back-pressures (BPs) coupled with multiple in situ 3D SMT to acquire high-resolution scans of the specimens at different axial strains. Specimens tested under low BP exhibited a large pore air volume change, which was not detected by the pump system that represents standard volume measurement. The increase in air volume caused a significant reduction in the degree of saturation leading to a possible transition from saturated to partially saturated constitutive behaviour. Undrained experiments exhibited a significant volume change contrary to the assumption of negligible volumetric strain for saturated undrained experiments. Air bubbles within the shear band for drained and undrained low-BP specimens showed opposite capillary pressure responses, increased for drained and decreased for undrained cases, due to the variation in the mechanism by which each of the two experiments predominantly counters the volume expansion within the shear band.
The proximity of the order of magnitude molecular mean free path and pore sizes in nanopores leads to remarkable interactions between molecules and walls. In such systems, the thermodynamic property and the transport behavior of fluids deviate from those at bulk conditions. Molecular dynamics (MD) simulation may be used to investigate the effects of confinement on fluid physics in nanosize pores. However, translating these subscale observations into a larger scale of interest remains a challenging task. In this work, we propose a modified extension of Peng-Robinson equation of state (PR-EOS) motivated and guided by MD simulation results. The shift of critical properties, i.e., pressure and temperature, are evaluated independently. A temperature dependent parameter is introduced to account for the existence of capillary pressure in the two-phase region. Here, this formulation is capable of capturing the shift in critical properties as well as density phase diagrams under various confinement scenarios. We incorporate the proposed EOS in the lattice Boltzmann method (LBM) to study coupled confined phase behavior and transport of methane in nano-size slit pores. Adsorptive strength is determined such that the resulting density ratio matches that from MD simulation. Results indicate that the use of the proposed EOS leads to smaller fluid velocities compared to PR-EOS. This effect is due to stronger interactions between fluid particles under confinement. Also, results show that transport characteristics are impacted by the pressure of the system: in systems with a relatively low pressure, transport seems to be dominated by Knudsen diffusion; however, at a higher pressure, the contribution from viscous flow increases as the pore widens while the influence of Knudsen and surface diffusion diminishes.
A potential geologic target for CO 2 storage should ensure secure containment of injected CO 2 . Traditionally, this objective has been achieved by targeting reservoirs with overlying seals-regionally extensive, low permeability units that have been proven capable of retaining buoyant fluid accumulations over geologic time. However, considering that the amount of CO 2 is limited by a decadal injection period, vertical migration of CO 2 can be effectively halted by a composite system of discontinuous shale/silt/mudstone barriers in bedded sedimentary rocks. Here, we studied the impact of depositional architectures in a composite confining system on CO 2 migration and confinement at reservoir scale. We stochastically generated lithologically heterogeneous reservoir models containing discontinuous barriers consistent with statistical distributions of net-sand-to-gross-shale ratio (NTG) and horizontal correlation lengths derived from well log data and observations of producing hydrocarbon fields in Southern Louisiana. We then performed an extensive suite of reservoir simulations of CO 2 injection and post-injection to evaluate the sensitivity of CO 2 plume migration and pressure response of the composite system to a series of geologic and fluid parameters including the lateral continuity of barriers, NTG, permeability anisotropy within the sand body, and capillary pressure contrast between the sand and shale facies. The results indicate that lateral continuity of barriers and NTG are the dominant controls on CO 2 plume geometry and pressure build-up in the reservoir, while the impact of NTG is particularly pronounced. The significance of intraformational barriers becomes apparent as they facilitate the local capillary trapping of CO 2 . Those barriers improve the pore space occupancy by promoting a more dispersed shape of the plume and ultimately retard the buoyancy-driven upward migration of the plume post injection.
A number of tasks for the heat pipe simulation tool Sockeye were identified based on customer needs and included both user documentation and code development: (a) Implement saturation-dependence in the capillary pressure model; (b) Provide an analytic utility to analyze the frozen startup limit; (c) Demonstrate custom friction closures; (d) Develop examples for different types of heat transfer sources and boundary conditions; (e) Demonstrate the capillary and sonic limits; and (f) Provide discretization details in the theory manual. These tasks were all completed, and results are given where applicable.
Ultrafast water transport in carbon nanotubes (CNTs) has drawn a great deal of attention in a number of applications, such as water desalination, power generation, and biomolecule detection. With the recent experimental advances in water filling of isolated CNTs, the Lucas-Washburn theory for capillary rise in tubes needs to be revisited for a better understanding of the physics and dynamics of water filling in CNTs. Here, the Lucas-Washburn theory is corrected for the hydrodynamic entrance effects as well as the variation of capillary pressure and hydrodynamic properties with the radius and length of CNTs. Due to the large slippage in CNTs, inclusion of the entrance effects is important particularly for the initial stages of filling where a L ∝ t scaling, as opposed to L 2 ∝ t , is observed in our molecular dynamics (MD) simulations. The corrected Lucas-Washburn theory is shown to predict the water filling dynamics in CNTs as observed in MD simulations. With the corrected theory, we achieve a better understanding of capillary rise and water filling dynamics in CNTs.
Two-phase fluid flow properties underlie quantitative prediction of water and gas movement, but constraining these properties typically requires multiple time-consuming laboratory methods. The estimation of two-phase flow properties (van Genuchten parameters, porosity, and intrinsic permeability) is illustrated in cores of vitric nonwelded volcanic tuff using Bayesian parameter estimation that fits numerical models to observations from spontaneous imbibition experiments. The uniqueness and correlation of the estimated parameters is explored using different modeling assumptions and subsets of the observed data. The resulting estimation process is sensitive to both moisture retention and relative permeability functions, thereby offering a comprehensive method for constraining both functions. The data collected during this relatively simple laboratory experiment, used in conjunction with a numerical model and a global optimizer, result in a viable approach for augmenting more traditional capillary pressure data obtained from hanging water column, membrane plate extractor, or mercury intrusion methods. This method may be useful when imbibition rather than drainage parameters are sought, when larger samples (e.g., including heterogeneity or fractures) need to be tested that cannot be accommodated in more traditional methods, or when in educational laboratory settings.
Accessibility of pore space in unconventional reservoirs is an important factor influencing both methane storage capacity and the kinetics of methane desorption. The determination of open (accessible) versus closed (inaccessible) porosity is therefore vital for the prediction of gas production potential. This study investigates accessibility of pores to methane in overmature middle Devonian Marcellus shale samples (cut parallel and perpendicular to bedding) using small and ultra-small neutron scattering (SANS and USANS) with contrast matching (CM), supplemented by other complementary techniques, such as mercury injection capillary pressure (MICP) and low pressure gas (N2 and CO2) adsorption. Our results demonstrate that for the samples studied, only about 6 % of pores with diameter 25-500 nm are accessible to methane. The accessibility fraction for pores larger than 500 nm is 35%. For nanopores smaller than 25 nm, pore accessibility could not be quantitatively determined due to increased methane density and condensation effects in confinement. Our observations indicate that methane penetrates the accessible small mesopores and micropores down to at least 1 nm in diameter, and the density of confined deuterated methane (CD4) is 0.68 g/cm3 for pores of diameter 25 nm and gradually increases with the decreased pore size. Moreover, elevated gas pressure causes formation of additional high-density methane nano-clusters. These clusters have a form of slightly anisotropic polydisperse discs oriented along bedding plane, about 1-12 nm in diameter and with average thickness of 3.6 nm. Utilizing samples cut parallel and perpendicular to the bedding, this study also briefly addresses anisotropy of pores. Based on the iso-size intensity ratio R, our SANS and USANS results demonstrate anisotropy in the out-of-bedding direction and suggest that the degree of anisotropy depends on the pore size. Specifically, for pore diameters ~2.5 to 250 nm, the extent of anisotropy is smaller than for pores ~500 nm to 6 µm in diameter. Finally, comparison of pore size distribution results calculated from SANS/USANS to those obtained using MICP shows good agreement at low pressures, but large difference at pressures above 1000 bar. Here, this discrepancy requires further testing; it is possible that the high mercury pressure used in MICP alters the mesopore and micropore structure of shales.
Within the Haynesville-Bossier Shale complex, the Bossier Shale has not been extensively studied by either industry and academia, despite it being an unconventional gas reservoir and a potential caprock for carbon storage in the underlaying Haynesville Shale. The lack of knowledge of the complex pore structures and fluid-rock interactions hinders the effective extraction of gas and the characterization of fluid reservoirs and sealing capacity. Integrated experimental studies of pore structure and fluid-rock interactions were conducted in seven Bossier Shale core samples collected in eastern Texas. Petrographic, geochemical, and petrophysical properties such as mineral composition, organic richness, thermal maturity, porosity, pore/pore throat diameter distribution, water-accessible pores, liquid water imbibition, and water vapor adsorption were characterized using complementary approaches of thin-section petrography, scanning electron microscopy, X-ray diffraction, total organic matter, pyrolysis, mercury intrusion porosimetry, nuclear magnetic resonance, (Ultra-) small angle X-rays scattering as well as small angle neutron scattering with deuterated liquids and contrast variation. Further, the results show that the thermally mature Bossier Shales are composed of mixed argillaceous mudstone, mixed mudstone, and mixed carbonate mudstone. The shale contains both organic and inorganic pores, with porosities of 3.24-9.37 %, pore-to-throat ratios of 1.65 to 19.4, and water-accessible pores accounting for 28.7-72.6 % of total pores. Approaches of liquid water imbibition and water vapor adsorption, with and without direct contact of water with shale samples, indicate that liquid water first enters the nano-sized pores under high capillary pressures, and water vapor adsorption is mainly controlled by both clay minerals and pores with diameters less than 10 nm. These findings contribute to a better understanding of pore structures and water-shale interactions and their controlling factors in the Bossier Shale.
Integration of petrophysical and geological information is critical to simulation of subsurface carbon storage (GCS). In this sense, two depositional facies were identified from the core description and well-log interpretation, namely massive (MS) and cross-bedded (CB) facies groups. Additionally, pore-scale characteristics were studied by a combination of techniques, e.g. Nuclear Magnetic Resonance (NMR) and mercury intrusion capillary pressure (MICP). Scanning electron microscope (SEM) and petrographic analyses show that the pore structure is dominantly controlled by the depositional environment and dolomite cementation. NMR-T 2 distributions of MS and CB facies show triple and quadruple modes, respectively. In addition, MICP of high- and low-permeability MS facies samples, and their CB facies group mixtures were collected. The MS sample pore-throat size distribution is uni-modal, while the triple-modal characteristic of the mixtures indicates heterogeneous pore structures at the sub-core scale for CB facies. The reliably estimates of porosity and permeability for both facies groups via NMR techniques and the MLR (Multiple Linear Regression) approach demonstrate the applicability of these techniques to eolian sandstone. Moreover, irreducible water saturation via the T 2 -cutoff method correlates strongly with T 2LM instead of porosity. Finally, the rock quality index and flow zone indicator were calculated based on Combinable Magnetic Resonance (CMR) log interpretations. This provides direct connection to properties measured in the well. Four flow units were classified for both facies groups. Finally, results show that better reservoir quality with significant heterogeneities is observed in the CB facies. This study highlights the importance integrating multiscale petrophysical properties including facies, pore architecture and diagenesis analysis with core- to log-scale property characterization. The results herein validate our reservoir characterization and flow unit classification in eolian reservoirs.
Cycling of lithium-ion batteries containing Ni-rich NMC cathodes at high voltage involves intense gas generation. From the safety standpoint, it is critical to understand how different gas species respond to changes of upper cut-off voltages. In this manuscript, we introduce a novel experimental set up for real-time analysis of gas generation in prismatic pouch cells. In a typical experiment, a lithium-ion pouch cell is directly connected to a quadruple mass spectrometer by glass capillary. Pressure difference helps move the generated gases to the mass spectrometer column for full analysis. The gaseous species are probed during both formation cycle and aging cycles in Li-ion pouch cells comprising NMC811 cathode and graphite anode. The gases are generated upon the creation of the solid electrolyte interphase on graphite during the first formation charge. Besides ethylene gas, C 2 H 4 , whose occurrence is related to the decomposition of ethylene carbonate, a minor solvent decomposition happened at a lower voltage is demonstrated to be associated with residue water on the electrode. For aging cycles, three upper cut-off voltages have been selected to investigate the gas evolution behavior upon charging to 4.2, 4.4 and 4.6 V. Higher upper cut-off voltages do not affect the amount of O 2 measured as a part of the generated oxygen could have reacted with electrolytes. Quite similar and low amounts of O 2 remain in the cells upon charging to the selected cut-off voltages. On the other hand, the generation of CO 2 has been found to be very sensitive to upper cut-off voltage. During similar aging cycles, 6167 nmol $g_{NMC}^{-1}$ of CO 2 is generated at 4.6 V, versus 1650 nmol $g_{NMC}^{-1}$ of CO 2 at 4.4 V, versus 91 nmol $g_{NMC}^{-1}$ of CO 2 at 4.2 V.
Proton exchange membrane fuel cells (PEMFC) with cathodes using platinum group metal-free (PGM-free) catalysts could significantly reduce costs, but the lower volumetric oxygen reduction reaction (ORR) activity requires thick electrodes that suffer from liquid water flooding and increased oxygen transport resistance. To address these challenges, we developed a 3D gas diffusion electrode (GDE) architecture to enhance liquid water removal through the diffusion media and reduce cathode saturation. The cathode features a uniform catalyst layer adjacent to the membrane for high ORR activity and then pillars of the hydrophilic catalyst layer that pass through the microporous layer (MPL), providing a low capillary pressure barrier pathway to the carbon fiber paper layer and channel. The non-planar cathode also increases the interfacial area between the catalyst layer and hydrophobic MPL for a greater fraction of the cathode with high O 2 concentration. Our studies included parametric experimental study of the pillar density to identify the optimum pitch between pillars. Our measurements show significant improvements in the mass transport region of the polarization curve with 3D structured electrodes leading to an 8% increase in maximum current density, 19% increase in maximum power density, and 16% increase in current density at 0.67 V with air.
Molecular diffusion plays an important role in oil and gas migration in tight shale formations. However, there are insufficient reference data in the literature to specify the diffusion coefficients within porous media. This study aims at calculating diffusion coefficients of shale gas, shale condensate, and shale oil at reservoir conditions with CO 2 injection for EOR/EGR. The large nano-confinement effects including large gas-oil capillary pressure and critical property shifts could alter the phase behaviors. This paper estimates the diffusivities of shale fluids in nanometer-scale shale rock from two perspectives: 1) examining the shift of diffusivity caused by nanopore confinement effects from phase change (phase composition and fluid property) perspective, and 2) calculating the effective diffusion coefficient in porous media by incorporating rock intrinsic properties (porosity and tortuosity factor). The tortuosity is obtained by using tortuosity-porosity relations as well as the measured tortuosity of shale from 3D imaging techniques. The results indicated that nano-confinement effects could affect the diffusion coefficient through altering the phase properties, such as phase compositions and densities. Compared to bulk phase diffusivity, the effective diffusion coefficient in porous shale rock is reduced by 10 2 to 10 4 times as porosity decreases from 0.1 to 0.03.