SEARCH · Engineering Papers
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.
The topics covered include the following: (1) CO2 laser kinetics modeling; (2) gas lifetimes in pulsed CO2 lasers; (3) frequency chirp and laser pulse spectral analysis; (4) LAWS A' Design Study; and (5) discharge circuit components for LAWS. The appendices include LAWS Memos, computer modeling of pulsed CO2 lasers for lidar applications, discharge circuit considerations for pulsed CO2 lidars, and presentation made at the Code RC Review.
A 3-dimensional kinetic model has been developed to study the dynamics of the storm time ring current in a dipole magnetic field. In this paper, the ring current model is extended to include a realistic, time-varying magnetic field model. The magnetic field is expressed as the cross product of the gradients of two Euler potentials and the bounce-averaged particle drifts are calculated in the Euler potential coordinates. A dipolarization event is modeled by collapsing a tail-like magnetosphere to a dipole-like configuration. Our model is able to simulate the sudden enhancements in the ring current ion fluxes and the corresponding ionospheric precipitation during the substorm expansion.
Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].
δ-Valerolactone (DVL) is a five-carbon (C5) cyclic ester that can undergo ring-opening polymerization to yield high-performance, biocompatible polyesters. But current market prices of C5 chemicals like DVL are very high due to poor availability of C5 feedstock in petroleum. Herein, we demonstrate a novel route to DVL synthesis via dehydrogenation of biomass-derived 2-hydroxytetrahydropyran (HTHP) over Cu/SiO 2 without the use of toxic reagents. Since HTHP exists in thermal equilibrium with 3,4-dihydropyran (DHP) via dehydration, and with 2,2’-oxybis(tetrahydropyran) and 5-(tetrahydropyran-2-yloxy)pentanal via acetalization, we have also determined the thermochemistry (ΔH rxn and ΔG rxn ) of each competing reaction using density functional theory (DFT) calculations at the M06-2X/cc-pVTZ level. Finally, by developing a kinetic model of all 8 reactions involved, we have achieved 84% selectivity to DVL at 150°C in a packed bed reactor for over 72 hours of time on stream.
Anthropogenic perturbations to the nitrogen cycle, primarily through use of synthetic fertilizers, is driving an unprecedented increase in the emission of nitrous oxide (N 2 O), a potent greenhouse gas and an ozone depleting substance, causing urgency in identifying the sources and sinks of N 2 O. Microbial denitrification is a primary contributor to biotic production of N 2 O in anoxic regions of soil, marine systems, and wastewater treatment facilities. Here, through comprehensive genome analysis, we show that pathway partitioning is a ubiquitous mechanism of complete denitrification within microbial communities. We have investigated mechanisms and consequences of process partitioning of denitrification through detailed physiological characterization and kinetic modeling of a synthetic community of Rhodanobacter thiooxydans FW510-R12 and Acidovorax sp. GW101-3H11. We have discovered that these two bacterial isolates, from a heavily nitrate (NO 3 − ) contaminated superfund site, complete denitrification through the exchange of nitrite (NO 2 − ) and nitric oxide (NO). The process partitioning of denitrification and other processes, including amino acid metabolism, contribute to increased cooperativity within this denitrifying community. We demonstrate that certain contexts, such as high NO 3 − , cause unbalanced growth of community members, due to differences in their substrate utilization kinetics. The altered growth characteristics of community members drives accumulation of toxic NO 2 − , which disrupts denitrification causing N 2 O off gassing.
Millions of tons of waste polyester plastics, including polyethylene terephthalate (PET), polytrimethylene terephthalate (PTT), and polybutylene terephthalate (PBT), end up in the environment as soil and water contaminants. Recent advances in enzymatic polyester degradation have motivated researchers toward biorecycling of plastic wastes. Leaf-branch compost cutinase (LCC) enzymes have been proven to be effective in the biodegradation of PET. This study focuses on enzymatic depolymerization of waste PET, PTT, and PBT materials by using an ICCG variant of LCC (LCC ICCG ) produced from Escherichia coli BL21(DE3). The degradation efficiency of the polyesters was determined by the monomer terephthalic acid (TPA) released from the depolymerization reaction. It was found that the most efficient depolymerization was achieved for PET, followed by PTT and PBT. A kinetic model based on Langmuir adsorption and the Michaelis-Menten equation was developed to describe the enzymatic depolymerization of PET, PTT, and PBT with various enzyme and substrate loadings. The model simulation results revealed that the LCCICCG enzyme loading should be linearly increased as the work capacity of the polyester substrate increases. A specific enzyme loading of 0.91 mg/g PET is suggested to achieve 90% depolymerization of PET within three days. Furthermore, the experimental data and model simulation results can be used to help further engineer the enzyme and process to achieve a complete biodegradation of polyester wastes at large scale.
Robust synthesis of vertically aligned carbon nanotubes (VACNT) at large scale is required to accelerate deployment of numerous cutting-edge devices to emerging commercial applications. To address this need, this work demonstrates that the structural characteristics of single-walled CNTs (SWCNTs) produced in a growth regime dominated by bulk diffusion of the gaseous carbon precursor are remarkably invariant over a broad range of process conditions (precursor concentration increased up to 30-fold, catalyst substrate area from 1 cm 2 to 180 cm 2 , growth pressure from 20 to 790 mbar, and gas flowrates up to 8-fold). A simple growth kinetics model is derived that accounts for both reactant bulk diffusion and competing byproduct formation in the presence of excess hydrogen, which quantitatively reproduces the experimental data in the entire isothermal parameter space. The model enables critical predictions for process scale-up optimization including a potential 6-fold increase in CNT production rate, a ~90% carbon conversion efficiency in select conditions, and elimination of reaction rate decay observed at high pressures by using a hydrogen-free growth environment. This model-guided opportunity to substantially improve synthesis throughput and efficiency, in conjunction with the remarkably invariant CNT properties, is appealing for reliable VACNT device fabrication at both small and industrial scales.
A detailed chemical kinetic model for oxidation of methylamine has been developed, based on theoretical work and a critical evaluation of data from the literature. In this work, the rate coefficients for the reactions of CH 3 NH 2 + O 2 → CH 2 NH 2 / CH 3 NH + HO 2 , CH 3 NH 2 + H → CH 3 + NH 3 , CH 3 NH → CH 2 NH 2 , and CH 3 NH + O 2 → CH 2 NH + HO 2 were calculated from ab initio theory. The mechanism was validated against experimental results from batch reactors, flow reactors, shock tubes, and premixed flames. The model predicts satisfactorily explosion limits for CH 3 NH 2 and its oxidation in a flow reactor. However, oxidation in the presence of nitric oxide, which strongly promotes reaction at lower temperatures, is only described qualitatively. Furthermore, calculated flame speeds are higher than reported experimental values; the model does not capture the inhibiting effect of the NH 2 group in CH 3 NH 2 compared to CH 4 . More work is desirable to confirm the products of the CH 3 NH + NO reaction and to look into possible pathways to NH 3 in methylamine oxidation.
Esters are often produced as unwanted byproducts during the catalytic upgrading of ethanol to diesel fuel precursors through Guerbet coupling. Removal of esters from the product stream is important to prevent the loss of downstream catalyst activity from ester-derived carboxylic acids. In this work, we studied ester hydrogenolysis to the parent alcohols as a viable route for enhanced diesel fuel production. Specifically, we investigated the reduction of hexyl acetate in butanol over ZrO 2 -supported Ni, Co, Cu, Rh, Pd, and Pt catalysts, where Cu/ZrO 2 was the most selective catalyst for the hydrogenolysis of hexyl acetate into hexanol and ethanol. Thermodynamic analysis reveals that a 90% alcohol yield can be obtained at 200 °C, 30 bar, and a relatively high H 2 :hexyl acetate molar ratio of 480:1. Experimentally, an alcohol yield of 88% yield was obtained with a 10 wt % Cu/ZrO 2 catalyst at these conditions with a residence time of 5.4 h kg cat kmol gas –1 . Catalytic tests on the support revealed that ZrO 2 catalyzes the transesterification reaction between hexyl acetate and butanol. However, only the Cu sites can catalyze the hydrogenolysis of the esters into the final alcohols. We developed a kinetic model for our experimental results, which shows that the transesterification and hydrogenolysis reactions run at two different timescales, the former being 10 times faster than the latter. Data regression has been used to develop a model to predict the mole fraction distribution of ester hydrogenolysis products over a wide range of contact times. Cu/ZrO 2 loses half its catalytic activity after 80 h of time on stream. Modeling of deactivation data reveals that the ZrO 2 support conserves a residual activity due to external active sites, while active sites over the Cu surface deactivate at different rates. Furthermore, the catalytic conversion of esters into their parent alcohols is relevant to the production of surrogate liquid fuels since alcohols can be bimolecularly dehydrated to produce a blend of ethers with diesel fuel-like properties.
In a circular economy, organic waste streams are valuable resources for sustainably producing chemicals and fuels. This study investigates a new methane-arrested anaerobic digestion (MAAD) process that converts high-strength cheese whey and brewery wastewater into carboxylic acids. The process was developed and optimized under various bench-scale semi-continuous (fed-batch) operating conditions (e.g., retention time, organic loading rate, pH, and feed/harvest frequency). The MAAD responses to various control-failure scenarios were also systematically investigated. The highest total acid productivity was 26g/(L liq ·d) with a substrate conversion of 0.79g COD digested /g COD fed at a hydraulic retention time (HRT) of approximately 2 d in a 14-L digester. The most stable conditions for digester operation (HRT 3 d at pH 6.0 and 40°C) were selected for process scale-up to 100gal (~380 L). Semi-continuous, pilot-scale MAAD successfully produced a total acid concentration of 40.6+/-1.1 g/L with 8.1% acetic acid, 45.1% butyric acid, and 44.5% lactic acid. The links between wastewater characteristics, operation mode, digester scale, and microbial community structure were statistically analyzed. The results show genera Sporolactobacillus and Clostridium positively correlate with butyric acid production (Pearson correlation coefficient>0.5). Moreover, four kinetic models were developed and fit to batch MAAD experimental datasets (R 2 >95%) and were successfully applied to predict the total acid production in both bench-and pilot-scale semi-continuous MAADs. In conclusion, this study shows MAAD has the potential for industrial-scale applications and is a robust platform to valorize low- or negative-value waste streams into high-value bioproducts.
Increasing peak cylinder pressures and operating temperatures of turbo-charged internal combustion engines (ICEs) engines present new challenges for existing materials employed in engine exhaust valves. Oxidation induced degradation of current Ni based alloys for these components will be a life-limiting mechanism while strict cost margins in the automotive transportation industry further limit the choice of suitable candidate materials. Metallic diffusion aluminide coatings can provide a cost-effective way to improve the oxidation resistance of underlying materials. In this study, the high temperature oxidation behavior of diffusion nickel aluminide coatings on two Ni based alloys (commercially available alloy 31V and a developmental alloy) during cyclic oxidation behavior in air+10% H 2 O at 900 °C was investigated. A complete depletion of the beneficial Al-rich β-(NiFe)Al phase was observed in the coating on alloy 31V while a significant fraction of this phase was retained in the coating on the second alloy. A coupled thermodynamic-kinetic model showed that the disappearance of the Al-rich phase in the coating on 31V was mainly due to the combined effect of steeper chemical potential gradients of Al and Ti between the coating and the substrate and their faster diffusivities compared to the coated developmental alloy. Finally, a higher content of Fe was shown to support the retention of β-(NiFe)Al phase in the coating on the developmental alloy while the presence of Ti in the alloy was shown to be detrimental for long-term coating performance.
Here, this work studies the interaction between colliding plasma jets to understand regimes in which jet merging results in shock formation vs regimes in which the shock formation is mitigated due to the collisionless interpenetration of the jets. A kinetic model is required for this study because fluid models will always produce a shock upon the collision of plasma jets. The continuum-kinetic, Vlasov–Maxwell–Dougherty model with one velocity dimension is used to accurately capture shock heating, along with a novel coupling with a moment equation to evolve perpendicular temperature for computational efficiency. As a result, this relatively inexpensive simulation can be used for detailed scans of the parameter space toward predictions of shocked vs shock-mitigated regimes, which is of interest for several fusion concepts such as plasma-jet-driven magneto-inertial fusion, high-energy-density plasmas, astrophysical phenomena, and other laboratory plasmas. The initial results obtained using this approach are in agreement with the preliminary outcomes of the plasma liner experiment.
We derive axisymmetric equilibrium equations in the context of the hybrid Vlasov model with kinetic ions and massless fluid electrons, assuming isothermal electrons and deformed Maxwellian distribution functions for the kinetic ions. The equilibrium system comprises a Grad–Shafranov partial differential equation and an integral equation. These equations can be utilized to calculate the equilibrium magnetic field and ion distribution function, respectively, for given particle density or given ion and electron toroidal current density profiles. The resulting solutions describe states characterized by toroidal plasma rotation and toroidal electric current density. Additionally, due to the presence of fluid electrons, these equilibria also exhibit a poloidal current density component. This is in contrast to the fully kinetic Vlasov model, where axisymmetric Jeans equilibria can only accommodate toroidal currents and flows, given the absence of a third integral of the microscopic motion.
Ethanol is an attractive oxygenate increasingly used for blending with petroleum-derived gasoline yielding beneficial combustion and emissions behavior for a range of internal combustion engine schemes, including stoichiometric spark-ignition and low temperature combustion (LTC). As such, it is important to fundamentally understand the autoignition behavior of gasoline/ethanol blends. This work utilizes a rapid compression machine (RCM) and a homogeneous charge compression ignition (HCCI) engine to experimentally quantify changes in fuel reactivity, through ignition delay times and preliminary heat release, for blends of 0 to 30% vol./vol. into a full boiling range research gasoline (FACE-F). Diluted/stoichiometric and undiluted/fuel-lean conditions are explored covering a wide range of compressed temperatures and pressures relevant to conventional and advanced, gasoline combustion engines. Detailed chemical kinetic modeling is undertaken using a recently updated gasoline surrogate model in conjunction with a five-component surrogate to model the RCM experiments and provide chemical insight into the perturbative effects of ethanol on the autoignition process. The diluted/stoichiometric RCM measurements reveal that within the low-temperature regime ethanol retards first-stage and main ignition delay times, and suppresses both the rates and extents of low-temperature heat release (LTHR), while within the intermediate-temperature regime ethanol only causes slight changes. Good agreement of ignition delay time and preliminary heat release prediction is found between model and experimental results. Sensitivity and flux analyses further show that ethanol blending effects are dominated by the competition between the H-atom abstraction from ethanol and other fuel components by OH radical at low temperatures and by HO 2 radical at intermediate temperatures. These findings are consistent across both fuel loading conditions explored in this study. In addition, when HCCI engine experiments are mapped onto undiluted/lean RCM measurements under a constant combustion phasing scenario, good correspondence between the two apparatuses is observed for LTHR and start of high-temperature heat release. Finally, the current study highlights the importance of characterizing LTHR in predicting fuel behaviors in high-boost/low-temperature engines, and demonstrates that RCM experiments can provide an alternative, and more-efficient avenue for such characterization.
In this paper, we study the kinetics of polyvinyl chloride (PVC) decomposition using a combination of experimental and computational approaches. Here, we develop a simplified kinetic model that contains only two steps: dehydrochlorination of PVC and further decomposition of the PVC residue. The model is consistent with density functional theory (DFT) calculations and experimental data. Dehydrochlorination is an autocatalytic reaction that starts with a tertiary chloride (Cl) and generates hydrogen chloride (HCl) and benzene as the main products. Benzene and HCl formation rates showed similar trends, indicating that HCl likely catalyzes a homolytic carbon-carbon (C-C) bond cleavage, which gives rise to benzene and an aliphatic fragment. We characterized the structure of dehydrochlorinated PVC (PVC residue) by using thermal gravimetric analysis (TGA), Fourier-transform infrared spectroscopy (FTIR), and nuclear magnetic resonance spectroscopy (NMR). FTIR and NMR results indicate that the PVC residue contains 20% quaternary carbon content, indicating a high concentration of cross-linked molecules. We predict that the most probable structure in the cross-linked centers of the PVC residue is cyclohexadiene, which is supported by DFT calculations, FTIR, and NMR.
Although trivalent manganese (Mn(III)) species have been recognized as crucial intermediates in the degradation of organic contaminants by Mn oxides, quantitative research on their specific roles remains scarce. Here, our study investigated the degradation processes of an organic pollutant, Bisphenol A (BPA), by dissolved Mn(III) and Mn(III)-bearing oxides, and elucidated the differences of the underlying mechanisms and reaction pathways between several Mn oxides and dissolved Mn(III). Our results indicated that BPA degradation rates with Mn(III)-bearing oxides alone follow the order: δ-MnO 2 >> γ-MnOOH > Mn 3 O 4 . Adding pyrophosphate (PP) significantly enhanced BPA degradation by promoting the formation of Mn(III)-PP complexes and exposing more reactive sites, achieved through destabilizing the crystal structure and mitigating of Mn(II) readsorption, particularly in γ-MnOOH and Mn 3 O 4 . Our kinetic model revealed that heterogeneous degradation by Mn oxides is the predominant reaction pathway, accounting for 61.4 %, 87.8 %, and 73.8 % of the total degraded BPA for δ-MnO 2 , γ-MnOOH, and Mn 3 O 4 , respectively, even in the presence of significant amount of dissolved Mn(III) intermediates due to high PP concentrations. These results offer mechanistic details on BPA degradation by Mn oxides and the influence of ligand concentration, providing helpful insights for optimizing degradation strategies of organic pollutants.
The Sabatier reaction (CO 2 + 4H 2 → CH 4 + 2H 2 O) is gaining renewed interest due to its potential to reduce energy carrier storage costs, serve as a feedstock for various organic chemicals, and supply in-situ propellant and life-support resources for long-duration Mars missions. This study demonstrates that combining a modest 2 mA electric field with H 2 feed modulation markedly elevates the CO 2 hydrogenation activity of 2 wt% Ru/CeO 2 catalyst. CO 2 conversion reaches 88 % and 93 % with a CH 4 yield of 83 % and 89 % at 350 °C and 450 °C, respectively. A simple lumped kinetic model reveals that the combined external perturbations not only shift the reaction mechanism but also redistribute key surface-adsorbed intermediates such as hydrogen adatoms and hydrogen-activated CO 2 among the Ru clusters, Ru/CeO 2 interface, and ceria surface. The electric field accelerates the conversion of adsorbed CO 2 to the hydrogenated CO 2 species on Ru and boosts CH 4 formation rate constant, while simultaneously suppresses the formation of undesired, non-reactive surface intermediates. Degree-of-rate-control analysis pinpoints proton migration across the metal-support interface as the decisive lever under these coupled perturbations. In conclusion, these findings establish that rational pairing of metal-support design with well-tuned electric fields and feed oscillations can unlock unprecedented Sabatier rates, guiding the development of next-generation reactors for efficient CO 2 to CH 4 conversion.
This work is an extension of a previous ring current decay model. In the previous work, a two-dimensional kinetic model was constructed to study the temporal variations of the equatorially mirroring ring current ions, considering charge exchange and Coulomb drag losses along drift paths in a magnetic dipole field. In this work, particles with arbitrary pitch angle are considered. By bounce averaging the kinetic equation of the phase space density, information along magnetic field lines can be inferred from the equator. The three-dimensional model is used to simulate the recovery phase of a model great magnetic storm, similar to that which occurred in early February 1986. The initial distribution of ring current ions (at the minimum Dst) is extrapolated to all local times from AMPTE/CCE spacecraft observations on the dawn and dusk sides of the inner magnetosphere spanning the L value range L = 2.25 to 6.75. Observations by AMPTE/CCE of ring current distributions over subsequent orbits during the storm recovery phase are compared to model outputs. In general, the calculated ion fluxes are consistent with observations, except for H+ fluxes at tens of keV, which are always over-estimated. A newly-invented visualization idea, designated as a chromogram, is used to display the spatial and energy dependence of the ring current ion differential flux. Important features of storm-time ring current, such as day-night asymmetry during injection and drift hole on the dayside at low energies (less than 10 keV), are manifested in the chromogram representation. The pitch angle distribution is well fit by the function, j(sub o)(1+Ay(exp n)), where y is sine of the equatorial pitch angle. The evolution of the index n is a combined effect of charge exchange loss and particle drift. At low energies (less than 30 keV), both drift dispersion and charge exchange are important in determining n.