Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “isotherm”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The energy saving potential of thermo-responsive desiccants for air dehumidification

The desiccant wheel is a promising technology for energy-efficient humidity control. However, its overall efficiency—hindered by the desiccant materials’ properties—remains low because traditional desiccants (e.g., silica gels) have a single isotherm regardless of their adsorption temperature. Thermo-responsive materials have been proposed to break this fixed affinity to water vapor, with drastically different adsorption isotherms depending on temperature. Its potential for improving dehumidification efficiency, however, has not been addressed. In this paper, we model the potential of a thermo-responsive interpenetrating polymer network (IPN) desiccant with temperature-dependent adsorption isotherms for energy-efficient dehumidification via a validated transient desiccant wheel model for humidity control of buildings. Thermo-responsive desiccants can improve the dehumidification performance due to their thermo-responsive switchable hydrophilicity below/above the lower critical solution temperature. Our analysis shows that thermo-responsive IPN desiccants can potentially reduce energy consumption by up to 30% compared to silica gels. Finally, the savings depend strongly on the critical temperature of the thermo-responsive desiccant and should be higher when the inlet temperatures are expected to be higher.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Influence of alloy solidification path on melt pool behavior in additive manufacturing

Numerical models used to study transport phenomena in laser-powder bed fusion processes often rely on assumptions and simplifications to reduce their computational expense. One common simplification is in the description of latent heat evolution during the solid-liquid phase change (i.e., the solidification pathway), justified by the fact that the mushy zone thickness is similar to the numerical grid spacing used for continuum transport models. The lack of resolution of transport phenomena in the mushy zone motivates the use of computationally convenient solidification paths such as linear or sigmoidal relationships over pathways derived from fundamental solidification theory such as equilibrium or Scheil models. In the present work, an uncertainty quantification (UQ) framework is used to analyze the influence of solidification pathway selection on the solidification dynamics and melt pool geometries in laser based additive manufacturing (AM) of IN625. Results show the solidification pathway has a quantifiable influence on the cooling rate at the liquidus isotherm, mushy zone thickness, and solidification time. Due to similarities in the latent heat evolution at the beginning of solidification, the equilibrium and Scheil models predict similar cooling rates near the liquidus isotherm, however the wider freezing range of Scheil leads to a wider mushy zone compared to equilibrium. The non-physical latent heat release profiles of sigmoidal and linear paths lead to significant overpredictions of cooling rates at the liquidus isotherm compared to equilibrium and Scheil. Finally, these results indicate that careful consideration should be given to the choice of solidification pathway to ensure reliable model predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Determination of the enthalpy of adsorption of hydrogen in activated carbon at room temperature

The development of high-performance materials for hydrogen storage by adsorption requires detailed understanding of the adsorbate-adsorbent interactions, e.g., the enthalpy of adsorption ΔH, which measures the interaction strength. The determination of ΔH for a weakly adsorbing gas such as hydrogen in a carbonaceous porous material is difficult experimentally, normally requiring measuring two cryogenic adsorption isotherms. Here we demonstrate a calculation of ΔH based on ca. room temperature adsorption isotherms at 273 K and 296 K using the Clausius-Clapeyron equation. This requires an estimation of the volume of the adsorbed film (~40%, ~12% of the total pore volume at 77 K, 296 K, respectively) obtained from fits of the excess adsorption isotherms to an Ono-Kondo model with the auxiliary use of a fixed point corresponding to the saturation film density (estimated as 100 ± 20 g/L) which appears to be remarkably sample and temperature independent, i.e., a property of the adsorbate. The calculated room temperature enthalpy of adsorption ΔH = 8.3 ± 0.4 kJ/mol is in excellent agreement with the low-coverage cryogenic determination of ΔH. The methodology hereby proposed facilitates reliable calculations of the enthalpy of adsorption at room temperatures for weakly-adsorbing gases.

08 HYDROGEN↗
DOE OSTI2025

Experimental observation of nonlinear relation between pressure and water flux is consistent with the solution-diffusion model

In several recent studies, it has been proposed that the fundamental understanding of penetrant transport in dense polymer membranes occurring via the solution-diffusion model, which has been the generally accepted theoretical framework for describing penetrant transport in such materials for the past several decades, is flawed. An alternate mechanistic framework based on the idea of two-phase flow in a porous medium (i.e., pore-flow) has been broadly advanced instead, with proponents of this approach claiming that the pore-flow theoretical framework provides the necessary mechanistic insight to design novel polymeric membrane materials for emerging applications. In this study, we show experimental results for hydraulic permeation of water that are entirely consistent with the solution-diffusion theory, without modification, for three dense polymeric membranes: crosslinked poly(ethylene glycol diacrylate) (XLPEGDA), Nafion 117 ionomer in the sodium counterion form (Nafion 117-Na), and cellulose acetate (CA). By measuring water flux at transmembrane pressures up to 240 bar, we observe a nonlinear relationship between the transmembrane pressure (TMP) and water flux, J w , for XLPEGDA and Nafion 117-Na, while this relationship is linear for CA. We demonstrate that the behavior of these three materials is described via the solution-diffusion model. According to the solution-diffusion model, flux is, to a good approximation, proportional to the transmembrane concentration difference induced by the pressure difference across the membrane, rather than to TMP itself. Water sorption isotherms are reported for all three materials. They further justify the nonlinear relationship between TMP and J w observed in XLPEGDA and Nafion 117-Na, emphasizing that the nonlinearity in the flux/TMP relationship stems from nonlinearities in the sorption isotherm with pressure. Additionally, the relationship between water flux and TMP can be predicted, a priori, with no adjustable parameters when a predictive model for the diffusion coefficient of water is employed in conjunction with the experimental water sorption isotherms in the solution-diffusion model. Furthermore, our results demonstrate the validity of the solution-diffusion model to describe transport of penetrants in dense polymer membranes, while highlighting the sensitivity of the solution-diffusion model to the many physical and mathematical simplifications commonly applied to the theory in literature.

materials↗
DOE OSTI2020

The Effect of η Phase Precipitates on the Creep Behavior of Alloy 263 and Variants

Phase stability is an important design parameter for Ni-based superalloys to be used in future advanced ultra-supercritical (AUSC) power plants as exposure times in this type of environment are considerable. In this investigation, microstructures based on candidate alloy 263 were obtained with varying amounts of η precipitates using isothermal exposure at 800°C for times ranging from 1,000 h to 10,000 h. The effect of η phase stability on the creep properties was determined using creep specimens isothermally aged at 800°C for 8 h, 3,000 h, 5,000 h and 10,000 h prior to creep screening. The creep life was found to exponentially decrease with increasing density of η phase while the elongation to failure was found to increase. Furthermore, the minimum creep rate was related to the density of η phase; a relationship that did not depend on the alloy formulation. Modification to the Ti and Al concentrations slowed down the γ' to η transformation while doubling the γ' fraction after standard heat treatment. By modifying the Ti and Al content, and thereby improving γ' stability over η, the creep lives of specimens isothermally aged for up to 5,000 h were greater than that of the nominal alloy in its standard aged condition.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The formation and evolution of Ni 2 Cr precipitates in Ni–Cr model alloys as a function of stoichiometry characterized by synchrotron x-ray diffraction

Ni-based alloys containing significant amounts of Cr, may precipitate new phases during long-term service at elevated temperatures. In Ni–Cr model alloys, the formation of Ni 2 Cr has been found to impact the material properties, including the strength and ductility after isothermal aging below the critical temperature of 590°C. In this work, we quantify the formation and evolution of long-range ordering in four Ni–Cr binary model alloys (Ni/Cr = 1.8, 2.0, 2.2, 2.4) after isothermal aging up to 10,000h at temperatures between 373 and 475°C. The alloys were characterized by hardness testing and synchrotron-based x-ray diffraction to quantify the impact of Ni 2 Cr phase fraction and precipitate size on mechanical properties. After 500h of isothermal aging, the formation of Ni 2 Cr was detected in all four alloys at 475°C. In the stochiometric alloy samples (Ni/Cr = 2.0), the formation of Ni 2 Cr was found after 500 and 3000h aging at 418 and 373°C, respectively. We found that the matrix lattice contraction and Ni 2 Cr phase fraction both saturate after early aging times. This is in stark contrast to the hardness and Ni 2 Cr precipitate size that both continue to increase with increasing aging time. Our results highlight that changes in hardness correlate linearly with Ni 2 Cr precipitate size rather than phase fraction. This important structure-property relationship can potentially help define Ni–Cr-based component lifetimes directly through an understanding of how Ni 2 Cr formation impacts strength and ductility. Here, we find that a precipitation hardening model for critical resolved shear stress with weakly coupled dislocations shows good agreement with the material property changes quantified from experimental measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Unraveling the multi-step crystallization mechanism of polytetrafluoroethylene, modified polytetrafluoroethylene, and their nanocomposites with boron nitride nanobarbs: Experimental insights and theoretical analysis

The non-isothermal crystallization behavior and kinetics of polytetrafluoroethylene (PTFE) composites with boron nitride nanobarb (BNNB), a new generation nanostructure with unique surface morphology and mechanical “barbs” have been analyzed, understanding these properties is essential for their high-end applications as thermal interface materials (TIM) for microwave, 5G and microelectronic devices. The analysis of the crystallization parameters includes crystallization onset, peak and end temperatures, crystallization half-life and overall crystallinity of PTFE, modified PTFE and their BNNB composites. The results were further analyzed using theoretical models such as the combined Avrami-Ozawa model. It was found that BNNB supports crystallization in the modified PTFE but shows minimal effect on the crystallization of PTFE. Due to the limitation of the classical theoretical models used above in fully characterizing the multi-step crystallization process of PTFE, an in-depth analysis using the model-free advanced isoconversional computation was used to characterize the PTFE crystallization based on the evolution of activation energy with fractional crystallinity and for the first time with temperature. Three kinetic regions were identified in the crystallization mechanism. Here, this study investigated the molecular organization and microstructural evolution of PTFE, modified PTFE and their composites during non-isothermal crystallization using advanced X-ray scattering measurements. An insight into the changes undergone by the material's microstructural units including crystallite size and morphology, lamellar thickness and lamellar interfacial layer thickness, and crystallographic phase dynamics during non-isothermal cooling from the melt, was provided here in this work. The effect of copolymer modification of PTFE and the inclusion of pristine and functionalized BNNB (a thermally conductive and electrically insulating ceramic) are both new investigations that provide valuable knowledge for the development of materials with strong matrix-nanofiller interaction and guidance for optimizing sintering and cooling cycles, two key steps in PTFE processing that largely affect the material microstructural features. Overall, the result of the three-part investigation demonstrates that BNNB supports crystallization in the modified PTFE up to 20 wt% concentration and at low and high cooling rates typically used in the industrial processing of PTFE.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Experimental and Theoretical Thermokinetic Analysis of the Na 3 VO 4 –CO 2 Reaction System under Different Physicochemical Conditions

Sodium vanadate (Na 3 VO 4 ) was synthesized, structural and microstructurally characterized as well as tested as possible carbon dioxide (CO 2 ) captor through thermodynamic calculations in addition to thermogravimetric dynamic and isothermal analyses. Ceramic structural characterization evidenced the formation of the Na 3 VO 4 crystal phase, while microstructural features evidenced dense agglomerates with a poor specific surface area, although synthesis temperature (600 °C) was not as high as that for other sodium ceramics. Na 3 VO 4 was investigated for the CO 2 capture process through dynamic and isothermal experiments, in the presence or absence of oxygen. All these experiments showed that Na 3 VO 4 possesses interesting CO 2 capture properties in a specific temperature range (520 and 600 °C). Isothermal products characterization allowed to elucidate the reaction pathway, implying the Na 4 V 2 O 7 formation, as possible reactive intermediate. In fact, reaction path evolution was supported with theoretical thermodynamic data. Moreover, the kinetic analysis, performed using the Avrami-Erofeev mathematical model and subsequent thermodynamic data obtention, probed the positive influence of the temperature, in both CO 2 superficial and bulk chemical sorption processes, while oxygen addition did not enhance the superficial process. Additionally, as it would be expected, kinetics were enhanced as a function of the CO 2 concentration.

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

Confinement-mediated phase behavior of hydrocarbon fluids: Insights from Monte Carlo simulations

The phase behavior of hydrocarbon fluids confined in porous media has been discussed to deviate significantly from that in the bulk environment due to the existence of sub-10nm pores. Though experiments and simulations have measured the bubble/dew points and sorption isotherms of hydrocarbons confined in both natural and synthetic nanopores, the confinement effects in terms of the strength of fluid-pore interactions tuned by surface wettability and chemistry have received comparably less discussion. More significantly, the underlying physics of confinement-induced phenomena remain obfuscated. In this work, we studied the phase behavior and capillary condensation of n-hexane to understand the effects of confinement at the molecular level. To systematically investigate the pore effects, we constructed two types of wall confinements; one is a structureless virtual wall described by the Steele potential and the other one is an all-atom amorphous silica structure with surface modified by hydroxyl groups. Our numerical results demonstrated the importance of fluid-pore interaction, pore size, and pore morphology effects in mediating the pressure-volume-temperature (PVT) properties of hydrocarbons. The most remarkable finding of this work was that the saturation pressure predicted from the van der Waals-type adsorption isothermal loop could be elevated or suppressed relative to the bulk phase, as illustrated in the graphical abstract. As the surface energy (i.e., fluid-pore interaction) decreased, the isothermal vapor pressure increased, indicating a greater preference for the fluid to exist in the vapor state. Sufficient reduction of the fluid-pore interactions could even elevate the vapor pressure above that of the bulk fluid.

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

Understanding Pore Filling Processes and Adsorption/Desorption Hysteresis in Nanoporous Metal–Organic Frameworks: Insights from Grand Canonical Monte Carlo Simulations and Free Energy Calculations

Grand canonical Monte Carlo (GCMC) simulations were used to investigate pore filling and hysteresis in nanoporous metal-organic frameworks (MOFs). Adsorption and desorption isotherms were calculated for argon at 87 K in 1866 MOFs from the CoRE MOF database and for short n-alkanes in selected MOFs, keeping the adsorbent structure rigid. Analysis of the molecular configurations showed two different mechanisms and origins of hysteresis: one involving a transition of the adsorbate arrangement in the pores similar to a gas-to-liquid transition associated with a large change in the loading and one more similar to a liquid-to-solid transition associated with a relatively small change in the loading. Our GCMC simulations in MOFs with diverse pore topologies indicate exceptions to an empirical relationship for the minimum diameter of a cylindical pore required for hysteresis as a function of the adsorbate diameter and reduced temperature. The simulations reveal some structures where isotherms exhibit two steps in the adsorption branch and only one step in the desorption branch. Hysteresis loops with a different number of adsorption and desorption steps are not common. Here, to better understand why hysteresis is observed in the GCMC simulations, the concept of the transition probability for observing a step in the adsorption isotherm at a given pressure in a GCMC simulation is introduced. We used two different methods to calculate the transition probabilities and find that these yield comparable results. Furthermore, the transition probability provides a measure for the length of GCMC simulations to yield reliable results.

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

Crystallization Mechanism in Spark Plasma Sintered Bulk Metallic Glass Analyzed using Small Angle Neutron Scattering

Understanding the thermal stability of metallic glasses is critical to determining their safe temperatures of service. In this paper, the crystallization mechanism in spark plasma sintered Fe 48 Cr 15 Mo 14 Y 2 C15B6 metallic glass is established by analyzing the crystal size distribution using x-ray diffraction, transmission electron microscopy and in-situ small angle neutron scattering. Isothermal annealing at 700 °C and 725 °C for 100 min resulted in the formation of (Fe,Cr) 23 C 6 crystals, measured from transmission electron micrographs, to be from 10 to 30 nm. The small angle neutron scattering intensity measured in-situ, over a Q-range of 0.02 to 0.3 Å -1 , during isothermal annealing of the sintered samples, confirmed the presence of (Fe,Cr) 23 C 6 crystals. The measured scattering intensity, fitted by the maximum entropy model, over the Q-range of 0.02 to 0.06 Å -1 , revealed that the crystals had radii ranging from 3 to 18 nm. The total volume fraction of crystals were estimated to be 0.13 and 0.22 upon isothermal annealing at 700 °C and 725 °C for 100 min respectively. The mechanism of crystallization in this spark plasma sintered iron based metallic glass was established to be from pre-existing nuclei as confirmed by Avrami exponents of 0.25 ± 0.01 and 0.39 ± 0.01 at the aforesaid temperatures.

42 ENGINEERING↗
DOE OSTI2020

Dynamic behavior of reversible oxygen migration in irradiated-annealed high temperature superconducting wires

Large-scale extrusion-based additive manufacturing of high-performance thermoplastic composites like fiber reinforced polyphenylene sulfide (PPS) is well-suited for tooling applications to lower manufacturing costs and lead times. Autoclave tooling requires good mechanical performance at temperatures even above the glass transition temperature (T g ). In this work, the authors have investigated a post-process isothermal annealing technique to improve the mechanical properties of various grades of carbon fiber reinforced PPS components printed on the Big Area Additive Manufacturing (BAAM) system. Since PPS is a semi-crystalline polymer, crystallinity can change during annealing and affect the mechanical properties of the part. In addition, isothermal annealing can also lead to solid-state structural changes in the form of thermal and oxidative branching and/or crosslinking reactions in some grades of PPS which can alter the crystallization process. This work reports the effect of annealing on dynamic mechanical properties of BAAM printed components, along with studies to determine the effect of annealing on crystallinity and the occurrence of oxidative reactions. Results showed that isothermal annealing at 250 °C for 18 h improved the storage modulus of all selected grades (neat and reinforced) of PPS at temperatures above T g . Annealing led to an overall increase in the degree of crystallinity, with secondary crystallization taking place. Although oxidative structural changes were observed to occur more on the surface of the PPS parts, they primarily influenced the size of crystals formed and did not significantly alter the degree of crystallinity at various regions within the sample.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Quantifying uncertainty in analysis of shockless dynamic compression experiments on platinum. I. Inverse Lagrangian analysis

Absolute measurements of solid-material compressibility by magnetically driven shockless dynamic compression experiments to multi-megabar pressures have the potential to greatly improve the accuracy and precision of pressure calibration standards for use in diamond anvil cell experiments. Here, to this end, we apply characteristics-based inverse Lagrangian analysis (ILA) to 11 sets of ramp-compression data on pure platinum (Pt) metal and then reduce the resulting weighted-mean stress–strain curve to the principal isentrope and room-temperature isotherm using simple models for yield stress and Grüneisen parameter. We introduce several improvements to methods for ILA and quasi-isentrope reduction, the latter including calculation of corrections in wave speed instead of stress and pressure to render results largely independent of initial yield stress while enforcing thermodynamic consistency near zero pressure. More importantly, we quantify in detail the propagation of experimental uncertainty through ILA and model uncertainty through quasi-isentrope reduction, considering all potential sources of error except the electrode and window material models used in ILA. Compared to previous approaches, we find larger uncertainty in longitudinal stress. Monte Carlo analysis demonstrates that uncertainty in the yield-stress model constitutes by far the largest contribution to uncertainty in quasi-isentrope reduction corrections. We present a new room-temperature isotherm for Pt up to 444 GPa, with 1-sigma uncertainty at that pressure of just under ±1.2%; the latter is about a factor of three smaller than uncertainty previously reported for multi-megabar ramp-compression experiments on Pt. The result is well represented by a Vinet-form compression curve with (isothermal) bulk modulus K 0 = 270.3 ± 3.8 GPa, pressure derivative K$^{'}_{0}$= 5.66 ± 0.10, and correlation coefficient $R_{K_{0},K^{'}_{0}}$= –0.843.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Thermal Stability Testing of LANA.75

Savannah River Site has used the metal hydride LaNi 4.25 Al 0.75 (LANA.75) in the Tritium Facilities for over two decades. LANA.75 beds store significant quantities of tritium but have a limited service life due to the radiolytic decay of tritium to 3He within the metal matrix. It has been shown that the isotherm performance of a tritium-aged LANA.75 sample can be restored by heating under vacuum. Additional investigation is needed to ensure there are no unexpected changes to the hydride before this technique is employed in full-scale beds in the Tritium Facilities. In addition, it is necessary to verify the regenerable behavior and thermal stability of LANA.75 on a small scale prior to it being implemented on a large scale. In this work, a non-tritiated bench-scale LANA.75 sample was held at 750°C under vacuum for 200 h to simulate exposure to multiple restoration evolutions. Hydride isotherm performance, chemical composition, crystallinity, and morphology are compared between the pre-restorative and post-restorative testing samples. No significant changes were observed in composition or crystallinity. Comparison of pre-anneal and post-anneal isotherms showed that performance improved rather than deteriorated during the evolution. Scanning electron microscopy analysis showed small growths on the particle surface after exposure to regeneration conditions. Additional testing will be required to determine the cause of these growths.

08 HYDROGEN↗
DOE OSTI2021

Discovery of two bright high-redshift gravitationally lensed quasars revealed by Gaia

We present the discovery and preliminary characterisation of two high-redshift gravitationally lensed quasar systems in Gaia Data Release 2 (DR2). Candidates with multiple close-separation Gaia detections and quasar-like colours in WISE, Pan-STARRS, and DES are selected for follow-up spectroscopy with the New Technology Telescope. We confirm DES J215028.71-465251.3 as a $z$ = 4.130 ± 0.006 asymmetric, doubly imaged lensed quasar system and model the lensing mass distribution as a singular isothermal sphere. The system has an Einstein radius of 1.202 ± 0.005 arcsec and a predicted time delay of ~122.0 d between the quasar images, assuming a lensing galaxy redshift of $z$ = 0.5, making this a priority system for future optical monitoring. We confirm PS J042913.17+142840.9 as a $z$ = 3.866 ± 0.003 four-image quasar system in a cusp configuration, lensed by two foreground galaxies. The system is well modelled using a singular isothermal ellipsoid for the primary lens and a singular isothermal sphere for the secondary lens with Einstein radii 0.704 ± 0.006 and 0.241 ± 0.030 arcsec, respectively. A maximum predicted time delay of 9.6 d is calculated, assuming lensing galaxy redshifts of $z$ = 1.0. Furthermore, PS J042913.17+142840.9 exhibits a large flux ratio anomaly, up to a factor of 2.66 ± 0.37 in i band, that varies across optical and near-infrared wavelengths. We discuss LSST and its implications for future high-redshift lens searches and outline an extension to the search using supervised machine learning techniques.

Astronomy & Astrophysics↗
DOE OSTI2024

Critical behavior in monoclinic Cr 3 ⁢Te 4

High-quality monoclinic Cr 3 ⁢Te 4 crystals are synthesized by the chemical vapor transport method and confirmed by synchrotron powder x-ray diffraction. Critical behaviors of Cr 3 ⁢Te 4 are studied by bulk isothermal magnetization measurements around its paramagnetic to ferromagnetic phase transition. From these data, we evaluate the modified Arrott plot, the Kouvel-Fisher plot, and critical isotherms, and estimate critical exponent values of β~0.3827, γ~1.2119, and δ~4.16, and the Curie temperature T C ≈321 K. After the scaling, the isothermal magnetization curves below and above the critical temperatures collapse into two independent universal branches which signifies the reliability of our critical exponent estimation. Finally, the estimated critical exponent values of Cr 3 ⁢Te 4 do not match well with theoretically calculated values of any three-dimensional models.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Pore characterization of the Marcellus Shale by nitrogen adsorption and prediction of its gas storage capacity

In a shale gas reservoir, pore characterization is an important factor used to determine gas storage capacity. However, the nanometer (nm)-scale pore system in shale is difficult to explore by traditional optical, scanning electron microscopy, or even nuclear magnetic resonance well logging. We have investigated the pore structure and storage capacity of the Marcellus Shale through integration of petrophysical analysis from laboratory and well-logging data and nitrogen adsorption. The isotherm of Marcellus Shale is a composite isotherm, which has features of type I, type II, and type IV isotherms with type H4 of the hysteresis loop, suggesting slit-like pores developed in the Marcellus Shale. Quantitative analysis of pore volumes from the nitrogen adsorption indicates that density porosity may be more properly used to approximate the shale porosity and estimate the shale gas volume. In addition, the specific surface area, micropore, and mesopore volumes have a positive relationship with the kerogen volume and total organic content (TOC). By using the Langmuir and Brunauer-Emmet-Teller models, the simulated result indicates that the higher adsorbed quantity of the Marcellus Shale could be the result of the increase of micropore volume contributed, by the increase of kerogen or TOC content. The proposed equations rapidly compute TOC, a key parameter to predict gas storage capacity in overmature shale such as the Marcellus Shale.

Geochemistry & Geophysics↗
DOE OSTI2023

Temperature-induced densification in compressed basaltic glass revealed by in-situ ultrasonic measurements

Abstract Acoustic velocities of a model basalt glass (64 mol% CaMgSi2O6 + 36 mol% CaAl2Si2O8) were measured along different pressure-temperature (P-T) paths. One set of experiments involved isothermal compression-decompression cycles, performed at temperatures of 300, 641, 823, and 1006 K and pressures up to 12.2 GPa. The other set of experiments involved constant-load heating-cooling cycles at temperatures up to 823 K and pressures up to 7.5 GPa. Both sets of experiments were performed in a multi-anvil apparatus using a synchrotron-based ultrasonic technique. Our results show that the glass compressed isothermally at 300 K (cold-compression) displays anomalously decreasing compressional (VP) and shear (VS) wave velocities with increasing pressure until ~8 GPa. Beyond 8 GPa, both VP and VS start to increase sharply with pressure and irreversible densification of the glass occurred, producing large hysteresis loops of velocities upon decompression. However, for the glass compressed isothermally at increasingly higher temperatures (hot-compression), the velocity minima gradually shift to lower pressures. At temperature close to the glass transition temperature Tg, the velocity minima disappear completely, displaying a monotonic increase of velocities during compression and higher VP and VS during decompression. In addition, constant-load heating-cooling experiments show that velocities generally decrease slightly with increasing temperature, but start to increase once heated above a threshold temperature (~650 K). During cooling the velocities increase almost linearly with decreasing temperature, resulting in higher velocities (~1.5–2.5% higher) when returned to 300 K. This implies that a temperature-induced densification may have occurred in the glass at high pressures. Raman spectra on recovered samples show that the hot-compressed and high-P heated glasses contain distinctly densified and depolymerized structural signatures compared to the initial glass and the cold-compressed glass below the velocity transition pressure PT (~8 GPa). Such densification may be attributed to the breaking of bridging oxygen bonds and compaction in the intermediate-range structure. Our results demonstrate that temperature can facilitate glass densification at high pressures and point out the importance of P-T history in understanding the elastic properties of silicate glasses. Comparison with melt velocity suggests that hot-compressed glasses may better resemble the pressure dependence of velocity of silicate melts than cold-compressed glasses, but still show significantly higher velocities than melts. If the abnormal acoustic behaviors of cold-compressed glasses were used to constrain melt fractions in the mantle low-velocity regions, the melt fractions needed to explain a given velocity reduction would be significantly underestimated at high pressures.

Geochemistry & Geophysics↗