National Institute for Research and Development in Electrical Engineering ICPE-CA

Sigla ICPE-CA MINISTRY OF RESEARCH AND INNOVATION Sigla MEC romana english



Advanced methodology for the kinetic analysis of complex heterogeneous processes with application in prediction of thermal behavior of materials and their thermal lifetime,
Acronym: MET-AV


Project Registration Code: PN-III-P4-ID-PCE-2016-0088

Contract number: 112PCE/2017

Financing: Public budget

Programme: Program 4 - Fundamental and Border Research

Research domain: 4 - Eco-Nano-Technology and Advanced Materials

Project type: Exploratory Research Projects(PCE)

Total amount of the contract: 850.000 lei
Of which, by financing sources:
Source 1 – public budget: 850.000 lei
Source 2 – own budget: 0 lei

Contract duration: 12.07.2017 - 31.12.2018)

Contracting Authority: Executive Agency for Higher Education, Research, Development and Innovation Funding (UEFISCDI)

Contractor: NATIONAL INSTITUTE FOR RESEARCH AND DEVELOPMENT IN ELECTRICAL ENGINEERING ICPE-CA BUCHAREST (INCDIE ICPE-CA)

Partner P1: NATIONAL INSTITUTE FOR RESEARCH AND DEVELOPMENT IN MICROTECHNOLOGIES IMT BUCHAREST (IMT)

Abstract:
Research teams:
Project objectives:
Project stages:
Expected results:
Contact:
Stage I:
Stage II


Abstract:
The critical analysis of the actual used methods of kinetic analysis of heterogeneous processes by using thermal analysis data (thermogravimetry, differential thermal analysis, differential scanning calorimetry) and some physical-chemical methods (IR spectroscopy, X-Ray diffraction, etc.) will be performed, putting in evidence the following open problems that will be tried solving: improving the accuracy of isoconversional methods by suggestion of an original method, assessment of errors in the evaluation of activation energy by model-free nonlinear methods, elaboration of a general algorithm for evaluation of the kinetic scheme (mechanism) and corresponding kinetic parameters for a complex heterogeneous process. The theoretical results will be checked for the following experimental data: non-isothermal and isothermal data corresponding to decomposition and/or thermo-oxidation of some polymers and polymeric materials, and the thermal decomposition of calcium carbonate. The general algorithm for evaluation of kinetic scheme and corresponding kinetic parameters of decomposition or/and thermo-oxidation of a material will be used for improving the procedures of prediction of thermal lifetime of materials that is an important parameter in designing of devices and equipments, or the prediction of duration of decomposition in industrial conditions of a raw material in order to obtain a compound or intermediate in chemical synthesis. A procedure for rapid prediction of thermal lifetime of a material by using thermal analysis methods and based on kinetic analysis of heterogeneous processes, will be elaborated. In comparison with Standardized method (IEC 216), this procedure has the following advantages: it could be applied for materials that exhibit of single-step thermal deterioration process, and also to materials that exhibit a complex thermal deterioration process; a less time consuming.

Research teams:
Coordinator CO - INCDIE ICPE-CA Bucharest
1. PhD Chim Budrugeac Petru – IDT I, Project manager

2. PhD Chim. Cucoş Andrei – CS II, Key person
3. PhD. Chim. Ştefănescu Carmen – CSII, Research team member
4. PhD Ing. Sbârcea Beatrice Gabriela - Research team member
5. Ing. Mitrea Sorina Mitrea – IDTI, Research team member
6. Ing. Chiose Ileana Laura - Research team member
7. Masterand Dascălu Radu-Cristian - Research team member

Project objectives:
The main objectives of the project are:
a. The use of thermal analysis methods (thermogravimetry (TG or TGA), differential thermal analysis (DTA), differential scanning calorimetry (DSC)) and some physical-chemical methods (IR spectroscopy (FTIR), X-Ray diffraction (XRD), etc.) for investigation of non-isothermal decomposition and/or thermo-oxidative degradation of solid compounds or materials
b. The improving of the methods for kinetic analysis of complex heterogeneous processes by suggestion and checking of an original isoconversional method, solving the problem of assessment of errors in activation energy evaluated by isoconversional non-linear methods, elaboration and checking for simulated and experimental data of a general algorithm (methodology) for evaluation of kinetic scheme and corresponding kinetic parameters of decomposition or/and thermo-oxidation of a material.
c. The use the general algorithm (methodology) for evaluation of kinetic scheme and corresponding kinetic parameters of decomposition or/and thermo-oxidation of a material for improving the procedures of prediction of thermal lifetime of materials or the prediction of decomposition duration in industrial conditions of a raw material in order to obtain a compound or intermediate in chemical synthesis.
d. Elaboration of a procedure for rapid prediction of thermal lifetime of a material by using thermal analysis methods, which is also applicable for prediction of the duration of decomposition of a raw material.


Project stages:
Stage I
Stage name: Methods of kinetic analysis of experimental data obtained by thermal analysis of materials characterized structurally by physical-chemical methods (IR spectroscopy, X-ray diffraction, etc.)
Activity type: Fundamental research
Stage duration: 12.07.2017 – 15.12.2017

Stage II
Stage name: DSSC achievement. Solve a problem regarding non-linear isoconversional methods and elaboration of a general algorithm for kinetic analysis of non-isothermal and isothermal data corresponding to heterogeneous processes

Activity type: Fundamental research
Partners involved: CO, P1
Stage duration: 16.12.2018 – 14.12.2018

Stage III
Stage name: Improving the procedures for prediction of the thermal life of the materials and the duration of decomposition of a raw material under industrial conditions
Activity type: Fundamental research
Stage duration: 15.12.2018 – 31.12.2019


Expected results:
• Documentary study;
• Physical-chemical characterization of materials, including the determination of their thermal and thermo-oxidative stabilities
• The improving of methods of kinetic analysis of heterogeneous processes (processes with participation of a solid compound or material), including the elaboration and checking of a general algorithm (methodology) for evaluation of kinetic scheme and corresponding kinetic parameters of decomposition or/and thermo-oxidation of a material ;
• The use of the general algorithm (methodology) for evaluation of kinetic scheme and corresponding kinetic parameters of decomposition or/and thermo-oxidation of a material for improving the procedures of prediction of thermal lifetime of materials or the prediction of duration of decomposition in industrial conditions of a raw material
• Scientific works communicated at national and international scientific events;
• Scientific papers submitted for publication in ISI rated journals.


Contact:
Coordinator CO: INCDIE ICPE-CA Bucharest
Project manager: PhD Petru BUDRUGEAC

Phone: 0040213467231 / extension 129, e-mail: petru.budrugeac@icpe-ca.ro


Stage I

Abstract

A critical documentary study concerning the methods of kinetic analysis of thermo-analytical data corresponding to complex heterogeneous processes has been presented. The open problems concerning the assessment of kinetic scheme and corresponding kinetic parameters have been put in evidence.
It has been suggested a simple and precise isoconversional method applicable for evaluation of activation energy dependence on the conversion degree. For the validation of the method the data obtained by thermal analysis of calcium carbonate, PVC and HDPE were used. There was a very good concordance between the activation energy values determined by the simple method proposed with those obtained by the "advanced isoconversion" and "incremental iterative" methods, which are much more laborious.
For calcium carbonate, high density polyethylene (HDPE) and vinyl polychloride (PVC), structural characterization was determined by X-ray diffraction and X-ray photoelectron spectroscopy (XPS).
Pentru materialele carbonatului de calciu, polietilenei de înaltă densitate (HDPE) şi policlorurii de vinil (PVC) s-au determinat caracterizarea structurală prin difracţie de raze X şi spectroscopiei fotoelectronice de raze X (XPS). Also, for HDPE and PVC, TG / DTG + DTA or DSC simultaneous thermal analyses coupled with FTIR were performed both under inert atmosphere and oxidizing atmosphere, in order to determine the volatile products as result of progressive heating. The data obtained were used for the validation of the proposed isoconversion method, but will be used also in the next stages of the present project which have as main objectives the elaboration of a general algorithm of kinetic analysis of the complex heterogeneous processes and the elaboration of a procedure for rapid determination of the lifetime thermal properties of materials and the duration of heterogeneous processes.


Dissemination of the results

• Submition of the following paper to Thermochimica Acta (ISI journal): „A simple and precise differential incremental isoconversional method to kinetic analysis of heterogeneous processes under arbitrary temperature programs”, author Petru Budrugeac
• Oral communication „General algorithm for evaluation of the kinetic scheme and corresponding kinetic parameters of heterogeneous processes. Applications.” author Petru Budrugeac, presented at international conference “4th Central and Eastern European Conference on Thermal Analysis and Calorimetry (CEEC-TAC4)” , Chişinău – Republica Moldova, 28 August 2017 - 31 August 2017
• Poster „Using of TG-FTIR method for studying the thermal decomposition of some electro-insulating materials”, authors Andrei Cucos, Petru Budrugeac, presented at international conference “4th Central and Eastern European Conference on Thermal Analysis and Calorimetry (CEEC-TAC4)” , Chişinău – Republica Moldova, 28 August 2017 - 31 August 2017
• Invited conference “Non-isothermal model-free prediction for assessment of conversion vs. time curves for complex processes under arbitrary temperature programs – advantages and limitations”, author P. Budrugeac”, presented at international conference “13th Mediterranean Conference on Calorimetry and Thermal Analysis (Medicta 2017)”, Loano – Italy, 24 September 2017 – 27 September2017



Stage II

Stage name: Solve a problem regarding non-linear isoconversional methods and elaboration of a general algorithm for kinetic analysis of non-isothermal and isothermal data corresponding to heterogeneous processes
Activity type: Fundamental research
Stage duration: 16.12.2018 – 14.12.2018
Activity 2.1. Solving the problem of error determination of activation energy determination by non-linear isoconversion methods
An original method for evaluation of the error in activation energy determinated by non-linear isoconversional (model-free) methods was elaborated and checked for some non-isothermal data. For this purpose, several linear and nonlinear isoconversional methods were applied for non-isothermal data corresponding to crystallization of (GeS2)0.3(Sb2S3)0.7 (4 constant heating rates), decomposition of ammonium perchlorate (6 constant heating rates), decomposition of poly(vinyl chloride) (PVC) (5 constant heating rates) and simulated data (12 constant heating rates). For each set of non-isothermal data, the linear and nonlinear isoconversional methods were applied for evaluation of activation energy (E). It has been considered some pairs “linear isoconversional method + nonlinear isoconversional method” (“differential isoconversional method suggested by Friedman + nonlinear differential method”; each “integral pair” corresponds to a certain approximation of the temperature integral). The comparison of the errors of activation energy) evaluated by a linear method (ΔLE) with the Fisher confidence intervals evaluated by the corresponding nonlinear method ΔFE shows that the relation ΔLE = b x ΔFE + c x (ΔFE)2 (b and c are parameters that not depend on the pair of isoconversional method) exhibits a very good accuracy of ΔLE vs. ΔFE fitting when ΔFE is evaluated for confidence level of 95%. Consequently, the errors in E evaluation by nonlinear methods could be determined by using the above equation and ΔFE values determined for the confidence level of 95%. This procedure was also checked for non-isothermal data corresponding to thermal decomposition of HDPE, which were not used for evaluation of b and c parameters.
Activity 2.2. Determination of thermal degradation curves of a kind of epoxy resin at minimum 3 constant heating rates
A “composite material based on epoxy resin (D010S+D110S cured with phtalic anhydride)” was characterized by X-ray fluorescence spectroscopy (XRF), X-ray diffraction (XRD), and differential scanning calorimetry (DSC). The thermal behavior of this material, used as electro-insulating material, was performed by the following thermal analysis methods: differential scanning calorimetry (DSC), and simultaneous thermogravimetry (TG) + differential thermogravitry (DTG) + differential thermal analysis (DTA).
The heating curves TG + DTG + DTA were recorded using STA 490C apparatus produced by Netzsch – Germany, in the following conditions: temperature range 25oC – 1000oC, heating rates of 2.99 K.min-1, 4.98 K.min-1, 10.04 K.min-1, 15.42 K.min-1, and 21.80 K.min-1, nitrogen flow (30 ml.min-1), and Pt-Rh crucible. It was obtained that the material exhibits two successive global processes characterized by mass losses, and peaks in DTG and DTA curves. The parameters of these processes depend on the heating rates.
Activity 2.3. Determination of thermal degradation isotherms of a kind of epoxy resin.
Four quasi-isothermal temperature programs (quasi-365; quasi-370; quasi-375 and quasi-380, where 365, 370, 375 and 380 are the temperatures (oC) of isothermal part of temperature programs) in which the thermal degradations of composite material based on epoxy resin were determined by analyzing the non-isothermal data obtained in step 2.2. For each temperature program, the dependence of mass loss vs. time was recorded.
Activity 2.4. Applying the non-linear regression method and some physical-chemical analysis methods for finding the kinetic scheme and the kinetic parameters corresponding to the heterogeneous process of decomposition of a kind of epoxy resin
The applied kinetic analysis algorithm of data obtained by thermal analysis involves the following successive steps:
(1) application for data recorded at steady heating rates of isoconversional methods to determine the dependence of apparent activation energy on the degree of conversion;
(2) processing data recorded at constant heating rates by using a non-linear regression program (we used ‘‘Netzsch Thermokinetics” program) to determine the kinetic scheme and the corresponding kinetic parameters;
(3) verification of the results obtained in step (2) for the "conversion degree vs. time or temperature” curves recorded at certain temperature programs, other than those used for kinetic parameter evaluation.
The first two steps of this procedure (algorithm) was applied for all non-isothermal data recorded at constant heating rates and whole process of thermal degradation of “composite material based on epoxy resin (D010S+D110S cured with phtalic anhydride)”. It was obtained that this process is complex one and consists in five successive reactions. For each step, the values of corresponding kinetic parameters were determined. In the third step of algorithm, the kinetic scheme and kinetic parameters thus determined were used to calculate the quasi-isothermal degradation TG curves that were compared to the corresponding experimental curves. It follows that there is a satisfactory concordance between calculated and experimental TG curves. Thus, the mechanism of the thermal degradation process and the corresponding kinetic parameters can be used for predictions corresponding to arbitrary temperature programs.
Activity 2.5. Dissemination of results through scientific communications presented at Conferences and / or Symposia and by sending articles for publication in ISI journals

I. Published paper
1. Petru Budrugeac, A simple and precise differential incremental isoconversional method to kinetic analysis of heterogeneous processes under arbitrary temperature programs, Thermochimica Acta, 661 (2018) 116-123, doi.org/10.1016/j.tca.2018.01.025
Abstract

A simple and precise linear differential incremental isoconversional (model-free) method for kinetic analysis of heterogeneous processes using thermo-analytical data recorded at arbitrary temperature programs has been suggested. This method has been applied to simulated data corresponding to a single-step process and to a complex process consisting in two consecutive reactions, and experimental data corresponding to thermal degradation of high density polyethylene (HDPE). The obtained results are consistent with those from other isoconversional methods (advanced isoconversional method suggested by Vyazovkin, iterative method suggested by Budrugeac and differential isoconversional method suggested by Friedman).

2. P.Budrugeac, Estimating errors in the determination of activation energy by nonlinear methods applied for thermoanalytical measurements performed under constant heating rates, Thermochimica Acta, 670 (2018) 1 - 6, doi.org/10.1016/j.tca.2018.09.020
Abstract

Several linear and nonlinear isoconversional methods have been applied for following non-isothermal thermoanalytical data: simulated data for two consecutive first order reactions (12 heating rates), crystallization of (GeS2)0.3(Sb2S3)0.7 (4 heating rates), decomposition of ammonium perchlorate (6 heating rates) and decomposition of poly(vinyl chloride) (PVC) (5 heating rates). It has been considered some pairs “linear isoconversional method + nonlinear isoconversional method”. The “differential pair” is “differential isoconversional method suggested by Friedman + nonlinear differential method”, while each “integral pair” corresponds to a certain approximation of the temperature integral. The values of activation energy (E), error of E obtained by linear method (ΔLE) and Fischer confidence interval obtained by nonlinear method ΔFE applying the procedure suggested by Vyazovkin and Wight have been determined for each pair of methods, several conversion degrees, and the confidence levels of 68.27%, 80%, 90% and 95%. It has been obtained that, for a certain pair of methods, (a) ΔFE values are substantially greater than ΔLE values, and (b) the values of E determined by linear method are identical with those determined by the nonlinear method. The statement (a) is explained by the procedure for (ΔFE) evaluation in which it is assumed that ΔFE correspond to maximum value of Fischer distribution function. According statement (b) it is expected that is a relationship between ΔLE and ΔFE. Both statements suggest that the error in E determined by a nonlinear isoconversional method is equal with ΔLE. Satisfactory fittings of ΔLE vs. ΔFE have be obtained for the relationships: (1) ΔLE = a x ΔFE and (2) ΔLE = b x ΔFE + c x (ΔFE)2, here a, b and c are parameters which depend on the confidence limit. These relations have been also checked for high density polyethylene (HDPE) decomposition data that were not used for their derivations. For all considered data, the best accuracy of fitting of ΔLE vs. ΔFE has been obtained for equation (2) and ΔFE determined for confidence level of 95%. It has been conclude that the evaluation of error in E determined by a nonlinear isoconversional method involves the following two successive steps: the determination of ΔFE for confidence level of 95%, and the application of relation

II. Communications presented at at “12th European Symposium on Thermal Analysis and Calorimetry (ESTAC12)”, which happened from 27th to 30th of August 2018, in Brasov –Romania
1. Petru Budrugeac, Critical study concerning the use of sinusoidal modulated thermogravimetric data (MTG) for evaluation of activation energy of heterogeneous processes. (Oral presentation)
2. Andrei Cucos, Petru Budrugeac, Thermal behaviour of some electro-insulating materials studied by coupled TG-FTIR technique (Poster)
3. Carmen Paraschiv, Gabriela Hristea, Marius Lungulescu, Beatrice-Gabriela Sbarcea, Virgil Marinescu, Investigation of thermal stability and morphology of zinc oxide-graphene oxide heterostructures obtained by hydrothermal synthesis (Poster)
4. Radu-Cristian Dascălu, Andrei Cucoş, Cătălin Maxim, Delia-Laura Popescu, Marius Andruh, Comparative Thermal Studies of Polynuclear Coordination Compounds with Organostannic Knots (Poster)

III. Communication presented at „PRIOCHEM - Priorităţile Chimiei pentru o Dezvoltare Durabilă”, 10 – 12 October, Bucureşti.
5. Radu-Cristian Dascălu, Andrei Cucoş, Cătălin Maxim, Delia-Laura Popescu, Marius Andruh, Synthesis, characterization and thermal behaviour of polynuclear coordination compounds with organostannic knots (Poster)
Each published article and communication is accompanied by the following Acknowledgements: The work was supported by the Romanian “Ministry of Research and Innovation –Executive Agency for Higher Education, Research, Development and Innovation Funding, UEFISCDI” research project MET-AV, PN-III-P4-ID-PCE Nr. 112/2017.