COMPARISON OF NUMERICAL SIMULATION RESULTS FOR TRANSPORT AND THERMODYNAMIC PROPERTIES OF THE SOLID FUELS COMBUSTION PRODUCTS WITH EXPERIMENTAL DATA

Abstract: Comparison of calculation results of transport properties of the solid fuels combustion products was made with known experimental data. Calculation was made by means of the modified program TETRAN developed in G.M. Krzhizhanovsky Power Engineering Institute. The calculation was spent with chemical reactions and phase transformations occurring during combustion. Also ionization of products of solid fuels combustion products at high temperatures was taken into account. In the capacity of fuels various Russian coals and some other solid fuels were considered. As a result of density, viscosity and heat conductivity calculation of a gas phase of solid fuels combustion products the data has been obtained in a range of temperatures 500-20000 . This comparison has shown good convergence of calculation results with experiment.


INTRODUCTION
Products of solid fuels combustion is difficult multicomponent system.Calculation of structure and properties of such system is an interesting and important problem.Similar researches can help to predict a course of reactions at phase and chemical transformations.The method and algorithm for calculation of solid fuels combustion products concentration and for calculation of transport properties is resulted lower and in previous papers.On the basis of this method, by means of modernized program TETRAN, the data on structure of various coal combustion products has been obtained and comparison of transport properties calculation results with known experimental data for pure gases, and for mixes is spent.

METHOD AND ALGORITHM
The system consisting from n components is considered.A current index of component number from n we designate ( 1,2,3 ) i i n , current indexes of phase number from r we designate , J E ( , 1,2,3 ) r J E . From a set of all components n we can allocate k linearly independent components.A current index of linearly independent components from k we designate ( 1,2,3 ) j j k .Any of i components can be expressed through linearly independent components j by means of set of stoichiometric factors ij ) .It is possible to form a rectangular matrix ij ) of these stoichiometric x Shmelkov Yuriy, Russia, Moscow, YurezSml@mail.ru;Samujlov Eugueny, Russia, Moscow, factors.If we designate a chemical symbol as , where ( ) A j E is basic components of ( ) j E , somehow distributed on r phases.For a finding the solid fuels combustion products structure the system of chemical thermodynamic equations for numbers of moles N i J all a component is uses.This system looking like: Types of the condensed phases are set preliminary, and their number is defined by Gibbs phase rule r k w 2 . The system of the nonlinear algebraic equations as a result of transformations is reduced to system of the iterative linear equations concerning logarithmic amendments for numbers of moles of linearly independent components which number is equal to number of various atoms in system, and full numbers of moles of the various condensed phases.Number of these equations much less than full number of the chemical thermodynamic equations.Algorithm of the linear iterative equations system decision is universal for change of quantity or structure of the condensed phases.Calculation is conducted, beginning from high temperatures where the concentrations of system atoms are set as initial approximation.Calculation is conducted simultaneously for two pairs independent parameters: T V , and T P, .In case of occurrence of problems with process convergence on one of pairs of variables the program automatically continues to consider and give out correct results only on variables T V , or T P, .After several steps on temperature «the bad point» manages and further calculation again on two pairs of variables is made.The program is provided by base of the thermodynamic given individual substances formed as a result of various reactions more than for 3000 substances.
The transport phenomenon in the gas phase in most cases satisfactorily describes Chapman-Enskog theory developed in the physical kinetics of gases for neutral gas [1], and further developed in the works [2,3].
For multi-component mixtures the expression for the coefficient of dynamic viscosity may be written as: . The kinematic viscosity coefficient and density shall be calculated in the following proportions: At heat conductivity calculation the big role have internal degrees of freedom and inelastic collisions of particles, including the collisions leading to chemical transformations.If inelastic collisions are frequent and probability of chemical reactions is great enough, so all components of gas are in local balance, with the big accuracy it is possible to use following expressions for forward heat conductivity of heavy components H O , internal heat conductivity int O and reaction heat conductivity R O : More details about method to calculate the composition and thermal properties is described in the article [4].

RESULTS
In this work comparison of calculation results of dynamic viscosity and heat conductivity with experimental data for nitrogen, water vapor and air is resulted.Data for comparison was received from the site of NIST (http://webbook.nist.gov/chemistry/)and from the book [5].In Fig. 1 comparison of the nitrogen dynamic viscosity with experiment is shown.The maximum calculation error makes 11 %.

Figure 1 :
Figure 1: Comparison of dynamic viscosity with experimental data for nitrogen In Fig. 2 comparison of the nitrogen heat conductivity with experimental data is shown.The maximum calculation error makes 12 %.

Figure 2 :
Figure 2: Comparison of heat conductivity with experimental data for nitrogen

Figure 3 :
Figure 3: Comparison of dynamic viscosity with experimental data for water vapor In Fig. 4 comparison of vapor heat conductivity with experiment is shown.The maximum calculation error makes 10 %.

Figure 4 :
Figure 4: Comparison of heat conductivity with experimental data for water vapor In Fig. 5 comparison of the air dynamic viscosity with experiment is shown.Data for comparison was received from the book [6].The calculation error makes 8,5 %.

Figure 5 :
Figure 5: Comparison of kinematic viscosity with experimental data for air

Figure 6 :
Figure 6: Comparison of heat conductivity with experimental data for air

Figure 7 :
Figure 7: Kinematic and dynamic viscosities of coal TROK I

Figure 8 :Figure 9 :
Figure 8: Kinematic and dynamic viscosities of coal CC2CCPOK I