Modelling Bidispersive Local Thermal Non-Equilibrium Flow

: In this work, we present a system of equations which describes non-isothermal ﬂow in a bidispersive porous medium under conditions of local thermal non-equilibrium. The porous medium consists of macro pores, and in the solid skeleton are cracks or ﬁssures which give rise to micro pores. The temperatures in the solid skeleton and in the ﬂuids in the macro and micro pores are all allowed to be independent. After presenting the general model, we derive a result of universal stability, which guarantees exponential decay of the solution for all initial data. We further present a concrete example by specializing the model to the problem of thermal convection in a layer heated from below.

In addition, another area which has separately attracted much attention is flow in a bidispersive (or double porosity) material.Here, the usual macro pores are present, but in addition, the solid skeleton contains cracks or fissures which give rise to micro pores.This is studied in detail in, for example, Falsaperla et al. [20], Nield [21], Nield and Bejan [1], Nield and Kuznetsov [9,10], Straughan [19,22].
The goal of this paper is to combine both approaches and develop a theory for a porous material with a double porosity structure, but one which allows for different temperatures in the solid skeleton, the fluid in the macro pores, and the fluid in the micro pores.
It is very important to realize that there are situations where the local thermal non-equilibrium theory predicts very different temperature evolution profiles from what is observed with a single temperature; see, for example, the very interesting analyses of Rees and Bassom [15] and Rees et al. [16].Furthermore, in a very useful work, Rees [13,14] has shown how one may derive estimates for the interaction coefficients for the temperatures in a local thermal non-equilibrium theory, thus rendering this theory to be of practical value.In addition, David et al. [23], Homand-Etienne and Houpert [24], and Siratovich et al. [25] show where the inclusion of thermal effects can lead to stress-induced micro cracking in rocks such as granite, and so we believe that the inclusion of local thermal non-equilibrium effects in a bidispersive porous medium will be useful in real life.
We thus consider a porous body which has a porosity φ, i.e., the ratio of the volume of the macro pores to the total volume of the saturated porous material.The solid skeleton contains much smaller micro pores which may be cracks, or may even be due to a man-made structure (e.g., the picture in Nield and Kuznetsov [10], p. 3069).The micro pores give rise to a porosity , which is defined as being the ratio of the volume occupied by the micro pores to the volume of the porous body which remains after the macro pores are removed.This means that the fraction of volume occupied by the micro pores is (1 − φ) and the fraction occupied by the solid skeleton is It is worth pointing out that perhaps the major reason for the recent interest in double porosity materials is due to the many applications in real engineering and geophysics situations.For example, landslides and land movement due to thermal gradients (e.g., Hammond and Barr [26], Montrasio et al. [27]); provision of clean and safe drinking water (e.g., Ghasemizadeh et al. [28], Zuber and Motyka [29]); the controversial area of hydraulic fracturing for natural gas (e.g., Huang et al. [30], Kim and Moridis [31]); and many other applications may be found in Straughan [19].

Basic Model
As we have written in the introduction and reiterate now, the porous body is composed of a solid skeleton, macro pores, and micro pores.So, besides a macroporosity φ, there is a microporosity such that (1 − φ) and (1 − )(1 − φ) denote the fractions of the body occupied by the micro pores and the solid skeleton, respectively.We follow Nield and Kuznetsov [10] and denote by U f i and U p i the pore-averaged velocities in the macro and micro pores.In addition, Nield and Kuznetsov [10] denote by T f and T p the temperature in the macro and micro pores.The solid skeleton does not move, but we denote the temperature there by T s .
The momentum and continuity equations in the macro and micro pores are given by Nield and Kuznetsov [10] for a Brinkman porous body.We follow these writers, but employ a Darcy theory and so omit the Brinkman terms (cf.Straughan [19]).Thus, the momentum and continuity equations for the macro pores are: where µ is the dynamic viscosity, K f , p f are the permeability and the pressure in the macro pores, g i is the gravity vector, α is the coefficient of thermal expansion in the fluid, and ρ 0 is a constant.The term D = φ + (1 − φ) , ζ is an interaction coefficient, and standard indicial notation is used throughout with subscript , i denoting ∂/∂x i .
In writing (1), we have followed Nield and Kuznetsov [10] and used a Boussinesq approximation to write the buoyancy term as linear in the weighted temperatures so that the buoyancy term has the form where T = φT f + (1 − φ) T p /D, T 0 is a reference temperature, and in (1) the constant terms in (3) have been absorbed in p f .In a similar manner, the momentum and continuity equations in the micro pores have the form To write the balance of energy equations for the temperature fields, we let V f i and V p i be the actual velocities in the macro and micro pores.These are connected to the pore-averaged velocities by Thus, the equations for the temperature fields T s , T f , and T p are where and where (ρc) s , (ρc) f , (ρc) p are the products of the density and the specific heat at constant pressure in the solid, in the fluid in the macro pores, and in the fluid in the micro pores, respectively.The terms κ s , κ f , and κ p are thermal conductivities in the solid, and in the fluid in the macro and micro pores, respectively.We denote by s, f , and p the solid, the macro pores, and the micro pores.The terms h, s 1 , and s 2 are interaction coefficients, and we have here assumed that the interactions are linear in the temperature differences.The governing equations are (1), (2), and ( 4)-( 8).These equations hold in a bounded regular domain Ω ⊂ R 3 , for t > 0.
Remark 1.As we believe this is the first article to develop a local thermal non-equilibrium theory for a bidispersive porous material, we have followed the lead of Banu and Rees [3] and Nield and Kuznetsov [10] in ignoring fluid acceleration terms in Equations ( 1) and (4).For a Darcy porous material where we expect relatively slow flow, we believe that this is realistic.We could easily include acceleration (inertia) terms in (1) and (4), and very little would change in the rest of the article.

Remark 2.
In this work, we treat the thermal coupling terms h, s 1 and s 2 as constants.This is in line with the treatment for a local thermal non-equilibrium single porosity model by Banu and Rees [3], and also the interaction coefficients are treated as constant by Nield and Kuznetsov [10] in their development of a bidispersive porous model.For the present scenario, it is conceivable that the thermal interaction coefficients may depend on the solution, but this would lead to an extremely complicated model.In connection with this, we point out that Franchi et al. [32] study continuous dependence on the interaction parameters in a single-temperature bidispersive porous medium.We believe such a continuous dependence (or structural stability) result could be established for the model presented here, but the calculations will be technically involved and this is deferred to a later article.

Universal Stability
In this section we define the governing equations on Ω × {t > 0} with conditions on the boundary for given functions T 1 , T 2 , T 3 , where n is the unit outward normal to Γ.
We let Ū f i , Ūp i , p f , pp , Ts , T f , and Tp be a basic solution to P. We wish to analyze the stability of this basic solution.Hence, we define perturbations u f i , u p i , π f , π p , θ s , θ f , and θ p by We now non-dimensionalize the governing equations for the perturbations u f i , u p i , π f , π p , θ s , θ f , and θ p .We choose length and time scales L and T = φ(ρc) f /h, and we define non-dimensional numbers Λ h , Λ p , and Λ s by where U and T are velocity and temperature scales.The quantities L, U, and T are specified exactly when one considers a particular problem of thermal flow, as for example in Section 4. Hence, the non-dimensional perturbation equations may be shown to have the form From the definition of the perturbations and the conditions (9), the boundary conditions are The initial conditions are of form Our goal is to now develop a result of universal stability in the sense of Serrin [33].That is, we find conditions on the base flow to ensure u f i , u p i , θ s , θ f , and θ p decay exponentially in a suitable measure for all initial data (i.e., a global stability estimate).
Let • and (•, •) be the norm and inner product on L 2 (Ω).To achieve the above aim, we commence by multiplying Equation (11a) by u f i and we integrate over Ω.Further, multiply (11c) by u p i and integrate over Ω. Use integration by parts and the boundary conditions (12), and then add the resulting equations to obtain We use the arithmetic-geometric mean inequality on the right hand side of this identity, and one may thus show where κ 1 , κ 2 are given by The next step is to multiply (11e) by θ s , (11f) by θ f , and (11g) by θ p , integrate each over Ω, use integration by parts and ( 12), and then add the results to find Denote by I f and and then we may employ the arithmetic-geometric mean inequality on the right of ( 15) to show that We now employ (14) to estimate the right hand side of ( 17), and thus we arrive at where where we have also employed Poincaré's inequality and λ 1 is the first eigenvalue in the membrane problem for Ω.
From inequality (18), we see that a condition for universal stability is that the coefficients of the last two terms are positive.If we return to the definition of κ 1 and κ 2 , then we find that a sufficient condition for universal stability is that It is likely that I f and I p are measurable quantities in a given situation and the numbers R 1 andR 2 may be calculated.Thus, this universal stability estimate should be of use.However, we do not expect it to be optimum, and for a particular flow situation it is better to return to the precise equations at hand.In the next section we present the appropriate equations for thermal convection in a plane layer.

Thermal Convection
We now specialize Equation (11) to the case of thermal convection.Suppose that the saturated porous medium is contained in the layer with the temperatures where T L , T U are constants with T L > T U .The velocity boundary conditions are The steady solution in whose stability we are interested has form is the temperature gradient.
If u f i , u p i , π f , π p , θ s , θ f , θ p denote perturbations to the steady solution, then one verifies that these quantities satisfy the equations where g is the gravity constant, k = (0, 0, 1), w f = u f 3 , and w p = u p 3 .We non-dimensionalize (22) with the time, length, pressure, and velocity scales (ρc) f d 2 /κ f , d, dζU, κ f /(ρc) f d; the temperature scale is and the Rayleigh number R a is , and H by Then, (22) may be rewritten in the non-dimensional form where K SV = 1 κ s /φκ f and K PV = 2 κ p /φκ f , and In addition, we suppose u f i , u p i , π f , π p , θ s , θ f , θ p satisfy a plane tiling periodicity in the (x, y) directions (cf.Straughan [19], p. 189) for a typical hexagonal cell structure.
We do not attempt a complete analysis of the nonlinear stability and linear instability of system (24) here.While Straughan [22] has numerically verified that the linear instability results of Nield and Kuznetsov [10]-for the case where θ s is neglected and the Brinkman equations are considered-are close to the nonlinear stability ones, the same is not true when the Darcy equations are taken into consideration (see Straughan [19]).Since we might expect oscillatory instabilities with θ s ≡ 0, we believe a careful instability analysis of (24) will require substantial numerical computation.This, together with a careful nonlinear energy stability analysis, will be a lengthy article, and we defer this to the future.However, we may commence a nonlinear energy stability here to give a view as to how things will behave.
Let the period cell of the solution be V, and let • and (•, •) denote the norm and inner product on L 2 (V).To develop a fully nonlinear energy stability analysis, we multiply (24a) by u f i , (24b) by u p i , and integrate each over V.We likewise multiply (24c) by θ s , (24d) by θ f , (24e) by θ p , and integrate each over V.After addition of the resulting identities and use of the boundary conditions, one may show that where The global nonlinear stability threshold follows by computing max H I/D over a suitable space H (cf. Straughan [19]), and is then compared to the equivalent linear instability threshold.

Conclusions
We have produced a system of equations for a saturated double porosity body under conditions of local thermal non-equilibrium.This is effectively a generalization of the bidispersive theory of Nield and Kuznetsov [9,10], and the local thermal non-equilibrium theory for a single porosity material of Nield [7,8] and Banu and Rees [3].
We have shown that one may derive a universal stability estimate when a general base flow is considered.This is in the spirit of the original universal stability estimate for the Navier-Stokes equations by Serrin [33].When one requires accurate quantitative stability estimates, then one must address a specific flow problem as we show in Section 4 when we introduce thermal convection.