Home

Staff
     
Consulting
               
Mathematica Tools

Customer Quotes
        
Contact

Older Mathematica Items

 

 

 

 

 

Scientific Arts

Differential Equations for Chemical Kinetics


Example of Differential Equations for Chemical Kinetics

Chemical reactions of a wide variety can be modeled with coupled (often nonlinear) differential equations. These describe the time evolution of the concentrations of the various chemical species: reactants, intermediaries, catalysts, and products. Such problems are quite simple to set up and solve with Mathematica. The function NDSolve can numerically integrate the differential equations that arise. The resulting concentrations can be plotted as a function of time and also be used to accurately compute the expected concentration of the molecular species.

Reaction with an intermediate

Here is the reaction, [Graphics:Images/chemicalkinetics_gr_1.gif] where the mechanism of the reaction proceeds through an intermediate species [Graphics:Images/chemicalkinetics_gr_2.gif].

[Graphics:Images/chemicalkinetics_gr_3.gif]

The rate equations corresponding to this reaction are

[Graphics:Images/chemicalkinetics_gr_4.gif]
[Graphics:Images/chemicalkinetics_gr_5.gif]
[Graphics:Images/chemicalkinetics_gr_6.gif]

where [Graphics:Images/chemicalkinetics_gr_7.gif], [Graphics:Images/chemicalkinetics_gr_8.gif], and [Graphics:Images/chemicalkinetics_gr_9.gif] represent the concentrations of the species [Graphics:Images/chemicalkinetics_gr_10.gif], [Graphics:Images/chemicalkinetics_gr_11.gif], and [Graphics:Images/chemicalkinetics_gr_12.gif] at the time [Graphics:Images/chemicalkinetics_gr_13.gif]. The parameters [Graphics:Images/chemicalkinetics_gr_14.gif], [Graphics:Images/chemicalkinetics_gr_15.gif], and [Graphics:Images/chemicalkinetics_gr_16.gif] (the "rate constants") are considered to be constant here. This might be the case if the reaction proceeds in contact with a heat bath. If the temperature of the mixture of reactants is allowed to be variable, then there needs to be an additional set of equations supplied which model the temperature variation of [Graphics:Images/chemicalkinetics_gr_17.gif], [Graphics:Images/chemicalkinetics_gr_18.gif], and [Graphics:Images/chemicalkinetics_gr_19.gif] with time.

Although these equations do not have an explicit analytic solution, they can be numerically integrated using the built in numerical differential equation solver, NDSolve. Here is the solution for a particular choice of the rate constants [Graphics:Images/chemicalkinetics_gr_20.gif], [Graphics:Images/chemicalkinetics_gr_21.gif], and [Graphics:Images/chemicalkinetics_gr_22.gif]. In this example the initial concentrations of [Graphics:Images/chemicalkinetics_gr_23.gif] and [Graphics:Images/chemicalkinetics_gr_24.gif] are equal and that of the intermediary [Graphics:Images/chemicalkinetics_gr_25.gif] is zero.

[Graphics:Images/chemicalkinetics_gr_26.gif]
[Graphics:Images/chemicalkinetics_gr_27.gif]

The result is expressed in terms of a set of replacement rules that give the functions as InterpolatingFunction objects. The notation [Graphics:Images/chemicalkinetics_gr_28.gif], for example, is a shorthand for the numerical information needed by the interpolation algorithm to reproduce the solution.

Here is plot of the result (we multiply the concentration of [Graphics:Images/chemicalkinetics_gr_29.gif] by a factor of 50 so that it can be seen in the graph on the same scale as [Graphics:Images/chemicalkinetics_gr_30.gif] and [Graphics:Images/chemicalkinetics_gr_31.gif]). The colors of the curves for [Graphics:Images/chemicalkinetics_gr_32.gif], [Graphics:Images/chemicalkinetics_gr_33.gif], and [Graphics:Images/chemicalkinetics_gr_34.gif], are black, red, and green respectively.

[Graphics:Images/chemicalkinetics_gr_35.gif]

[Graphics:Images/chemicalkinetics_gr_36.gif]

For this example, as the reaction proceeds, the intermediary is created and then settles down to an approximately steady state value.

Reaction with Catalysts

[Graphics:Images/chemicalkinetics_gr_37.gif]

The rate equations corresponding to this reaction are

[Graphics:Images/chemicalkinetics_gr_38.gif]
[Graphics:Images/chemicalkinetics_gr_39.gif]
[Graphics:Images/chemicalkinetics_gr_40.gif]
[Graphics:Images/chemicalkinetics_gr_41.gif]

Here is a particular case of the integration of these differential equations. The initial concentrations of [Graphics:Images/chemicalkinetics_gr_42.gif] and [Graphics:Images/chemicalkinetics_gr_43.gif] are equal and the rate constants [Graphics:Images/chemicalkinetics_gr_44.gif] and [Graphics:Images/chemicalkinetics_gr_45.gif]chosen as shown.

[Graphics:Images/chemicalkinetics_gr_46.gif]
[Graphics:Images/chemicalkinetics_gr_47.gif]

Here is plot of the result. The colors of the curves for [Graphics:Images/chemicalkinetics_gr_48.gif], [Graphics:Images/chemicalkinetics_gr_49.gif], [Graphics:Images/chemicalkinetics_gr_50.gif], and [Graphics:Images/chemicalkinetics_gr_51.gif] are black, red, green, and blue respectively.

[Graphics:Images/chemicalkinetics_gr_52.gif]

[Graphics:Images/chemicalkinetics_gr_53.gif]

It can be seen that, as the reaction proceeds towards completion, the concentrations of the catalysts [Graphics:Images/chemicalkinetics_gr_54.gif] and [Graphics:Images/chemicalkinetics_gr_55.gif] eventually return to their initial values.


For further information on our services send email to info@scientificarts.com .
Contents of this web site Copyright © 1999-2011, Scientific Arts, LLC.

s