Chemical Reaction and Internal Heating Effects on the Double Diffusive Convection in Porous Membrane Enclosures Soaked with Maxwell Fluid

In this paper, the joint impact of the interior heating and chemical reaction on the double diffusive convective flow in porous membrane enclosures soaked by a non-Newtonian Maxwell fluid is investigated applying linear and nonlinear stability techniques. The porous enclosures are square, slender and rectangular. Using the linear stability analysis, the expression for the critical thermal Rayleigh–Darcy number, above which the convective movement occurs, is derived analytically in terms of associated physical parameters. A nonlinear stability examination reliant on the Fourier double series is executed to calculate the convective heat and mass transports of the arrangement. It is observed that the pattern of convective activity is oscillatory only in the occurrence of a relaxation parameter and the threshold value of the relaxation parameter for the occurrence of the oscillatory pattern depends on the other physical parameters. The onset of convective instability accelerates with the increasing chemical reacting parameter, the interior heating parameter, the solute Rayleigh–Darcy number, the Lewis number, the Vadasz number, and the relaxation parameter, while it delays with the heat capacity ratio. The convective heat and mass transfers increase with the solute Rayleigh–Darcy number, the Vadasz number, the relaxation parameter, and the aspect ratio (for rectangular enclosure), while it decreases with the heat capacity ratio and the aspect ratio (for slender enclosure). Additionally, the convective heat transfer enhances with the interior heating parameter, while the convective mass transfer enhances with the chemical reacting parameter and the Lewis number. The effects of Vadasz number, heat capacity ratio, and relaxation parameter are witnessed only on the oscillatory pattern of convection and unsteady convective heat and mass transfers. Further, some existing literature results are compared with the current findings.


Introduction
Double diffusive convective motion is encouraged by two components whose densities are different. These density gradients result in two different diffusion rates which are very crucial for this type of convection. Convection that occurs due to concentration and temperature gradients has gained popularity in previous years because of its realistic applications in science and engineering. Some of the vital sectors of relevances in engineering contain foodstuff and chemical processes [1][2][3][4], casting of metals [5][6][7], crystal growth [8], petroleum production [9][10][11], biomechanical and geosciences applications [12][13][14][15]. The problem of the double diffusive convective progress in a porous layer was explored by Poulikakos [16]. He derived the boundaries describing the areas of direct and overstable styles in terms of critical parameters. The anisotropic effect on the double diffusive convective movement in a porous surrounding substance was explored by Gaikwad et al. [17]. They considered the impacts of anisotropy parameters, solute Rayleigh number, the Soret, and the Dufour factors on the stationary and oscillatory convections as well as on the warmth and mass carrying. Kuznetsov and Nield [18] investigated the heterogeneity effect on the beginning of double diffusive convective movement analytically exploiting linear stability theory. The impact of throughflow on the onset of double-diffusive motion in a permeable layer was explored by Kiran [19] and Shivakumara and Khalili [20]. They found that the consequence of throughflow has either to become stable or to unsettle the arrangement. It depends on the direction of throughflow. Javaheri et al. [21] deliberated the double diffusive convective flow as a purpose for the geological congregation of carbon dioxide. The exponential discrepancy of boundary conditions on the unsteady double diffusive natural convective progress inside the porous enclosures was investigated by Al-Mudhaf [22]. Altawallbeh et al. [23] calculated the power of internal heating on binary instability in a permeable layer. They obtained the critical Rayleigh number and wave number for the stationary and oscillatory mode of convections exploiting the linear hypothesis. Malashetty and Biradar [24] examined the consequence of chemical feedback on the binary diffusive instability in a porous medium layer. They found that the chemical feedback may have a stabilizing or destabilizing influence on the stability criterion. An outstanding review of the studies associated with the double-diffusive instability has been provided by Nield and Bejan [13], Yadav [25], Mojtabi and Charrier-Mojtabi [26], and Vafai [27].
It is known that fluids that take place in the majority of the above-mentioned applications and in nature exhibit non-Newtonian fluid features [28][29][30][31][32][33][34][35][36][37][38]. The problem of double-diffusive instability in a porous medium layer drenched by non-Newtonian fluid taking an Oldroyd model was explored by Malashetty and Swamy [39]. They derived the arrival conditions for stationary, oscillatory, and finite amplitude convective motions analytically. They found a contrast amid the progression of thermal transmission, solute transmission, and viscoelasticity to establish the convective motion as an oscillatory pattern. Malashetty et al. [40] and Kumar and Bhadauria [41] extended this problem with the thermal non-equilibrium effect. The linear stability examination of a Maxwell liquid with double-diffusive motion was undertaken by Awad et al. [42]. They demonstrate that the critical Darcy-Rayleigh number reduces with the relaxation time. Wang and Tan [43] inspected the dual-diffusive instability for non-Newtonian liquid in an absorbent medium considering an amalgamation of Maxwell and Darcy models. They observed that the Soret outcome and relaxation time sped up the onset of convection. The consequence of interior heating on the beginning of double diffusive instability in a coupled stress non-Newtonian liquid flooded porous layer was scrutinized by Gaikwad and Kouser [44]. They found that the arrival of both stationary and oscillatory convective movement is increased by the interior Rayleigh number. Gaikwad and Dhanraj [45] scrutinized the combined weight of anisotropic and interior heating on the binary flow in a non-Newtonian Maxwell liquid flooded permeable layer. They observed that the inner Rayleigh number, mechanical anisotropy factor, and relaxation parameter sped up the start of convective activity, while the thermal anisotropy factor delayed it. Very recently, the impact of larger frequency pulsation in the gravity force on the double-diffusive convective activity with non-Newtonian viscoelastic liquid-filled porous matrix was inspected by Zhao et al. [46].
From the literature inspection, clearly no effort has been made to determine the mutual impact of chemical response and interior heating on the beginning of double diffusive convective motion in a permeable layer flooded by a non-Newtonian fluid. However, there are several practical applications in enhanced oil recovery systems (such as during polymer-flooding practices), a packed bed reactor, chemical privilege equipment, food manufacturing, metal casting procedures, and geophysical arrangements where the porous material may offer its source of warmth and the chemical reaction can take place among the chemical species in the porous substance and the non-Newtonian fluid. Therefore, the present effort intends to inspect the mutual effect of the interior heating and chemical reaction on the double diffusive convective motion in permeable enclosures flooded by a non-Newtonian fluid applying linear and nonlinear stability techniques. To model the non-Newtonian behavior of the fluid, the Maxwell model is used. This is a rate variety of non-Newtonian fluid models in which stress relaxation is known. The Maxwell fluid model acceptably describes the flow behavior of non-Newtonian fluids consisting of a substructure, for instance lubricants with polymer additions, electro-rheological fluids, liquid crystals, blood, and suspension fluids [47][48][49][50]. This work is presented as follows: In Section 1, an introduction of the problem under investigation is provided. In Section 2, the mathematical formulation of the problem is presented. The perturbation equations are obtained in Section 3. In Section 4, the conditions for the start of binary convective flow are derived. The convective heat and mass transports are derived in Section 5. Section 6 presents key results and discussion. At last, this work is completed with a conclusion in Section 7.

Mathematical Formulation
The system examined is a non-Newtonian Maxwell fluid-saturated porous cavity of length H x and width H z with a third dimension infinitely extended so that the fluid flow and heat and mass transport can be taken as two dimensional [29,51], as demonstrated in Figure 1.
diffusive convective motion in a permeable layer flooded by a non-Newtonian fluid. However, there are several practical applications in enhanced oil recovery systems (such as during polymer-flooding practices), a packed bed reactor, chemical privilege equipment, food manufacturing, metal casting procedures, and geophysical arrangements where the porous material may offer its source of warmth and the chemical reaction can take place among the chemical species in the porous substance and the non-Newtonian fluid. Therefore, the present effort intends to inspect the mutual effect of the interior heating and chemical reaction on the double diffusive convective motion in permeable enclosures flooded by a non-Newtonian fluid applying linear and nonlinear stability techniques. To model the non-Newtonian behavior of the fluid, the Maxwell model is used. This is a rate variety of non-Newtonian fluid models in which stress relaxation is known. The Maxwell fluid model acceptably describes the flow behavior of non-Newtonian fluids consisting of a substructure, for instance lubricants with polymer additions, electro-rheological fluids, liquid crystals, blood, and suspension fluids [47][48][49][50]. This work is presented as follows: In Section 1, an introduction of the problem under investigation is provided. In Section 2, the mathematical formulation of the problem is presented. The perturbation equations are obtained in Section 3. In Section 4, the conditions for the start of binary convective flow are derived. The convective heat and mass transports are derived in Section 5. Section 6 presents key results and discussion. At last, this work is completed with a conclusion in Section 7.

Mathematical Formulation
The system examined is a non-Newtonian Maxwell fluid-saturated porous cavity of length x H and width z H with a third dimension infinitely extended so that the fluid flow and heat and mass transport can be taken as two dimensional [29,51], as demonstrated in Figure 1. It is assumed that the porous cavity is subjected to an internal heat supply of strength S and first-order chemical response of rate K R among a chemical genus in the permeable matrix and the Maxwell fluid. The temperatures θ L and θ U , and solute concentrations φ L and φ U , are consistently forced on the bottom and top boundary walls such that θ L > θ U Membranes 2022, 12, 338 4 of 26 and φ L > φ U , while other boundary walls are accepted to be insulated and impermeable to mass transmission. It is specified that the porous cavity is homogeneous and in local thermal balance with Maxwell fluid.
The continuity equation for the Maxwell fluid in a porous cavity is [52][53][54]: Here, V D is Darcy's velocity of the Maxwell fluid, ∇ = i x ∂ ∂x + i z ∂ ∂z and, i x and i z are unit vectors in x and z paths.
The momentum equation by taking the Darcy rule for the porous medium, the Maxwell model for non-Newtonian fluid, and the Boussinesq approximation for density disparity with temperature and solute concentration in the buoyancy force is [55][56][57]: where, τ represents the time, θ symbolizes the temperature, φ represents the solute concentration, ρ 0 symbolizes the density at θ U , K represents the permeability of the porous medium, γ 1 symbolizes the stress relaxation feature time constant, µ symbolizes the viscosity of Maxwell fluid, P symbolizes the pressure, β θ and β φ symbolize the thermal and solute expansion coefficients, respectively. The energy equation for a Maxwell fluid in a heat-generating permeable cavity is [31,58,59]: Here, S symbolizes the power of the interior heating, k E symbolizes the effectual thermal conductivity of the porous medium, (ρc) and (ρc) E symbolize the heat capacities of the Maxwell fluid and effectual permeable medium, correspondingly.
The conservation equation for the solute with the chemical reaction of rate K R in the porous matrix is [60] ∂ ∂τ Here, D S is the solutal diffusivity and ε is the porosity of the permeable medium. On eliminating the pressure term and using the stream function χ as u = −∂χ/∂z and w = ∂χ/∂x, the Equations (1)-(4) can be written as: ∂φ ∂τ The boundary situations are: For no-dimensional examination, the dimensionless variables are described as: where α E = k E (ρc) . Then, the dimensionless forms of Equations (5)-(8) become: Here,

Basic Condition
For a time-free calm solution of Equations (10)- (12), it is assumed that the temperature and solute distributions for the basic solution are: On solving the Equations (11)-(13) for the basic solution, we have: In the lack of interior heating and chemical reaction, Equations (15) and (16) give: Equation (17) is the one found by Kuznetsov and Nield [61] for the case of pure fluid.

Perturbation Equation
Now, perturbation on the basic condition is imposed as: where,χ ,θ , andφ are the perturbed variables on their basic estimates. Replacing Equation (18) into Equations (10)-(13), we have:

Linear Stability Consideration
In this segment, the thresholds of the marginal and oscillatory type of convective flows are determined applying linear theory. Now, it is supposed that the amplitudes of the perturbation are extremely small and expressed as [13,29,60,62]: where λ and σ represent the wavenumber and enlargement rate of disturbances, correspondingly. On using Equation (23) in Equations (19)- (22) and avoiding the nonlinear terms with perturbed variables, we have: Here, d dz ≡ D. To find an approximate solution to the system of Equations (24)- (27), the Galerkin routine is utilized [63][64][65][66][67][68]. The trial functions (fulfilling the boundary circumstances) are picked as: where, E, F, and G are unidentified coefficients. On applying Equation (28) into Equations (24)- (26) and for the non-singular solution, we have: where, . Now, from Equation (25), we have: where,

Marginal Pattern of Convection
The marginal pattern of convection can happen, if σ = 0. Thus, Equation (30) gives the marginal thermal Rayleigh-Darcy number R M DT as: The threshold of the R M DT signifying the beginning of marginal convection occurs at λ c where λ c = √ a satisfies the equation: In the nonattendance of chemical response ( K RN → 0 ), Equations (33) and (34) turn into: For the case of a single component ( R DS → 0 ), Equations (35) and (36) give: The Equations (37) and (38) are the same as found by Yadav and Maqhusi [29] for the lack of viscosity variation.
For the case of one component ( R DS → 0 ), Equations (39) and (40) offer the critical marginal thermal Rayleigh-Darcy number R M DT,c as: This coincides with the conclusion of Nield and Kuznetsov [70] for the convection in a rectangular box. Additionally, Equation (41) agrees with the experimental results obtained by Horton and Rogers [71] and Katto and Masuoka [72].

Oscillatory Pattern of Convection
The oscillatory pattern of convective motion occurs when σ = 0 and N 2 = 0. Then, Equation (30) proposes the oscillatory thermal Rayleigh-Darcy number R Os DT as: Additionally, from Equation (30), the rate of oscillation σ satisfies the following dispersion relation: Here, From Equation (43), it is clear that the oscillatory manner of convective flow is conceivable only if σ 2 > 0. The analytical appearance for R Os DT specified by Equation (42) is minimized respecting the wave number λ numerically when σ 2 > 0 for diverse estimates of involved physical parameters to identify their impacts on the start of the oscillatory type of convective motion.

Weak Nonlinear Stability Investigation
The linear stability inquiry offers the threshold for the start of convective movement of the Maxwell fluid in terms of R DT,c but does not calculate the convective heat and mass transport. To get these extra details, here the weak nonlinear stability analysis is used as [17,44,56,60,73]: where, B 11 (τ), C 11 (τ), C 02 (τ), D 11 (τ), and D 02 (τ) are undetermined amplitudes and to be determined. On replacing Equations (44)- (46) into Equations (19)- (21), we obtain the nonlinear stability equations as: where . The above nonlinear equations are not fit for analytical inspection for the time-reliant variables. Thus, we solved it numerically utilizing the Runge-Kutta-Fehlberg technique (RKF45). The results are also validated with ODE45 solver in MATLAB. For the initial state, we select B 11 = 1, C 11 = 0, C 02 = 0, D 11 = 0, and D 02 = 0.

Convective Heat and Mass Transports
Convective heat and mass transfers play a very vital role in detecting the convective motion in more early stages. Heat and mass transfers can be calculated in spans of Nusselt number Nu and Sherwood number Sh, individually, and described as [17,44,56,60,73]: On applying Equations (15), (16), (45) and (46) into Equations (62) and (63), we have:

Results and Discussion
The impact of interior heating and chemical reaction on the beginning of dual diffusive convective motion and the convective heat and mass transfers in non-Newtonian Maxwell fluid-saturated permeable square (A = 1), slender (A < 1), and rectangular (A > 1) enclosures were explored. Employing the linear stability philosophy, the criteria for the start of the marginal and oscillatory pattern of convective motions were derived analytically in terms of R DT,c , which is the function of physical parameters K RN   Rayleigh-Darcy number DS R (Figure 3iii), the aspect ratio A (Figure 3iv tion parameter  (Figure 4i), the Lewis number Le (Figure 4ii), the heat m (Figure 4iii), and the Vadasz number Va (Figure 4iv). From these plo that the neutral curves are linked in a topological way. This establishes t stability of the arrangement is specified in the span of , DT c R , which at low   Figure 2i, it is recognized that the power of the basic temperature allocation increases with accumulating interior heating parameters S N and the profile of basic temperature allocation alters from linear to nonlinear with S N . This happened because the increasing S N provides more warming to the Maxwell fluid layer, which enhances the strength of buoyancy force and as a result, more disturbances are seen in the system. From Figure 2ii, we observed a similar result on the allocation of basic state solute concentration φ b with increasing K RN and Le. This happened because an increase in the strength of K RN and Le creates more disturbances in the system.    (Figure 3i), the interior heating parameter S N (Figure 3ii), the solute Rayleigh-Darcy number R DS (Figure 3iii), the aspect ratio A (Figure 3iv), the relaxation parameter γ (Figure 4i), the Lewis number Le (Figure 4ii), the heat capacity ratio m (Figure 4iii), and the Vadasz number Va (Figure 4iv). From these plots, it is found that the neutral curves are linked in a topological way. This establishes that the linear stability of the arrangement is specified in the span of R DT,c , which at lower values the arrangement is stable and at R DT somewhat greater than R DT,c , convective activity starts. From these figures, it is also established that by increasing the values of K RN , S N ,R DS , γ, Le, and Va, the estimate of R DT,c tends to lessen, i.e., the arrangement goes to destabilize, while m has a stabilizing impact on the stability of the structure. The aspect ratio A does not affect R DT,c . The marginal pattern of the convective motion is found to be free with m, Va, and γ. Figure 5 illustrates the variations in R DT,c , λ c , and σ as a function of γ for varied estimates of the chemical reacting parameter K RN and the interior heating parameter S N . From Figure 5i,ii, we found that with a boost in the values of K RN , S N , and γ, the critical thermal Rayleigh-Darcy number R DT,c diminishes. This shows that the chemical reacting parameter K RN , the interior heating parameter S N , and the relaxation parameter γ speed up the beginning of convective activities. This is due to the fact that increasing the chemical reacting parameter K RN and the interior heating parameter S N creates more disturbances by increasing the energy supply to the system. Further, the destabilizing effect of the relaxation parameter γ is found because thermal diffusivity of the system increases with γ (from the definition of γ). From Figure 5iii,iv, it is noted that the critical wave number λ c decreases with K RN , S N , and γ. This illustrates that the magnitude of convection cells increases with K RN , S N , and γ. From Figure 5v,vi, we detected that the frequency of oscillations σ decreases with accumulating K RN and S N , while an opposite result is seen with γ. Figure 6 exhibits the impact of R DS and Le on the stability of the scheme. From Figure 6i,ii, it is noticed that an improvement in the estimation of R DS and Le is to speed up the marginal and oscillatory patterns of convective motions.    This is on the ground that the disturbance to the arrangement increases with DS R and Le . Furthermore, an increase in the estimate of Le increases the threshold estimate of  at which the pattern of instability is amended. From Figure 6iii,iv, it is found that the critical wave number  c declines with DS R and Le for a marginal pattern of convection, while a reverse result is seen for an oscillatory pattern of convection. From This is on the ground that the disturbance to the arrangement increases with R DS and Le. Furthermore, an increase in the estimate of Le increases the threshold estimate of γ at which the pattern of instability is amended. From Figure 6iii,iv, it is found that the critical wave number λ c declines with R DS and Le for a marginal pattern of convection, while a reverse result is seen for an oscillatory pattern of convection. From Figure 6v,vi, it is established that the frequency of oscillations σ shrinks with R DS and Le. Le .  The impacts of the heat capacity ratio m and the Vadasz number Va on R DT,c , λ c , and σ are exposed in Figure 7. From Figure 7i,ii, it is found that R DT,c surges with m for oscillatory pattern of convection, while this result is opposite with Va. This appears that the outcome of m delays the start of an oscillatory pattern of convection. This is due to fact that the energy restoring capacity of arrangement enhances with escalating the heat capacity ratio m.
Further, the Vadasz number Va advances the start of the oscillatory type of movement. A similar result of the Vadasz number Va on the system was also observed by Kumar and Bhadauria [41] and Malashetty and Biradar [76]. From Figure 7iii,iv, it is proven that the critical wave number λ c increases with m and Va. The frequency of oscillations σ diminishes with m, while this result is opposite with Va as found from Figure 7v,vi. From Figure 7, it is also noticed that increasing Va decreases the threshold estimate of γ at which the pattern of instability alters, while this result is reversed with m.  Further, the Vadasz number Va advances the start of the oscillatory type of movement. A similar result of the Vadasz number Va on the system was also observed by Kumar and Bhadauria [41] and Malashetty and Biradar [76]. From Figure 7iii,iv, it is proven that the critical wave number  c increases with m and Va . The frequency of oscillations  diminishes with m , while this result is opposite with Va as found from Figure 7v,vi. From Figure 7, it is also noticed that increasing Va decreases the threshold estimate of  at which the pattern of instability alters, while this result is reversed with  Figure 8 demonstrates the power of the aspect ratio A on the stability of the system. From this graph, it is found that λ c reduces with escalating A. This illustrates that the dimension of convective cells enhances with increasing the aspect ratio A. From Figure 8, it is also found that A has no control on R DT,c and σ. that the dimension of convective cells enhances with increasing the aspect ratio A Figure 8, it is also found that A has no control on , DT c R and  .   To observe the effect of K RN , S N , R DS , Le, and A on steady-state heat and mass spreads, the Nusselt number Nu and Sherwood number Sh are plotted in Figure 9 as a function R DT for diverse values of these parameters. From this figure, we recognize that if R DT increases from one to five or six times of R DT,c , the heat and mass transport increased significantly, and if R DT increases more, it remains moderately constant. Membranes 2022, 12, x FOR PEER REVIEW 23 of 30  The convective steady mass transfer increased with increasing K RN , R DS , Le, and A for the slender enclosure (A < 1), while it decreased with A for rectangular enclosure A. From Figure 9, it is also found that increases in the values of S N , R DS , and A for the slender enclosure (A < 1) amplify the convective heat transport in the scheme, while for a rectangular enclosure (A > 1), convective heat transport decreases with A.
A for the slender enclosure ( ) , while it decreased with A for rectangular enclosure ( )  1 A . From Figure 9, it is also found that increases in the values of N S , DS R , and A for the slender enclosure ( )  1 A amplify the convective heat transport in the scheme, while for a rectangular enclosure ( ) , convective heat transport decreases with A .  The nonlinear unsteady ordinary differential Equations (47)-(51) are solved numerically by applying the RKF45 method with realistic initial circumstances. The achieved outcomes are presented in Figures 10 and 11. It is noticed that the greatest increase in Nusselt number Nu and Sherwood number Sh appear near to the opening time; it reveals the spatial progress of increasing frequency. Lastly, the oscillations reach the steady situation for a sufficiently large amount of time. From Figures 10 and 11, it is also found that increases in the estimates of K RN , S N , R DS , R DT , Le, Va,A for the slender enclosure (A < 1) and γ enhance the convective unsteady mass transmission in the arrangement, while it decreases with A for rectangular enclosure (A > 1) and m. The convective unsteady heat transportation in the structure increases with increasing S N , R DS , R DT , Va,A for the slender enclosure (A < 1) and γ, whereas it reduces with A for the rectangular enclosure (A > 1) and m. Table 1. Comparison of R DT,c , λ c , and σ for different values of γ, K RN , and S N at Le = 2, m = 1.5, A = 1, R DS = 3, S N = 0.5, K RN = 0.5, and Va = 3.   In Figure 12, a comparison is made between the obtained results and the experimental results offered by Serkitjis [77] for normal fluid. From this figure, it is noticed that the heat spread observed by the experiment are very close to our outcomes if R DT ≥ R DT,c , whereas for R DT < R DT,c , the experimental results are slightly lower, here R DT,c = 4π 2 .
Nusselt number Nu and Sherwood number Sh appear near to the opening time; it reveals the spatial progress of increasing frequency. Lastly, the oscillations reach the steady situation for a sufficiently large amount of time. From Figures 10 and 11, it is also found that increases in the estimates of  In Figure 12, a comparison is made between the obtained results and the experimental results offered by Serkitjis [77] for normal fluid. From this figure, it is noticed that the heat spread observed by the experiment are very close to our outcomes if

Conclusions
The double diffusive convective progress in a Maxwell fluid occupying the interior heat-generating and chemical reacting porous enclosures was investigated analytically as well as numerically taking linear and weak nonlinear stability theories. The consid-

Conclusions
The double diffusive convective progress in a Maxwell fluid occupying the interior heat-generating and chemical reacting porous enclosures was investigated analytically as well as numerically taking linear and weak nonlinear stability theories. The considered enclosures were rectangular (A > 1), square (A = 1), and slender (A < 1). The impact of various parameters on the onset of the convective flow, and on the convective heat and mass transports of the system were achieved. The important conclusions are as follows:

•
Increasing the chemical reacting parameter K RN , the interior heating parameter S N , the solute Rayleigh-Darcy number R DS , the relaxation parameter γ, the Lewis number Le, and the Vadasz number Va accelerates the onset of double diffusive convective motion, while it delays with increasing the heat capacity ratio m.

•
The dimension of convective cells enhances by increasing the aspect ratio A, the chemical reacting parameter K RN , the interior heating parameter S N , and the relaxation parameter γ, while it decreases with the heat capacity ratio m.

•
Increasing R DS and Le enhances the size of marginal convective cells, while this result is opposite for oscillatory convection.

•
The convective mass transfer in the system is augmented with increasing K RN , S N , R DS , R DT , Le, Va, A (for the slender enclosure) and γ, while it reduces with A (for rectangular enclosure) and m.

•
The convective heat transport in the system is enhanced with increasing S N , R DS , R DT , Va, A (for the slender enclosure) and γ, whereas it reduces with A (for rectangular enclosure) and m.

•
The marginal pattern of the convective motion and steady heat and mass transport are observed to be free with m, Va, and γ.