Abstract

The present paper has enormous applications including fields of crystal growth, material processing, spacecraft, underground spread of chemical contaminants, petroleum reservoirs, waste dispersal and fertilizer migration in saturated soil, alloy solidification, and many more. The importance of double-diffusive convection has been recognized in various engineering applications, and it has thoroughly been investigated theoretically. In the presence of a constant heat source/sink on both layers, the double component-magneto-Marangoni-convection in a composite layer is examined. Due to heat, this combined layer is enclosed by adiabatic boundaries and exposed to basic temperature gradients. For the system of ordinary differential equations, the thermal surface-tension-driven (Marangoni) number, which also happens to be the Eigenvalue, is solved in closed form. For basic temperature gradients, the eigenvalues, or thermal Marangoni numbers are determined in closed form for lower rigid and higher free boundaries with surface tension, depending on both temperature and concentration. The effect of several parameters on the eigenvalue against depth ratio was studied. It is observed that the inverted parabolic temperature gradient on double component magneto-Marangoni-convection in a combined layer is the most stable of the three temperature gradients and for larger values of depth ratios, all physical parameters are nominal for the porous region dominant combined system.

1. Introduction

Owing to its use in geophysics, astrophysics, engineering, and technology, MHD flow research is important. The influence of a magnetic field on the distribution of the heat by a heat source/sink where the liquid can produce and absorb hot radiation is influential when it comes to confined use and high operating temperatures. Hill [1] studied double-diffusive convection in a fluid-saturated porous medium with an internal heat source depending on concentration. Using linear stability analysis, Liu and Umavathi [2] investigated the double-diffusive convection of sparsely packed micropolar fluid-saturated porous media. Using the Runge–Kutta technique, Bhuvanavijaya and Mallikarjuna [3] investigated the effect of nonuniform heat sources on convective heat and mass transport over a vertical plate in a rotating system embedded in a non-Darcy porous medium. In a layered system, Köllner and Boeck [4] demonstrated a two-dimensional simulation of solutal Rayleigh–Bénard–Marangoni convection. Alizadeh et al. [5] investigate the possibility of Marangoni convection causing instability in porous medium statistically. They discovered that porous media’s damping forces postpone the commencement of convection. In the context of the Soret effect, Gaikwad and Dhanraj [6] used linear and weakly nonlinear stability calculations to investigate the beginning of Darcy–Brinkman double-diffusive convection in a binary viscoelastic fluid-saturated porous layer. Kumar et al. [7] explored the topic of triple-diffusive convection and the impact of the solute Rayleigh numbers, concentration Rayleigh number, Lewis number and modified diffusivity ratio on system stability. Mebarek–Oudina [8] investigated the numerical modeling of hydrodynamic stability in a vertical annulus with varying lengths of heat source. He discovered that as the heat source length ratio increases, the critical Rayleigh number decreases. Awasthi et al. [9] conduct a linear stability study for the beginning of triple-diffusive convection in the presence of an internal heat source in a Maxwell fluid-saturated porous layer. It has been discovered that Lewis numbers have a destabilizing impact, whereas solute Rayleigh numbers have a stabilizing effect. Using the shooting technique, Reddy and Dinesh [10] investigated coupled highly nonlinear partial differential equations for the effects of internal heat generation, combined effects of Soret and Dufour, and changeable fluid characteristics. Basha et al. [11] use the Crank-Nicolson technique to investigate the effects of Lorentz force on the fluid transport parameters of a chemically reactive nanofluid with two types of geometries. They discovered that the Hartman number has a significant impact on fluid flow and heat transfer properties. Athira et al. [12] investigated the effects of nonlinear convection and induced magnetic fields in a viscous fluid flow over a porous plate when chemical reactions and heat sources/sinks were present. K. H. Ahmed and W. M. Ahmed [13] used the finite volume approach to investigate the natural convection in a nanofluid including water-copper nanoparticles. Manjunatha and Sumithra [14] studied the triple diffusive convection in the absence of a heat source. Sumithra et al. [15, 16] discuss single components, heat sources, and with/without magnetic fields in the presence of temperature gradients. They developed the closed form of the thermal Marangoni number solution.

Recently, Hemanthkumar et al. [17] looked into the impact of a horizontal pressure gradient on the beginning of Soret-driven thermosolutal porous convection. Using normal modes, Jayalatha and Suma [18] investigated the effect of heat sources on instability in rotating viscoelastic liquids. Supraja and Raju [19] give a numerical investigation of heat transmission by the natural convection of a fluid inside a square cavity with two inner bodies. Using the incompressible smoothed particle hydrodynamics approach, Raizah and Aly [20] investigated the suspension of a nanoencapsulated phase change material during a double-diffusive convection flow inside a compound cavity. Double component convection in rectangular hollows with varied factor ratios is investigated numerically by Yang and Zhao [21]. The effects of magnetic field and heat source on double-diffusive convection two layer system studied by Manjunatha and Sumithra [22]. The effects of waveforms of internal heat source modulation on triple diffusive convection in viscoelastic liquids were investigated by Arshika et al. [23]. Kumar Pundir et al. [24] investigated the linear stability of a fluid layer with nanoparticles in the presence of a solute gradient. Using linear and nonlinear stability approaches, Yadav et al. [25, 26] investigated the combined influence of heating and chemical reaction on the double-diffusive convective in a non-Newtonian Maxwell fluid and Kuvshiniski type saturated porous layer. They determined that oscillatory convection is only possible when the solute Rayleigh–Darcy number is negative and that it is also influenced by other physical factors. Animasaun et al. [27] studied the Ratio of Momentum Diffusivity to Thermal Diffusivity. Sumithra et al. [28] evaluated the role of thermal diffusion on the commencement of double-diffusive convection in a composite system and the influence of the Soret effect on the system’s stability. For adiabatic-adiabatic and adiabatic-isothermal boundary conditions, Manjunatha et al. [29] investigated the effect of a constant heat source/sink on double-diffusive convection in a combined layer system in the absence of a magnetic field. Some interesting and recent papers on nanofluid flows are [3032]. Convective stability analysis has been studied by a few authors for porous layers [3339].

The current paper investigates the effect of fundamental temperature gradients on double component-magneto-Marangoni-convection (DCMC) in a combined layer with a constant heat source (sink) at both levels. The Eigen value problem was solved by using the exact method and the effect of several parameters on the eigenvalue against the depth ratio was studied. Numerous applications in astronomy, engineering, geophysics, climatology, and crystal formation [40] will surely benefit from this work.

2. Research Methodology

Consider a two-component, subject to the Boussinesq approximation, electrically conducting fluid-saturated, incompressible, sparsely packed porous layer of thickness underlying a two-component fluid layer of thickness with an imposed magnetic field intensity in the vertical Z-direction and with heat sources and , respectively. Let and are the temperature and concentration difference between lower and upper boundaries. The constant temperature and concentration & . At the upper free region, surface tension force acts which varies linearly with temperature and concentration respectively in the form , where , , is the unperturbed value. The porous layer’s lower surface is hard, while the fluid layer’s upper surface is free, with surface tension effects depending on temperature and concentration as shown in Figure 1. The rectangular coordinate system’s origin is defined at the interface of the fluid-porous composite system and the magnetic field is imposed along z–axis upwards. The Darcy model predicts that at the interface, the velocity, shear stress, normal stress, heat, heat flux, mass, and mass flux are all continuous.

The basic equations are (see Venkatachalappa et al. [41] and Shivakumara et al. [42]).

2.1. Fluid Layer: Region-1

The fundamental equations are

2.2. Porous Layer: Region-2

where subscript “f” and “p” denotes the fluid and porous layer, are the velocities, is the density, are the viscosities, are the hydromagnetic pressures, are the temperatures, are the thermal diffusivities, , are the heat sources, are the salinity diffusivities, are the concentrations, are the magnetic permeability’s, are the magnetic viscosities, is the magnetic field, is the porosity, is the permeability of the porous medium, is the ratio heat capacities and is the specific heat.

The goal of this research is to see if a quiescent state can withstand tiny perturbations superimposed on the basic state.Region-1Region-2

Temperature and concentration distributions in the basic state are derived, respectively. where the subscript “b” denotes the basic state and & are the nondimensional basic temperature gradients in region-1 and region-2, respectively.

At the interface and note that and .

To investigate the stability of the basic state, we apply infinite perturbations to region-1 and region-2, respectivelywhere are velocities, are the pressures, are the temperatures, are the salinities and is the magnetic field respectively are perturbed quantities for region-1 and region-2. Substituting (16) and (17) in (1)–(12) and omitting the primes for simplicity, linearizing the equations, and taking curl twice on the momentum equations in both the regions to eliminate pressures. In the dimensionless formulation, scales for length, time, velocity, temperature, species concentration, magnetic field, and temperature gradient, respectively, , , , , , , in the region-1 and , , , , , , are the corresponding quantities in the region-2. The scales for the two regions are determined separately (refer F. Chen and C. F. Chen [43], Chen [44] and Nield [45]), so that each region is of unit depth such that and .

The following perturbed equations areRegion-1Region-2

Solutions for dependent variables are achieved in regions 1 and 2 utilizing normal mode analysis of the dimensionless equations for perturbed variables such as

With and such that , here , are horizontal wave numbers and & are wave numbers in x and y direction, respectively, and are frequencies in region-1 and region-2, respectively. Introducing (26) and (27) in (18)–(25), the following equations are obtained.Region-1: Region-2:

In the above-given expressions, , are the Prandtl numbers, , are the Chandrasekhar numbers and , are the diffusivity ratios, and is the square root of Darcy number.

The linear stability analysis technique and the concept of exchange of stability are applicable for the current scenario, and the study is confined to stationary convection and take and (see Shivakumara et al. [42], Nield [45] and Sumithra [46]). The eigenvalue problem is transformed into the following equations:Region-1Region-2

In the above equations , , , , are, namely, the modified internal Rayleigh number, the internal Rayleigh number, the vertical velocity, the temperature distribution, and concentration distribution, respectively, for region-1 and , , , , are the similar quantities in the region-2. Because the horizontal wave numbers for the combined layers must be the same so that we have here is the depth ratio.

3. Boundary Conditions

The following boundary conditions are nondimensional and are subject to the normal mode analysis:

The velocity conditions are

The conditions for adiabatic-adiabatic boundary conditions are

The salinity conditions are

In the above-given expressions, is the solute diffusivity ratio, is the thermal ratio, is the viscosity ratio, is the thermal Marangoni number (tMn), is the solute Marangoni number , and is the surface tension.

4. Methods of Solution

The eigenvalue problem with an eigenvalue of is formed by the boundary conditions of (36) and the equations (30)–(35). We address this problem using the exact technique procedure, which produces acceptable results when dealing with difficulties of this nature (see Manjunatha et al. [29]).

4.1. Velocity Profile

Using the velocity conditions (37), the solutions of and are obtained by solving (30)–(33) and appropriately written as follows:where

(Refer Appendices).

4.2. Salinity Profile

From (32)–(35), the & we get using conditions (39), as follows:

5. Thermal Marangoni Number

The Eigen value problem is solved for the thermal Marangoni number using the exact technique. The three temperature profiles are considered and the Marangoni numbers are calculated analytically.

5.1. Linear Temperature Profile Model

Consider the linear temperature profile

Introducing equation (43) into equations (31)–(34), we obtained using (38) as follows:

For the linear model, the thermal Marangoni number from (36), is as follows:where

5.2. Parabolic Temperature Profile Model

Consider the (see Sparrow et al. [47]) parabolic temperature profile of the form

The temperature distributions are determined using the temperature boundary conditions (38) by substituting (47) into (31)–(34), as follows:

For the parabolic model, the thermal Marangoni number from (36), is as follows:where

5.3. Inverted Parabolic Temperature Profile Model

Consider the inverted parabolic temperature profile of the form

The temperature distributions are determined using the temperature boundary conditions (38) by substituting equations (51) into equations (31)–(34) as follows:

For an inverse parabolic temperature profile model, the thermal Marangoni number of (36) has the form:where

6. Analysis and Discussion of Results

The present study aims at solving exactly the problem of Marangoni convection in presence of a heat source and temperature gradients. The exact method provides useful results and also general basic temperature profiles can be readily treated with minimum mathematical computations. The graphs are plotted using MATHEMATICA version 11. The thermal Marangoni numbers (tMn’s), linear profile , parabolic profile , and inverted parabolic profile . These t & are drawn versus the depth ratio . The Chandrasekhar number, modified internal Rayleigh numbers for region-1 & region-2, solute Marangoni number, and diffusivity ratio impacts on the tMn are demonstrated in the following figures for various given parameters are , , , , , and .

The implications of the Chandrasekhar number for linear, parabolic, and inverted parabolic temperature profiles are shown in Figure 2. The curves for all three profiles are slightly diverging, showing that bigger depth ratio values, i.e., for the Porous layer Dominant Composite System (PDCS), are more prominent. The value of raises the tMn’s for all three temperature gradients for a fixed depth ratio, hence the Double Component-Magneto-Marangoni-Convection (DCMC) may be advanced by decreasing the values of as well as the possibility of the system becoming unstable. A modified internal Rayleigh number for fluid region influence on the tMn’s are shown in Figure 3. The tMn’s increases as the value increases, indicating that the DCMC can be delayed by increasing the value of from sink to source. Hence, the heat absorption renders the system steady. This parameter’s influence is also stronger for larger depth ratios, implying that the modified internal Rayleigh number is useful for PDCS. The effect of the modified internal Rayleigh number for the porous region on the tMn’s is seen in Figure 4 for all three temperature gradients. The Marangoni number grows as the value increases from sink to source, hence the DCMC can be postponed by increasing the values of the modified internal Rayleigh number. As a result, heat absorption stabilizes the system. This parameter is also effective for a modest value of the depth ratios for the three temperature gradients, as shown in the figure. For all three temperature profiles shown in Figure 5, the influence of sMn on the Marangoni number is comparable. The tMn for DCMC grows as the value increases, hence the DCMC can be delayed by raising the value of . The diverging curves show that the effect of the sMn is more intense for larger depth ratios, i.e., for PDCS. This demonstrates that there is a second diffusing component, at which point the DCMC influences the stability of the system. The tMn for DCMC is versus for different values of the diffusivity ratio , as shown in Figure 6. According to the diverging curves, the diffusivity ratio is only important for larger values of depth ratio, i.e., for PDCS. Also, for all three temperature profiles, the consequence is a drop in tMn’s. As a result, the diffusivity ratio has the effect of preplanning the DCMC and therefore destabilizing the composite layer system.

7. Conclusion

The influence of various physical parameters on the onset of double-component magneto convection are shown graphically in section 6 based on which the following conclusions may be drawn. The core findings from this study can be given as follows:(1)In DCMC of the composite layer, the inverse parabolic temperature gradient is the most stable of the three temperature gradients.(2)For combined systems where large values of the depth ratio, i.e., the porous layer predominates, all physical parameters apply.(3)One can postpone the DCMC in a combined layer by increasing the values of Chandrasekhar number, modified internal Rayleigh numbers, and solute Marangoni number for the set of physical parameters chosen in the study, and one can prepone the DCMC in a composite layer by increasing the values of diffusivity ratio.(4)By managing these parameters, stability-demanding conditions such as crystal growth, and material processing can be effectively achieved.(5)In the process of manufacturing of permanent magnetic material, the destabilizing effects under normal gravity conditions are more efficient. Thus, by increasing the diffusivity ratio of the porous medium, a high-quality permanent magnetic material can be made.

Appendix

Nomenclature

, :Horizontal wave numbers
:Concentrations
:Reference concentration
:Thickness in fluid & porous layer
:Depth ratio
:Magnetic field
:Applied magnetic field
:Permeability of the porous medium
:Ratio of heat capacities
:Solute Marangoni number
:Thermal Marangoni number
:Frequencies
:Total pressure in region-1&2
:Prandtl number
:Chandrasekhar number
:Modified internal Rayleigh number for fluid & porous region
:Internal Rayleigh number for fluid & porous region
:Salinity profile for region-1 and region-2
:Solutal ratio
:Thermal ratio
:Reference temperature
:Velocity vectors
:Velocity profile
:Constant heat sources
:Laplacian operator
:Temperature profiles
:Fluid density
:Porous parameter
:Porosity
:Thermal diffusivities
:Solute diffusivities
:Magnetic viscosities
:Viscosity ratio
:Surface tension
:Fluid viscosity
:Effective viscosity of the fluid in the porous layer
:Diffusivity ratio
:Magnetic permeability
:Temperature distributions in fluid & porous layer
:Subscript denotes fluid and porous layer.

Data Availability

The raw data supporting the conclusions of this article will be made available by the corresponding author without undue reservation.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The article was supported by the Princess Nourah Bint Abdulrahman University Researchers Supporting Project number (Grant no. PNURSP2023R59), Princess Nourah Bint Abdulrahman University, Riyadh, Saudi Arabia.