Next Article in Journal
Reliability Optimization of Hybrid Systems Driven by Constraint Importance Measure Considering Different Cost Functions
Next Article in Special Issue
Fully Nonlinear Evolution of Free-Surface Waves with Constant Vorticity under Horizontal Electric Fields
Previous Article in Journal
Probability Spaces Identifying Ordinal and Cardinal Utilities in Problems of an Economic Nature: New Issues and Perspectives
Previous Article in Special Issue
Start-Up Multilayer Electro-Osmotic Flow of Maxwell Fluids through an Annular Microchannel under Hydrodynamic Slip Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of LTNE on Two-Component Convective Instability in a Composite System with Thermal Gradient and Heat Source

by
Varalakshmi K. Balaji
1,
Manjunatha Narayanappa
1,*,
Ramalingam Udhayakumar
2,
Ghada AlNemer
3,*,
Sumithra Ramakrishna
4 and
Gangadharaih Yeliyur Honnappa
5
1
Department of Mathematics, School of Applied Sciences, REVA University, Bengaluru 560064, India
2
Department of Mathematics, School of Advanced Sciences, Vellore Institute of Technology, Vellore 632014, India
3
Department of Mathematical Sciences, College of Science, Princess Nourah bint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia
4
Department of Mathematics, Nrupathunga University, Bengaluru 560064, India
5
Department of Mathematics, RV Institute of Technology & Management, Bengaluru 560076, India
*
Authors to whom correspondence should be addressed.
Mathematics 2023, 11(20), 4282; https://doi.org/10.3390/math11204282
Submission received: 11 September 2023 / Revised: 10 October 2023 / Accepted: 10 October 2023 / Published: 13 October 2023
(This article belongs to the Special Issue Numerical and Analytical Study of Fluid Dynamics)

Abstract

:
An analytical study is conducted to examine the influence of thermal gradients and heat sources on the onset of two-component Rayleigh–Bènard (TCRB) convection using the Darcy model. The study takes into account the effects of local thermal non-equilibrium (LTNE), thermal profiles, and heat sources. The composite structure is horizontally constrained by adiabatic stiff boundaries, and the resulting solution to the problem is obtained using the perturbation approach. The various physical parameters have been thoroughly examined, revealing that the fluid layer exhibits dominance in the two-layer configuration. It has been observed that the parabolic profile demonstrates greater stability in comparison to the step function. Conversely, in the setup where the porous layer dominates, the step function plays a crucial role in maintaining stability. The porous layer, model (iv), exhibits greater stability in the predominant combined structure, while the linear configuration is characterized by higher instability.

1. Introduction

Double-component convection is very important for studying the evolution of systems with various causes of density fluctuations. Such systems occur in the sun, the oceans, and in the mantle of the Earth. Modeling of two-component convection in a composite layer can be used in many situations, including the storage of nuclear waste, thermal insulation, grain storage, chemical dispersion through water-saturated soil, and soil pollution. However, one more crucial use is for modeling the boundary conditions at interfaces. This matter has garnered significant theoretical interest due to its importance in understanding interface boundary conditions, as well as its aforementioned practical implications. The role of LTNE is significant in inflows with pronounced temperature differences between the solid phase and fluid phase. Wang et al. [1] studied solar air receiver and radiation effects using the LTNE model. Internal heat, using a local thermal non-equilibrium model, was analyzed analytically by Altawallbeh et al. [2,3] in a Maxwell fluid-porous layer with and without a magnetic field. Abidin et al. [4] conducted a study to examine the effects of viscosity in a binary fluid layer. In a two-component system, the closed form of the solution to Marangoni convection with a magnetic field and salinity gradients was studied by Komala and Sumithra [5]. Kuznetsov and Nield [6] studied the LTNE effects for a nanofluid in the presence of a porous domain using the Darcy model. Thirupathi Thumma and Mishra [7] studied nanofluid flow over a stretching sheet in the presence of a heat source/sink using the domain decomposition method. Heat transmission in a porous medium under LTNE conditions was studied by Shashi Prabha et al. [8]. Double-component convective flow over a chemically reactive plate with a heat source/sink was investigated numerically by Kannan and Pullepu [9]. Using the LTNE model, Astanina et al. [10] studied the heat source effect for porous cavities. Hema et al. [11] used the LTNE model to investigate how heat flow affects two-component convection. Pulkit et al. [12] analyzed the phenomenon of double-diffusive convection in a rotating couple-stress ferromagnetic fluid-porous medium, taking into consideration the influence of both gravitational and magnetic fields. Shukla and Gupta [13] studied three-component convection in nanofluids using the LTNE model. Two-component convection in a combined setup with a heated and salted subsurface was investigated by Mahajan and Tripathi [14]. Meften and Ali [15] studied two-component convection with variable viscosity. Capone et al. [16,17] investigated the convection in a rotating, viscous, anisotropic porous medium in the LTNE regime. Using the Brinkmann–Forchheimer model, Meften et al. [18] examined the LTNE effects for dual-component convection with rotation. Safia et al. [19] investigated the impact of magnetic field, temperature, and concentration on convection in a nanofluid. Double diffusion of a nanofluid was investigated by Tahar et al. [20], who focused on magneto-natural convection. A stability assessment of Brinkman–Bènard convection was studied by Siddabasappa et al. [21]. Thermosolutal convection in an LTNE porous medium was examined by Noon and Haddad [22]. Corcione and Quintino [23] examined the dual-component effects for water-based nanofluids using Rayleigh–Bènard convection. Gangadharaih et al. [24] studied the gravity fluctuation and throughflow for an anisotropic porous layer. Ahmed [25] conducted a numerical study on convection in a fluid layer within a square-packed bed enclosure using the LTNE model. Atul and Anand Kumar [26] studied the effects of the solute boundary conditions and the heat source on two-component convection. Yellamma et al. [27] investigated the third-component effect for a combined system with a magnetic field, heat source, and thermal profiles. They obtained the thermal Marangoni number analytically.
Inspired by the aforementioned literature review, the current study examines the effects of LTNE in a composite layer with six thermal gradients on the onset of SCRB convection in a two-layer configuration with an incompressible horizontal fluid flow. The problem has been solved by using the perturbation technique and it is noted that the parabolic profile demonstrates greater stability in comparison to the step function. Conversely, in the setup where the porous layer dominates, the step function plays a crucial role in maintaining stability. This work undoubtedly holds tremendous potential for a multitude of applications in the fields of pure crystal growth formation, as elucidated by Rudolph et al. [28].

2. Materials and Methods

Assume a continuous horizontal layer of incompressible fluid with a thickness d f L and a densely packed region-II d p L of the same fluid that lies behind region-I and is heated continuously by heat sources Q f L and Q p L The area above region-1 is expected to be affected by surface tension effects and concentration as functions of temperature and concentration, while the area below region-II is assumed to be stable. The coordinate system used is Cartesian, with the z-axis pointing upward and the origin located at the interface between the fluid and porous layers, as shown in Figure 1. Both the solid and liquid phases are assumed to be in the LTNE, and it is thought that a solid–fluid field model can characterize the temperature differences between the solid and liquid phases separately for the porous layer.
For the composite system Darcy model, the basic equations assume the Boussinesq hypothesis (see Sumithra and Shyamala [29]).
Fluid layer (region-I):
f L   ·   q f L = 0
q f L t f L + ( q f L   · f L ) q f L = f L P f L ρ 0 + μ f L ρ 0 f L 2 q f L                                                                                           [ 1 α f L ( T f L T f L 0 ) + α f L s ( C f L C f L 0 ) ] g   k ^
T f L t + ( q f L · f L ) T f L = κ f L f L 2 T f L + Q f L
C f L t + ( q f L   ·   f L ) C f L = κ f L s f L 2 C f L
The porous layer (region-II) (see Nield and Bejan [30]):
p L ·   q p L = 0
ϕ 1 q p L t p L + ϕ 2 ( q p L · p L ) q p L = p L P p L ρ 0 + μ p L K   ρ 0 q p L                                 [ 1 α f p L ( T f p L T p L 0 ) α s p L ( T s p L T p L 0 ) + α p L ( C p L C p L 0 ) ] g   k ^
( ρ c p ) f p L [ ϕ T f p L t p L + ( q p L   ·   p L ) T f p L ] = ϕ κ f p L p L 2 T f p L + h ( T s p L T f p L ) + Q p L
( 1 ϕ ) ( ρ c p ) s p L T s p L t p L = ( 1 ϕ ) κ s p L p L 2 T s p L h ( T s p L T f p L )
ϕ C p L t p L + ( q p L   .   p L ) C p L = κ p L s   p L 2 C p L
For the two regions, the basic state is expressed as follows:
[ q f L , P f L , T f L , T s f L , C f L ]   = [ 0 , P f L b ( z f L ) , T f L b ( z f L ) , T s f L b ( z f L ) , C f L b ( z f L ) ] , d f L Δ T   d T f L b d z f L = F f L ( z f L )
[ q p L , P p L , T f p L , T s p L , C p L ] = [ 0 , P p L b ( z p L ) , T f p L b ( z p L ) , T s p L b ( z p L ) , C p L b ( z p L ) ] , d p L Δ T   d T p L b d z p L = F p L ( z p L )
The basic state salinity and temperature distributions are, respectively, found to be:
T f L b ( z f L ) = Q f L z f L ( z f L d f L ) 2 κ f L + ( T U T 0 d f L ) F f L ( z f L ) + T 0       0 z f L d f L
T f p L b ( z p L ) = T s p L b ( z p L ) = Q p L z p L ( z p L + d p L ) 2 ϕ κ p L + ( T L T 0 d p L ) F p L ( z p L ) + T 0       0 z p L d p L
C f L b ( z f L ) = ( C U C 0 d f L ) z f L + C 0       0 z f L d f L
C p L b ( z p L ) = ( C L C 0 d p L ) z p L + C 0       0 z p L d p L
Here ,   T 0 = κ f L d p L T u + κ p L d f L T L κ f L d p L + κ p L d f L + d f L d p L ( Q f L d f L + Q p L d p L ) 2 ( κ f L d p L + κ p L d f L )
C 0 = κ f L s C U d p L + κ p L s C L d f L κ p L s d f L + κ f L s d p L .
where the b symbolizes the basic state, and F f L ( z f L ) and F p L ( z p L ) are the non-dimensionalized thermal gradient that fulfils the situation 0 1 F f L ( z f L )     d z f L = 0 1 F p L ( z p L )     d z p L = 1 .
To study the basic solution’s stability, the following mild disturbances are introduced.
[ q f L , P f L , T f L , C f L ] = [ 0 , P f L b ( z f L ) , T f L b ( z f L ) , C f L b ( z f L ) ] + [ q f L , P f L , θ , S ]
[ q p L , P p L , T f p L , T s p L , C p L ] = [ 0 , P p L b ( z p L ) , T p L b ( z p L ) , T s p L b ( z p L ) , C p L s b ( z p L ) ]                                                                                                                                                                                                 + [ q p L , P p L , θ p L , θ s p L , S p L ]
The prime symbol represents a perturbation from the corresponding values at equilibrium. The physical values in Equations (1)–(9) are now augmented with the introduction of (12) and (13), and subsequently linearized based on convection. The pressure term is subsequently eliminated from Equations (2) and (6) through the application of the curl operation twice, while retaining only the upright factor.
The variables are non-dimensionalized and normal mode study is applied (see Yellamma et al. [31]); the corresponding distinctive quantities are found to be:
In 0 z f L 1 ,
( D f L 2 a f L 2 ) 2 W f L = R f L a f L 2 θ f L R f L s a f L 2 S f L
( D f L 2 a f L 2 ) 2 θ f L + [ F f L ( z f L ) + R I ( 2 z f L 1 ) ] W f L = 0
τ f L s ( D f L 2 a f L 2 ) 2 S f L + W f L = 0
In 0 z p L 1 ,
( D p L 2 a p L 2 ) 2 W p L = R f p L β 2 a p L 2 θ f p L τ s p L R s p L β 2 a p L 2 θ s p L + τ p L s R p L s β 2 a p L 2 S p L
ϕ ( D p L 2 a p L 2 ) θ f p L + [ F p L ( z p L ) + R I M ( 2 z p L + 1 ) ] w p L = H ( θ s p L θ f p L )
( 1 ϕ ) ( D p L 2 a p L 2 ) θ s p L = H ( θ s p L θ f p L )
τ p L s ( D p L 2 a p L 2 ) 2 S p L + W p L = 0
where R f L = g α f L ( T 0 T U ) d f L 3 ν f L κ f L ,   R p L = g α p L ( T L T 0 ) d p L 3 ν p L κ p L ,   R f L s = g α f L s ( C 0 C U ) d f L 3 ν f L κ f L s ,   R p L s = g α p L s ( C L C 0 ) d p L 3 ν p L κ p L s ,   R s p L = g α s p L ( T L T 0 ) d p L 3 ν p L κ s p L are, respectively, the thermal, solute, and solid phase Rayleigh numbers, R I = Q f L d f L 2 κ f L , R I M = Q p L d p L 2 κ f p L ,   R I = R I 2 ( T 0 T U ) ,   R I M = R I M 2 ( T L T 0 ) ,   β 2 = K d p L 2 = D a , H = h d p L 3 κ f p L are the internal and corrected Rayleigh numbers, the Darcy number, and the inter-phase heat transfer coefficient, respectively; α f p L ,   α s p L are the fluid and solid phase thermal expansion coefficients; κ f p L ,   κ s p L are the fluid phase and solid phase thermal diffusivities; and κ f L s , κ p L s are the fluid and porous layer solute diffusivities (see nomenclature).
After non-dimensionalization, normal mode expansion is carried out on the proper boundary conditions.
At z f L = d f L , the situations are
W f L ( 1 ) = 0 , D f L W f L ( 1 ) = 0 , D f L θ f L ( 1 ) = 0 , D f L S f L ( 1 ) = 0
At z p L = d p L , the situations are
W p L ( 0 ) = 0 ,   D p L W p L ( 0 ) = 0 ,   D p L θ f p L ( 0 ) = 0 ,   D p L θ s p L ( 0 ) = 0 ,   D p L S p L ( 0 ) = 0
At z f L = 0 , z p L = 1 , the situations are
T ^ W f L ( 0 ) = d ^ 2 W p L ( 1 ) ,   T ^ D f L W f L ( 0 ) = d ^ μ ^ D p L W p L ( 1 ) ,
θ f L ( 0 ) = T ^ d ^ 2 θ f p L ( 1 ) ,   θ f L ( 0 ) = T ^ d ^ 2 θ s p L ( 1 ) ,
D f L θ f L ( 0 ) = d ^ 2 D p L   θ f p L ( 1 ) ,   D f L θ f L ( 0 ) = d ^ 2 D p L θ s p L ( 1 ) ,
S f L ( 0 ) = S ^   S p L ( 1 ) ,   D f L S f L ( 0 ) = D p L S p L ( 1 ) ,
T ^ d ^ β 2 D f L 3 W f L ( 0 ) = D p L W p L ( 1 ) + μ ^   β 2 D p L 3 W p L ( 1 ) ,
T ^ D f L 2 W f L ( 0 ) = μ ^ D p L 2 W p L ( 1 ) .
where T ^ = ( T L T 0 ) ( T 0 T U ) ,   S ^ = C L C 0 C 0 C U ,   d ^ = d p L d f L , and μ ^ = μ p L μ f L are the thermal, salinity, depth, and viscosity ratio, respectively.

3. Method of Solution

The task at hand involves addressing the eigenvalue problem within the context of arbitrary boundary conditions related to temperature and concentration. This will be accomplished by employing the standard perturbation method, with the wave number a f L 2 , a p L 2 serving as the perturbation parameter. The results of such a study can be used to assess the reliability of the method, and they can also be used to justify the method’s application to the solution of problems involving convective instability in which the critical stability parameter must be determined analytically.
The powers a f L 2 , a p L 2 have been added to the dependent variables in both levels, as seen below:
[ W f L θ f L S f L ] = j = 0 a f L 2 j [ W f L j θ f L j S f L j ]   and   [ W p L θ f p L θ s p L S p L ] = j = 0 a p L 2 j [ W p L j θ f p L j θ s p L j S p L j ]
When solving a f L 2 , a p L 2 using the preceding equation in (14)–(20), we obtain zero-order expressions (see Appendix A).
The first-order equations for a f L 2 , a p L 2 are as follows:
For region-I,
D f L 4 W f L 1 R f L T ^ + R f L s S ^ = 0
D f L 2 θ f L 1 T ^ + [ F f L ( z f L ) + R I ( 2 z f L 1 ) ] W f L 1 = 0
τ f L s D f L 2 S f L 1 τ f L s S ^ + W f L 1 = 0
For region-II,
D p L 2 W p L 1 + R f p L β 2 + τ s p L R s p L β 2 τ p L s R p L s β 2 = 0
( ϕ D p L 2 H ) θ f p L 1 + H θ s p L 1 + [ F p L ( z p L ) + R I M ( 2 z p L + 1 ) ] W p L 1 ϕ = 0
[ ( 1 ϕ ) D p L 2 H ] θ s p L 1 + H θ p L 1 ( 1 ϕ ) = 0
τ p L s D p L 2 S p L τ p L s + W p L = 0
The related boundary conditions are as follows:
W f L 1 ( 1 ) = 0 , D f L 1 W f L 1 ( 1 ) = 0 , D f L 1 θ f L 1 ( 1 ) = 0 , D f L 1 S f L 1 ( 1 ) = 0 ,
T ^ W f L 1 ( 0 ) = d ^ 2 W p L 1 ( 1 ) ,   T ^ D f L W f L 1 ( 0 ) = d ^ μ ^ D p L W p L 1 ( 1 ) ,
θ f L 1 ( 0 ) = T ^ d ^ 2 θ p L 1 ( 1 ) ,   θ f L 1 ( 0 ) = T ^ d ^ 2 θ s p L 1 ( 1 ) , S f L 1 ( 0 ) = S ^ d ^ 2 S p L 1 ( 1 ) ,
D f L θ f L 1 ( 0 ) = d ^ 2 D p L θ f p L 1 ( 1 ) , D f L θ f L 1 ( 0 ) = d ^ 2 D p L θ s p L 1 ( 1 ) , D f L S f L 1 ( 0 ) = d ^ 2 D p L S p L 1 ( 1 ) ,
T ^ d ^ β 2 D f L 3 W f L 1 ( 0 ) = D p L W p L 1 ( 1 ) + μ ^ β 2 D p L 3 W p L 1 ( 1 ) ,   T ^ D f L 2 W f L 1 ( 0 ) = μ ^ D p L 2 W p L 1 ( 1 ) ,
W p L 1 ( 0 ) = 0 ,   D p L W p L 1 ( 0 ) = 0 ,   D p L θ f p L 1 ( 0 ) = 0 ,   D p L θ s p L 1 ( 0 ) = 0 , D p L S p L 1 ( 0 ) = 0 .
The solutions of W f L 1 and W p L 1 are as follows:
W f L 1   ( z f L ) = δ f L 1 + δ f L 2   z f L + δ f L 3   z f L 2 + δ f L 4     z f L 3 + ( R f L T ^ R f L s S ^ ) 24 z f L 4
W p L 1 ( z p L ) = δ p L 5 + δ p L 6   z p L β 2 ( R f p L + R s p L   τ s p L R p L s   τ p L s ) 2 z p L 2
where δ f L 1 , δ f L 2 , δ f L 3 , δ f L 4 , δ f L 5 , δ f L 6 are the constants obtained utilizing velocity boundary conditions; for the form they are in, see Appendix A.

4. Condition of Solvability

The solvability condition is derived by solving the differential equation with the necessary boundary conditions concerning heat.
d ^ 2 ( 1 + R I M + τ f L s ) 0 1 W p L 1 ( z p L )   · F p L ( z p L )   d z p L + 2 d ^ 2 R I M 0 1 W p L 1 ( z p L ) · F p L ( z p L )   d z p L + ( 1 R I + τ p L s ) 0 1 W f L 1 ( z f L )   ·   F f L ( z f L )   d z f L + 2 R I 0 1 W f L 1 ( z f L ) · F f L ( z f L )   d z f L                                                                                                                                                             = ( d ^ 2 + T ^ ) + ( d ^ 2 + S ^ ) τ f L s τ p L s
Computing R C for different profiles can be achieved by substituting expressions W f L 1 ( z f L ) and W p L 1 ( z p L ) from Equations (29) and (30) into Equation (31).

5. Thermal Gradients

The following thermal gradients are considered for the present study (refer to Table 1; see Rudraiah et al. [32], Vasseur and Robillard [33], and Shivakumara [34]).

5.1. Model (i): Linear Thermal Gradient

The R C is attained for the linear model F f L ( z f L ) = 1 and F p L ( z p L ) = 1 by utilizing (31):
R C = ( d ^ 2 + T ^ ) + ( d ^ 2 + S ^ ) τ f L s τ p L s d ^ 4 κ p L s 2 τ p L s Ψ 1 R p L s { Φ 11 + Φ 12 } T ^ { Φ 13 + Φ 14 } + [ α f p L   κ f p L 2 + α s p L κ s p L 2 τ s p L ] d ^ 4   Ψ 1
where Ψ 1 = Φ 1 + Φ 2 + Φ 3 + Φ 4 + Φ 5 + Φ 6 + Φ 7 + Φ 8 + Φ 9 + Φ 10 ,   Φ 1 = d ^ 2 β 2 ( 1 + R I M + τ f L s ) 6 ,   Φ 2 = d ^ 2 β 2 R I M 4 ,   Φ 3 = d ^ 2 β 2 ( 1 R I + τ p L s ) 2 ,   Φ 4 = d ^ β 2 μ ^ ( 1 R I + τ p L s ) 2 T ^ ,   Φ 5 = β 2 μ ^ ( 1 R I + τ p L s ) 6 T ^ ,   Φ 6 = ( 1 R I + τ p L s ) 24 d ^ T ^ ,   Φ 7 = d ^ 2 β 2 R I 2 T ^ ,   Φ 8 = 2 d ^ β 2 μ ^ R I 3 T ^ ,   Φ 9 = β 2 μ ^ R I 4 T ^ ,   Φ 10 = R I 15 d ^ T ^ ,   Φ 11 = S ^ ( 1 R I + τ p L s ) 120 ,   Φ 12 = 2 S ^ R I 144 ,   Φ 13 = ( 1 R I + τ p L s ) 120 ,   Φ 14 = 2 R I 144 .

5.2. Model (ii): Parabolic Thermal Gradient

The R C is attained for the parabolic model F f L ( z f L ) = 2 z f L ,   F p L ( z p L ) = 2 z p L by utilizing (31):
R C = ( d ^ 2 + T ^ ) + ( d ^ 2 + S ^ ) τ f L s τ p L s d ^ 4 κ p L s 2 τ p L s Ψ 2 R p L s { Ν 11 + Ν 12 } T ^ { Ν 13 + Ν 14 } + [ α f p L   κ f p L 2 + α s p L κ s p L 2 τ s p L ] d ^ 4   Ψ 2
where Ψ 2 = Ν 1 + Ν 2 + Ν 3 + Ν 4 + Ν 5 + Ν 6 + Ν 7 + Ν 8 + Ν 9 + Ν 10 ,   Ν 1 = d ^ 2 β 2 ( 1 + R I M + τ f L s ) 4 , Ν 2 = 2 d ^ 2 β 2 R I M 5 , Ν 3 = d ^ 2 β 2 ( 1 R I + τ p L s ) 2 T ^ , Ν 4 = 2 d ^ μ ^ β 2 ( 1 R I + τ p L s ) 3 T ^ , Ν 5 = β 2 μ ^ ( 1 R I + τ p L s ) 4 T ^ , Ν 6 = ( 1 R I + τ p L s ) 15 d ^ T ^ , Ν 7 = 3 d ^ 2 β 2 R I 3 T ^ , Ν 8 = d ^ β 2 μ ^ R I T ^ ,   Ν 9 = 2 β 2 μ ^ R I 5 T ^ , Ν 10 = R I 9 d ^ T ^ , Ν 11 = S ^ ( 1 R I + τ p L s ) 72 , Ν 12 = S ^ R I 42 , Ν 13 = ( 1 R I + τ p L s ) 72 , Ν 14 = R I 42 .

5.3. Model (iii): Inverted Parabolic Thermal Gradient

The R C is attained for this model F f L ( z f L ) = 2 ( 1 z f L ) and F p L ( z p L ) = 2 ( 1 z p L ) by utilizing (31):
R C = ( d ^ 2 + T ^ ) + ( d ^ 2 + S ^ ) τ f L s τ p L s d ^ 4 κ p L s 2 τ p L s Ψ 3 R p L s { A 11 + A 12 } T ^ { A 13 + A 14 } + [ α f p L   κ f p L 2 + α s p L κ s p L 2 τ s p L ] d ^ 4     Ψ 3
where Ψ 3 = A 1 + A 2 + A 3 + A 4 + A 5 + A 6 + A 7 + A 8 + A 9 + A 10 ,   A 1 = d ^ 2 β 2 ( 1 + R I M + τ f L s ) 12 , A 2 = d ^ 2 β 2 R I M 10 , A 3 = d ^ 2 β 2 ( 1 R I + τ p L s ) 2 T ^ , A 4 = d ^ β 2 μ ^ ( 1 R I + τ p L s ) 3 T ^ A 5 = β 2 μ ^ ( 1 R I + τ p L s ) 12 T ^ , A 6 = ( 1 R I + τ p L s ) 60 d ^ T ^ , A 7 = d ^ 2 β 2 R I 3 T ^ , A 8 = d ^ β 2 μ ^ R I 3 T ^ , A 9 = β 2 μ ^ R I 10 T ^ , A 10 = R I 45 d ^ T ^ , A 11 = S ^ ( 1 R I + τ p L s ) 360 , A 12 = S ^ R I 252 , A 13 = ( 1 R I + τ p L s ) 360 , A 14 = R I 252 .

5.4. Model (iv): Piecewise Linear Gradient Heated from Below

F f L ( z f L ) = { ε f L 1 , 0 z f L ε f L 0 , ε f L z f L 1   and   F p L ( z p L ) = { ε p L 1 , 0 z p L ε p L 0 , ε p L z p L 1
The Rc is attained for this model by utilizing (35) in (31) and is calculated as follows:
R C = ( d ^ 2 + T ^ ) + ( d ^ 2 + S ^ ) τ f L s τ p L s d ^ 4 κ p L s 2 τ p L s Ψ 4 R p L s { B 11 + B 12 } T ^ { B 13 + B 14 } + [ α f p L   κ f p L 2 + α s p L κ s p L 2 τ s p L ] d ^ 4     Ψ 4
where Ψ 4 = B 1 + B 2 + B 3 + B 4 + B 5 + B 6 + B 7 + B 8 + B 9 + B 10 ,   B 1 = d ^ 2 β 2 ( 1 + R I M + τ f L s ) ε p L 2 6 , B 2 = d ^ 2 β 2 R I M ε p L 3 4 , B 3 = d ^ 2 β 2 ( 1 R I + τ p L s ) 2 T ^ , B 4 = d ^ β 2 μ f L ( 1 + R I + τ p L s ) ε f L 2 T ^ B 5 = β 2 μ ^ ( 1 R I + τ p L s ) ε f L 2 6 T ^ , B 6 = ( 1 R I + τ p L s ) ε f L 3 24 d ^ T ^ , B 7 = d ^ 2 β 2 R I ε f L 2 T ^ , B 8 = 2 d ^ β 2 μ ^ R I ε f L 2 3 T ^ , B 9 = β 2 μ ^ R I ε f L 3 4 T ^ , B 10 = R I ε f L 4 15 d ^ T ^ , B 11 = S ^ ( 1 R I + τ p L s ) ε f L 4 120 , B 12 = S ^ R I ε f L 5 72 , B 13 = ( 1 R I + τ p L s ) ε f L 4 120 , B 14 = R I ε f L 5 72 .

5.5. Model (v): Piecewise Linear Gradient Cooled from Above

F f L ( z f L ) = { 0 , 0 z f L ( 1 ε f L ) ε f L 1 , ( 1 ε f L ) z f L 1   and   F p L ( z p L ) = { 0 , 0 z p L ( 1 ε p L ) ε p L 1 , ( 1 ε p L ) z p L 1
The Rc is attained for this model by utilizing (37) in (31) and is calculated as follows:
R C = ( d ^ 2 + T ^ ) + ( d ^ 2 + S ^ ) τ f L s τ p L s d ^ 4 κ p L s 2 τ p L s Ψ 5 R p L s { 11 + 12 } T ^ { 13 + 14 } + [ α f p L   κ f p L 2 + α s p L κ s p L 2 τ s p L ] d ^ 4     Ψ 5
where Ψ 5 = 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 9 + 10 ,   1 = d ^ 2 β 2 ( 1 + R I M + τ f L s ) 2 ε p L ( 1 ( 1 ε f L ) 3 ) 3 ,   2 = d ^ 2 β 2 R I M ε p L ( 1 ( 1 ε p L ) 4 ) 4 ,   3 = d ^ 2 β 2 ( 1 R I + τ p L s ) 2 ,     4 = d ^ β 2 μ ^ ( 1 + R I + τ p L s ) 2 T ^ ε f L ,   R 5 = β 2 μ ^ 1 R I + τ p L s 1 1 ε f L 3 6 T ^ ε f L ,   R 6 = 1 R I + τ p L s 1 1 ε f L 4 24 d ^ T ^ ε f L ,   R 7 = d ^ 2 β 2 R I 1 1 ε f L 2 2 ε f L T ^ ,   R 8 = 2 d ^ β 2 μ ^ R I 1 1 ε f L 3 3 T ^ ε f L ,   9 = β 2 μ ^ R I ( 1 ( 1 ε f L ) 4 ) 4 T ^ ε f L ,   10 = R I ( 1 ( 1 ε f L ) 5 ) 15 d ^ T ^ ε f L ,   11 = S ^ ( 1 R I + τ p L s ) ( 1 ( 1 ε f L ) 5 ) 120 ε f L ,   12 = S ^ R I ( 1 ( 1 ε f L ) 6 ) 72 ε f L ,   13 = ( 1 R I + τ p L s ) ( 1 ( 1 ε f L ) 5 ) 120 ε f L ,   14 = R I ( 1 ( 1 ε f L ) 6 ) 72 ε f L .

5.6. Model (vi): Step Function

The basic temperature in this profile drops rapidly by a certain amount Δ T f L at z f L = ε f L and Δ T p L at z p L = ε p L , otherwise it is uniform. Accordingly,
F f L ( z f L ) = δ ( z f L ε f L )   and   F p L ( z p L ) = δ ( z p L ε p L )
The R C is attained for this model by utilizing (39) in (31) and is calculated as follows:
R C = ( d ^ 2 + T ^ ) + ( d ^ 2 + S ^ ) τ f L s τ p L s d ^ 4 κ p L s 2 τ p L s Ψ 6 R p L s { 25 + 26 } T ^ { 27 + 28 } + [ α f p L   κ f p L 2 + α s p L κ s p L 2 τ s p L ] d ^ 4   Ψ 6
where Ψ 6 = 15 + 16 + 17 + 18 + 19 + 20 + 21 + 22 + 23 + 24 ,   15 = d ^ 2 β 2 ( 1 + R I M + τ f L s ) ε p L 2 2 ,   16 = d ^ 2 β 2 R I M ε p L 3 4 ,   17 = d ^ 2 β 2 ( 1 R I + τ p L s ) 2 T ^ ,   18 = d ^ β 2 μ ^ ( 1 + R I + τ p L s ) ε f L T ^ ,   19 = β 2 μ ^ ( 1 R I + τ p L s ) ε f L 2 2 T ^ ,   20 = ( 1 R I + τ p L s ) ε f L 3 6 d ^ T ^ ,   21 = d ^ 2 β 2 R I ε f L T ^ ,   22 = 2 d ^ β 2 μ ^ R I ε f L 2 T ^ ,   23 = β 2 μ ^ R I ε f L 3 T ^ ,   24 = R I ε f L 4 6 d ^ T ^ ,   25 = S ^ ( 1 R I + τ p L s ) ε f L 4 24 ,   26 = S ^ R I ε f L 5 12 ,   27 = ( 1 R I + τ p L s ) ε f L 4 24 ,   28 = R I ε f L 5 12 .

6. Results and Discussion

The theoretical investigation focuses on the initiation of TCRB convection in a combined structure. This system comprises a region-I positioned above a region-II that is saturated with the same fluid. The analysis is conducted in accordance with the LTNE model. The analytical calculations and graphical representations are performed utilizing the advanced computational software known as MATHEMATICA, version 11. The eigenvalue issue is resolved by the utilization of the perturbation method, resulting in the derivation of an analytical expression for the critical Rayleigh number. This expression is obtained for six distinct thermal gradients, as indicated in Table 1. The graph depicts the relationship between the critical Rayleigh number ( R C ) and the depth ratio ( d ^ ) for all the cases examined in the current study, specifically in the context of adiabatic rigid boundaries. Also, there is a noticeable pattern in the profiles’ stability. The range of values under consideration in the present analysis is as follows: α s p L = 2 ,   κ f p L = 0.2 ,   κ s p L = 0.3 ,   τ f L s = 0.2 ,   τ p L s = 0.3 ,   R f L s = 10 ,   R p L s = 10 ,   T ^ = 1.0 ,   R I = 0.1 ,   R I M = 1.0 ,   β = 1 ,   κ p L s = 0.5 ,   κ f L s = 0.3 and α f p L = 1 .
The impact of the thermal expansion ratio α f p L on R C is illustrated in Figure 2. The curves exhibit divergence when considering the values α f p L = 1 ,   5 , and 10. The figure clearly demonstrates a positive correlation between α f p L and R C , indicating that an increase in α f p L results in a corresponding increase in R C . Consequently, the setup can achieve stability, thereby delaying the onset of TCRB convection. Based on the graph, it can be observed that the linear profile exhibits greater stability compared to the parabolic profile. In comparison, PLHB exhibits greater stability, while SF demonstrates a higher degree of instability. The impact of the parameter is minimal in the parabolic profiles, while it is more pronounced in the linear profiles. According to the data presented in Figure 3, an increase in the variable solid phase expansion ratio α s p L leads to a corresponding increase in the R C . Consequently, this adjustment in the setup has the potential to enhance the stability and delay the onset of TCRB convection. Based on the graph, it can be observed that the linear profile exhibits the highest level of stability, while the parabolic profile demonstrates instability. In a similar vein, it can be observed that PLHB exhibits a higher degree of stability compared to SF, which is relatively less stable. The impact of the parameter is minimal in parabolic profiles, whereas it is more significant in linear profiles. It is observed that nearly equivalent outcomes are achieved for the variables α f p L and α s p L . It is evident from Figure 4 that the growth in the porous parameters β leads to a rise in the R C . Consequently, this allows for the stabilization of the setup, resulting in a delayed onset of TCRB convection. As a result, the system is postponed. The graphs indicate that the linear and PLHB models exhibit greater stability, while the parabolic and SF models demonstrate higher levels of instability. Additionally, the impact of the porous parameter is significantly more pronounced for the linear and PLHB models, moderately significant for the inverted parabolic and PLCA models, and detrimental for the parabolic and SF thermal gradient models. Figure 5 and Figure 6 exhibit the effects of the fluid and solid phase thermal diffusivity ratios, respectively, κ f p L and κ s p L . As the ratios increase, R C increases. As a result, the setup is stable, delaying the onset of TCRB convection. Hence the system is postponed. The effects of the corrected internal Rayleigh numbers R I and R I M are illustrated in Figure 7 and Figure 8. The system experiences destabilization as the R I parameter increases, resulting in the earlier occurrence of TCRB convection for higher values of this parameter. The TCRB convection effectively regulates smaller values of the R I . It is observed in Figure 8 that there is an inverse response from the porous layer R I M . This enhances the system stability by increasing the value of R C . The observed curves exhibit a divergence pattern, suggesting that the sensitivity arises when the porous layer predominates over the composite layer across all thermal gradients. Therefore, the TCRB convection is delayed. The influence of the solute Rayleigh numbers R f L s and R p L s over the R C are depicted in Figure 9 when R f L s = R p L s = 10 , 20 , 30 , and 30 . With the increase in R f L s and R p L s , R C decreases, destabilizing the setup. TCRB convection, hence, starts up quickly. The observed phenomenon could potentially arise from the existence of a secondary constituent.

7. Conclusions

The primary focus of this theoretical investigation is the initiation of the two-component Rayleigh–Bènard (TCRB) convection in a combined structure. The system consists of a region-I situated above a region-II that is fully saturated with the same fluid. The analysis is performed following the LTNE model. The present endeavor necessitates the examination of the situation within the framework of arbitrary boundary conditions for salinity and temperature. The eigenvalue problem is solved using the perturbation method, which leads to the derivation of an analytical expression for the critical Rayleigh number. The aforementioned expression is derived from six unique thermal gradients, as outlined in Table 1. The present study has yielded the following outcomes.
(i)
The onset of TCRB convection is supported by a corrected internal Rayleigh number in region-I and the solute Rayleigh number. By increasing the system’s instability, the six profiles speed up the beginning of TCRB convection with the corrected internal Rayleigh number in region-I and the solute Rayleigh number.
(ii)
When applied to each of the six profiles, the thermal ratio, the thermal diffusivity ratio, the porous parameter, and the corrected internal Rayleigh number all work together to postpone the onset of TCRB convection.
(iii)
In a combined structure arrangement with region-I on top, the parabolic profile is more stable than the step function, but in a setup with region-II on top, the step function takes the lead.
(iv)
The PLHB structure exhibits the highest stability among the combined structures in region-II, whereas the linear structure demonstrates the lowest stability.
(v)
The commencement of TCRB convection in a composite fluid and porous layer system may be efficiently controlled by selecting the proper parameters.

Author Contributions

Conceptualization, V.K.B., M.N., S.R. and R.U.; methodology, V.K.B. and M.N.; validation, G.A., M.N., R.U. and G.Y.H.; formal analysis, G.Y.H. and M.N.; investigation, R.U. and G.A.; writing—original draft preparation, G.Y.H., M.N., S.R., V.K.B. and G.A.; writing—review and editing, G.A.; visualization, all the authors; supervision, M.N. and G.A.; project administration, M.N. and R.U.; funding acquisition, G.A. All authors have read and agreed to the published version of the manuscript.

Funding

Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R45), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Data Availability Statement

Not applicable.

Acknowledgments

Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R45), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

English LettersGreek Letters
a f L , a p L Wave numbers ρ 0 Reference density, kg/m3
C f L , C p L Concentrations, mol/m3 ρ f L , ρ p L Fluid and porous layer density, kg/m3
C 0 Reference concentration, mol/m3 μ f L , μ p L Fluid viscosity and porous viscosity,
kg/ms
C p Specific heat capacity,
J/kg∙K
Δ T Temperature difference, K
D f L , D p L D f L = d d z f L , D p L = d d z p L Δ C Concentrate difference, mol/m3
D a Darcy number κ f L , κ p L Thermal diffusivities in fluid and porous layer, m2/s
d f L , d p L Thickness, m κ f L s , κ p L s Solute diffusivity ratio
d ^ Depth ratio κ p L , κ s p L Thermal diffusivity of fluid and solid phase in porous layer, m2/s
g Gravity, m/s2 α f L , α p L Thermal expansion coefficients in fluid and porous layer, 1/K
H Scaled interface heat transfer coefficient, W/mK α f p L , α s p L Thermal expansion coefficients in fluid and solid phase in porous layer, 1/K
h Inter-phase heat transfer, W/m2/K α f L s , α p L s Solute thermal expansion coefficients in fluid and porous layer, 1/K
K Permeability, H/m τ s p L Inter-phase thermal diffusivity ratio, m2/s
k ^ (0, 0, 1) τ f L , τ p L Thermal diffusivity ratio, m2/s
P f L , P p L Pressure, kg m−1s−2 τ f L s , τ p L s Solute thermal diffusivity ratios, m2/s
q f L , q p L Velocity vectors, m/s θ f L , θ p L Amplitude of perturbed temperature, K
R C Critical Rayleigh number β Porous parameter
R f L s , R s p L Solute Rayleigh numbers ϕ Porosity
R f L , R p L Rayleigh numbers μ ^ Viscosity ratio
R I , R I M Internal Rayleigh numbersSubscripts
R I , R I M Corrected internal Rayleigh numbers b Basic state
S ^ Solute ratio f L Fluid layer
S f L , S p L Amplitude of perturbed concentration p L Porous layer
T f L , T p L , T s p L Temperatures, K f p L Fluid phase porous layer
T 0 Interface temperature, K s p L Solid phase porous layer
T ^ Thermal ratio f L s Fluid layer salinity
W f L , W p L Dimensionless vertical velocities p L s Porous layer salinity
F f L ( z f L ) , F p L ( z p L ) Fluid and porous layer temperature gradient

Appendix A

The zero-order expressions are as follows:
For region-I,
D f L 4 W f L 0 = 0
D f L 2 θ f L 0 + [ F f L ( z f L ) + R I ( 2 z f L 1 ) ] W f L 0 = 0
τ f L s D f L 2 S f L 0 + W f L 0 = 0
For region-II,
D p L 2 W p L 0 = 0
ϕ D p L 2 θ f p L 0 + [ F p L ( z p L ) + R I M ( 2 z p L + 1 ) ] w p L 0 + H ( θ s p L 0 θ f p L 0 ) = 0
( 1 ϕ ) D p L 2 θ s p L 0 + H ( θ s p L 0 θ p L 0 ) = 0
τ p L s D p L 2 S p L 0 + W p L 0 = 0
The zero-order boundary conditions are as follows:
W f L 0 ( 1 ) = 0 , D f L W f L 0 ( 1 ) = 0 , D f L θ f L 0 ( 1 ) = 0 ,       D f L S f L 0 ( 1 ) = 0 , T ^ W f L 0 ( 0 ) = W p L 0 ( 1 ) ,
T ^ d ^ D f L W f L 0 ( 0 ) = μ ^   W p L 0 ( 1 ) , θ f L 0 ( 0 ) = T ^ θ f p L 0 ( 1 ) , θ f L 0 ( 0 ) = T ^ θ s p L 0 ( 1 ) , S f L 0 ( 0 ) = S ^ S p L 0 ( 1 ) ,
D f L S f L 0 ( 0 ) = D p L S p L 0 ( 1 ) , D f L θ f L 0 ( 0 ) = D p L θ p L 0 ( 1 ) , D f L θ f L 0 ( 0 ) = D p L θ s p L 0 ( 1 ) ,
T ^ d ^ 2 D f L 2 W f L 0 ( 0 ) = μ ^ D p L 2 W p L 0 ( 1 ) , T ^ d ^ 3 β 2 D f L 3 W f L 0 ( 0 ) = D p L W f L 0 ( 1 ) + μ ^ β 2 D p L 3 W p L 0 ( 1 ) ,
W p L 0 ( 0 ) = 0 , D p L W p L 0 ( 0 ) = 0 , D p L θ f p L 0 ( 0 ) = 0 , D p L θ s p L 0 ( 0 ) = 0 , D p L S p L 0 ( 0 ) = 0 .
The zeroth-order equations have solutions which are given by
W f L 0 ( z ) = 0 , θ f L 0 ( z ) = T ^ , W p L 0 ( z p L ) = 0 , θ f p L 0 ( z p L ) = 1 , θ s p L 0 ( z p L ) = 1 .
δ f L 1 = d ^ 2 β 2 ( R f p L + R s p L   τ s p L R p L s   τ p L s ) 2 T ^ , δ f L 2 = μ ^ d ^ β 2 ( R f p L + R s p L   τ s p L R p L s   τ p L s ) T ^ ,
δ f L 3 = μ ^ β 2 ( R f p L + R s p L   τ s p L R p L s     τ p L s ) 2 T ^ , δ f L 4 = ( R f p L + R s p L   τ s p L R p L s   τ p L s ) 2 T ^ ,
δ p L 5 = 0 ,     δ p L 6 = 0 .

References

  1. Wang, P.; Vafai, K.; Liu, D.Y. Analysis of radiative effect under local thermal non-equilibrium conditions in porous media-application to a solar air receiver. Numer. Heat Transf. Part A Appl. 2014, 65, 931–948. [Google Scholar] [CrossRef]
  2. Altawallbeh, A.A.; Hashim, I.; Tawalbeh, A.A. Thermal non-equilibrium double diffusive convection in a Maxwell fluid with internal heat source. J. Phys. Conf. Ser. 2018, 1132, 012027. [Google Scholar] [CrossRef]
  3. Altawallbeh, A.A.; Hashim, I.; Bhadauria, B.S. Magneto-double diffusive convection in a viscoelastic fluid saturated porous layer with internal heat source. AIP Conf. Proc. 2019, 2116, 030015. [Google Scholar] [CrossRef]
  4. Abidin, N.H.; Mokhtar, N.F.; Majid, Z.A. Onset of Bènard-Marangoni instabilities in a double diffusive binary fluid layer with temperature dependent viscosity. Numer. Algebra Control Optim. 2019, 9, 413–421. [Google Scholar] [CrossRef]
  5. Komala, B.; Sumithra, R. Effects of non-uniform salinity gradients on the onset of double diffusive magneto—Marangoni convection in a composite layer. Int. J. Adv. Sci. Technol. 2019, 28, 874–885. [Google Scholar]
  6. Kuznetsov, A.V.; Nield, D.A. Effect of local thermal non-equilibrium on the onset of convection in a porous medium layer saturated by a nanofluid. Transp. Porous Med. 2010, 83, 425–436. [Google Scholar] [CrossRef]
  7. Thumma, T.; Mishra, S.R. Effect of nonuniform heat source/sink, and viscous and Joule dissipation on 3D Eyring–Powell nanofluid flow over a stretching sheet. J. Comput. Des. Eng. 2020, 7, 412–426. [Google Scholar] [CrossRef]
  8. Prabha, G.S.; Bharathi, M.C.; Kudenatti, R.B. Heat transfer through mixed convection boundary layer in a porous medium: LTNE analysis. Appl. Therm. Eng. 2020, 179, 115705. [Google Scholar] [CrossRef]
  9. Kannan, R.M.; Pullepu, B. Numerical solutions of double diffusive convective flow past a chemical reactive vertically inclined infinite plate with heat source/sink. Adv. Math. Sci. J. 2020, 9, 1623–1636. [Google Scholar] [CrossRef]
  10. Astanina, M.S.; Sheremet, M.; Umavathi, C.J. Unsteady natural convection in a partially porous cavity having a heat-generating source using local thermal non-equilibrium model. Int. J. Numer. Methods Heat Fluid Flow 2019, 29, 1902–1919. [Google Scholar] [CrossRef]
  11. Hema, M.; Shivakumara, I.S.; Ravisha, M. Double diffusive LTNE porous convection with Cattaneo effects in the solid. Heat Transf. 2020, 49, 3613–3629. [Google Scholar] [CrossRef]
  12. Nadian, P.K.; Pundir, R.; Pundir, S.K. Study of double–diffusive convection in a rotating couple stress ferromagnetic fluid in the presence of varying gravitational field and horizontal magnetic field saturating in a porous medium. J. Math. Comput. Sci. 2021, 11, 1784–1809. [Google Scholar]
  13. Shukla, S.; Gupta, U. LTNE effects on triple-diffusive convection in nanofluids. ASME J. Heat. Transf. 2022, 144, 092501. [Google Scholar] [CrossRef]
  14. Mahajan, A.; Tripathi, V.K. Stability of a chemically reacting double-diffusive fluid layer in a porous medium. Heat Transf. 2021, 50, 6148–6163. [Google Scholar] [CrossRef]
  15. Meften, G.A.; Ali, A.H. Continuous dependence for double diffusive convection in a Brinkman model with variable viscosity. Acta Univ. Sapientiae Math. 2022, 14, 125–146. [Google Scholar] [CrossRef]
  16. Capone, F.; Gianfrani, J.A. Natural convection in a fluid saturating an anisotropic porous medium in LTNE: Effect of depth-dependent viscosity. Acta Mech. 2022, 233, 4535–4548. [Google Scholar] [CrossRef]
  17. Capone, F.; Gianfrani, J.A. Thermal convection for a Darcy–Brinkman rotating anisotropic porous layer in local thermal non-equilibrium. Ric. Mat. 2022, 71, 227–243. [Google Scholar] [CrossRef]
  18. Abed Meften, G.; Ali, A.H.; Al-Ghafri, K.S.; Awrejcewicz, J.; Bazighifan, O. Nonlinear stability and linear instability of double-diffusive convection in a rotating with LTNE Effects and symmetric properties: Brinkmann-Forchheimer model. Symmetry 2022, 14, 565. [Google Scholar] [CrossRef]
  19. Akram, S.; Athar, M.; Saeed, K.; Umair, M.Y. Nanomaterials effects on induced magnetic field and double-diffusivity convection on peristaltic transport on Prandtl nanofluids in inclined asymmetric channel. Nanomater. Nanotechnol. 2022, 12, 18479804211048630. [Google Scholar] [CrossRef]
  20. Tayebi, T.; Dahmane, F.; Jamshed, W.; Chamkha, A.J.; El Din, S.M.; Raizah, Z. Double diffusive magneto-natural convection of nano fluid in an enclosure equipped with a wavy porous cylinder in the local thermal non-equilibrium situation. Case Stud. Therm. Eng. 2023, 43, 102785. [Google Scholar] [CrossRef]
  21. Siddabasappa, C.; Siddehwar, P.G.; Mallikarjunaih, S.M. Analytical study of Brinkman-Bènard convection in a bidisperse porous medium: Linear and weakly nonlinear study. Therm. Sci. Eng. Prog. 2023, 39, 101696. [Google Scholar] [CrossRef]
  22. Noon, N.J.; Haddad, S.A. Stability Analysis of double diffusive convection in local thermal non-equilibrium porous medium with Internal heat source and reaction effects. J. Non-Equilib. Thermodyn. 2023, 48, 25–39. [Google Scholar] [CrossRef]
  23. Corcione, M.; Quintino, A. Double-Diffusive effects on the onset of Rayleigh-Benard convection of water-based nanofluids. Appl. Sci. 2022, 12, 8485. [Google Scholar] [CrossRef]
  24. Gangadharaiah, Y.H.; Manjunatha, N.; Udhayakumar, R.; Almarri, B.; Elshenhab, A.M.; Honnappa, N. Darcy–Brinkman double diffusive convection in an anisotropic porous layer with gravity fluctuation and throughflow. Mathematics 2023, 11, 1287. [Google Scholar]
  25. Mehdy, A.N. Double diffusive free convection in a packed bed square enclosure by using local thermal non-equilibrium model. J. Eng. 2023, 18, 121–136. [Google Scholar] [CrossRef]
  26. Jakhar, A.; Kumar, A. Instability analysis of double diffusive convection under time dependent solute boundary conditions in the presence of internal heat generator. Phys. Fluids 2023, 35, 077101. [Google Scholar] [CrossRef]
  27. Yellamma; Narayanappa, M.; Udhayakumar, R.; Almarri, B.; Ramakrishna, S.; Elshenhab, A.M. The impact of heat source and temperature gradient on Brinkman–Bènard triple-diffusive magneto-Marangoni convection in a two-layer system. Symmetry 2023, 15, 644. [Google Scholar] [CrossRef]
  28. Rudolph, P.; Wang, W.; Tsukamoto, K.; Wu, D. Transport phenomena of crystal growth-heat and mass transfer. AIP Conf. Proc. 2010, 1270, 107–132. [Google Scholar]
  29. Sumithra, R.; Venkatraman, S. Local thermal non-equilibrium dominant Darcy-Rayleigh-Bènard-magneto-Marangoni convection in a composite layer. J. Mines Met. Fuels 2022, 70, 1–12. [Google Scholar] [CrossRef]
  30. Nield, D.A.; Bejan, A. Convection in Porous Media, 3rd ed.; Springer: New York, NY, USA, 2006. [Google Scholar]
  31. Manjunatha, N.; Khan, U.; Elattar, S.; Eldin, S.M.; Chohan, J.S.; Sumithra, R.; Sarada, K. Onset of triple-diffusive convective stability in the presence of a heat source and temperature gradients: An exact method. AIMS Math. 2023, 8, 13432–13453. [Google Scholar]
  32. Rudraiah, N.; Veerappa, B.; Balachandra, R.S. Effects of nonuniform thermal gradient and adiabatic boundaries on convection in porous media. J. Heat Transf. 1980, 102, 254. [Google Scholar] [CrossRef]
  33. Vasseur, P.; Robillard, L. The Brinkman model for natural convection in a porous layer: Effects of nonuniform thermal gradient. Int. J. Heat Mass. Transf. 1993, 36, 4199. [Google Scholar] [CrossRef]
  34. Shivakumara, I.S. Onset of convection in a couple-stress fluid-saturated porous medium: Effects of non-uniform temperature gradients. Arch. Appl. Mech. 2010, 80, 949–957. [Google Scholar] [CrossRef]
Figure 1. Model of the problem.
Figure 1. Model of the problem.
Mathematics 11 04282 g001
Figure 2. The effect of the fluid phase thermal expansion ratio α f p L = 1 , 5 , 10 on R C . (a) Linear, parabolic and inverted parabolic, (b) Heating from below, cooling from above and step function.
Figure 2. The effect of the fluid phase thermal expansion ratio α f p L = 1 , 5 , 10 on R C . (a) Linear, parabolic and inverted parabolic, (b) Heating from below, cooling from above and step function.
Mathematics 11 04282 g002
Figure 3. Influence of the solid phase thermal expansion ratio α s p L = 2 , 6 , 10 on R C . (a) Linear, parabolic and inverted parabolic, (b) Heating from below, cooling from above and step function.
Figure 3. Influence of the solid phase thermal expansion ratio α s p L = 2 , 6 , 10 on R C . (a) Linear, parabolic and inverted parabolic, (b) Heating from below, cooling from above and step function.
Mathematics 11 04282 g003
Figure 4. Impact of the porous parameter β = 0.8 , 0.9 , 1.0 on R C . (a) Linear, parabolic and inverted parabolic, (b) Heating from below, cooling from above and step function.
Figure 4. Impact of the porous parameter β = 0.8 , 0.9 , 1.0 on R C . (a) Linear, parabolic and inverted parabolic, (b) Heating from below, cooling from above and step function.
Mathematics 11 04282 g004
Figure 5. The effect of the fluid phase thermal diffusivity κ f p L = 0.2 , 0.4 , 0.5 , on R C . (a) Linear, parabolic and inverted parabolic, (b) Heating from below, cooling from above and step function.
Figure 5. The effect of the fluid phase thermal diffusivity κ f p L = 0.2 , 0.4 , 0.5 , on R C . (a) Linear, parabolic and inverted parabolic, (b) Heating from below, cooling from above and step function.
Mathematics 11 04282 g005
Figure 6. The effect of the variation of the solid phase thermal diffusivity ratio κ s p L = 0.3 , 0.5 , 0.7 , on R C . (a) Linear, parabolic and inverted parabolic, (b) Heating from below, cooling from above and step function.
Figure 6. The effect of the variation of the solid phase thermal diffusivity ratio κ s p L = 0.3 , 0.5 , 0.7 , on R C . (a) Linear, parabolic and inverted parabolic, (b) Heating from below, cooling from above and step function.
Mathematics 11 04282 g006
Figure 7. Influence of corrected internal Rayleigh number in region-I R I = 0.1 , 0.5 , 1.0 on R C . (a) Linear, parabolic and inverted parabolic, (b) Heating from below, cooling from above and step function.
Figure 7. Influence of corrected internal Rayleigh number in region-I R I = 0.1 , 0.5 , 1.0 on R C . (a) Linear, parabolic and inverted parabolic, (b) Heating from below, cooling from above and step function.
Mathematics 11 04282 g007
Figure 8. The influence of corrected internal Rayleigh number in region-II R I M = 0.1 , 0.5 , 1.0 on R C . (a) Linear, parabolic and inverted parabolic, (b) Heating from below, cooling from above and step function.
Figure 8. The influence of corrected internal Rayleigh number in region-II R I M = 0.1 , 0.5 , 1.0 on R C . (a) Linear, parabolic and inverted parabolic, (b) Heating from below, cooling from above and step function.
Mathematics 11 04282 g008
Figure 9. Impacts of solute Rayleigh numbers R f L s = R p L s = 10 , 20 , 30 on R C . (a) Linear, parabolic and inverted parabolic, (b) Heating from below, cooling from above and step function.
Figure 9. Impacts of solute Rayleigh numbers R f L s = R p L s = 10 , 20 , 30 on R C . (a) Linear, parabolic and inverted parabolic, (b) Heating from below, cooling from above and step function.
Mathematics 11 04282 g009
Table 1. Thermal gradients for fluid and porous layer.
Table 1. Thermal gradients for fluid and porous layer.
ModelThermal GradientsFluid Layer (Region-I)Porous Layer (Region-II)
Model (i)Linear F f L ( z f L ) = 1 F p L ( z p L ) = 1
Model (ii)Parabolic F f L ( z f L ) = 2 z f L F p L ( z p L ) = 2 z p L
Model (iii)Inverted parabolic F f L ( z f L ) = 2 ( 1 z f L ) F p L ( z p L ) = 2 ( 1 z p L )
Model (iv)Piecewise linear gradient heated from below (PLHB) F f L ( z f L ) = { ε f L 1 , 0 z f L ε f L 0 , ε f L z f L 1 F p L ( z p L ) = { ε p L 1 , 0 z p L ε p L 0 , ε p L z p L 1
Model (v)Piecewise linear gradient cooled from above (PLCA) F f L ( z f L ) = { 0 , 0 z f L ( 1 ε f L ) ε f L 1 , ( 1 ε f L ) z f L 1 F p L ( z p L ) = { 0 , 0 z p L ( 1 ε p L ) ε p L 1 , ( 1 ε p L ) z p L 1
Model (vi)Step function (SF) F f L ( z f L ) = δ ( z f L ε f L ) F p L ( z p L ) = δ ( z p L ε p L )
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Balaji, V.K.; Narayanappa, M.; Udhayakumar, R.; AlNemer, G.; Ramakrishna, S.; Honnappa, G.Y. Effects of LTNE on Two-Component Convective Instability in a Composite System with Thermal Gradient and Heat Source. Mathematics 2023, 11, 4282. https://doi.org/10.3390/math11204282

AMA Style

Balaji VK, Narayanappa M, Udhayakumar R, AlNemer G, Ramakrishna S, Honnappa GY. Effects of LTNE on Two-Component Convective Instability in a Composite System with Thermal Gradient and Heat Source. Mathematics. 2023; 11(20):4282. https://doi.org/10.3390/math11204282

Chicago/Turabian Style

Balaji, Varalakshmi K., Manjunatha Narayanappa, Ramalingam Udhayakumar, Ghada AlNemer, Sumithra Ramakrishna, and Gangadharaih Yeliyur Honnappa. 2023. "Effects of LTNE on Two-Component Convective Instability in a Composite System with Thermal Gradient and Heat Source" Mathematics 11, no. 20: 4282. https://doi.org/10.3390/math11204282

APA Style

Balaji, V. K., Narayanappa, M., Udhayakumar, R., AlNemer, G., Ramakrishna, S., & Honnappa, G. Y. (2023). Effects of LTNE on Two-Component Convective Instability in a Composite System with Thermal Gradient and Heat Source. Mathematics, 11(20), 4282. https://doi.org/10.3390/math11204282

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop