Reaction computer programs

The text listed below provides more details on how the potentiometric titration data may be used to calculate equilibrium constants. This text provides a number of examples and includes a discussion of several computer programs that have been developed to model equilibrium reactions. 

When the kinetics are unknown, still-useful information can be obtained by finding equilibrium compositions at fixed temperature or adiabatically, or at some specified approach to the adiabatic temperature, say within 25°C (45°F) of it. Such calculations require only an input of the components of the feed and produc ts and their thermodynamic properties, not their stoichiometric relations, and are based on Gibbs energy minimization. Computer programs appear, for instance, in Smith and Missen Chemical Reaction Equilibrium Analysis Theory and Algorithms, Wiley, 1982), but the problem often is laborious enough to warrant use of one of the several available commercial services and their data banks. Several simpler cases with specified stoichiometries are solved by Walas Phase Equilibiia in Chemical Engineering, Butterworths, 1985). 

The holistic thermodynamic approach based on material (charge, concentration and electron) balances is a firm and valuable tool for a choice of the best a priori conditions of chemical analyses performed in electrolytic systems. Such an approach has been already presented in a series of papers issued in recent years, see [1-4] and references cited therein. In this communication, the approach will be exemplified with electrolytic systems, with special emphasis put on the complex systems where all particular types (acid-base, redox, complexation and precipitation) of chemical equilibria occur in parallel and/or sequentially. All attainable physicochemical knowledge can be involved in calculations and none simplifying assumptions are needed. All analytical prescriptions can be followed. The approach enables all possible (from thermodynamic viewpoint) reactions to be included and all effects resulting from activation barrier(s) and incomplete set of equilibrium data presumed can be tested. The problems involved are presented on some examples of analytical systems considered lately, concerning potentiometric titrations in complex titrand + titrant systems. All calculations were done with use of iterative computer programs MATLAB and DELPHI. 

Two standard estimation methods for heat of reaction and CART are Chetah 7.2 and NASA CET 89. Chetah Version 7.2 is a computer program capable of predicting both thermochemical properties and certain reactive chemical hazards of pure chemicals, mixtures or reactions. Available from ASTM, Chetah 7.2 uses Benson s method of group additivity to estimate ideal gas heat of formation and heat of decomposition. NASA CET 89 is a computer program that calculates the adiabatic decomposition temperature (maximum attainable temperature in a chemical system) and the equilibrium decomposition products formed at that temperature. It is capable of calculating CART values for any combination of materials, including reactants, products, solvents, etc. Melhem and Shanley (1997) describe the use of CART values in thermal hazard analysis. 

A copy (5V4 inch floppy disk) of a menu-driven computer program to calculate Gibbs free energy of formation and change in Gibbs free energy for reactions (including random access data file of compound coefficients) is available for a nominal fee. For details, contact C. L. Yaws, Dept, of Chem. Eng. Lamar University, P.O. Box 10053, Beaumont, Texas 77710, USA. 

Equations 5-110, 5-112, 5-113, and 5-114 are first order differential equations and the Runge-Kutta fourth order numerical method is used to determine the concentrations of A, B, C, and D, with time, with a time increment h = At = 0.5 min for a period of 10 minutes. The computer program BATCH57 determines the concentration profiles at an interval of 0.5 min for 10 minutes. Table 5-6 gives the results of the computer program and Figure 5-16 shows the concentration profiles of A, B, C, and D from the start of the batch reaction to the final time of 10 minutes. 

The computer program PLUG51 used Equation 5-334 to determine the conversions and the compositions of the components. Table 5-13 illustrates the results of the computer program, and Figure 5-32 shows the plots of the rates of each reaction as a function of V/F. In both instances, the rates decrease toward zero as V/F increases. Figure 5-33 shows the plot of the total conversion versus V/F. 

Using the same values of the kinetic parameters as in Type 1, and given C o = 0-1 mo 1/1, it is possible to solve Equation 6-155 with Equations 6-127 and 6-128 simultaneously to determine the fractional conversion X. A computer program was developed to determine the fractional conversion for different values of (-iz) and a temperature range of 260-500 K. Eigure 6-30 shows the reaction profile from the computer results. 

The reaction flow-charts of Part Two, and indeed all chemical formulae which appear in this book, were generated by computer. The program used for these drawings was ChemDraw adapted for the Macintosh personal computer by Mr. Stewart Rubenstein of these Laboratories from the molecular graphics computer program developed by our group at Harvard in the 1960 s (E. J. Corey and W. T. Wipke, Science, 1969,166, 178-192) and subsequently refined. 

Huff, J. E., Computer Programming For Runaway Reaction Venting, I. Chem. E. Symposium Series 85, 109-129, 1984. 

Changes in free energy and the equilibrium constants for Reactions 1, 2, 3, and 4 are quite sensitive to temperature (Figures 2 and 3). These equilibrium constants were used to calculate the composition of the exit gas from the methanator by solving the coupled equilibrium relationships of Reactions 1 and 2 and mass conservation relationships by a Newton-Raphson technique it was assumed that carbon was not formed. Features of the computer program used were as follows (a) any pressure and temperature may be specified (b) an inert gas may be present (c) after... 

The reaction enthalpy and reaction entropy were derived from the curves comparing with data for the all-or-none model140 using the computer program of Rosenbrock139 (Table 6). 

The mechanism of radiative transfer in flares was found to depend on compn, flare diameter and pressure (Ref 69). The flare efficiency calcn is complicated by the drop-off in intensity at increasing altitudes and at very large diameters owing to the lower reaction temps (Ref 11, p 13) and the narrowing of the spectral emittance band (Ref 35). The prediction of the light output in terms of compn and pressure (ie, altitude) is now possible using a computer program which computes the equilibrium thermodynamic properties and the luminance (Ref 104) Flare Formulations... 

A new chapter (5) on reaction intermediates develops a number of methods for trapping them and characterizing their reactivity. The use of kinetic probes is also presented. The same chapter presents the Runge-Kutta and Gear methods for simulating concentration-time profiles for complex reaction schemes. Numerical methods now assume greater importance, since useful computer programs are available. The treatment of pH profiles in Chapter 6 is much more detailed. 

Proton transfer reactions, 143-144, 144 activation energy, 149,164 all-atom models for, 146-148 Cys 25-His 159 in papain, 140-143 computer program for EVB calculations, 150-151... 

Radial distribution function, 79 computer program for calculating, 96-106 Rate constant, see Rate of reaction Rate of reaction ... 

Solution The computer programs used for the consecutive reaction examples can be used. All that is needed is to modify the subroutine Reactor. Results are shown in Table 6.5. 

A general methodology has been developed for the treatment of NMR data of polymer mixtures. The methodology is based on reaction probability models and computer optimization methods, resulting in a family of computer programs called MIXCO. The use of MIXCO programs enabled three components to be resolved from NMR tacticity data of fractionated polybutylene. 

The pre-gel model calculates the weight average molecular weight of the reaction mixture, while the post-gel model calculates the weight of the sol fraction and the effective crosslink density. A simple computer program using the derived expressions has been written in BASIC and runs on IBM-PC compatible computers. The importance of secondary reactions on cure in typical coatings is discussed. 

Contemporary s Tithetic chemists know detailed information about molecular structures and use sophisticated computer programs to simulate a s Tithesis before trying it in the laboratory. Nevertheless, designing a chemical synthesis requires creativity and a thorough understanding of molecular structure and reactivity. No matter how complex, every chemical synthesis is built on the principles and concepts of general chemistry. One such principle is that quantitative relationships connect the amounts of materials consumed and the amounts of products formed in a chemical reaction. We can use these relationships to calculate the amounts of materials needed to make a desired amount of product and to analyze the efficiency of a chemical synthesis. The quantitative description of chemical reactions is the focus of Chapter 4. 

More complicated reactions schemes, including first-order reversible consecutive processes and competitive consecutive reactions, are considered in a textbook by Irwin [89]. Professor Irwin s textbook also includes computer programs written in the BASIC language. These programs can be used to fit data to the models described. 

Gasteiger, J., and Jochum. C. EROS — A Computer Program for Generating Sequences of Reactions. 74,93-126 (1978). 

---