# Suppression of Chaos in Porous Media Convection under Multifrequency Gravitational Modulation

Authors: Karam Allali
Journal: Advances in Mathematical Physics
 (2018)

---

## Abstract

Suppression of chaos in porous media convection under multifrequency gravitational modulation is investigated in this paper. For this purpose, a two-dimensional rectangular fluid-saturated porous layer heated from below subjected to a vertical gravitational modulation will be considered. The model consists of nonlinear heat equation coupled with a system of equations describing the motion under Darcy law. The time-dependent gravitational modulation is assumed to be with two frequenciesσ1 and σ2. A spectral method of solution is used in order to reduce the problem to a system of four ordinary differential equations. The system is solved numerically by using the fifth- and a sixth-order Runge-Kutta-Verner method. Oscillating and chaotic convection regimes are observed. It was shown that chaos can be suppressed by appropriate tuning of the frequencies’ ratio η=σ2/σ1.


---

## Body

## 1. Introduction

Several studies have been devoted to investigating the effect of a periodic gravitational modulation on the convective instability in a fluid layer. Indeed, the modulation of gravity can lead to stabilizing or destabilizing effect on the dynamics and the convective properties of the fluid [1, 2]. In practice, the gravitational modulation can be achieved by oscillating, in the vertical direction, the fluid container that is already subjected to a constant gravitational field. The modulation of gravity, with some specific amplitudes and frequencies, can contribute to the improvement of product performances such as solidification process [3, 4], frontal polymerization [5], or crystal growth [6].The effect of gravitational modulation on natural convection in a porous layer was studied by Govender [7]. The linear stability analysis method was used to reduce the problem to a Mathieu equation and it was shown that increasing the frequency of vibration stabilizes the fluid convection. Elhajjar et al. [8] studied the influence of small-amplitude and high-frequency vertical vibrations on the Soret-driven convection flow. Both the direct numerical simulations and linear stability analysis were used in order to show that vibrations delay the transition from unicellular to bicellular flow regime. More recently, Vadasz et al. [9, 10] studied the periodic and chaotic natural convection in a porous layer subject to vertical vibrations. It was shown that periodic and chaotic solutions alternate when Rayleigh number increases. All those previous works consider the effect of one frequency gravitational modulation. However, adding another frequency or more can change considerably the onset of convection regime. For instance, it was shown that the gravitational modulation with two incommensurate frequencies produces a convective stabilizing or destabilizing effect depending on the frequencies’ ratio [11, 12]. In the absence of gravitational modulation, it was shown that chaos and thermal convection instabilities can be controlled or even suppressed by tuning of the temperature boundary values [13, 14].Motivated by the recent works of Vadasz et al. [9, 10], the aim of this present paper is to continue the investigation of the influence of gravitational modulation on convective instability in porous layer heated from below, by assuming that the modulation is with two frequencies. To this end, we will consider that the system containing the fluid and the porous matrix are subjected to two-frequency vertical gravitational modulation. This external excitation causes a time-dependent acceleration given by g+b(t), where g is the gravity acceleration and b(t)=λ1sin⁡ν1t+λ2sin⁡(ν2t), where λ1 and λ2 and ν1 and ν2 are the amplitudes and the frequencies of the vibration, respectively. We will examine the effect of the two frequencies’ ratio on the convective instability of the fluid flow and whether the chaos regime can be suppressed by tuning of the second frequency. Indeed, in many situations, chaos presents a troublesome phenomenon that may lead to physical system damage, thermal explosion due to irregular temperature oscillations, deadly epilepsy, and cardiac arrhythmias [15–18]. Accordingly, the need for chaos suppression is very important. Among those situations, one will be presented in this paper and we will study the effectiveness of introducing a second frequency gravitational modulation in chaos suppression. The need to add such second frequency gravitational modulation is particularly interesting when the first frequency is constrained to take some specific values or belongs to a narrow interval for which the chaos persists. In this paper, we will study the importance of adding a second frequency gravitational modulation in controlling chaos for both problems dealing with low or high Prandtl number fluids.The paper is organized as follows: the next section introduces the model, while Section3 deals with the numerical method used in our simulations. Results and discussions are provided in Section 4. The last section concludes the work.


## 2. Mathematical Formulation

### 2.1. The Mathematical Model

We consider a two-dimensional rectangular fluid-saturated porous layer subjected to a vertical gravitational modulation, as shown in Figure1. We assume that the fluid is incompressible and is heated from bellow. The model of such a process can be described by the energy equation coupled with the hydrodynamics equations under the Boussinesq-Darcy approximation:(1)∂T∂t+v·∇T=κΔT,∂v∂t+μKv+∇p=gβT-TC1+λ1sin⁡ν1t+λ2sin⁡ν2tγ,∇·v=0,where T denotes the temperature, v=(u,v) is the velocity, p is the pressure, κ is the coefficient of thermal diffusivity, μ is the kinematic viscosity, g is the gravity, γ is a unit vector in the downward vertical direction (in the sense of gravity), TC is the cold wall temperature, and K is the permeability. We consider the following boundary conditions:(2)x=0,l×l,∂T∂x=u=0,y=0,T=TC,v=0,y=l,T=TH,v=0;here TH is the hot wall temperature; l×l and l are the horizontal and the vertical dimensions of the porous medium. The boundary conditions for the temperature correspond to the adiabatic (no-flux) condition at the lateral boundaries and fixed temperature at the lower and upper boundaries. For the velocity, its normal component at the boundary is zero. This means that the fluid does not intersect the boundary. In what follows, we will denote by DTc=TH-TC the characteristic temperature difference.Figure 1
A two-dimensional rectangular porous domain heated from below subjected to two-frequency gravitational modulation.


### 2.2. The Dimensionless Model

To obtain the dimensionless model, we introduce new spatial variablesx′=x/l and y′=y/l, time t′=κ/l2t, velocity l/κv, pressure K/μκp, and the frequencies σi=l2/κνi (i=1,2). Denoting θ=T-TC/DTc and keeping for convenience the same notation for the other variables, we obtain(3)∂θ∂t+u∂θ∂x+v∂θ∂y=∂2θ∂x2+∂2θ∂y2,(4)χ∂u∂t+u+∂p∂x=0,(5)χ∂v∂t+v+∂p∂y=-Raθ1+λ1sin⁡σ1t+λ2sin⁡σ2t,(6)∂u∂x+∂v∂y=0.Here Ra=gβDTcKl/κμ is the Rayleigh number and (u,v) is the velocity vector. The parameter χ=1/PrD stands for the inverse of Darcy-Prandtl number with PrD=Pr/Da; Pr=μ/κ is the Prandtl number and Da=K/l2 is the Darcy number.The system of (3)–(6) is supplemented by the free surface boundary conditions for the velocity, zero boundary condition for the temperature at the upper and lower boundaries of the rectangular domain, and adiabatic boundary condition at the side boundaries:(7)x=0,l,∂θ∂x=u=0,y=0,θ=0,v=0,y=1,θ=1,v=0.


## 2.1. The Mathematical Model

We consider a two-dimensional rectangular fluid-saturated porous layer subjected to a vertical gravitational modulation, as shown in Figure1. We assume that the fluid is incompressible and is heated from bellow. The model of such a process can be described by the energy equation coupled with the hydrodynamics equations under the Boussinesq-Darcy approximation:(1)∂T∂t+v·∇T=κΔT,∂v∂t+μKv+∇p=gβT-TC1+λ1sin⁡ν1t+λ2sin⁡ν2tγ,∇·v=0,where T denotes the temperature, v=(u,v) is the velocity, p is the pressure, κ is the coefficient of thermal diffusivity, μ is the kinematic viscosity, g is the gravity, γ is a unit vector in the downward vertical direction (in the sense of gravity), TC is the cold wall temperature, and K is the permeability. We consider the following boundary conditions:(2)x=0,l×l,∂T∂x=u=0,y=0,T=TC,v=0,y=l,T=TH,v=0;here TH is the hot wall temperature; l×l and l are the horizontal and the vertical dimensions of the porous medium. The boundary conditions for the temperature correspond to the adiabatic (no-flux) condition at the lateral boundaries and fixed temperature at the lower and upper boundaries. For the velocity, its normal component at the boundary is zero. This means that the fluid does not intersect the boundary. In what follows, we will denote by DTc=TH-TC the characteristic temperature difference.Figure 1
A two-dimensional rectangular porous domain heated from below subjected to two-frequency gravitational modulation.


## 2.2. The Dimensionless Model

To obtain the dimensionless model, we introduce new spatial variablesx′=x/l and y′=y/l, time t′=κ/l2t, velocity l/κv, pressure K/μκp, and the frequencies σi=l2/κνi (i=1,2). Denoting θ=T-TC/DTc and keeping for convenience the same notation for the other variables, we obtain(3)∂θ∂t+u∂θ∂x+v∂θ∂y=∂2θ∂x2+∂2θ∂y2,(4)χ∂u∂t+u+∂p∂x=0,(5)χ∂v∂t+v+∂p∂y=-Raθ1+λ1sin⁡σ1t+λ2sin⁡σ2t,(6)∂u∂x+∂v∂y=0.Here Ra=gβDTcKl/κμ is the Rayleigh number and (u,v) is the velocity vector. The parameter χ=1/PrD stands for the inverse of Darcy-Prandtl number with PrD=Pr/Da; Pr=μ/κ is the Prandtl number and Da=K/l2 is the Darcy number.The system of (3)–(6) is supplemented by the free surface boundary conditions for the velocity, zero boundary condition for the temperature at the upper and lower boundaries of the rectangular domain, and adiabatic boundary condition at the side boundaries:(7)x=0,l,∂θ∂x=u=0,y=0,θ=0,v=0,y=1,θ=1,v=0.


## 3. The Spectral Method of Solution

In order to perform the mathematical analysis of our problem, we first introduce, due to the incompressibility property of the fluid, the stream functionψ defined by(8)
u
v
=
∂
ψ
∂
y
-
∂
ψ
∂
x
.We apply the rotational operator to the system of (4)-(5) in order to eliminate the pressure; the system becomes(9)∂θ∂t+∂ψ∂y∂θ∂x-∂ψ∂x∂θ∂y=∂2θ∂x2+∂2θ∂y2,(10)χ∂Δψ∂t+Δψ+Ra∂θ∂x1+λ1sin⁡σ1t+λ2sin⁡σ2t=0.To apply the spectral method [19, 20], we will look for a solution of the problem in the form of the basic solution added to a variation one. More precisely, we will separate the stream function and temperature into a basic conduction part and variation convection one in the form(11)ψ=ψB+ψV=ψV,θ=θB+θV=1-y+θV,where ψB and θB are the basic solution of problem (9)-(10). The variation part of the solution is chosen as follows: (12)ψV=A11sin⁡πxlsin⁡πy,θV=B11cos⁡πxlsin⁡πy+B02sin⁡2πy.These two representations are equivalent to a Galerkin expansion of the solution in x and y directions, truncated when (i+j)=2, where i and j are the Galerkin summation indices in x-direction and in y-direction, respectively. By substituting the form of solutions (11) in (10), we will have(13)-χddt+1π21l2+1A11sin⁡πxlsin⁡πy-Ra1+λ1sin⁡σ1t+λ2sin⁡σ2tπlB11sin⁡πxlsin⁡πy=0;by setting (14)ϱ=l2l2+1,we will have(15)π2χϱdA11dt+π2ϱA11+πRal1+λ1sin⁡σ1t+λ2sin⁡σ2tB11sin⁡πxlsin⁡πy=0.Again, by substituting the form of solutions (11) in (9), we will have(16)dB11dtcos⁡πxlsin⁡πy+dB02dtsin⁡2πy-A11B11π2lsin2πxlsin⁡πycos⁡πy+A11πlcos⁡πxlsin⁡πy-A11B11π2lcos2πxlsin⁡πycos⁡πy-A11B022π2lcos⁡πxlsin⁡πycos⁡2πy=-B11π2cos⁡πxlcos⁡πy-4B02π2sin⁡2πy-B11π2l2cos⁡πxlsin⁡πy,which means(17)dB11dtcos⁡πxlsin⁡πy+dB02dtsin⁡2πy-A11B11π22lsin⁡2πy+B11π2ϱcos⁡πxlsin⁡πy-A11B02π2lcos⁡πxlsin⁡3πy+A11πlcos⁡πxlsin⁡πy+4B02π2sin⁡2πy+A11B02π2lcos⁡πxlsin⁡πy=0;then(18)dB11dt+π2ϱB11+πϱA11+π2ϱA11B02cos⁡πxlsin⁡πy+dB02dt-π22lA11B11+4π2B02sin⁡2πy-π2lcos⁡πxlsin⁡3πyA11B02=0.By multiplying (15) by sin⁡(πx/l)sin⁡(πy) and integrating over the space domain (0,1)×(0,l), we have (19)l4π2χϱdA11dt+π2ϱA11+πRal1+λ1sin⁡σ1t+λ2sin⁡σ2tB11=0.Again, multiplying (18) by cos⁡(πx/l)sin⁡(πy) and integrating over the space domain (0,1)×(0,l), we will have (20)l4dB11dt+π2ϱB11+πϱA11+π2ϱA11B02=0.Also, multiplying (18) by sin⁡(2πy) and integrating over the space domain (0,1)×(0,l), it implies that (21)l2dB02dt-π22lA11B11+4π2B02=0,so (22)dA11dt=-1χA11+Raπς1+λ1sin⁡σ1t+λ2sin⁡σ2tB11,dB11dt=-π2l2+1l21πςA11+B11+1ςA11B02,dB02dt=-π2l2+1l24l2l2+1B02-12ςA11B11.The first equilibrium point of this system is (23)A11f=0,B11f=0,B02f=0,which is expected to be stable when Ra/π2ς2<1. However, for Ra/π2ς2>1, the following second equilibrium point is expected to be stable: (24)A11∗=±l8Ra-8π2ς2π2l2+π2=±22ϱRaπ1-π2ς2Ra,B11∗=±lπςRa8Ra-8π2ς2π2l2+π2=±22ϱςRa1-π2ς2Ra,B02∗=-Ra+π2ς2πRa=1ππ2ς2Ra-1.By assuming that (25)1R=π2ς2Ra,we will have (26)A11∗=±22ϱςR-1,B11∗=±22ϱπR-1R,B02∗=-1πRR-1.We set (27)X=-2ϱςR-1-1A11∗,Y=22ϱπR-1R-1B11∗,Z=-1πRR-1-1B02∗;we obtain (28)X˙=-ϱχπ2X-1+λ1sin⁡σ1t+λ2sin⁡σ2tY,Y˙=-Y+RX-R-1XZ,Z˙=4ϱXY-Z.By setting (29)α=ϱχπ2,W=σ1t,we will obtain the following system describing the problem dynamics:(30)X˙=-αX-Y+αλ1sin⁡W+λ2sin⁡ηWY,Y˙=-Y+RX-R-1XZ,Z˙=4ϱXY-Z,W˙=σ1,with the two modulation frequencies’ ratio given by (31)η=σ2σ1.


## 4. Numerical Simulations

The system of (30) will be solved numerically using the fifth- and a sixth-order Runge-Kutta-Verner method. The used initial conditions are given as follows: (32)X0=Y0=Z0=0.9,W0=0,and we will use the following parameters for all the numerical simulations: (33)ϱ=0.5,λ1=1,λ2=1.In order to compare our results with those in the work of Vadasz et al. [9], we will use the following parameter describing natural convection with low Prandtl number, that is, α=5. Meanwhile, to compare our results with the recent work taking into account high Prandtl number [10], we will fix α=500. The other control parameters R, σ1, and η will be varied adequately in order to check the impact of the two-frequency vibration on the behaviour of solution.Figure2 shows the projection of the solution on X-Y and X-Z planes for α=5, R=25, σ1=0, and η=0. These parameters correspond to the absence of gravitational modulation case. The plots show chaotic behaviour of the solution, which are similar to those obtained in Vadasz et al. [9]. Figure 3 shows the behaviour of solution when the porous layer is subjected to one-frequency gravitational modulation; all the parameters in this figure are similar to the previous one, except σ1=4.5. From this figure, it can be seen that the chaos still persists and the one-frequency gravitational modulation has no effect in controlling the chaos, since the leading Lyapunov exponent remains positive. Figure 4 shows a periodic convective regime, which means that suppression of chaos can be gained by changing one control parameter and making it nonzero: η=6. One can conclude that the second gravitational modulation has a significant effect in chaos suppression. This result is confirmed by plotting the dynamics of the problem’s Lyapunov exponents (see Figure 5). From this figure, we clearly see that all the problem’s Lyapunov exponents are negatives, which support the periodicity behaviour of the convective regime. Hence, for the problems dealing with low-Prandtl number fluids, one can expect the suppression of chaos by introducing a second gravitational modulation frequency.Figure 2
The solution on a phase diagram forα=5, R=25, σ1=0, and η=0. Y versus X (a) and Z versus X (b).
(a) (b)Figure 3
The solution on a phase diagram forα=5, R=25, σ1=4.5, and η=0. Y versus X (a) and Z versus X (b).
(a) (b)Figure 4
The solution on a phase diagram forα=5, R=25, σ1=4.5, and η=6. Y versus X (a) and Z versus X (b).
(a) (b)Figure 5
The dynamics of the problem’s Lyapunov exponents forα=5, R=25, σ1=4.5, and η=6.Figure6 shows two projections from the phase diagram on both X-Y and X-Z planes for α=500, R=75, σ1=6.2, and η=0. This case corresponds to one-frequency gravitational modulation and the plots are similar to those obtained in Vadasz et al. [10]; this case corresponds to high-Prandtl number fluids. It can be seen from these two figures that the convective regime is chaotic. Figure 7 shows the behaviour of solution when the second modulation is introduced: η=3; the chaotic regime is completely suppressed for a second frequency greater only three times than the first one. It was found from the Lyapunov exponents’ dynamics (not shown) that all the exponents are negatives, which confirms the chaos suppression. One can conclude that suppression of chaos can be achieved by looking only for small natural integer frequencies’ ratio, which is more convenient to be reproduced experimentally.Figure 6
The solution on a phase diagram forα=500, R=75, σ1=6.2, and η=0. Y versus X (a) and Z versus X (b).
(a) (b)Figure 7
The solution on a phase diagram forα=500, R=75, σ1=6.2, and η=3. Y versus X (a) and Z versus X (b).
(a) (b)


## 5. Conclusion

In this paper, we have studied the effect of multifrequency gravitational modulation on chaos suppression in the problem dealing with natural convection in porous layer heated from below. For this purpose, we have assumed that the system containing the fluid and the porous matrix are subjected to two-frequency vertical gravitational modulation. The considered model consists of heat equation coupled with the hydrodynamics equations under Darcy law. The mathematical analysis of the problem was performed by using the spectral method in order to reduce the problem to a system of four ordinary differential equations describing the dynamics of solution. The fifth- and a sixth-order Runge-Kutta-Verner method was used to solve numerically the reduced system of differential equations. Numerical simulations were conducted for both cases of low- and high-Prandtl number fluids. It was shown that the chaotic regime can be totally suppressed by choosing adequately the second gravitational modulation frequency. This situation happens when the frequencies’ ratio is equal to six in the case of low Prandtl number or is equal to three in the case of high Prandtl number, which can be more convenient to be reproduced experimentally, since one can look only for small natural integer frequencies’ ratio.



---

*Source: 1764182-2018-04-03.xml*