Abstract

A mathematical model of the dynamics of the self-ignition of a reaction-diffusion system is studied in this paper. An approximate analytical method (modified Adomian decomposition method) is used to solve nonlinear differential equations under steady-state condition. Analytical expressions for concentrations of the gas reactant and the temperature have been derived for Lewis number (Le) and parameters , , and . Furthermore, in this work, the numerical simulation of the problem is also reported using MATLAB program. An agreement between analytical and numerical results is noted.

1. Introduction

Nonlinear dynamical phenomena in combustion process are an active area of experimental and theoretical research. Mathematical models that describe this phenomenon can be considered as nonlinear dynamical systems. The dynamic characterization of such models was developed by Continillo et al. [1]. The detailed numerical simulation of autoignition of coal stockpiles leads to the observation of steady regimes. To better investigate this phenomenon, two simplified distributed-parameter models were discussed which incorporate heat conduction, mass diffusion, and one-step Arrhenius exothermic chemical reaction. Both model equations were solved with straight forward finite-difference schemes [2]. The problem of spontaneous ignition of coal stockpiles is challenging for safety implications and for its theoretical complexity: a spontaneous combustion reaction takes place in a bed of solid fuel, while flow, driven by natural convection generated by the onset of temperature gradients within the pile, occurs. Coal stockpiles self-ignite when reaction of coal with oxygen present in the atmosphere generates heat, that is, not efficiently removed toward the external ambient [3]. Continillo et al. [4, 5] have analyzed the self-combustion of coal piles in the absence of natural convection. The three main phenomena in the self-ignition of coal stockpile are convection, reaction, and diffusion. On the other hand, Continillo et al. [6] studied the dynamic behavior of a two-dimensional coal pile also by accounting for natural convection. As part of a comprehensive study of self-heating of coal stockpiles, a simple mathematical model has been developed. To our knowledge, no rigorous analytical expressions of gas reactant () and temperature () have been derived for all possible values of parameters under steady-state conditions. The purpose of this paper is to derive approximate analytical expressions for gas reactant concentration and temperature using the modified Adomian decomposition method.

2. Mathematical Formulation of the Boundary Value Problem

The nonlinear differential equations are those of a distributed-parameter dynamic model of heterogeneous reaction in a one-dimensional layer. The gaseous reactant diffuses through the reacting medium and a first order one-step exothermic chemical reaction takes place. The reaction rate depends on the temperature through the Arrhenius exponential. The Arrhenius rate equation is a mathematical expression which relates the rate constant of a chemical reaction to the exponential value of the temperature. The model nonlinear equations in dimensionless form are [1] where is the concentration of the gas reactant, is the temperature, Le is the Lewis number (the ratio between mass and heat diffusivities),    is the dimensionless heat of reaction, is the thermal Thiele modulus (the ratio of the time scale of the limiting transport mechanism to the time scale of intrinsic reaction kinetics [7]), and is the dimensionless activation energy (minimum amount of energy between reactant molecules for effective collisions between them). The boundary conditions are Under steady-state condition, the equations become with boundary conditions

3. Analytical Solution of Nonlinear Dynamics of a Self-Igniting Reaction-Diffusion System under Steady-State Condition Using Modified Adomian Decomposition Method

In the recent years, much attention is devoted to the application of the Adomian decomposition method to the solution of various scientific models [8]. An efficient modification of the standard Adomian decomposition method for solving initial value problem in the second order partial differential equation yields the MADM. The MADM without linearization, perturbation, transformation, or discretization gives an analytical solution in terms of a rapidly convergent infinite power series with easily computable terms. The results show that the rate of convergence of modified Adomian decomposition method is higher than standard Adomian decomposition method [9–13]. Using this method (see Appendix A), we obtain an approximate analytical expression of concentration of gas reactant () and temperature () (see Appendix B) as follows:

4. Numerical Simulation

Nonlinear diffusion equation (3) for the boundary condition (4) is also solved numerically. We have used the function pdex1 in MATLAB software to solve the initial-boundary value problems for the nonlinear differential equations numerically. This numerical solution is compared with our analytical results in Figures 1–4. Upon comparison, it gives a satisfactory agreement for all values of the dimensionless parameters  ,  , and . The MATLAB program is also given in Algorithm 1.

function pdex4
m = 0;
x = linspace(0,1);
t = linspace(0,100000);
sol = pdepe(m,@pdex4pde,@pdex4ic,@pdex4bc,x,t);
u1 = sol(:,:,1);
u2 = sol(:,:,2);
figure
plot(x,u1(end,:))
title(β€˜u1(x,t)')
xlabel(β€˜Distance x')
ylabel(β€˜u1(x,2)')
%––––––––––––––––––––––––––––––––––
figure
plot(x,u2(end,:))
title(β€˜u2(x,t)')
xlabel(β€˜Distance x')
ylabel(β€˜u2(x,2)')
%––––––––––––––––––––––––––––––––––
function c,f,s = pdex4pde(x,t,u,DuDx)
c = 1; 1 ;
f = 1; 1 .* DuDx;
l = 0.233;
p = 35000;
b = 4.287;
g = 13.6;
F = βˆ’p*u(1)*exp(βˆ’g/u(2))/l;
F1 = b*p*u(1)*exp(βˆ’g/u(2));
s = F; F1 ;
%––––––––––––––––––––––––––––––––––
function u0 = pdex4ic(x);
u0 = 1; 1 ;
%––––––––––––––––––––––––––––––––––
function pl,ql,pr,qr = pdex4bc(xl,ul,xr,ur,t)
pl = ul(1)βˆ’1; ul(2)βˆ’1 ;
ql = 0; 0 ;
pr = 0; ur(2)βˆ’1 ;
qr = 1; 0 ;

5. Discussion

Equations (5) and (6) represent the simplest form of approximate analytical expressions for the concentration of gas reactant and temperature for all values of parameters ,  ,  , and . The Thiele number (thermal) is the ratio of layer thickness (L) and thermal diffusivity (). Equation (5) represents the new approximate analytical expression of concentration of gas reactant. The numerical solution is compared with the analytical results in Figures 1–4. These figures represent the analytical and numerical concentration profiles of gas reactant and temperature for different values of parameters  , , and . Figure 1 represents the dimensionless concentration versus dimensionless spatial coordinate for . From the figure it is inferred that the value of decreases when the value of or layer thickness increases. Figure 2 illustrates the dimensionless concentration versus dimensionless spatial coordinate for and we infer that the dimensionless temperature increases with the increase in values of layer thickness. Figure 3 represents the dimensionless concentration versus dimensionless spatial coordinate for values of . From the figure it is inferred that the value of decreases when the value of dimensionless activation energy () increases. Figure 4 illustrates the dimensionless temperature versus dimensionless spatial coordinate for the values of ; the dimensionless temperature increases with increase in values of dimensionless activation energy (). From Figures 2 and 4 it is evident that the maximum value for dimensionless concentration is 1 and that temperature attains its maximum value when the spatial coordinate  . Figure 5 confirms the results given by Figures 1 to 4.

Our analytical results are compared with the numerical results for the dimensionless concentration in Table 1. The maximum relative error between our analytical results and simulation results for the concentration is 1.3%. Also in Table 2, our analytical results are compared with the numerical results for the dimensionless temperature . Satisfactory agreement is noted. The maximum relative error in this case is 0.4%.

6. Conclusion

The system of steady-state nonlinear differential equations in different dynamic modes for the concentration of gas reactant and temperature has been solved analytically. The model investigated the influence of parameters over the temperature and concentration of gas reactant in the dynamic mode. The approximate analytical expression for the steady-state concentration of gas reactant and temperature is obtained using the modified Adomian decomposition method. A satisfactory agreement with the numerical result is noted. The analytical results will be useful to characterize the model predictions for the various values of parameters Thiele number (), Lewis number (Le), layer thickness (L), and thermal diffusivity .

Appendices

A. Basic Concept of Modified Adomian Decomposition Method

Consider the nonlinear differential equation in the form with initial condition where is a nonlinear real function, is the given function, and and are constants. We propose the new differential operator as follows: So, problem (A.1) can be written as The inverse operator is therefore considered a twofold integral operator as follows: By operating on (A.4), we have The modified Adomian decomposition method introduces solution and the nonlinear function by infinite series where the components of the solution will be determined recurrently and the Adomian polynomials of   are evaluated using the formula By substituting (A.7) and (A.8) into (A.6), By using the modified Adomian decomposition method, the components can be determined as which gives From (A.9) and (A.12), we can determine the components , and hence the series solution of in (A.7) can be immediately obtained.

B. Analytical Solution of Nonlinear Reaction-Diffusion Equation (3) Using Modified Adomian Decomposition Method

In this appendix we derive the general solution of nonlinear equation (3) by using modified Adomian decomposition method. We write (3) in the operator form as Applying the inverse operator on both sides of (B.1) and (B.2) yields where ,  ,  , and are constants of integration.

Let where where In view of (B.5), (B.7), and (B.8), (B.3) gives We identify the zeroth component as Using initial condition (4), we get = 0 and = 1; and the remaining components are as the recurrence relation where are the Adomian polynomials of . We can find as follows: In view of (B.6), (B.9), and (B.10), (B.4) gives We identify the zeroth component as Using initial condition (4), we get = 0 and = 1; and the remaining components are as the recurrence relation where are the Adomian polynomials of . We can find as follows: Hence solving (B.14) and (B.19) we get the following solutions: Adding (B.13), (B.21), and (B.23) we get (5) in the text. Adding (B.18), (B.22), and (B.24) we get (6) in the text.

C.

See Algorithm 1.

Notations

Symbols
:Concentration of the free stream
:Specific heat
:Mass diffusivity
:Activation energy
:Preexponential factor
:layer thickness
Le:Lewis number,
:Gas reactant
:Dimensionless temperature
:Temperature
:Free stream temperature
:Dimensionless time
:Time
:Dimensionless concentration .
Greek letters
:Thermal diffusivity
:Dimensionless heat reaction
:Dimensionless activation energy,
:Enthalpy of reaction
:Density
:Thiele number, .

Conflict of Interests

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

Acknowledgments

This work is supported by the University Grant Commission (UGC) Minor Project no. F. MRP-4122/12 (MRP/UGC-SERO), Hyderabad, Government of India. The authors are thankful to Shri. S. Natanagopal, secretary at Madura College Board, and Dr. R. Murali, principal at Madura College (autonomous), Madurai, Tamil Nadu, India, for their constant encouragement.