Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gas property variation”

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

Sensitivity of time-dependent density functional theory to initial conditions

Time-dependent density-functional theory is mathematically formulated through nonlinear coupled time-dependent three-dimensional partial differential equations, and it is natural to expect a strong sensitivity of its solutions to variations of the initial conditions, akin to the butterfly effect ubiquitous in classical dynamics. Since the Schrödinger equation for an interacting many-body system is, however, linear and mathematically the exact equations of the density-functional theory reproduce the corresponding one-body properties, it would follow that the Lyapunov exponents are also vanishing within a density-functional theory framework. Whether for realistic implementations of the time-dependent density-functional theory the question of the absence of the butterfly effect and whether the dynamics provided is indeed a predictable theory was never discussed. At the same time, since the time-dependent density-functional theory is a unique tool allowing us to study the nonequilibrium dynamics of strongly interacting many-fermion systems, the question of predictability of this theoretical framework is of paramount importance. Here our analysis, for a number of quantum superfluid many-body systems (unitary Fermi gas, nuclear fission, and heavy-ion collisions) with a classical equivalent number of degrees of freedom O(10 10 ) and larger, suggests that its maximum Lyapunov exponents are negligible for all practical purposes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Smart Manufacturing Pathways for Industrial Decarbonization and Thermal Process Intensification

Rapid decarbonization is fast becoming the primary environmental and sustainability initiative for many economic sectors. Industry consumes more than 30 % of all primary energy in the United States and accounts for nearly 25 % of all greenhouse gas (GHG) emissions. More than 70 % of energy consumed by the industrial sector is related to thermal processes, which are also the largest contributors of carbon emissions, overwhelmingly due to the combustion of fossil fuels. Thermal process intensification (TPI) seeks to dramatically improve the energy performance of thermal systems through technology pillars focusing on alternative energy sources and processes, supplemental technologies, and waste heat management. The impacts of TPI have significant overlap with the goals of industrial decarbonization (ID) that seeks to phase out all GHG emissions from industrial activities. Emerging supplemental technologies such as smart manufacturing (SM) and the industrial internet of things (IoT) enable significant opportunities for the optimization of manufacturing processes. Combining strategies for TPI and ID with SM and IoT can open and enhance existing opportunities for saving time and energy via approaches such as tighter control of temperature zones, better adjustment of thermal systems for variations in production levels and feedstock properties, and increased process throughput. Data collected by smart processes will also enable new advanced solutions such as digital twins and machine learning algorithms to further improve thermal system savings. Herein, this paper examines the individual pathways of TPI, ID, and SM and how the combination of all three can accelerate energy and GHG reductions.

42 ENGINEERING↗
DOE OSTI2022

Stresses within rare-earth doped yttria-stabilized zirconia thermal barrier coatings from in-situ synchrotron X-ray diffraction at high temperatures

There is a growing interest for smart coatings that can be integrated into turbine engines for in-situ temperature measurements or health monitoring. The addition of rare-earth dopants into standard thermal barrier ceramic top coat materials is used to obtain luminescent coatings that enable spectral measurements, for real-time temperature or health monitoring. The thermomechanical performance and durability of such novel coating compositions in extreme environments still remains to be evaluated. Consequently, the ability to manufacture sensor coatings which present suitable thermal properties needs to be demonstrated. For this study, highly luminescent erbium and europium doped yttria-stabilized zirconia and state-of-the-art yttria-stabilized zirconia coatings manufactured by air plasma spray were characterized to determine the effects of the embedded rareearth dopants on coating internal strain and stress and to quantify and compare their high temperature response using synchrotron X-ray diffraction. In-situ depth-resolved strain measurements were performed at 15 μm intervals along the depth of the coatings to evaluate materials response at key locations, specifically at layer interfaces. In-plane stress was calculated for the coatings and a finite element model was implemented to supplement the results and enable further predictions. The results show that the sensor coatings that were manufactured in this work revealed only minor variations in the strain response of sensor coatings under a typical thermal cycle and especially at temperatures closer to that of gas turbine operating conditions, compared to state-of-the-art coatings. This work demonstrates the viability of manufacturing rare-earth doped yttria-stabilized zirconia coatings that provide beneficial spectroscopic monitoring capabilities while having minimal impact on the thermomechanical response of the thermal barrier coatings.

Air plasma spray↗
DOE OSTI2025

High-Temperature Gas Sensor Materials with Properties Predicted via First-Principles Calculations with Machine Learning Modeling and Experimental Corroboration

Understanding the temperature dependence of functional properties of sensing materials is vital for their applications in combustion environments. The electron-phonon coupling that derives the electronic structure change with temperatures is a key property of interest as it affects other sensing responses. Herein, we first assess the temperature dependence of band gap renormalization in sensing materials by employing Allen-Heine-Cardona (AHC) theory with density functional theory (DFT) simulations corroborated with experimental observation. As the AHC calculations are impractical for high-throughput screening of materials, we employ data-driven Gaussian process regression to predict the parameters employed in the O’Donnell empirical model from a set of physical features. To mitigate the reliability issues arising from the small size of the dataset, we apply a Bayesian technique to improve the generalizability of the data-driven models as well as to quantify the uncertainty associated with theoretical predictions. These models capture well the overall trend of the O’Donnell parameters with respect to a reduced feature set obtained by transforming the available physical features. Quantifying the associated uncertainty helps us understand the reliability of the predictions and, therefore, the variation of bandgap as a function of temperature for other novel materials. The predicted candidates from machine learning models are further validated by experiments and DFT calculations.

bandgap renormalization↗
DOE OSTI2024

Identification of Suitable Vacuum Gas Oils as Plasticizers Using HT-GC × GC-HRMS

Vacuum gas oils (VGOs) have long been heralded as effective plasticizers due to their high boiling-points and lubricating properties. One of these VGOs, HyVac Oil 93050, is a paraffinic plasticizer that is currently utilized in polymer–plasticizer explosive formulations, while this oil effectively increases the elasticity and decreases the sensitivity for current formulations. We seek to identify alternative VGOs with similar density, viscosity, molecular composition, and impurities for future formulations. In this study, 18 VGOs, including HyVac Oil 93050, were initially evaluated for their density and viscosity. The oils were then ranked based on physical characteristics and analyzed for molecular composition using high-temperature comprehensive two-dimensional gas chromatography with high-resolution time-of-flight mass spectrometry (HT-GC × GC-HRMS). Further, the HT-GC × GC-HRMS chromatograms of the top 10 most similar VGOs to HyVac Oil 93050, according to density, were compared utilizing the Pearson correlation coefficient. Three of the oils with densities similar to those of HyVac Oil 93050 had a Pearson correlation within the lot-to-lot variation of HyVac Oil 93050. For the top 3 candidates, Pearson correlation was then utilized as a feature selection technique to discover significant chemical differences. Positive chemical ionization (PCI) and negative chemical ionization (NCI) HT-GC × GC-HRMS chromatograms were also evaluated and aided in the discovery and identification of several key compounds including additives that acted as stabilizers for the VGOs. After this further chemical analysis, two of the original 17 potential VGOs were determined to be physically and chemically similar to HyVac Oil 93050.

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

Elevated temperature contact creep and friction of nickel-based superalloys using machine learning assisted finite element analysis

Nickel-based superalloys with superior thermochemical, mechanical, and tribological properties are highly utilized for critical components in several high temperature applications such as gas turbines and nuclear reactors. Inconel 617, in particular, is considered as one of the main candidate superalloys for tribo-components in very-high-temperature gas-cooled nuclear reactors. Recent findings indicate that this alloy grows unique surface oxide especially in a high-temperature helium environment with distinctive wear, friction, and contact properties. Here, this study investigates the high temperature contact area evolution and frictional behavior of Inconel 617 using finite element simulation and provides predictive models for the contact and friction performance at different normal loads, dwell times, and temperatures. High temperature helium-aged Inconel 617 top surface properties (up to 600 °C) are utilized along with a single asperity-based deformable elastic-plastic contact model under combined normal and tangential loading. Machine learning is used to assist the finite element results and to predict friction coefficient as well as contact area evolution. While a small difference is observed in the instantaneous friction coefficient (no dwell time) for all temperatures, friction coefficient increases considerably with dwell time. This shows that the effect of contact creep for longer dwell times significantly dominants the effect of high temperature variation in basic mechanical parameters such as modulus and yield strength. It is found that increasing temperature and dwell times lead to the friction coefficient increase, yet the dominance of dwell time effects decreases at higher temperatures and loads. While the analysis is presented for Inconel 617, the methodology is easy to be generalized and can be applied to other HT alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Magnetic and electron transport properties of Co 2 Si nanomagnets

In this study, magnetotransport and ferromagnetism in thin films of Co 2 Si nanoclusters are investigated experimentally and theoretically. The nanoclusters are fabricated by an inert-gas condensation-type cluster-deposition method and have an average size of 11.3 nm. Unlike the bulk Co 2 Si that exhibits a very weak net magnetic moment only below 10 K, the nanoclusters exhibit room-temperature ferromagnetism with a substantial saturation magnetization. Key features of the system are its closeness to the Stoner transition, magnetic moments induced by spin polarization starting from surface atoms, and nonuniaxial anisotropy associated with the orthorhombic crystal structure of Co 2 Si . A method is introduced to determine the effective anisotropy using the experimental magnetization data of this complex system and its relationship with the two lowest-order nonuniaxial anisotropy constants. On decreasing temperature from 300 K, the nanoclusters show electron-transport properties unusual for a ferromagnetic metal, including an increase of Hall resistivity and a nonmonotonic change of negative magnetoresistance with a peak at around 100 K. The underlying physics is explained on the basis of the large polarization of surface spins and variation in the degree of their misalignments due to temperature-dependent effective anisotropy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermal and solubility effects on fault leakage during geologic carbon storage

Geologic carbon storage (GCS) is a promising method for reducing anthropogenic CO 2 emissions to the atmosphere. To safely deploy GCS in the field, it is necessary to assess risks and the effect of uncertainty on safe storage. The effect of uncertainty can be quantified using batches of simulations, but the high computational costs of high-resolution simulations necessitate use of reduced-order models (ROMs). Previous work involves ROMs for quantifying the risk of different potential leakage paths from storage reservoirs to shallow formations. However, previous studies on development of fault-leakage ROMs have limited numbers of uncertain parameters and do not explicitly examine impacts of CO 2 solubility and thermal stresses on fault reactivation, which can generate high-permeability pathways and compromise CO 2 storage. Here, we analyze an ensemble of simulations considering CO 2 leakage from a storage reservoir to a shallow aquifer through a fault while varying a number of uncertain parameters related to thermo-hydro-mechanical properties and CO 2 injection. We show the effects of solubility on: free-phase CO 2 -leakage rates, brine-leakage rates, and poroelastic fault destabilization. We find that CO 2 solubility is more important for estimating free-phase CO 2 -leakage rates compared to brine-leakage rates or poroelastic fault destabilization. We also find that thermal stresses and overpressures have different spatial distributions within the fault, indicating that the spatial variability of overpressures due to variation in flow parameters does not necessarily make the spatial variability of thermal stresses negligible. We suggest the use of the CO 2 phase-change path as a variable in future fault-leakage ROMs.

03 NATURAL GAS↗
DOE OSTI2021

Dynamic deuterium retention properties of tungsten measured using laser-induced breakdown spectroscopy

Dynamic deuterium (D) retention properties of tungsten (W) are investigated under steady-state D plasma exposure in the PISCES-A linear plasma device. In contrast to static retention, dynamic retention is quickly released from the material at the termination of the incident plasma flux, and thus in-operando laser-induced breakdown spectroscopy (LIBS) measurements have been conducted during steady-state D plasma exposure. A procedure is, first, established to extract the dynamic retention component from the in-operando LIBS D I 656.1 nm line intensity, which can contain contributions from dynamic and static retention as well as background D/D 2 gas, excited by both steady-state and laser-induced plasmas. Using the developed procedure, the dynamic D retention in W is systematically examined while scanning the following plasma exposure parameters: incident ion energy, E i , sample temperature, T s , incident ion flux, Γ i . No clear E i dependence is seen in the range of E i ~ 45–175 eV, as expected from a small variation in the average implantation depth (~3-6 nm) of D in W. The dynamic retention is found to monotonously decrease with increasing T s from 348 to 573 K, while the total static retention is reported to peak at T s ~ 500–600 K. It is revealed that the dynamic retention linearly increases with increasing Γ i , and then saturates at Γ i ≥ 0.75 × 10 21 m –2 s –1 . Possible physical mechanisms for the observed dependence of the dynamic retention on T s and Γ i are discussed. Furthermore, sequential pure He, followed by pure D, plasma exposures show that the dynamic D retention is not strongly affected by He bubbles in the near-surface region.

super-saturated surface layer↗
DOE OSTI2021

Spectroscopic and Photometric Monitoring of a Poorly Known Highly Luminous OH/IR Star: IRAS 18278+0931

We present the time-dependent properties of a poorly known OH/IR star, IRAS 18278+0931 (hereafter IRAS 18+09), toward the Ophiuchus constellation. We have carried out long-term optical/near-infrared photometric and spectroscopic observations to study the object. From optical R- and I-band light curves, the period of IRAS 18+09 is estimated to be 575 ± 30 days and the variability amplitudes range from ΔR ∼ 4.0 mag to ΔI ∼ 3.5 mag. From the standard period–luminosity relations, the distance (D) to the object, 4.0 ± 1.3 kpc, is estimated. Applying this distance in the radiative transfer model, the spectral energy distribution is constructed from multiwavelength photometric and IRAS-LRS spectral data, which provide the luminosity, optical depth, and gas mass-loss rate of the object to be 9600 ± 500 L {sub ⊙}, 9.1 ± 0.6 at 0.55 μm, and 1.0 × 10{sup −6} M {sub ⊙} yr{sup −1}, respectively. The current mass of the object is inferred to be in the range 1.0−1.5 M {sub ⊙} assuming solar metallicity. Notably, the temporal variation of atomic and molecular features (e.g., TiO, Na i, Ca i, CO, H{sub 2}O) over the pulsation cycle of the OH/IR star illustrates the sensitivity of the spectral features to the dynamical atmosphere as observed in pulsating AGB stars.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The chemical properties of the Milky Way’s on-bar and off-bar regions: evidence for inhomogeneous star formation history in the bulge

ABSTRACT Numerous studies of integrated starlight, stellar counts, and kinematics have confirmed that the Milky Way is a barred galaxy. However, far fewer studies have investigated the bar’s stellar population properties, which carry valuable independent information regarding the bar’s formation history. Here, we conduct a detailed analysis of chemical abundance distributions ([Fe/H] and [Mg/Fe]) in the on-bar and off-bar regions to study the azimuthal variation of star formation history (SFH) in the inner Galaxy. We find that the on-bar and off-bar stars at Galactocentric radii 3 kpc < rGC < 5 kpc have remarkably consistent [Fe/H] and [Mg/Fe] distribution functions and [Mg/Fe]–[Fe/H] relation, suggesting a common SFH shared by the long bar and the disc. In contrast, the bar and disc at smaller radii (2 kpc < rGC < 3 kpc) show noticeable differences, with relatively more very metal-rich ($\rm [Fe/H] \sim 0.4$) stars but fewer solar abundance stars in the bar. Given the three-phase star formation history proposed for the inner Galaxy in Lian et al., these differences could be explained by the off-bar disc having experienced either a faster early quenching process or recent metal-poor gas accretion. Vertical variations of the abundance distributions at small rGC suggest a wider vertical distribution of low-α stars in the bar, which may serve as chemical evidence for vertical heating through the bar buckling process. The lack of such vertical variations outside the bulge may then suggest a lack of vertical heating in the long bar.

Lian, Jianhui↗
DOE OSTI2021

Cometary Activity Begins at Kuiper Belt Distances: Evidence from C/2017 K2

We study the development of activity in the incoming long-period comet C/2017 K2 over the heliocentric distance range 9 ≲ r {sub H} ≲ 16 au. The comet continues to be characterized by a coma of submillimeter-sized and larger particles ejected at low velocity. In a fixed co-moving volume around the nucleus we find that the scattering cross section of the coma, C, is related to the heliocentric distance by a power law, C∝r{sub H}{sup -s}, with heliocentric index s = 1.14 ± 0.05. This dependence is significantly weaker than the r {sub H} {sup -2} variation of the insolation as a result of two effects. These are, first, the heliocentric dependence of the dust velocity and, second, a lag effect due to very slow-moving particles ejected long before the observations were taken. A Monte Carlo model of the photometry shows that dust production beginning at r {sub H} ~ 35 au is needed to match the measured heliocentric index, with only a slight dependence on the particle size distribution. Mass-loss rates in dust at 10 au are of order 10{sup 3} kg s{sup -1}, while loss rates in gas may be much smaller, depending on the unknown dust to gas ratio. Consequently, the ratio of the nongravitational acceleration to the local solar gravity, α', may, depending on the nucleus size, attain values of ~10{sup -7} ≲ α' ≲ 10{sup -5}, comparable to values found in short-period comets at much smaller distances. Nongravitational acceleration in C/2017 K2 and similarly distant comets, while presently unmeasured, may limit the accuracy with which we can infer the properties of the Oort cloud from the orbits of long-period comets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Heat transfer and pseudo phase transition for low-Reynolds, mixed-convection channel flow in the supercritical thermodynamic regime

Fluids at supercritical thermodynamic conditions are inherently complex due to the large variations in thermodynamic and transport properties. Multiple recent numerical and experimental investigations indicate ongoing interest for these fluids, especially supercritical CO and supercritical water, for a variety of applications. The potential heat transfer benefits in this regime are directly influenced by the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. Limited studies have been conducted for low-Reynolds/intermediate-Rayleigh numbers mixed-convection channel flow of supercritical water, which is the focus of this study. To investigate the thermally driven hydrodynamic instabilities in this regime, we use a high-order fully-implicit numerical method. As expected, in this channel configuration and flow regime, a competition arises between the forced and the natural convection. Buoyancy, driven by the thermal gradient at the wall, pushes the fluid upwards in the channel, while inertial forces drive the fluid from inlet to outlet. We observe significant differences in the heat transfer coefficient between the top and bottom walls of the channel flow, due to natural convection, and unstable plumes. We further quantify the buoyancy forces along the channel flow through evaluation of the Richardson number. Simulations show that with a higher negative Richardson number, there are more unstable structures along the channel, resulting in enhanced mixing. Furthermore, we demonstrate the existence of gas-like, liquid-like, and the possibility of pseudo-two-phase regimes in the supercritical fluid. The pseudo-two-phase regime occurs when the range of fluid temperatures is crossing the pseudo-critical line, also known as the Widom line. The density of the fluid varies significantly near the Widom line, and further influences the natural convection and heat transfer in the flow.

42 ENGINEERING↗
DOE OSTI2023

Effect of material composition of diamond field emission array cathodes on quality of transversely shaped beams

Diamond field emitter array field emission cathodes (DFEA FECs) are attractive for the next generation of injectors due to their ability to produce transversely shaped beams without the need for complex masking or laser schemes. However, reliability of this cathode technology remains a challenging issue as principal mechanisms guiding and allowing for output beam shaping remained poorly understood. Here, this paper reports the results of testing two DFEA FECs with the same pattern and emitter tip geometry. Although both cathodes were able to sustain gradients of 44 MV/m and produce maximum output integral charge of 0.5 nC per radio frequency pulse, their emission patterns were different. One cathode did not produce a shaped beam, while the other one did. This difference was explained by the extent of the local variations of the diamond material across the arrays as discovered by spatially resolved Raman spectroscopy. The main practical takeaways were (i) tip sharpness was not a prerequisite for producing a shaped beam and instead (ii) material characteristics resulting in different cathode ballast resistance affected emission spatial uniformity across the array and hence the beam shaping.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Investigating Performance of Geothermal Reservoirs With Bi-Directional Heat Transfer for Seasonal Thermal Storage Using Thermal-Hydraulic Reservoir Simulations

This paper presents initial simulation results of deep subsurface seasonal thermal storage of waste heat from a natural gas-fired power plant during the summer months for providing heating during winter months at the Cornell University main campus in Ithaca, New York. An injection temperature profile is calculated with the MEnU tool and numerical reservoir simulations are conducted with the simulator FALCON. Both a porous medium and large single fracture geothermal reservoir are considered at about 3 km depth with initial reservoir temperature of 75 degrees C. Simulations indicate that a much smaller reservoir - which is technically easier and less risky to develop - is sufficient when operating the reservoir for seasonal thermal storage instead of for continuous heat production. For example, storing 13 GWh of heat seasonally in a 100 m thick porous medium based on the Galway formation requires about 65 m of well spacing to limit the yearly variation in production temperature to 10 degrees C. Operating this reservoir for continuous heat production would require a well spacing of 365 m to limit the thermal decline to 10 degrees C after 10 years. Neglecting wellbore heat losses and assuming uniform and homogeneous properties, a round-trip thermal efficiency of 100% can be obtained if the yearly averaged injection temperature matches the initial reservoir temperature.

district heating↗
DOE OSTI2025

Summer Aerosol and Trace Gas Observations in Houston, Texas Using an Adaptable Mobile Facility

An aerosol container featuring a shared inlet system was deployed to Houston, Texas in July 2022, enabling direct, high-time-resolution in situ measurements of aerosols and trace gases. The internal rack system and floorplan was designed for adaptable modularity to elucidate aerosol physicochemical processes at fine scales. The design allowed for the deployment of a core instrument suite and additional customized research grade instruments. A heterogeneous mixture of aerosols was observed during three regimes: (1) intermittent black carbon (BC) and diurnal variations in aerosol chemical composition, (2) observed particle growth associated with SO 2 , (3) transported supermicron dust. The high variability of observed particles and gases in high time resolution indicated a complex urban area with multiple local and regional sources and processes. Particle growth rates of 7–16 nm/hr were observed for submicron particles during periods when SO 2 was >0.5 ppbv. Two periods of multi-day long-range transport events of dust from the African Sahara were observed in the supermicron and submicron particle modes with total mass concentrations up to 30 μg m −3 . Aerosol scattering angstrom exponents and extinction coefficients (B ext ) increased with humidity as a function of particle composition. The measurements demonstrate collaborative capabilities that can be used to increase observations of aerosol processing, microphysical and optical properties, internal mixing state, and supermicron aerosol that are not parameterized or missing in global Earth energy system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Assessing the Effects of Temperature and Oxygen Vacancy on Band Gap Renormalization in LaCrO 3-δ : First-Principles and Experimental Corroboration [plus supplemental information]

Understanding the temperature dependence of functional properties in high-temperature gas sensors is vital for applications in combustion environments. Temperature effect on the electronic structure due to electron-phonon coupling is a key property of interest as this influences other responses of sensors. In this work, we assess the impact of temperature on band gap renormalization of pristine and oxygen-vacant LaCrO 3–δ perovskite employing Allen–Heine–Cardona theory with first-principles simulations and corroborate with experimental observation. Antiferromagnetic cubic LaCrO 3 shows a direct ground-state band gap of 2.62 eV that is reduced by over 1 eV due to the presence of oxygen vacancies, which can form endothermically. We find excellent agreement in temperature-dependent band gap shift in LaCrO 3 between theory and an in-house experiment, proving that the theory can adequately predict renormalization on the band gap in a magnetic system. Band gaps in cubic LaCrO 3–δ are found to monotonically narrow by 1.13 eV in pristine and by around 0.62 eV in oxygen-vacant structures as temperature increases from 0 to 1500 K. The predicted band gap variations are rationalized using an analytical model. Furthermore, the experimental zero-temperature band gaps are extracted from the model fits that can provide useful insights on the simulated band gaps.

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

On the pressure and temperature dependence of adsorption densities and other thermodynamic properties in gas shales

Black shale source rocks have become a key natural gas resource in the US and China. Unlike in conventional hydrocarbon reservoirs, a significant fraction of the gas in shales is adsorbed onto the inordinate specific surface areas (SSA) of micropores inside organic matter with densities exceeding those of bulk gas by an order of magnitude. Yet, the distribution of adsorbed versus bulk gas remains poorly understood. Experiments measure excess adsorption isotherms, which are corrected to absolute adsorption amounts by a factor that involves the density, &#x3C1; abs , of the adsorption layer(s). Constant values of &#x3C1; abs are widely used in the literature but, we argue, are inconsistent with the pressure- and temperature-dependence of adsorption behavior. In this work, we propose a more reliable approach that assumes a constant volume of each adsorption layer that can be determined from the SSA of the substrate. We measure nine high-pressure methane excess adsorption isotherms up to 15 MPa for three shale samples from China and Germany at temperatures of 65, 75, and 95 °C, and low-pressure carbon dioxide and nitrogen isotherms to estimate the SSA. We also investigate another 24 isotherms at 35, 50, and 60° C or 8 Chinese shale samples and 3 isotherms at 60, 90, and 120° C for an Argentinian shale sample from the literature. A new variation of the multi-layer Ono-Kondo lattice model is introduced to derive the temperature and pressure dependent adsorption layer densities, which are subsequently used to obtain the absolute adsorption isotherms of interest. Finally, this practical methodology has the potential to significantly improve our estimates of shale gas-in-place resources, and adsorption modeling more broadly.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗