Modeling of Rotary Kiln in Cement Industry

Cement production is a highly energy-intensive process, and the rotary kiln is the most important part of the process. Having a comprehensive model of the kiln in order to reduce manufacturing costs, better performance can be created. In this paper, the influence processes in a simulated cement rotary kiln and operating parameters on the output of the study were to develop and validate the systems using the same batch. The followings were examined: solid phase, gas and coating temperature change in a rotary kiln.


Introduction
In today's world, powering along the supply of raw materials, cement production is the most important factor.Cement production process is extremely time-consuming [1] and the rotary kiln is the most important part of the process, which is considered as the heart of it.Rotary kiln due to its complex nature (physical and chemical reactions simultaneously) and various parameters affecting the quality of research on the dynamics of complex, is very necessary [2].Ongoing maintenance and operation of the rotary kiln are vital.It is very expensive, because that would stop the kiln production.After having a comprehensive model of the process, we can create a better performance of the kiln, and the model can be used for operator training by experienced operators expert in the industry used [3], and eventually a reduced energy consumption in the cement industry.

Description of Cement Rotary Kiln
A simple system for cement rotary kiln is shown in Figure 1.Rotary kiln is basically a cylinder of length 50 -120 m (modern kilns) and 3.5 to 6.5 feet in diameter.
During kiln mainly on how long it takes until the temperature of the raw material clinker temperature of about 1370˚C is, it depends.The modern design of the plant, a cyclone preheater to increase the temperature of the solid input materials are used to shorten the length of the kiln.
Raw materials to the kiln consist of calcium carbonate (CaCO 3 ), silica (SiO 2 ), shale (Al 2 O 3 ) and iron oxide (Fe 2 O 3 ) feeds.These compounds are very fine powder and mixed with cement to form [4].

Energy Balance
For the balance of power in the rotary kiln, the average temperature of the walls was used in this model.The following equations for the gas, and solid walls, respectively [5] written.
According to the equation: where i m  (i mass flow rate in kg/s), C pi (i specific heat capacity in kJ/(kg•˚C)), ΔT i (i temperature change according to ˚C), Q (heat in kJ) is, and ṁ i replaced by the following equation: That ρ i (i density in kg•m 3 ), v i (i speed in m/s), A i (i surface area in m 2 ).Therefore the temperature in the gas phase transport equations are: where T w (wall temperature by ˚C), T g (gas temperature in ˚C), T s (solid temperature in ˚C), β 1 (heat transfer coefficient between the wall and the gas in W/˚C), β 2 (transfer coefficient between the solid and the gas temperature in W/˚C), Q comb (combustion heat in W), v g (gas velocity in m/s), ρ g (average gas density of 0.85 kg/m 3 ), ∂T g /∂z (gas temperature changeover the elements in ˚C /m), C pg (gas specific heat capacity of 1173.8 kJ/kg/˚C) and A g (of gas per m 2 ) and v g obtained from the following equations: ( ) That f m  (fuel mass flow rate) and a m  (air mass flow rate) is: That r 1 (radius of the kiln according tom) and p (surrounded by a solid angle 3π/2) and the solid phase are: That β 3 (heat transfer coefficient between the wall and solid in terms of W/˚C), Q c (heat of reaction in terms of W), v s (solid velocity in m/s), ρ s (solid density of 890 kg/m 3 ), ∂T s /∂z (solid temperature changes over the elements in ˚C/m), C ps (specific heat capacity of solid 1089.97 kJ/kg/˚C) is, A s (solid surface in m 2 ) and v s calcu- lated from the following equations: ( ) Wall: That β 4 (heat transfer coefficient between the wall and the environment in terms of W/˚C) and to calculate the heat of reaction is: where ) h o fraction of the radiation 0.0757, f 1 (conductivity gas to the wall), f 2 (conductivity solid to gas), f 3 (conductivity wall of the gas), ε w = 0.9, ε s = 0.8, ε g = 0.1, respectively coefficient gas, solid wall and h (radiation coefficient between solid and wall W/m 2 /˚C) is calculated from the following equation: That f 4 (conductivity of the wall) and To calculate the heat transfer coefficient between the wall and the environment (external air), Spang of Equation ( 16) is used.This is a perfect description of the heat transfer coefficients between the shell and the environment.Since ƒ 4 the kiln was fixed, β 4 not sensitive to the conditions in the longitudinal direction.To increase the accuracy of the model, the heat transfer coefficient of the outer shell of the total displacement csh a h → and radiation Rsh a h → was considered as follows: ( ) That Re (Reynolds number), Gr (Grashof number), k a (air conductivity), Pr (Prandtl number), D (shell diameter in meters) is: where C in Equation ( 18) is obtained as follows: where T a (ambient temperature on ˚C), T sh (shell temperatures in ˚C) and ε sh = 0.5 (emission shell) and σ = 5.6697 × 10 −8 W/m 2 •˚C 4 .As to the equation ( 16) deduced csh a h → is temperature dependent.To Pr, Re and Gr can be strongly influenced by environmental conditions.Therefore, to improve the model, first, the temperature of the kiln shell, the shell temperature scanner must be registered.
To calculate the film temperature (T f according ˚C) fluid inside the kiln elements are considered in the following equation can be used [6]: In each of the elements considered by the kiln shell surface temperature T sh in ˚C, T a is the ambient temperature in ˚C.So much csh a h → can be calculated at any point in the kiln.After calculating the temperature of the kiln is calculated film dimensionless groups.Heat transfer in the kiln body and Reynolds number Grash of What ratio are studied.If the Reynolds number is greater than the number Grashf ( where Nu forced (Nusselt forced displacement) and: where Nu free (free convection Nusselt number) and Ra (Riley number) [7]:

Flame Model
Plug flow of cement rotary kiln flame model is used [8].In this model, the flame is divided into n cut to size.The thickness of each slice (step-size) mathematical equations solver that will be discussed later.
The number of sections is calculated as follows: Simulation step-size Flame length overall, F L , the equation Beer is obtained: The AF * (stoichiometric air-fuel ratio), ρ e (equivalent to gas density), ρ cp (density of combustion products), ρ se (solid density published), d 0 (equivalent diameter of the burner) is.Gorog, cement kiln burner Coaxial type equipped [8]: ( ) The f m  (fuel mass flow rate in kg/s) and sa m  (secondary air mass flow rate in kg/s) can be recorded from the primary data cement kiln.The pa m  (primary air mass flow rate in kg/s) is calculated by knowing sa m  and 10% excess air to the burner, is calculated.
G F and also G Pa in Equation ( 26) to the flow rate of fuel and primary air.Now Q comb in Equation ( 3) is calculated by the following expression: where T sa (depending on the secondary air temperature ˚C), LHV (lower heating value fuel in J/kg), T pa (primary air temperature in ˚C), T ref (adjusted reference temperature ˚C) is.

Rotary Kiln Reaction Kinetic Rates
Rates of chemical reactions rotary kiln using Arrhenius equation can be expressed as follows [9]: The solid temperature T s by ˚C, and R = 8.314 J/mol•˚C gas constant, A i (i based on the frequency factor 1/h) and E i (i activation energy in J/k mol) and k i (Arrhenius constant reaction i by 1/h) is.
In which:

CaO,C S,C S,C A,C AF, H O,CO i =
And by chemical reactions in cement rotary kiln, we have: In these expressions M i (i molecular mass in kg/m 3 ) and R i (kinetic rates i), C i (level i) is.

Mass Balance
Mass imbalances, Equations ( 30) and (31) for 10 solid compound involved in the reactions occurring in cement kilns and other gas (CO 2 ) is written: R i (kinetic rates of the solid phase) and R j (gas phase kinetic rates) are presented in the previous section.Kinetic parameters are presented in the following Table 1.

Modeling Assumptions
These assumptions are as follows: • Internal and external diameter of the kiln is assumed.
• The specific and reaction heat were independent of temperature and they were constant along the axial direction.
• Conduction in gases and solids in the axial direction of the wall was ignored.
• Displacement and diffusion coefficients are independent of temperature and location.
• Height and speed of the solids in the kiln section is assumed.
• solid material conveyed by the exhaust gas flow is not included in the model • Arrhenius reaction rate specified by law.
• Average amount of coating conductivity 0.74 W/m 2 /˚C was considered.
• The conductivity of the refractory could be estimated by Equation (32) which was correlated from the experimental data, given by the refractory manufacturer for the magnetite-fried brick type is: • Conductivity of the metallic shell 45 W/m 2 /˚C is considered close to the alloy carbon steel.
• Since conductivity changes during the refractory lining is not possible, the starting amount was fixed.This assumption can be a source of grave error in the prediction model.
• The numbers of points scanned shell temperature kiln for a full rotation every 40 elements.At any point in the calculation of the axial position of an average value of all points in between.

Cement Rotary Kiln Studied
Length, diameter and outdoor cement kiln 70, respectively, 4.1 and 4.5 meters.Average thickness of 15 cm refractory magnetite, which can be considered to be uniform throughout the cooking area.Tilt cylinder to facilitate axial displacement of the solid bed, moving toward the discharge end in a state where the hot gases circulating in the opposite direction, was 4%.

Shell Temperature Scanner
In Figure 2, the temperature measured by the scanner shell is provided.As can be downloaded from the possibility of severe changes in temperature measurement scale, there is a scanner.The reason for this phenomenon is the thickness of the coating.Since the coating thickness of 10 to 15 inches of cooking area is ideal for protection from fire in all areas [10].To keep the coating thickness, shell temperature should be between 190˚C and 220˚C.Temperature below a solid barrier to the movement of the kiln and the refractory layer is higher than a disadvantage.

Simulation
Mass and energy equations are integrated with a set of differential equations and algebraic models are flame by software MATLAB R2014a (ver.8.3) is solved.Equations for the particular type of cooking area Hlknndh ODE15S are used.After solving the model, the wall temperature profile to the application submitted by Hlknndh ODE is solved.To solve the model equations, ODE during kiln divided into 70 elements.The size of each element of the program was considered a meter.The generated two-dimensional matrix, the temperature of the gas phase, solid coating in the oven during the show.Resistance between the wall and the surface layer is shown in Figure 3.

Calculate the Energy Loss
Heat flow in cylindrical coordinates (without heat) as follows [9]: For one-dimensional mode with the following assumptions, Equation (33) can be simplified: Assumptions 1) Steady-state heat transfer through the layers of the wall of the kiln (Steady-State) was considered.
2) The heat conductivity in the z direction was ignored.Thus, according to Figure 3, the following boundary conditions must be defined.
Coating layer: at , , & at , At each stage (ΔZ), in the inner wall of the kiln temperature (T w ) to find.Equations ( 1) to (28) are solved simultaneously in the program.Then, using Equations (38), (37) and (36), Q pass , T b , T c , respectively, can be calculated.The simulation code, both inside the kiln wall temperature (T w ) and the temperature of the kiln shell (T sh ) as a two-dimensional matrix in terms of the length of the kiln is considered.We can conclude that the thermal  resistance Equation ( 35) to (38) can be used to adapt the matrix elements, developed.Therefore, the heat loss from the body of the kiln is calculated (see Table 2).

Results and Validate
Now using equations results in mind, we can change the temperature of the solid phases, gas and coating in a rotary kiln were investigated (Figure 4).With the inlet and outlet temperature of the gas and solid modeling results with experimental results, we can infer that the model has good functionality.
Also according to the model equations, we can obtain the amount of waste in the rotary kiln (Figure 5).

Discussion and Conclusions
Since the temperature changes of the gas phase, solid coating of cement kilns was essential; therefore, an integrated model was developed in the cement kiln.First, a mathematical model of steady state (Steady-State) to estimate the inner wall surface of the rotary kiln temperature profiles was formulated.Then, by calculating the  temperature profile in the kiln and temperature profiles measured by the outside, the heat dissipation of body heat transfer resistance in the adjacent layers of the cylinder was estimated by the model.A comparison of model and data was collected from an industrial kiln, confirming that the feature is good for heat dissipation body temperature changes and the different phases.In addition, it is concluded that the heat loss from the body in the rotary kiln sintering temperature is greater than the thickness of the coating layer of the coating material, which does not create a problem.

Figure 1 .
Figure 1.A simple view of a cement rotary kiln.
μ (kinematic viscosity in terms of kg/m•s), ρ (density in kg/m 3 ), u (velocity in m/s), d (diameter in m), g (acceleration due to gravity in m/s 2 times with 9.8), T surf (depending on the surface temperature ˚C), T a (depending on the ambient temperature ˚C), L (characteristic length in m), ϑ (dynamic viscosity in m 2 /s), k a (air conductivity in W/m•˚C), β (coefficient of gas expansion in 1/˚C) [6] and on the free movement:

Figure 2 .
Figure 2. The temperature in the kiln shell scanner studied.
above equation with boundary conditions (34), the following equation is obtained:

Figure 4 .
Figure 4.The temperature changes of the gas phase, the phase of the study materials and coating over 70 m rotary kiln.

Figure 5 .
Figure 5.The total amount of waste in the body of the 70-meter rotary kiln studied.

Table 1 .
Parameters of reaction kinetics of cement rotary kiln.

Table 2 .
Rotary kiln operating conditions studied.