Search NASA⌕ Search

SEARCH · Search NASA

Results for “Enthalpy method”

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 145 records · Page 8

Analysis and prediction of reaction kinetics using the degree of rate control

“Degree of rate control” (DRC) analysis provides a quantitative approach for analysing the kinetics of multi-step reaction mechanisms that has been widely applied to both heterogeneous and homogeneous catalysis research, as well as electrocatalysis. The DRC of any given transition state or intermediate is defined as a partial derivative such that it approximately equals the fractional increase in net rate to the product of interest per differential decrease in its standard-state free energy for that species (÷RT), holding constant the standard-state free energies of all other transition states and intermediates. Even very complex mechanisms usually have only a few species with non-zero DRCs and are thus the species whose interactions with the catalyst most strongly affect the net rate. These key DRC values thus offer a simple and intuitive route to optimize catalyst materials, especially with the assistance of computational methods like density functional theory (DFT). These high-DRC species are also the species whose energetics must be most accurately measured or calculated to achieve an accurate kinetic model for any reaction mechanism. In simple cases with a single “rate-determining step”, the DRC for its transition state (TS) is + 1. Catalyst-bound intermediates, on the other hand, often have negative DRCs equal to a small integer times their fractional occupancy of catalyst sites. The apparent activation energy equals a weighted average of the standard-state enthalpies (relative to reactants) of all the species (intermediates, transition states and products) in the reaction mechanism, each weighted by its DRC (+RT). It has been shown that the apparent transfer coefficient in electrocatalysis, an inverted form of the Tafel slope, is a weighted average of the number of electrons transferred to generate each intermediate or product species in the mechanism, weighted again by the DRC. Quantitative analysis of kinetic isotope effects (KIEs, or the ratio of net rates for different reactant isotopes) in complex mechanisms has shown that the logarithm of the KIE equals the weighted average over all species in the mechanism of the difference between the two isotopes in their standard-state free energies (÷RT), again weighted by the DRC. The reaction orders with respect to fluid-phase concentrations of reactants, products and intermediates have also been proven to be directly related to DRCs. Thus, there are numerous experimental observables which equate to short linear combinations of DRCs, so that combinations of experimental measurements might provide access to DRC values. Since its invention for transition states in 1994 and its generalization to include intermediates in 2009, the DRC has thus far mainly been calculated for microkinetic models based either entirely on DFT or on DFT with the key energies (i.e., those for high-DRC species) being fine-tuned to match experiments. The relationships summarized in this work provide new opportunities for using experiments earlier in the development and optimization of microkinetic models that require input from computational catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Combined First-Principles and Experimental Investigation into the Reactivity of Codeposited Chromium–Carbon under Pressure

High-pressure synthesis in the diamond anvil cell suffers from the lack of a general approach for the control of precursor stoichiometry and homogeneity. Here, we present results from a new method we have developed that uses magnetron cosputtering to prepare stoichiometrically precise and atomically mixed amorphous films of Cr:C. Laser-heated diamond anvil cell experiments carried out on a flake of this sample at pressures between 13.5 and 24.3 GPa lead to the observation of Cr 3 C (Pnma) over the entire pressure range–in good agreement with our in-house theoretical predictions–but also reveal two other metastable phases that were not expected: a novel monoclinic chromium carbide phase and the NaCl-type CrC (Fm3̅m) phase. The unexpected stability of CrC is investigated by using first-principles methods, revealing a large stabilizing effect tied to substoichiometry at the carbon site. These results offer an important case study into the current limitations of crystal structure prediction methods with regard to phase complexity and bolster the growing need for advanced theoretical approaches that can more completely survey experimentally unexplored phase space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Two dimensional topology optimization of heat exchangers with the volume fraction method

We perform a comparison study of two topology optimizations methods applied to the design of a two fluids heat exchanger modeled with a coupled thermal-flow problem. The flow follows an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation without internal generation and high Peclet number. To keep the two fluid phases separated, we solve two Stokes-Brinkman equations, where the Brinkman term models the other phase as a solid. These two velocity fields are then fed to the heat transfer equation. Our goal is to maximize the enthalpy at the cold outlet while constraining the pressure drop. We first solve the design modeling the solid and fluid phases with a volume fraction variable. A SIMP-like penalization in the Brinkman term drives the optimization to a discrete design. The cost and constraint function derivatives are calculated with the library pyadjoint and the optimization is performed by IPOPT. We present optimized designs in two dimensions and discuss the influence of the parameters.

Beck, VictorA.↗

Fundamentals of SrCoO 3-δ Based Oxygen Deficient Perovskites as Cathodes for Solid Oxide Fuel Cells

The local structure and oxygen stoichiometry in oxide systems have a profound impact on oxygen electrocatalysis encountered in metal-air batteries and solid oxide fuel cells (SOFCs). However, this knowledge is often obtained under conditions different from the real working conditions of the material, resulting in misinterpretation and misunderstanding. This PhD dissertation aims to obtain the structure and oxygen-stoichiometry information of a class of perovskite oxides under their real working conditions in solid oxide fuel cells. Several perovskite oxides were selected for the study: Sr 0.9 Y 0.1 CoO 3-δ (SYC10), Sr 0.9 Y 0.3 CoO 3-δ (SYC30), SrCo 0.9 Nb 0.1 O 3-δ (SCN10) and SrCo 0.9 Ta0.1O 3-δ (SCT10). The local crystal structure and oxygen stoichiometry of these materials were systematically characterized with in-situ neutron diffraction (ND). The oxygen stoichiometry was also measured by thermogravimetric analysis (TGA) and iodometric titration methods. To establish the correlation among structure, oxygen stoichiometry and oxygen reduction reaction (ORR) activity, high temperature electrochemical impedance spectroscopy, oxygen permeation and single cell testing were also performed. In chapter 2, SYC10 is described to have a more symmetrical structure and $V_o^{··}$ distribution, and higher $V_o^{··}$ concentration than SYC30. Molecular orbital energy analysis based on the local structure of the ORR-active Co1-polyhedra indicates SYC30 has a higher Fermi level relative to O-2p energy level in the active Co1-polyhedra, and thus a higher motional enthalpy for $V_o^{··}$ migration. However, the ORR activity of Y-doped SrCoO 3-δ (SCO) are found uncompetitive with other popular catalysts. In chapter 3 and 4, systematic structural studies are presented for Nb and Ta-doped SCO, namely SCN10 and SCT10, respectively. The methodology of unrevealing the structure-activity relationships in chapter 2 is also found applicable to SYC10, SCN10 and SCT10 as presented in chapter 5. The electrochemical tests indicate that SCT10 is the best catalyst due to its excellent thermal and electrochemical stability. In chapter 6, SCT10 is evaluated in a real solid oxide fuel cell in the form of nanoscaled SCT10 layer coating on the commercial (La 0.6 Sr 0.4 ) 0.95 Co 0.2 Fe 0.8 O 3-δ (LSCF). The bilayer structured cathode shows excellent properties of coarsening-resistant, Sr-segregation-free and high ORR-activity.

Yang, Tianrang↗

Assessing entropy for catalytic processes at complex reactive interfaces

When chemical reactions are accelerated by a catalyst, entropy differences between reactants and their transient intermediates can be the driving force behind the promotion or inhibition of desired and parasitic chemical pathways. Understanding and controlling catalytic processes therefore requires both a fundamental and practicable understanding of entropy in addition to enthalpy. In unstructured media such as the vapor phase equilibrated with sparsely covered surfaces, entropy can be adequately accounted for by well-established approaches based on translational, rotational, and harmonic vibrational partition functions. However, these approximations become inadequate in more complex condensed phase environments, e.g., solid liquid interfaces of confined reaction spaces. In this chapter, we provide an overview of the state-of-art in the computational quantification of entropy and its known ramifications on catalysis. The fundamental roles of thermodynamics and kinetics in catalysis are covered in enough detail to appreciate and contextualize the computational methods employed to compute chemically accurate estimates of entropy. These methods are discussed in appropriate detail and range from the ubiquitous harmonic oscillator approximation where entropy unrelated to high frequency oscillations is typically underestimated, to enhanced free energy sampling with molecular dynamics where the desired accuracy must be weighed against the associated computational cost of obtaining it. The rising importance of machine learning and artificial intelligence in accelerating methodological progress in this field is touched upon, as well. Finally, applications, successes, and pitfalls of using these methods are provided to showcase past and present accomplishments while clarifying where improvements in both understanding and methodology are still needed.

Kollias, Loukas↗

Electronic states and transitions of PrO and PrO + probed by threshold ionization spectroscopy and spin–orbit multiconfiguration perturbation theory

The precise ionization energy of praseodymium oxide (PrO) seeded in supersonic molecular beams is measured with mass-analyzed threshold ionization (MATI) spectroscopy. Here, a total of 33 spin–orbit (SO) states of PrO and 23 SO states of PrO + are predicted by second-order multiconfigurational quasi-degenerate perturbation (MCQDPT2) theory. Electronic transitions from four low-energy SO levels of the neutral molecule to the ground state of the singly charged cation are identified by combining the MATI spectroscopic measurements with the MCQDPT2 calculations. The precise ionization energy is used to reassess the ionization energies and the reaction enthalpies of the Pr + O → PrO + + e – chemi-ionization reaction reported in the literature. An empirical formula that uses atomic electronic parameters is proposed to predict the ionization energies of lanthanide monoxides, and the empirical calculations match well with available precise experimental measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High accuracy barrier heights, enthalpies, and rate coefficients for chemical reactions

Abstract Quantitative chemical reaction data, including activation energies and reaction rates, are crucial for developing detailed kinetic mechanisms and accurately predicting reaction outcomes. However, such data are often difficult to find, and high-quality datasets are especially rare. Here, we use CCSD(T)-F12a/cc-pVDZ-F12// ω B97X-D3/def2-TZVP to obtain high-quality single point calculations for nearly 22,000 unique stable species and transition states. We report the results from these quantum chemistry calculations and extract the barrier heights and reaction enthalpies to create a kinetics dataset of nearly 12,000 gas-phase reactions. These reactions involve H, C, N, and O, contain up to seven heavy atoms, and have cleaned atom-mapped SMILES. Our higher-accuracy coupled-cluster barrier heights differ significantly (RMSE of ∼5 kcal mol −1 ) relative to those calculated at ω B97X-D3/def2-TZVP. We also report accurate transition state theory rate coefficients $${k}_{\infty }(T)$$ k ∞ ( T ) between 300 K and 2000 K and the corresponding Arrhenius parameters for a subset of rigid reactions. We believe this data will accelerate development of automated and reliable methods for quantitative reaction prediction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Inorganic Halide Double Perovskites with Optoelectronic Properties Modulated by Sublattice Mixing

All-inorganic halide double perovskites have emerged as a promising class of materials that are potentially more stable and less toxic than lead-containing hybrid organic-inorganic perovskite optoelectronic materials. In this work, 311 cesium chloride double perovskites (Cs 2 BB'Cl 6 ) were selected from a set of 903 compounds as likely being stable on the basis of a statistically learned tolerance factor (t) for perovskite stability. First-principles calculations on these 311 double perovskites were then performed to assess their stability and identify candidates with band gaps appropriate for optoelectronic applications. We predict that 261 of the 311 Cs 2 BB'Cl 6 compounds are likely synthesizable on the basis of a thermodynamic analysis of their decomposition to competing compounds (decomposition enthalpy <0.05 eV/atom). Of these 261 likely synthesizable compounds, 47 contain no toxic elements and have direct or nearly direct (within 100 meV) band gaps between 1 and 3 eV, as computed with hybrid density functional theory (HSE06). Within this set, we identify the triple-alkali perovskites Cs 2 [Alk] + [TM] 3+ Cl 6 , where Alk is a group 1 alkali cation and TM is a transition-metal cation, as a class of Cs 2 BB'Cl 6 double perovskites with remarkable optical properties, including large and tunable exciton binding energies as computed by the GW-Bethe-Salpeter equation (GW-BSE) method. We attribute the unusual electronic structure of these compounds to the mixing of the Alk-Cl and TM-Cl sublattices, leading to materials with small band gaps, large exciton binding energies, and absorption spectra that are strongly influenced by the identity of the transition metal. The role of the double-perovskite structure in enabling these unique properties is probed through an analysis of the electronic structures and chemical bonding of these compounds in comparison with other transition-metal and alkali transition-metal halides.

14 SOLAR ENERGY↗

An ab initio molecular dynamics investigation of the thermophysical properties of molten NaCl-MgCl 2

Molten salts have many applications in the nuclear and solar energy industries for thermal storage and heat transfer applications. However, there is a knowledge gap in molten salt thermophysical properties which hinders the technical readiness level of molten salt applications, especially in the nuclear industry. A common method of investigating new materials is through ab initio Molecular Dynamics (AIMD) simulations which is an effective tool to investigate structural and thermophysical properties at realistic temperatures. NaCl-MgCl 2 is an inexpensive salt that is a good candidate for use as a heat transfer medium in solar power applications or in the secondary loop of a nuclear reactor. In this article, the thermophysical properties of NaCl-MgCl 2 are computed via AIMD calculations to supplement the limited experimental studies in the literature. Here a wide range of compositions and temperatures for the pseudo-binary NaCl-MgCl 2 were used to calculate the density, heat capacity, compressibility, enthalpy of mixing, and volumetric thermal expansion coefficient. AIMD is shown to accurately model the densities of molten NaCl-MgCl 2 as there is good agreement with the available literature. This work observed a transition to a monotonic increase of the density with respect to MgCl 2 composition occurring above 1100 K. The heat capacity values increase uniformly with respect to concentration of MgCl 2 at a rate of 2.85 J/mol-K per 10 mol% of MgCl 2 . Select thermophysical properties are fit to a Redlich-Kister expansion for utilization in multiphysics simulations.

36 MATERIALS SCIENCE↗

A review of cladding failure thresholds in RIA conditions based on transient reactor test data and the need for continued testing

Transient reactor experiments on light water reactor (LWR) fuel pins had been conducted since the beginning of the nuclear era to help determine core coolability and cladding failure thresholds. During one such test in November of 1993 at the CABRI transient test reactor on a test involving a high burnup fuel rod with a corroded Zircaloy-4 cladding it was first observed that cladding failures could occur prior to a departure from nucleate boiling (pre-DNB) at lower-than-expected peak radial average enthalpies. Thirteen additional tests would be performed in the CABRI reactor over the next decade on fuel rods with higher burnups [1]. Additionally, a larger testing program at the NSRR reactor in Japan with high burnup fuel would uncover a similar trend of pre-DNB ruptures in high burnup test rods at lower-than-expected peak enthalpies [2]. Numerous out of pile testing programs involving a variety of innovative mechanical testing techniques have been employed in an attempt to better quantify the failure thresholds of corroded zirconium alloy cladding in these rapid heating and loading conditions [3][4][5]. While previously interim guidance had been issued, in June of 2020 the NRC officially published updated regulatory guidance to account of these pre-DNB failures in Regulatory Guide 1.236 [6]. This paper will present an independent review of the publicly available transient reactor test database on higher burnup LWR pins conducted at the CABRI and NSRR reactors as well as review of a selection of published out of pile mechanical testing methods. The purpose of the review is to determine how well the new regulatory limits are supported by experimental data. The review will identify if additional transient reactor tests could provide additional support for the NRC guidance or identify the need for revisions. The evaluation will consider how far the existing database can be extrapolated when considering low hydrogen zirconium alloy claddings (with and without protective coatings) containing very high burnup (> 70 MWd/kgU) UO2 fuel pellets. Finally, the authors will suggest how out of pile mechanical tests can be used in conjunction with a limited number of transient reactor tests to develop cladding specific failure thresholds in RIA type transients.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

An Implementation of Replica Exchange with Dynamical Scaling for Efficient Large-Scale Simulations

An implementation of the replica exchange with dynamical scaling (REDS) method in the commonly used molecular dynamics program GROMACS is presented. REDS is a replica exchange method that requires fewer replicas than conventional replica exchange while still providing data over a range of temperatures and can be used in either constant volume or constant pressure ensembles. Details for running REDS simulations are given, and an application to the human islet amyloid polypeptide (hIAPP) 11-25 fragment shows that the model efficiently samples conformational space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molten Flux Synthesis of Plutonium (IV) Silicates K 2 PuSi 2 O 7 and Rb 2 PuSi 6 O 15

Few plutonium containing silicates have been reported in the literature and only one other tetravalent plutonium silicate structure has been reported using the flux growth method. With the goal of understanding actinide incorporation into silicates, we describe herein two new plutonium (IV) silicates, Rb 2 PuSi 6 O 15 and K 2 PuSi 2 O 7 , synthesized using a mixed alkali chloride/alkali fluoride flux growth. Both structures crystallize in the C2/c space group with lattice parameters a=24.2980(10) Å, b=7.1466(3) Å, c=17.2025(8) Å and $\blacksquare$96.6560(10)° for Rb 2 PuSi 6 O 15 , and a=9.9478(2) Å, b=5.59820 (10) Å, c=13.1848(2) Å, and $\blacksquare$105.7400(10)° for K 2 PuSi 2 O 7 . Comparisons to the synthesis of other reported alkali actinide silicates (Ce, U, and Th) are made. Raman spectra and DFT calculations of formation enthalpies and vibrational modes provide confirmation of structure refinements and enhanced product characterization for comparison with other materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Model for the Solid–Liquid Interfacial Free Energy at High Pressures

The free energy involved in the formation of an interface between two phases (e.g., a solid–liquid interface) is referred to as the interfacial free energy. For the case of solidification, the interfacial free energy dictates the height of the energy barrier required to nucleate stable clusters of the newly forming solid phase and is essential for producing an accurate solidification kinetics model using classical nucleation theory (CNT)-based methods. While various methods have been proposed for modeling the interfacial free energy for solid–liquid interfaces in prior literature, many of these formulations involve making restrictive assumptions or approximations, such as the system being at or near equilibrium (i.e., the system temperature is approximately equal to the melt temperature) or that the system is at pressures close to atmospheric. However, these approximations and assumptions may break down in highly non-equilibrium situations, such as in dynamic-compression experiments where metastable liquids that are undercooled by hundreds of kelvin or over-pressurized by several gigapascals or more are formed before eventually solidifying. We derive a solid–liquid interfacial free-energy model for such high-pressure conditions by considering the enthalpies of interactions between pairs of atoms or molecules. We also consider the contribution of interface roughness (disordering) by incorporating a multilayer interface model known as the Temkin n-layer model. Our formulation is applicable to a diverse variety of materials, and we demonstrate it by developing models specifically for two different materials: water and gallium. We apply our interfacial free-energy formulation to CNT-based kinetics simulations of several suites of dynamic-compression experiments that cause liquid water to solidify to the high-pressure solid polymorph ice VII and have found good agreement to the observed kinetics with only minor empirical fitting.

36 MATERIALS SCIENCE↗

Phase equilibria of advanced technology uranium silicide-based nuclear fuel

The phases in uranium-silicide binary system were evaluated in regards to their stabilities, phase boundaries, crystal structures, and phase transitions. The results from this study were used in combination with a well assessed literature to optimize the U-Si phase diagram using the CALPHAD method. A thermodynamic database was developed, which could be used to guide nuclear fuel fabrication, could be incorporated into other nuclear fuel thermodynamic databases, or could be used to generate data required by fuel performance codes to model fuel behavior in normal or off-normal reactor operations. The U 3 Si 2 and U 3 Si 5 phases were modeled using the Compound Energy Formalism model with 3 sublattices to account for the variation in composition. The crystal structure used for the USi phase was the tetragonal with an I4/mmm space. Above 450°C, the U 3 Si 5 phase was modeled. The composition of the USi 2 phase was adjusted to USi 1.84 . The calculated invariant reactions and the enthalpy of formation for the stoichiometric phases were in agreement with experimental data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Electronic properties of ThC − and ThC

The present study investigates the properties of low-lying states of ThC − and ThC using correlated wave function theories. To this end, we employed multireference calculations and various coupled-cluster approaches in combination with large correlation-consistent basis sets. These methods were applied to examine potential energy curves (PECs), electron configurations, energetics, spectroscopic constants, and spin–orbit coupling effects for 9 states of ThC − and 18 states of ThC. The ground states of ThC − and ThC were identified as single-reference 2 Σ + 1/2 (I) and $^3Σ^+_{0^-} (\textrm{I})$, respectively. Electron detachment from the 7s orbital of the Th center in ThC − [ 2 Σ + (I); 1σ 2 2σ1 3 σ2 1 π 4 ] yields ThC [ 3 Σ + (I); 1σ 2 2σ 1 3σ 1 1π 4 ]. The calculated adiabatic detachment energy (ADE) and vertical detachment energy (VDE) for this process are 1.591 and 1.604 eV, respectively. Furthermore, the dissociation energy (D 0 ) of ThC [$^3Σ^+_{0^-} (\textrm{I})$] is predicted to be 5.099 eV. The standard enthalpy of formation, $ΔH^{°}_{\textrm{f}}$ (298 K), of ThC is estimated to be 822.52 ± 6 kJ mol −1 .

74 ATOMIC AND MOLECULAR PHYSICS↗

Thermodynamic properties of Ca–Pb electrodes determined by electromotive force measurements

Thermodynamic properties of Ca–Pb alloys are investigated by electromotive force (emf) measurements to determine equilibrium cell potentials and phase properties for their application in energy storage systems such as liquid metal batteries. Using the electrochemical cell Ca(s) | CaF 2 (s) | Ca(in Pb) at 700–1060 K, cell emf is measured for thirteen Ca–Pb alloys at mole fractions, x Ca = 0.06–0.80. At 873 K, the equilibrium potentials of liquid Ca–Pb alloys are 0.57–0.62 V versus Ca and the activity values are as low as a Ca = 6.2 × 10 -8 at x Ca = 0.06. In addition, the emf values as a function of temperature provide partial molar quantities (entropy and enthalpy) as well as phase transitions which are corroborated by determining transition temperatures and phase constituents using differential scanning calorimetry (DSC) and powder X-ray diffraction (XRD). This study establishes the fundamental data necessary for the design of Pb-containing liquid metal electrodes through the integration of electrochemical, thermal, and structural properties of Ca–Pb electrodes.

25 ENERGY STORAGE↗

Designing the Platinum Catalyst Layer for Improved Performance and Durability in Anion Exchange Membrane Water Electrolysis

To lower the cost of hydrogen produced by anion exchange membrane water electrolysis (AEMWE), it is critical to reduce the use of platinum group metal (PGM) catalysts within the device. While iridium has been successfully replaced with PGM-free catalysts at the anode, platinum-based (Pt) cathode catalysts are still required to meet the activity and durability targets. This study investigates the impact of commercial Pt/C catalyst loading, ionomer type and content, and electrode fabrication method on the cathode catalyst layer properties and AEMWE performance with the aim of determining the feasibility of reduced Pt loadings. While increased Pt loading is found to improve beginning-of-life performance, the effects are minimal above 0.6 mg/cm 2 . Ink characterization shows that ionomer type and content affect the ink stability, particle size, and percent of unbound ionomer, which further impact the homogeneity of the sprayed catalyst layers. The 5% PiperION cathode exhibited the highest performance, which may be attributed to a balance between the small particle size and the low proportion of unbound ionomer, minimizing kinetic and transport losses. Theoretical calculations show that the ionomers interact differently with the Pt surface, causing different surface charges and water adsorption strength and activating different mechanisms for hydrogen evolution. Pt-PiperION lowered the enthalpy of water-splitting by 0.1 eV compared to Pt alone and allowed for equal site access between adsorbed H* and OH* (both adsorbed at circa −2.2 eV). Although catalyst-coated membrane (CCM) fabrication techniques are desirable for scale-up, no performance enhancement is observed compared with the catalyst-coated substrate approach. Durability, as measured by degradation rates, Pt loss, and catalyst layer restructuring, was found to improve with increased Pt loadings, higher ionomer content, and CCM architectures. These findings provide important insight into the significant role of the cathode in AEMWE and strategies for maintaining the performance with low Pt loading or PGM-free catalysts.

08 HYDROGEN↗

Reexamining supercritical gas adsorption theories in nano-porous shales under geological conditions

Adsorption is an important phenomenon in surface chemistry, especially for nano-porous shales. In the shale gas-in-place, adsorbed gas could contribute up to 85%. However, adsorption is hard to quantitatively characterize due to different adsorption mechanisms, patterns, surfaces, and pore sizes. Moreover, key thermodynamic parameters, such as the enthalpy of adsorption, are challenging to determine due to uncertainties in adsorbed gas densities used in constructing absolute isotherms. In this paper, we revisit the Brunauer, Emmett, and Teller (BET) model, analytically simplify Ono-Kondo (OK) models for subsurface shales, and compare commonly used mono- and multilayer adsorption models (e.g., Langmuir, supercritical Dubinin-Radushkevich (SDR), supercritical BET (SBET), and simplified OK (OKs) models) with recently proposed pressure-dependent adsorption densities to develop a practical and reliable methodology that can be used in the supercritical state, typical for subsurface black shale conditions.Three independent data sets were used for nitrogen and methane adsorption isotherms at different temperatures. Furthermore, we demonstrate that adsorption predicted by the SDR model is comparable or lower than that predicted by the SBET model, but higher than the amounts predicted by Langmuir and OK models. The nitrogen BET method tends to underestimate the accessible SSA for methane. Measurement of isosteric heat of adsorption is suggested instead of the experimental fitting method, due to the significant difference between the calculated results by two commonly used methods. Experimental fitting and simulation methods are also briefly reviewed to guide future research on shale gas adsorption.

58 GEOSCIENCES↗