Thermal Convective Instabilities and Chaos in a Rotating Hybrid Nanofluid Layer with Cattaneo–Christov Heat Flux Model

Laboratoire de la Mécanique des Fluides, de la Dynamique Non-Linéaire et de la Modélisation des Systèmes Biologiques (LMFDNMSB), Institut de Mathématiques et de Sciences Physiques (IMSP), Université d’Abomey-Calavi (UAC), Godomey, Benin Département de Physique, FAST-Natitingou, Université Nationale des Sciences, Technologies, Ingénierie et Mathématiques (UNSTIM) d’Abomey, Abomey, Benin Laboratoire de Physique et Applications du Centre Universitaire de Natitingou, Université Nationale des Sciences, Technologiques, Ingénierie et Mathématiques (UNSTIM) d’Abomey, Abomey, Benin Department of Industrial Computer and Electrical Engineering, Université Nationale des Sciences, Technologiques, Ingénierie et Mathématiques (UNSTIM) d’Abomey, Abomey, Benin Département de Physique, ENS-Natitingou, Université Nationale des Sciences, Technologies, Ingénierie et Mathématiques (UNSTIM) d’Abomey, Abomey, Benin


Introduction
In 1995, Choi introduced, at the Argonne National Laboratory of the University of Chicago in the U.S., the concept of nanofluid [1].
is new generation of fluids consists of dispersing nanoparticles (assemblies of a few hundred to a few thousand atoms, leading to an object with at least one dimension of thousands of atoms, leading to an object of which at least one of the dimensions is of nanometric size) in a base liquid (water, oil, ethylene glycol, toluene). e use of these nanofluids in some industrial operations that involve heat transfer by convection is a promising alternative solution to improve thermal performance. us, the problem of natural convection in a nanofluid layer heated from below has been studied by several researchers [2][3][4][5][6] with the Fourier law. e flow and heat transport of nanomaterial with quadratic radiative heat flux and aggregation kinematics of nanoparticles reported by Mahanthesh [7] revealed that the suspension of the nanoparticles increases the thermal conductivity and, thus, improves the temperature and reduces the heat flux at the plate. e Rayleigh-Bénard convection in nanofluid submerged with dust particles was investigated by Shalini and Mahanthesh [8]. ey pointed out that the inclusion of nano and dust particles reduces the Rayleigh number while the rotation postpones the onset of convection and stabilizes the system. Ahuja and Sharma [9] conducted a comprehensive review of the instability of Rayleigh-Bénard convection in nanofluids by summarizing in their work the studies related to the instability of a horizontal nanofluid layer under the impact of various parameters such as rotation magnetic field, Hall currents, and LTNE (local thermal nonequilibrium) effects in porous and nonporous media. e thermal convection in a rotating fluid layer provides a system to study hydrodynamic instabilities, pattern formation, and spatiotemporal chaos in nonlinear dynamical systems with many practical applications in engineering, such as food processing, chemical processes, solidification, centrifugal casting of metals, and rotating machines [10].
To compensate for the defect and all the disadvantages of mono nanofluids, it is essential to combine several nanoparticles to prepare a hybrid nanofluid [11]. Natural magnetohydrodynamic convection in a triangular cavity filled with a hybrid (copper-alumina)/water nanofluid with localized heating from below and internal heat has been investigated by Rashad et al. [12]. ey concluded that a hybrid nanofluid composed of equal amounts of copper and alumina nanoparticles dispersed in water has no significant improvement on the average Nusselt number compared to the mono nanofluid and that the effect of increasing the hybrid nanoparticles becomes significant in the case where natural convection is very low. Aladin et al. [13] also studied the significant effects of suction and magnetic field on a moving plate containing a hybrid (copper-alumina)/water nanofluid. ey have proved that the hybrid nanofluid gives better results than the mono nanofluid.
According to Myson and Mahanthesh [14], the hybrid nanofluid delays the convection and will further enhance the heat transfer rate, but the Casson parameter advances the convection while reducing the heat transfer rate. Mackolil and Mahanthesh [15] illustrated the optimization of heat transfer in the thermal Marangoni and nonlinear convective flow of a hybrid nanomaterial with sensitivity analysis. It is shown that the hybrid nanomaterial possesses enhanced thermal fields for nanoparticle volume fractions less than 0.02. e sensitivity computation of nonlinear convective heat transfer in hybrid nanomaterial between two concentric cylinders with irregular heat sources was also made by riveni and Mahanthesh [16].
Given its advantages and industrial applications, especially in chemical reactions, biological systems, crystal production, petroleum reservoir modeling, and packed-bed catalytic filtration, chaotic convection in a hybrid nanofluid layer should receive considerable attention due to the performance of nanofluids. Jawdat et al. [17], Moaddy et al. [18,19], and Bhardwaj and Chawla [20] all contributed well to nonlinear dynamical analysis of the thermal convection in a horizontal nanofluid layer heated from below in the presence or absence of a magnetic field. ey studied the effect of nanoparticles on chaotic convection in a layer of fluid (water) heated from below and noticed that the stability region can be increased by using nanofluids and that the onset of chaotic convection can be delayed under the influence of nanoparticles. Also, variations in temperature and magnetic field strength cause the system to transition from a steady state to chaos and back to a steady state. e case of hybrid nanofluid was first presented by Dèdèwanou et al. [21] with the Fourier law. ey found that the copper nanofluid makes it possible to quickly switch from chaotic to periodic regimes compared to the alumina nanofluid, and the use of hybrid nanoparticles allows further control of the chaos in the base fluid by expanding the convective flow and reducing the chaotic flow.
Furthermore, Maxwell and Cattaneo modified Fourier's law by taking into account the aspect of thermal relaxation time in the propagation of heat [22]. In order to eliminate the heat flow and thus obtain a single equation for the temperature field, Christov [23] proposed a generalization of the material-invariant Maxwell-Cattaneo law, in which the relaxation time of the heat flow is given by the convex Oldroyd upper derivative.
is new law was used by Straughan [24] to study thermal convection in an ordinary fluid. He concluded that the thermal relaxation time is significant if the Cattaneo number is large enough, and the convection mechanism changes from stationary to oscillatory convection with narrower cells. Indeed, some researchers used the Cattaneo-Christov model to appreciate the effects of temporal relaxation on the thermal behavior of a nanofluid [25][26][27][28][29][30][31]. Alebraheem and Ramzan [26] have studied the heat and mass transfer of Casson nanofluid flow containing gyrotactic microorganisms past a swirling cylinder by considering the Cattaneo-Christov heat flux model. According to their numerical solution of the subject system, which is framed via the bvp4c technique of MATLAB software, the concentration of the fluid is reduced owing to the increase in values of the brownian motion parameter and local Reynolds number, but the diminishing density of microorganisms is perceived for mounting estimates of the bioconvection Péclet number. Multiple slip impacts in the MHD hybrid nanofluid flow with Cattaneo-Christov heat flux and autocatalytic chemical reaction were investigated by Gul et al. [32]. ey found that the fluid temperature is diminishing function of the thermal slips parameters but increased for high estimates of the heat source-sink and nanoparticle volume concentration parameters while entropy number augmented for higher thermal relaxation parameter and Reynolds number. Lu et al. [33] have also studied a thin film flow of nanofluid comprising carbon 2 Complexity nanotubes influenced by Cattaneo-Christov heat flux and entropy generation. ey showed that the velocity and temperature distributions increase as the solid volume fraction escalates. Recently, a three-dimensional flow of gold-silver/engine oil hybrid nanofluid owing to a rotating disk of variable thickness with Cattaneo-Christov heat flux has been addressed by Zhang et al. [34]. ey proved that the performance of the hybrid nanofluid is far better than the common nanofluid according to the surface temperature and heat transfer rate. is same remark is made from the results obtained for the model-based comparative study of magnetohydrodynamics unsteady hybrid nanofluid flow between two infinite parallel plates with particle shape effects [35]. Considering hybrid nanofluid Yamada-Ota and Xue flow models in a rotating channel with the modified Fourier law, it is observed that the velocity profile decreases for the higher rotation parameter while it increases for the escalated slip parameter, but the fluid concentration and temperature are on the decline for higher surface catalyzed reaction and thermal relaxation parameters respectively [36]. Ramdan et al. [37] have analyzed the hydrodynamic and heat characteristics of the three-dimensional flow of a steady, laminar, and incompressible convective graphene-copper oxide/water and graphene-silver/water hybrid nanofluids (used as a solar energy absorber) with varied particle shapes in a porous medium. eir study revealed that the rotational parameter has declined the velocity profiles but enhanced the temperature profiles, and the decline effect is significant in the case of graphene-copper oxide/water whereas the enhancement effect of temperature is significant for graphene-silver/water. A comparative analysis of magnetized partially ionized copper, copper oxide-water, and kerosene oil nanofluid flow with Cattaneo-Christov heat flux was made by Abid et al. [38]. ey noted the greater effective thermal conductivity for copper-water partially ionized nanofluid as compared to other given partially ionized nanofluids (copper-kerosene oil, copper oxide-water/kerosene oil partially ionized nanofluids). Ramzan et al. [39] developed a mathematical model for the nanofluid flow containing carbon nanotubes with ethylene glycol as a base fluid in a rotating channel with an upper permeable wall by adding the Cattaneo-Christov heat flux's impact with thermal stratification. e displacement of the lower plate at variable velocity, caused by the rotation of the fluid, produces forced convection with rotation and centripetal impact. Nevertheless, the upper plate is porous. Chu et al. [40] investigated a numerical solution for MHD Maxwell nanofluid with gyrotactic microorganisms, a higher-order chemical reaction in the presence of variable source/sink, and Newtonian heating in a rotating flow on a deformable surface and noted that on incrementing the conjugate heat parameter and thermal relaxation time, the rate of heat transfer augments, but the rate of heat transfer decreases on varying the fluid relaxation time.
Despite all the above, the study of nanoparticles requires more attention due to their industrial uses. After inspecting the scientific literature, we noted that no work has yet addressed the chaotic aspect of thermal convection in hybrid nanofluids, taking into account the thermal relaxation time, although this would be very useful in some applications like petroleum reservoir modeling, chemical reactions, thermal transport in biological tissue, and surgical operations. Nevertheless, Layek and Pati [41] studied the effects of thermal lag on the onset of convection, its bifurcations, and the chaos of a horizontal layer of the heated Boussinesq fluid underneath via a five-dimensional nonlinear system. A comparative study of the five-dimensional system obtained for the case of a hybrid nanofluid was made by Dèdèwanou et al. [42]. erefore, the objective of the present paper is to investigate the effects of hybrid nanoparticles on the occurrence of thermal convection instability and chaos in a rotating fluid layer heated from below with the Cattaneo-Christov heat flux model. A specific objective is to determine the analytical expression for the stationary Rayleigh number that can be used to study the nonlinear dynamics of thermal convection in rotating hybrid nanofluid flow in the presence of thermal relaxation time.
is work aims to study the different transition regimes as a function of the thermophysical properties of nanofluids and then to show the effects of hybrid nanoparticles, Taylor number, and Cattaneo number on the chaotic behavior of natural convection in a basic fluid such as water via dynamical systems.
In the next section, the thermal convection in an infinite horizontal rotating hybrid nanofluid layer with the hyperbolic Cattaneo-Christov heat flux is outlined. Section 3 discusses the theory of conduction, stationary convection, and oscillatory convection, where we generalize and simplify the expression for the Rayleigh number by deriving a number of new analytical results. In order to reduce the set of equations governing the dynamic behavior of thermal convection in the hybrid (alumina-copper)/water nanofluid, discretized models in four and six dimensions are developed in Section 4 using the Galerkin expansion. We have studied the nature of the nonlinear dynamics of the obtained dynamical systems and determined the fixed points by analyzing the stability of the stationary solutions. ese analyses have allowed us to justify the influence of the hybrid alumina-copper nanoparticle, the Cattaneo number, and the Taylor number on number on the transition from chaos to periodicity and vice versa in the fluid. In Section 5, we present the different simulations performed, and the results obtained are discussed. e conclusions are drawn in Section 6.

Problem Formulation.
We consider an infinite rectangular cavity with two horizontal walls maintained at different temperatures. is cavity, heated from below with a thermal relaxation time, is filled with a hybrid nanofluid (water and nanoparticles), subject to gravity acting downwards and to rotation. In order to develop our numerical model, it is necessary to adopt certain assumptions, namely, the flow is assumed to be permanent and incompressible, the mixture is assumed to be homogeneous, single-phase, Newtonian, the nanoparticles are spherical, and the mass transfer between the particles and the fluid is negligible. e Cartesian coordinate system used is such that the y-axis follows the horizontal and the vertical z-axis is collinear with gravity. e geometry of our problem is presented in Figure 1.

Governing Equations.
In this section, we have studied the equations governing the dynamic and thermal fluxes with boundary conditions and nondimensional parameters characterizing the thermal convection in a rotating hybrid nanofluid layer in the presence of thermal relaxation time.
Taking into account the listed assumptions and using the hybrid nanofluid model proposed in Section 3, the equations governing the conservation of mass, momentum, energy, and heat flux for a laminar flow of the hybrid nanofluid are written in their dimensional form respectively as follows [10,42]: In equation (2), e i � (0, 0, 1) is the unit vector and ∇ 2 � (z 2 /zx 2 ) + (z 2 /zy 2 ) + (z 2 /zz 2 ) is the Laplacian operator. e use of the Boussinesq approximation allows us to define the density as a function of temperature as e density, thermal expansion coefficient, heat capacity, dynamic viscosity, thermal conductivity of the dynamic viscosity, and thermal conductivity of the hybrid nanofluid, respectively, are defined as follows: to where e thermophysical properties of nanoparticles and water used in this work are summarized in Table 1 (see [13]).

Linear Stability Analysis
We consider a classical Rayleigh-Bénard problem of linear stability of convective rolls in a horizontal fluid layer with unconstrained boundary conditions. us, the temperature boundary conditions are T � T c at z � 0 and T � T 0 at z � 1 with T c > T 0 . As for the velocity, its component along the z axis is zero at the boundaries.

Steady-State Solutions.
A time-independent quiescent solution of equations (1)-(4) with temperature and heat flux varying in the z direction only, is obtained by reducing equations (2)-(4) to en, the steady state solutions are given with the temperature gradient χ and the profile of P 0 defined by

Simplified Set of Equations.
In order to simplify the parametric representation of this physical problem and to find the characteristic properties of the system, it is necessary to recast the flow equations. us, the following normalized quantities are introduced: Considering small perturbations on the basic solutions as follows: and neglecting the products of the primed quantities, we obtain the following dimensional equations: where the dimensionless parameters are defined by Now we eliminate the pressure from the nondimensional equations (14)-(17) by taking the curl-curl of equation (15), the divergence of equation (15), the inner product of any vector equation with e i ;, and denoting the divergence of the heat flux Q � (zQ i /zx i ), to obtain, after dropping the tilde notation for brevity, the equations which describe the evolution of the conduction steady state perturbations in a conveniently simplified form with four variables such as the z-component of the velocity field, the vorticity ζ � zv/zx − zu/zy, the heat flux Q, and the temperature T. An evolutionary equation for the vorticity can be obtained from the equation of motion by taking curl, then the dot product with e 3 for the vertical component.
us, eliminating the pressure and introducing the vorticity in equation (2) allows us to obtain the following: From equations (19) and (22), we obtain the following: Complexity

Normal Modes and Analytical Solution.
e linear stability of the conduction solutions is studied by writing the perturbations in separable form and assuming an exponential time dependence with the plane tiling function satisfying where W, Θ, Φ are eigenfunctions. e substitution of equation (24) into the differential equations (20), (21), and (23) leads to where D � d/dt, D 2 � d 2 /dt 2 and Equations (26)- (28), which represent the starting point for analytical and numerical calculations on thermal convective instability are used to study the occurrence of stationary and oscillatory convection in nanofluids. ey are equivalent to those obtained by Straughan [24] and Bissell [43,44] in the case of an ordinary fluid not subjected to a Coriolis force, i.e. when φ 1 � φ 2 � 0 and Ta � 0. e disappearance of tangent shear stresses at the free surface and the conservation of the mass equation allow us to obtain the boundary conditions of the free surface defined by In order to obtain an approximate solution of equations (26)-(28), we used the Galerkin weighted residual method by choosing the test function written as which fulfill the conditions at the borders mentioned in equation (30).
By substituting the test functions defined in equation (31) into equations (26)- (28) and performing some integrations, we obtain the following matrix equation: where J � D 2 − κ 2 . For this matrix equation (32) to admit a nontrivial solution, the Rayleigh number of the base fluid must be in the following form: 6 Complexity For the following, let σ � jω with j 2 � − 1 and ω the real frequency.
us, the expression for the Rayleigh number defined in equation (33) becomes with Since the Rayleigh number is a real and positive physical quantity, then for the expression equation (33) to exist, ω must be zero or Δ 2 � 0.

Stationary Convection.
According to the stability exchange principle for the stationary case, the stability margin is characterized by the frequency equal to zero. is condition allows to obtain from the expression equation (33), the Rayleigh number of the base fluid of the stationary convection expressed as follows: If c 1 � c 2 � c 3 � c 4 � 1 and the cavity is not rotating, i.e. φ 1 � φ 2 � 0 and Ta � 0, equation (37) is equivalent to the classical Rayleigh number of stationary convection in ordinary fluids [42].
In the absence of the rotation force, we have e absolute critical Rayleigh number of the hybrid nanofluid in this case is defined as with the corresponding wavenumber We note that the Rayleigh number obtained for stationary convection in nanofluids is not a function of the Prandtl number or the Cattaneo number of the base fluid. us, the same results can be associated with the more usual Fourier law for mono nanofluids or hybrid nanofluids.

Oscillatory
Convection. Now, we study the effects of the Cattaneo number, the nanoparticles, and the rotation on the oscillating convection. In this case, we must have ω ≠ 0 and Δ 2 � 0. erefore, the Rayleigh number of the base fluid for oscillatory convection is given by e corresponding oscillatory frequency ω must verify the following equations: For oscillatory instability to be possible, the value of ω 2 generated by equation (42) must be positive.

Dynamical System Analysis
In order to explore how thermal relaxation time, hybrid nanoparticles, and rotation affect the nonlinear stability of the onset of thermal convection in a horizontal layer of ordinary fluid such as water, we reduce the problem to the classical case of two-dimensional convective rolls in a fluid layer with unconstrained horizontal boundaries. us, we assume that all physical quantities are independent of x. We consider the early stages of nonlinear convection when the basic structure of the convective rolls is still determined by the dynamic behavior of the linearized solution. e real components of the fluid velocity are expressed in terms of partial derivatives of the stream function as follows: Eliminating the pressure from equation (2) and introducing the expressions of the stream function defined in equations (43) into the resulting equation and equations (3) and (4), we obtain with the appropriate dimensionless variables [10,42], the following new equations: 1 e ratio of the Rayleigh number of the hybrid nanofluid to that of the base fluid gives the effective Rayleigh number of the hybrid nanofluid as a function of the thermophysical properties and Rayleigh number of the heat transfer fluid defined as Similarly, the Cattaneo number of the hybrid nanofluid is defined as follows:

Reduced Set of Equations.
e solution of the coupled nonlinear system of partial differential equations (44) and (45) will be obtained by representing the current function and the temperature using the Galerkin expansion in the following form [42]: is representation is equivalent to a Galerkin expansion of the solution in the y and z directions, truncated when i + j � 2, where i is the Galerkin summation index in the y direction and j is the Galerkin summation index in the z direction. Substituting equations (49) into equations (44) and (45), multiplying the equations by the orthogonal eigenfunctions corresponding to equations (44), and integrating over the domain and wavelength of the convection cell in the vertical and horizontal directions respectively, i.e., π/κ 0 dy 1 0 dz(.), we obtain a set of three ordinary differential 8 Complexity equations for the time evolution of the second-order amplitudes expressed by After introducing new variables of amplitudes defined as and the expressions in equations 50, we obtain the following system: erefore, we can reduce the amplitude equations of system (53) to a system of first-order nonlinear equations by introducing the amplitudes us, we obtain the six-dimensional system by describing the nonlinear dynamic behavior of thermal convection in mono or hybrid nanofluids, presented as follows: where the dot (.) denote the time derivative d/dt. When φ 1 � φ 2 � 0, φ 1 ≠ 0 and φ 2 � 0, φ 1 � 0 and φ 2 � 0, system equation (63) corresponds to the base fluid, aluminawater nanofluid, copper-water nanofluid, respectively.
In the absence of the thermal relaxation time, (54) and (46) are reduced to the following system: Lorenz [45] has investigated the nonlinear analysis of convection in pure fluid confined in a nonporous cavity by using the Fourier law. His nonlinear dynamic system has been analyzed and solved for Pr f � 10, so that there are convection cells in the domain and that the boundary conditions are satisfied [45,46]. Bissell [43] analyzed the oscillatory convection with the Cattaneo-Christov hyperbolic heat-flow model and included the effects owing to Prandtl number, which in some circumstances can be used as a control parameter. He showed that the Cattaneo threshold can be conceived equivalently as a Prandtl threshold, so that system bifurcations could potentially be triggered by varying the Prandtl number. For small values of Cattaneo number, a five-dimensional nonlinear system obtained by Layek and Pati [41] undergoes a subcritical transition to chaos similar to the Lorenz system but undergoes a period-doubling transition to chaos when Pr f � 5 and C f � 0.001. For increasing values of Prandtl number, he found that the fine-structure of the period-doubling cascade is interrupted and that this is due to the generation of internal noise that fastens the transitional process. With the critical value of the wavenumber corresponding to the convection threshold, the expression of the Rayleigh number of the base fluid defined in equation (53) gives: Ra fc � 27π 4 ϵ/4.

Dissipation Effect.
e nonlinear dynamical system (54) has the reflection symmetry (U, Y, P) ⟶ − (U, Y, P) and We note that ∇ → . ϑ → < 0 whatever the values of Pr f , c 5 , c 6 and δ. en system (54) is dissipative and its solutions are bounded in phase space. erefore, if a set of initial points in phase space occupies the region ϑ(0) at t � 0, then after some time, t, the end points of the corresponding trajectories will fill a volume (57)

Equilibrium Points and eir Stability.
In this section, the nature of the nonlinear dynamics of systems (54) and (55) is determined around the fixed points by analyzing the stability of stationary solutions. e hybrid nanofluid is confined in a nonporous cavity so that there are convection cells in the domain and the boundary conditions are satisfied.

e
Case of C f � 0. Considering the general form of system (54) defined by _ X � F(X s ) and the equilibrium (stationary or fixed) points X s defined by F(X s ) � 0, we obtained three fixed equilibrium points of the system, including the first one is the stationary solution and the other two are the convection solutions. e linear stability of the points can be obtained by linearizing the nonlinear dynamical system equation (55). us, the resulting Jacobian matrix is as follows: 10 Complexity Using test solutions of the form exp(ξt), ξ ∈ C, the stability of the fixed point corresponding to the conduction solution is controlled by the roots of the following characteristic polynomial equation: is equation (61) generates four eigenvalues. e first one given by ξ � − λc 7 is always negative, and the other three are the solutions of the following equation: From this equation, the fixed point of the stationary solution is stable if and only if c 5 c 6 ϵ R f < c 7 (T f + c 2 5 ). us, the critical value of the rescaled Rayleigh number of the base fluid at which the fixed point of the stationary solution loses its stability and that of the convection solution takes over is expressed by is expression is a function of the thermophysical properties of the hybrid nanofluid and the rescaled Taylor number, so the transition from conduction to stationary convection depends on the volume fraction of the nanoparticles and the effect of rotation as shown in Figure 2.
Data analysis of the curves constructed in Figure 2 shows that when the value of Taylor number is less than about 0.33, 0.315, 0.30, and 0.293 for φ 1 � φ 2 � 0.01, 0.02, 0.03, and 0.04, respectively, R fc1 decreases but increases for higher values of T f . Taking T f � 0 for example, we found R fc1 � 1.2 like Gupta [10] for the ordinary fluids (φ � 0). But when φ 1 � φ 2 � 0.01, 0.02, 0.03, 0.04; we have R fc1 ≈ 1.176, 1.157, 1.142, 1.13. us, it is then possible to reduce or increase conduction in a heat transfer fluid using hybrid nanoparticles under the effect of rotation. Using the same test solutions of the form exp(ξt), ξ ∈ C, the stability of the fixed points corresponding to the convection solution is controlled by the roots of the following characteristic polynomial equation: is equation is solved numerically for different values of the parameters to study the stability of the fixed points of the convection solutions.

e
Case of C f ≠ 0. In the presence of the thermal relaxation time, the elimination of a quadratic factor, which is not associated with the beginning of the instability, Complexity 11 allowed us to find the characteristic polynomial equation of the fixed point corresponding to the immobile support at the origin whose roots control its stability, which is presented as follows: (65) When ξ � 0, a stability exchange occurs and stationary convection takes over. e corresponding critical rescaled Rayleigh number of the base fluid from which this phenomenon is observed is equivalent to equation (63).

Results and Discussion
We performed numerical simulations to investigate the influence of hybrid nanoparticles and rotation on the dynamic behavior of thermal convection in a base fluid (water) in the presence of thermal relaxation time. Using free boundary conditions, we determined the analytical expressions of Rayleigh numbers of the base fluid for stationary and oscillatory convection as a function of the thermophysical properties of the hybrid nanofluid. We observe that the stationary Rayleigh number of the base fluid does not depend on the Prandtl number and the Cattaneo number. Figure 3 shows the variation of the stationary Rayleigh number of the base fluid as a function of wavelength for different values of the volume fraction of the hybrid nanoparticles (alumina-copper) with a fixed value of Taylor number. From these plotted curves, we find that the stationary Rayleigh number increases with the value of the volume fraction of hybrid nanoparticles. us, the addition of the hybrid nanoparticles (alumina-copper) to the base fluid (water) subjected to the rotation stabilizes the stationary convection. Figure 4 shows the variation of the stationary Rayleigh number of the base fluid as a function of the wavelength for different values of the Taylor number with a fixed value of the volume fraction of the hybrid nanoparticles (alumina-copper).
From these plotted curves, it can be seen that the stationary Rayleigh number increases with an increasing Taylor number.
us, the rotation stabilizes the stationary convection in the hybrid nanofluid. e fourth-order Runge-Kutta method, the polynomial companion matrix, and the standard eigenvalue solver of the Lapack method are used to numerically solve systems equations (54) and (55). We took the initial conditions U(0) � Y(0) � 0.8, Z(0) � 0.92195, P(0) � 0.8, S(0) � 0.92195 and V(0) � 0.8. In order to guarantee the results, our different numerical simulations are compared with the results obtained by Dèdèwanou et al. [42] and Gupta [10]. We present in Figures 5-8, the bifurcation diagrams representing the minima and maxima of the posttransient regimes of the solutions of the amplitude Z(t) as a function of function of R f when thethermal relaxation time is zero using Pr f � 10 and λ � 8/3. ese diagrams show that system equation (55) can have chaotic, periodic, or multiperiodic behavior depending on the parameter values chosen. By comparing the diagrams in Figures 5-7, we notice that, for T f � 0.2, when the volume fraction of the hybrid nanoparticles increases, the      To confirm this prediction of the bifurcation diagrams, we constructed in Figure 9, the chaotic behavior of the system in the base fluid case by choosing R f � 80 and T f � 0.5. As shown in Figure 10, we set the values of the rescaled Taylor and Rayleigh numbers by varying the volume fraction of the hybrid nanoparticles to construct the in-plane phase spaces with their corresponding time evolutions. For φ 1 � φ 2 � 0, 0.02 and 0.04, the system is chaotic in period 2 and period 1, respectively. erefore, by analyzing these curves shown in Figure 10, it can be deduced that the addition of the hybrid nanoparticles in a heat transfer fluid makes the convection periodic. On the other hand, in the case of Figure 11, we have fixed the values of the volume fraction of hybrid nanoparticles and the rescaled Rayleigh number of the base fluid by varying the value of the rescaled Taylor number. e analysis of these curves shows that the system leaves the chaotic regime to the periodic regime when the value of the rescaled Taylor number decreases. erefore, increasing the rescaled Taylor number increases the periodicity of the system. Furthermore , we present in Figures 12-15, the bifurcation diagrams representing the minima and maxima of the posttransient regimes of the solutions of the amplitude Z(t) as a function of R f when the thermal relaxation time exists using Pr f � 5 like Layek and Pati [41]. ese diagrams show that system 63 can also have chaotic, periodic, or multiperiodic behavior depending on the parameter values chosen. Comparing the diagrams in Figures 12-14, it can be seen that, for T f � 0.2 and C f � 0.001, increasing the volume fraction of the nanoparticles hybrid nanoparticles decreases the domain of chaotic behavior with the increase of the Rayleigh number values of the base fluid. On the other hand, the comparison of the diagrams in Figures 13 and 15 show that, for φ 1 � φ 2 � 0.02 and T f � 0.2, increasing the values of the Cattaneo number largely increases the range of chaotic behavior with increasing values of the Rayleigh number of the base fluid. Referring to Figure 14, for φ 1 � φ 2 � 0.04, T f � 0.2 and C f � 0.001, the system is chaotic for 14 < R f < 16.
In Figures 16 and 17, we have constructed the phase spaces in the X − Z plane for different values of the control parameters of system equation (54). When we set T f � 0.2 and C f � 0.001 (see Figure 16), we notice in the base fluid case (φ 1 � φ 2 � 0) that the system is in period 4 for R f � 166. For φ 1 � φ 2 � 0.02, the system is in period 2 and for φ 1 � φ 2 � 0.04, the system has quasi-chaotic behavior. For φ 1 � φ 2 � 0.02 and R f � 100 fixed, the system is in period 1 for T f � 0.2 and C f � 0.001. On the other hand, for T f � 0.2 and C f � 0.003, T f � 0.2 and C f � 0.005, T f � 0.3 and C f � 0.005, T f � 0.5 and C f � 0.001, and T f � 1.7 and C f � 0.005, the system is chaotic.