ORIGINAL RESEARCH article
Sec. Colloidal Materials and Interfaces
Volume 10 - 2023 | https://doi.org/10.3389/fmats.2023.1158644
Darcy–Benard–Oldroyd convection in anisotropic porous layer subject to internal heat generation
- 1Department of Mathematics, Government College (Autonomous), Kalaburagi, India
- 2Department of Mathematics, School of Applied Sciences, REVA University, Bangalore, India
- 3Department of Mathematical Sciences, Faculty of Science and Technology, Universiti Kebangsaan Malaysia, Bangi, Selangor, Malaysia
- 4Department of Mathematics and Social Sciences, Sukkur IBA University, Sukkur, Sindh, Pakistan
- 5Department of Mathematics, Sharnbasva University, Kalaburagi, India
- 6Mechanical Engineering, Future University in Egypt, New Cairo, Egypt
- 7Department of Mathematics, College of Science, Qassim University, Buraydah, Saudi Arabia
- 8University Centre for Research and Development and Department of Mechanical Engineering, Chandigarh University, Mohali, Punjab, India
An anisotropic horizontal porous layer saturated with viscoelastic liquids of the Oldroyd-B type is explored to determine how the internal heat source affects thermal convection. As a momentum equation, a modified Darcy–Oldroyd model is used that takes into account the anisotropy of the porous layer. The energy equation is formulated in such a way that the influence of internal heat sources and anisotropy in thermal diffusivity on the stability criterion may be easily identified. The effects of anisotropy, viscoelasticity, and internal heat generation on the onset of thermal convection are investigated using linear stability analysis. It is understood that convection begins via an oscillatory mode instead of a stationary mode because viscous relaxation, thermal diffusions, and internal heat generation mechanisms compete with one another. Both steady and unsteady finite-amplitude convections are studied using nonlinear stability analysis with the truncated Fourier series method. The effect of different governing parameters on the system’s stability and on convective heat transfer is studied. The present investigation has been significantly validated by the recovery of several prior results as special situations. The findings presented in this work are anticipated to have significant implications for a number of real-world applications, including modeling of oil reservoirs, crude oil extraction, crystal growth, the pharmaceutical and medical industries, and the use of geothermal energy, among others.
Numerous real-world applications exist for the theoretical and practical analysis of convective heat transfer in porous media. Geothermal energy systems, hydrocarbon reserves, nuclear reactors, medicine, and the chemical industry are among the many examples. Ingham and Pop (2005), Nield and Bejan (2006), Vafai (2005), and Vadasz (2008) have documented the developments in porous medium thermal convection. An essential component of rheological research is the study of the viscoelastic properties of the asthenospheric and mantle components of the Earth (see Lowrie, 2020). Snow systems and the rheology of food transport involve viscoelastic liquids saturated in porous media. Less research has been conducted on Rayleigh–Benard convection in liquids with viscoelastic properties than on thermal convection in a Newtonian fluid (refer to Li and Khayat, 2005).
Due to their coagulated viscosity, polymeric liquids are unaffected by flow and turbulence, whereas Newtonian fluids are affected. Moreover, oscillatory convection is believed to exist as viscoelastic fluids are characterized by their elasticity. According to Kim et al. (2003), the system is unstable due to its elasticity but stable due to its porosity. The oscillatory mode is appropriate for studying convection because convection reduced to supercritical and stable bifurcation forms does not vary with elasticity. However, few authors have investigated oscillatory convection in viscoelastic liquids, and many researchers have not addressed a comparable problem in porous media. Thermal convection in viscoelastic fluids was studied by O’Connell and Budiansky (1977), Griffiths (1987), Rudraiah et al. (1989), Yoon et al. (2004), Malashetty et al. (2006a), and Swamy et al. (2012).
When a porous medium exhibits anisotropy in its mechanical and thermal properties, the flow behavior is significantly altered. The assumption underlying Darcy’s model is that the fluid flow is sufficiently slow for inertial effects to be negligible, and that the fluid complies with the laws of continuum mechanics. If the porous medium is anisotropic, then the medium’s permeability will depend on the flow direction. In other words, the permeability of the porous medium will vary in different orientations. Darcy’s law must be modified to account for the anisotropy of the medium in this instance. This can be accomplished by introducing a permeability tensor that describes the medium’s permeability in various orientations. Numerous applications in the actual world involve anisotropic porous media. Anisotropy results from the incorrect orientation of a solid matrix or the asymmetries of the natural porous medium; it is also a characteristic of porous synthetic materials, such as the fibrous material utilized for insulation and pelleting, both of which are beneficial to chemical engineering processes.
McKibbin et al. (1985) and Storesletten (1998) have exhaustively documented the work pertaining to convection in anisotropic porous media. Govender (2006) described the effects of anisotropy on thermal convection in a porous layer. Malashetty et al. (2006b)presented a study on the effect of a time-periodic modulated temperature field on the stability of a viscoelastic fluid saturated within an anisotropic porous layer. Saravanan and Purusothaman (2009) studied non-Darcian effects in anisotropic porous media. Many pertinent studies on anisotropy have been conducted (Malashetty and Swamy, 2007a; Malashetty and Swamy, 2007b; Malashetty et al., 2009; Sivakumar and Saravanan, 2009; Malashetty and Swamy, 2010; Agarwal et al., 2011; Malashetty et al., 2011; Srivastava et al., 2011; Chandra and Satyamurty, 2012; Swamy et al., 2013; Swamy et al., 2014; Swamy, 2017; Swamy et al., 2019; Capone et al., 2020).
Understanding thermal convection in a porous stratum with internal heat generation is crucial in a variety of natural and artificial systems. The interior of the Earth provides an example of internal heat production. Due to the disintegration of radioactive isotopes and residual heat from the planet’s formation, the Earth’s core is thought to be the source of significant internal heating. This heat production is responsible for the elevated temperature of the Earth’s interior and is one of the driving forces behind plate tectonics and volcanism. Electronic devices, such as computer processors, are another example of internal heat generation. A computer processor’s electronic components generate heat due to the passage of electric current through them. If the generated heat is not effectively dissipated, the processor may malfunction or even fail. To prevent this, electronic devices are frequently equipped with cooling systems, such as heat sinks and fans, to dissipate heat and maintain secure component temperatures. Chemical reactions can also generate heat internally. Due to the exothermic nature of the combustion reaction, for example, heat is produced when fuel is consumed. This internal heat production can be utilized to generate electricity in power facilities or to heat industrial processes.
Internal heat generation results in a nonlinear temperature gradient. A temperature gradient is the temperature change over a specified distance. A nonlinear temperature gradient indicates that the change in temperature over a given distance is not constant. In other terms, the change in temperature is not linear. This refers to the fact that the system or material in question generates its own heat rather than receiving it from an external source. Therefore, thermal convection occurs even though the temperature difference between the lower and upper surfaces is insufficient for convection to begin. Consequently, the production of internal heat is an additional essential mechanism for regulating the onset of convection in the porous layer. Numerous researchers have conducted extensive research on thermal convection in absorbent layers with internal heat generation [refer to (Thirlby, 1970; Mahabaleshwar et al., 2017; Ahmed and Rashed, 2019; Yadav et al., 2021; Enagi et al., 2022; Raju et al., 2022; Upadhya et al., 2022)].
Now, let us see the physical mechanisms involved in the present problem. As a result of the internal heating, a temperature gradient is generated within the porous layer, with the temperature being higher near the heat source and decreasing toward the top surface of the layer. This temperature gradient creates density differences within the viscoelastic liquid, causing it to flow from the heat source toward the top surface. The fluid flow of the viscoelastic liquid is influenced by the anisotropy of the porous layer. If the porous layer is more permeable in one direction than in other directions, the fluid flow will be preferentially directed along the more permeable direction, resulting in anisotropic convection. The properties of the porous medium, such as its porosity, permeability, and tortuosity, also affect the rate of fluid flow and the way in which heat is transferred through the medium. For example, a higher porosity of the porous layer can increase the flow rate of the viscoelastic liquid, while a higher tortuosity can hinder the flow, leading to a slower rate of fluid motion. The viscoelasticity of the liquid can also influence the way in which fluid flow and heat transfer occur within the porous layer. A viscoelastic liquid has both viscous and elastic properties, which can lead to nonlinear behavior in fluid flow and heat transfer. For instance, the elasticity of the liquid can cause it to exhibit oscillatory behavior in response to the temperature gradient, resulting in oscillatory fluid flow and heat transfer. The heating of the bottom surface by an external means can also affect the flow of the viscoelastic liquid. For instance, uniformly heating the bottom surface can generate a temperature gradient that interacts with the temperature gradient generated by the internal heat source, resulting in a more complex flow pattern.
A stability analysis of an Oldroyd-B fluid in a porous medium with the combined influence of anisotropy and the internal heat source is rare. Such investigations are still greatly desired. The primary objective of this study is to determine how the viscoelastic parameters, internal heat generation coefficient, and anisotropy parameters impact the onset criterion of thermal convection and the heat transfer across the layer.
2 Mathematical formulation
A shallow horizontal anisotropic porous layer saturated with Oldroyd-B liquid is considered. The surfaces held at
Because the fluid is at rest in the basic state, we can determine its mass, pressure, and temperature by
The stability of the system is studied by imposing the perturbations on the basic state.
where prime represents the quantity in the perturbed state. We use Eqs. 5–8 in the governing Eqs. 1–4 and eliminate
3 Galerkin weighted-residual technique
In the linear stability theory, the imposed perturbations are anticipated to be infinitesimal, and hence, the nonlinear terms in Eqs. 9–11 can be ignored. The normal mode analysis is used to solve the resulting eigenvalue problem. Thus, it is supposed that the infinitesimal perturbations are the periodic waves of form
where the real variable k is the wavenumber, and the complex variable
In natural convection, the spatial variation of temperature and fluid flow patterns can take on different wavelengths, depending on factors such as the geometry of the system, the temperature differences within the fluid, and the properties of the fluid itself. The wavelength of these patterns can be quantified in terms of the wavenumber, and the behavior of natural convection can be analyzed in terms of the relationship between the wavenumber and other parameters such as the Rayleigh number.
According to the Galerkin method, we choose
where A1 and B1 are constants, and
The trial functions that satisfy the boundary conditions (15) are obviously
On using these, one can estimate X1, X2, X4, X5 and X6 value and then substituting into Eq. 20 to get
Now, we discuss the behavior of Eq. 19 with the nonzero growth rate, that is,
The real part of Eq. 19 then symbolizes the expression for the oscillatory Rayleigh number:
4 Weak nonlinear theory
Nonlinear stability analysis is preferred to measure convection amplitudes and heat transfer. This facilitates comprehension of the physical mechanism with a little mathematical labor. This is a basic step toward grasping the full nonlinearity of the problem. Because the perturbations are assumed to be of finite amplitude, it is reasonable to represent them in the form of a limited Fourier series, as follows:
The finite amplitudes of A and B subscripts for unsteady nonlinear convection are to be assessed by the dynamics of the system. Using Eqs. 24 and 25 in Eqs. 9 and 10 and comparing the coefficients of like terms, the following fourth-order Lorenz system of autonomous nonlinear differential equations is obtained:
There is no suitable analytical method to obtain a closed-form solution of the aforementioned system. Thus, a competent numerical technique is recommended. Although it may not be possible to make precise quantitative predictions, there are several qualitative insights that can be gleaned from the available data or theoretical models. As system of Eq. 26 is homogeneously bounded with time, it retains numerous features of the full problem. For
This implies that the system is constrained. In dynamical systems theory, a system is said to be constrained if it is subject to certain limitations or conditions. These constraints can arise due to physical, mathematical, or other reasons. For example, a mechanical system may be constrained by rigid walls or other physical barriers that restrict the motion of its components. When a system is constrained, the possible states that it can occupy are limited to a subset of its phase space. The phase space is the space of all possible states of the system, and it is often represented as a high-dimensional space in which each dimension corresponds to a particular variable or parameter of the system. Because the system is constrained, the phase space paths are drawn toward a set of measure zero or a fixed point. A set of measure zero is a subset of the phase space that has zero volume or probability. In other words, it is a set of states that is extremely unlikely to be reached by the system. A fixed point, on the other hand, is a state of the system that does not change over time. When the phase space paths are drawn toward a set of measure zero or a fixed point, this can result in volume shrinkage in the phase space. This means that the volume of the phase space that is accessible to the system becomes smaller over time as the system is forced to occupy a smaller subset of the phase space due to the constraints. This is revealed by Eq. 28 through the fact that if a set of preliminary points in space fills a volume
The larger values of the Darcy–Prandtl number and strain-retardation number and smaller values of the Deborah number and heat-capacities ratio are used to emphasize the exponential deterioration of volume with time.
Upon switching from qualitative exploration, we now look into the existence of an analytical solution. As the finite amplitude instability can be well explored analytically using the truncated model shown in Eqs. 24 and 25, a closed-form solution of Eq. 26 is used for the steady case. The following expression is obtained by Eq. 26 after removing all terms except
Rather than merely defining the onset criterion, the impact of the Rayleigh number can be swiftly documented by analyzing its effect on heat transport. In the study of convection, determining the quantity of heat transported past the layer is of utmost importance. Because the basic state is immobile, heat transfer in this state is limited to conduction. However, as the Rayleigh number exceeds the threshold, convection develops. The Nusselt number is used to describe the convection-heat transport throughout the stratum.
Substituting the value of
This can be further simplified as
This analysis is valid for Rayleigh numbers close to their threshold value. By including more terms in the Fourier expansion, one can anticipate better outcomes. The Runge–Kutta–Gill method is used to clarify the unsteady Eq. 26. The calculated amplitude values for various time intervals are then used to estimate Nu as a function of t.
5 Results and discussion
The onset of convection refers to the point at which fluid motion due to buoyancy forces begins to occur in a fluid that is initially at rest. This occurs when a fluid is subjected to a temperature difference that is large enough to cause density variations within the fluid, which in turn generate buoyancy forces that drive fluid motion. The onset of convection can be characterized by a critical value of a dimensionless parameter called the Rayleigh number. The Rayleigh number represents the ratio of buoyancy forces to viscous forces in the fluid. For a fluid layer that is initially at rest, the onset of convection occurs when the Rayleigh number exceeds a critical value that is specific to the geometry and boundary conditions of the system. Above this critical value, the buoyancy forces overcome the viscous forces and initiate fluid motion. The critical Rayleigh number can be determined through theoretical analysis or experimental observation. The onset of convection is an important phenomenon in many natural and industrial systems, including geophysical flows, crystal growth, and industrial heat transfer processes. Understanding the onset of convection is important for predicting the behavior of these systems and for designing efficient heat transfer systems.
The primary objective of investigating a convection problem is to determine the smallest possible Rayleigh number that demonstrates convection. Fine-tuning the parameters that define the Rayleigh number is advantageous for deferring or accelerating convective motions. This study examines the interaction between variable permeability, thermal diffusivity, and internal heat generation. The point
In Figure 1B, the effect of the mechanical anisotropy parameter, which signifies heterogeneity in the permeability of the porous layer, is expressed. The values of both
One of the main objectives of the present paper is to analyze the significance of controlling the convection by the anisotropic nature of the porous layer. The detailed behavior of
The impact of
To explore the activity of
This study helped us to measure the amplitude of convective motions and the amount of heat transfer. The Nusselt number (Nu) that signifies the extent of convective heat transfer is calculated. The control of Rayleigh number over Nu is presented in Figures 3A–C. The value of Nu is found to be 1 at the onset, while as the
From Figures 3A,B, Nu is found to upsurge for the higher
In unsteady finite-amplitude analysis, the fourth-order Lorenz model has been solved numerically using the Runge–Kutta–Gill method. The amplitudes are obtained as the functions of t and then substituted into the expression of Nu. In general, the Nusselt number is a function of several parameters, including the fluid velocity, temperature, and physical properties of the fluid, and it can vary with time during transient heat transfer processes. The behavior of the heat transfer coefficient with respect to time is displayed through the curves in the (Nu, t) plane, as shown in Figures 4A–D. One can observe from these figures that at the onset of natural convection, heat transfer initially occurs through conduction alone, and Nu has a value of 1, which corresponds to purely conductive heat transfer. As the fluid flow starts to develop, convective heat transfer becomes dominant, and the Nusselt number becomes sensitive to time. The behavior of the Nusselt number during this transient phase is oscillatory, with the value of Nu fluctuating about a mean value. This means that the heat transfer coefficient, which is related to the Nusselt number, varies periodically in time, with some oscillations. However, as time progresses, the oscillatory behavior of the Nusselt number starts to decay, and after a short period of time, the heat transfer process reaches a steady state. At this point, the Nusselt number reaches a mean value that is similar to the value of Nu obtained in the steady-state, finite-amplitude case. Thus, in general, these figures display that the behavior of the Nusselt number during a transient heat transfer process is initially oscillatory and sensitive to time but eventually reaches a steady-state mean value. This behavior is a result of the interplay between conduction and convection heat transfer mechanisms during the transient phase of the process.
Figures 4A, B depict a considerable enhancement in the value of the heat transfer coefficient with the Prandtl number and the heat capacities ratio. It also shows that the sensitivity of Nu with respect to t increases with increasing Pr and
The aforementioned results support the subsequent conclusions. Due to a competition between thermal diffusion, viscoelastic relaxation, anisotropy in permeability, thermal diffusivity, and internal heat generation, the conduction state degenerated into convective motions via the oscillatory mode. Viscoelasticity, the heat capacity ratio, and the Prandtl number have no effect on stationary convection. Increasing anisotropy in permeability, the coefficient of internal heat generation, and stress-relaxation parameters are associated with an early onset. It has been discovered that anisotropy in thermal diffusivity and heat capacity ratios delays convection. Strain-retardation time reinforces oscillatory case stability. The range of retardation parameter values within which oscillatory convection occurs is determined by the magnitude of the relaxation time. Beyond this range, oscillatory convection ceases, and stationary mode instability is then established. The Prandtl number indicates the effect of destabilization on oscillatory convection. The coefficient of heat transfer increases with the Rayleigh number, the internal Rayleigh number, the Deborah number, the Prandtl number, anisotropy in permeability, and the heat capacity ratio. Increasing values of the strain-retardation parameter and anisotropy in thermal diffusivity indicate an inverse trend. It is possible to promote or inhibit convection in a given system by adjusting the relevant parameters in accordance with practical application. In other words, the onset and intensity of convection can be manipulated by adjusting the system’s controlling parameters, such as temperature gradients, fluid properties, and geometrical factors. Several of the prior results were obtained through the use of special cases. This provides substantial support for the results of the current investigation (Eswaramoorthi et al., 2023).
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
MS, BH, AS: conceptualization, methodology, software, formal analysis, and writing—original draft. AH: writing—original draft, data curation, investigation, visualization, and validation. UK: conceptualization, writing—original draft, writing—review and editing, supervision, and resources. RK: validation, investigation, writing—review and editing, and formal analysis. VK: writing—review and editing, data curation, validation, and resources. All authors contributed to the article and approved the submitted version.
This work was partially funded by the research center of the Future University in Egypt 2023.
The authors would like to express their gratitude to the reviewers for their constructive remarks and insightful suggestions, which have significantly improved the current work.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Agarwal, S., Bhadauria, B. S., and Siddheshwar, P. G. (2011). Thermal instability of a nanofluid saturating a rotating anisotropic porous medium. Spec. Top. Rev. Porous Media 2 (1), 53–64. doi:10.1615/specialtopicsrevporousmedia.v2.i1.60
Ahmed, S. E., and Rashed, Z. Z. (2019). MHD natural convection in a heat generating porous medium-filled wavy enclosures using Buongiorno's nanofluid model. Case Stud. Therm. Eng. 14, 100430. doi:10.1016/j.csite.2019.100430
Eswaramoorthi, S., Loganathan, K., Faisal, M., Botmart, T., and Shah, N. A. (2023). Analytical and numerical investigation of Darcy-Forchheimer flow of a nonlinear-radiative non- Newtonian fluid over a Riga plate with entropy optimization. Ain Shams Eng. J. 14 (3), 101887. doi:10.1016/j.asej.2022.101887
Mahabaleshwar, U. S., Basavaraja, D., Wang, S., Lorenzini, G., and Lorenzini, E. (2017). Convection in a porous medium with variable internal heat source and variable gravity. Int. J. Heat. Mass Transf. 111, 651–656. doi:10.1016/j.ijheatmasstransfer.2017.04.030
Malashetty, M. S., Shivakumara, I. S., Kulkarni, S., and Swamy, M. (2006a). Convective instability of Oldroyd-B fluid saturated porous layer heated from below using a thermal non-equilibrium model. Transp. Porous Med. 62, 123–139. doi:10.1007/s11242-005-1893-0
Malashetty, M. S., Siddheshwar, P. G., and Swamy, M. (2006b). Effect of thermal modulation on the onset of convection in a viscoelastic fluid saturated porous layer. Transp. Porous Media 62, 55–79. doi:10.1007/s11242-005-4507-y
Malashetty, M. S., Swamy, M. S., and Sidram, W. (2011). Double diffusive convection in a rotating anisotropic porous layer saturated with viscoelastic fluid. Int. J. Therm. Sci. 50, 1757–1769. doi:10.1016/j.ijthermalsci.2011.04.006
Malashetty, M. S., and Swamy, M. (2007b). The effect of rotation on the onset of convection in a horizontal anisotropic porous layer. Int. J. Therm. Sci. 46, 1023–1032. doi:10.1016/j.ijthermalsci.2006.12.007
Malashetty, M. S., Tan, W. C., and Swamy, M. (2009). The onset of double diffusive convection in a binary viscoelastic fluid saturated anisotropic porous layer. Phys. Fluids 21 (8), 084101. doi:10.1063/1.3194288
McKibbin, R. (1985). “Thermal convection in layered and anisotropic porous media: A review,” in Convective flows in porous media. Editors R. A. Wooding, and I. White (Wellington, NZ: Department of Scientific and Industrial Research), 113–127.
Enagi, N. K., Chavaraddi, K. B., Kulkarni, S., and Ramesh, G. K. (2022). Effect of maximum density and internal heating on the stability of rotating fluid saturated porous layer using LTNE model. Heliyon 8, e09620. doi:10.1016/j.heliyon.2022.e09620
Raju, C. S. K., Ahammad, N. A., Sajjan, K., Shah, N. A., Yook, S. J., and Kumar, M. D. (2022). Nonlinear movements of axisymmetric ternary hybrid nanofluids in a thermally radiated expanding or contracting permeable Darcy Walls with different shapes and densities: Simple linear regression. Int. Comm. Heat. Mass 135, 106110. doi:10.1016/j.icheatmasstransfer.2022.106110
Saravanan, S., and Purusothaman, A. (2009). Floquet instability of a gravity modulated Rayleigh–Benard problem in an anisotropic porous medium. Int. J. Therm. Sci. 48, 2085–2091. doi:10.1016/j.ijthermalsci.2009.04.001
Srivastava, A. K., Bhadauria, B. S., and Kumar, J. (2011). Magnetoconvection in an anisotropic porous layer using thermal nonequilibrium model. Spec. Top. Rev. Porous Media 2 (1), 1–10. doi:10.1615/specialtopicsrevporousmedia.v2.i1.10
Swamy, M. S. (2017). Combined effect of thermal modulation and AC electric field on the onset of electrothermoconvection in anisotropic porous layer. Am. J. Heat Transf. 4 (3), 95–114. doi:10.7726/ajhmt.2017.1011
Swamy, M. S., Naduvinamani, N. B., and Sidram, W. (2012). Onset of Darcy-Brinkman convection in a binary viscoelastic fluid saturated porous layer. Transp. Porous Med. 94, 339–357. doi:10.1007/s11242-012-0008-y
Swamy, M. S., Patil, S., and Pallavi, S. P. (2019). Soret and dufour effect induced double- diffusive reaction-convection in anisotropic porous layer. J. Nanofluids 8, 1329–1337. doi:10.1166/jon.2019.1688
Swamy, M. S., Shivakumara, I. S., and Naduvinamani, N. B. (2014). Effect of gravity modulation on electrothermal convection in dielectric fluid saturated anisotropic porous layer. J. Heat. Transf. 136, 032601.
Swamy, M. S., Shivakumara, I. S., and Sidram, W. (2013). The onset of convection in a gravity-modulated viscoelastic fluid-saturated anisotropic porous layer. Spec. Top. Rev. Porous Media 4 (1), 69–80. doi:10.1615/specialtopicsrevporousmedia.v4.i1.70
Upadhya, S. M., Raju, S. V. S. R., Raju, C. S. K., Shah, N. A., and Chung, J. D. (2022). Importance of entropy generation on Casson, Micropolar and Hybrid magneto-nanofluids in a suspension of cross diffusion. Chin. J. Phys. 77 (2022), 1080–1101. doi:10.1016/j.cjph.2021.10.016
Yadav, D., Mahabaleshwar, U. S., Wakif, A., and Chand, R. (2021). Significance of the inconstant viscosity and internal heat generation on the occurrence of Darcy Brinkman convective motion in a couple-stress fluid saturated porous medium an analytical solution. Int. Comm. Heat. Mass Transf. 122, 105165. doi:10.1016/j.icheatmasstransfer.2021.105165
Keywords: thermal convection, internal heat generation, Rayleigh number, Oldroyd-B model, viscoelastic fluid
Citation: Swamy MS, Hanumagowda BN, Khan U, Vidyashree K, Hassan AM, Mohammed Saeed A and Kumar R (2023) Darcy–Benard–Oldroyd convection in anisotropic porous layer subject to internal heat generation. Front. Mater. 10:1158644. doi: 10.3389/fmats.2023.1158644
Received: 04 February 2023; Accepted: 01 June 2023;
Published: 15 June 2023.
Edited by:Mustafa Inc, Firat University, Türkiye
Reviewed by:Francesco Costanzo, The Pennsylvania State University (PSU), United States
Jawali C. Umavathi, Gulbarga University, India
Copyright © 2023 Swamy, Hanumagowda, Khan, Vidyashree, Hassan, Mohammed Saeed and Kumar. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ahmed M. Hassan, email@example.com