Asymptotic behaviour for convection with anomalous diffusion

We investigate the fully nonlinear model for convection in a Darcy porous material where the diffusion is of anomalous type as recently proposed by Barletta. The fully nonlinear model is analysed but we allow for variable gravity or penetrative convection effects which result in spatially dependent coefficients. This spatial dependence usually requires numerical solution even in the linearized case. In this work we demonstrate that regardless of the size of the Rayleigh number the perturbation solution will decay exponentially in time for the superdiffusion case. In addition we establish a similar result for convection in a bidisperse porous medium where both macro and micro porosity effects are present. Moreover, we demonstrate a similar result for thermosolutal convection.


Introduction
Anomalous or fractional heat and mass diffusion defines deviations from the traditional diffusion model, which describes the movement of molecules as a random walk with a constant diffusion coefficient.In anomalous diffusion, the mean squared displacement of molecules grows with time slower or faster than linearly.With a slower growth, the regime is subdiffusive, meaning molecules spread out more slowly than expected.When the growth is faster than linear, the process is superdiffusive, meaning molecules spread out more rapidly than in standard diffusion.A thorough review on this topic can be found in Henry et al. [1].In recent work, Barletta [2] has proposed a model for convective motion of a fluid in a saturated porous medium where the diffusion coefficient is time dependent but may be of subdiffusion or superdiffusion type.He showed, for the linearized theory, that regardless of the size of the Rayleigh number the solution to the perturbation system may grow in the short term, but eventually the perturbation will decay for superdiffusion, and will grow indefinitely for subdiffusion.The object of this paper is to analyse the completely nonlinear problem, and we show that exponential decay holds in this case for superdiffusion.We establish this result in a saturated porous material but also allow for variable gravity effects, cf.Straughan [3], or for penetrative convection driven by an internal heat source, cf.Straughan [4, p. 343].This extension is important because mathematically the coefficients in the perturbation equations often depend on the vertical coordinate z.For such spatial dependence it is usually not possible to obtain an analytical solution by a normal mode procedure, cf.Barletta [5,6], and one must resort to numerical techniques.We additionally establish nonlinear decay of the solution when the porous medium is bidisperse, i.e. the medium has macro pores, but the solid skeleton possesses fissures or cracks, necessitating the inclusion of micropores, and hence a bidisperse or double porosity structure.Furthermore, we prove a similar asymptotic decay result when thermal and salt effects are present, i.e. in the case of thermosolutal or double diffusive convection The extension to spatially dependent coefficients is essential for real life applications.Also, consideration of double diffusive flow in a bidisperse porous medium is a topic of immense importance in current everyday life.For example, this phenomenon is proving important in renewable energy research, especially in solar pond technology, see Dineshkumar and Raja [7], Wang et al. [8].A particularly important application of double diffusion in a bidisperse porous medium is to magma flow in a volcano, see e.g.Vieira et al. [9], Allocca et al. [10], Toy et al. [11], Singh et al. [12], Bagdassarov and Fradkov [13], De Campos et al. [14].Given the recent seismic and volcanic activity in the Campi Flegrei region near Pozzuoli, an understanding of this scenario, and the potential effects of anomalous diffusion, is of vital importance.
We now present the equations for convection models and deliver a fully nonlinear analysis of asymptotic solution behaviour.

Thermal convection in a Darcy porous material
The subject of thermal convection in a Darcy porous material is investigated by Barletta [2] who allows for anomalous thermal diffusion by introducing a statistically motivated time dependent diffusion coefficient leading to a diffusion of form ϕD r rt r−1 ∆C, where C may be a concentration or a temperature field.The coefficients ϕ, D r , r are positive constants being porosity, constant diffusion coefficient, and an exponent.The variable t denotes time and ∆ is the three -dimensional Laplace The basic solution of Barletta [2] is the classical one of Chandrasekhar [15], where C (or T in the case of temperature) is linear in the vertical coordinate z and the velocity is 0.
We commence with the non-dimensional perturbation equations of Barletta [2, eqs.(36) -( 38)], which employ a Boussinesq approximation, cf.Barletta [16], and have form (in our notation) where u i , π, θ are the velocity perturbation, pressure perturbation, and temperature perturbation, w = u 3 , k = (0, 0, 1), and Ra is the Rayleigh number.It is convenient to rescale these equations and employ the parameter R = √ Ra and we write (1) in the form ( Throughout we investigate the superdiffusion case where r > 1. Barletta [2] essentially linearizes (2) and employs a normal mode analysis, cf.Barletta [5,6], to derive very interesting novel behaviour.In the interests of encompassing greater physical behaviour we allow for the effect of variable gravity, cf.Straughan [3], and for penetrative convection driven by an internal heat source, cf.Straughan [4, p. 343], and the perturbation equations (2) are replaced by where g, h are bounded functions of z, k is a positive constant and r > 1.The domain for equations ( 3) is (x, y) ∈ R 2 , {z ∈ (0, 1)} and t > 0. The boundary conditions are together with u i , π, θ being periodic in x, y.The periodicity ensures cellular structure of the convection cells as explained in detail by Chandrasekhar [15, pp. 43-52].Suppose for constants c 1 and c 2 .Then let V be a periodic cell for the perturbation solution to (3).Let (•, •) and • be the inner product and norm on L 2 (V ).
Multiply (3) 1 by u i and integrate over V , and multiply (3) 3 by θ and integrate over V .After integration by parts and use of the boundary conditions one may show and We next add ( 6) and ( 7) and then use the arithmetic -geometric mean inequality on the result to obtain d dt The function θ satisfies Poincaré's inequality ∇θ 2 ≥ π 2 θ 2 and then from (8) we derive Now employ an integrating factor and integrate (9) to see that It follows from ( 10) that as t → ∞, θ(t) tends to zero very rapidly no matter how large θ(0) is.From ( 6) one shows and so decay of θ(t) also implies decay of u(t) .

Remark
The exponential decay is obtained regardless of the size of the Rayleigh number Ra = R 2 .The result (10) applies also to the Barletta model where g = h = 1.Observe that we here employ no linearization and show decay for a solution to the fully nonlinear equations.

Convection in a bidisperse porous material
For a single temperature, T , in the macro and micro pores, equations for thermal convection in a bidisperse Darcy porous material are given by Gentile and Straughan [17].If we adopt the anomalous diffusion term of Barletta [2], then the non-dimensional perturbation equations for thermal convection in a bidisperse porous medium may be shown to be, cf.Gentile and Straughan [17], Straughan [18], Here u f i , u p i , π f , π p and θ are nonlinear perturbations to the velocity in the macropores, velocity in the micropores, pressure in the macropores, pressure in the micropores, and temperature, respectively, with w f = u f 3 , w p = u p 3 .The parameters ξ and K r are an interaction coefficient and the relative permeability K r = K f /K p , where K f and K p are the permeabilities in the macro and micro pores.
Equations ( 11) hold on {(x, y) ∈ R 2 }, {z ∈ (0, 1)}, t > 0, and the boundary conditions are together with u f i , u p i , π f , π p , θ being periodic in x, y.To derive an asymptotic behaviour result we multiply (11) 1 by u f i , (11) 3 by u p i and (11) 5 by θ and integrate each over V using integration by parts and the boundary conditions.After addition of the equations for u f i and u p i this leads to and d dt We now add equations ( 13) and ( 14) and we use the arithmetic -geometric mean inequality on the (θ, w f + w p ) terms to arrive at Now employ Poincaré's inequality on the last term and use an integrating factor as before to obtain This establishes that since r > 1, θ(t) decays rapidly as t → ∞, for any Ra and θ(0) .From ( 13) one shows Hence, from ( 16) decay of θ(t) guarantees decay of both u f (t) and u p (t) .

Double diffusive porous convection
In the case of double diffusive convection in a Darcy porous material when the layer is heated below and salted above or below, the basic solution is as in e.g.Straughan [4, pages 238-241], and the perturbation equations are given as, Straughan [4,eqs. (14.20)- (14.22)].If we started at the outset with the anomalous diffusion of Barletta [2] for the temperature of form k 1 t r−1 ∆θ and for salt as k 2 t s−1 ∆φ, for constants k 1 , k 2 and r > 1, s > 1, where φ is now the salt perturbation, then the non-dimensional perturbation equations are 0 where Ra = R 2 is the Rayleigh number, C = R 2 s is the salt Rayleigh number and Le is a constant known as the Lewis number.
We allow for different exponents of anomalous diffusion and without loss of generality we here assume s > r > 1.We consider the case of the minus sign in (18) 4 which corresponds to salting from above.The analysis in the other case of the plus sign is easier and we omit details.
To derive an asymptotic result we note that for t ∈ (0, 1), t r−1 > t s−1 .From (18) 1 we may obtain From (18) 3 and (18) 4 one finds d dt We add (19) and (20) and use the arithmetic -geometric mean inequality and Poincaré's inequality to obtain d dt Now put A = max{4R 2 , 4R 2 s }, B = min{2k 1 , 2k 2 }, and from (21) we may obtain for t ∈ (0, 1], where F (t) = θ 2 + Le φ 2 .Integrate ( 22) with an intergrating factor to obtain We now apply a similar argument to the above on the time interval (1, ∞), observing then that t s−1 > t r−1 > 1.In this case we obtain instead of (22), This is integrated with an integrating factor to arrive at where G(0) is as defined.Thus, Inequality (24) shows that as t → ∞, θ(t) and φ(t) both decay rapidly regardless of the size of Ra, C, or the initial data.By using the arithmetic -geometric mean inequality in (19) one shows Then from ( 24) and (25), u likewise decays as t → ∞.

Conclusions
We have extended the interesting result of Barletta [2] for the asymptotic behaviour of the solution to convection in a Darcy porous material for a superdiffusion model to the fully nonlinear case and we have shown that the perturbation velocity and temperature will always decay to zero, at least in L 2 norm.We have shown that this result may be extended to other convection in porous media scenarios.In particular, we allow for effects such as variable gravity, or convection with internal heat source.We also established an asymptotic decay result in the important problems of bidisperse convection and double diffusive convection.