Research Article | Open Access
Instability of Vertical Constant Through Flows in Binary Mixtures in Porous Media with Large Pores
A binary mixture saturating a horizontal porous layer, with large pores and uniformly heated from below, is considered. The instability of a vertical fluid motion (throughflow) when the layer is salted by one salt (either from above or from below) is analyzed. Ultimately boundedness of solutions is proved, via the existence of positively invariant and attractive sets (i.e. absorbing sets). The critical Rayleigh numbers at which steady or oscillatory instability occurs are recovered. Sufficient conditions guaranteeing that a secondary steady motion or a secondary oscillatory motion can be observed after the loss of stability are found. When the layer is salted from above, a condition guaranteeing the occurrence of “cold” instability is determined. Finally, the influence of the velocity module on the increasing/decreasing of the instability thresholds is investigated.
Convection in fluid mixtures saturating porous media has attracted—in the past as nowadays—the attention of many scientists due to its practical applications like, for example, in geothermal energy exploitation, extraction of oil from underground reservoirs, ground-water pollution, underground flows movement, thermal engineering, crystal growth, polymer engineering, and ceramic processing. For this reason, several studies have been addressed to this topic [1–10].
The models describing the fluid motion in porous media are reaction-diffusion dynamical systems of P.D.Es, which, as it is well known, play an important role in the modeling and studying of many phenomena , for instance, [11, 12] and references . Several geophysical and technological applications involve nonisothermal flow of fluids through porous media called throughflow (i.e., there is flow across the porous medium and the basic flows nonquiescent) which affects the stability of the system significantly. Precisely, in some situations, a vertical motion (throughflow) is observable in a horizontal porous layer heated from below and the problem to determine until this motion is stable is of fundamental importance especially in applications involving cloud physics, hydrological/geophysical studies, seabed hydrodynamics, subterranean pollution, and many industrial and technological processes [1, 13–16]. This is because, after the loss of stability, a secondary motion arises and this motion can be steady or oscillatory.
The effect of vertical throughflow on convective instability, in either the fluid or the porous layer, has been extensively discussed by several investigators. In the fluid layer, the problem is of interest because of the possibility of controlling the convective instability by adjusting the throughflow . In the porous layer, the “in situ processing of energy resources” such as coal, oil shale, or geothermal energy often involves the throughflow in the porous medium. The importance of buoyancy-driven instability in such systems may become significant when precise processing is required . Besides, the throughflow effect in such situations may be of attention due to the opportunity of controlling the convective instability by regulating the throughflow in accumulation to the gravity. Throughflows play an important role also in the directional solidification of concentrated alloys, in which a mushy zone exists and regarded as a porous layer . As demonstrated in several studies, this problem also has a very strong role in mineralization and ore body formation in fluid-saturated porous rock masses. For example, in  a hydrothermal system consisting of a horizontal porous layer with upward throughflow has been considered. This horizontal layer of a fluid-saturated porous medium may undergo constant temperature and overpressure at its top, whereas it may undergo constant vertical temperature gradient and eventually constant injection of mass flux at its bottom. The overpressure can result from the presence of impermeable seals in geophysics, while it can be induced by surface structures or human activities in geoenvironmental engineering. The vertical temperature gradient may be generated by either geothermal sources in geophysics or buried heat-generating waste in geoenvironmental engineering. The eventual mass flux may be generated by dehydration reactions or by compaction of underlying soft layers or reservoirs. In this situation, it has been studied the destabilizing effect caused by the increase in upward throughflow. Rapid developments in modern technologies, during the recent past, have posed challenges in studying convective instability problems in different fluid dynamical systems. In fact the effect of vertical throughflow on the convection has been studied in many cases [13–25] (for example, the effect in a rectangular box is considered in ; the effect combined with a magnetic field in ; stability analysis when the density is quadratic in temperature in ; the effect with an inclined temperature gradient in ). The stabilizing effect of rotation and Brinkman law (holding for large pores) on a vertical constant throughflow has been analyzed in , while in [24, 25] the stability analysis of a vertical constant throughflow in binary mixtures saturating horizontal porous layers has been performed. In the previous papers, only the stability of a throughflow has been investigated, neglecting the question to determine the critical Rayleigh thermal numbers at which instability sets in and how instability arises.
In the present paper, we will focus on the instability analysis of a vertical constant throughflow in a horizontal porous layer, with large pores, uniformly heated from below and uniformly salted by one chemical either from above or below. In particular we determine the critical Rayleigh thermal numbers at which instability occurs and investigate for the kind of secondary motion arising.
Section 2 deals with the introduction of the mathematical problem. The onset of instability is investigated in Section 3. In particular, the critical Rayleigh thermal numbers at which steady or oscillatory instability can occur are determined and sufficient conditions guaranteeing that—after the loss of stability—a secondary steady or oscillatory motion arises have been determined. Further, in the case of a layer salted from above, the onset of “cold” instability (i.e., instability independent of the temperature gradient) has been investigated. The influence of the module of the vertical throughflow velocity, on the instability thresholds, is analyzed in Section 4. The paper ends with Section 5 in which the obtained results are highlighted.
2. Statement of the Problem
Let us consider a fluid moving mixture saturating a horizontal porous layer of depth , uniformly heated from below, with large pores, salted by a chemical specie (salt) either from below or from above. Denoting by an orthogonal frame of reference with fundamental unit vectors ( pointing vertically upwards), the set of governing equations for the fluid motion in the Darcy-Oberbeck-Boussinesq scheme, according to the Brinkman law, is given by where is the fluid mixture density and , , , , , , , , .
The subscripts and refer to the solid matrix and fluid, respectively, and . To (1) we append the boundary conditions where are positive constants such that . On considering the following dimensionless variables where according to the layer is salted from below or above, system (1), omitting all the primes, reduces to where and the boundary conditions (3) become A vertical constant throughflow solution of (5)-(7) is given by where is a constant. Setting system (5) becomes where , , and To (10) we append the initial-boundary conditions with . In the sequel, as usually is done in convection problems in layers, we suppose that (i)the perturbations are periodic in the and directions of periods and , respectively;(ii) is the periodicity cell;(iii) belong to and can be expanded in Fourier series, uniformly convergent in , together with all their first derivatives and second spatial derivatives.
Let us denote by(i) the norm in ;(ii) where .
3. Instability Analysis
In this section we study the instability of the throughflow solution (8). To this end, let us neglect the nonlinear terms in (10), and let us denote by the solution of the linearized version of the system (10)-(12). On taking the third component of the double curl of (10)1, one obtains that the linear system governing the evolution of is under the boundary conditions On looking for solutions of normal modes type, in view of periodicity in the and directions, one has with . Setting it easily turns out that where . Substituting (16) in (14), one obtains under the boundary conditionsIn order to determine an approximate solution of (19)-(20), let us employ an order-1 Galerkin weighted residuals method . We choose as trial functions (satisfying the boundary conditions (20)) where are constants. Substituting (21) in (19), we obtain three residuals. Making these residuals orthogonal to the trial functions over the range , one obtains a system of three linear algebraic equations in the three unknown . The vanishing of the determinant of the coefficient matrix, gives the instability threshold with . When the vertical constant throughflow (8) loses its stability, a secondary motion (steady or oscillatory) arises.
From (16), the critical Rayleigh number at which a secondary steady motion arises (say ) is obtained on substituting in (22) and is given by It turns out that (see  for details) with Hence the critical Rayleigh thermal number at which steady instability can occur is The critical Rayleigh number at which a secondary oscillatory motion arises (say ) is obtained on substituting , in (22). Setting the imaginary part equal to zero it follows that Substituting in the real part of the right-hand side of (22), it follows that Then—in view of (24)—the critical Rayleigh thermal number at which oscillatory instability can occur is We remark that (27) excludes a oscillatory motion not only when but also when We have tested numerically that the order-1 Galerkin approximation gives good estimates of the critical Rayleigh thermal numbers if is small compared to unity, as it usually happens for problems of this type. Hence, from now on, we will assume that is small compared to the unity.
Now, limiting the analysis to the case in which the layer is salted from above (), we investigate for sufficient conditions ensuring the loss of stability of the vertical throughflow, independently of the thermal Rayleigh number, i.e., we look for the onset of “cold convection” . Setting the following theorem holds.
Theorem 2. In the case , if either or holds true, then cold instability arises.
Proof. Let us first observe that or implies that the throughflow (8) is unstable. Hence the proof follows since either (32) or (33)—guaranteeing, respectively, that the second right-hand side of (26) and (29) is negative—ensures instability independent of .
According to the definition given in , the critical solute Rayleigh number for the onset of cold convection is .
3.1. Onset of Secondary Motions
In order to establish if a secondary steady or oscillatory motion arises when (8) loses its stability, we have to compare the critical Rayleigh thermal numbers for the onset of steady and oscillatory instability. Since the order-1 Galerkin method gives a good approximation of the critical instability thresholds and has to capture the physics of the problem, we assume that the comparison relation between the exact critical instability thresholds is the same of that one between and . The following results hold true.
Theorem 3. If either or then instability occurs via a steady state at .
Theorem 4. If then instability occurs via an oscillatory state at .
Proof. The proof follows since (37) guarantees that .
4. Stabilizing/Destabilizing Effect of the Throughflow on the Onset of Instability
In this section we investigate for the influence of on the instability thresholds (26) and (29). We say that the throughflow has a stabilizing effect if it inhibits the onset of instability (i.e., increases and ). Vice versa, we say that the throughflow has a destabilizing effect.
Lemma 5. If or the throughflow has a stabilizing effect on the onset of steady instability. If the throughflow has a destabilizing effect on the onset of steady instability.
Lemma 6. If (38) or the throughflow has a stabilizing effect on the onset of oscillatory instability. If the throughflow has a destabilizing effect on the onset of oscillatory instability.
In this paper we analyze a vertical fluid motion of a binary mixture saturating a horizontal porous layer with large pores, uniformly heated from below and salted by one salt (either from above or from below). In particular,(1)the definitely boundedness of solutions (existence of absorbing sets) has been recalled;(2)the critical Rayleigh numbers at which steady or oscillatory instability occurs have been recovered;(3)sufficient conditions guaranteeing that a steady or oscillatory secondary motion sets in after the loss of stability have been found.(4)the onset of “cold” instability, possible only when the layer is salted from above, has been analyzed;(5)the stabilizing/destabilizing effect of the vertical throughflow has been investigated.
We conclude this section by showing some numerical simulations on the performed analysis. Figure 1 shows the behavior of (solid line) and (dashed line) as functions of with the other parameters set as =1, =4.8, =80000, =3.2, =. With this choice . According to (43), and decrease with . From Figure 1 it turns out that there exists a critical value of , , such that the motion arising when the throughflow is no longer observable is oscillatory or steady according to or , respectively. This result is in agreement with the results performed in Section 3 since and guarantee, respectively, that (35) and (37) hold. Finally, in order to investigate the influence of the Darcy number (and hence of the Brinkman viscosity), let us observe that the critical wave number given in (25) is a decreasing function of , while given in (24) is an increasing function of (see Figure 2). Hence, as one is expected, and , respectively, given by (26) and (29), are increasing functions of ; i.e., the Brinkman viscosity has a stabilizing effect (see Figure 3).
No data were used to support this study.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
This paper has been performed under the auspices of the G.N.F.M. of INdAM.
- D. A. Nield and A. Bejan, Convection in Porous Media, Springer, New York, NY, USA, 2017.
- S. Rionero and L. Vergori, “Long-time behaviour of fluid motions in porous media according to the Brinkman model,” Acta Mechanica, vol. 210, no. 3-4, pp. 221–240, 2010.
- A. M. Elaiw, A. A. Bakr, M. A. Alghamdi, and F. S. Ibrahim, “Effect of variable viscosity on vortex instability of non-darcy mixed convection boundary layer flow adjacent to a nonisothermal horizontal surface in a porous medium,” Mathematical Problems in Engineering, vol. 2012, Article ID 691802, 14 pages, 2012.
- F. Ali, N. A. Sheikh, M. Saqib, and I. Khan, “Unsteady MHD flow of second-grade fluid over an oscillating vertical plate with isothermal temperature in a porous medium with heat and mass transfer by using the Laplace transform technique,” Journal of Porous Media, vol. 20, no. 8, pp. 671–690, 2017.
- F. Ali, M. Gohar, I. Khan, N. A. Sheikh, S. A. A. Jan, and M. Saqib, “Magnetite molybdenum disulphide nanofluid of grade two: a generalized model with caputo-fabrizio derivative,” Microfluidics and Nanofluidics, vol. 10, 2018.
- M. Saqib, F. Ali, I. Khan, N. A. Sheikh, and A. Khan, “Entropy generation in different types of fractionalized nanofluids,” Arabian Journal for Science and Engineering, pp. 1–10, 2018.
- F. Ali, N. A. Sheikh, I. Khan, and M. Saqib, “Influence of a porous medium on the hydromagnetic free convection flow of micropolar fluid with radiative heat flux,” Journal of Porous Media, vol. 21, no. 2, pp. 123–144, 2018.
- M. Saqib, I. Khan, and S. Shafie, “Application of Atangana-Baleanu fractional derivative to MHD channel flow of CMC-based-CNTs nanofluid through a porous medium,” Chaos, Solitons & Fractals, vol. 116, pp. 79–85, 2018.
- M. Saqib, F. Ali, I. Khan, N. A. Sheikh, and S. B. Shafie, “Convection in ethylene glycol-based molybdenum disulfide nanofluid: Atangana–Baleanu fractional derivative approach,” Journal of Thermal Analysis and Calorimetry, pp. 1–10, 2018.
- N. A. Sheikh, F. Ali, I. Khan, M. Saqib, and A. Khan, “MHD flow of micropolar fluid over an oscillating vertical plate embedded in porous media with constant temperature and concentration,” Mathematical Problems in Engineering, vol. 2017, Article ID 9402964, 20 pages, 2017.
- S. Rionero and I. Torcicollo, “On the dynamics of a nonlinear reaction–diffusion duopoly model,” International Journal of Non-Linear Mechanics, vol. 99, pp. 105–111, 2018.
- S. Rionero and I. Torcicollo, “Stability of a continuous reaction-diffusion Cournot-Kopel duopoly game model,” Acta Applicandae Mathematicae, vol. 132, pp. 505–513, 2014.
- P. Kiran, “Throughflow and non-uniform heating effects on double diffusive oscillatory convection in a porous medium,” Ain Shams Engineering Journal, vol. 7, no. 1, pp. 453–462, 2016.
- C. Zhao, B. E. Hobbs, and H. B. Mühlhaus, “Theoretical and numerical analyses of convective instability in porous media with upward throughflow,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 23, no. 7, pp. 629–646, 1999.
- D. A. Nield and A. V. Kuznetsov, “Onset of convection in a porous medium with strong vertical throughflow,” Transport in Porous Media, vol. 90, no. 3, pp. 883–888, 2011.
- D. A. Nield, “Throughflow effects in the Rayleigh-Benard convective instability problem,” Journal of Fluid Mechanics, vol. 185, pp. 353–360, 1987.
- D. A. Nield and A. V. Kuznetsov, “The effect of vertical throughflow on the onset of convection in a porous medium in a rectangular box,” Transport in Porous Media, vol. 90, no. 3, pp. 993–1000, 2011.
- Y. N. Murty, “Effect of throughflow and magnetic field on Bénard convection in micropolar fluids,” Acta Mechanica, vol. 150, no. 1-2, pp. 11–21, 2001.
- A. A. Hill, S. Rionero, and B. Straughan, “Global stability for penetrative convection with throughflow in a porous material,” IMA Journal of Applied Mathematics, vol. 72, no. 5, pp. 635–643, 2007.
- D. A. Nield, “Convection in a porous medium with inclined temperature gradient and vertical throughflow,” International Journal of Heat and Mass Transfer, vol. 41, no. 1, pp. 241–243, 1998.
- R. De Luca, “Global nonlinear stability and “cold convection instability” of non-constant porous throughflows, 2D in vertical planes,” Ricerche di Matematica, vol. 64, no. 1, pp. 99–113, 2015.
- F. Capone and R. De Luca, “Coincidence between linear and global nonlinear stability of non-constant throughflows via the Rionero “auxiliary system method”,” Meccanica, vol. 49, no. 9, pp. 2025–2036, 2014.
- F. Capone and R. De Luca, “On the stability-instability of vertical throughflows in double diffusive mixtures saturating rotating porous layers with large pores,” Ricerche di Matematica, vol. 63, no. 1, pp. 119–148, 2014.
- F. Capone, R. De Luca, and I. Torcicollo, “Longtime behavior of vertical throughflows for binary mixtures in porous layers,” International Journal of Non-Linear Mechanics, vol. 52, pp. 1–7, 2013.
- F. Capone, V. De Cataldis, R. De Luca, and I. Torcicollo, “On the stability of vertical constant throughflows for binary mixtures in porous layers,” International Journal of Non-Linear Mechanics, vol. 59, pp. 1–8, 2014.
- R. Temam, Infinite-Dimensional Dynamical Systems in Mechanics and Physics, vol. 68 of Applied Mathematical Sciences, Springer, New York, NY, USA, 2nd edition, 1997.
- S. Rionero, “'Cold convection’ in porous layers salted from above,” Meccanica, vol. 49, no. 9, pp. 2061–2068, 2014.
Copyright © 2019 Florinda Capone et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.