Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shale”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A multiscale packed-bed reactor model for sustainable ethylene production via chemical looping oxidative coupling of methane

The rising global warming concerns and shale gas discovery have prompted research in the direction of greenhouse gas (GHG), such as methane, reduction and conversion. Oxidative coupling of methane (OCM) offers a pathway to low carbon-intense valorization of methane while producing ethylene, a chemical regarded as central to the petrochemical industry. Even after decades of OCM discovery, researchers keep understanding the process and underlying chemical reactions in a pursuit to achieve industrial viability for OCM. Here, in general, OCM suffers from low C 2 selectivity, yield and reactor temperature runaways due to highly exothermic nature of its reactions. Computational Fluid Dynamics (CFD) tools help analyze spatial gradients within the reactor to deeply understand the diffusion of species, mass and heat transfer phenomena. Furthermore, challenges associated with scaling up such as hot spot formation and parametric sensitivity can be addressed without having to expend on costly experiments. The current paper presents a multiscale packed-bed reactor CFD model coupled with a chemical kinetic model for the chemical looping OCM. The CFD model includes two scales i.e., macroscale for catalyst bed and microscale for individual pellets. Moreover, a chemical kinetic model based on 10 gas-phase reactions is integrated with the CFD model. An additional surface reaction for the formation of gas-phase oxygen from catalyst surface is added to account for the absence of feed oxygen. The model is calibrated against experimental results. The calibrated model captures trends in CH 4 conversion, C 2 selectivity and C 2 yield within a ± 4.35 % range across a temperature range of 700-900 °C. Moreover, model fidelity is evaluated by varying key computational parameters such as mesh resolution and time step size. The model is also verified by varying the inlet methane concentration and the gas hourly space velocity (GHSV) and comparing the results with literature. A sensitivity analysis and scale-up of the current model is undergoing.

Chemical looping↗
DOE OSTI2022

Measurement and modeling of methane diffusion in hydrocarbon mixtures

Methane (CH 4 ) dissolution and diffusive mass transfer in liquid hydrocarbon mixtures is of key interest in the context of enhanced oil recovery from tight shales. In this paper, we have studied CH 4 dissolution and diffusion in normal alkane mixtures at a temperature of 50 ⁰C and at pressures of ~8 MPa. For the measurement of CH 4 dissolution/diffusion in bulk liquid hydrocarbon mixtures, we have utilized a high pressure and temperature, constant-volume diffusion (CV-D) setup. Here, we have studied CH 4 diffusion in mixtures with three long-chain normal alkanes: decane (C 10 ), dodecane (C 12 ) and hexadecane (C 16 ). We have measured CH 4 solubility and diffusion in its binary mixtures with each of these three normal alkanes, as well as in ternary and quaternary mixtures. During the experiments, the swelling of the liquid mixture due to CH 4 dissolution was measured in situ via a cathetometer and was subsequently integrated into the data analysis. A key conclusion from this study is that the solubility and transport properties of the multicomponent mixtures can be predicted accurately from binary mixture measurements using an appropriate Equation of State (EOS) and Wilke’s simplification of the classical Maxwell-Stefan (MS) diffusion theory. This observation can facilitate accurate prediction of diffusive mass transfer in more complex liquid hydrocarbon mixtures.

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

Bridging adsorption behavior of confined CH 4 -CO 2 binary mixtures across scales

An accurate understanding of the competitive adsorption of CH 4 -CO 2 binary mixtures in nano-confined systems is critical for engineering CO 2 storage in shale gas reservoirs. Due to difficulties in making reliable experimental observations in nano-scale, atomistic simulations (ASs), such as the Grand Canonical Monte Carlo (GCMC) method, provide a viable approach to studying the adsorption behavior of confined fluids. ASs are, however, limited in the size of the compositional domain due to the high computational cost. This work proposes a framework that combines AS and the lattice Boltzmann (LB) method to bridge the physics of confined fluids across scales. The Peng–Robinson equation of state (PR-EoS) produces fugacity coefficients, which serve as input for conducting multi-component GCMC simulations. These GCMC simulations explore the competitive adsorption behavior of CH 4 -CO 2 in nano-slits at various composition, pressure, and channel-width conditions. Both components generate adsorption layers with high densities near the walls with CO 2 preferentially adsorbing compared to CH 4 on the organic walls of carbon sheets. At the mesoscale, a pseudopotential model represents the intermolecular forces in multi-component, multiple-relaxation-time LB simulations. The LB simulations are in good agreement with the GCMC results, allowing us to obtain values for tunable LB parameters. We then extend the use of LB to simulate adsorption behavior in complex networks with nano-sized channels. The phase behavior and fluid properties in the complex geometries of nano-channels differ from nano-slits and bulk systems. Furthermore, the bridging of physics from GCMC (microscale) to LB (mesoscale) via the macroscale PR-EoS connects the adsorption behavior of binary systems across scales.

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

Experimental and theoretical determinations of hydrogen isotopic equilibrium in the system CH 4 —H 2 —H 2 O from 3 to 200°C

The stable isotopic composition of methane (CH 4 ) is commonly used to fingerprint natural gas origins. Over the past 50 years, there have been numerous proposals that both microbial and thermogenic CH 4 can form in or later attain hydrogen isotopic equilibrium with water (H 2 O) and carbon isotopic equilibrium with carbon dioxide (CO 2 ). Evaluation of such proposals requires knowledge of the equilibrium fractionation factors between CH 4 and H 2 O or CO 2 at the temperatures where microbial and thermogenic CH 4 form in or are found in the environment, which is generally less than 200°C. Experimental determinations of these fractionation factors are only available above 200°C, requiring extrapolation of these results beyond the calibrated range or the use of theoretical calculations at lower temperatures. Here, we provide a calibration of the equilibrium hydrogen isotopic fractionation factor for CH 4 and hydrogen gas (H 2 ) ( D α CH4(g)–H2(g) ) based on experiments using γ-Al 2 O 3 and Ni catalysts from 3 to 200°C. Results were regressed as a 2 nd order polynomial of 1000 × ln D α CH4(g)–H2(g) vs. 1/T (K -1 ) yielding: 1000 × l n D α C H 4 ( g ) - H 2 ( g ) = 3.5317 × 10 7 T 2 + 2.7749 × 10 5 T - 179.48 We combine this calibration with previous experimental determinations of hydrogen isotope equilibrium between H 2 , H 2 O(g), and H 2 O(l) and we provide an interpolatable experimental calibration of 1000 × ln D α CH4(g)–H2O(l) from 3 to 200°C. Our resulting 4th order polynomial is the following equation: 1000 × l n D α C H 4 ( g ) - H 2 O l = - 7.9443 × 10 12 T 4 + 8.7772 × 10 10 T 3 - 3.4973 × 10 8 T 2 + 5.4398 × 10 5 T - 382.05 At 3°C, the value from our calibration differs by 93‰ relative to what would be calculated based on the extrapolation of the only experimental calibration currently available to temperatures below its calibrated range (lowest temperature of 200°C; Horibe and Craig, 1995). We additionally provide new theoretical estimates of hydrogen isotopic equilibrium between CH 4 (g), H 2 (g), and H 2 O(g) and carbon isotopic equilibrium between CH 4 (g) and CO 2 (g) using Path Integral Monte Carlo (PIMC) calculations. Our PIMC calculations for hydrogen isotopic equilibrium between CH 4 and H 2 agree 1:1 with our experiments. Finally, we compile carbon and hydrogen isotopic measurements of CH 4 , CO 2 , and H 2 O from various environmental systems and compare observed differences between carbon and hydrogen isotopes to those expected based on isotopic equilibrium. We find that isotopic compositions of some microbial gases from marine sedimentary, coalbed, and shale environments are consistent with those expected for CH 4 H 2 O(l) hydrogen and CH 4 CO 2 carbon isotopic equilibrium. In contrast, microbial terrestrial and pure culture gases are not consistent with both CH 4 H 2 O(l) hydrogen and CH 4 CO 2 carbon isotopic equilibrium. Overall, these results are explained qualitatively using previously developed conceptual models that link free energy gradients available to microorganisms to the degree that their enzymes can promote isotope-exchange reactions between CH 4 , CO 2 , and H 2 O.

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

Experimental and numerical investigation of fracture conductivity between non-smooth rock surfaces with and without proppant

The enhancement of fracture conductivity is vital for the efficient recovery of subsurface resources, such as geothermal energy and petroleum hydrocarbons. Proppants, granular materials injected into hydraulic fractures to maintain their conductivity, have been studied primarily in the context of smooth fractures (i.e., fractures between smooth rock surfaces). However, non-smooth fractures (i.e., fractures between rough rock surfaces) are common in geoenergy reservoirs and thus require further investigations. In this study, we conducted laboratory measurements of fracture conductivity on shale slabs with non-smooth surfaces and carried out numerical simulation using the lattice Boltzmann (LB) method, which aimed to investigate the conductivity of non-smooth fractures with and without proppants placement. When ceramic proppant with an areal concentration of 2 lb/ft 2 was placed in the fracture, the conductivity was enhanced by roughly 3-8 times compared to fractures without proppant. In fractures with proppant, gas-measured conductivity was higher than that measured with water due to proppant embedment caused by water. The experiments demonstrate the advantages of using proppant in fractures, even if the rock surface roughness can provide certain fracture conductivity via the self-propping mechanism. For fractures without proppants, high rock surface roughness is not necessarily favorable for enhancing fracture conductivity because the self-propping mechanism requires shear slip along the fracture surface. If there is no shear slip, high rock surface roughness can cause a detrimental effect on the fracture conductivity due to the interlocking effect. Utilizing advanced experimental equipment and LB modeling, this research explores the interplays between proppant placement, fracture geometry, and stress conditions to develop a comprehensive understanding of the productivity in non-smooth fractures. Further, the outcomes of this investigation indicate the importance of creating fractures with surface roughness during hydraulic fracturing and will contribute to the development of more efficient stimulation techniques for subsurface energy extraction.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Soft interface instability and gas flow channeling in low-permeability deformable media

Understanding gas percolation through a clay layer or a shale formation is of great importance for the development of a geologic repository for nuclear waste disposal, a subsurface system for gas storage, and an engineering approach for hydrocarbon extraction from unconventional reservoirs. Gas injection experiments have revealed complex dynamic behaviours of gas percolation through water saturated compacted bentonite, characterized by a high breakthrough pressure, rapid breakthrough, a pressure/stress decay after the breakthrough, a relatively high migration rate, high-frequency periodic/nonperiodic variations in flow rate, stepwise rate reductions during relaxation, and low gas saturation over the whole process, all indicating channelling nature of the processes. Using linear stability analyses, we show that this channelling can autonomously emerge from the instability of the deformable interface between the injected gas and the compacted bentonite matrix driven by local stress concentration, pore dilation, and hydrologic gradient. Channel patterns formed would possess a fractal geometry. We further show that, once a percolating channel is established, the gas injected would percolate through the channel in a chain of gas bubbles, also due to the interface instability, resulting in periodic/chaotic variations in gas flow rate. Our work provides a unified explanation for key features observed for gas percolation in low-permeability deformable media. The work also suggests a possibility of designing an engineered barrier system for a nuclear waste repository that can have controllable gas release while limit water transport.

Bentonite↗
DOE OSTI2023

Prediction and demonstration of periodic tensile cracking in rate-dependent porous cement

Periodic dilation bands are shown to occur when cylinder-shaped saturated cement specimens are subjected to rapid decompression, with the spacing between the bands scaling as a power law of the decompression time. This behavior is predicted by theory, which is specified to this scenario based on prior work that was more general and applied to compaction. Specifically, the formation of the bands is shown to coincide with periodicity that naturally arises in the effective stress when the material follows a Terzaghi-type consolidation law, although in reverse for dilation and with material deformation described by a rate-dependent viscoplastic law. When strain rate is proportional to effective stress to a power that is greater than one, periodic regions of tensile effective stress arise through the poromechanical fluid-solid coupling. The theoretically-predicted exponent of the resulting power-law relationship between dilation band spacing and unloading time successfully brackets the experimental results, for which it is observed that the exponent is slightly stronger than a square-root relationship (0.28 to 0.67) at testing temperatures of 20 °C and 90 °C. Further, these results demonstrate that rapid-depressurization leads to periodic fracturing that could, on the one hand, be detrimental to the isolation provided by cement used to seal wellbores in the petroleum industry. On the other hand, the dilation bands could also be favorable to production if generated in low-permeability reservoir rocks such as shales that are targeted for petroleum production or granites that are targeted for geothermal energy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Evaluation of the economic implications of varied pressure drawdown strategies generated using a real-time, rapid predictive, multi-fidelity model for unconventional oil and gas wells

Experience has suggested that pressure maintenance in hydraulically fractured reservoirs via lower, more sustained production drawdowns may offer improved cumulative recovery and overall resource extraction efficiency compared to more rapid drawdown approaches aimed at generating high initial production. However, given the inherent variability of oil and natural gas markets, operators pursue production strategies that maximize profitability over resource extraction efficiency. This study focuses on evaluating the implications of contrasting pressure drawdown strategies on the long-term production and resulting economics for a real, producing unconventional gas well in the Marcellus Shale of the Appalachian Basin using a techno-economic analysis approach. Our research combines elements of well-specific horizontal well design, production forecasting, equipment sizing and capital cost estimation, operating cost estimation, and revenue and tax calculations. Gas production forecast outlook scenarios were generated under varying pressure drawdowns using two approaches: 1) a novel physics-informed machine learning workflow and 2) traditional reservoir simulation. A discounted cash flow model was used to evaluate the resulting economic implications for each drawdown scenario—generating output for exploring the coupled effect of factors like the timing and volume of gas production, prevailing economic and market conditions for natural gas, and overall estimated ultimate recovery on profitability metrics such as internal rate of return and net present value. Results show that there is potential to maximize the cumulative gas produced in the specific case study well by employing a lower pressure drawdown. Conversely, the greatest profitability is achieved using rapid drawdown as signified by a small, specific subset of our outlook scenarios. On an averaging basis, we find that the combinations of highest cumulative producing and most profitable scenarios occur under lower drawdowns with long (>40 years) producing timeframes, but require higher relative gas price and lower discounting considerations. Further, the machine learning predictive outlooking capability proved effective for enabling rapid generation of a multitude of scenario forecasts. As a result, a variety of prominent example cases could be generated to strike the balance of greater productivity and economic return given their associated producing features and economic conditions when compared to similar producing scenarios—critical insight that offers improved decision support for unconventional oil and gas operations.

42 ENGINEERING↗
DOE OSTI2025

A review on advances in oxidative coupling of methane (OCM) for industrial use and prospects of CO 2 –H 2 O splitting integration

The discovery of shale gas reserves has encouraged the development of direct methods for methane conversion into valuable chemicals, offering an alternative to indirect approaches that involve an energy-intensive and intermittent syngas production step, leading to high CO 2 emissions. Amongst the direct methods, the oxidative coupling of methane (OCM) is a potential pathway to reduce CO 2 emissions and can produce commodity chemicals such as ethylene, a chemical regarded as central to the petrochemical industry. Even though OCM has been studied for over four decades, the technology still has not found commercial application. Amongst the challenges regarding industrial deployment of OCM, the most significant one is the requirement of a high ethylene yield of 30 % which is currently reported to be around 20 %. Moreover, the highly exothermic nature of the process and controlling the carbon selectivity over oxides of carbon (COx) is the heart of the problem. Numerous researchers have presented promising results in terms of catalysts, reactor designs and feeding strategies for OCM. However, due to lack of inclusiveness in the results, none of the combination of catalysts, reactors and system optimizations has been able to bring about its industrial viability. The current paper presents an extensive review of the noteworthy attempts to achieve industrial targets for OCM. Moreover, a comprehensive criteria is presented which highlights the desired end state for the industrial deployment of OCM technology. Furthermore, the criteria is based on literature survey and a comparison with industrially deployed ethylene production plants i.e., naphtha or ethane steam cracker plants. Finally, a novel integration technology is presented which includes a combination of OCM and CO 2 -H 2 O splitting in a chemical looping reactor design to enable efficient energy utilization and minimal heat losses to the environment.

CO2 Splitting↗
DOE OSTI2026

Coupling geophysical, geological, geochemical and mineralogical assessments to examine preferential contaminant transport pathways in interbedded fractured bedrock

This study shows that a multi-faceted approach, combining borehole geophysical logging and surface seismic P-wave first-arrival tomography with confirmatory coring, well installation, and chemical and mineralogical analysis, is effective for identifying difficult-to-locate preferential contaminant transport pathways in deeper fractured bedrock. Seismic tomography detected porous 10–20 m wide elongated fractured conduits that allow acidic groundwater contaminated with uranium (U) and nitrate (NO 3 - ) to migrate within interbedded shale-limestone bedrock over 1000 m from a former disposal facility (S-3 Ponds site) located at the DOE Y-12 National Security Complex in Tennessee (USA). Conventional drilling techniques would easily miss these conduits because they are oriented parallel with fractured bedding planes. Synchrotron analysis of aquifer solids revealed that > 95 % of the U is hexavalent (U VI ). This uranyl (UO 2 2+ ) species is coordinated with carbonate, iron oxide, silicate and phosphate minerals within cemented white to yellow precipitates, which contain U concentrations as high as ∼21.6 % by weight fraction. Identifying the presence of these mineral phases, enables a further understanding of the potential effectiveness of remediation actions. The combination of methodologies presented here can also be applied to other explorations, such as the detection of water supply.

Fate and transport↗
DOE OSTI2020

Towards real-time forecasting of natural gas production by harnessing graph theory for stochastic discrete fracture networks

In this work, we compare hydrocarbon production curves obtained from a graph-based reduced-order model with the high-fidelity Discrete Fracture Network (DFN) predictions for a fracture network created using data from a real shale site. We observe that the bounds for the high fidelity DFN model lie within the bounds for the reduced order model, implying that the reduced-order model provides a conservative estimate. Moreover, we found that except for first-passage times and late arriving mass, the production curves from the reduced-order model predict transport accurately. However, it is to be noted that the results are inspite of trading a three-dimensional geometry for a reduced system in the form of a graph, one that is 500–1000 times faster in terms of computational efficiency (for this particular application). In addition, we also compare the production curves for large drawdown and small drawdown using our graph approach. The reduced-order model is successful in showing that the long term productivity is higher in case of small drawdown although the initial productivity is higher for large drawdown. Thus, this reduced-order model offers great potential in uncertainty quantification for production, as well as in providing operators with information to make real-time decisions for optimal production.

03 NATURAL GAS↗
DOE OSTI2021

Characterizing mineralization on low carbon steel exposed to aerated and degassed synthetic hydraulic fracture fluids

This study investigates the mineralization occurring at the interface between injected hydraulic fracturing fluids and steel production casing during the shut-in period of unconventional shale wells. Mild steel coupons were reacted with synthetic hydraulic fracturing fluids with varying base fluid and dissolved gas concentrations. Experiments were carried out in an autoclave system set at 50 °C and 2000 PSI for 2 and 14 days. Physical and mineralogical changes to the steel coupons’ surface were characterized using an FEI Quanta 600 field emission gun (FEG) scanning electron microscope (SEM) and a PANalytical XPert Pro diffractometer. Barite precipitation occurred ubiquitously during the first 48 hours in all experimental scenarios. Sulfate radicals created from degradation and reduction of ammonium persulfate breakers facilitate barite formation. Concomitant barite and corrosion related precipitates on the sample surfaces suggest iron hydroxy salts act as nucleation sites for mineral scale growth. These results suggest hydraulic fracturing fluids may have deleterious effects on the wellbore before production of the well begins.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Critical mineral source potential from oil & gas produced waters in the United States

The volume of produced water, a by-product of oil & gas operations and other energy processes, has been growing across the United States (U.S.) along with the need to manage or recycle this wastewater. Produced water contains many naturally occurring elements of varying concentrations, including critical minerals which are essential to the clean energy transition. However, the current understanding of critical mineral concentrations in produced water and the associated volumes across the U.S. is limited. This study has assessed available databases and literature to gain insight into the presence and concentration of five high priority critical minerals, namely cobalt, lithium, magnesium, manganese, and nickel. The U.S. Geological Survey's National Produced Waters Geochemical Database was the main data source used for determining average critical mineral concentrations in produced water from the major oil and gas reservoirs in the U.S. The volumes of produced water for these major reservoirs were coupled with these concentrations to provide insights into where critical minerals are likely to have high abundance and therefore more recovery options. The analysis indicated the highest recovery potential for lithium and magnesium from produced water in the Permian basin and the Marcellus shale region. However, these assessments should be considered conservative due to the limited availability of reliable concentration data. Finally, it is expected more critical mineral recovery options could emerge with comprehensive characterization data from more recent and representative sources of produced water.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Comparison of direct and CO 2 -oxidative dehydrogenation of propane

Propane dehydrogenation to propylene has received increasing attention due to the fast growth in propylene demand and the exploration of shale gas containing propane. Direct dehydrogenation of propane (DDHP) offers high propylene selectivity but is limited by quick deactivation due to coke formation. CO 2 -assisted oxidative dehydrogenation of propane (CO 2 -ODHP) can consume CO 2 and meanwhile reduce coke deposition via the reverse Boudouard reaction. Here in the current review, direct and CO 2 -assisted dehydrogenation of propane has been compared from three aspects: reaction mechanisms, catalyst compositions, and CO 2 footprint analysis. As the average CO 2 emissions from electricity generation decrease due to the utilization of renewable energy, CO 2 -ODHP potentially leads to a net-negative CO 2 footprint while DDHP cannot.

58 GEOSCIENCES↗
DOE OSTI2025

Autocorrelation and Multifractal Detrended Fluctuation Analyses Reveal Superdiffusive Mass Transport in Solvent-Filled Nanoporous Media

Fluorescence fluctuation spectroscopy experiments were conducted to better understand the complex mass transport dynamics of organic molecules in liquid-filled nanoporous media. Anodic aluminum oxide (AAO) membranes incorporating 10 and 20 nm diameter cylindrical pores were employed as model materials. Nile red (NR) dye was used as a fluorescent tracer. The dye was dissolved separately in ethanol and toluene at a concentration of 20 nM and used to fill the membrane nanopores. Confocal fluorescence microscopy was employed to capture photon intensity time series data reflecting apparent diffusion of the dye within the pores. Autocorrelation of these data revealed that NR diffusion within the membranes occurred over a broad range of time scales. The autocorrelation decays were fit to a model for one-dimensional diffusion incorporating both fast and slow components having apparent diffusion coefficients, D f and D s , differing by a factor of ∼100. The fast mechanism was attributed to hindered bulk-like diffusion in the central pore cavity, while slow diffusion likely involved absorption of the dye to the pore surfaces. Unfortunately, important evidence of diffusion anomalies is lost in the broad autocorrelation decays obtained. The method of multifractal detrended fluctuation analysis (MF-DFA) was applied to the same data as a means to overcome this limitation. MF-DFA revealed that time series acquired from within the nanopores were multifractal and exhibited evidence of anomalous superdiffusion, likely resulting from the participation of a desorption-mediated diffusion mechanism. Monte Carlo simulations of time series modeling desorption-mediated diffusion in cylindrical nanopores provided support for this assignment. Here, the new knowledge gained affords an improved understanding of hydrocarbon dynamics within nanoporous oil and gas shales.

Diffusion↗
DOE OSTI2022

Autothermal Chemical Looping Oxidative Dehydrogenation of Ethane: Redox Catalyst Performance, Longevity, and Process Analysis

Energy-efficient upgrading of stranded ethane from shale gas to olefins holds the promise of increasing the supply of useful chemical feedstocks while reducing flaring-based CO 2 emissions. Previously, we reported a modular ethane-to-liquids (M-ETL) system based on a chemical looping oxidative dehydrogenation (CL-ODH) scheme. In this article, we present long-term (>1200 h) results of Li 2 CO 3 -promoted La 0.8 Sr 0.2 FeO 3 corresponding to ~4125 CL-ODH cycles in a large laboratory-scale packed bed reactor. Temperature monitoring along the bed length confirmed the exothermicity of both the oxidative dehydrogenation and regeneration steps, enabling an autothermal operation. Product gas analysis indicated that the redox catalyst maintains a high C 2+ selectivity (~90%) and ethane conversion (~67%) at 735 °C after continuous cycling of >1000 h. As a result, product distributions and heats of reactions were used to update our M-ETL process model for revised techno-economic analysis, demonstrating that the current system is economically viable with relatively low required selling prices across a wide range of operating scenarios.

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

Hierarchical Galloaluminosilicate MFI Catalysts for Ethane Nonoxidative Dehydroaromatization

This study presented an approach to overcome diffusion limitation and imbalance between Brønsted and Lewis acidities that affect conversion and aromatic yield in alkane nonoxidative dehydroaromatization (DHA) conversion. Three different Ga-promoted ZSM-5 catalysts were synthesized for the ethane DHA reaction. The hierarchical galloaluminosilicate Meso-GaZSM-5 was prepared by the solid-state crystallization technique wherein Ga was introduced through the in situ synthesis process. The microporous galloaluminosilicate Micro-GaZSM-5 was synthesized by the hydrothermal technique wherein Ga was introduced in situ during zeolite formation. The hierarchically structured Ga-Meso-ZSM-5 was synthesized by the solid-state crystallization method, but the Ga was introduced by the incipient wetness method. Additional Pt was incorporated by the incipient wetness technique for all three catalysts to enhance the DHA activity. The solid-state crystallization technique created the hierarchical structure without using a template. With the in situ technique, Ga species replaced the Al framework, reducing the strong Brønsted acidity. Meanwhile, they created the extraframework GaO + , introducing the strong Lewis acid sites into catalysts. Due to better synergistic interaction between the oxidation state of metals, acidic properties, and catalyst porosity, Pt-Meso-GaZSM-5 showed higher C 2 H 6 conversion and aromatic selectivity. The improvement on diffusivity not only reduces the effect of channel blockage but also improves the utilization of active sites. As a result, this work sheds light on the direct, nonoxidative conversion of shale gas without going through the syngas route.

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

Multiscale Modeling of a Direct Nonoxidative Methane Dehydroaromatization Reactor with a Validated Model for Catalyst Deactivation

Due to the recent boom in shale gas production, aromatics production using direct nonoxidative methane dehydroaromatization (DHA) is being investigated extensively. However, due to rapid coke formation, catalysts in the nonoxidative methane DHA reactors get deactivated, which is one of the critical issues for the commercial success of the methane DHA process. In this paper, a model for catalyst deactivation is developed. Rate models for other DHA reactions are developed by considering the decrease in the catalyst activity with time. Due to the very fast coke formation rate on the fresh catalyst, there is coke formation immediately upon the introduction of the feed. Therefore, an algorithm is developed for estimation of the initial state of the reactor and the kinetic parameters by coupling an iterative direct substitution approach with an optimization approach. Transient experimental data from an in-house reactor are first reconciled and then used for developing the kinetic model including the coke formation model. Using the rate model, a dynamic, heterogeneous, multiscale reactor model with embedded heating is developed. Here, the model couples the catalyst pellet level model with a reactor level model. Impacts of temperature, L/D ratio, and scheduling of reactors on variability in conversion and yield with time are studied.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗