Next Article in Journal
Test and Mesoscopic Analysis of Chloride Ion Diffusion of High-Performance-Concrete with Fly Ash and Silica Fume
Next Article in Special Issue
Electroplating on Al6082 Aluminium: A New Green and Sustainable Approach
Previous Article in Journal
Degradation and Stabilization of Resin-Dentine Interfaces in Polymeric Dental Adhesives: An Updated Review
Previous Article in Special Issue
Optimisation of Thiourea Concentration in a Decorative Copper Plating Acid Bath Based on Methanesulfonic Electrolyte
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Throughflow and Coriolis Force on the Onset of Double-Diffusive Convection with Internal Heat Source

by
Kallu Vetty Muhammed Rafeek
1,
Gudala Janardhana Reddy
1,*,
Ravi Ragoju
2,
Gundlapally Shiva Kumar Reddy
2 and
Mikhail A. Sheremet
3,4,*
1
Laboratory on Computational Fluid Dynamics, Department of Mathematics, Central University of Karnataka, Kalaburagi 585367, Karnataka, India
2
Department of Applied Sciences, National Institute of Technology Goa, Ponda 403401, Goa, India
3
Laboratory on Convective Heat and Mass Transfer, Tomsk State University, 634050 Tomsk, Russia
4
Butakov Research Center, National Research Tomsk Polytechnic University, 634050 Tomsk, Russia
*
Authors to whom correspondence should be addressed.
Coatings 2022, 12(8), 1096; https://doi.org/10.3390/coatings12081096
Submission received: 10 April 2022 / Revised: 15 May 2022 / Accepted: 28 July 2022 / Published: 2 August 2022

Abstract

:
The present research examines the joint influence of throughflow and Coriolis force on the onset of double-diffusive convection with an internal heat source modelled by Darcy’s law. The BVP4C routine in MATLAB R2020a is used to solve the eigenvalue problem numerically. Critical Rayleigh numbers are obtained for designated values of governing parameters. The effect of the internal heat source parameter, Taylor number, Darcy number, and Peclet number on the system’s stability is investigated. The internal heat source parameter has a destabilizing influence on the system, according to our findings. The reason behind this observation is that the presence of an internal heat source in the porous medium may cause more molecular diffusion inside the medium. The Taylor number, on the other hand, stabilizes the system for both upward and downward throughflow because rotation introduces vorticity into the fluid. Thus, the fluid moves with higher velocity in horizontal planes. The velocity of the fluid perpendicular to the planes reduces as a result of this motion. Thus, the onset of convection is delayed.

1. Introduction

Convection with an internal heat source has been of interest to many researchers, since it has extensive applications in astrophysics [1], combustion modelling [2], geophysics [3], the miniaturization of electronic components [4] and thermal ignition [5]. Non-linear temperature circulation in the system is steered by the presence of an internal heat source, and convection may occur. Tritton and Zarraga [6] conducted the first experimental research on internally heated thermal convection, which was followed by theoretical studies by Roberts [7] and Thirby [8]. The origins of this phenomenon were established by Tveitereid and Palm [9]. Takashima ([10,11]) explored convection in a biased fluid layer with internal heat. They discovered that an inner heat source has a significant impact on the stability of the onset of convection. Natural convection affected by the heat source with the outcome of heat source distribution was tested by Tasaka and Takeda [12]. The analysis of free convection in a porous medium plays a substantial role in many areas, such as geophysical systems, the Earth’s oceans, magma chambers, the petroleum industry and many engineering applications. Nield and Bejan [13] provided a brief review of this topic.
There are many studies related to geophysical sciences and technological applications which contain non-isothermal motion of liquids, known as throughflow. With the height of the fluid, this flow transforms the basic state temperature equation from linear to non-linear, altering the system’s stability. Jones and Persichetti [14] discussed the thermal instability in packed beds with throughflow. Wooding [15] studied the Rayleigh instability of a thermal boundary layer flow in saturated porous medium, in which he showed that the layer is stable provided that the Rayleigh number for the system does not exceed a critical positive value and that the wave number of the critical neutral disturbance is finite. Homsy and Sherwood [16] investigated the linear and energy theory on thermal instability in porous media with throughflow. They showed that the fluid can lose stability through either a buoyantly driven mode or through a continuous analogue of the Saffman–Taylor mode. Sutton [17] and Shivakumara [18] explored the impact of throughflow on the onset of convection in a horizontal porous layer.
The effect of throughflow and internal heat generation on the onset of convection in a porous material was examined by Khalili and Shivakumara [19]. They concluded that, if the boundaries are of the same type, throughflow destabilizes the system, but this is not true when internal heat source is absent. In a porous material with a tilted temperature drop and vertical throughflow, Brevdo [20] described 3D absolute and convective instability at the onset of convection. He deduced the fact that for marginally supercritical values of the vertical Rayleigh number, the destabilization has the character of absolute instability in all of the cases in which the horizontal Rayleigh number is zero or the Peclet number is zero. Shivakumara and Sureshkumar [21] investigated convective instability in non-Newtonian fluids in vertical throughflow, and concluded that throughflow has an essential influence depending on the nature of the borders and fluid flow directions. Yadav [22] scrutinized the impacts of throughflow and a varied gravity field on the onset of convective flow in a porous medium layer numerically by employing the higher order Galerkin method, and showed that both the throughflow and gravity variation parameters postpone the onset of convective motion. Later, many researchers such as Kiran [23], Bhadauria and Kiran [24], Kiran [25], Shinkumara and Nanjundappa [26], Reza and Gupta [27], Nield and Kuznetsov [28], Yadav [29] and Kiran and Bhadauria [30] investigated the effect of throughflow with different external effects.
Convection driven by the internal heating of porous material was investigated by Gasser and Kazimi [31] and Tveitereid [32]. Yadav et al. [33] performed linear analysis and used the Galerkin method to explain the onset of convection in rotating porous media with an inner heater. Mahabaleshwar et al. [34] analyzed the convection heat transfer in a porous zone with modulated gravity and an inner heater. Some interesting results can be found in [35,36,37].
Riahi [38] investigated nonlinear convection in a porous region with an inner heater and discovered that the non-uniform internal heat source could reduce or enhance the ideal Rayleigh value and cell size. In the case of linear and nonlinear stability, Rionero and Straughan [39] derived a critical Rayleigh value. The internal heat-generating porous medium in vertical cavities was investigated by Du and Bilgen [40]. Hewitt et al. [41] overworked the underlying theory, mechanism, correlations and methodologies of heat transfer. Choi et al. [42] investigated a variety of characteristics of convection flow in porous medium caused and prolonged by a constant inner heater. Brinkman convection in a rotating porous zone filled with a nanofluid with an inner heater was investigated by Yadav et al. [43]. An internal heater’s influence on the onset of convection in a porous medium filled with a nanosuspension was studied by Yadav et al. [44]. Barletta et al. [45] studied the influence of viscous dissipation in the porous material. The stability of mixed thermal convection in a porous zone was discussed by Sphaier et al. [46] using the generalized integral transform technique. The effects of inner thermal production and throughflow on convective instability in an anisotropic porous medium were studied by Vanishree [47].
The present paper aims to study the thermal convection stability with an internal heater in porous material. The plan of this research is as follows. In Section 2, we describe the considered problem. In Section 3, we discuss the basic state. In Section 4, the linear instability analysis is performed. The method of solution is described in Section 5. The numerical outcomes and discussions are presented in Section 6. The research ends with conclusions.

2. Governing Equations

We consider heat conducting liquid in a porous zone placed between two infinitely parallel horizontal plates at z = 0 and z = L which are set to rotate at a fixed angular velocity Ω ¯ = Ω e ¯ z . The upper and lower bounding surfaces of the layer are assumed to be stress-free. The z-axis is oriented upward, so that g = g e ¯ z where g’ is the modulus of g’ and e ¯ z is the unit vector along the z-direction in Figure 1. Physical properties of the fluid are assumed to be constant, except density in the buoyancy term, so that the Boussinesq approximation is valid. The control Oberbeck–Boussinesq equations are [48]:
u = 0
u = P + D a 2 u + R a θ e ¯ z + T a u × e ¯ z
θ t + u · θ = 2 θ + Q
with the following boundary conditions
u = 0 ,   θ = 1   on   z = 0 u = θ = 0   on   z = 1 .
where u is the velocity, t is the time, θ is the temperature, Da is the Darcy number, Ra is the Rayleigh number, Ta is the Taylor number and Q is the parameter of inner heater, defined as follows:
D a = µ e µ K L 2 , R a = ρ 0 g β Δ T L 3 α µ , Q = Q L 2 k Δ T , T a = 2 ρ 0 Ω 0 K μ ϕ 2 ,   P = K α µ P ¯ + ρ 0 g z .  
The rescaling used to obtain Equations (1)–(3) with conditions (4) is
x , y , z = x L , y L , z L ,   u , v , w = L u α , L v α , L w α ,   t = α L 2 A t ,   α = k ρ c p f ,
where asterisks refer to dimensional quantities, L is the channel height, A is the ratio of volumetric thermal capacity of liquid filled porous material to the fluid, i.e., A = ρ c m ρ c p f , α is the thermal diffusivity, k is the thermal conductivity, K is the permeability, ρ0 is the mean flow density, μ is the dynamic viscosity, ϕ is the porosity, ω is the angular velocity, μe is the effective viscosity, Q′ > 0 is a (fixed) inner heater and β is the thermal expansion coefficient.

3. Basic State

The basic steady motion of Equations (1)–(3) is defined by a uniform throughflow,
u b = 0 ,   v b = 0 ,   w b = P e
where b stands for the basic state, P e = w 0 L α is the Peclet number and w0 is the prescribed vertical throughflow velocity. The Peclet number is positive for upward throughflow and negative for downward throughflow. By substituting Equation (6) into Equation (3), we obtain the basic temperature profile as
θ b = e z P e e P e   1 e P e + z + e z P e z e P e 1 Q 1 e P e P e  

4. Linear Stability Study

The perturbation of the basic state can be defined as
u = u b + ε u ,   θ = θ b + ε θ ,   P = P b + ε P .
where ε 1 . By substituting Equation (8) into Equations (1)–(4) and by neglecting terms O ε 2 or higher, we obtain the linearized governing equations as follows:
u = P ¯ + D a 2 u + R a θ e ¯ z + T a u × e ¯ z
θ t + F z w + P e θ z = 2 θ
w = θ = 0   at   z = 0 ,   1
where F z = Q P e + e z P e P e + Q 1 e P e   . By taking the third components of the curl of (9) and the double curl of (9), we obtain
1 D a 2 ω z T a w z = 0
1 D a 2 2 w R a h 2 θ + T a ω z z = 0
where ω z = × u e ^ z . By removing ωz from Equations (12) and (13), one obtains
1 D a 2 2 2 w R a 1 D a 2 h 2 θ + T a 2 w z 2 = 0
Normal modes can be defined by the perturbations
w , θ = W z ,   θ z e i m x + l y + σ t
where q = (m, l, 0) is the wave vector, with q = m 2 + l 2 expressing the wave number, and σ is a complex characteristic, where its real part, σr, is the raising rate of instability and its imaginary part, σi, is the angular frequency. Substituting the above expression into Equations (10) and (14), we obtain:
1 D a D 2 q 2 2 D 2 q 2 W + R a 1 D a D 2 q 2 q 2 θ + T a D 2 W = 0
D 2 q 2 P e D i ω θ W F z = 0
W = D 2 W = D 4 W = θ = 0   on   z = 0 ,   1 .
The principle of exchange of stabilities is used. In other words, the marginal stability condition can be found for stationary modes. Hence, Equations (16)–(18) become
1 D a D 2 q 2 2 D 2 q 2 W + R a 1 D a D 2 q 2 q 2 θ + T a D 2 W = 0
D 2 q 2 P e D θ W F z = 0
W = D 2 W = D 4 W = θ = 0   on   z = 0 ,   1 .

5. Solution Methodology

The eigenvalue problem, defined by Equations (19)–(21), is worked out by employing the BVP4C routine in MATLAB R2020a. To achieve a non-trivial solution to the eigenvalue problem, the normalizing condition w 0 = 1 is considered. We compute the eigenvalue Ra using this normalization condition. The critical Rayleigh and wave numbers are acquired by using index-linked instruction in MATLAB R2020a. The comparative and conclusive patience were taken as 10−6 and 10−10 independently to achieve higher-order exactness.

6. Results and Discussion

This section contains the numerical results and discussions. A numerical study of the eigenvalue problem corresponding to a convection problem with the uniform internal heat source and throughflow was performed in this paper. The non-dimensional parameters governing the onset of convection are the Rayleigh number, Ra, inner heater parameter, Q, Taylor number, Ta, Darcy number, Da, and Peclet number, Pe. The BVP4C routine in Matlab R2020a is used to work out the eigenvalue problem for linear stability analysis.
We verified our results with those found in the literature to validate our analysis method. In the absence of throughflow and rotation, the current problem reduces to that of Gasser and Kazimi [31] for Darcy porous media. Table 1 demonstrates that our numerical results are really close to Gasser and Kazimi’s critical external Rayleigh number. Furthermore, the current numerical results are validated by comparing them to those found by Barletta et al. [48] for Q = 0, Ta = 0 and Da = 0 (see Table 2).
Table 3 illustrates the Rayleigh number critical values for various Q and Pe, for constant Da = 0.01, Da = 0.1, and Ta = 50. A graphical representation of these values is given in Figure 2. The results for upward throughflow are shown in Figure 2. The Rac reduces with a growth of Q in this figure, suggesting that the internal heat source parameter has a destabilizing influence on the system. The reason behind this observation is that the presence of an internal heat source in the porous medium may cause more molecular diffusion inside the medium.
Table 4 shows the Rac for different Q and Pe, along with fixed Da = 0.01, Da = 0.1 and Ta = 50. Figure 3 illustrates a visual behavior of these values. The results for downward throughflow are shown in Figure 3. As can be seen in this figure, the critical Rayleigh number reduces as Q increases, implying that Q has a destabilizing influence on the system.
The behavior of Rac versus the Taylor number is shown in Table 5 and Figure 4. Figure 4 shows that as the Taylor number rises, the critical Rayleigh number rises as well, indicating the stabilizing impact of the Taylor number. This can be explained as follows: rotation introduces vorticity into the fluid. Thus the fluid moves in horizontal planes with higher velocity. On account of this motion, the velocity of the fluid perpendicular to the planes reduces. Thus, the onset of convection is delayed.
Table 6 and Figure 5 show the dependence of Rac on Ta for various Pe. Only the results for downward throughflow are shown in this part. The critical Rayleigh number rises with the increase in Ta, as shown in Figure 5, indicating that the Taylor number has a stabilizing impact.

7. Conclusions

This study has examined the linear instability of rotation convection in a porous zone with an inner heater. The behavior of various parameters such as the inner heater coefficient Q, critical Rayleigh number Rac, Peclet number Pe, Taylor number Ta, and Darcy number Da has been analyzed. The following are the most important findings from the linear instability:
-
In the absence of throughflow and rotation, the Rac and the wave number for the Darcy porous medium match with those found in the literature and reported by Gasser and Kazimi [31].
-
In the absence of an inner heater and rotation, the critical values of Ra for the Darcy porous medium are identical to those discovered by Barletta et al. [48].
-
The system is destabilized by the internal heat source parameter.
-
The Taylor number has a stabilizing impact on the considered unit for both upward and downward throughflows.

Author Contributions

Conceptualization, G.J.R.; methodology, K.V.M.R., G.J.R. and R.R.; software, R.R. and G.J.R.; writing—original draft preparation, K.V.M.R., G.J.R., R.R., G.S.K.R. and M.A.S.; writing—review and editing, K.V.M.R., G.J.R., R.R., G.S.K.R. and M.A.S.; visualization, K.V.M.R., G.J.R. and R.R. All authors have read and agreed to the published version of the manuscript.

Funding

The first author wishes to express his gratitude to the Department of Science and Technology, Government of India for the granting of DST-Inspire Fellowship (IF190169) (Grant no. C/4084/IFD/2020-21) and to the Central University of Karnataka for providing the research facilities. This research of the fifth author was supported by the Tomsk Polytechnic University Development Programme (Priority-2030-NIP/EB-002-0000-2022).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data can be found in this paper.

Acknowledgments

The authors are grateful to editor and anonymous referees for their suggestions and comments that helped in the enhancement of the clarity and depth of the interpretation of the results in the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tritton, D.J. Internally heated convection in the atmosphere of Venus and in the laboratory. Nature 1975, 257, 110–112. [Google Scholar] [CrossRef]
  2. Hurley, M.J.; Gottuk, D.T.; Hall, J.R., Jr.; Harada, K.; Kuligowski, E.D.; Puchovsky, M.; Watts, J.M., Jr.; Wieczorek, C.J. (Eds.) SFPE Handbook of Fire Protection Engineering; Springer: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
  3. Mckenzie, D.P.; Roberts, J.M.; Weiss, N.O. Convection in the Earth’s mantle: Towards a numerical simulation. J. Fluid Mech. 1974, 62, 465–538. [Google Scholar] [CrossRef]
  4. Chadwick, M.L.; Webb, B.W.; Heaton, H.S. Natural convection from two-dimensional discrete heat sources in a rectangular enclosure. Int. J. Heat Mass Transf. 1991, 34, 1679–1693. [Google Scholar] [CrossRef]
  5. May, H.O. A numerical study on natural convection in an inclined square enclosure containing internal heat sources. Int. J. Heat Mass Transf. 1991, 34, 919–928. [Google Scholar] [CrossRef]
  6. Tritton, D.J.; Zarraga, M.N. Convection in horizontal layers with internal heat generation. Experiments. J. Fluid Mech. 1967, 30, 21–31. [Google Scholar] [CrossRef]
  7. Roberts, P.H. Convection in horizontal layers with internal heat generation. Theory. J. Fluid Mech. 1967, 30, 33–49. [Google Scholar] [CrossRef]
  8. Thirlby, R. Convection in an internally heated layer. J. Fluid Mech. 1970, 44, 673–693. [Google Scholar] [CrossRef]
  9. Tveitereid, M.; Palm, E. Convection due to internal heat sources. J. Fluid Mech. 1976, 76, 481–499. [Google Scholar] [CrossRef] [Green Version]
  10. Takashima, M. The stability of natural convection in an inclined fluid layer with internal heat generation. J. Phys. Soc. Jpn. 1989, 58, 4431–4440. [Google Scholar] [CrossRef]
  11. Takashima, M. The stability of natural convection in an inclined fluid layer with internal heat generation. II. J. Phys. Soc. Jpn. 1991, 60, 455–465. [Google Scholar] [CrossRef]
  12. Tasaka, Y.; Takeda, Y. Effects of heat source distribution on natural convection induced by internal heating. Int. J. Heat Mass Transf. 2005, 48, 1164–1174. [Google Scholar] [CrossRef]
  13. Nield, D.A.; Bejan, A. Convection in Porous Media, 3rd ed.; Springer: New York, NY, USA, 2013. [Google Scholar]
  14. Jones, M.C.; Persichetti, J.M. Convective instability in packed beds with throughflow. AIChE J. 1986, 32, 1555–1557. [Google Scholar] [CrossRef]
  15. Wooding, R.A. Rayleigh instability of a thermal boundary layer in flow through a porous medium. J. Fluid Mech. 1960, 9, 183–192. [Google Scholar] [CrossRef]
  16. Homsy, G.M.; Sherwood, A.E. Convective instabilities in porous media with throughflow. AIChE J. 1976, 22, 168–174. [Google Scholar] [CrossRef]
  17. Sutton, F.M. Onset of convection in a porous channel with net throughflow. Phys. Fluids 1970, 13, 1931. [Google Scholar] [CrossRef]
  18. Shivakumara, I.S. Effects of throughflow on convection in porous media. Proc. Seventh Asian Congr. Fluid Mech. 1997, 2, 557–560. [Google Scholar]
  19. Khalili, A.; Shivakumara, I.S. Onset of convection in a porous layer with net throughflow and internal heat generation. Phys. Fluids 1998, 10, 315. [Google Scholar] [CrossRef]
  20. Brevdo, L. Three-dimensional absolute and convective instabilities at the onset of convection in a porous medium with inclined temperature gradient and vertical throughflow. J. Fluid Mech. 2009, 641, 475–487. [Google Scholar] [CrossRef]
  21. Shivakumara, I.S.; Sureshkumar, S. Convective instabilities in a viscoelastic fluid-saturated porous medium with throughflow. J. Geophys. Eng. 2007, 4, 104–115. [Google Scholar] [CrossRef]
  22. Yadav, D. Numerical investigation of the combined impact of variable gravity field and throughflow on the onset of convective motion in a porous medium layer. Int. Commun. Heat Mass Transf. 2019, 108, 104274. [Google Scholar] [CrossRef]
  23. Kiran, P. Throughflow and gravity modulation effects on heat transport in a porous medium. J. Appl. Fluid Mech. 2016, 93, 1105–1113. [Google Scholar] [CrossRef]
  24. Bhadauria, S.; Kiran, P. Chaotic and oscillatory magneto-convection in a binary viscoelastic fluid under g-jitter. Int. J. Heat Mass Transf. 2015, 84, 610–624. [Google Scholar] [CrossRef]
  25. Kiran, P. Throughow and g-jitter effects on binary fluid saturated porous medium. Appl. Math. Mech. 2015, 36, 1285–1304. [Google Scholar] [CrossRef]
  26. Shinkumara, I.S.; Nanjundappa, C.E. Effects of quadratic drag and throughflow on double diffusive convection in a porous layer. Int. Commun. Heat Mass Transf. 2006, 33, 357–363. [Google Scholar] [CrossRef] [Green Version]
  27. Reza, M.; Gupta, A.S. Magnetohydrodynamic thermal instability in a conducting fluid layer with throughflow. Int. J. Non-Linear Mech. 2012, 47, 616–625. [Google Scholar] [CrossRef]
  28. Nield, D.A.; Kuznetsov, A.V. The effects of combined horizontal and vertical heterogeneity on the onset of convection in a porous medium with horizontal throughflow. Int. J. Heat Mass Transf. 2011, 54, 5595–5601. [Google Scholar] [CrossRef]
  29. Yadav, D. The onset of Darcy-Brinkman convection in a porous medium layer with vertical throughflow and variable gravity field effects. Heat Transf. 2020, 49, 3161–3173. [Google Scholar] [CrossRef]
  30. Kiran, P.; Bhadauria, S. Throughflow and rotational effects on oscillatory convection with modulation. Nonlinear Stud. 2016, 23, 439–455. [Google Scholar]
  31. Gasser, R.D.; Kazimi, M.S. Onset of convection in a porous medium with internal heat generation. J. Heat Transf. 1976, 98, 49–54. [Google Scholar] [CrossRef]
  32. Tveitereid, M. Thermal convection in a horizontal porous layer with internal heat sources. Int. J. Heat Mass Transf. 1977, 20, 1045–1050. [Google Scholar] [CrossRef]
  33. Yadav, D.; Wang, J.; Jinho, L. Onset of Darcy-Brinkman convection in a rotating porous layer induced by purely internal heating. J. Porous Media 2017, 20, 691–706. [Google Scholar] [CrossRef]
  34. Mahabaleshwar, U.S.; Basavaraja, D.; Wang, S.; Lorenzini, G.; Lorenzini, E. Convection in a porous medium with variable internal heat source and variable gravity. Int. J. Heat Mass Transf. 2017, 111, 651–656. [Google Scholar] [CrossRef]
  35. Gholamalizadeh, E.; Pahlevanzadeh, F.; Ghani, K.; Karimipour, A.; Nguyen, T.K.; Safaei, M.R. Simulation of water/FMWCNT nanofluid forced convection in a microchannel filled with porous material under slip velocity and temperature jump boundary conditions. Int. J. Numer. Methods Heat Fluid Flow 2020, 30, 2329–2349. [Google Scholar] [CrossRef]
  36. Heydari, A.; Akbari, O.A.; Safaei, M.R.; Derakhshani, M.; Alrashed, A.A.A.A.; Mashayekhi, R.; Ahmadi Sheikh Shabani, G.; Zarringhalam, M.; Nguyen, T.K. The effect of attack angle of triangular ribs on heat transfer of nanofluids in a microchannel. J. Therm. Anal. Calorim. 2018, 131, 2893–2912. [Google Scholar] [CrossRef]
  37. Maleki, H.; Safaei, M.R.; Togun, H.; Dahari, M. Heat transfer and fluid flow of pseudo-plastic nanofluid over a moving permeable plate with viscous dissipation and heat absorption/generation. J. Therm. Anal. Calorim. 2019, 135, 1643–1654. [Google Scholar] [CrossRef]
  38. Riahi, N. Non-linear convection in a porous medium with internal heat sources. Int. J. Non-Linear Mech. 1984, 19, 469–478. [Google Scholar] [CrossRef]
  39. Rionero, S.; Straughan, B. Convection in a porous medium with internal heat source and variable gravity effects. Int. J. Eng. Sci. 1990, 28, 497–503. [Google Scholar] [CrossRef]
  40. Du, Z.G.; Bilgen, E. Natural convection in vertical cavities with internal heat generating porous medium. Heat Mass Transf. 1992, 27, 149–155. [Google Scholar] [CrossRef]
  41. Hewitt, G.F.; Shires, G.L.; Bott, T.R. Process Heat Transfer; CRC Press: Boca Raton, FL, USA, 1994. [Google Scholar]
  42. Choi, E.; Chamkha, A.; Nandakumar, K. A bifurcation study of natural convection in porous media with internal heat sources: The non-Darcy effects. Int. J. Heat Mass Transf. 1998, 41, 383–392. [Google Scholar] [CrossRef]
  43. Yadav, D.; Lee, J.; Cho, H.H. Brinkman convection induced by purely internal heating in a rotating porous medium layer saturated by a nanofluid. Powder Technol. 2015, 286, 592–601. [Google Scholar] [CrossRef]
  44. Yadav, D.; Bhargava, R.; Agrawal, G.S. Boundary and internal heat source effects on the onset of Darcy–Brinkman convection in a porous layer saturated by nanofluid. Int. J. Therm. Sci. 2012, 60, 244–254. [Google Scholar] [CrossRef]
  45. Barletta, A.; Celli, M.; Rees, D.A. The onset of convection in a porous layer induced by viscous dissipation: A linear stability analysis. Int. J. Heat Mass Transf. 2009, 52, 337–344. [Google Scholar] [CrossRef]
  46. Sphaier, L.A.; Barletta, A.; Celli, M. Unstable mixed convection in a heated inclined porous channel. J. Fluid Mech. 2015, 778, 428–450. [Google Scholar] [CrossRef]
  47. Vanishree, R.K. Effects of throughflow and internal heat generation on a thermo convective instability in an anisotropic porous medium. J. Appl. Fluid Mech. 2014, 7, 581–590. [Google Scholar]
  48. Barletta, A.; Rossi di Schio, E.; Storesletten, L. Convective roll instabilities of vertical throughflow with viscous dissipation in a horizontal porous layer. Transp. Porous Media 2010, 81, 461–477. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic of the problem.
Figure 1. Schematic of the problem.
Coatings 12 01096 g001
Figure 2. Dependence of Rac on Q for (a) Da = 0.01 and (b) Da = 0.1.
Figure 2. Dependence of Rac on Q for (a) Da = 0.01 and (b) Da = 0.1.
Coatings 12 01096 g002
Figure 3. Dependence of Rac on Q for (a) Da = 0.01 and (b) Da = 0.1.
Figure 3. Dependence of Rac on Q for (a) Da = 0.01 and (b) Da = 0.1.
Coatings 12 01096 g003
Figure 4. Dependence of Rac on Ta for (a) Da = 0.01 and (b) Da = 0.1.
Figure 4. Dependence of Rac on Ta for (a) Da = 0.01 and (b) Da = 0.1.
Coatings 12 01096 g004
Figure 5. Dependence of Rac on Ta for (a) Da = 0.01 and (b) Da = 0.1.
Figure 5. Dependence of Rac on Ta for (a) Da = 0.01 and (b) Da = 0.1.
Coatings 12 01096 g005
Table 1. Critical Rayleigh number values for Ta = 0; Da = 0 and Pe = 0 compared with [31].
Table 1. Critical Rayleigh number values for Ta = 0; Da = 0 and Pe = 0 compared with [31].
QGasser and Kazimi [31]Present Study
039.4839.4788
534.5934.5953
1027.0227.0162
1521.4521.4436
2017.6317.6267
2514.9214.9165
3012.9112.9117
4010.1610.1606
508.378.3690
607.117.1121
805.475.4670
1004.444.4391
Table 2. Critical Rayleigh number values for Q = 0; Ta = 0 and Da = 0 compared with [48].
Table 2. Critical Rayleigh number values for Q = 0; Ta = 0 and Da = 0 compared with [48].
PeBarletta et al. [48]Present StudyPeBarletta et al. [48]Present Study
−0.00139.478439.478420.00139.478439.47842
−0.0139.478639.478560.0139.478639.47856
−0.139.492439.492370.139.492439.49237
−140.875140.87507140.875140.87507
−245.077645.07761245.077645.07761
−352.068452.06842352.068452.06842
−461.666461.66642461.666461.66642
−573.414673.41456573.414673.41456
−686.619286.61920686.619286.61920
−7100.581100.580857100.581100.58085
−8114.833114.832608114.833114.83260
−9129.167129.166859129.167129.16685
−10143.518143.5184910143.518143.51849
−15215.283215.2828015215.283215.28280
Table 3. Critical Rayleigh number for upward throughflow fixed at Ta = 50.
Table 3. Critical Rayleigh number for upward throughflow fixed at Ta = 50.
QDa = 0.01Da = 0.1
Pe = 0.001Pe = 0.01Pe = 0.1Pe = 0.001Pe = 0.01Pe = 0.1
1470.7949470.8283471.2445360.2181360.2767360.9591
2456.8040456.8651457.5531353.1931353.3036354.5031
3436.5979436.6802437.5732342.6116342.7646344.3875
4413.2066413.3041414.3438329.7012329.8865331.8297
5389.0000389.1083390.2508315.5820315.7903317.9600
6365.4518365.5675366.7803301.1062301.3299303.6501
7343.3180343.4387344.6974286.8434 287.0766289.4889
8322.8993323.0231324.3092273.1351 273.3735275.8350
9304.2439304.3692305.6686260.1618 260.4023262.8818
10287.2738287.3996288.7016247.9998 248.2401250.7157
11271.8562271.9817273.2784236.6618 236.9005239.3566
12257.8414257.9660259.2515226.1242 226.3602228.7861
13245.0812245.2044246.4745216.3440 216.5764218.9648
14233.4375233.5590234.8107207.2690 207.4974209.8434
15222.7858222.9053224.1365198.8445 199.0685201.3693
16213.0155213.1329214.3422191.0165 191.2360193.4896
17204.0293204.1446205.3312183.7335 183.9484186.1542
18195.7422195.8553197.0187176.9479 177.1581179.3200
19188.0798188.1907189.3307170.6157 170.8213172.9315
20180.9771181.0857182.2025164.6970 164.8982166.9615
Table 4. Critical Rayleigh number for downward throughflow fixed at Ta = 50.
Table 4. Critical Rayleigh number for downward throughflow fixed at Ta = 50.
QDa = 0.01Da = 0.1
Pe = −0.001Pe = −0.01Pe = −0.1Pe = −0.001Pe = −0.01Pe = −0.1
1470.7877470.7561470.5231360.2053 360.1488359.6803
2456.7906456.7312456.2149353.1688 353.0604352.0715
3436.5798436.4992435.7636342.5778 342.4269341.0094
4413.1850413.0890412.1932329.6602 329.4769327.7328
5388.9761388.8691387.8591315.5359 315.3296313.3517
6365.4262365.3117364.2219301.0567 300.8349298.6994
7343.2913343.1717342.0281286.7918 286.5603284.3257
8322.8720322.7493321.5716273.0823 272.8456270.5550
9304.2162304.0919302.8952260.1085 259.8696257.5545
10287.2460287.1211285.9168247.9465 247.7077245.3905
11271.8285271.7039270.5004236.6089 236.3718234.0684
12257.8139257.6902256.4938226.0720 225.8375223.5590
13245.0539244.9316243.7470216.2925 216.0615213.8155
14233.4106233.2900232.1205207.2184 206.9914204.7831
15222.7593222.6405221.4885198.7949 198.5722196.4048
16212.9894212.8727211.7397190.9678 190.7496188.6251
17204.0038203.8892202.7763183.6859 183.4722181.3917
18195.7171195.6047194.5126176.9012 176.6922174.6559
19188.0552187.9450186.8739170.5701 170.3656168.3734
20180.9530180.8450179.7950164.6525164.4525162.5039
Table 5. Critical Rayleigh number for upward throughflow fixed at Q = 2.
Table 5. Critical Rayleigh number for upward throughflow fixed at Q = 2.
TaDa = 0.01Da = 0.1
Pe = 0.001Pe = 0.01Pe = 0.1Pe = 0.001Pe = 0.01Pe = 0.1
5111.6703111.7037112.0815148.2322 148.3047149.0853
10162.3775162.4197162.8861180.1151 180.1960181.0679
15207.2615207.3097207.8364207.7227 207.8099208.7501
20248.5374248.5894249.1605232.6185 232.7106233.7058
25287.2261287.2787287.8852255.5851255.6815256.7225
30323.9144323.9714324.6015277.0832 277.1830278.2636
35358.9829359.0413359.6915297.4115297.5145298.6298
40392.6975392.7571393.4230316.7779 316.8836318.0299
45425.2544425.3149425.9932335.3336335.4418336.6159
50456.8040456.8651457.5531353.1931 353.3036354.5031
55487.4646487.5261488.2217370.4458370.5583371.7811
60517.3322517.3939518.0952387.1627 387.2771388.5215
65546.4853546.5471547.2528403.4019403.5180404.7824
70574.9901575.0518575.7605419.2113 419.3290420.6121
75602.9020602.9637603.6743434.6310434.7502436.0508
Table 6. Critical Rayleigh number for downward throughflow fixed at Q = 2.
Table 6. Critical Rayleigh number for downward throughflow fixed at Q = 2.
TaDa = 0.01Da = 0.1
Pe = −0.001Pe = −0.01Pe = −0.1Pe = −0.001Pe = −0.01Pe = −0.1
5111.6629111.6289111.3183148.2162148.1449147.4869
10162.3681162.3262161.9452180.0972180.0177179.2843
15207.2509207.2037206.7770207.7035207.6178206.8289
20248.5259248.4750248.0171232.5982232.5076231.6749
25287.2114287.1590286.6808255.5639255.4693254.6004
30323.9019323.8463323.3516277.0612276.9631276.0635
35358.9700358.9130358.4081297.3889297.2878296.3616
40392.6844392.6263392.1149316.7546316.6508315.7013
45425.2412425.1823424.6673335.3097335.2036334.2333
50456.7906456.7312456.2149353.1688353.0604352.0715
55487.4512487.3915486.8757370.4210370.3107369.3048
60517.3187517.2588516.7450387.1375387.0254386.0041
65546.4718546.4120545.9013403.3764403.2626402.2272
70574.9766574.9167574.4103419.1854419.0700418.0217
75602.8885602.8289602.3276434.6048434.4880433.4276
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rafeek, K.V.M.; Reddy, G.J.; Ragoju, R.; Reddy, G.S.K.; Sheremet, M.A. Impact of Throughflow and Coriolis Force on the Onset of Double-Diffusive Convection with Internal Heat Source. Coatings 2022, 12, 1096. https://doi.org/10.3390/coatings12081096

AMA Style

Rafeek KVM, Reddy GJ, Ragoju R, Reddy GSK, Sheremet MA. Impact of Throughflow and Coriolis Force on the Onset of Double-Diffusive Convection with Internal Heat Source. Coatings. 2022; 12(8):1096. https://doi.org/10.3390/coatings12081096

Chicago/Turabian Style

Rafeek, Kallu Vetty Muhammed, Gudala Janardhana Reddy, Ravi Ragoju, Gundlapally Shiva Kumar Reddy, and Mikhail A. Sheremet. 2022. "Impact of Throughflow and Coriolis Force on the Onset of Double-Diffusive Convection with Internal Heat Source" Coatings 12, no. 8: 1096. https://doi.org/10.3390/coatings12081096

APA Style

Rafeek, K. V. M., Reddy, G. J., Ragoju, R., Reddy, G. S. K., & Sheremet, M. A. (2022). Impact of Throughflow and Coriolis Force on the Onset of Double-Diffusive Convection with Internal Heat Source. Coatings, 12(8), 1096. https://doi.org/10.3390/coatings12081096

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