Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gas phase”

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

Mass transport in nanoarray monolithic catalysts: An experimental-theory study

Reducing the mass transfer resistance globally of a catalyst is a key to enhancing the catalytic reaction kinetics and fully utilizing the catalyst activity. Despite the success in tailoring the external mass transfer in the widely studied washcoat monoliths, the internal mass transfer resistance is difficult to be reduced due to the requirement of increasing macroporosity while maintaining high specific surface area and mechanical stability. Therefore, nanostructured array-based monolithic catalysts (nanoarray catalysts) have been developed in the past decade as a promising class of structured catalysts that may complement or substitute washcoat catalysts. This work fundamentally elucidates the enhanced mass transport properties of the nanoarray monolithic catalysts by a combination of experimental measurements and theoretical modeling. Using a low-dimensional model, the relative contributions of resistances were quantified in terms of chemical kinetics, internal and external mass transfers based on a probe model of C2H4 oxidation over the TiO2 supported Pt-based monolithic catalysts. The nanoarray catalysts displayed a lower internal mass transfer resistance than the washcoat counterparts as a result of the high macroporosity and small thickness of nanoarray layers. Finally, the nanoarray configuration provides a new pathway towards designing high-performance monolithic reactors and catalysts with low internal diffusion limitations for various gas phase reactions.

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

Process intensification of CO 2 capture by low-aqueous solvent

Low-energy solvent-based CO 2 absorption processes have drawn attention as a next-generation post-combustion CO 2 capture technology to reduce CO2 emissions from fossil fuel– or biomass-fired power generation and industrial flue gas streams. A low-aqueous (or water-lean) solvent process may substantially reduce the thermal energy consumption for solvent regeneration. Low-aqueous solvent–based processes are thermally sensitive, requiring a delicate temperature control within the absorber because of the fast exothermic amine-CO 2 reaction and low heat capacity organic diluent. This reaction may result in heat accumulation in a packed absorption column and undesirable CO 2 desorption occurring as the solvent moves through the column, reducing the solvent’s CO 2 capture efficiency if its temperature is not controlled. Using a 3D printed intensified packing device, enhanced heat and mass transfer were demonstrated in an amine-CO 2 scrubbing process using low-aqueous solvent. The multifunctional intensified device facilitates contact of the reactive solvent and gas phases in a single stage and heat removal by a cooling fluid flowing through channels in the interior of the corrugated plates of the device. These functionalities led to effective thermal management along the column via intrastage cooling and significant improvement in CO 2 uptake under a wide range of operating conditions. Intrastage cooling effectively reduced the solvent average temperature along the column by ~10 °C and, as a result, the solvent’s capture efficiency improved by up to 25%.

42 ENGINEERING↗
DOE OSTI2023

Numerical modeling of chemical looping oxidative dehydrogenation of ethane in parallel packed beds

Chemical looping oxidative dehydrogenation (CL-ODH) of ethane has the potential to be a highly efficient alternative to steam cracking for ethylene production. Accurate reactor modeling is of critical importance to efficiently scale up and optimize this new technology. This study reports a one-dimensional, heterogeneous packed bed model to simulate the CL-ODH of ethane to ethylene with a Na 2 MoO 4 -promoted CaTi 0.1 Mn 0.9 O 3 redox catalyst. Here, the overall reaction kinetics was well-described by coupling the gas-phase steam cracking of ethane with the reduction kinetics of the redox catalyst by H 2 and C 2 H 4 . The impact of H 2 on the formation rate of CO 2 byproduct from C 2 H 4 conversion was also thoroughly investigated to validate the applicability of the kinetic model under operational environments. The temperature variation within the different CL-ODH steps and the temperature distribution along the bed were also carefully considered. The accuracy of the model was validated by experiments conducted in a large lab-scale packed bed reactor (200 g catalyst loading), with an average deviation of 2.8% in terms of ethane conversion and ethylene yield. The model was subsequently used to optimize the operating parameters of the CL-ODH reactor, indicating that up to 63.7% single-pass C 2 +olefin yield can be achieved with the current redox catalyst bed whereas further optimization of the redox catalyst to inhibit C 2 H 4 activation can result in 69.4% single-pass C 2 +yield while maintaining low CO 2 selectivity.

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

Theoretical and kinetic modeling study of H 2 S pyrolysis

Hydrogen sulfide pyrolysis was investigated theoretically and through chemical kinetic modeling. Reactions on the SHH potential energy surface, primarily S + H 2 (+Ar) ⇌ H 2 S (+ Ar) (R1) and S + H 2 ⇌ SH + H (R6b) were characterized by ab initio calculations. Results for k 1 were in good agreement with experiment, but deviated strongly below 2000 K from values previously used in modeling. Collider efficiencies for H 2 S, S 2 , and N 2 compared to Ar were calculated for R1. Hydrogen sulfide decomposition experiments reported in literature were re-examined in terms of an updated detailed chemical kinetic model. Concentration profiles for the atomic S at high temperature in shock tubes supported the present value of k 1 and served to constrain the rate constants for reaction of S with SH and H 2 S. To explain results from batch and flow reactors, conducted at high H 2 S concentrations in the 900–1400 K range, a very fast rate constant was required for HSS + H ⇌ SH + SH. Under dilute conditions, the gas-phase chemistry was too slow to compete and the decomposition of H 2 S was controlled by loss on the reactor surface.

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

Statistical distributions for transient transport

Here, this paper introduces the use of statistical distributions based on transport differential equations for clear distinction of transport modes within transient kinetic experiments. More specifically, novel techniques are developed for the transient data obtained through the Temporal Analysis of Products (TAP) reactor and are applicable to experiments where pulse response into a gas flow is used. The methodology allows distinguishing between two domains of diffusion transport in heterogeneous catalytic systems, i.e., Knudsen and non-Knudsen diffusion, using statistical fingerprints, and finding the transition domain. Two distribution parameters were obtained that directly result in coefficients that correspond to the concentration and the rate of transport. Using a linear relationship between the rate and concentration coefficients, Knudsen diffusion is revealed when the rate of transport is constant and non-Knudsen diffusion is confirmed when the rate of transport coefficient is a function of the concentration coefficient. As a result, accurate transport information can be extracted from experimental data even in the presence of comprising instrument drift or noise particularly when analyzing higher pressure pulse responses with complex transport. This enables more direct investigation of experiments influenced by gas-phase reactions.

TAP reactor↗
DOE OSTI2024

PFAS remediation: Evaluating the infrared spectra of complex gaseous mixtures to determine the efficacy of thermal decomposition of PFAS

Due to their widespread production and known environmental contamination, the need for the detection and remediation of per- and polyfluoroalkyl substances (PFAS) has grown quickly. While destructive thermal treatment of PFAS at low temperatures (e.g., 200 to 500oC) is of interest due to lower energy and infrastructure requirements, the range of possible degradation products remains underexplored. To better understand the low temperature decomposition of PFAS species, we have coupled gas-phase infrared spectroscopy with a multivariate curve resolution (MCR) analysis and a database of high-resolution PFAS infrared reference spectra to detect and quantify a complex mixture resulting from potassium perfluorooctanesulfonate (PFOS-K) decomposition. Nine prevalent decomposition products (namely smaller perfluorocarbon species) are identified and quantified.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Three-dimensional thermo-mechanical simulations of heterogeneous solid propellants

Here in this work, we present a numerical framework that describes thermo-mechanical deformations in a burning heterogeneous solid propellant. These deformations are quasi-static at time scales associated with combustion, and the resulting thermo-mechanical formulation is discretized on a Cartesian grid using a hypoelastic law. We use a weak form of Chorin-type projection method to deal with large difference in shear modulus of the constituent materials. Extending our previous two-dimensional work, grid convergence studies for a three-dimensional propellant configuration are presented for the stress, velocity, and reference map components. Finally, simulations are carried out for a random propellant pack that is coupled to a gas phase, and we present results for the pack undergoing combustion, with and without deformations.

simulations↗
DOE OSTI2024

Mesoscale modeling and semi-analytical approach for the microstructure-aware effective thermal conductivity of porous polygranular materials

Here we established a comprehensive modeling approach for investigating the microstructure-aware effective thermal conductivity ($κ_{eff}$) for porous microstructures containing solid particles and gaseous pores. Our approach combines the mesoscale computational modeling framework and the semi-analytical method, allowing for efficient prediction of $κ_{eff}$ for realistic porous microstructures, while considering complicated microstructural thermal conduction pathways effectively in the prediction. We used the diffuse-interface mesoscale computational model to generate extensive simulated $κ_{eff}$ data for realistic digital representations of microstructures with wide ranges of porosity ($f_p$), thermal conductivity of the gas phase ($κ_g$), and thermal conductivity of the solid phase ($κ_s$). From the simulated data, we identified two property variation regimes for $κ_{eff}$: (1) a slow $κ_{eff}$ increase for $κ_s ~ κ_g$; and (2) a faster $κ_{eff}$ increase for $κ_s \gg κ_g$. To capture the key features of the relationship between the microstructure and $κ_{eff}$, we derived a semi-analytical model by introducing structure and intensification factors. The two new factors incorporate the calibrated effective contribution of the solid volume with $κ_s$ and additional interfacial effects into the prediction of $κ_{eff}$, respectively, allowing for consideration of parallel, serial, and interfacial conduction mechanisms effectively. Using the selected simulation data, we quantified key model parameters within the semi-analytical model and verified that the parameterized model exhibits excellent agreement with simulated $κ_{eff}$ for the entire range of the parameter space.

36 MATERIALS SCIENCE↗
DOE OSTI2021

TUMME: Tsinghua University Minnesota Master Equation program

We report that TUMME is a program for assembling and solving master equations for gas-phase chemical kinetics based on chemically significant eigenmodes. TUMME has interfaces to the Gaussian, Polyrate, and/or MSTor output files that allow the master equation code to obtain the microcanonical flux coefficients needed for the coefficient matrix of the master equation. The flux coefficients for reactions with barriers can be calculated by multi-structural variational transition state theory with small-curvature tunneling (MS-VTST/SCT) or by simpler approximations to this such as conventional transition state theory without tunneling (also called RRKM theory). The flux coefficients for barrierless reactions are provided by a hard-sphere model. TUMME is written in double precision with Python 3; quadruple and octuple precision are also available for some subtasks in C++. The Python code can run in serial or parallel (MP or MPI), and the C++ code can run on a single processor or on multiple processors with OpenMP.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

MSTor 2023 : A new version of the computer code for multistructural torsional anharmonicity, now with automatic torsional identification using redundant internal coordinates

The MSTor program is a computer program for calculating partition functions and thermodynamic functions of complex gas-phase molecules with multiple torsions; the multi-structural approximation with torsional anharmonicity (MS-T) can be used based on either a coupled torsional potential or an uncoupled torsional potential. The program can also carry out calculations in the multiple-structure local harmonic (MS-LH) approximation or multi-structural local quasiharmonic (MS-LQ) approximation or by the dual-level MS-T method. Furthermore, the program package includes eight utility codes that can be used as stand-alone programs to help the user to generate the input files for the code or to generate comparison results. The 2023 version of MSTor includes a new capability, namely automatic identification of torsional modes with redundant internal coordinates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

QMMM 2023: A program for combined quantum mechanical and molecular mechanical modeling and simulations

Combined quantum mechanical and molecular mechanical (QM/MM) methods play an important role in multiscale modeling and simulations. QMMM 2023 is a general-purpose program for single-point calculations, geometry optimizations, transition-state optimizations, and molecular dynamics (MD) at the QM/MM level. It calls a QM package and an MM package to perform the required single-level calculations and combines them into a QM/MM energy by a variety of schemes. QMMM 2023 supports GAMESS-US, Gaussian, and ORCA as QM packages and Tinker as the MM package. Four types of treatments are available for embedding the QM subsystem in the MM environment: mechanical embedding with gas-phase calculations of the QM region, electronic embedding that allows polarization of the QM region by the MM environment, polarizable embedding for mutual polarization of the QM and MM regions, and flexible embedding for both mutual polarization and partial charge transfer between the QM and MM regions. Boundaries between QM and MM regions that pass through covalent bonds can be treated by several methods, including the redistributed charge (RC) scheme, redistributed charge and dipole (RCD) scheme, balanced-RC scheme, balanced-RCD scheme, screened charge scheme that takes account of charge penetration effects, and smeared charge scheme that delocalizes the MM charges near the QM–MM boundary. Geometry optimization can be done using the optimizer implemented in QMMM 2023 or the Berny optimizer in Gaussian through external calls to Gaussian. Molecular dynamics simulations can be performed at the pure-MM level, pure-QM level, fixed-partitioning QM/MM level, and adaptive-partitioning QM/MM level. As a result, the adaptive-partitioning treatments permit on-the-fly relocation of the QM–MM boundary by dynamically reclassifying atoms or groups into the QM or MM subsystems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Electronic structure of reduced CeO 2 (111) surfaces interacting with hydrogen as revealed through electron energy loss spectroscopy in comparison with theoretical investigations

Based on both, ab-initio cluster calculations as well as periodic density functional based loss function calculations, we have assigned the origins of the valence electron excitation regime in electron energy loss spectra of well-ordered ceria films in (111) orientation at various states of reduction as well as after exposing the reduced films to hydrogen from the gas phase (Li et al, Angew. Chem. Int. Ed., 58 (2019) 14686–14693). The explicit calculation of intensity distributions using the dipole approximation allow us to draw conclusions about the nature of oxygen vacancies, which occur upon reduction and how those interact with hydrogen to form hydride species. Herein, we also reports a brief discussion of vibrational excitations of hydrogen loaded ceria and corroborates previous interpretations, based on inelastic neutron scattering.

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

Deep root activity overprints weathering of petrogenic organic carbon in shale

The oxidation of organic carbon in sedimentary bedrock (petrogenic OC, OC petro ) is increasingly recognized as a potential source of CO 2 to the atmosphere. Recent studies provide evidence for the mobilization and oxidation of OC petro in sedimentary bedrock during rock weathering. However, the mechanisms and rates remain uncertain, particularly where overlying soils and vegetation drive contemporaneous oxidation of recently fixed organic carbon. Here, in this study, we quantify OC petro weathering across a 16 m shale depth profile in a steep, rapidly eroding forested hillslope in the Northern California Coast Ranges. We report solid and gas phase radiocarbon and stable isotope analyses of samples extracted from specialized in-situ samplers, and a supporting laboratory incubation experiment of the shale regolith. OC petro is removed from the weathered bedrock at a rate of approximately 0.12 gC/m 3 yr, which is orders of magnitude lower than the rate of OC petro oxidation we achieved in the laboratory with crushed samples (557.1 gC/m 3 yr). This disparity occurs despite high O 2(g) content across the depth profile, indicating that physical accessibility of OC petro can regulate oxidative weathering. There is no direct radiocarbon evidence of OC petro oxidation in CO 2(g) across the upper 13 m of the weathering profile during both wet and dry seasons. Instead, vadose zone CO 2(g) production at the site is dominated by respiration of recently fixed carbon associated with deep rooting. OC petro is clearly mobilized across the vadose zone during weathering in this rapidly eroding, oxygen-rich, biologically dynamic hillslope, but at rates far below what can be measured given the contribution of root-derived CO 2(g) .

58 GEOSCIENCES↗
DOE OSTI2025

Zigzag flow reactor for weekly thermochemical energy storage

This paper describes theoretical models and experimental performance of a novel Zigzag Flow Reactor (ZFR) for weekly thermochemical energy storage. The ZFR reduces redox-active metal oxide (MO x ) particles at high temperature (up to ~1100 °C) under inert gas sweep. A physical model demonstrates the approach to process equilibrium by minimizing the associated exergy destruction in a finite number of reaction steps, establishing the thermodynamic requirements for a practical reactor. The model results show several cost-relevant parameter tradeoffs, and the tradeoff analysis implies a cost-optimized set of boundary conditions. Numerical models and prototypes show that the ZFR enables significant gas phase homogenization while simultaneously enabling a customizable MO x residence time in the reactor, both key requirements for approaching an equilibrium process. A scaling model demonstrates the simplicity and affordability of sizing the ZFR to grid-scale levels, with fabrication costs at least five times lower than previously proposed scalable reactor concepts. As a result, a laboratory ZFR prototype achieved an energy storage density of ~90 Wh/kg with CaAl 0.2 Mn 0.8 O 3-δ as the MO x , at temperatures of ~850 °C in >10 h of total runtime.

Thermochemical energy storage↗
DOE OSTI2020

Assessment of various tar and soot treatment methods and a priori analysis of the steady laminar flamelet model for use in coal combustion simulation

Soot formed in the fuel-rich regions of coal flames is known to have a significant impact on the radiative heat transfer characteristics of the flame. It is therefore important to account for the effect of soot radiation when simulating coal combustion systems. Several recent coal combustion studies have used acetylene as a surrogate for devolatilized tar and neglect soot formation altogether, while others have utilized empirical models for tar and soot or even chemical mechanisms which include tar and soot as species. In this work, we investigate several methods of treating coal-derived tar and soot and the impact each treatment has on predicted furnace temperature, particle mass loss, and other quantities of interest. Furthermore, we perform a priori comparisons of simulation results to calculations from a steady laminar flamelet model in order to assess the accuracy of a flamelet approach for each tar and soot treatment considered. Data used as a basis for the study is generated through One-Dimensional Turbulence (ODT) simulation of a coal combustion system in which high- fidelity models are used for gas phase chemistry and devolatilization, where char combustion is neglected.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Tracking natural CO 2 migration through a sandstone aquifer using Sr, U and C isotopes: Chimayó, New Mexico, USA

The geochemical and isotopic characteristics of groundwaters in Chimayó, New Mexico, reflect processes that affect water quality in the Tesuque Aquifer, which overlies a leaking natural CO 2 source in a structurally complex region. In this study, select isotopes (δ 13 C, 87 Sr/ 86 Sr, 234 U/ 238 U) are applied to groundwaters to better understand CO 2 transport mechanisms and related water-rock interactions that impact Chimayó groundwater. Carbon stable isotope ratios of dissolved inorganic carbon (DIC; δ 13 CDIC = -15.10‰ to 4.50‰) identify a distinct source of upward-migrating CO 2 that interacts with the groundwater. Additionally, groundwater 87 Sr/ 86 Sr compositions (0.7098 to 0.7154) reflect intrusion of varying amounts of saline water associated with the high CO 2 source, while 234 U/ 238 U ratios corroborate presence of deep groundwater and suggest impacts from distinctive natural uranium sources are affecting groundwater. Previous work proposed two CO 2 transport mechanisms at the site: (1) dissolved in deep brine that underlies the aquifer and (2) in the gas phase; this study uses isotope mixing models to identify wells that are affected by these CO 2 transport mechanisms and demonstrates that both transport mechanisms are associated with impaired groundwaters. Overall, this study demonstrates that applying multiple isotope systems (δ 13 CDIC, 87 Sr/ 86 Sr, 234 U/ 238 U) is a dynamic tool for identifying and measuring the impact of CO 2 leakage from a sequestration site.

58 GEOSCIENCES↗
DOE OSTI2020

Temperature statistics in a radiatively heated particle-laden turbulent square duct flow

Radiation absorption by preferentially concentrated particles in a turbulent square duct flow is studied experimentally. The particle-laden flow is exposed to near-infrared radiation, and the gas phase temperature statistics are measured along the wall bisector of the duct. It is found that the instantaneous temperature fluctuations are comparable to the overall mean temperature rise. The temperature statistics at the duct centerline and near the wall are qualitatively different. The former reflects preferential concentration in isotropic flows while the latter displays evidence of particle clustering into streamwise elongated streaks. Comparison of the experimental data to a simplified heat transfer model suggests that the Lagrangian evolution of particle clusters and voids, and turbulent mixing in the vicinity of particle clusters, are important. This work was motivated by particle solar receiver technology, but the findings are also relevant to systems where there is localized heat release or mass transfer from disperse particles or droplets. It shows that obtaining Lagrangian histories of particle trajectories is an important next step towards understanding thermal transport phenomena in particle-laden turbulent flows.

42 ENGINEERING↗
DOE OSTI2021

Mathematical modeling of novel porous transport layer architectures for proton exchange membrane electrolysis cells

Thin foil based porous transport layers (PTLs) that contain highly structured pore arrays have shown promise as anode PTLs in proton exchange membrane electrolysis cells. These novel PTLs, fabricated with advanced manufacturing techniques, produce thin, tunable, multifunctional layers with reduced flow and interfacial resistances and high thermal and electric conductivities. To further optimize their design, it is important to understand their fundamental impact on the transport of protons, electrons, and liquid/vapor mixtures in the electrode. In this work, we develop a two-dimensional multiphysics model to simulate the coupled electrochemistry and multiphase transport in an electrolysis cell operated with the novel PTL architecture. The results show that larger pores improve access of water to the anode catalyst layer, which is beneficial for both the oxygen evolution reaction and membrane hydration. Larger pore sizes also improve oxygen gas transport from the catalyst layer, because generated oxygen gas is forced to travel in-plane through the anode catalyst layer until it reaches a pore opening that is connected to a channel. The discussed results confirm that the proposed thin foil based PTLs are fundamentally different from conventional PTLs, such as felts or layered meshes. The model developed in this work also provides generalizable insight into fundamental PEMEC phenomena, such as the competition between liquid and gas phase transport, membrane hydration and water management, and nonuniform electrochemical reactions, which are processes relevant to all PEMEC designs.

25 ENERGY STORAGE↗