Abstract

The objects of this investigation were the comparative kinetic analysis of the NiO and NiCl2 reduction by hydrogen during an induction period and elimination of the calcination during the synthesis of Ni/Al2O3 catalysts. The effect of temperature and time on NiO and NiCl2 reduction degrees was studied. Avrami I equation was selected as the most favorable kinetic model and used to determine activation energy of the NiO and NiCl2 reduction for the investigated temperature range (623–923 K) and time intervals (1–5 minutes). The investigation enabled reaching conclusions about the reaction ability and rate of the reduction processes. Afterward, Ni/Al2O3 catalysts were obtained by using oxide and chloride precursor for Ni. The catalysts were supported on alumina-based foam and prepared via aerosol route. Properties of the samples before and after low-temperature hydrogen reduction (633 K) were compared. Obtained results indicated that the synthesis of Ni/Al2O3 catalysts can be more efficient if chloride precursor for Ni is directly reduced by hydrogen during the synthesis process, without the calcination step. In addition, Ni-Pd/Al2O3 catalysts with different metal content were prepared by using chloride precursors. Lower reduction temperature was utilized and the chlorides were almost completely reduced at 533 K.

1. Introduction

Nowadays, nickel is one of the most technologically important metals due to its physical and chemical properties. It is widely used as the alloying element, particularly in demanding corrosion-resistant applications [13]. The production of fine Ni powders, which are used in electronics industry and for the polycomponent catalysts manufacturing, is a great challenge in the field of new materials [46]. A widely used method for Ni powder production is reduction of NiO by hydrogen:

Two main characteristics in the kinetics of metal oxide reduction are the existence of an induction period and the possibility for autocatalysis [7, 8]. On the basis of both theoretical and experimental analyses of the NiO reduction by hydrogen, Delmon [9] determined number of potential Ni nuclei and the other parameters which confirmed the law of the new phase nuclei formation. Eliminating the influence of the diffusion controlled stage, Delmon reduced the spherical NiO particles (diameter from 0.7 to 0.8 μm). The reduction was carried out in a temperature range from 406 to 498 K, with both the homogeneous and the heterogeneous nucleation of Ni on the added foreign particles. The foreign particles, as the Ni crystallization nuclei, caused positive effect on the NiO reduction rate. The induction period of reduction process, which occurred with the homogeneous nucleation, lasted 2.2 minutes, while in the process occurring with addition of foreign nuclei, it lasted 1.6 minutes. Richardson et al. [10] studied mechanism of the reduction of porous bulk NiO at elevated temperatures (448–573 K). Their results show a direct transformation without accumulation of any intermediate phase. Without addition of water to the reducing gas (H2), reduction had the following stages: (1) an induction period, where the initial NiO reduction occurred and Ni clusters were formed, (2) increase of the reduction rate with the enlargement of Ni clusters, and (3) a pseudo-first-order (excess H2) process, where NiO disappeared, Ni was formed instead, and the reduction slowed at a fractional conversion of about 0.8 [10]. A number of research groups greatly contributed to investigation of kinetics and mechanism of the NiO reduction by hydrogen. Activation energy values for the NiO reduction were 20 and 19 kJ/mol for powder and pellet, respectively [11], and 90 kJ/mol for Goro NiO [12]. Plascencia and Utigard [7] determined activation energy for the reduction of two commercial nickel oxides. The values were 32 and 30 kJ/mol for Tokyo and Sinter 75 sample, respectively. Janković [13] reduced NiO prepared by sol-gel method. Experimental results were compared with 15 theoretical reaction models. He proposed sequential-reaction mechanism and average apparent activation energy value was 101 kJ/mol. Considering other metal oxides, Piotrowski et al. [14] obtained activation energy values of 28 and 123 kJ/mol when Fe2O3 was reduced by H2 and CO, respectively. According to that research group, the processes were initially chemically controlled [14]. According to Go et al. [15], the water splitting reaction with Fe2O3, ZnFe2O4, and MnFe2O4 was controlled by the product-layer diffusion. Activation energy values ranged from 57 to 110 kJ/mol.

Another method for the synthesis of Ni powder is reduction of NiCl2 by hydrogen:

In comparison to the NiO hydrogen reduction, the kinetics of the NiCl2 reduction by the same reducing agent has been the subject of a smaller number of studies. Suh et al. [16] investigated the kinetics of gas phase reduction of NiCl2. The aim was to produce Ni particles in hydrogen atmosphere using a tubular furnace reactor. They obtained the activation energy value of 104 kJ/mol. Stopić et al. [17] conducted kinetic analysis of NiCl2 reduction, with 0.1 wt% of Pd, Cu, or Ni and without additives. The activation energy of NiCl2 reduction, at temperatures from 573 to 823 K, was 54 kJ/mol. The presence of additives caused a decrease in the activation energy at lower temperatures (533 to 653 K): 33 kJ/mol for Ni(Pd); 50 kJ/mol for Ni(Cu) and Ni(Ni). The influence of Pd was explained by dissociation and hydrogen spillover, where reactive molecular hydrogen was formed and intensified the reduction process. Good solubility of Cu in Ni decreased energy bonds in the Ni lattice and added Ni(HCOO)2 formed artificial Ni nuclei. Effect of Pd to a decrease in particle size and improved uniformity of particle shape was the most dominant. Stopić et al. [18] synthesized fine nanosized Ni powder by ultrasonic spray and hydrogen reduction pyrolysis. Addition of Pd and Cu enabled complete NiCl2 reduction at 1173 K in dynamic conditions. When Pd was added, ideally spherical nonagglomerated Ni particles were obtained. The presented method was proven to be efficient for the production of fine nanocrystalline Fe-Ni alloy particles from corresponding metal chlorides [19]. Pérez-Hernández et al. [20] successfully produced Ni-based catalysts by impregnation and they used NiCl2 as a precursor for Ni. That preparation method included calcination at 773 K for 2 h in air and reduction by hydrogen for 2 h at 773 K [20]. However, the calcination may not be necessary during the synthesis of catalysts. In the earlier work of authors [21], addition of Pd, Cu, and Fe positively influenced the reduction process of NiCl2 in hydrogen atmosphere, at low temperatures (533 to 653 K). The highest reduction degree of 58 wt% was achieved with Pd addition at 653 K for 24 minutes. It was concluded that obtained results can particularly be significant for the synthesis of Ni-based catalysts. Reduction of chloride precursors can enable economic and technological benefits in the catalysts production, by elimination of the calcination step and reduction at low temperatures. According to Juan-Juan et al. [22], the calcination (773 K or 973 K for 12 h) can be avoided. Ni/Al2O3 catalysts prepared by impregnation and obtained by the direct hydrogen reduction of nickel salt at 773 K or 973 K for 2 h had good performance in reforming of CH4 with CO2. These authors noted that the calcination caused formation of hardly reducible NiAl2O4 phase in the catalysts [22]. It is known that modification with other active components, especially noble metals (Pd, Pt), improves catalytic properties of Ni-based catalysts [23, 24].

Catalytic activity can be enhanced if active components are supported on reticulated ceramic foams [25]. Using of alumina-based foams, synthesized as described in the previous research of authors [26], can contribute to energy savings in the complete process of catalyst manufacturing. The foams with improved mechanical properties were prepared by polymer replication. In comparison to the current synthesis processes, lower temperature was utilized during the sintering step (1673 K) [26]. Authors developed synthesis method for Ni-based catalysts supported on alumina-based foams in the earlier work [27]. That investigation included dependence of the catalyst properties from synthesis method and the addition of Pd, Cu, and Fe. The foams were impregnated with ultrasonically aerosolized solutions of metal chlorides. Oxides and chlorides on the foam surface were obtained by calcination and drying, respectively. Then, the samples were reduced by hydrogen at low temperatures (533 and 633 K). Ni/Al2O3 catalyst prepared by using chloride precursor had Ni coating without agglomerates and cracks, while those obtained from oxide precursor had residual agglomerates. All of the samples with chlorides reached higher reduction degree than those with oxides. The catalyst with 0.1 wt% of Pd, prepared from chlorides, was nearly completely reduced at 533 K [27].

This paper aims to compare kinetics of NiO and NiCl2 reduction by hydrogen during an induction period and to consider the possibility of direct reduction of chloride precursor for Ni during the catalysts synthesis. First, kinetic analysis was made in order to reach conclusions about reduction mechanism. By testing kinetic equations for solid state reactions [28], it was concluded that the experimental results were best fitted by Avrami I equation. Activation energy of NiO and NiCl2 reduction was determined from slopes of corresponding Arrhenius plots and by using Arrhenius equation. Afterwards, Ni/Al2O3 catalysts with 20 wt% Ni loading, supported on alumina-based foam, were prepared by using aerosol route, as previously described [27]. Ni coating was obtained from oxide and chloride precursor. For comparison purposes, the following catalysts were also synthesized: Ni/Al2O3 with 5 wt% Ni loading and Ni-Pd/Al2O3 with 5 wt% and 20 wt% Ni loading. For synthesis of these catalysts, only chloride precursors for Ni and Pd were used. Properties of the samples before and after hydrogen reduction were examined. Elimination of the calcination step and the reduction at low temperature were proposed.

2. Materials and Methods

For the reduction of NiO and NiCl2, the used materials included NiO powder and (MERCK, pro analysis). The NiCl2 powder was prepared by dehydration of for 1 h at 473 K and grinding in a mortar. The experimental investigations of NiO and NiCl2 reduction were conducted in a tubular quartz reactor set in an electric resistance furnace. Silica vessel carriers contained 10 g of each powder. The samples were reduced in static conditions, from 1 to 5 minutes at temperatures from 623 to 923 K (hydrogen flow rate 18 L/h).

Apparatus for the reduction experiments is illustrated in Figure 1 and presented in the previous research of authors [21]. Hydrogen was supplied from the high-pressure container linked to the columns system for purification (the silica gel, the calcium chloride, and the copper powder columns) through the wash bottle with H2SO4. The electric furnace with special construction provided a simple furnace opening. Automatic measuring and regulation of temperature during the reduction process was enabled. The reaction tube and silica carrier for NiO and NiCl2 samples ensured high temperatures operation.

A flow of nitrogen was injected into the reaction area in order to remove the residual oxygen and attain a neutral atmosphere. After the nitrogen inlet was stopped, hydrogen was introduced through the wash bottle (7) and the rotameter (8) into the silica tube containing the silica vessel (10) with the NiO or NiCl2 sample. Then, the silica tube was put into the electric furnace, previously heated to the selected temperature. The reaction time was measured from that moment. During the reduction, gaseous products were led away through the water-cooled condenser (12) and the wash bottle (14) to the gas spout.

After completion of the desired reduction time, the silica carrier was taken away from the tube and cooled down in a desiccator to room temperature. The reduction degrees were determined on the basis of the sample mass losing during the reduction.

The Avrami (I) equation was applied to data obtained for time-temperature dependence of reduction degrees. Then, activation energy values were determined from slopes of Arrhenius plots and by using the Arrhenius equation.

Reduction of NiO and NiCl2 was carried out at temperatures from 623 to 923 K in order to investigate the reaction ability and rate of the reduction processes at low temperatures. Obtained results were the basis for further research, where Ni/Al2O3 catalysts were synthesized by using oxide and chloride precursor for Ni. The aim was to select precursor that offers the following two possibilities during the catalyst synthesis: elimination of the calcination step and reduction of Ni/Al2O3 catalysts at low temperature (633 K).

For the synthesis of Ni/Al2O3 catalysts, the following materials were used: (MERCK, pro analysis) and alumina-based foam that was produced as described in our earlier work [26]. The catalysts with 20 wt% of Ni were obtained according to previously presented method [27]. First, aqueous solution of NiCl2 with metallic ion concentration of 0.15 mol/L was prepared. Then, the foams were placed in a tubular quartz reactor and impregnated with ultrasonically aerosolized chloride solution by using two different procedures. In the first procedure, impregnation at 773 K and calcination at 773 K were carried out for 1 h, which resulted in formation of oxide precursor for Ni on the foam surface. In the second procedure, chloride precursor was formed by impregnation at 473 K and drying at that temperature for 1 h. After the calcination and drying, the samples were reduced by hydrogen at 633 K for 1.5 h (static conditions, H2 flow rate: 20 L/h). Apparatus for the impregnation experiments was presented in the previous research of authors [27]. Subsequently, the following catalysts were prepared: Ni/Al2O3 and Ni-Pd/Al2O3 with 5 wt% Ni loading and Ni-Pd/Al2O3 with 20 wt% Ni loading. With respect to Ni, modified catalysts contained 0.1 wt% of added Pd. Those samples were synthesized according to the second procedure in order to obtain chlorides on the foam surface. The same materials were used for the synthesis of all the samples, except that the additional material was PdCl2 (MERCK, pro analysis). Lower amount of chloride solutions (0.15 mol/L metallic ion concentration) was used to decrease metal content. After the drying, the samples were reduced in hydrogen flow at 533 and 633 K for 1.5 h (static conditions, H2 flow rate: 20 L/h).

SEM analysis of the samples, before and after the reduction, was conducted in order to compare their properties.

3. Results and Discussion

The influence of reduction temperature and time on the NiO and NiCl2 reduction degrees is illustrated in Figure 2.

As expected, reduction degrees of NiO and NiCl2 increased with an increase in time and temperature. The peak reduction degree values reached 72 and 81 wt% for NiO and NiCl2, respectively, at 923 K for 5 minutes. Regarding the NiO reduction (Figure 2(a)), obtained results clearly show a pronounced induction period of 2 minutes, which decreases with an increase in temperature. In comparison with NiO, it is clearly visible that the induction period of NiCl2 was less pronounced and lasted about 1 minute (Figure 2(b)). This is the most probable explanation for slightly higher reduction degree of NiCl2, compared to NiO, at the same time and temperature. Induction period is an initiation of the reduction process, which occurs before actual reduction reaction is started. During induction period, crystallization centers of new phase (nuclei) start to develop and then begin to grow. After that, the reaction accelerates [7, 8]. In general, induction period lasts longer at lower reduction temperatures.

Kinetic analysis of the NiO and NiCl2 reduction by hydrogen was made on the basis of the experimental results shown in Figure 2, by using appropriate kinetic equations [28] (Table 1).

By testing these equations, it was concluded that the experimental results of the temperature influence on the NiO and NiCl2 reduction degrees were best fitted by Avrami I equation: where is a rate constant with unit of reciprocal time and is reduction degree (wt%).

Linearization of the experimental results by Avrami I equation is presented in Figure 3.

Equation (3) was developed assuming the cylindrical shape of particles, the accidental nucleation of new phase, and the spreading of the reaction front (surface, area) in the beginning stage of reaction. The excellent linear correlations between and time for various temperatures indicate that, during the induction period, the reduction of NiO and NiCl2 by hydrogen is controlled by chemical reaction.

The relationship between the basic kinetic parameters of a process (rate constant— and activation energy—) is usually described by Arrhenius equation:

The rate constant for various temperatures was determined from the slope of the curves presented in Figure 3 and obtained values were used for the Arrhenius plots. The activation energy values for the NiO and NiCl2 reduction during the induction period were calculated from the slopes of the plots, using Arrhenius equation, as shown in Figure 4.

The activation energy values for the temperature interval from 623 to 923 K were 29 and 28 KJ/mol for NiO and NiCl2 reduction, respectively (Figures 4(a) and 4(b)). Although the obtained activation energies are not typical for the kinetic region, time intervals of the NiO and NiCl2 reduction were short.

Considering that the experimental results were best fitted by using Avrami I equation (3), it was concluded that the reduction processes were controlled by chemical reaction. Piotrowski et al. [14] calculated activation energy values of 28 and 123 kJ/mol for the reduction of Fe2O3 by H2 and CO, respectively. The research group has reported that the processes were initially chemically controlled. For the reduction with both reducing agents, the same mechanism was suggested. Go et al. [15] investigated the water splitting reaction with Fe2O3, ZnFe2O4, and MnFe2O4. The reaction was controlled by the product-layer diffusion and activation energy values varied from 57 to 110 kJ/mol.

In the present research, slightly lower activation energy of the NiCl2 reduction indicates that the NiO reduction rate is slower compared to NiCl2. In addition, obtained results show that NiO phase is more stable and less reducible than NiCl2.

After the comparative kinetic analysis, it was concluded that using of chloride precursor is promising for the synthesis of Ni powder and Ni-based catalysts, enabling lower energy consumption.

In the further course of research, authors synthesized Ni/Al2O3 catalysts with 20 wt% of Ni. SEM micrographs of oxide precursor (after calcination at 773 K for 1 h) and Ni coating (after reduction at 633 K for 1.5 h) on the foam surface are presented in Figure 5.

Figure 6 shows SEM micrographs of chloride precursor (obtained by drying at 473 K for 1 h) and Ni coating (after reduction at 633 K for 1.5 h) on the foam surface.

The calcined sample, shown on Figure 5(a), had sponge-like agglomerates and uneven particle distribution. This morphology was most probably obtained because of the combined complex mechanisms of rapid water evaporation from the aerosol droplets, dehydration of , and oxidation of NiCl2 at 773 K. Due to cell morphology of the foam, its surface was not completely covered with oxide precursor. After the reduction at 633 K, a Ni coating completely covered the foam due to complex mechanism of mass transport. However, large residual sponge-like agglomerates were noted (Figure 5(b)). Calcination treatment increases stability of oxides, which is probable explanation for the residual agglomerates after the reduction. In contrast to that, SEM of the dried sample shows complete coverage of the foam surface with a chloride precursor layer (Figure 6(a)). For the synthesis of this sample, lower temperature was applied during the impregnation and drying. Evaporation of water from the aerosol droplets was slower, which resulted in this morphology of the chloride layer. The chloride layer had cracks and inhomogeneous thickness. The cracks appeared due to the evaporation of water from the aerosol and dehydration of at 473 K. The main reason for inhomogeneous distribution of chloride particles could be the foam cell morphology. Although the crust was cracked and had uneven thickness, relatively smooth Ni coating was formed after the reduction (Figure 6(b)). That coating covered the whole foam surface. In addition, cracks and agglomerates were not noted. The complete coverage of the foam with Ni coating and disappearance of cracks can be explained by complex mechanism of mass transport during the reduction. Mass transport took place when Ni atoms gained sufficiently high diffusion mobility at elevated temperature. Contacts between Ni particles were formed in cracks of the chloride layer and the growth of contact surfaces occurred. The morphologies of Ni coatings on the foam surface, obtained after the reduction of oxide and chloride precursors, indicate that the chloride phase has higher reducibility than the oxide phase. That is consistent with the comparative kinetic analysis results and it was confirmed in our earlier work [27]. Reduction degree of the sample obtained from chloride precursor (35.2 wt%) was approximately twice higher than that of the sample obtained from oxide (15.4 wt%) at 633 K [27]. It was concluded that the Ni/Al2O3 catalysts can be synthesized by the direct reduction of chloride precursor on the foam surface, without the calcination step.

The further course of research included preparation of Ni/Al2O3 and Ni-Pd/Al2O3 catalysts with 5 wt% Ni loading by using only chloride precursors for metals. These samples were synthesized in order to investigate influence of Ni loading and the addition of Pd on properties of the catalysts. For comparison purposes, Ni-Pd/Al2O3 catalyst with 20 wt% Ni loading was also prepared.

SEM micrographs of chloride precursor (obtained by drying at 473 K for 1 h) and Ni particles (after reduction at 633 K for 1.5 h) on the foam surface are shown in Figure 7.

SEM analysis of the dried sample clearly showed presence of fine, submicron-sized, and relatively unevenly dispersed chloride particles on the foam surface (Figure 7(a)). Lower amount of the chloride solution and the foam cell morphology led to formation of this microstructure. After the reduction at 633 K, submicron-sized, isolated, and island-like Ni particles were obtained (Figure 7(b)). Reduction degree of Ni/Al2O3 sample with 5 wt% Ni loading significantly increased and it reached a value of 67.7 wt% at 633 K. Considering Ni/Al2O3 sample with 20 wt% Ni loading, 35.2 wt% of chloride precursor was reduced at the same temperature [27]. The morphology obtained with lower chloride content resulted in higher reduction degree because fine isolated chloride particles had good exposure to hydrogen stream.

High reducibility of the catalysts at low temperatures (533 and 633 K) is achieved by modification with only 0.1 wt% of Pd. Active Ni and Pd particles can be obtained by the direct reduction of the corresponding chlorides, supported on a ceramic support. These chlorides are nearly completely reduced at 533 K: reduction degree of Ni-Pd/Al2O3 sample with 20 wt% Ni loading, and 0.1 wt% of Pd reaches a value of 98.2 wt% [27]. In the present investigation, properties of synthesized Ni-Pd/Al2O3 catalysts with different metal loading were compared. Metals were obtained by reduction of corresponding chlorides and lower reduction temperature (533 K) was utilized.

SEM micrographs of Ni-Pd/Al2O3 catalysts with 5 wt% and 20 wt% Ni loading (after reduction at 533 K for 1.5 h) are presented in Figure 8.

After the SEM analysis of reduced Ni-Pd/Al2O3 sample with 5 wt% of Ni, fine, submicron-sized, and isolated island-like particles were noticed on the foam surface (Figure 8(a)). Higher Ni loading (20 wt%) resulted in formation of relatively smooth Ni-Pd coating over the whole foam surface (Figure 8(b)). Ni-Pd/Al2O3 catalysts (reduced at 533 K) and Ni/Al2O3 catalysts (reduced at 633 K) had similar morphologies for the same metal loading. Chloride precursors for Ni and Pd were almost completely reduced at 533 K. Reduction degree of Ni-Pd/Al2O3 with 5 wt% Ni loading reached a value of 99.4 wt%. This value was only slightly higher than the reduction degree of Ni-Pd/Al2O3 with 20 wt% Ni loading (98.2 wt%) [27]. Considering that similar reduction degrees were obtained for different metal loading, Ni-Pd/Al2O3 catalysts with 20 wt% of Ni could be more suitable for catalytic processes because the metal coating completely covers the foam surface.

Presented results indicate that the calcination step can successfully be eliminated during the production of Ni-based catalysts. In addition to low-temperature reduction of the catalysts, it is possible to achieve further energy savings by using alumina-based ceramic foams that are sintered at low temperatures.

4. Conclusions

Analysis and study of the experimental results showed that the NiO and NiCl2 reduction by hydrogen occurred in the temperature range of 623 to 923 K. The NiCl2 and NiO reduction degrees reached peak values of 81 and 72 wt%, respectively, at 923 K for 5 minutes. Considering both systems, the duration of the induction period decreased with the increase of temperature. Under the same reduction conditions, slightly higher NiCl2 reduction degree was obtained due to less pronounced induction period (about 1 minute).

The experimental results of the temperature influence on the reduction degrees were best fitted by Avrami I equation, which is applied when the rate of a process depends on the rate of chemical reaction and when the reaction front spreads until the maximum reaction rate is achieved. The activation energy values for the induction period of the NiO and NiCl2 reduction were 29 and 28 KJ/mol, respectively, in the temperature interval from 623 to 923 K. Obtained results indicated that the rate of NiO reduction was slower compared to NiCl2.

It was concluded that the calcination step could successfully be eliminated by direct reduction of chloride precursor in the manufacturing of Ni/Al2O3 catalysts. Ni coating, supported on alumina-based foam, was synthesized from oxide and chloride precursor. Higher reducibility of the chloride phase was confirmed after the examination of prepared Ni/Al2O3 catalysts with 20 wt% Ni loading. Reduction degree of chloride on the foam surface reached approximately twice higher value (35.2 wt%) than the oxide (15.4 wt%) at 633 K. In contrast to Ni coating obtained from oxide precursor, when chloride was reduced at 633 K, relatively smooth coating without residual agglomerates and cracks was formed. When Ni/Al2O3 catalyst with lower Ni loading (5 wt%) was obtained by using chloride precursor, significantly higher reduction degree (67.7 wt%) was achieved at 633 K. Submicron-sized isolated Ni particles were formed on the foam surface. Additional energy savings in the production of monolithic Ni-based catalysts can be achieved by means of utilizing lower temperatures in each step of production. Alumina-based foams with enhanced properties can be sintered at low temperatures. Metal chlorides can be deposited on the foam surface by aerosol assisted routes and almost completely reduced at very low temperature, 533 K, if modified with addition of only 0.1 wt% of Pd. Desired metal loading is obtained by using different amount of chloride precursor solutions. Ni-Pd/Al2O3 catalysts with 5 wt% and 20 wt% of Ni were prepared by using chloride precursors for Ni and Pd. After reduction at 533 K, these catalysts had similar microstructures to Ni/Al2O3 catalysts with the same metal content, which were obtained by reduction of chloride precursor at 633 K. Reduction degree of Ni-Pd/Al2O3 catalysts reached a value of 99.4 wt% and 98.2 wt% for 5 and 20 wt% Ni loading, respectively.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This paper was done with the financial support of the Ministry of Education, Science and Technological Development of the Republic of Serbia and it is a result of Projects no. 34033 and no. 34023.