1. Introduction
The convection of a binary mixture, kept in a porous medium, provides a rich analysis for the study of non-equilibrium phenomena. The antagonistic relationship between the flow-driven by buoyancy effects, namely the density difference caused by the temperature and concentration variations, is often the novel aspect of convection studies. Such phenomena are generally referred to as the natural convection of combined heat and mass transfer, thermohaline or double-diffusive convection. In many engineering applications and problems related to seawater and mantle in the crust of Earth, such double-diffusive convective flow often occurs. Some more examples include the distribution of chemical pollutants in water-saturated soil, grain storage facilities and the migration of moisture through the air contained in insulated fibers. The oscillatory instability is a characteristic feature of double-diffusive convection. In different branches of fluid mechanics and condensed matter physics, the convective systems which exhibit oscillatory convection as a first instability show the interesting physical phenomena. In simple fluids, only the oscillatory instability can occur as secondary bifurcation. Since higher bifurcations are more difficult to describe theoretically, it will be useful to have a real physical system that shows oscillatory instability as the first bifurcation.
The oscillatory instability depends on the stabilizing effect of imposed concentration gradient, such a gradient can also be developed in response to the applied temperature gradient. The first instability with the spatial structure of the convection rolls vary from the stationary bifurcation to the oscillatory bifurcation as the control parameters of the system keeps changing. Specifically, when the superfluid is heated from below in horizontal convective system, it becomes unstable with the oscillatory convection [
1]. It is interesting to observe that oscillatory instability occurs along with the concentration perturbations, which stabilizes the system and the temperature perturbations destabilize the convective system. At a certain critical value of temperature gradient, the relaxation time of the temperature fluctuations is much larger than that of the concentration fluctuations.
Further research concerning the normal fluid mixture subjected to the vertical temperature gradient was investigated by Steinberg [
2] using the linear stability analysis. He pointed out the existence of oscillatory instabilities when the mixture is heated from below. Later, the influence of two-fluid effects near the onset of convection in
–
mixture was investigated by Steinberg and Brand [
3]. They observed that the potential chemical perturbations reduced the stability of system. Similar flow–visualization experiments in the ethanol-water mixture were performed by Kolodner et al. [
4]. They have shown that the conducting state became unstable due to the existence of oscillatory convection. In this convective system, apart from the stability, the occurrence of forward and backward bifurcations over the particular range of were shown by Guenter and Ingo [
5] using heat transport measurements. In the ethanol-water and
–
mixtures, Barbara et al. [
6] experimentally presented the data for the critical Rayleigh number and neutral frequency as a function of separation ratio. Based on the sign and magnitude of Soret coefficient, Brand and Steinberg [
7] showed that when a layer of two miscible fluids embedded in a porous medium was heated from below or above, any one of two bifurcations will occur, namely, stationary or oscillatory instabilities. Later, Steinberg and Brand [
8] illustrated a similar result for a layer of two miscible fluids in a porous medium, that is heated from below or above and with a rapid chemical reaction. Further, the amplitude equation was derived by Brand and Steinberg [
9] near the threshold for both stationary and oscillatory instabilities. In a simple fluid, they have shown when heating is conducted from below for the Rayleigh-Bénard convection, the first instability is always stationary. Later, depending on the importance of occurrence of parameters in these systems, Steinberg and Brand [
10] predicted the existence of forwarding bifurcation for stationary convection, as well as forwarding and inverse bifurcations for oscillatory convection.
It can be observed from the above studies that the natural convection due to the double diffusion alone has been widely studied and well-reported, while very few studies have been devoted to the double-diffusive convection in porous media. Thermohaline convection is an example of double-diffusive convection. Double-diffusive convection in porous media is the focus of an extensive study due to its significance in predicting the movement of groundwater in aquifers and the method of energy extraction from geothermal reservoirs. Brand and Hohenberg [
11] derived an amplitude equation using the weakly nonlinear analysis for the double-diffusive convection and discussed the experimental realizability in alcohol-water and binary-fluid (
–
) mixtures. For such systems, Rehberg and Ahlers [
12] elaborated on the occurrence of steady or oscillatory bifurcations based on the mean temperature. The thresholds for finite-amplitude, oscillatory and monotonic convective instabilities are calculated both analytically and numerically by Bahloul et al. [
13] in terms of the governing parameters of these systems. Using a Galerkin method, Bahloul et al. [
14] investigated the nonlinear state of ideal straight rolls arising in these systems. They observed that the stability regions of these rolls are restricted by the cross roll, zigzag and oscillatory instabilities.
For a horizontal layer saturated in a porous medium, the oscillatory instability may be possible when a strongly stabilizing solute gradient is opposed by a destabilizing thermal gradient (Nield [
15]). An order of magnitude predictions for the total heat and mass transfer rates in such systems and their valid domains are shown to be in agreement with the results produced by discrete numerical experiments (Trevisan and Bejan [
16,
17]). Near the onset, when the negative separation ratio (
) decreases, the traveling wave bifurcation branch is shifted towards the lower values of
R without qualitative changes as long as the Lewis number (
) is small and the normalized porosity is large (Augustin et al. [
18]). The heat and concentration flux values were obtained experimentally by Griffiths [
19] through an interface between two-fluid layers and with different temperature values and salt concentrations. He showed that the differences in horizontal properties in each convecting layer lead to an increased heat flow through the two-layer device due to the bending of interface. Hounsou et al. [
20] studied the criterion of the appearance of stationary convection in a porous horizontal layer saturated by a binary mixture of ferro-fluids, which is heated from below. The effect of Darcy number, a ratio of viscosities, the magnetic and binary parameters are illustrated by them on the appearance and size of the convection cells. Recently Hu and Zhang [
21] discovered a new oscillatory instability with a positive separation ratio. This ratio showed the interactions among several factors and provided the analytical criteria to predict the onset of Rayleigh -Bénard convection under the ideal boundary conditions. Rashmi and Murthy [
22] have established the nonlinear stability theory for a uniform flow along a Brinkman porous layer with two impermeable and isothermal horizontal boundaries. The cross-diffusion effect is also considered as a factor contributing to the convective instability in the medium.
The nonlinear theory of Rayleigh-Bénard convection was carried out by Kuo [
23] using the Fourier analysis of perturbations which are valid for all the Rayleigh numbers
(critical stationary convection Rayleigh number) and he compared his theoretical results with the experimental observations. For the convection in porous media, the finite amplitude convection was studied by Palm et al. [
24] and their theoretical results showed very good agreement with the experiments. Rameshwar et al. [
25,
26] and Rawoof Syeed et al. [
27], have investigated the steady-state finite amplitude cellular convection for the double diffusive physical models. It can be noted that the research carried by Rameshwar et al. [
26] is based on stationary convection with the eigenvalue,
and the present research work is based on oscillatory convection with eigenvalue,
. The dynamic behaviors of stationary convection and oscillatory convection are different in nature. For stationary convection (steady) model, the approximate solutions until
(
is the expansion parameter) are computed, but in the present oscillatory convection (unsteady) model, due to the high complexity of analytical expressions of the eigenfunctions, we are able compute the approximate solutions until
.
Much less attention is paid to the nonlinear studies of the considered complex physical model, while near the onset of oscillatory convection, only weakly nonlinear studies exist in the literature. Thus, the objectives of the present study were as follows:
Investigate theoretically the linear and nonlinear behaviour of the Rayleigh-Bénard convection system of a superfluid mixture, –, kept in a porous medium;
Obtain the critical values of control parameters such as R, from the linear stability analysis for the onset of stationary and oscillatory convection;
Find the R values for different bifurcations viz. pitchfork bifurcation, Hopf bifurcation, Takens–Bogdanov bifurcation and co-dimension two bifurcation;
Solve the nonlinear partial differential equations using the perturbation method proposed by Kuo [
23], until the
and obtain the approximate solutions (eigenfunctions) to analyse the time-dependent nonlinear behavior of the convective system;
Find the local () and time–dependent averaged () Nusselt numbers on the hot wall to understand the development of heat flow and rate of heat transfer, respectively;
Obtain the cellular pattern of the fluid flow, hot regions (isotherms), concentration regions (isoconcentrations) from the eigenfunctions related to streamfunction, temperature and concentration, respectively;
Study the heatline patterns of the flow by using the heatfunction;
Study the effect of physical parameters on the entropy generation.
The manuscript is arranged, further, into nine sections.
Section 2 lists the hydrodynamic governing equations related to the double-diffusive flow in a porous medium under the Boussinesq approximation. For this mathematical model, in
Section 3, the linear stability analysis is presented. In
Section 4, the method of solution for this model by using the perturbation method is elaborated. The approximate solutions for the field variables are also presented in this section. In
Section 5, the approximate analytical expression for
is given. In
Section 6, the streamlines and isotherm patterns of the flow are presented.
Section 7 and
Section 8 deal with the heatline patterns of the flow by using the heat function and the entropy generation, respectively. Finally, the conclusions of the present study are given in
Section 9.
2. Mathematical Formulation
The binary mixture
–
, placed in a sparsely packed porous medium of thickness
d was considered. This layer is parallel to the horizontal
-plane with a very large horizontal extension that is kept in the gravitational field
. The layer has its interfaces at the vertical coordinates
and
. A static temperature difference across the layer was assumed to be imposed. The temperature on the bottom plate
was assumed to be
and on the top plate
the temperature is
. The basic equations were taken according to the Boussinesq approximation and the Darcy model. Based on the above assumptions, the momentum equation is given by [
7]:
In above equation, the quantity
‘a’ measures the microscale inertial effects such as thermal and particle diffusion anisotropies. The inertial term
is associated with
‘a’. In the present study, the term
‘a’ is neglected since the contribution of anisotropy is assumed to be very small [
17,
28]. The other notations are explained in the Nomenclature. The other hydrodynamic basic equations for the conservation of mass, the concentration, and the temperature are given by [
7,
29]:
where
is the Soret coefficient,
for liquids; additionally,
. Near the onset of convection, Equations (
1)–(
4) give rise to the static solution as
where
is the applied temperature gradient. The Equations (
1)–(
4) are non-dimensionalized using the scales listed below,
The perturbed velocity, temperature and concentration quantities are given by
Using Equations (
5)–(
7), the basic Equations (
1)–(
4) are transformed into the perturbed dimensionless equations. For convenience, the asterisk is removed from
,
and
. The perturbed dimensionless equations are given by
where
The control physical parameters in Equations (
8)–(
10) are given by
R =
which accounts for buoyancy effects,
d relates to the porous medium,
, is the separation ratio due to the Soret effect and
is the ratio of mass diffusivity to thermal diffusivity. As
decreases, the thermal diffusivity increases and accordingly the mass diffusivity decreases. In the linear part of Equations (
8)–(
10),
and
C are eliminated to obtain the following form:
and
The horizontal planes
and
are assumed to be kept at constant boundary conditions and hence
and also the normal component of the velocity should vanish on the horizontal planes, i.e.,
The conditions (
12) and (
13) are applicable for any combination of top–bottom boundaries such as free–free or rigid–rigid or rigid–free or free–rigid. The stress-free boundary conditions [
30] are considered and are given by
Since the selected physical system is a double-diffusive system, it is unstable to either stationary convection or oscillatory convection at the onset.
6. Distortion of Streamlines and Isotherms
The fluid flow displayed using the stream function
is obtained from the
u and
w velocity components. For two-dimensional flow, the relationship between
and velocity components is given by [
34]
which yields the following single equation and by utilizing the accompanying limit conditions for
and removing the arbitrary integration constants (Kuo [
23]):
The family of streamlines gives a clear view of the entire flow domain and its important characteristics. Whereas the points with uniform temperature connected with lines that are called as the isotherms. The snapshots of the flow field and the hot regions close to the beginning of oscillatory convection are shown in terms of streamlines and isotherms, respectively, for distinct values of R, , , and t.
The combined heat and mass transfer phenomena are extended to the convective heat transfer principle of heatlines along with the graphical presentation of convective mass transfer. A representative set of streamlines, isotherms and isoconcentration lines is presented in
Figure 8,
Figure 9,
Figure 10,
Figure 11 and
Figure 12 for the variation in
R,
,
,
t and
. The circulation strength is assumed to be positive in the anticlockwise direction and negative in the clockwise direction. The snapshot of streamlines, isotherms and isoconcentration lines for increasing values of
R with
,
,
=
and
are displayed in
Figure 8a–l. The fluid flow is represented at the beginning by cellular patterns (
Figure 8a). The cell lying in
(
Figure 8a) has the absolute maximum and minimum circulation strength as
and
, respectively. From
Figure 8b, near the start of oscillatory convection, isotherms depict the effect on the temperature field with absolute maximum and minimum values
and
, respectively. Isotherms are observed to be nearly horizontal parallel lines, suggesting that conduction carries much of the heat transfer.
Figure 8c is also a plot of the isoconcentration lines that are calculated when
, and the absolute maximum and minimum values of these isoconcentration lines are
and
, respectively.
Figure 8d–f illustrate streamlines, isotherms and isoconcentration lines, respectively, for
for the same values of other physical parameters as those defined in
Figure 8a–c. The temperature, concentration and gravitational buoyancy forces act together and change the fluid flow pattern, as shown by the streamlines (
Figure 8d). The bicellular patterns become slightly deformed (
Figure 8d) and, for a cell lying between
, the absolute maximum and absolute minimum of circulation strength are
and
, respectively. The isotherms in
Figure 8e unambiguously indicate the presence of a convective mode of heat flow inside the fluid with absolute maximum and minimum values as
and
, respectively.
Figure 8f plots the isoconcentration lines, indicating the absolute maximum and minimum values as
and
, respectively.
Figure 8g, is plotted for
. As the thermal buoyancy increases, two small vortices, namely
B and
, near the top right and left bottom plates are generated with the equal circulation strength of
. The original cell becomes more deformed as
R changes from
to
(
Figure 8g). For a cell lying between
, the absolute maximum circulation strength is
and the minimum is
. The isotherms in
Figure 8h exhibit, inside the fluid layer, the convective mode of heat transfer creates the oval-shaped circular patterns in the middle region with absolute maximum and minimum values of
and
, respectively.
Figure 8i is a plot of the isoconcentration lines and their the absolute maximum and minimum values are
and
, respectively.
Figure 8j, is plotted for
. When
R is increased from
to
, due to the increased thermal buoyancy, the basic cell encounters more deformation, as shown in
Figure 8j. Additionally, the two small generated vortices,
B and
, are grown in size and shape with the increased circulation strength of
. Thus, the basic cell divided into two more new vortices, namely
A and
near the top and bottom of the basic cell with the same circulation strength of
. The isotherms that are shown in
Figure 8k demonstrate the strong convective heat flow in the fluid layer, with the absolute maximum and minimum values of
and
, respectively, which produces the circular oval patterns in the middle of the domain.
Figure 8l plots the isoconcentration lines, with the absolute maximum and minimum values of
and
, respectively. From
Figure 8a–l, it is noted that when
R is increased from
to
, oscillatory flow moves very quickly and the transfer of mass increases due to thermal buoyancy.
Figure 9 illustrates the streamlines, isotherms and isoconcentrations for distinct values of
with a fixed set of other parameters, say,
,
,
and
.
Figure 9a–c are plotted for
=
. In
Figure 9a the primary concentric roll patterns of streamlines are tilted towards the right side by retaining the flow orientation in the considered range of
. The absolute maximum and minimum values of circulation strength are
and
, respectively.
Figure 9d–f are plotted for
=
, with the same fixed set of parameters and contour lines defined in
Figure 9a–c. In
Figure 9d, the flow pattern in the range
becomes more deformed by creating the two additional vortices
B and
at the top and bottom limits, respectively, and on either side of the primary cell. The circulation strength of the two vortices,
B and
is
. In the area,
, circulation strength has absolute maximum and minimum values of
and
, respectively. From
Figure 9a,d, we can observe that the streamline deformation is increased when
values decrease. This implies that the effect of an increase in
values stabilizes the convective system. The flow of heat transfer is depicted using the isotherms and is shown in
Figure 9b,e for
=
and
, respectively. The isotherms appeared as plumes in the range of
, as shown in
Figure 9b. The boundary layer thickness increased when
decreased from
to
. The isoconcentrations (concentration field) are plotted in
Figure 9c,f for
2 and
, respectively. The isoconcentrations occurred in the form of vortices in the range of
. The size of vortex structures increased as
decreased. The boundary layer of the concentration region increases as
decreases from
to
. For a particular
value, the distribution of the concentration field is higher than the temperature field due to Soret effect.
The distributions of velocity, temperature and concentration fields are plotted in
Figure 10a–f for distinct values of
with
,
,
and
. In
Figure 10a, the streamlines are plotted for
. In the range of
the basic primary cell has been deformed and two new vortices, namely,
B and
, are generated at the top left and right bottoms, respectively, of the basic primary cell with a circulation strength of
. The absolute maximum circulation strength of the basic cell, namely A, is 0.036. In
Figure 10d, the streamlines are plotted for the decreased value of
i.e.,
=
, with the same parameters defined to plot the contour lines as shown in
Figure 10a. The generated two vortices, namely,
B and
in
Figure 10a become enlarged with the increased circulation strength of
, as shown in
Figure 10d. The primary cell is split into two vortices, namely,
A and
, with the absolute maximum circulation strength of the basic cell
. As
decreases, the thermal diffusivity increases and accordingly, the mass diffusivity decreases. As
decreases from
to
, the streamlines become more deformed. This implies that the decreasing
destabilizes the convective system.
The flow of heat transfer is depicted by using the isotherms and is shown in
Figure 10b,e for
=
and
, respectively. It is observed that the isotherms in the hot regions get enlarged as
decreases. Hence, we can conclude that, with the decreasing values
, the system becomes unstable.
Figure 10c is a plot of the isoconcentration lines for
=
and the absolute maximum and minimum values of isoconcentrations are
and
, respectively.
Figure 10f is a plot of the isoconcentration for
=
, and the absolute maximum and minimum values of isoconcentrations are
and
, respectively.
Figure 10c,f show
and
, respectively. The size of the vortex patterns decreases as
decreases. Finally, as
decreases, the thickness of thermal boundary layer increases, and the boundary layer thickness of the concentration field decreases. Hence, the increasing values of
stabilize the system.
Figure 11 illustrates the pattern of streamlines, isotherms and isoconcentration lines for distinct values of
t with
,
,
and
.
Figure 11a represents the pattern of streamlines for
and have the absolute maximum and absolute minimum values of circulation strength
and
, respectively, in the range of
. The streamlines are of parallel concentric roll patterns, thus indicating the stable fluid flow at
.
Figure 11d shows the pattern of streamlines for
t =
for the same values of other physical parameters that are considered in
Figure 11a. The absolute maximum and absolute minimum values of circulation strength are
and
, respectively. The stable primary cells have been deformed and two new vortices are generated inside the primary cell, namely,
B and
, with an absolute minimum value of circulation strength as
. The isotherms plotted in
Figure 11b,e are for the same physical parameters that are considered in
Figure 11a,d, respectively.
Figure 11b shows the isotherms that appeared in the form of plumes in the middle region of the temperature field. As
t increases from
to
, the size of the plume structure increases vertically. This result suggests that the heat transfer is enhanced as time increases for the fixed values of the physical parameters.
Figure 11c is a plot of the isoconcentration lines for
t =
and the absolute maximum and minimum values of isoconcentrations are
and
, respectively.
Figure 11f plots the isoconcentration lines for the same values of physical parameters as those defined in
Figure 11c and for
t =
. The absolute maximum and minimum values of isoconcentrations are
and
, respectively. As
t increases from
to
, the convective mass transport of the fluid increases. Finally, the boundary layer thickness of both the temperature and concentration fields decrease as time increases.
Figure 12 shows the streamline patterns, isotherms and isoconcentration lines for distinct fixed values of
with
,
,
and
.
Figure 12a, illustrates the streamlines for
=
. The primary basic cell has been deformed and is pushed into the middle of the region by the two generated vortices, namely,
B and
near the top and bottom boundaries, respectively, in the considered range of
with the circulation strength of
. In
Figure 12d, the streamlines are drawn for the same values of physical parameters as defined in
Figure 12a and
=
. In the considered range of
, the two vortices
B and
, shown in
Figure 12a, increased in size with the circulation strength of
(
Figure 12d). The basic cell is split into two new vortices, namely,
A and
, with a maximum circulation strength of
. This implies that the decreasing
destabilizes the convective system. The flow of heat transfer is depicted through isotherms and is shown in
Figure 12b,e for
=
and
, respectively. It can be observed that, for the decreased values of
, the hot regions increased. Hence, as
decreases, the system becomes unstable.
Figure 12c is a plot of the isoconcentration lines for
=
, the absolute maximum and minimum values of isoconcentrations are
and
, respectively. For
=
, in
Figure 12f, the plot of the isoconcentrations are shown with the absolute maximum and minimum values are
and
, respectively. This figure reveals that, as
decreases from
to
, mass transfer decreases due to less permeability of the porous medium.
Topology of Flow
Fluid flow topology is explored based on the Euler number,
. As described by Jana et al. [
35], on the surface,
is the aggregate of the Poincare indices of the critical points. The invariance of topology relation is defined by
where the number of elliptic, hyperbolic and parabolic points are denoted by
,
and
, respectively [
36,
37]. The vorticity contours for
,
= −0.01 and
= 0.06 are shown in
Figure 13a. For a fixed interval of
X, topological rule given in Equation (
49) is satisfied with
and
. A similar study has been carried out for vorticity contours with
depicted in
Figure 13b. For
, Equation (
49) is also satisfied by the vorticity contours as exhibited in
Figure 13b by
= 8,
= 8 and
= 0. The vorticity contours appeared in the form of vortices (elliptical shape) for
, and they are known as primary contours. Each pair of consecutive vortices has the counter-rotating form. The value of vorticity on a line is constant. The value of vorticity decreases as it moves away from the center of a primary cell, i.e., the circulation of fluid is less near the center. When
R increased from
to
, two new vortices showed their presence at the top and bottom walls and on either side of the primary vortex. Due to the presence of these new vortices, the primary vortex became deformed and, inside the primary vorticity contour, two new vortices emerged. Thus, the strength of vorticity and the velocity of the fluid increases as
R increases.
7. Heatfunction
The convective heat transfer phenomenon is represented by heatlines in the present study. Morega and Bejan [
38] successfully used the concept of heatlines to visualize the fluid flow. Various researchers [
39,
40,
41,
42,
43,
44,
45,
46] enlarged this concept to a wide range applications of natural convective systems. The nondimensional heatfunction,
is defined as
where
and
denotes the static temperature. Here
denotes the arbitrary reference temperature. It should be noticed that different values of
lead to different slopes of
[
40,
41]. The
defined in Equations (
50) and (
51) is similar to the one defined by Kimura and Bejan [
42] but with opposite signs. Eliminating the temperature gradients from the above Equations (
50) and (
51) by cross differentiation, the following Poisson-type equation is obtained:
The following boundary conditions for
are obtained by integrating the Equations (
50) and (
51):
The integral expressions that exist in the boundary conditions (
53)–(
56) represent the average Nusselt numbers on the corresponding boundaries. The total heat transfer across the system in terms of
is plotted graphically in
Figure 14,
Figure 15,
Figure 16,
Figure 17 and
Figure 18.
Figure 14a–d show the pattern of heatlines for
,
,
and
respectively, for
,
,
and
. From
Figure 14a, the heatline contours within the domain are found to be normal to
and
lines due to conduction dominant heat transfer. The absolute maximum and absolute minimum values of heatlines are
and
, respectively, in the considered range
. In the neighborhood of the
line, the heatlines are vertically parallel. The curvature at the central part of each heatline increases with
X and the heatlines are similar to that of parabolic structure. It represents the nonlinear propagation of heat transfer that occurs at
. Hence the transition takes place from the conduction to the convection state at
.
Figure 14b shows the heatlines for
and the absolute maximum and minimum values of heatlines are
and
, respectively. The heatlines shown in
Figure 14a are converted to curves and closed paths as depicted in
Figure 14b when
increased to
.
Figure 14c illustrates the heatlines for
and have the absolute maximum and minimum
and
, respectively.
Figure 14d shows the heatlines for
, which have the absolute maximum and minimum
and
, respectively. This figure shows that some of the heatlines shown in
Figure 14c are divided into closed paths (vortices) and curves. It can also be noted from these figures that the size of closed paths increases as
R increases. Thus, from
Figure 14a–d, we can conclude that the heat transfer intensity increases with
R.
Looking at
Figure 15a,b, the influence of
is analyzed on the heat transfer for fixed values of other physical parameters,
,
,
and
. For
, the heatlines drawn in
Figure 15a, have the absolute maximum and the absolute minimum values
and
, respectively, in the assumed range
.
Figure 15b shows the heatlines for
and the absolute maximum and absolute minimum values are
and
, respectively.
Figure 15a shows the pattern of heatlines, namely, the combination of curves and vortices with vortices located in the central part of the considered domain.
Figure 15b shows that the deformation of the curved lines and the size of the vortices are reduced compared to those lines shown in
Figure 15a. Thus,
Figure 15a,b show the decreasing intensity of heat transfer as
increases from −0.02 to −0.01.
Figure 16a,b illustrate the visualization of heat flow via heatlines for different
and for the fixed values of other physical parameters,
,
,
and
.
Figure 16a illustrates heatlines for
and has the absolute maximum and the absolute minimum
and
, respectively.
Figure 16b shows heatlines for
and have the absolute maximum and the absolute minimum
and
, respectively, in the considered range of
X. The heat flow patterns of
Figure 16a,b are similar to those in
Figure 15a,b, respectively. The deformation of the curved lines and the size of vortices decreases as
increases from 0.05 to 0.06. Thus, as
increases the thermal diffusivity decreases.
The snapshot of heatlines for fixed
,
,
and
for different
t are shown in
Figure 17a,b.
Figure 17a shows the pattern of heatlines for
, which have a structure similar to that of heatlines for
with the absolute maximum and absolute minimum being
and
, respectively.
Figure 17b illustrates the heatlines for
. The absolute maximum and the absolute minimum are
and
, respectively, in the considered range of
X. Some heatlines shown in
Figure 17a are converted into curves and closed paths as
t increases. Thus, as
t increases the intensity of heat flow also increases.
Figure 18a,b show the visualization of heat flow via heatlines for different values of
and for fixed values of other parameters
,
,
and
.
Figure 18a illustrates heatlines for
and have the absolute maximum and absolute minimum
and
, respectively.
Figure 18b illustrates heatlines for
with the absolute maximum and absolute minimum
and
, respectively, in the considered range
.
Figure 18a,b show that, as
decreases, the intensity of the heat flow increases.
8. Entropy Generation
It is important to understand the efficient use of energy resources, together with the minimum degradation of energy, i.e., the minimum generation of entropy should be efficient during heat transfer. In the present work, entropy generation is related to the irreversible existence of heat transfer and viscous effects within the fluid and the fluid-solid interfaces. The non-dimensional form of the local entropy generation rate (
) is calculated by the addition of the heat transfer irreversibility (
) and the fluid friction irreversibility (
) as discussed in [
47,
48,
49] and is expressed as:
where
w and
T are the known non-dimensional velocity and temperature functions given in
Section 2 and
) is the irreversibility distribution ratio, which is initially minimized numerically (fixed as
) and the physical parameter
was fixed as
. Al-Hadhrami et al. [
50] proposed three different models, namely, Darcy, Darcy–Freshman and Brinkman–Brinkman, for viscous dissipation flow through porous media. Hooman and Gurgenci [
51] compared these three models. They concluded that these models were effectively showed the same results only for small values of
. The viscous dissipation model proposed by Baytas [
52] is employed to obtain Equation (
57).
Figure 19a–f show the results of entropy generation for different
R. At
, the pattern of
is the combination of vortices and curves. The vortices are converted into the wavy pattern as
R increases. At
, the pattern of
appeared in the form of vortices for
(
Figure 19b). For
, the above pattern is converted into the combination of elliptical and deformed vortices. Each deformed vortex encloses five other vortices (
Figure 19d). When
R is increased to
, the patterns become more deformed, as shown in
Figure 19f. It is observed that, for
(
Figure 19a,b), the maximum entropy generations due to
and
are 0.03098 and 0.00032, respectively. The maximum entropy generation due to
is 3.69005 when
and 18.6827 when
(
Figure 19c,e). Similarly, the maximum entropy generation due to
is 0.5214 when
and 4.86893 when
, i.e., for small values of
R, the
and
are small and they increase with
R. This is because, at low
R, the viscous force dominates and the heat transfer is mainly due to conduction. As already mentioned, the
is the sum of the two contributions, namely,
and
. Thus, it can be concluded that
increases as
R increases.