Study of local thermal nonequilibrium in porous media due to temperature sudden change and heat generation/Vietinio siluminio nepusiausvirumo priklausomybes nuo staigiu temperaturos pokyciu ir silumos generavimo poringoje aplinkoje tyrimas.

Kayhani, M.H. ; Abbasi, A.O. ; Sadi, M. 等

1. Introduction

Much effort has been devoted recently to determine conditions which guarantee that the assumption of local thermal equilibrium (LTE) is accurate when modeling of heat transfer in porous media. When it is accurate, then the thermal field is well-approximated by a single thermal energy equation. An excellent review of conductive effects in a stagnant porous medium may be found in Cheng and Hsu [1]. In their chapter these authors consider periodic media and their aim is to determine the effective thermal conductivity of the combined medium in the terms of the conductivities of the constituent phases. Therefore Cheng and Hsu provide important information for those wishing to use a single temperature field to model a two-phase saturated porous medium, or equivalently a composite solid consisting of two different constituents. In other circumstances, local thermal nonequilibrium (LTNE) prevails and it is necessary to employ two energy equations, one for each phase. The first papers which used two different temperature fields presented by Anzelius [2] and Schumann [3], and they were both published about eighty years ago. In their presented energy equations, we see that diffusion and advective (u dT/dx) terms have been neglected in the work of Anzelius. The numerical study by Combarnous [4] predated by a couple of decades further work on fully nonlinear convection using this model. Nakayama et al. [5] have proposed the nonthermal equilibrium two-energy equations model for conduction and convection, in which the two-energy equations for the individual phases at constant porosity are combined together and solved analytically. Neild and Bejan [6] stated the simplest equations which are generally regarded as modeling unsteady heat transfer in a saturated porous medium where LTE does not apply.

Great heat generation is one of the efficacious reasons to create LTNE condition between phases, (e.g., in the fluid phase this factor is appeared as chemical reaction). In the absent of fluid flow, Rees [7] determines both analytical and numerical formulae for interfacial heat transfer coefficient h in the porous media, when a uniform heat generation in fluid phase is produced and uphold LTNE condition. Nonthermal equilibrium heat transfer in the stagnant porous medium with variable porosity is analyzed by Nazari and Kowsari [8], where heat generation takes place within the solid phase. They use from energy equations of the solid and fluid phases with the assumption of steady-state and one-dimensional heat conduction.

Temperature sudden change is another effective factor for LTNE condition. When the temperature at the bounding surface changes significantly with respect to time in each phase, the local volumes of the solid and fluid phase can not react quickly and thereupon two equations are used to model the fluid and solid phases separately. In the two-field model, the energy equations are coupled by means of terms which account for the heat lost or gained from the other phase. Nouri-Borujerdi et al. [9] inspect the effect of LTNE on the evolution of the stagnant temperature field in a semiinfinite porous medium and then conduction takes place more rapidly in one phase than in the other, although local thermal equilibrium is always approached as time increases and in continuance Kayhani et al. [10] studied the effect of LTNE on a two-dimensional porous media under a step temperature change on the boundaries.

In summarize, special situations for LTNE to occur in the porous media are: 1) the great heat generation is happened in each phases (Chemical Reactions), 2) boundary temperature change suddenly along time, 3) hot fluid is injected in the cold porous media and 4) phases have different specific heat capacities and thermal conductivities.

In some of the LTNE applications such as fruit drying technology [11], heat transfer in biological tissues [12] and thermal analysis of the porous burner that is reviewed by Mujeebu et al. [13], temperature sudden change is governed and in some others, such as chemical catalyst [14] and nuclear reactors [15], heat generation affect on heat transfer process as well as previous condition.

In the present paper, we assume a two-temperature model for conduction in a stagnant porous medium that is saturated with the incompressible fluid which temperature change suddenly in boundary x = 0 and simultaneously uniform significant heat generation takes place within the solid phase. Also we consider how the heat generation in solid alters the behavior of temperature gradients in the different values of conductivity and thermal diffusivity ratio. The present paper is in continuance performed research projects by Nouri-Borujerdi et al. [9] and Kayhani et al. [10]. These projects followed from the similarity solution with complicated calculations and without presence of heat generation term. We simplify the solution method using direct numerical method and inspect the effect of heat generation as well as temperature sudden change in porous media.

2. Model development

As shown in Fig. 1, consider semiinfinite porous media that is saturated with a stagnant incompressible fluid. Using one-dimensional heat conduction in the porous media, the nonthermal equilibrium energy equations of the fluid and solid phases are as follows [6, 16]:

[epsilon])[([rho]c).sub.f][partial derivative][T.sub.f]/[partial derivative]t = [epsilon][nabla]([k.sub.f][nabla][T.sub.f])+h([T.sub.s] - [T.sub.f) (1)

(1 - [epsilon])[([rho]c).sub.s][partial derivative][T.sub.s]/[partial derivative]t = (1 - [epsilon]) ([nabla][([k.sub.s][nabla][T.sub.f]).sub.s] + q'") - h ([T.sub.s] - [T.sub.f) (2)

The subscripts f and s denote fluid and solid phases respectively. The quantities [epsilon], [rho] and c are the porosity, density and specific heat capacity and [rho] is the uniform heat generation per unit solid volume. The last term in energy equations represent the coupled heat transfer between the two phases because of the existing temperature difference. Many of scientists attempted to determine suitable values of h have generally relied upon averaging methods, and various assumptions then need to be made about closure; (see Rees [17]). Some of these formulations for determination of h, yield a zero value for h when Re = 0, which implies that there is no transfer of heat between the separate phases when the porous medium is stagnant. However, some of others yield nonzero values for h in the absence of flow, but the resulting expressions are independent of the conductivity of the solid phase. Based on presented models for h, we assume interfacial heat transfer coefficient as follows:

h = [epsilon][k.sub.f]/[L.sup.2] (3)

By introducing the following dimensionless parameters:

[??] = x/L (4a)

[??] = [[alpha].sub.f]/[L.sup.2] t (4b)

[[theta].sub.fluid] = [T.sub.f] - [T.sub.ref]/[q.sup.m.sub.f] [L.sup.2]/[k.sub.f] (4c)

[[phi].sub.solid] = [T.sub.s] - [T.sub.ref]/[q.sup.m.sub.f][L.sup.2][k.sub.f] (4d)

The governing Eqs. (1) and (2) can be nondimensionalized as:

Fluid [right arrow] [[theta].sub.t] = [[theta].sub.xx] + H ([phi] - 0) (5)

Solid [right arrow] [alpha][[phi].sub.xx] + H [gamma](theta] - [phi]) + 1 (6)

[FIGURE 1 OMITTED]

Where that the nondimension parameters are:

[alpha] = [[alpha].sub.f]/[[alpha].sub.s] diffusivity ratio (7)

H = h[L.sup.2]/[epsilon][k.sub.f] interfacial heat transfer coefficient (8)

[gamma] = [epsilon][k.sub.f]/(1 - [epsilon][k.sub.s] conductivity ratio (9)

The above parameters are constant and do not vary with temperature. This assumption helps us to verify the changes of phases together as well as simplification. Nouri-Borujerdi et al. [9] eliminate parameter H from energy equations by using natural coordinates provided by Carslaw and Jaeger [18]. But in this paper, by using the natural length scale in relation (3), parameter H will be equal to unit.

Assuming a high temperature sudden change at the boundary x = 0, LTNE condition between fluid and solid phases is possible and heat generation term amplify it. According to the Fig. 1, initial and boundary conditions can be simply defined as the following form:

t = 0 [right arrow] [theta]( x = 0 ) = [phi]( x = 0 ) = 1 (10a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (10b)

Boundary conditions (10b) showed a sudden change of temperature in the boundary of semiinfinite domain which induced LTNE between the phases.

3. Solution method

There are many numerical methods for solving ordinary differential equations which each of these methods have certain accuracy. In order to solve Eqs. (5) and (6), two different numerical methods have been used. For the calculation of second order derivatives, we use from compact finite difference and about time from forth order Runge-Kutta methods. Compact finite difference is the useful method to discrete domain with high accuracy. Basic of this method is very simple and similar to finite difference method but with the less error. For example in derivation, we use from backward, forward and central operators but in the compact finite difference method we mix them and use from an operator for derivation. This method was completed by Hirsh [19] and Lele [20] generalized it.

At early and late times we check the results using perturbation method. We determine the power series solution of Eqs. (5) and (6). At the suitable order it is possible to proceed easily analytically. We have to use from the other numerical solution (Shooting Method) at the end of analytical procedure.

4. Solution at early and late times

In this section we are going to examine results of numerical method in the special case at early and late times. In this case we assume no heat generation occurred in the solid phase. Fig. 2 shows the result of numerical method without heat generation when just temperature sudden change is the reason of heat transfer in porous media.

[FIGURE 2 OMITTED]

In this figure, at early times when x [right arrow] 0, we use from power series solution as follow:

[theta](x,t) = [[theta].sub.0] (x) + t[[theta].sub.l](x)+ [t.sup.2][[theta].sub.2](x)+ ... (11)

[phi](x,t) = [[phi].sub.0] (x)+ t [[phi].sub.1] (x) + t [[phi].sub.2] (x) + ... (12)

At [theta](t) by derivation than time and place and

then situation in the Eqs. (5) and (6), we obtain simplified equations as below:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (13)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (14)

So boundary equations (10b) change according to the power series solution as below:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (15)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (16)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (17)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (18)

We solve Eqs. (13) and (14) with boundary Eqs. (15) to (18) using from shooting method. Result of numerical solution of Eqs. (5) and (6) at the boundary x = 0 and early time t= 0.001 is equal result of analytical solution of Eqs. (13) and (14). This result for [DELTA]x = 0.15 is 1.1281.

We also repeat above progress for late time. In Fig. 2, at late times when x [right arrow] [infinity], we use from other power series solution:

[theta](x,t) = [[theta].sub.0](x) + 1/t [[theta].sub.1](x)+ 1/[t.sup.2] [[theta].sub.2](x) + ... (19)

[phi](x,t) = [[phi].sub.0] (x) + 1/t [[phi].sub.1] (x) + 1/[t.sup.2] [[phi].sub.2](x)+ ... (20)

At [theta](t) by derivation than time and place and then situation in the Eqs. (5) and (6), we obtain simplified equations as below:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (21)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (22)

We solve Eqs. (21) and (22) with previous boundary equations using from shooting method. The result of numerical solution of Eqs. (5) and (6) at the boundary x = 0 and late time t = 1000 is equal result of analytical solution of Eqs. (21) and (22). This result for [DELTA]x = 0.15 is 0.0333. Thereupon we can compare values of temperature gradients using results of Compact numerical method and perturbation method at early and late times. Compact numerical method is reliable and expandable to different cases.

5. Result and discussion

In this section we are going to present the complete description about behaviour of temperature gradients in the different conditions. Diffusivity ratio [alpha] and porosity-modified conductivity ratio [gamma] are two important parameters in this section. In the previous sections, we certified accuracy of results for temperature gradients using Compact numerical method.

Due to the same initial conditions in both phases, graphs in the Fig. 3 have the same start points. In the presence of heat generation, all of the graphs attain the steady state in the negative value of temperature gradients. In the through time domain, initial condition of phases is constant in the first node and therefore temperature of second node can get to higher temperature values than the first node. Finally, temperature gradients of fluid and solid phases are stabilized in the negative values. Temperature gradient difference is created from the difference between the solid and fluid thermal diffusivities. In the Fig. 3, a, with progress in time, this difference becomes greater. So, the rate of difference quicken by the existence of generation term in the solid phase. As we know, [alpha] is the ratio of fluid diffusivity to solid diffusivity and when [alpha] < 1:

[[alpha].sub.f]/[[alpha].sub.s] < 1 [right arrow] [k.sub.f]/[([rho]c).sub.f]<[k.sub.s]/[([rho]c).sub.s] [right arrow] [([rho]c).sub.f] (23)

Unequal Eq. (23) states that specific heat of solid phase is less than fluid phase. This means that solid phase will be heated and cooled too early rather than fluid phase. So when [alpha] < 1, the solid curve is placed lower than the fluid curve. According to the presence of heat generation term in this case, the decrease rate of temperature gradient of solid phase rather than fluid phase becomes greater and there is no contact point between solid and fluid graphs.

According to the unequal Eq. (23), in the Fig. 3, b when [alpha] > 1, specific heat capacity of solid phase is more than fluid phase and solid phase will be heated and cooled too late and the solid curve is placed higher than the fluid curve, but warming late value of solid phase decrease. In the Fig. 3, b, the difference between temperature gradients is less than Fig. 3, a. For the influence of continuous heat generation, this difference dwindles and finally curves obtain a coincidence point. At very late times, graph slope decreases to zero and the temperature gradient in both phases moves towards the infinite with a constant value. Similar slops zero at the late times represent that heat generation effect is counteracted. When the graphs slope is zero, heat transfer happens between phases yet, but it should be noted that due to the lack of heat generation effect, the amount of heat transfer always remains constant.

[FIGURE 3 OMITTED]

As it was stated previously, solid specific heat for [alpha] > 1 and [alpha] < 1, is larger and smaller than the fluid phase specific heat capacity respectively. The important point in Fig. 4 is that as the amount of [alpha] increase, time to counteract heat generation effect increases. Fluid phase graphs are stabilized in the negative value of temperature gradient with the constant slop zero and different diffusivity ratios [alpha] don't affect on this value. Solid phase has the similar status too.

[FIGURE 4 OMITTED]

We find out that different diffusivity ratios [alpha] do not affect on final values of temperature gradient in the phases, but according to the Fig. 5, a and b, different conductivity ratios change final values of temperature gradient in the phases and final distance between two phases. According to the relation (9), (1 - [epsilon])[k.sub.s] becomes greater than [epsilon][k.sub.f] with decreasing [gamma] or porosity decreases. With decreasing void volumes in the porous media, solid phase contribution in the heat generating increases more than before. This result is reverse with increasing [gamma] in the heat transfer process.

Fig. 6 depicts temperature gradients for [alpha] = 4 and different values of [gamma]. Heat generation effect is counteracted in earlier times when [gamma] increases.

[FIGURE 5 OMITTED]

[FIGURE 6 OMITTED]

For [gamma] = 1 and different values of [alpha], differences of temperature gradient between solid and fluid phases are depicted in Fig. 7. As we see, for [alpha] > 1, as a increases, the difference between solid and fluid phases will increase and vice versa, for [alpha] < 1 it will decrease. Better to say, by going away from the boundary [alpha] = 1 , the difference of the solid and fluid phases becomes more. Moreover, with reduction of [alpha], whether this coefficient is greater or less than unity causes that the point of maximum differences will tend to early times. Based on previous explanation, validating such argument is quite easy. With reduction of a and tending towards zero, fluid specific heat capacity tends towards infinity, although solid specific heat capacity will tend towards zero.

[FIGURE 7 OMITTED]

The differences of temperature gradient between solid and fluid phases are depicted in Fig. 8 for [alpha] = 4 and different values of [gamma]. The graphs state that when [gamma] increases, the value of the difference decreases and temperature gradient of the phases are more coinciding. Contact point of the phases tend to earlier time when [gamma] increases.

[FIGURE 8 OMITTED]

6. Conclusion

In this note, we have considered the local thermal nonequilibrium due to the temperature sudden change and great heat generation in porous media. Energy equations presented by Neild and Bejan [6] and Kaviany [16] were used as governing equations. After nondimensionalising, new parameters such as diffusivity ratio, scaled interfacial heat transfer coefficient and porosity-modified conductivity ratio were defined. Governing equations are solved numerically using Compact method and results valid using perturbation method. The effect of defined parameters on the behaviour of solid and fluid phases in porous media is investigated and results are presented in the form of various graphs. The results showed that after the passage of time, temperature gradients in both phases reached to a negative fixed amount and remain constant. The diffusivity ratio affected the behaviour and positioning of temperature gradient in the both phases. Also, the effect of porosity-modified conductivity ratio was discussed to reach equilibrium conditions. In two final figures, the difference amount of temperature gradient for all states between both phases was presented.

Acknowledgment

A.O. Abbasi author would like to thank associate professor M.H. Kayhani, M. Nazari and expresses his gratitude to mechanical engineering department in Shahrood University of technology for spiritually and Iranian Fuel Conservation Company (IFCO) for financially supporting this work.

Received August 30, 2010

Accepted January 17, 2011

References

[1.] Cheng, P.; Hsu, C.T. 1998. In: Ingham, D.B.and Pop, I. (Eds), Transport Phenomena in Porous Media I, Pergamon Title: 57-76.

[2.] Anzelius, A. 1926. Uber Erwarmung vermittels durchstromender Medien, Zeit. Schrift fur angewandte Mathematik und Mechanik 6: 291-294.

[3.] Schumann, T.E.F. 1929. Heat transfer: A liquid flowing through a porous prism, J. Franklin Instit. 208: 405-416.

[4.] Combarnous, M. 1972. Description du transfert de chaleur par convection naturelle dans une couche poreuse horizontale a laide d'un coefficient de transfert solide-fluide, C.R. Acad. Sci. Paris A. 275: 1375-1378.

[5.] Nakayama, A.; Kawakawa, F.; Sugiyama, M. and Xu G. 2001. A two energy equation model for conduction and convection in porous media, Int. J. Heat Mass Transfer 44: 4375-4379.

[6.] Neild, D.A. and A. Bejan. 2006. Convection in Porous Media. 3rd Edn., Springer, USA..

[7.] Rees, D.A.S. 2009. Microscopic modeling of the two-temperature model for conduction in heterogeneous media: Three dimensional media. In: Proceedings of the 4th International Conference on Applications of Porous Media, (ICAPM2009), University of CALGARY Turkey, August 10-12, 2009: 42-48. http://icapm2009.yeditepe.edu.tr/.

[8.] Nazari, M. and Kowsary, F. Analytical solution of nonequilibrium heat conduction in porous medium incorporating a variable porosity model with heat generation. -J. Heat Transfer., 2009, vol.131, issue 1, 014503.

[9.] Nouri-Borujerdi, A.; Noghreabadi, A.R. and Rees, D.A.S. 2007. The effect of local thermal nonequilibrium on impulsive conduction in porous media, Int J. Heat and Mass Transfer 50: 3244-3249.

[10.] Kayhani, M.H.; Maghrebi, M.J.; Sadi, M. and Nazari, M. 2009. The effect of local thermal nonequilibrium on a 2D symmetric porous. In: Proceedings of the 4th International Conference on Applications of Porous Media, (ICAPM2009), University of CALGARY Turkey, August 10-12, 2009: 85-90, http://icapm2009.yeditepe.edu.tr/.

[11.] Singh, R.P. and Medina, A.G. 1998. Food Properties and Computer-Aided Engineering of Food Processing System. -Boston, London: Kluwer Publisher.

[12.] Fung, Y.C. and Tang, H.T. 1975. Longitudinal dispersion of tracer particles in the blood flowing in a pulmonary alveolar sheet. Trans ASME, J. Applied Mech. 42: 536-540.

[13.] Mujeebu, M.A.; Abdullah, M.Z.; Abu-Bakar, M.Z.; Mohamad, A.A. and Abdullah M.K. 2009. Applications of porous media combustion technology-A review, Applied Energy. 86: 1365-1375.

[14.] Whitaker, S. 1986. Heat Transfer in Catalytic Packed Bed Reactors. Handbook of Heat and Mass Transfer. Houston: Gulf Publishing Co.

[15.] Lloyd, G.; Kim, K.J.; Razani, A. and Feldman, K.T. 1998. Thermal conductivity measurements of metal hydride compacts developed for high power reactors, J. Thermo-physics and Heat Transfer 12: 132-137.

[16.] Kaviany, M. 1995. Principle of Heat Transfer in Porous Media. 2nd Edn. -Berlin: Springer.

[17.] Rees, D.A.S. 2010. Microscopic modeling of the two-temperature model for conduction in heterogeneous media, J. Porous Media 13: 125-143.

[18.] Carslaw, H.S. and Jaeger, J.C. 1986. Conduction of Heat in Solids. Oxford University Press.

[19.] Hirsh, R.S. 1975. Higher order accurate difference solutions of fluid mechanics problems by a compact differenceing technique, J. Comput. Phys. 19: 90-109.

[20.] Lele, S.K. 1992. Compact finite difference schemes with spectral-like resolution, J. Comput. Phys. 103: 16-42.

M. H. Kayhani*, A. O. Abbasi**, M. Sadi***

* Mechanical Engineering Department, Shahrood University of Technology, P.O. 3619995161 Shahrood, Iran, E-mail: m.kayhani@yahoo.com

** Mechanical Engineering Department, Shahrood University of Technology, P.O. 3619995161 Shahrood, Iran, E-mail: ata.o.abbasi@gmail.com

*** Mechanical Engineering Department, Shahrood University of Technology, P.O. 3619995161 Shahrood, Iran, E-mail: m.sadi@gmail.com