Reaction Rate Equation Research Articles

Long-short-term memory model for quantitative analysis of corrosion products in concrete

A quantitative analysis model is proposed to describe the content evolution of corrosion products in concrete structures subjected to sulfate or chloride attack. The model is inspired by the idea of the Long-Short-Term Memory (LSTM) algorithm in the machine learning field. For establishing the relationship between engineering data and theoretical computations, the training samples of the model are generated by solving the rate equations of the sulfate corrosion reactions in concrete with wide-range initial conditions and validated by the existing experimental data. The Pearson correlation coefficient is employed to determine the model features. The model performance is comprehensively evaluated and discussed. The results show that the proposed model has enough accuracy and feasibility. It could effectively predict the contents of the corrosion products during the sulfate attack by several input values rather than the initial contents of all chemical constituents. The proposed model builds a bridge between experimental methods and theoretical predictions, adequately inheriting their advantages.

Kinetic models for the methanation of COx gases to produce methane: A critical analysis

An exhaustive review of the available literature regarding the kinetic models for the methanation of COx gases to produce natural gas (CH4) is reported. Since the oldest and recent reported kinetic studies about this topic were discussed and analyzed. The characteristics of the experimental studies such as the reaction system, operating conditions, catalyst used, and the experimental procedure, as well as the derived reaction rate equations, reaction rate coefficients and activation energies are highlighted. The results of the various kinetic models are also discussed in order to identify appropriate approaches that can be used as reference to develop new kinetic models. Most of the works reported values of reaction order of 0.2–1 and activation energies of 75–110 kJ/mol in a temperature range of 220–450 °C. Average absolute error (ε) values were calculated to evaluate the accuracy of the predicted data, which resulted to be in the range of 0.15 %–8.08 % showing generally a good agreement.

Homogeneous and heterogeneous ionic liquids catalyze CO2 cycloaddition reaction

The synthesis of cyclic carbonates, which are essential solvents for lithium-ion batteries, through the cycloaddition of carbon dioxide (CO2) and epoxides, and further ester exchange with methanol to form linear carbonates has always been one of the effective ways for CO2 resource utilization. In this work, the kinetics and reactor technology selection of ionic liquids as catalysts for CO2 conversion were analyzed from the perspective of process systems engineering. Firstly, ethanol (EtOH) was used as the solvent, and 1‑butyl‑3-methylimidazolium bromide ([Bmim] Br) was used to catalyze the CO2 cycloaddition reaction in a batch reactor. The effect of reaction conditions on the kinetic performance of synthesizing propylene carbonate (PC) was studied. Secondly, a supported polyelectrolyte liquid membrane was designed and prepared, and a continuous ionic liquid membrane reactor was established for the continuous catalytic reaction of CO2 and epichlorohydrin (PO). The performance of the continuous catalytic reaction was studied and compared with batch reaction experiments. Based on the results of kinetic experiments, the reaction order of the reaction rate equation and the activation energy in the reaction system were determined. Finally, under the premise of achieving high PC yield goals and optimizing process conditions, a dual zone dynamic modeling was conducted on the batch reactor to guide the comprehensive evaluation of technology, economy, and environment under uncertain conditions, laying the foundation for the pilot scale clean production of high-value chemicals synthesized by chemical conversion of CO2.

Analysis of Student Misconception Using Three-Tier Multiple-Choice Diagnostic Test on Reaction Rate Material

Misconception or conceptual errors in learning are errors in understanding a meaning. This descriptive research aims to analyze the level of students' misconceptions on the reaction rate material and analyzing misconceptions on each concept in the reaction rate material. The reseach method used is descriptive research with a quantitative approach. The sample in this research are 30 students of Grade XI students of SMAN 2 Percut Sei Tuan. The instrument used in this research is three-tier multiple-choice diagnostic test, with total 20 questions.. The results showed that the average of students’ misconception of SMAN 2 Percut Sei Tuan who experience misconception is 61,6%. Students’ misconception occurs in each concepts of reaction rate material. Percentages of students' misconceptions on reaction rate and reaction rate expressions have percentage 54.25%, in concept reaction order have percentage 62.87%, in concept factors affecting reaction rate have percentage 60.58%, in concept collision theory have percentage 73.70%, and in concept reaction rate equation have percentage 52.58%.

Analysis of Chemical Kinetics of Multistep Reactions by Mean Reaction Time.

Although methods to derive the rate equation from a kinetic model have been known for over a century, it remains mathematically challenging to derive the rate equation for complex reactions involving multiple steps, as the derivation requires a solution for simultaneous differential equations. Furthermore, the derived kinetic equations are often difficult to intuitively understand. Here, we report a radically different approach to analyze chemical kinetics using the mean reaction time, the average of the time required for the completion of a chemical reaction. This mathematical description of chemical kinetics provides an intuitive understanding of how individual steps in a multistep reaction contributes to the overall kinetics. We demonstrate here the step-by-step approach to derive the formula for the mean reaction time from kinetic models. Surprisingly, one can derive the formula for the mean reaction time without solving simultaneous differential equations or using a steady-state approximation. Being an ensemble-averaged value, the mean reaction time cannot describe the distribution of the reaction time of the individual chemical entities. However, the mean reaction time reveals invaluable insights on how the energy levels of the ground states and the transition states affect the kinetics of the multistep reaction. As proof of principle, we apply the mean reaction time to enzyme kinetics and demonstrate that one can readily derive the expressions of the kinetic parameters (k cat/K M and k cat) even without deriving the Michaelis-Menten equation.

Exploring tubular steady‐state laminar flow reactors with orthogonal collocation

AbstractThe laminar flow reactor (LFR) is one of the most comprehensive problems in chemical reaction engineering, as its modeling involves mass, heat, and momentum conservation equations coupled with chemical reaction rate equations. It is relatively easy to grasp its basic operation, but the solution of the problem is far from trivial. Although there are analytical solutions that simplify the problem, students and tutors must use efficient numerical strategies to appropriately solve some LFR problems. In the present investigation, we solve several different cases of two‐dimensional cylindrical LFR, beginning with the isothermal case, as a benchmark. Calculations were performed using the Chebyshev orthogonal collocation technique by custom scripts in Scilab, Mathematica©, and Matlab®, and were compared with solutions available from the ECRE Version of COMSOL®, which was used as the reference, as well as available analytical solutions in the literature. After analyzing the case of the isothermal LFR with Newtonian fluid, we explored non‐Newtonian fluids, including Carreau and Bingham fluids, whose LFR results are not available in the preceding literature. For Newtonian fluids, besides (a) the isothermal case, we also explore other three nonisothermal design cases: (b) cooling jacket with a fixed heat exchange coefficient at the wall, (c) adiabatic operation, (d) nonisothermal LFR with isothermal wall. Through a performance criterion, the different operation models are compared and we show that nonisothermal design cases perform better than the isothermal case. Computation times, for different scenarios, are quite short and taking 25 nodal points or more suffice for an accurate and timely solution of the more complex problem Case (b). The numerical modeling approach is useful from a pedagogical standpoint, as one can compare numerical results with classical assumptions, and progress from a more restrictive conceptual model (segregated flow) to an array of different kinds of operation with the LFR in which one needs to consider diffusion and temperature distribution.

Modeling of Lithium-Flow Batteries with Aqueous Electrolyte

Flow batteries with aqueous electrolyte at the cathode are safe, low-cost and ecofriendly energy storage systems. When the aqueous cathode is coupled with a Li-metal anode with organic electrolyte, these batteries have extremely high energy densities (over 4,000 Wh/Kg), which makes them attractive not only for grid applications (where flow batteries are often used) but also in maritime applications such long endurance surface vessel systems and wave gliders [1-4]. In this presentation, we derive the mathematical model for lithium-flow batteries and use this model to predict the performance of the batteries under different discharge conditions, cathode microstructures, water flow velocities, etc. The results are then compared to experimental data already available in the literature. The mathematical model is based on coupling the flow equations (Navier-Stokes) with electrochemical equations for reaction rates and using average microscopic models for the cathode. For instance, Fig. 1 shows the voltage of a flow battery as a function of time for different water velocities in a flow battery simulated by solving the above model. For low water velocities, the water is not able to provide enough oxygen in the cathode and the voltage of the battery collapses after a few seconds. However, as long as the water velocity is above approx. 1.2 cm/s, the battery will have enough oxygen to operate continuously; in this case, the higher the flow of the water, the higher the voltage of the battery in the stationary state will be because of lower overpotential losses. In these simulations, the anode contained Li-metal. The cathode had a length of 1 cm (measured in the direction of water flow) and a thickness of 100 μm. The cell was discharged at a small rate of 0.01 mA/cm2. The incoming water had an oxygen concentration of 5 mg/l (which is equal to the oxygen concentration in seawater) and was pumped through the cathode at different velocities. A detailed derivation of the mathematical model and more simulation results and predictions will be presented at the meeting.

A theory of phonon-induced friction on molecular adsorbates

In this manuscript, we provide a general theory for how surface phonons couple to molecular adsorbates. Our theory maps the extended dynamics of a surface's atomic vibrational motions to a generalized Langevin equation, and by doing so captures these dynamics in a single quantity: the non-Markovian friction. The different frequency components of this friction are the phonon modes of the surface slab weighted by their coupling to the adsorbate degrees of freedom. Using this formalism, we demonstrate that physisorbed species couple primarily to acoustic phonons while chemisorbed species couple to dispersionless local vibrations. We subsequently derive equations for phonon-adjusted reaction rates using transition state theory and demonstrate that these corrections improve agreement with experimental results for CO desorption rates from Pt(111).

A clean metallurgical process for vanadium precipitation from vanadium-rich solutions

Ammonium salt vanadium precipitation process has excellent performance in industrial production, but its shortcomings such as high cost and environmental impact are not conducive to sustainable development. In this paper, ethylenediamine was selected as the precipitant, and vanadate was precipitated from acidic simulated vanadium-rich solutions by the ammonium salt precipitation process. The high purity V2O5 product was obtained after calcination. Under optimal conditions, the precipitation rate of vanadium reached 94.8 %, and the purity of vanadium was as high as 99.1 %. The vanadium precipitation product was hexamethyl ammonium, and the vanadium product was V2O5. It can be derived by the differential method that the reaction order was n = 2.4, the reaction rate equation was CV-1.4=kt+a, and the apparent activation energy was 27.8097 kJ/mol, which can be classified as the diffusion-controlled reaction. This innovative process enables a savings of 5050 yuan per ton of V2O5. On this basis, a complete V2O5 production process with zero waste discharge was designed.

ELECTROSYNTHESIS OF Cu2O IN THE Cu(s)|CuSO4(aq)|Fe(s) REACTOR AND DEGRADATION POTENT FOR METHYLENE BLUE

Research on Cu synthesis2O uses the electrolysis method in the Cu(s)|CuSO4(aq)|Fe(s) system and degradation tests on methylene blue. Cu2O is a blackish-red solid which has properties as a semiconductor with a band gap of around 2-2.2 eV. These transition metal oxides are widely used as catalysts, solar cell components, and photodetectors. This research aims to synthesize Cu2O with variations in electrolysis time, determine the characteristics of Cu2O using FTIR and XRD, and test the degradation of methylene blue. The method used is electrolysis with a fixed potential and time variations. Research stages include solution preparation, determination of decomposition potential, electrolysis of samples under various conditions, and analysis. The results of determining the decomposition potential are used as the application potential value in the copper electrolysis process, namely 3.25 V. Sample electrolysis was carried out at a potential of 3.25 V, solution pH 4.0 in a volume of 100 mL at 30, 60, 90, and 120 respectively. and 150 minutes. The products obtained were analyzed using Fourier Transform Infra-Red (FTIR) and X-ray diffraction (XRD). The results of FTIR analysis show that before calcination there is a Cu-O functional group with a wave number of 477.23 cm-1 and after calcination, there is a Cu-O functional group with a wave number of 538.28 cm-1 which indicates the formation of oxide. The results of XRD analysis show that the value of 2θ before calcination is 36.553; 50,481 ; 62.855 which shows the peak of Cu and after calcination 29.928; 36,027 ; 61.568 shows the Cu2O peak. The highest photocatalytic activity test results for Cu2O against methylene blue were 85.96%. The rate of degradation of methylene blue by Cu2O follows the second order reaction rate equation with the R2 value and reaction rate constant being 0.9418 and 0.00769572 L mol-1 s-1, respectively.

Expansion displacement mechanics model of concrete under seawater corrosion

In this study, the variation in the expansion displacement of concrete samples with different water-cement ratios under five corrosion solutions (single sulfate salt and coupled sulfate-chloride salt) is explored. The expansion displacement evolution of these concrete samples under sulfate corrosion (single salt corrosion) and sulfate-chloride corrosion (double salt corrosion) is comprehensively examined. The results reveal that the continuous accumulation of corrosion damage eventually manifests in the form of expansion displacement. Based on the experimental results and the chemical reaction rate equation of the delayed ettringite formation and Friedel’s salt generation, an evolution model of expansion force is established. According to this model and the Weibull distribution law of damage, a expansion displacement mechanics model is proposed to predict the expansion displacement behavior of concrete under sulfate corrosion as well as combined sulfate-chloride corrosion.

Synthesis of Propylene Glycol Methyl Ether Acetate: Reaction Kinetics and Process Simulation Using Heterogeneous Catalyst

Propylene glycol methyl ether acetate (PGMEA) serves as a crucial solvent in semiconductor and display material processes, demanding high purity and low acidity. Despite its significance, its conventional synthesis method using homogeneous catalysts requires extensive purification. Our study explores the use of Amberlyst-15, a stable solid catalyst, to streamline this process. Through batch reactions with a 1:1 reactant ratio at various temperatures and modeling using an integrated reaction rate equation, we obtained kinetic parameters. These parameters were used to predict the kinetics under different reactant ratios and different catalyst amounts, and the predictions match well with experimental results, especially when we used the catalyst amount scaled by the amount of the limiting reactant (PGME) rather than the total amount of the reactants. This highlights the importance of reporting kinetic parameters with proper scaling for catalyst used. Furthermore, we integrated these parameters into process simulations to determine the length of a plug flow reactor (PFR), constructed a PFR system, and confirmed that the simulation results matched well with experimental data obtained from the PFR system. Our findings suggest Amberlyst-15’s potential in simplifying PGMEA synthesis, promising advancements in industrial applications.

Liquid Hot Water (LHW) and Hydrothermal Carbonization (HTC) of Coffee Berry Waste: Kinetics, Catalysis, and Optimization for the Synthesis of Platform Chemicals

Colombia is the world’s leading producer of mildly washed arabica coffee and produces 12.6 million bags of green coffee, but at the same time, 784,000 tons of waste biomass are dumped in open fields, of which only 5% is recovered or used. The objective of this project was to evaluate the production of platform chemicals from these coffee wastes for sustainable resource management. To achieve this, biomass characterization was carried out using proximate analysis, ultimate analysis, and structural analysis. Hydrothermal valorization was carried out at a temperature range of 120–180 °C (LHW) and 180–260 °C (HTC) for one hour. The platform chemicals obtained were quantified by HPLC-RI and monitored by pH and conductivity, and the solid fraction was characterized by monitoring the functional groups in IR spectroscopy and elemental analysis. Hydrolysis processes were obtained at 150 °C, production of platform chemicals at 180 °C, and maximum concentration at 180 °C-4 h; over 200 °C, degradation of the products in the liquid fraction starts to take place. Homogeneous basic and acid catalysts were used to improve the yields of the reaction. The kinetics of the hydrolysis of lignocellulosic structures to sugars were also analyzed and described, and reaction orders of 1 (LHW), 3 (HTC), and their respective reaction rate equations were reported.

The Impact of Effective Matter Mixing Based on Three-dimensional Hydrodynamical Models on the Molecule Formation in the Ejecta of SN 1987A

To investigate the impact of matter mixing on the formation of molecules in the ejecta of SN 1987A, time-dependent rate equations for chemical reactions are solved for one-zone and one-dimensional (1D) ejecta models of SN 1987A. The latter models are based on the 1D profiles obtained by angle-averaging of the three-dimensional (3D) hydrodynamical models, which effectively reflect the 3D matter mixing; the impact is demonstrated, for the first time, based on 3D hydrodynamical models. The distributions of initial seed atoms and radioactive 56Ni influenced by the mixing could affect the formation of molecules. By comparing the calculations for spherical cases and for several specified directions in the bipolar-like explosions in the 3D hydrodynamical models, the impact is discussed. The decay of 56Ni, practically 56Co at later phases, could heat the gas and delay the molecule formation. Additionally, Compton electrons produced by the decay could ionize atoms and molecules and could destroy molecules. Several chemical reactions involved with ions such as H+ and He+ could also destroy molecules. The mixing of 56Ni plays a nonnegligible role in both the formation and destruction of molecules through the processes above. The destructive processes of carbon monoxide and silicon monoxide due to the decay of 56Ni generally reduce the amounts. However, if the molecule formation is sufficiently delayed under a certain condition, the decay of 56Ni could locally increase the amounts through a sequence of reactions.

A parameter determination method of powder burn model based on Levenberg-Marquardt optimization

ABSTRACT The powder burn model is a mathematical model used to characterize the combustion behavior of energetic materials. In this paper, a parameter determination method for the powder burn model is proposed based on the Levenberg-Marquardt (L-M) optimization algorithm, establishing connections between internal ballistics parameters and powder burn model parameters. Taking a 4/1 single-base propellant as an example, the detailed fitting process is presented, and error analysis is conducted. The fittings of the gas state equation and reaction rate equation of powder burn model converge at the 15th and 22nd iterations, respectively. Both fittings displayed a high level of accuracy, with an R-squared value of 0.999. The comparison with the results of the closed bomb experiments demonstrates excellent agreement between pressure-time curves calculated by the powder burn model and the internal ballistics model, confirming their suitability for describing the combustion behavior of propellants. This method provides a foundation for numerically simulating the internal ballistics combustion behavior of energetic materials.

Correct Use of the Bhatia and Perlmutter Random Pore Model under Nonisothermal Conditions.

The classic random pore model (RPM) of Bhatia and Perlmutter [A random pore model for fluid-solid reactions: I. Isothermal, kinetic control. AIChE J.1980, 26, 379-386], which has been used in the kinetic analysis of numerous gas-solid processes involving porous materials under reaction control, was originally derived only for isothermal conditions. This has not prevented many researchers from using it in nonisothermal conditions, a frequent procedure whose validity has not been demonstrated until now. In this work, different questions are answered: Are the isothermal equations of the RPM valid under nonisothermal conditions? And, related to this question, is the commonly used conversion function correct? This work provides the scientific basis that was missing until now for the use of the RPM in nonisothermal conditions, as well as the correct conversion function for such use over the entire particle size range. The said RPM conversion function has been obtained from the exact formulation of the reaction rate equation, and as it should have been originally, it has the grain model function as one of its particular solutions. In this way, the gap that prevented the use of the random pore model for very low particle size values has finally been closed.

Kinetics of Homogeneous Reaction of Potassium Methoxide Based on K2CO3 Catalyst in Transesterification of RBDPO to Biodiesel

Biodiesel production is generally catalyzed by potassium methylate or sodium methylate catalysts based on KOH and NaOH and these catalysts are still imported. The search for a cheap and effective catalyst continues to be carried out by researchers. One of the catalyst support materials currently in use involves impregnating K2CO3 with various substances, resulting in a heterogeneous catalyst. In this study, it was tried to use K2CO3 dissolved in methanol to produce a homogeneous potassium methylate catalyst. Potassium methylate-based homogeneous catalyst K2CO3-methanol is proven to have a very high function in the transesterification reaction of Refined Bleached Deodorized Palm Oil (RBDPO) into biodiesel, this is evidenced by the use of a catalyst percentage of 2% w and 30% w methanol to the weight of RBDPO resulting in an acid content in biodiesel of only 0.12% and a total glycerol of 0.124% in reaction time 3 hours, with the purity of the methyl ester in biodiesel reaching 98.80%. Meanwhile, for the calculation of homogeneous reaction kinetics, a reaction rate equation is produced where the order of the RBDPO transesterification reaction is order 2 (two) and the reaction rate constant is 0.0044.

BioModME for building and simulating dynamic computational models of complex biological systems.

Molecular mechanisms of biological functions and disease processes are exceptionally complex, and our ability to interrogate and understand relationships is becoming increasingly dependent on the use of computational modeling. We have developed "BioModME," a standalone R-based web application package, providing an intuitive and comprehensive graphical user interface to help investigators build, solve, visualize, and analyze computational models of complex biological systems. Some important features of the application package include multi-region system modeling, custom reaction rate laws and equations, unit conversion, model parameter estimation utilizing experimental data, and import and export of model information in the Systems Biology Matkup Language format. The users can also export models to MATLAB, R, and Python languages and the equations to LaTeX and Mathematical Markup Language formats. Other important features include an online model development platform, multi-modality visualization tool, and efficient numerical solvers for differential-algebraic equations and optimization. All relevant software information including documentation and tutorials can be found at https://mcw.marquette.edu/biomedical-engineering/computational-systems-biology-lab/biomodme.php. Deployed software can be accessed at https://biomodme.ctsi.mcw.edu/. Source code is freely available for download at https://github.com/MCWComputationalBiologyLab/BioModME.

First-principles-based microkinetic modeling of methanol steam reforming over Cu(111) and Cu(211): structure sensitive activity and selectivity.

The development of hydrogen energy is widely recognized as a key approach to addressing the energy and carbon emission challenges. Methanol steam reforming is a promising hydrogen production scheme that can provide high-purity hydrogen. In this work, we studied the primary reaction mechanisms of methanol steam reforming over the Cu(111) and Cu(211) surfaces using density functional theory (DFT) calculations and microkinetic simulations. A detailed kinetic perspective on the reaction mechanism, which is often overlooked in previous research that relies solely on DFT calculations, is provided in the current work. Our findings reveal that under typical experimental conditions, the dominant mechanism on the Cu(111) surface is the methyl formate mechanism, while the H2COO dehydrogenation mechanism is dominant on Cu(211). The activity over the Cu(111) surface was slightly higher than that over Cu(211). Based on the degree of rate control analysis results, a reaction rate equation was derived to quantitatively explain the trend of activity under different operating conditions. It was also found that CO2 selectivity was significantly higher over Cu(211) than over the Cu(111) surface. Furthermore, based on the Wulff construction scheme, copper nanoparticle models with different sizes were constructed, and a detailed structure sensitivity study was executed. This comprehensive investigation sheds light on the mechanisms of methanol steam reforming reactions over the Cu(111) and Cu(211) surfaces, providing essential insights for the design of high-performance catalysts for hydrogen production.

First principle based rate equation (1pRE) for reduction kinetics of Fe2O3 with syngas in chemical looping

Chemical looping combustion (CLC) of solid fuel is a promising technology with inherent CO2 separation and low energy penalty for CO2 capture. Syngas is main intermediate species of solid fuel conversion in CLC, the reduction kinetics of oxygen carriers with syngas play a crucial role in CLC systems. However, the current research obtains the reaction rate constants by fitting the apparent models with the experimental data, and cannot explain the reduction kinetics behavior from a microscopic level. It remains a challenge to compute the reduction kinetics of oxygen carriers with syngas directly from first-principles density functional theory (DFT) without fitting experimental data. This study proposes a first-principle-based rate equation (1pRE) theory and integrates it into the random pore model (RPM) to predict the kinetics of Fe2O3 reduction by syngas in CLC. The developed 1pRE theory utilizes DFT calculations to search for reaction pathways and energy barriers of elementary reactions. Then the DFT data are introduced into the statistical mechanics partition function and transition state theory (TST) to calculate the reaction rate constants. Microkinetic rate equations of elementary reactions occurring at the surface scale are developed to describe the change of surface coverage of different surface species. The 1pRE theory is integrated into the RPM to account for the influence of particle-scale structural changes on the overall conversion rate during the reduction process. The theory can predict the reduction kinetics of oxygen carriers without fitting experimental data and establishes a connection between microscopic insights and macroscopic phenomena. The accuracy was validated by experimental data of Fe2O3 oxygen carriers obtained from the thermogravimetric analyzer (TGA) in the atmosphere of syngas. The developed 1pRE predicts the reduction kinetics of oxygen carriers accurately and can be used to optimize the design of oxygen carrier materials and the scale up of CLC reactors.