1.4: Kinetic Analysis Based on Thermogravimetry Data
- Page ID
- 222824
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)4.1 Dehydration of Calcium Oxalate Monohydrate
Let us consider the dehydration as the reaction that occurs by the scheme
, that is,
Figure 4.1 shows the TG curves of the dehydration of calcium oxalate monohydrate.

The CaC2O4 · H2O dehydration was studied on a Netzsch TG 209 F 3 Tarsus thermo-microbalance. The experiment was run at three heating rates: 5, 7.5, and 10 K/min. Three measurements were taken at each heating rate, other conditions being identical. Standard aluminum crucibles without lids were used as holders. The process was carried out in a dry air flow at a rate of 200 mL/min. The initial reagent was freshly precipitated calcium oxalate with a particle size of 15–20 μm. Weighed portions of the reagent were 5–6 mg for each heating rate.
4.2 Computational Procedure. Solution of the Inverse and Direct Kinetic Problems. Quasi-One-Stage Process
As follows from Figure 4.1, the dehydration in the given temperature range can be considered a quasi-one-stage reaction at all heating rates used. The experimental data obtained on NETZSCH equipment are processed with the NETZSCH Proteus software.
For further work with the experimental data using the NETZSCH Thermokinetics software, it is necessary to export data from the Proteus program in the tabulated form (measured signal as a function of temperature or time) as an ASCII file. To do this, the user must select the desired curve and click the Extras Export data button in the Proteus toolbar. The user must enter the lower and upper limits of the data range to be exported. To correctly specify the limits, the derivative of the selected curve is used. The left- and right-hand limits are chosen in the ranges where the derivative becomes zero (Figure 4.2).

Remember that the derivative of the selected curve can be obtained by clicking the corresponding icon in the NETZSCH Proteus program window.











4.3 Analysis of Computation Results
Let us consider the computation results obtained by the linear regression method for the CaC2O4 · H2O dehydration (Figure 4.19).
Figure 4.19 presents the Arrhenius parameters and the form and characteristics of the function best fitting the experimental results (from the statistical viewpoint). For the reaction under consideration, the best fitting function is the Prout–Tompkins equation with autocatalysis (the Bna code), which is indicated at the top left of the table. However, before discussing the meaning of the results obtained, let us consider the F-test: Fit Quality window (Figure 4.20).

4.3.1 The F-test: Fit quality and F-test
Step significance windows present the statistical analysis of the fit quality for different models. This allows us to determine using the statistical methods which of the models provides the best fit for the experimental data.
To perform such an analysis, Fisher’s exact test is used. In general, Fisher’s test is a variance ratio which makes it possible to verify whether the difference between two independent estimates of the variance of some data samples is significant. To do this, the ratio of these two variances is compared with the corresponding tabulated value of the Fisher distribution for a given number of degrees of freedom and significance level. If the ratio of two variances exceeds the corresponding theoretical Fisher test value, the difference between the variances is significant.
In the Thermokinetics software, Fisher’s test is used for comparing the fit qualities ensured by different models. The best-fit model, that is, the model with the minimal sum of squared deviations, is taken as a reference (conventionally denoted as model 1). Then, each model is compared to the reference model. If the Fisher test value does not exceed the critical value, the difference between current model 2 and reference model 1 is insignificant. There is no reason to then believe that model 1 provides a more adequate description of the experiment in comparison to model 2.
The Fexp value is estimated by means of Fisher’s test:
\[F_{e x p}=\frac{L S Q_{1} / f_{1}}{L S Q_{2} / f_{2}} \label{4.2}\]
The Fexp value is compared with the Fisher distribution Fcrit(0.95) for the significance level of 0.95 and the corresponding number of degrees of freedom.
Figure 4.21).

In the Const. column, the option ‘false’ is set for the parameters that should be varied and the option ‘true’ is chosen for the parameters that remain constant. The three columns to the right of this column are intended for imposing constraints on the selected values. The computation results are presented in Figure 4.25.


The following conclusions can be drawn from Table 4.1: first, the autocatalysis parameters for the Bna and CnB functions are almost zero, that is, all is reduced to the Fn function. Second, the error of the Arrhenius parameters for Fn are minimal. Hence, the calcium oxalate dehydration is better described by the nth-order function. The reaction order can be considered to be 1/3, that is, the process is described by the “contracting sphere equation.” This means that the sample consists of spherical particles of the same size and that the dehydration is a homothetic process, that is, the particles during decomposition undergo a self-similar decrease in size. Such a mechanism is inherent in thermolysis of inorganic crystal hydrates. Thus, the problem of CaC2O4·H2O dehydration macrokinetics can be thought of as solved.
Function code | log A | E, kJ/mol | Reaction order n | log K cat 1 | Exp a1 |
---|---|---|---|---|---|
Bna | 6.8±0.6 | 73±6 | 0.34±0.25 | — | 0.06±0.14 |
CnB | 6.8±1 | 74±8 | 0.41±0.48 | -0.45±1.9 | — |
Fn | 7.0±0.1 | 75.2±0.8 | 0.34±0.25 | — | — |
The project created in the NETZSCH Thermokinetics software is saved by clicking on the common button (Figure 4.27).

