Next Article in Journal
Carnosic Acid and Carnosol: Analytical Methods for Their Determination in Plants, Foods and Biological Samples
Next Article in Special Issue
Properties of Glass-Ceramics Prepared from Industrial Multi-Wastes
Previous Article in Journal
Chitosan and Metal Oxide Functionalized Chitosan as Efficient Sensors for Lead (II) Detection in Wastewater
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Calculation and 3-D Imaging of the Arrhenius Temperature Integral

1
The State Key Laboratory of Refractories and Metallurgy, Wuhan University of Science and Technology, Wuhan 430081, China
2
School of Chemical Engineering, University of New South Wales, Sydney, NSW 2052, Australia
3
Key Laboratory for Ferrous Metallurgy and Resources Utilization of Ministry of Education, Wuhan University of Science and Technology, Wuhan 430081, China
*
Author to whom correspondence should be addressed.
Separations 2023, 10(9), 480; https://doi.org/10.3390/separations10090480
Submission received: 28 July 2023 / Revised: 23 August 2023 / Accepted: 30 August 2023 / Published: 31 August 2023

Abstract

:
The Arrhenius temperature integral is typically used in non-isothermal kinetic analysis, which is widely applied in gas–solid reactions in separation processes. In previous studies, researchers provided various methods to solve the temperature integral, but the error usually became significant when the value of x (x = Ea/RT) was too large or too small. In this paper, we present a new series method and design a computer program to calculate the temperature integral. According to the precise calculation of the temperature integral, we first reveal the relationship among the integral, the temperature, and the activation energy, and we find an interesting phenomenon in which the 3-D image of the temperature integral is of self-similarity according to fractal theory. The work is useful for mechanism and theoretical studies of non-isothermal kinetics.

1. Introduction

Non-isothermal chemical kinetics are very important for gas–solid reactions in separation processes. As a result, they have been extensively introduced in disciplines, such as chemical engineering [1,2], mineral technology [3], metallurgical engineering [4], materials science [5], biomass energy [6,7], etc. The reaction rate of a solid-involved reaction is now generally accepted to be
d α d t = k f α
where α denotes the conversion ratio of the concerned substance at time t, f(α) is the mechanism function, and k is the rate constant expressed by the Arrhenius equation:
k = A exp E a / R T
where A, T, and Ea are the pre-exponential factor, temperature, and activation energy of the reaction, respectively.
For a linear heating process, the heating rate β is usually a function of time,
β = d T d t
Combining all of the above equations, re-arranging them, and integrating Equation (1), we obtain
g α = α 0 α d α f α = A β T 0 T exp E a R T d T
where T0 represents the initial reaction temperature and is set at 298 K in this study, and the initial reaction extent α0 usually equals zero.
The Arrhenius temperature integral is also known as the integral of Boltzmann factor Ψ(T), which is derived from Equation (4),
Ψ T = T 0 T exp E a R T d T
where another variable x could be introduced as
x = E a R T
Then, the temperature integral could be calculated as
Ψ T = E a R x E a R T 0 e x x 2   d x
where the upper limit of Equation (7) is usually considered to be ∞ regarding T0 in the denominator, which is equal to 0.
Then, the expression of the temperature integral as shown in Equation (7) is extensively regarded as the most concise and the best one for calculation.
The solution of the temperature integral is usually the basic demand for a non-isothermal kinetic study [8]. However, it cannot be analytically solved. With the solution of the temperature integral, constant activation energy and the corresponding kinetic mechanism could be determined in model-based (single-step) reactions, and variable activation energies could be calculated in model-free iso-conversional (multi-step) reactions.
The model-based non-isothermal analysis was dominant for a long time because of its ability to determine the kinetic triplet (pre-exponential factor, activation energy, and reaction model) [9,10]. However, these methods suffer from several problems, among which is their inability to uniquely determine the reaction model. This issue led to the decline of these methods in favor over the latest two decades of iso-conversional methods that evaluate kinetics without model assumptions [11]. Regardless of the kind of method, the evaluation of the temperature integral is an inevitable issue. As a result, there is a bottleneck in kinetic analysis, namely the calculation of the temperature integral.
The temperature integral was first encountered by scientists in the early thermogravimetric analyses [9,10,11,12,13]. Reich and Levi [14] made mistakes since they did not recognize that it is a transcendental function. Although different efforts have been made to solve this puzzle over the past six decades [15,16,17,18,19,20,21,22,23,24,25,26,27,28], troubles still occasionally arise since ‘there is no true value for the integral’, as stated by Tang [27]. The temperature integral has been playing an ‘enigmatic’ role in non-isothermal kinetic analysis. That is, it has appeared to be a necessary evil to be addressed [25].
Flynn [25] and Órfão [17] provided excellent reviews of the temperature integral. According to Flynn’s categorization, these solutions of the temperature integral are classified into three categories: series solutions, complex approximations [29,30], and simple approximations. It is given accurately for x greater than 2 by the rational fraction expression. Truncating the number of terms in the fraction expression, we can obtain one, two, three, and four “rational approximations”. Other series solutions, such as Schlömilch expansion, Bernouli number expansion, and Asymptotic expansion, are all available for the calculation limitation of x in different ranges [13,25].
In model-based non-isothermal analysis, simple and complex approximations were preferred and truncated from those series solutions because the truncated solutions could give the temperature integral an analytical expression and make the whole non-isothermal kinetics analytically solvable. For example, truncating the two-term asymptotic approximation, we can obtain the famous linear plot of ln(g(α)/T2) against 1/T, resulting in a slope of Ea/R. In the other words, it is exactly the famous Coats–Redfern method [12]. In addition to the truncations, some modified approximations were also proposed to improve the calculation precision [8,17,18,19,20,26,27,28].
The numerical calculation of the temperature integral is usually carried out by series methods, and model-free iso-conversional non-isothermal kinetics are usually required. Although the precision of the numerical calculation is better today, there are still backward calculations because of the limitations of the series method itself. Numerical calculation based on the quadrature method was occasionally mentioned in previous papers [24,28]. However, the residual error usually increases with the rise in the upper limit temperature of the integration, or it increases with the decrease in x [25]. As a disquieting situation, some published papers provided the wrong quadrature results as standard integral values, such as in Ref. [24].
There were hundreds of papers focused on the solution of the temperature integral once it was suggested that a relative error of 10% for the temperature integral could be acceptable for many calculations. However, more accurate values are still necessary for scientific purposes, such as distinguishing between reaction mechanisms [25]. Furthermore, a more accurate solution possible today with computer technology. All in all, accurate values and numerical calculation methods of the temperature integral are always attractive in non-isothermal kinetic analysis.

2. Theory

So far, most of the solutions for the temperature integral have been based on the concise Equation (7). However, the upper limit of the integration may be far from ∞, resulting in the corresponding inaccurate mathematical equations of series methods.
At the same time, all approximations of the temperature integral take a detour; errors and troubles cannot be eradicated. The only way to obtain the exact numerical solution of the integral is to develop an appropriate series solution, which easily reaches convergence. Here, we provide a new series solution based on Taylor’s expansion:
Ψ T 0 + Δ T = = n Ψ T 0 + Ψ T 0 Δ T + 1 2 ! Ψ T 0 Δ T 2 + + 1 n ! Ψ n T 0 Δ T n
where ΔT is the difference between T and T0, T = T0 + ΔT, and n denotes the order of Taylor expansion.
Then, the exact derivative value at temperature T0 could be calculated using Equations (7)–(9).
( n = 1 )             Ψ T = exp E a R T = g T
( n = 2 )             Ψ T = E a R × T 2 g T
( n >   2 )     Ψ n T = E a R × i = 0 n 1 C n 1 i g i T T 2 n i 1
The n-th-order derivative expression shown as Equation (9) is established according to Leibniz formula [31]. With the help of the general form of the n-th-order derivative equation, we designed a computer program to calculate the integral of the Boltzmann factor with Equation (6). The numerical value can converge within a sufficiently small residual error when the derivative order n is large enough, and ΔT is small enough. Computing the results shows that the value usually converges with reasonably little error when n is greater than 4 and ΔT is no greater than 50 K. Consequently, a multi-step integral method is applied to compute the integral value when the difference between the upper and lower limit temperature is larger than 50 K. The corresponding error analysis of the present calculation method is also provided below.

3. Error Analysis

Table 1 shows the converging evolution of the integral value with increasing order of Taylor expansion. Most of the values converge within a relative error of 1‰ when the Taylor order is larger than 4. A relative error within 10−8 can be achieved when the Taylor order increases to 10. To reach high accuracy, a 3D image (Figure 1a) of the integral function is established based on 200 × 200 calculation nodes, and each node is computed with a Taylor order of 20.
As shown in Figure 1a–d, the integral values have differences in orders of magnitude from each other with different temperatures and activation energies. To control computing error, the relative residual was defined as follows:
δ = Ψ n + 1 , T , E a Ψ n , T , E a Ψ n + 1 , T , E a
where Ψ(n, T, Ea) denotes the integral of Boltzmann factor with upper limit temperature of T, activation energy of Ea, and Taylor order of n. The relative residual δ denotes a tolerance when it is presupposed as a maximum limitation that should not be exceeded.
Figure 1a–d demonstrates the minimum required Taylor order for calculating the temperature integral under error tolerances of 10−4, 10−8, 10−12, and 10−16, respectively. According to Figure 1a,b, the relative residual errors in Figure 1a are no greater than 10−4 under the domain of definition, and most of the values have a residual error of 10−8. The smaller that the minimum required Taylor order is, the more easily that the iteration reaches tolerance. In previous studies, it is simply concluded that the errors become significant when the value of x (x = Ea/RT) is too small or too large [25]. However, according to these figures, a different trend is discovered. Generally, the trend that is proved by these four figures is that the minimum required Taylor order rises with a decreasing upper limit or an increase in the activation energy when the upper limit temperature is greater than about 350 K, and the minimum required Taylor order rises with increasing the upper limit or increase in the activation energy when the temperature is between 298 and 350 K.
Furthermore, the relationships among T, Ea, and Taylor order n are more complicated than monotonic, and some singular points also exist in the plots. However, the maximum required Taylor order is no greater than 45 under a tolerance of 10−16, which proves that it is easier to approach convergence than the former methods [25].
Another series solution is established to perform a comparison. Using the Newton–Leibniz quadrature method as a usual idea to resolve it, we obtain:
Ψ T = T 0 T exp E a R T d T = = n i = 0 n m × exp E a R T 0 + i m
where
m = T T 0 n
However, we find that the Newton–Leibniz quadrature method was difficult to converge, as shown in Table 2. Moreover, it is not realistic to adopt this method in kinetic analysis because it always requires a massive amount of time to reach high precision. As shown in Table 1, the results converge within 8 accurate digits after the decimal point when the Taylor order n is only larger than 15.
Table 3 provides specific expressions of the commonly used approximate solutions for the temperature integral. The values of temperature integral calculated by MATLAB were used to check the accuracy of other approximation methods. Table 4 reports the temperature integral values calculated by different methods. Except for the present solution, the other values are calculated based on a published article [29], in which the original values are equivalent to Ψ/[T·exp(−Ea/(RT))]. In Table 4, the temperature integral values of the present method and Órfão III’s approximation are fully consistent with MATLAB’s values. Therefore, the new series solution can provide sufficient, accurate values. In contrast, other approximations provide unreliable values when the value of x tends toward minority.

4. Results and Discussion

The dependence of the temperature integral on Ea and T has always been a mysterious relationship. Here, we show the true value of the integral in Figure 2a. This image is drawn based on 200 × 200 calculation nodes, and each node is computed with the Taylor order of 20 to ensure the numerical value is accurate enough. It can be seen that the integral value declines by an order of magnitude with the growth of the activation energy, proving that the integral value is extremely sensitive to the activation energy. Therefore, the approximation method often provides an inaccurate value and fails to distinguish the kinetic mechanisms [26]. The integral value tends to be equivalent to the value of ΔT when Ea infinitely approaches 0 kJ/mol, while the integral value declines by an order of magnitude and tends to be 0 when Ea is greater than about 75 kJ/mol.
There is also an interesting phenomenon in that the profile of the 3D surface is relatively similar to others if the independent variables (temperature, activation energy) are arbitrarily truncated, such as in images (b), (c), and (d). Image (a) is drawn based on a 200 × 200 node network, and each node is computed with the Taylor order of 20 to ensure the numerical value is accurate enough. Independent variables Ea and T vary from 2 to 200,002 J/mol and 300 to 1200 K, respectively. In the aforementioned node network (200 × 200), mid nodes (Ea varies from 50,002 to 149,002 J/mol, and T varies from 529.5 to 975 K) are chosen to form a new node network (100 × 100); then, image (b) can be drawn. The node network of image (c) is 50 × 50, and its values of Ea and T change from 20,002 to 69,002 J/mol and from 300 to 520.5 K, respectively. Regarding image (d), its node network is 49 × 63, and corresponding ranges of Ea and T are 52,002 to 100,002 J/mol and 651 to 930 K, respectively. As we can see, the truncated images b, c, and d are similar to the original image a. This phenomenon indicates the self-similarity property of the figure according to fractal theory [38].

5. Conclusions

Because of the demand for accurate solutions related to non-isothermal kinetics, previous researchers have provided various methods to solve the temperature integral, but the error usually becomes significant when the value of x (x = Ea/RT) is too small or too large. A numerical method is proposed to obtain an accurate value of the temperature integral. Multi-step and 20th Taylor expansion are applied to ensure the high accuracy of results calculated by the new method.
The results show that, when the upper limit temperature is greater than about 350 K, the error rises with the decreasing upper limit or the increasing activation energy; when the temperature is between 298 and 350 K, the error rises with the increasing upper limit or the increasing activation energy.
The temperature integral function is very sensitive to the activation energy according to the error analysis and could be well displayed by the 3-D image. The reason why the solution is always divergent when the value of x (x = Ea/RT) is too small or too large, as widely encountered in previous study, can be well interpreted by the contour image of the minimum Taylor order required when using the present method within the iterative tolerances. More interesting, the 3-D image of the temperature integral demonstrates self-similarity according to fractal theory.

6. Resource Codes

Here we share our resource codes for scientific purposes. Anyone who directly uses/translates these codes to produce new papers must cite this paper. All codes below are written in Visual Basic language.
These codes are easily transplanted into a project because they are written in functions. To use these codes, one should first define the environmental temperature (T0, default as 298 K), the presupposed Taylor expansion order (s), the Arrhenius energy (Ea), and the upper limit temperature (T1) of the integral. If the precision worsens, it only requires a smaller iteration temperature (default as 50 K in the code) and a larger Taylor expansion order. Usually, an iteration temperature of 50 K and a Taylor expansion order of 10 are good enough to reach a rational precision.
 
Option Explicit
Public T0 As Double’ T0 is the initial reaction temperature, default as 298K
Public i As Integer, k As Integer
Public ddg() As Double
Public store() As Double ‘Dynamic array for storage
Public store2() As Double
Public sum As Double
 
Function Psai(s As Integer, Ea, t1) ‘ Main function: calculation at an arbitrary temperature T1
Dim hun As Integer, jj As Integer
hun = Int((t1 − T0)/50)
ReDim Preserve store2(hun)
Dim sum0 As Double
sum0 = 0
If hun = 0 Then
Psai = gT(s, Ea, t1, T0)
Else
For jj = 1 To hun
store2(jj − 1) = gT(s, Ea, T0 + jj * 50, T0 + (jj − 1) * 50)
sum0 = sum0 + store2(jj − 1)
Next jj
Psai = sum0 + gT(s, Ea, t1, T0 + hun * 50)
End If
End Function
 
Function Factorial(n) ‘The n-th Factorial
If n = 0 Then
Factorial = 1
Else
Factorial = 1
For i = 1 To n
Factorial = Factorial * i
Next
End If
End Function
 
Function dg(n As Integer, Ea, T) As Double ‘ n-th derivative of g(T) for successive calculation.
Erase ddg()
ReDim ddg(n) ‘The n-th item in the polynomial expression of n-th-order derivative
ReDim Preserve store(n)
sum = 0
ddg(0) = Exp(−Ea/(8.314 * T))
store(0) = ddg(0)
If n > 0 Then
ddg(1) = Ea * Exp(−Ea/(8.314 * T)) * T ^ −2/8.314
store(1) = ddg(1)
End If
If n = 0 Then
dg = ddg(0)
ElseIf n = 1 Then dg = ddg(1)
Else
Dim m As Integer
m = n − 1
For k = 0 To m Step 1
ddg(k + 1) = (−1) ^ (m − k) * (m + 1 − k) * (Factorial(m)/Factorial(k)) * store(k) * T ^ (k − m − 2)
sum = sum + ddg(k + 1)
Next k
store(n) = sum * Ea/8.314
dg = store(n)
End If
End Function
 
 
Function dg_X(n As Integer, Ea, T) ‘ n-th-order derivative of g(T), unconditional
ReDim Preserve store(n)
For i = 0 To n
store(i) = dg(i, Ea, T)
Next i
dg_X = store(n)
End Function
 
Function gT(n As Integer, Ea, T, T00) ‘ Calculation for temperature difference smaller than 50 K
Dim j As Integer
sum = 0
For j = 1 To n
sum = sum + dg_X(j − 1, Ea, T00) * (T − T00) ^ j/Factorial(j)
Next j
gT = sum
End Function

7. Standard Integral Values

To be compared with the other values of temperature integrals, the standard integral values calculated by the present method are reported in Table 5 below.

Author Contributions

Conceptualization: W.Z.; writing—original draft: W.Z., Q.Z. and K.L.; writing and editing: X.Y. and Y.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Key Technologies R & D Program of Hubei Province, grant number 2022BCA058 and The APC was funded by the same funder.

Data Availability Statement

All data are available by contacting the corresponding author.

Acknowledgments

Many thanks to Chengzhi Li, Juhua Zhang, and Guangqiang Li for their suggestions and to Zhengliang Xue for his support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gautam, R.; Vinu, R. Unraveling the interactions in fast co-pyrolysis of microalgae model compounds via pyrolysis-GC/MS and pyrolysis-FTIR techniques. React. Chem. Eng. 2019, 4, 278–297. [Google Scholar] [CrossRef]
  2. Nestler, F.; Müller, V.P.; Ouda, M.; Hadrich, M.J.; Schaadt, A.; Bajohr, S.; Kolb, T. A novel approach for kinetic measurements in exothermic fixed bed reactors: Advancements in non-isothermal bed conditions demonstrated for methanol synthesis. React. Chem. Eng. 2021, 6, 1092–1107. [Google Scholar] [CrossRef]
  3. Zhang, J.; Zhang, W.; Xue, Z. Oxidation Kinetics of Vanadium Slag Roasting in the Presence of Calcium Oxide. Miner. Process. Extr. Metall. Rev. 2017, 38, 265–273. [Google Scholar] [CrossRef]
  4. Wang, H.H.; Li, G.Q.; Yang, J.; Ma, J.; Khan, B.S. The behavior of phosphorus during reduction and carburization of high-phosphorus oolitic hematite with H2 and CH4. Metall. Mater. Trans. B Process Metall. Mater. Process. Sci. 2016, 47, 2571–2581. [Google Scholar] [CrossRef]
  5. Jiang, Y.H.; Liu, F.; Song, S.J. An extended analytical model for solid-state phase transformation upon continuous heating and cooling processes: Application in γ/α transformation. Acta Mater. 2012, 60, 3815–3829. [Google Scholar] [CrossRef]
  6. Fermoso, J.; Gil, M.V.; Pevida, C.; Pis, J.J.; Rubiera, F. Kinetic models comparison for non-isothermal steam gasification of coal–biomass blend chars. Chem. Eng. J. 2010, 161, 276–284. [Google Scholar] [CrossRef]
  7. Tai, H.-C.; Li, G.-C.; Huang, S.-J.; Jhu, C.-R.; Chung, J.-H.; Wang, B.Y.; Hsu, C.-S.; Brandmair, B.; Chung, D.-T.; Chen, H.M.; et al. Chemical distinctions between Stradivari’s maple and modern tonewood. Proc. Natl. Acad. Sci. USA 2017, 114, 27–32. [Google Scholar] [CrossRef]
  8. Aghili, A. Representation and evaluation of the Arrhenius and general temperature integrals by special functions. Thermochim. Acta 2021, 705, 179034. [Google Scholar] [CrossRef]
  9. Zhang, W.; Li, K.; Dong, J.H.; Li, C.Z.; Liu, A.H.; Zhang, J.H.; Xue, Z.L. Kinetic triplet from low-temperature carburization and carbon deposition reactions. J. Iron Steel Res. Int. 2022, 29, 1545–1558. [Google Scholar] [CrossRef]
  10. Li, K.; Zhang, W.; Fu, M.; Li, C.; Xue, Z. Discussion on Criterion of Determination of the Kinetic Parameters of the Linear Heating Reactions. Minerals 2022, 12, 81. [Google Scholar] [CrossRef]
  11. Li, K.; Gan, C.; Zhang, W.; Li, C.; Li, G. Validity of isothermal kinetic prediction by advanced isoconversional method. Chem. Phys. 2023, 567, 111801. [Google Scholar] [CrossRef]
  12. Coats, A.W.; Redfern, J.P. Kinetic parameters from thermogravimetric data. Nature 1964, 201, 68–69. [Google Scholar] [CrossRef]
  13. Doyle, C.D. Series approximation to the equation of thermogravimetric data. Nature 1965, 207, 290–291. [Google Scholar] [CrossRef]
  14. Reich, L.; Levi, D. Thermal stability indices for polymeric materials based on energy considerations. Die Makromol. Chemie. 1963, 66, 102–113. [Google Scholar] [CrossRef]
  15. Doyle, C.D. Integral methods of kinetic analysis of thermogravimetric data. Macromol. Chem. Phys. 1964, 80, 220–224. [Google Scholar] [CrossRef]
  16. Zsakó, J. Kinetic analysis of thermogravimetric data, VI—Some problems of deriving kinetic parameters from TG curves. J. Therm. Anal. Calorim. 1973, 5, 239–251. [Google Scholar] [CrossRef]
  17. Órfão, J.J.M. Review and evaluation of the approximations to the temperature integral. AIChE J. 2007, 53, 2905–2915. [Google Scholar] [CrossRef]
  18. Neglur, R.; Grooff, D.; Hosten, E.; Aucamp, M.; Liebenberg, W. Approximation-based integral versus differential isoconversional approaches to the evaluation of kinetic parameters from thermogravimetry: Kinetic analysis of the dehydration of a pharmaceutical hydrate. J. Therm. Anal. Calorim. 2016, 123, 2599–2610. [Google Scholar] [CrossRef]
  19. Han, J.; Liu, D.; Qin, L.; Chen, W.; Xing, F. A modified temperature integral approximation formula and its application in pyrolysis kinetic parameters of waste tire. Energy Sources Part A Recover. Util. Environ. Eff. 2018, 40, 220–226. [Google Scholar] [CrossRef]
  20. Lin, Y.; Tian, Y.; Xia, Y.; Fang, S.; Liao, Y.; Yu, Z.; Ma, X. General distributed activation energy model (G-DAEM) on co-pyrolysis kinetics of bagasse and sewage sludge. Bioresour. Technol. 2019, 273, 545–555. [Google Scholar] [CrossRef]
  21. Senum, G.I.; Yang, R.T. Rational approximations of the integral of the Arrhenius function. J. Therm. Anal. Calorim. 1977, 11, 445–447. [Google Scholar] [CrossRef]
  22. Zsakó, J. Kinetic analysis of thermogravimetric data—XIV. Three integral methods and their computer programs. J. Therm. Anal. Calorim. 1980, 19, 333–345. [Google Scholar] [CrossRef]
  23. Serra, R.; Nomen, R.; Sempere, J. The Non-Parametric Kinetics A New Method for the Kinetic Study of Thermoanalytical Data. J. Therm. Anal. Calorim. 1998, 52, 933–943. [Google Scholar] [CrossRef]
  24. Rao, V.K.; Bardon, M.F. Integral of the boltzmann factor-A new approximation. J. Therm. Anal. 1996, 46, 323–326. [Google Scholar] [CrossRef]
  25. Flynn, J.H. The ‘Temperature Integral’—Its use and abuse. Thermochim. Acta 1997, 300, 83–92. [Google Scholar] [CrossRef]
  26. Vyazovkin, S. Modification of the integral isoconversional method to account for variation in the activation energy. J. Comput. Chem. 2001, 22, 178–183. [Google Scholar] [CrossRef]
  27. Tang, W.; Liu, Y.; Zhang, H.; Wang, C. New approximate formula for Arrhenius temperature integral. Thermochim. Acta 2003, 408, 39–43. [Google Scholar] [CrossRef]
  28. Capela, J.M.V.; Capela, M.V.; Ribeiro, C.A. Approximations for the generalized temperature integral: A method based on quadrature rules. J. Therm. Anal. Calorim. 2009, 97, 521–524. [Google Scholar] [CrossRef]
  29. Budrugeac, P. Applicability of some approximations of the temperature integral used for heating processes to processes taking place on cooling. J. Therm. Anal. Calorim. 2016, 124, 479–485. [Google Scholar] [CrossRef]
  30. Chen, H.X.; Liu, N.A. Approximations for the temperature integral. J. Therm. Anal. Calorim. 2008, 92, 573–578. [Google Scholar] [CrossRef]
  31. Chu, W. Derivative inverse series relations and lagrange expansion formula. Int. J. Number Theory 2013, 9, 1001–1013. [Google Scholar] [CrossRef]
  32. Doyle, C.D. Estimating isothermal life from thermogravimetric data. J. Appl. Polym. Sci. 1962, 6, 639–642. [Google Scholar] [CrossRef]
  33. Gorbachev, V.M. A solution of the exponential integral in the non-isothermal kinetics for linear heating. J. Therm. Anal. Calorim. 1975, 8, 349–350. [Google Scholar] [CrossRef]
  34. Pérez-Maqueda, L.A.; Criado, J.M. The Accuracy of Senum and Yang’s Approximations to the Arrhenius Integral. J. Therm. Anal. Calorim. 2000, 60, 909–915. [Google Scholar] [CrossRef]
  35. Urbanovici, E.; Segal, E. Some problems concerning the temperature integral in non-isothermal kinetics: Part I. Generalities and some simple applications. Thermochim. Acta 1990, 168, 71–87. [Google Scholar] [CrossRef]
  36. Chen, H.; Liu, N. New procedure for derivation of approximations for temperature integral. AIChE J. 2006, 52, 4181–4185. [Google Scholar] [CrossRef]
  37. Ji, L.Q. New rational fraction approximating formulas for the temperature integral. J. Therm. Anal. Calorim. 2008, 91, 885–889. [Google Scholar] [CrossRef]
  38. Mandelbrot, B. How long is the coast of Britian? Statistical self-similarity and fractal dimension. Science 1967, 156, 636–638. [Google Scholar] [CrossRef]
Figure 1. The minimum Taylor order required when using present method within the iterative tolerances of (a) 10−4, (b) 10−8, (c) 10−12, and (d) 10−16.
Figure 1. The minimum Taylor order required when using present method within the iterative tolerances of (a) 10−4, (b) 10−8, (c) 10−12, and (d) 10−16.
Separations 10 00480 g001aSeparations 10 00480 g001b
Figure 2. Three-dimensional image of the integral function of Boltzmann factor.
Figure 2. Three-dimensional image of the integral function of Boltzmann factor.
Separations 10 00480 g002aSeparations 10 00480 g002bSeparations 10 00480 g002c
Table 1. The integral values of Boltzmann factor with varying Taylor order n.
Table 1. The integral values of Boltzmann factor with varying Taylor order n.
Taylor Order nEa = 2 J/molEa = 200 J/molEa = 20,000 J/molEa = 200,000 J/mol
T = 600 KT = 1200 KT = 600 KT = 1200 KT = 600 KT = 1200 KT = 600 KT = 1200 K
1301.82115578 901.64929838 284.65221331 867.71805582 1.47187736 41.45574855 1.11339239 × 10−176.64568260 × 10−8
2301.83262904 901.66613166 285.72844879 869.31399874 1.84315996 44.66857854 2.03501118 × 10−179.66472837 × 10−8
3301.83160096 901.66487181 285.63545060 869.19891467 1.88712698 44.76038220 3.15684847 × 10−171.05477139 × 10−7
4301.83170869 901.66499125 285.64482423 869.20938372 1.88871043 44.75981330 4.17892561 × 10−171.07227989 × 10−7
5301.83169624 901.66497815 285.64378407 869.20828368 1.88857776 44.75969041 4.88644170 × 10−171.07474026 × 10−7
δ4.12482 × 10−81.5 × 10−83.6 × 10−61.3 × 10−67 × 10−52.7 × 10−60.1447920.002289
6301.83169778 901.66497973 285.64390756 869.20841067 1.88857039 44.75968614 5.27169484 × 10−171.07498862 × 10−7
7301.83169758 901.66497953 285.64389214 869.20839504 1.88857151 44.75968706 5.44047391 × 10−171.07500633 × 10−7
8301.83169761 901.66497955 285.64389414 869.20839705 1.88857153 44.75968709 5.50078992 × 10−171.07500716 × 10−7
9301.83169761 901.66497955 285.64389388 869.20839679 1.88857152 44.75968708 5.51849690 × 10−171.07500718 × 10−7
10301.83169761 901.66497955 285.64389391 869.20839682 1.88857152 44.75968708 5.52277114 × 10−171.07500718 × 10−7
δ001.05 × 10−103.45 × 10−11000.0007740
11301.83169761 901.66497955 285.64389391 869.20839682 1.88857152 44.75968708 5.52361451 × 10−171.07500718 × 10−7
12301.83169761 901.66497955 285.64389391 869.20839682 1.88857152 44.75968708 5.52374846 × 10−171.07500718 × 10−7
13301.83169761 901.66497955 285.64389391 869.20839682 1.88857152 44.75968708 5.52376503 × 10−171.07500718 × 10−7
14301.83169761 901.66497955 285.64389391 869.20839682 1.88857152 44.75968708 5.52376652 × 10−171.07500718 × 10−7
15301.83169761 901.66497955 285.64389391 869.20839682 1.88857152 44.75968708 5.52376660 × 10−171.07500718 × 10−7
δ0000001.45 × 10−80
16301.83169761 901.66497955 285.64389391 869.20839682 1.88857152 44.759687085.52376660 × 10−171.07500718 × 10−7
17301.83169761 901.66497955 285.64389391 869.20839682 1.88857152 44.759687085.52376660 × 10−171.07500718 × 10−7
18301.83169761 901.66497955 285.64389391 869.20839682 1.88857152 44.759687085.52376660 × 10−171.07500718 × 10−7
19301.83169761 901.66497955 285.64389391 869.20839682 1.88857152 44.759687085.52376660 × 10−171.07500718 × 10−7
20301.83169761 901.66497955 285.64389391 869.20839682 1.88857152 44.759687085.52376660 × 10−171.07500718 × 10−7
δ00000000
Table 2. Results of the Newton–Leibniz quadrature method.
Table 2. Results of the Newton–Leibniz quadrature method.
Quadrature Node mEa = 2 J/molEa = 200 J/molEa = 20,000 J/molEa = 200,000 J/mol
T = 600 KT = 1200 KT = 600 KT = 1200 KT = 600 KT = 1200 KT = 600 KT = 1200 K
1 × 102304.84987233 910.68041913 288.48730649 877.78754503 1.91652897 45.37009283 6.12793756 × 10−171.16596876 × 10−7
1 × 103302.13351522 902.56652506 285.92824787 870.06645411 1.89135955 44.82059549 5.58241547 × 10−171.08390284 × 10−7
δ−0.008990585−0.00899−0.00895−0.00887−0.01331−0.01226−0.09772−0.07571
1 × 104301.86187937 901.75513412 285.67232943 869.29420397 1.88885024 44.76577660 5.52961379 × 10−171.07589474 × 10−7
δ−0.000899059−0.0009−0.0009−0.00089−0.00133−0.00122−0.00946−0.00739
1 × 105301.83471578 901.67399501 285.64673746 869.21697755 1.88859939 44.76029602 5.52435114 × 10−171.07509592 × 10−7
δ−8.99949 × 10−5−9 × 10−5−9 × 10−5−8.9 × 10−5−0.00013−0.00012−0.00095−0.00074
1 × 106301.83199942 901.66588110 285.64417826 869.20925489 1.88857430 44.75974797 5.52382505 × 10−171.07501606 × 10−7
δ−8.99949 × 10−6−9 × 10−6−9 × 10−6−8.9 × 10−6−1.3 × 10−5−1.2 × 10−5−9.5 × 10−5−7.4 × 10−5
Table 3. Expressions of the approximation for temperature integral.
Table 3. Expressions of the approximation for temperature integral.
ModelApproximationP(x)
DDoyle [32] exp ( 5.331 1.052 x )
CRCoats and Redfern [12] exp ( x ) x 2 1 2 x
G1Gorbachev (1st degree) [33] exp ( x ) x 1 x + 2
SY4Senum and Yang (4th degree) [21] exp ( x ) x x 3 + 18 x 2 + 86 x + 96 x 4 + 20 x 3 + 120 x 2 + 240 x + 120
PC8Perez-Maqueda and Criado (8th degree) [34] exp ( x ) x x 7 + 70 x 6 + 1886 x 5 + 24920 x 4 + 170136 x 3 + 577584 x 2 + 844560 x + 357120 x 8 + 72 x 7 + 2024 x 6 + 28560 x 5 + 216720 x 4 + 880320 x 3 + 1794240 x 2 + 1572480 x + 403200
US4Urbanovici and Segal IV [35] exp ( x ) x 2 ( 2 + 2 x 1 + 8 x + 4 x 2 )
CLChen and Liu [36] exp ( x ) x 2 3 x 2 + 16 x + 4 3 x 2 + 12 x + 30
O3Órfão III [17] exp ( x ) x 0.9999936 x 3 + 7.5739391 x 2 + 12.4648922 x + 3.6907232 x 4 + 9.5733223 x 3 + 25.6329561 x 2 + 21.099653 x + 3.95849
J3Ji III [37] exp ( x ) x 2 x 3 + 9.27052 x 2 + 16.79440 x + 1.20025 x 3 + 11.27052 x 2 + 33.33602 x + 24.21457
Table 4. The temperature integral value calculated by the present method and the other approximation methods.
Table 4. The temperature integral value calculated by the present method and the other approximation methods.
X = Ea/(RT) (T = 600 K)MTALABPresent
Solution (n = 10)
Numerical QuadratureDCRG1SY4PC8US4CLO3J3
17.8106 × 1017.8106 × 1018.9108 × 1016.6486 × 10−1−2.2086 × 1026.3661 × 1017.7575 × 1017.8076 × 1017.6156 × 1019.8401 × 1017.8106 × 1017.8362 × 101
22.1597 × 1012.1597 × 1012.2517 × 1016.2429 × 10−1−6.6416 × 1011.9419 × 1012.1575 × 1012.1597 × 1012.1429 × 1012.8311 × 1012.1597 × 1012.1601 × 101
41.9093 × 1001.9093 × 1001.9187 × 1001.7034 × 10−11.3648 × 1001.8221 × 1001.9092 × 1001.9093 × 1001.9059 × 1002.5173 × 1001.9093 × 1001.9094 × 100
61.9083 × 10−11.9083 × 10−11.9096 × 10−13.1556 × 10−21.6513 × 10−11.8579 × 10−11.9083 × 10−11.9083 × 10−11.9071 × 10−12.4537 × 10−11.9083 × 10−11.9083 × 10−1
82.0481 × 10−22.0481 × 10−22.0490 × 10−25.1394 × 10−31.8868 × 10−22.0126 × 10−22.0481 × 10−22.0481 × 10−22.0475 × 10−22.5632 × 10−22.0481 × 10−22.0481 × 10−2
102.2981 × 10−32.2981 × 10−32.2982 × 10−37.8368 × 10−42.1792 × 10−32.2700 × 10−32.2981 × 10−32.2981 × 10−32.2978 × 10−32.8087 × 10−32.2981 × 10−32.2981 × 10−3
122.6575 × 10−42.6575 × 10−42.6576 × 10−41.1470 × 10−42.5601 × 10−42.6332 × 10−42.6575 × 10−42.6575 × 10−42.6572 × 10−43.1836 × 10−42.6575 × 10−42.6575 × 10−4
143.1404 × 10−53.1404 × 10−53.1402 × 10−51.6322 × 10−53.0546 × 10−53.1182 × 10−53.1404 × 10−53.1404 × 10−53.1402 × 10−53.6997 × 10−53.1404 × 10−53.1404 × 10−5
163.7724 × 10−63.7724 × 10−63.7724 × 10−62.2751 × 10−63.6926 × 10−63.7512 × 10−63.7724 × 10−63.7724 × 10−63.7723 × 10−64.3821 × 10−63.7724 × 10−63.7724 × 10−6
184.5901 × 10−74.5901 × 10−74.5900 × 10−73.1218 × 10−74.5126 × 10−74.5690 × 10−74.5901 × 10−74.5901 × 10−74.5900 × 10−75.2684 × 10−74.5901 × 10−74.5901 × 10−7
205.6429 × 10−85.6429 × 10−85.6430 × 10−84.2306 × 10−85.5651 × 10−85.6213 × 10−85.6429 × 10−85.6429 × 10−85.6428 × 10−86.4106 × 10−85.6429 × 10−85.6429 × 10−8
Table 5. Standard values of temperature integral calculated by present method.
Table 5. Standard values of temperature integral calculated by present method.
T/KEa = 2 J/molEa = 20,002 J/molEa = 40,002 J/molEa = 60,002 J/molEa = 80,002 J/molEa = 100,002 J/molEa = 120,002 J/molEa = 140,002 J/molEa = 160,002 J/molEa = 180,002 J/molEa = 200,002 J/mol
3001.9984 × 1006.4068 × 10−42.0545 × 10−76.5897 × 10−112.1142 × 10−146.7844 × 10−182.1777 × 10−216.9916 × 10−252.2453 × 10−287.2121 × 10−322.3172 × 10−35
3181.9984 × 1018.1688 × 10−33.4100 × 10−61.4524 × 10−96.3035 × 10−132.7834 × 10−161.2483 × 10−195.6766 × 10−232.6130 × 10−261.2157 × 10−295.7090 × 10−33
3363.7971 × 1011.9709 × 10−21.0915 × 10−56.3972 × 10−93.9296 × 10−122.5059 × 10−151.6455 × 10−181.1055 × 10−217.5617 × 10−255.2465 × 10−283.6823 × 10−31
3545.5959 × 1013.6632 × 10−22.7010 × 10−52.1870 × 10−81.8958 × 10−111.7247 × 10−141.6239 × 10−171.5678 × 10−201.5425 × 10−231.5400 × 10−261.5556 × 10−29
3727.3947 × 1016.0526 × 10−25.9036 × 10−56.5170 × 10−87.8001 × 10−119.8417 × 10−141.2870 × 10−161.7265 × 10−192.3604 × 10−223.2746 × 10−254.5961 × 10−28
3909.1935 × 1019.3187 × 10−21.1877 × 10−41.7523 × 10−72.8220 × 10−104.7989 × 10−138.4604 × 10−161.5301 × 10−182.8199 × 10−215.2736 × 10−249.9776 × 10−27
4081.0992 × 10021.3659 × 10−12.2414 × 10−44.3258 × 10−79.1447 × 10−102.0423 × 10−124.7284 × 10−151.1229 × 10−172.7174 × 10−206.6724 × 10−231.6574 × 10−25
4261.2791 × 10021.9288 × 10−14.0116 × 10−49.9213 × 10−72.6921 × 10−97.7171 × 10−122.2931 × 10−146.9882 × 10−172.1700 × 10−196.8371 × 10−222.1792 × 10−24
4441.4590 × 10022.6432 × 10−16.8611 × 10−42.1336 × 10−67.2838 × 10−92.6265 × 10−119.8160 × 10−143.7621 × 10−161.4691 × 10−185.8205 × 10−212.3327 × 10−23
4621.6389 × 10023.5330 × 10−11.1278 × 10−34.3345 × 10−61.8290 × 10−88.1498 × 10−113.7631 × 10−131.7817 × 10−158.5949 × 10−184.2063 × 10−202.0823 × 10−22
4801.8189 × 10024.6228 × 10−11.7901 × 10−38.3718 × 10−64.2978 × 10−82.3293 × 10−101.3079 × 10−127.5303 × 10−154.4167 × 10−172.6281 × 10−191.5817 × 10−21
4981.9988 × 10025.9379 × 10−12.7540 × 10−31.5456 × 10−59.5179 × 10−86.1862 × 10−104.1652 × 10−122.8751 × 10−142.0216 × 10−161.4420 × 10−181.0404 × 10−20
5162.1787 × 10027.5038 × 10−14.1200 × 10−32.7401 × 10−51.9987 × 10−71.5384 × 10−91.2263 × 10−111.0021 × 10−138.3418 × 10−167.0436 × 10−186.0155 × 10−20
5342.3586 × 10029.3461 × 10−16.0104 × 10−34.6837 × 10−54.0012 × 10−73.6055 × 10−93.3646 × 10−113.2183 × 10−133.1354 × 10−153.0985 × 10−173.0969 × 10−19
5522.5385 × 10021.1490 × 1008.5707 × 10−37.7461 × 10−57.6707 × 10−78.0102 × 10−98.6610 × 10−119.5978 × 10−131.0833 × 10−141.2401 × 10−161.4358 × 10−18
5702.7184 × 10021.3962 × 1001.1972 × 10−21.2433 × 10−41.4140 × 10−61.6953 × 10−82.1043 × 10−102.6767 × 10−123.4675 × 10−144.5560 × 10−166.0542 × 10−18
5882.8984 × 10021.6786 × 1001.6410 × 10−21.9419 × 10−42.5151 × 10−63.4332 × 10−84.8509 × 10−107.0233 × 10−121.0355 × 10−131.5485 × 10−152.3418 × 10−17
6063.0783 × 10021.9987 × 1002.2112 × 10−22.9587 × 10−44.3305 × 10−66.6785 × 10−81.0659 × 10−91.7432 × 10−112.9029 × 10−134.9027 × 10−158.3734 × 10−17
6243.2582 × 10022.3588 × 1002.9328 × 10−24.4062 × 10−47.2372 × 10−61.2522 × 10−72.2419 × 10−94.1121 × 10−117.6803 × 10−131.4547 × 10−142.7864 × 10−16
6423.4382 × 10022.7613 × 1003.8340 × 10−26.4261 × 10−41.1769 × 10−52.2698 × 10−74.5293 × 10−99.2588 × 10−111.9271 × 10−124.0673 × 10−148.6810 × 10−16
6603.6181 × 10023.2084 × 1004.9459 × 10−29.1936 × 10−41.8663 × 10−53.9888 × 10−78.8190 × 10−91.9973 × 10−104.6052 × 10−121.0768 × 10−132.5458 × 10−15
6783.7980 × 10023.7022 × 1006.3022 × 10−21.2921 × 10−32.8918 × 10−56.8119 × 10−71.6597 × 10−84.1417 × 10−101.0522 × 10−112.7106 × 10−137.0607 × 10−15
6963.9780 × 10024.2448 × 1007.9394 × 10−21.7866 × 10−34.3858 × 10−51.1330 × 10−63.0267 × 10−88.2811 × 10−102.3065 × 10−116.5139 × 10−131.8601 × 10−14
7144.1579 × 10024.8382 × 1009.8971 × 10−22.4329 × 10−36.5210 × 10−51.8388 × 10−65.3615 × 10−81.6009 × 10−94.8659 × 10−111.4995 × 10−124.6725 × 10−14
7324.3378 × 10025.4841 × 1001.2217 × 10−13.2666 × 10−39.5191 × 10−52.9175 × 10−69.2446 × 10−82.9995 × 10−99.9066 × 10−113.3172 × 10−121.1231 × 10−13
7504.5178 × 10026.1844 × 1001.4944 × 10−14.3290 × 10−31.3660 × 10−44.5323 × 10−61.5545 × 10−75.4593 × 10−91.9514 × 10−107.0716 × 10−122.5910 × 10−13
7684.6977 × 10026.9407 × 1001.8124 × 10−15.6671 × 10−31.9292 × 10−46.9041 × 10−62.5538 × 10−79.6713 × 10−93.7276 × 10−101.4565 × 10−115.7539 × 10−13
7864.8777 × 10027.7547 × 1002.1808 × 10−17.3347 × 10−32.6845 × 10−41.0326 × 10−54.1051 × 10−71.6707 × 10−86.9194 × 10−102.9052 × 10−111.2332 × 10−12
8045.0576 × 10028.6277 × 1002.6045 × 10−19.3922 × 10−33.6840 × 10−41.5183 × 10−56.4661 × 10−72.8188 × 10−81.2505 × 10−95.6239 × 10−112.5568 × 10−12
8225.2376 × 10029.5613 × 1003.0890 × 10−11.1908 × 10−24.9903 × 10−42.1970 × 10−59.9934 × 10−74.6526 × 10−82.2043 × 10−91.0586 × 10−105.1393 × 10−12
8405.4175 × 10021.0557 × 1013.6398 × 10−11.4956 × 10−26.6780 × 10−43.1317 × 10−51.5172 × 10−67.5228 × 10−83.7955 × 10−91.9410 × 10−101.0035 × 10−11
8585.5975 × 10021.1615 × 1014.2624 × 10−11.8620 × 10−28.8349 × 10−44.4018 × 10−52.2653 × 10−61.1931 × 10−76.3933 × 10−93.4725 × 10−101.9066 × 10−11
8765.7774 × 10021.2737 × 1014.9628 × 10−12.2991 × 10−21.1564 × 10−36.1058 × 10−53.3297 × 10−61.8581 × 10−71.0549 × 10−86.0707 × 10−103.5313 × 10−11
8945.9574 × 10021.3925 × 1015.7468 × 10−12.8168 × 10−21.4983 × 10−38.3647 × 10−54.8224 × 10−62.8447 × 10−71.7073 × 10−81.0385 × 10−96.3849 × 10−11
9126.1373 × 10021.5179 × 1016.6205 × 10−13.4258 × 10−21.9229 × 10−31.1326 × 10−46.8879 × 10−64.2858 × 10−72.7130 × 10−81.7405 × 10−91.1286 × 10−10
9306.3173 × 10021.6500 × 1017.5901 × 10−14.1377 × 10−22.4458 × 10−31.5167 × 10−49.7100 × 10−66.3598 × 10−74.2375 × 10−82.8613 × 10−91.9529 × 10−10
9486.4972 × 10021.7888 × 1018.6618 × 10−14.9650 × 10−23.0845 × 10−32.0100 × 10−41.3520 × 10−59.3033 × 10−76.5120 × 10−84.6192 × 10−93.3118 × 10−10
9666.6772 × 10021.9345 × 1019.8419 × 10−15.9209 × 10−23.8590 × 10−32.6375 × 10−41.8606 × 10−51.3426 × 10−69.8548 × 10−87.3300 × 10−95.5104 × 10−10
9846.8571 × 10022.0872 × 1011.1137 × 1007.0196 × 10−24.7913 × 10−33.4288 × 10−42.5324 × 10−51.9130 × 10−61.4698 × 10−71.1444 × 10−89.0052 × 10−10
10027.0371 × 10022.2468 × 1011.2553 × 1008.2759 × 10−25.9062 × 10−34.4183 × 10−43.4108 × 10−52.6929 × 10−62.1623 × 10−71.7594 × 10−81.4468 × 10−9
10207.2170 × 10022.4134 × 1011.4096 × 1009.7058 × 10−27.2309 × 10−35.6459 × 10−44.5484 × 10−53.7473 × 10−63.1398 × 10−72.6657 × 10−82.2872 × 10−9
10387.3970 × 10022.5872 × 1011.5774 × 1001.1326 × 10−18.7955 × 10−37.1572 × 10−46.0085 × 10−55.1581 × 10−64.5032 × 10−73.9834 × 10−83.5609 × 10−9
10567.5770 × 10022.7680 × 1011.7592 × 1001.3153 × 10−11.0633 × 10−29.0047 × 10−47.8665 × 10−57.0268 × 10−66.3830 × 10−75.8746 × 10−85.4638 × 10−9
10747.7569 × 10022.9561 × 1011.9556 × 1001.5206 × 10−11.2779 × 10−21.1248 × 10−31.0212 × 10−49.4789 × 10−68.9471 × 10−78.5562 × 10−88.2686 × 10−9
10927.9369 × 10023.1513 × 1012.1674 × 1001.7504 × 10−11.5272 × 10−21.3954 × 10−31.3149 × 10−41.2667 × 10−51.2409 × 10−61.2315 × 10−71.2350 × 10−8
11108.1168 × 10023.3537 × 1012.3952 × 1002.0067 × 10−11.8156 × 10−21.7199 × 10−31.6801 × 10−41.6778 × 10−51.7036 × 10−61.7524 × 10−71.8216 × 10−8
11288.2968 × 10023.5634 × 1012.6395 × 1002.2914 × 10−12.1474 × 10−22.1067 × 10−32.1310 × 10−42.2035 × 10−52.3165 × 10−62.4671 × 10−72.6550 × 10−8
11468.4768 × 10023.7803 × 1012.9010 × 1002.6067 × 10−12.5276 × 10−22.5652 × 10−32.6840 × 10−42.8704 × 10−53.1210 × 10−63.4376 × 10−73.8259 × 10−8
11648.6567 × 10024.0045 × 1013.1804 × 1002.9548 × 10−12.9615 × 10−23.1059 × 10−33.3578 × 10−43.7104 × 10−54.1682 × 10−64.7432 × 10−75.4538 × 10−8
11828.8367 × 10024.2360 × 1013.4782 × 1003.3380 × 10−13.4544 × 10−23.7402 × 10−34.1741 × 10−44.7609 × 10−55.5203 × 10−66.4836 × 10−77.6941 × 10−8
12009.0166 × 10024.4748 × 1013.7950 × 1003.7584 × 10−14.0124 × 10−24.4808 × 10−35.1571 × 10−46.0659 × 10−57.2526 × 10−68.7835 × 10−71.0748 × 10−7
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, W.; Zheng, Q.; Yu, X.; Shen, Y.; Li, K. Numerical Calculation and 3-D Imaging of the Arrhenius Temperature Integral. Separations 2023, 10, 480. https://doi.org/10.3390/separations10090480

AMA Style

Zhang W, Zheng Q, Yu X, Shen Y, Li K. Numerical Calculation and 3-D Imaging of the Arrhenius Temperature Integral. Separations. 2023; 10(9):480. https://doi.org/10.3390/separations10090480

Chicago/Turabian Style

Zhang, Wei, Qiaoyu Zheng, Xiaobing Yu, Yansong Shen, and Kui Li. 2023. "Numerical Calculation and 3-D Imaging of the Arrhenius Temperature Integral" Separations 10, no. 9: 480. https://doi.org/10.3390/separations10090480

APA Style

Zhang, W., Zheng, Q., Yu, X., Shen, Y., & Li, K. (2023). Numerical Calculation and 3-D Imaging of the Arrhenius Temperature Integral. Separations, 10(9), 480. https://doi.org/10.3390/separations10090480

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop