Determination of Kinetic Parameters for Carbon Fiber and Glass Fiber Epoxy Composites Decomposition.
Abstract not provided.
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.
Abstract not provided.
Reactor temperature manipulation to increase product yields of chemical reactions is a known technique used in many industrial processes. In the case of exothermic chemical reactions, the well-known Le Chatelier’s principle predicts that a decrease in temperature will displace the chemical reaction toward the formation of products by increasing the value of the equilibrium constant. The reverse is true for endothermic reactions. Reactor temperature manipulation in an industrial system, however, affects the values of many variables, including physical properties, transport parameters, reaction kinetic parameters, etc. In the case of reactive absorption, some variables change with increasing temperatures due to solute absorption, while others change in such a way that the solute absorption rate decreases. For example, temperature drop increases product formation for exothermic reactions but reduces the value of transport parameters, leading to decreasing interfacial concentrations and absorption rates. Therefore, temperature manipulation strategies must be designed carefully to achieve the process goals. In this work, we theoretically study the use of temperature as a tool to increase CO 2 absorption by solvents in a semi-batch reactor. A computer code has been developed and validated using reported experimental data. Calculated results demonstrate an increase in absorbed CO 2 of more than 28% with respect to the highest temperature used. Despite high agitation and high gas flow rate, the system is mass transfer controlled at short times, becoming kinetically controlled as time increases. An operating strategy to decrease cooling energy costs is also proposed. This study reveals that reactor temperature manipulation can be an effective process to improve CO 2 absorption by solvents in two-phase semi-batch reactors.
In this study, the catalytic pyrolysis of low-density polyethylene (LDPE), polyethylene terephthalate (PET) and their mixture (1 : 1 wt/wt) with three zeolite catalysts (HZSM-5, H-beta, HY) was investigated using a thermogravimetric analyzer (TGA) and a Pyroprobe® micro-reactor coupled to a gas chromatograph mass spectrometer (Py-GC/MS). The TGA results demonstrated that during pyrolysis at 10 °C min –1 , on average, zeolite catalysts decreased the maximum decomposition temperature by 149 °C for LDPE while only decreasing by 8 °C for PET. The derivative thermogravimetric (DTG) curve evidenced interactions when the two polymers were catalytically co-pyrolyzed for all the three catalysts. A lumped nth order reaction scheme was able to accurately model both non-catalytic and catalytic pyrolysis and co-pyrolysis by using the least squares fitting approach for determining the kinetic parameters. The kinetic model was able to model well the interaction effects observed during catalytic co-pyrolysis of LDPE and PET with HZSM-5, H-beta, and HY (Fit% Wt% > 96%, Fit% DTG > 94%). Py-GC/MS experiments for the catalytic fast pyrolysis of LDPE revealed that HZSM-5 resulted in the highest selectivity to aromatic hydrocarbons (31.6%) and HY resulted in the highest selectivity to gasoline range C 5 –C 10 paraffins and olefins (40.9%). Catalytic fast pyrolysis of PET showed high selectivity to benzene for all the catalysts (>43%) and that HZSM-5 resulted in the highest selectivity to polyaromatic hydrocarbons (24.7%). Finally, the catalytic fast co-pyrolysis of LDPE and PET revealed interaction effects for all the three catalysts evidenced by a positive synergy% for alkylated benzenes (3–142%) and polyaromatics (105–187%) with a concomitant negative synergy% for benzene (24–36%) and C 5 –C 10 paraffins and olefins (27–53%).
Oxidative cleavage of alkenes and unsaturated fatty acids with hydrogen peroxide gives an efficient and sustainable process to obtain mono- and di-acids for polymers and lubricants with fewer safety risks and less environmental impact than processes that utilize ozone or other inorganic oxidizers (e.g., permanganate, dichromate, etc.). Guided by insight into the mechanisms for competing reaction pathways (i.e., epoxidation of alkene on W–(η2-O2) complexes vs. H2O2 decomposition) and the apparent kinetics derived from kinetic experiments, here, we postulate that W-based heterogeneous catalysts can provide high performance and stable operations at low H2O2 concentrations. Semi-batch reactors with continuous introduction of H2O2 solutions offer the means to maintain low H2O2 concentrations while providing sufficient quantities of H2O2 to satisfy the reaction stoichiometry. We derived simple kinetic model equations for the epoxidation, ring-opening, oxidative cleavage, and oxidation steps and fit theses equations to batch experimental data to obtain kinetic parameters. This kinetic model describes the concentration profiles of reactant, oxidant, and products well as shown by agreement with experimental data. Further predictions of the optimal H2O2 feed rate for semi-batch operation utilized by the proposed rate expressions and the reactor design equations suggest that low H2O2 feed rate increases selectivity towards oxidative cleavage products and selective use of H2O2 for oxidative cleavage pathway. Comparisons of oxidative cleavage of 4-octene in batch and semi-batch reactors show that semi-batch reactors with optimized molar feed rates of H2O2 increased oxidative cleavage product selectivities (76% to 99%; with an increase in butyric acid selectivity from 1% to 55%) and H2O2 selectivity (3% to 30%). In addition, semi-batch reaction conditions used avoid H2O2-mediated dissolution of W-atoms from the catalyst. Analysis of these findings suggest that solid oxide catalysts will be effective for continuous oxidative cleavage reactions if deployed within fixed-bed reactors that allow for distributed introduction of reactants and therefore low in situ concentrations of H2O2.
The prevalent use of organic materials in manufacturing is a fire safety concern, and motivates the need for predictive thermal decomposition models. A critical component of predictive modeling is numerical inference of kinetic parameters from bench scale data. Currently, an active area of computational pyrolysis research focuses on identifying efficient, robust methods for optimization. This paper demonstrates that kinetic parameter calibration problems can successfully be solved using classical gradient-based optimization. We explore calibration examples that exhibit characteristics of concern: high nonlinearity, high dimensionality, complicated schemes, overlapping reactions, noisy data, and poor initial guesses. The examples demonstrate that a simple, non-invasive change to the problem formulation can simultaneously avoid local minima, avoid computation of derivative matrices, achieve a computational efficiency speedup of 10x, and make optimization robust to perturbations of parameter components. Techniques from the mathematical optimization and inverse problem communities are employed. By re-examining gradient-based algorithms, we highlight opportunities to develop kinetic parameter calibration methods that should outperform current methods.
We undertook a detailed rheological investigation to evaluate the kinetic parameters of the forward and reverse Diels–Alder (DA) reactions of a model network cross-linked using a furan prepolymer and a common aromatic bismaleimide. At high temperature where the Winter–Chambon’s criterion of frequency-independence was more applicable, a multiwave technique permitted van’t Hoff analysis and calculation of the reaction thermodynamic parameters, specifically the enthalpy and entropy of the reaction: ΔH° = –38.3 ± 5.2 kJ mol –1 and ΔS° = –94.3 ± 13.4 J mol –1 . At mild temperatures where the G'–G" crossover point is experimentally convenient to measure gelation, isothermal tests were used to obtain reasonable fDA kinetic parameters from Eyring analysis such as the apparent activation enthalpy and entropy of ΔH$^{‡}_{fDA}$ = 76.8 ± 6.9 kJ mol –1 and ΔS$^{‡}_{fDA}$= –82.8 ± 22.2 J mol –1 K –1 . Comparable rheokinetic methods include cross-linking density measurements and stress relaxation tests to calculate effective kinetics, whereas the critical gel conversion was consistently applied here. As a result, rate data are fitted with the Arrhenius equation for comparison purposes and the Eyring equation to demonstrate its broader utility.
The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry-based transport code.
The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code. This manuscript serves as a single manual for two separate codes: VARI3D and PERSENT. The VARI3D code (VARIational 3D) is based upon the classic finite difference diffusion theory solver available in DIF3D. The PERSENT code (PERturbation and SENitivity for Transport) is based upon the variational nodal method employed in DIF3D termed VARIANT. The VARIANT solver was added to DIF3D in 1995 and has seen continued development and use for the last 18 years. Because VARI3D primarily uses deprecated coding practices, rather than incorporating the perturbation and sensitivity treatments for transport within VARI3D, a new coding development was built using modern Fortran coding. The primary purpose of this manual is to describe the theory behind PERSENT (and by convenience, that of VARI3D) and discuss the input and output of PERSENT along with giving potential users an idea of how to use it. While this manuscript does describe the input and output of VARI3D, the PERSENT code is intended to be the replacement capability of VARI3D as PERSENT can generate nearly identical (if not superior) diffusion theory results. In this manuscript, the relevant aspects of generalized perturbation theory and exact perturbation theory that apply to both VARI3D and PERSENT are covered. The input and output of VARI3D is displayed by excerpting several of the example problems. Similarly, the input and output of PERSENT is displayed along with tips on how best to use the code. Note that the input and output of the inhomogeneous solver wrapped around DIF3D (DIF3D_IFS) is also discussed as it is needed to carry out some of the sensitivities in PERSENT such as reaction rate ratios. This manuscript describes several perturbation and sensitivity problems, and the results computed using PERSENT. From these sections, potential users should find that PERSENT provides not only the typical tables of numbers desired in perturbation and sensitivity analysis work, but also can visually plot the result for a more thorough understanding of the space and energy distribution (Section 5). Overall, PERSENT is observed to produce accurate reactivity worths and sensitivities for the displayed set of test problems and clearly demonstrates the need to have a transport-based sensitivity capability as evident from the thousands of percent errors observed in the 21-group hexagonal fast reactor problem (covered in Section 7). The uncertainty calculation capability is described in Section 3 and demonstrated in Section 7.
Visible photon fluxes can influence the rate and selectivity of heterogeneously catalyzed reactions on metal nanoparticle surfaces. Models describing the influence of photon fluxes have typically introduced photon flux dependent apparent thermal kinetic parameters (reaction orders, activation energies, binding energies, etc.). This has relied on empirical fitting of reaction rate data, making mechanistic interpretations of how photon fluxes influence elementary step rates challenging and inconsistent with fundamental descriptions of photochemistry on metal surfaces developed from surface science studies. Using the CO adsorption–desorption quasi-equilibrium reaction on Pt/Al 2 O 3 catalysts as a model system, we measured steady state adsorbed CO (CO*) coverages under isothermal and isobaric (1 mbar CO) conditions as a function of temperature (473–573 K) and of 440 nm photon flux ((0.1–5.2) × 10 3 # hv Pt site –1 s –1 ) using in situ IR spectroscopy. Steady state CO* coverage on Pt was photon flux dependent with increasing photon flux causing decreasing coverage, consistent with photons driving CO* desorption rates faster than thermal CO* desorption rates. However, photon flux dependent CO* coverages were essentially temperature independent, inconsistent with models that describe photon effects using perturbations to apparent thermal kinetic parameters. Instead, 120 steady state CO* coverages as a function of temperature and photon flux are quantitatively described by a kinetic model in which the overall desorption rate is a summation of independent thermal and photon induced CO* desorption rates. Site-resolved analysis reveals distinct kinetic parameters for photon driven desorption of CO* from well-coordinated, under-coordinated, and highly under-coordinated Pt sites, with temperature-dependent apparent quantum efficiencies (AQE) consistent with temperature dependence of vibrational quanta distribution of adsorbed CO. The rigorous kinetic rate laws for independent photon and thermal driven pathways allow for predictive modeling of the influence of photon fluxes on the rates of CO* desorption under catalytic conditions. Further, the analysis provides evidence that steady state continuous wave photon fluxes can drive desorption/adsorption reactions on metal surfaces out of thermal equilibrium, reconciling surface science observations of molecular photodesorption with applied catalysis. The work establishes a general kinetic framework to be considered for photon driven processes on metals, and defines catalyst, reaction, and photon flux characteristic design principles for breaking Sabatier limitations.
The equivalent fundamental-mode source refers to a “source” which is identically distributed in space, energy, and angle as that of a fundamental-mode fission source distribution. It has previously been shown that the system fixed source multiplication ($M_{fs}$) relates to the effective multiplication factor ($\kappa_{eff}$) as $M_{fs} = \frac{g^*}{1-k_{eff}}$ The term $\mathcal{g}$* therefore relates the system source distribution to the system equivalent fundamental-mode source. As previously shown, the $\mathcal{g}$* term is approximately 1 for deeply subcritical systems and then diverges as the system multiplication increases. For a system with a point source located in the center of an assembly (such a sphere), $\mathcal{g}$* will be greater than 1; this happens because there is a higher probability of induced fission when all of the starter neutrons originate from the center of the sphere. For a uniform source, $\mathcal{g}$* is less than 1, but does not diverge from 1 as dramatically as seen by a point source. One reason that $\mathcal{g}$* is an important parameter is that it is used in measurement methods to infer reactor kinetics parameters ($β_{eff}$ being one example). This work will discuss multiple ways to simulate $\mathcal{g}$* using MCNP®6.2, which include use of traditional methods versus a single input file method that was recently published. A discussion related to methods to simulate $M_{fs}$ will be presented. Spherical systems of Highly Enriched Uranium (HEU) will be compared with previous works. In addition, results of plutonium including multiple reflector materials will be investigated. Simulated systems will include both parametric studies as well as configurations used in recent experiments performed at the National Criticality Experiments Research Center (NCERC). These studies will include subtleties associated with system geometry that have not been previously shown. Last, this work will investigate the relationship between $\mathcal{g}$* and reactor kinetics parameters (such as $β_{eff}$).
Prior engineering of the ethanologen Zymomonas mobilis has enabled it to metabolize xylose and to produce 2,3-butanediol (2,3-BDO) as a dominant fermentation product. When co-fermenting with xylose, glucose is preferentially utilized, even though xylose metabolism generates ATP more efficiently during 2,3-BDO production on a BDO-mol basis. To gain a deeper understanding of Z. mobilis metabolism, we first estimated the kinetic parameters of the glucose facilitator protein of Z. mobilis by fitting a kinetic uptake model, which shows that the maximum transport capacity of glucose is seven times higher than that of xylose, and glucose is six times more affinitive to the transporter than xylose. With these estimated kinetic parameters, we further compared the thermodynamic driving force and enzyme protein cost of glucose and xylose metabolism. It is found that, although 20% more ATP can be yielded stoichiometrically during xylose utilization, glucose metabolism is thermodynamically more favorable with 6% greater cumulative Gibbs free energy change, more economical with 37% less enzyme cost required at the initial stage and sustains the advantage of the thermodynamic driving force and protein cost through the fermentation process until glucose is exhausted. Glucose-6-phosphate dehydrogenase (g6pdh), glyceraldehyde-3-phosphate dehydrogenase (gapdh) and phosphoglycerate mutase (pgm) are identified as thermodynamic bottlenecks in glucose utilization pathway, as well as two more enzymes of xylose isomerase and ribulose-5-phosphate epimerase in xylose metabolism. Acetolactate synthase is found as potential engineering target for optimized protein cost supporting unit metabolic flux. Pathway analysis was then extended to the core stoichiometric matrix of Z. mobilis metabolism. Growth was simulated by dynamic flux balance analysis and the model was validated showing good agreement with experimental data. Dynamic FBA simulations suggest that a high agitation is preferable to increase 2,3-BDO productivity while a moderate agitation will benefit the 2,3-BDO titer. Taken together, this work provides thermodynamic and kinetic insights of Z. mobilis metabolism on dual substrates, and guidance of bioengineering efforts to increase hydrocarbon fuel production.
Oxidative cleavage of alkenes and unsaturated fatty acids with hydrogen peroxide gives an efficient and sustainable process to obtain mono- and di-acids for polymers and lubricants with fewer safety risks and less environmental impact than processes that utilize ozone or other inorganic oxidizers (e.g., permanganate, dichromate, etc.). Guided by insight into the mechanisms for competing reaction pathways (i.e., epoxidation of alkene on W–(η 2 -O 2 ) complexes vs. H 2 O 2 decomposition) and the apparent kinetics derived from kinetic experiments, here, we postulate that W-based heterogeneous catalysts can provide high performance and stable operations at low H 2 O 2 concentrations. Semi-batch reactors with continuous introduction of H 2 O 2 solutions offer the means to maintain low H 2 O 2 concentrations while providing sufficient quantities of H 2 O 2 to satisfy the reaction stoichiometry. We derived simple kinetic model equations for the epoxidation, ring-opening, oxidative cleavage, and oxidation steps and fit theses equations to batch experimental data to obtain kinetic parameters. Furthermore, this kinetic model describes the concentration profiles of reactant, oxidant, and products well as shown by agreement with experimental data. Further predictions of the optimal H 2 O 2 feed rate for semi-batch operation utilized by the proposed rate expressions and the reactor design equations suggest that low H 2 O 2 feed rate increases selectivity towards oxidative cleavage products and selective use of H 2 O 2 for oxidative cleavage pathway. Comparisons of oxidative cleavage of 4-octene in batch and semi-batch reactors show that semibatch reactors with optimized molar feed rates of H 2 O 2 increased oxidative cleavage product selectivities (76% to 99%; with an increase in butyric acid selectivity from 1% to 55%) and H 2 O 2 selectivity (3% to 30%). In addition, semibatch reaction conditions used avoid H 2 O 2 -mediated dissolution of W-atoms from the catalyst. Analysis of these findings suggest that solid oxide catalysts will be effective for continuous oxidative cleavage reactions if deployed within fixed-bed reactors that allow for distributed introduction of reactants and therefore low in situ concentrations of H 2 O 2 .
Wastewater is a misplaced resource well suited to recover nutrients, value-added chemicals, energy, and clean water. A photoelectrochemical device is proposed to transform wastewater nitrates to ammonia and nitrous oxide, coupled with water oxidation. Numerical models were developed to quantify the dependence of process efficiencies and nitrogen-removal rates on light absorber band gaps, electrocatalytic kinetic parameters, competing oxygen reduction and hydrogen evolution reactions, and the reacting nitrate species concentrations that affect the mass-transfer limited current densities. With a single light-absorber and state-of-the-art catalysts, optimal solar-to-chemical efficiencies of 7% and 10% and nitrogen-removal rates of 260 and 395 gN m -2 day -1 are predicted for nitrate reduction to ammonia and nitrous oxide, respectively. The influence of competing reactions on the performance depends on the nitrate concentration and band gap of the light absorber modeled. Oxygen reduction is more dominant than hydrogen evolution to compete with the nitrate reduction reaction, but it is mass-transfer limited. Even with kinetic parameters that enhanced the driving forces for the competing reactions, the performance is only minimally affected by these reactions for optimally selected band gaps and nitrate concentrations larger than 100 mM. Here, theoretically predicted peak nitrogen removal rates and specific energy intensities are competitive with reported estimates for electrochemical and Sharon-Anammox processes for ammonia recovery and nitrogen removal, respectively. This result, together with the added benefit of harnessing sunlight to produce value-added products, indicates promise in the photoelectrochemical approach as a tertiary pathway to recover nutrients and energy from wastewater nitrates.
Abstract Accurate, predictive reaction models are critical for the design and optimization of chemical looping combustion (CLC) reactors. The formulation and estimation of kinetic parameters for these reaction models using a first‐principles equation‐oriented (EO) approach is particularly beneficial as large amounts of experimental data spanning process‐relevant conditions can be used to estimate parameters in a computationally tractable way. This work demonstrates the application of a novel EO framework to develop reduction reaction kinetic models of an iron‐based CLC oxygen carrier (OC). An optimization problem is formulated to estimate kinetic parameters that provide the best fit to the experimental data. The model predicts the state of the OC with mean square error values of 2.5%–4.4% across the full range of validation data, including multiple reduction cycles.
Mean-field microkinetic models (MKMs) are appealing for their relatively facile construction, computational tractability, and high-throughput catalyst screening capabilities. As such, they will continue to be a valuable tool for materials design in heterogeneous catalysis even as the field aims to describe more complex systems. Numerous prior reports have provided the groundwork for constructing first-principles-based MKMs, including the analysis of strategies for incorporating lateral interactions into thermodynamic parameters (e.g., adsorption energies). Yet, there remains a need for concerted dialogue on methods for calculating and incorporating coverage-dependent kinetic parameters into MKMs. In this Perspective, we assess strategies for doing so, including the corresponding key physical implications and computational challenges. Here, we emphasize that decoupling thermodynamic and kinetic parameters within MKMs can violate thermodynamic consistency and risk unphysical solutions. For some reactions and catalyst materials, scaling relationships can predict coverage-dependent activation energies, but there are several exceptions evident in the literature, indicating that this approach is not universally applicable and that the field could benefit from research aimed at elucidating the limitations. Conducting high-coverage transition state searches is a rigorous but computationally costly strategy, and the effects of various methods for mitigating this cost on resulting energetics have yet to be broadly explored and validated. The goal of this Perspective is to generate discussion on and inspire focused research into the physical relevance of approaches for describing coverage-dependent reaction barriers in MKMs, including the development of computationally tractable methodologies, to advance the applicability of MKMs across diverse reaction chemistries and conditions.
Determining kinetic parameters such as activation energy (Ea), pre-exponential factor (A), rate constants, and reaction order (n) via thermal analysis techniques is important to material synthesis and fabrication, industrial production, biomedicals, energy storage, catalysis, etc. Various kinetic methods based on the non-isothermal techniques have been developed to obtain the kinetic parameters. In this review, we focus on introducing various kinetic methods in non-isothermal analysis and their applications in different reactions, including thermal decomposition, solid-state phase transformation, crystallization, thermal ignition, curing process, etc.
Previous models for H 2 and H 2 O outgassing from moisture exposed LiH specimens were based on temperature programmed desorption/decomposition (TPD) experiments and analyses that yielded kinetic descriptions which were specific to defined sets of initial conditions. Subsequently the inputs to the kinetic models were refined to correlate with infrared spectra of the adsorbed species. In this current enhancement, better agreement between prediction and experiment was obtained through: 1) refining the reaction mechanisms for some of the decomposition reactions; 2) developing the potential range of kinetic parameters by simultaneous fitting of all the TPD outgassing curves obtained at different heating rates; and, 3) limiting the scope of possible kinetic parameters to an optimized subset by matching the outgassing prediction curves at 3 different temperatures to the outcomes from the accurate but cumbersome, non-intuitive, and effort demanding isoconversional analysis.
Temperature-controlled, time-resolved picosecond electron pulse radiolysis was utilized to measure the rate of reaction between the solvated electron (eS–) and dissolved oxygen in n-dodecane solutions from 2.5 to 40 °C for the first time. At 20.0 °C, the reaction rate was determined to be k(eS– + O2) = (4.54 ± 0.21) × 1010 M-1 s-1, with an activation energy of Ea = 14.4 ± 1.3 kJ mol-1. These newly determined kinetic parameters are important for predicting and managing the effects of aerated environments on the degradation of organic solvents used in nuclear fuel reprocessing technologies.