Next Article in Journal
Spinel CoFe2O4 Nanoflakes: A Path to Enhance Energy Generation and Environmental Remediation Potential of Waste-Derived rGO
Next Article in Special Issue
Convective Heat Transfer in Magneto-Hydrodynamic Carreau Fluid with Temperature Dependent Viscosity and Thermal Conductivity
Previous Article in Journal
Low-Power Complementary Inverter Based on Graphene/Carbon-Nanotube and Graphene/MoS2 Barristors
Previous Article in Special Issue
The Impacts of Viscoelastic Behavior on Electrokinetic Energy Conversion for Jeffreys Fluid in Microtubes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Thermal Performance Analysis of an Al2O3-Water Nanofluid Flow in a Composite Microchannel

Centre for Advanced Mechanical and Green Technology, Faculty of Engineering and Technology, Multimedia University, Bukit Beruang, Melaka 75450, Malaysia
*
Author to whom correspondence should be addressed.
Nanomaterials 2022, 12(21), 3821; https://doi.org/10.3390/nano12213821
Submission received: 8 October 2022 / Revised: 16 October 2022 / Accepted: 20 October 2022 / Published: 28 October 2022
(This article belongs to the Special Issue Advances of Nanoscale Fluid Mechanics)

Abstract

:
Partial filling of porous medium insert in a channel alleviates the tremendous pressure drop associated with a porous medium saturated channel, and enhances heat transfer at an optimum fraction of porous medium filling. This study pioneered an investigation into the viscous dissipative forced convective heat transfer in a parallel-plate channel, partially occupied with a porous medium at the core, under local thermal non-equilibrium condition. Solving the thermal energy equation along the Darcy–Brinkman equation, new exact temperature fields and Nusselt number are presented under symmetrical isoflux thermal boundary condition. Noteworthy is the heat flux bifurcation at the interface between the clear fluid and porous medium driven by viscous dissipation, in cases where the combined hydrodynamic resistance to fluid flow and thermal resistance to fluid conduction is considerable in low Darcy number porous medium insert. However, viscous dissipation does not affect the qualitative variation of the Nusselt number with the fraction of porous medium filling. By using Al2O3-Water nanofluid as the working fluid in a uniformly heated microchannel, partially filled with an optimum volume fraction of porous medium, the heat transfer coefficient improves as compared to utilizing water. The accompanied viscous dissipation however has a more adverse impact on the heat transfer coefficient of nanofluids with an increasing Reynolds number.

1. Introduction

Forced convection heat transfer in porous medium is encountered in broad engineering applications, such as heat exchangers, oil recovery, electronics cooling, heat pipes, catalytic reactors, fuel cells, thermal insulation, geothermal systems, and many more [1]. To cut down the tremendous pressure drop required to drive a flow through a channel saturated with porous medium, partial filling of porous medium becomes an alternative to augment convection heat transfer while maintaining the pumping cost at an acceptable level [2,3,4].
Studies on a channel partially filled with porous medium are less common than ones which are fully filled. Past studies in references [3,5,6,7,8,9,10,11,12,13,14,15] developed the velocity and temperature distributions based on the local thermal equilibrium (LTE) assumption in a channel partially filled with porous medium. The interests span various configurations of porous medium insert [3,5,6,7,8,9,10,11,12,13,14,15], laminar and turbulent flow [6,7], thermal dispersion effect [9], effect of moving boundary condition [10], and pulsating flow [15]. The fluid flow and heat transfer at the interface region was analyzed for three different cases, namely two different porous medium, a fluid region, and a porous medium, and lastly, an impermeable medium and a porous medium [5]. The study derived velocity and temperature distribution for all three interface conditions, and found great agreement between their analytical and numerical solutions, by using the empirical hypothesis [8]. The non-Darcy effects were numerically studied on the fully developed forced convection in a parallel plate channel partially filled with a porous medium, by assuming Darcy–Brinkman–Forchheimer model in a porous medium [9]. The study indicated that convective heat transfer enhances the thermal dispersion effect, more significantly at smaller values of thickness ratio for the fluid layer to porous medium layer. An investigation on a Couette flow in a composite channel partially filled with a porous medium, was carried out by adopting the Brinkman–Forchheimer–extended Darcy equation, for the porous medium region [10]. The study obtained the velocity and temperature distributions by using the stress jump condition at the interface region [11] and showed that Nusselt number depends on the size of the gap between a porous medium and moving plate; Nusselt number increases with a decreasing gap size and Darcy number D a . A partially filled channel with porous medium placed at the core, imposed with symmetrical boundary heat flux condition, was analyzed both analytically and numerically in another study [12]. The use of porous layer for heat transfer enhancement and the accompanying pressure drop was evaluated. The study showed that the overall performance is more strongly affected by D a , but not so much by the thermal conductivity ratio of porous medium to fluid layer, as the ratio increases beyond two. Subsequently, a theoretical analysis was conducted on the heat transfer and fluid flow of a channel, with a porous medium inserted adjacent to the wall, assuming uniform heat flux thermal boundary condition [13]. The study perceived a good agreement between the analytical and numerical results, and showed that the maximum overall Nusselt number based on the definition of mean thermal conductivity of porous medium and clear fluid is attainable for a thermal conductivity ratio of one. Notably, the location of porous medium insert in a circular pipe influenced heat transfer rates substantially [14]. A study on the heat transfer for a pulsating flow in a circular pipe partially filled with porous medium was tackled numerically [15]. The study analyzed the effect of different parameters on convective heat transfer and determined the optimal thickness of the porous layer.
LTE assumption in the aforementioned studies is valid when there is a negligible temperature difference between solid and fluid phases. Local thermal non-equilibrium (LTNE) or two-equation model in a porous medium reflects the fluid and solid temperatures more accurately. However, investigations on the enhancement of heat transfer in a channel partially filled with porous medium under LTNE conditions remain relatively scarce. One area of interest lies in the modelling of temperature and heat flux at the interface between the porous medium and clear fluid. For a uniform heat flux boundary condition, two approaches were proposed to model the thermal boundary condition at the interface [16]. In the first approach, depending upon the effective conductivities and temperature gradient of solid and fluid phases, the total heat flux is divided between two phases, whereas an equal amount of heat flux received by individual phases at the boundary is modelled in the second approach. In another study, a jump condition was introduced at the interface between a porous medium and a homogenous fluid [17]. The study used the volume-average technique for thermal energy equations in solid and fluid phases and showed that a jump condition contains an excess in heat flux at an interface. A systematic analysis was presented to model different interfacial momentum and thermal boundary conditions between a porous medium and a fluid layer [18]. Variances among different models were compared and reviewed, revealing that for most of the practical applications, the variance in different models has the least impact on temperature and Nusselt number distributions. Another study adopted the two aforementioned approaches [16] to model heat transfer in porous medium, and termed the varied heat flux, and uniform heat flux interface models as model A, and model B respectively, in the channels partly occupied by porous medium [19]. Under LTNE condition, and with model A applied at the interface, the study obtained the exact temperature field for a parallel plate partially occupied by porous insert [4]. Nusselt number is affected by three conditions, namely channelling effect, fluid conduction from the interface, and solid conduction and internal heat exchange in a porous medium. The study also highlighted that an increased thickness of porous medium insert amplified the error in an LTE model and suggested the optimum fraction of porous material to enhance the heat transfer to be 0.8 under reasonable pressure drop. Amid growing interest to examine the validity of interface boundary conditions between two regions, a study reported the restrictions on the validity of thermal boundary conditions at the interface and presented the largest Nusselt number for model A [20]. The study further explored the heat flux bifurcation inside the composite channel by considering the inertial and thermal dispersion effects, and applying three different interfacial thermal conditions [21]. The study discerned the heat flux bifurcation phenomenon at the porous-fluid interface when the temperatures of the solid and fluid phases are not equal, while the inertial terms show weaker effects when the thermal dispersion effect is included and vice versa. The exact solutions were obtained for both fluid and solid phases by inserting the metallic foams on both sides of the tube under the LTNE model [22]. In order to predict the effectiveness of inserting porous material into the tube, an analytical investigation was performed on forced convection in a heated tube with a porous medium core, and a tube with a porous medium insert at the side, adopting model A as a thermal boundary condition under the LTE, and LTNE conditions [23]. The study took into account of thermal dispersion effects and found that the heat transfer performance of the tube with a porous medium core is higher in a lower range of pumping power than that of the tube with porous medium layer at the side, while the latter’s performance exceeds the former in a higher range of pumping power. More recently, the variations between the resulting temperature distribution and the Nusselt number for model A and model B were analyzed [24,25]. An analytical solution was developed to predict the validity of the LTE model [24]. Although a higher Nusselt number was obtained for model A, the heat transfer performance depends less on the interface model and D a when the volume fraction of porous medium is less than 0.9. An analytical study on a partially filled channel under LTNE condition, subjected to a constant heat flux, with the porous insert attached to both sides of channel walls [26], was subsequently looked into. The study adopted the stress jump and continuity condition at the interface and examined the effects of stress jump coefficient on the fluid flow and heat transfer using model A. It predicted an optimum volume fraction of 0.2 for clear fluid region in order to enhance the heat transfer. The Nusselt number depends more strongly on the thermal conductivity ratio than the porous medium fraction. The studies [4,22,23,26,27,28,29,30] showed that the more commonly used model A leads to a higher Nusselt number.
Viscous dissipation is a heat source owing to work done by viscous stresses. The effects of viscous dissipation were neglected in most analyses of forced convection in a channel partially filled with porous medium, exposed to LTNE conditions [4,20,23,24,26,27,31]. Studies [29,32] show that the internal heat source alters the temperature, Nusselt number, and entropy generation of the system significantly in either the solid or fluid phases in porous medium. With viscous dissipation incorporated in the energy equations, the mathematical model examined the heat transfer and entropy generation for a partially filled channel, subjected to constant heat flux at one wall, and an adiabatic condition on the other wall using the LTNE model [33]. The study reported the bifurcation of heat transfer and entropy generation as well as the possibility of having a temperature higher than that of the heated wall. A new analytical solution was next presented for a composite channel partly occupied by a porous material next to the wall [34]. The work highlighted heat flux bifurcation near the interface region and a more dominant viscous dissipation effect in the clear fluid region, intensified with an increasing D a . In another study, internal heat generation effects, within the fluid and solid phases of a porous medium under LTNE condition, for models A and B were analysed for a partially filled channel, along with the phenomenon of heat flux bifurcation at the porous-fluid interface [35]. The first and second law analyses were undertaken in a microreactor partially filled with porous medium insert, with internal heat generation [36]. Imposing asymmetrical temperatures at the boundary surface and assuming a thick outer wall, internal heat generations was found to significantly affect the system. In another study, three viscous dissipation models were studied, namely, the Darcy model, the form drag model, and the clear fluid compatible model, for a channel incorporated with a porous insert next to the wall [37]. The study presented an analytical solution, without involving the convection term, for the energy equation. The temperatures and heat transfer coefficients are affected more significantly by viscous dissipation with an increasing Da. The best heat transfer enhancement is achieved in the fluid-compatible viscous dissipation model. For both LTE and LTNE conditions, a recent study assessed the effect of internal heat sources on the thermal performance and reported that an increasing thickness of porous insert causes singularities in Nusselt number [38]. The study is however restricted to uniform heat sources in fluid and solid phases.
To this end, existing literature points to the fact that the LTNE model is able to present a more comprehensive temperature field, while viscous dissipation influences the temperature field more significantly in cases of large D a besides imposing heat flux bifurcations at the interface in a channel partially filled with porous medium adjacent to the wall. The configuration of the composite wall has also been highlighted to have an effect on the thermal performance. To the best of the authors’ knowledge, the effects of viscous dissipation on the thermal aspects of a composite channel, with a porous insert at its core, has not been investigated under LTNE condition despite its rather common configuration in a composite channel. This study hence intends to bridge the gap by focusing on the arising heat transfer characteristics manifested by viscous dissipation in the composite channel, uniformly heated at the boundaries. The results are subsequently applied to a composite microchannel, with Al2O3-Water nanofluid as the coolant.

2. Problem Formulation

This problem is formulated for forced convective heat transfer in a parallel plate channel, partially filled with porous medium at the core, as shown in Figure 1. The plates are separated by a distance 2 H , with a porous medium of thickness 2 h inserted at the core of the channel, subjected to isoflux thermal boundary condition. The flow is assumed to be hydrodynamically fully developed, steady, and laminar within the channel, while LTNE is assumed in the porous medium. Due to the symmetry of the problem, only half of the channel is solved in this study.

2.1. Problem Formulation

The velocity profiles presented in [24] are adopted in this study, obtained by solving the Navier Stokes equation and Brinkman-extended Darcy equation, in the clear fluid and porous medium, respectively. Interface continuity is assumed for the velocity and shear stress. The velocity distribution [24] is given as
U c = Y 2 2 + A 1 Y + A 2
U p = A 3 cosh A 0 Y + D a  
where A 0 A 3 are the constants listed in Appendix A.

2.2. Problem Formulation

The governing thermal energy equations for the respective clear fluid, and LTNE porous medium, are [34]:
Clear fluid:
k f 2 T c y 2 + μ f d u c d y 2 = ρ f c p , f u c T c   x      
Porous medium:
k s , e f f 2 T s y 2 h i a i T s T f = 0
k f , e f f 2 T f y 2 + h i a i T s T f + μ f K u p 2 + μ e f f d u p d y 2 = ρ f c p , f u p T f   x
where in the thermal energy for fluid in porous medium, Equation (5), a clear fluid compatible model [39] for viscous dissipation is incorporated.
In this study, the temperature at the interface between the clear fluid region and the porous medium is assumed to be constant. Based on the difference in effective thermal conductivities and temperature gradients between the two phases in a porous medium, heat flux is divided between clear fluid and porous medium [16].
The corresponding thermal boundary conditions are [16]
y = 0 , T f y = 0 , T s y = 0
y = h ,             T f = T s = T c ,             k f , e f f   T f y + k s , e f f   T s y = k f   T c y = q i n t
y = H ,             k f T c y = q w
For a fully developed condition subjected to constant heat flux, T x = T f x = T c x = T m x = constant   , and is computed from
T x = 1 ρ f c p , f u m H q w + μ f K 0 h u p 2 d y + μ f ϵ 0 h d u p d y 2 d y + μ f h H d u c d y 2 d y          
where u m = 1 H 0 h u p d y + h H u c d y .
By introducing
θ = k s , e f f T T i n t q w H ,   B i = h i a i H 2 k s , e f f ,   κ = k f , e f f k s , e f f , B r = μ f u m 2 H q w K , γ = q i n t q w , ϵ = k f , e f f k f
The resulting non-dimensional equations are
κ ϵ 2 θ c Y 2 = U c ω U m B r D a U m 2 d U c d Y 2
κ 2 θ f Y 2 + B i θ s θ f = U p ω U m B r U m 2 U p 2 B r D a U m 2 ϵ d U p d Y 2
2 θ s Y 2 B i θ s θ f = 0
where
ω = 1 + B r U m 2 0 η U p 2 d Y + B r D a U m 2 ϵ 0 η d U p d Y 2 d Y + B r D a U m 2 η 1 d U c d y 2 d Y
The corresponding thermal boundary conditions are
Y = 0 ,             θ f Y = 0 ,             θ s Y = 0
Y = η ,             θ f = θ s = θ c = 0 ,             κ ϵ θ c Y = κ θ f Y + θ s Y = γ
Y = 1 ,             κ ϵ θ c Y = 1 .  
The resulting Equations (18) and (19) are obtained are after some manipulations to decouple Equations (12) and (13).
4 θ f Y 4 B i 1 + κ κ 2 θ f Y 2 = U p ω B i κ U m ω κ U m d 2 U p d Y 2 B r B i κ U m 2 U p 2 + D a ϵ d U p d Y 2
+ B r κ U m 2 d 2 U p 2 d Y 2 + D a ϵ   d 2 ( d U p / d Y ) 2 d Y 2
4 θ s Y 4 B i 1 + κ κ 2 θ s Y 2 = U p ω B i κ U m B r B i κ U m 2 U p 2 + D a ϵ d U p d Y 2
To solve Equations (18) and (19), additional thermal boundary conditions are derived from Equations (12) an (13), giving
Y = 0 ,             3 θ f Y 3 = 0
Y = 0 ,             3 θ s Y 3 = 0
Y = η ,             2 θ f Y 2 = U p ω U m κ B r U p 2 U m 2 κ B r D a U m 2 ϵ κ d U p d Y 2 = φ
Y = η ,             2 θ s Y 2 = 0
where φ is a coefficient defined in Appendix A.

2.3. Analytical Solution

2.3.1. Temperature

The analytical solution for the temperatures in clear fluid and porous medium are derived as follows:
θ c = D 1 Y + D 2 Y 2 + D 3 Y 3 + D 4 Y 4 + D 5
θ f = E 1 Y 2 + E 2 cosh A 0 Y + E 3 cosh C 0 Y + E 4 cosh 2 A 0 Y + E 5
θ s = E 1 Y 2 + F 1 cosh A 0 Y + F 2 cosh C 0 Y + F 3 cosh 2 A 0 Y + F 4

2.3.2. Dimensional Temperature Variation with x and y

In the respective porous medium, and clear fluid region, the dimensional temperatures derived from Equations (24) and (25), are recast as
T c = θ c q w H k s , e f f + T c x x + T i n
T f = θ f q w H k s , e f f + T f x x + T i n

2.4. Problem Formulation

Nusselt number is defined as
N u = q w H k f T w T m
can be recast as
N u = 2 ϵ   κ ( θ w θ m , f )
where the bulk mean fluid temperature is obtained from
θ m , f = 1 U m 0 η U p θ f d Y + η 1 U c θ c d Y
The new Nusselt number expression is hence obtained as
N u = G 1 + 420 2 A 2 D 5 + A 3 E 2 + 2 D a E 5 + 12 D a A 0 3 E 5 η + A 3 E 2 + 2 D a E 4
3 sinh 2 A 0 η + 2 A 3 E 4 sinh 3 A 0 η 420 A 2 D 1 + A 1 D 5 η 2 140   ( 2 A 1 D 1 + 2 A 2 D 2
D 5 2 D a E 1 + 4 D a A 0 3 E 1 ) η 3 + 105 D 1 2 A 1 D 2 + A 2 D 3 η 4 + 84 [ D 2 2 ( A 1 D 3
+ A 2 D 4 ) ] η 5 + 70 D 3 2 A 1 D 4 η 6 + 60 D 4 A 3 E 2 + 2 D a E 4 η 7
where the coefficient G 1 is given in Appendix A.

2.5. Verification of Results

Nusselt number is compared with the literature [40,41], reduced to cases completely filled with clear fluid η = 0 , and fully filled with porous medium ( η = 1 ). In the porous medium, as κ , B i 0 , and D a ,   a porous medium approaches a clear fluid [40]. Table 1 shows that the results in this study concur excellently with the literature.

3. Results and Discussion

The optimum volume fraction of porous medium η used is dependent on the parameters κ, B i , and D a . For brevity, we fixed the thickness of the porous medium to be η = 0.9 , for all subsections except for Section 3.4, which used the optimum volume fraction η = 0.45 , in a microchannel.

3.1. Velcoity Profile

Velocity profile presented in [24] is graphed to facilitate the discussion in the ensuing sections. Figure 2 represents the maximum velocity obtained across the channel for various D a and porous medium thickness. It depicts the highest velocity achievable for the respective D a , at a specific thickness of porous medium. Beyond this thickness, velocity drops as porous medium insert fills up the channel.

3.2. Temperature Profile

Figure 3a–d depict the dimensionless temperature at specified B i and κ for D a = 10 4 , in a channel partially filled with porous medium up to η = 0.9 , with B r as a parameter. An increasing B r amplifies the temperature difference between solid and fluid in the porous region, as well as the clear fluid temperature. In Figure 3a, high κ, and low B i render a large temperature difference between fluid and solid phases, indicating weak interfacial heat transfer between two phases. Fluid and solid temperatures are lowered as B r increases. For low κ and low B i , the same trend is observed in Figure 3c, accompanied by a larger temperature drop due to viscous dissipation. When there is strong resistance to fluid conduction in porous medium, internal heat exchange from a higher solid temperature to a lower fluid temperature becomes significant; a small B i gives rise to a large temperature difference. For a low κ and a big B i , depicted in Figure 3d, heat flux bifurcation is observed near the interface with the oppositely signed temperature gradients in the fluid and solid phases. The phenomenon magnifies inside the porous medium as the viscous dissipation increases. Heat flux bifurcation is reinforced by dominant solid conduction when κ is small, and B i is large.
Table 2 delineates the role of viscous dissipation in heat flux bifurcation at the interface region at fixed B r , for a set of stipulated parameters. It shows that increasing B r and B i intensifies the phenomenon of heat flux bifurcation at the interface, while a decrease in D a further magnifies this phenomenon at the interface. As D a decreases, a more intense viscous dissipation resulting from a large velocity gradient in the interface region intensifies the heat flux bifurcation; manifested as the temperature fluctuation at the interface region. It is noteworthy that when the permeability of the porous medium deviates more from that of clear fluid, hence exerting larger resistances to the fluid and heat flow, heat flux bifurcation intensifies despite the fact that a composite channel with a lower D a gives a better overall convection heat transfer enhancement.

3.3. Nusselt Number

Figure 4 shows the N u variation in the channel with the porous medium thickness, having D a and B r as parameters. The results are compared to those without viscous dissipation in Ref. [24].
In Figure 4a, for κ and B i both specified as 0.1, N u increases with an increase in the thickness of the porous medium and reaches a maximum at the optimal porous medium thickness. Beyond the optimal thickness, N u drops more significantly as D a decreases. The surge in N u is obvious as the thickness of porous medium constrained by the optimal thickness increases. The trend is not unexpected in light of the maximum velocity plot in Figure 2, whereby the optimal thickness corresponds to one that gives the highest achievable velocity. N u decreases sharply toward a more packed porous medium due to a rather drastic drop in velocity, which ultimately reduces the heat transfer. The trend is alike in the presence of viscous dissipation, except for a reduced N u . It is worth mentioning that the viscous dissipation effect is most prominent at the optimal porous medium thickness. Figure 4b depicts the N u variation for κ = 0.1 and B i = 10 . The trend of variation for D a = 10 3 is different in this case, with the maximum N u achieved for a fully filled channel ( η = 1 ) . When κ is small, and B i is large for a larger D a , the convection effect in a clear fluid is not sufficient to dissipate heat from the wall, hence extending the diffusion of heat to porous medium substantially. Under circumstances where the thermal resistance to fluid conduction is much larger than the thermal resistance to solid conduction, heat transfer in porous medium is dominated by a different transport mechanism as compared to that in Figure 4a, giving rise to an intensified heat flux bifurcation like what is shown in Figure 4d as heat tends to diffuse through the solid in porous medium. Nonetheless, viscous dissipation effects tend to resemble the velocity variation in Figure 2 for smaller D a , due to enhanced convection.
Likewise, for Figure 4c and d, as κ increases, the same trend is obtained, as discussed in Figure 4a. Figure 4c,d depict that increasing κ lowers the N u , because less amount of heat flux is transferred to the solid phase in the porous region from the interface. The trend, however, is on the contrary to the results observed in [33], which demonstrated asymptotes in such a plot when the top wall is insulated while the bottom wall is heated when porous medium insert is placed next to a heated wall.
Figure 5 depicts the variation of N u against B r with D a as the parameter. Figure 5a,b are plotted for small κ values, with B i fixed at 0.1 and 10, respectively, and Figure 5c,d are plotted for higher κ values for the same B i . Figure 5a–d show the same pattern of variation whereby N u decreases with increasing B r for all D a more significantly for larger D a .   N u decreases sharply toward a larger D a and an increasing B r due to a rather drastic hike in the velocity gradient, caused by viscous effects, which ultimately reduces the heat transfer. The extent of B r effect commensurate with the velocity surge at a specified porous medium filling as shown in Figure 2. Figure 5 also depicted that κ and B i have a more obvious effect on N u when these combined effects impose significant resistance to fluid conduction in porous region. Otherwise, N u hinges largely on the convection effect inside the clear fluid only region.

3.4. Microchannel Partially Filled with Porous Medium

The theoretical results in this study are applied to capture the effects of viscous dissipation in a microchannel partially filled with silicon-based porous medium saturated with Al2O3-Water nanofluid, as depicted in Figure 6. The specifications of the microchannel were taken from [42] and stated in Table 3.
The thermophysical properties of Al2O3 nanoparticles, obtained from the correlations in [43] with a d n p of 40 nm are given in Table 4.

3.4.1. Heat Transfer Coefficient

The heat transfer coefficient for 0% to 4% volume fraction of Al2O3-Water nanofluid at R e = 25 and 50 are presented in Figure 7a,b respectively with and without viscous dissipation. Figure 7 shows an increasing heat transfer coefficient with the volume fraction of nanoparticles. For this particular specification, the highest heat transfer coefficient is obtained when the microchannel is partially filled with a porous medium, η = 0.45 due to the synergistic effect of a more uniform velocity distribution and temperature profile for D a = 0.1 . The heat transfer coefficients for water are lower than ones for nanofluid and are barely affected by viscous dissipation. As Re increases to R e = 50 in Figure 7b, marked effect of viscous dissipation is notable where the heat transfer coefficient is lowered significantly because of the overriding effect of the heat source over heat convection.

3.4.2. Temperature Field in a Microchannel

Two-dimensional temperature line plots are produced for a fluid phase in porous material and clear fluid region in the microchannel based on Equations (27) and (28) as depicted in Figure 8. Figure 8 is plotted for the same D a with a porous insert fraction, η = 0.45 , and the wall heat flux is specified as 1 × 10 4 W/m2, the isotherms are calculated for a particular Reynolds number ( R e ) , from which the mean flow velocity and B r are determined. A temperature of 300 K is assumed at X = 0 , the leading edge of lower wall. Figure 8a,b depicted the isotherms for 0% volume fraction of nanoparticles, whereas those of 4% volume fraction of nanoparticles were presented in Figure 8c,d for R e = 25 , indicating negligible influence of viscous dissipation. On the other hand, having Al2O3-Water nanofluid in Figure 8c,d, for the same R e , the fluid temperature elevates along the microchannel as the viscous dissipation effects are considered. Hence, nanofluid increases the temperature along the flow as viscous dissipation effect is considered. It is also noteworthy that the suspension of nanoparticles gives rise to a lower temperature when viscous dissipation is assumed negligible.
Table 5 compares the heat transfer coefficient for 0%, and 4% volume fraction of Al2O3-Water nanofluid at fixed R e for the partially filled porous medium, with η = 0.45 . Viscous dissipation causes a more significant deterioration in the heat transfer coefficient of the 4% volume fraction of Al2O3-Water nanofluid as compared to water at R e = 25 , and with an increased R e , R e = 50 , the heat transfer coefficient diminished by 6% in a composite microchannel having the 4% Al2O3-Water as coolant. Nanofluid is a better coolant than water but viscous dissipation is deemed to cause a more adverse effect on heat convection in a nanofluid flow as the R e and the corresponding Br go up in a microchannel.

4. Conclusions

New exact Nusselt numbers and temperature profiles in a channel partially filled with porous medium at the core are presented for a forced convective heat transfer, with viscous dissipation accounted for under local thermal non-equilibrium conditions. Viscous dissipation has the most significant effect on the temperature at the interface in cases where there is strong thermal resistance to fluid diffusion in porous material combined with strong hydrodynamics resistance to fluid flow, causing heat flux bifurcation at the interface. Viscous dissipation notably causes a more significant temperature decrease across the channel with an increasing Da and hence lowers the Nusselt number. The extent of Br effect commensurates with the highest velocity achievable at a specified porous medium filling, noting that convection in a clear fluid is not sufficient to dissipate the wall heat flux. The heat transfer performance deteriorates when the viscous dissipation effect is accounted for, more so when Da approaches 1, in line with the Nu trend. The two-dimensional isotherm plots in the microchannel highlighted the significance of viscous dissipation even for low Reynolds number flow. Maximum heat transfer coefficient is obtained with a certain fraction of porous medium in the channel. The 4% vol. fraction Al2O3-water is a better coolant than water in the microchannel in cases where the Reynolds number is low. It is noteworthy that viscous dissipation however has a more notable adverse impact on the heat transfer coefficient of the nanofluid, more so as Reynolds number increases. This study contributes to the heat transfer coefficient computation in a channel partially filled with a porous medium at the core, which is useful to forced convection in a partially filled microchannel where viscous dissipation is significant.

Author Contributions

Conceptualization, G.M.C. and M.F.B.; methodology, G.M.C. and M.F.B.; validation, G.M.C. and M.F.B.; formal analysis, G.M.C. and M.F.B.; investigation, M.F.B.; resources, G.M.C.; data curation, M.F.B.; writing—original draft preparation, M.F.B.; writing—review and editing, G.M.C., M.F.B. and C.P.T.; visualization, M.F.B.; supervision, M.F.B. and C.P.T.; project administration, G.M.C.; funding acquisition, G.M.C. All authors have read and agreed to the published version of the manuscript.

Funding

The authors would like to express their utmost gratitude to the Ministry of Higher Education, Malaysia for the financial support provided by the FRGS grant: FRGS/1/2019/TK03/MMU/02/1.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

φ = ω U m D a + A 3 cosh A 0 η B r D a + A 3 cosh A 0 η 2 A 0 2 A 3 2 B r D a sinh A 0 η 2 / ϵ κ U m 2
A 0 = ϵ / D a
A 1 = 2 D a η cosh A 0 η 1 + 2 D a + η 2 ϵ 3 / 2 sin h A 0 η 2 D a η cosh A 0 η 2 1 + η ϵ 3 / 2 sin h A 0 η
A 2 = D a 1 2 η cosh A 0 η + 2 D a + 1 + η η ϵ 3 / 2 sin h A 0 η 2 D a η cosh A 0 η 2 1 + η ϵ 3 / 2 sin h A 0 η
A 3 = D a 2 D a + 1 + η 2 2 D a cosh A 0 η 2 1 + η ϵ 3 / 2 sin h A 0 η
C 0 = B i 1 + κ κ
D 1 = 2 1 3 A 1 + 3 A 1 2 B r D a + U m ω 3 A 1 ω 6 A 2 ω + 6 U m ϵ 6 κ U m 2
D 2 = A 1 2 B r D a + A 2 ω U m ϵ 2 κ U m 2
D 3 = A 1 2 B r D a + ω U m ϵ 6 κ U m 2
D 4 = 2 B r D a + ω U m ϵ 24 κ U m 2
D 5 = η 24 κ U m 2 { 2 B r D a 4 + 6 A 1 2 2 + η + η 3 4 A 1 3 + η 2 + { ω [ 4
12 A 2 212 A 2 2 + η + η 3 4 A 1 3 + η 2 24 U m } U m } ϵ
E 1 = B i A 0 2 A 3 2 B r D a A 3 2 B r + 2 D a B r D a ω U m ϵ 4 C 0 κ ϵ U m 2
E 2 = A 3 A 0 2 B i 2 B r D a ω U m κ U m 2 A 0 4 A 0 2 C 0
E 3 = 1 2 C 0 2 ( 4 A 0 4 5 A 0 2 C 0 + C 0 2 ) κ 2 U m 4 ϵ { 4 A 0 2 C 0 C 0 A 3 A 0 2 B i ( 2 B r D a
ω U m ) cosh A 0 η + A 0 2 C 0 4 A 0 2 C 0 2 κ ϵ U m 2 A 6 φ A 3 2 4 A 0 2 B i B r A 0 2 D a + ϵ
+ C 0 C 2 cosh 2 A 6 η ] } sec h C 0 η }
E 4 = A 3 2 4 A 0 2 B i B r A 0 2 D a + ϵ 2 κ ϵ U m 2 16 A 0 4 4 A 0 2 C 0
E 5 = 1 16 κ 2 ϵ 2 U m 2 A 0 2 C 0 2 { 2 A 0 2 [ 4 κ ϵ U m 2 C 0 φ 2 B i { A 0 2 A 3 2 B r D a [ A 3 2 B r
+ 2 D a B r D a ω U m ] ϵ } + C 0 η 2 B i { A 0 2 A 3 2 B r D a A 3 2 B r + 2 D a B r D a ω U m ϵ ] }
+ 4 C 0 A 3 A 0 2 B i 2 B r D a ω U m cosh A 0 η + C 0 A 3 A 0 2 B i 2 B r D a ω U m
cosh 2 A 0 η }
F 1 = A 3 B i 2 B r D a + ω U m κ U m 2 A 0 4 A 0 2 C 0
F 2 = 1 2 κ 2 ϵ U m 4 C 0 2 ( 4 A 0 2 5 A 0 2 C 0 + C 0 2 ) { 4 A 0 2 C 0 C 0 A 3 B i 2 B r D a + ω U m
cosh A 0 η + A 0 2 C 0 [ 4 A 0 2 B i A 0 2 A 3 2 B r D a A 3 2 B r + 2 D a B r D a ω U m ϵ
C 0 B i A 0 2 A 3 2 B r D a A 3 2 B r + 2 D a B r D a ω U m ϵ + C 0 A 3 2 B i B r A 0 2 D a + ϵ
cosh 2 A 0 η } sec h C 0 η κ U m 2 }
F 3 = A 3 2 B i B r A 0 2 D a + ϵ 2 κ ϵ U m 2 16 A 0 4 4 A 0 2 C 0
F 4 = 2 A 0 2 2 + C 0 η 2 + 4 C 0 A 3 B i 2 B r D a + ω U m cosh A 0 η + C 0 cosh A 0 η 4 A 0 2 C 0 2
G 1 = 2 ϵ / { κ { D 1 + D 2 + D 3 + D 4 + D 5 1 U m { 1 840 [ 2 42 D 2 + 35 D 3 + 30 D 4 + 70 D 5
+ 7 [ 5 3 + 8 A 1 + 12 A 2 D 1 + 2 A 2 20 D 2 + 15 D 3 + 12 D 4 + 60 D 5 + A 1 ( 30 D 2 + 24 D 3 + 20 D 4
+ 60 D 5 ) + 1 / A 0 2 C 0 [ A 0 A 3 cosh C 0 η sinh A 0 η + 1 12 A 0 3 [ 24 A 0 A 3 E 1 η cosh A 0 η
+ 12 A 0 2 D a E 2 + 6 A 0 2 A 3 E 4 + 126 A 0 2 A 3 E 5 + 24 A 3 E 1 + 12 A 0 2 η 2 ] sinh A 0 η + 6 A 0 3 A 3 E 2 η
+ 1 / C 0 A 6 A 0 2 E 3 sinh C 0 η [ A 6 A 0 2 D a A 3 A 6 cosh A 0 η ]

References

  1. Bejan, A. Convection in Porous Media; Springer: Berlin/Heidelberg, Germany, 2006; ISBN 9780387290966. [Google Scholar]
  2. Vafai, K. Handbook of Porous Media, 3rd ed.; CRC Press: Boca Raton, FL, USA, 2015; ISBN 9781439885574. [Google Scholar]
  3. Mohamad, A.A. Heat Transfer Enhancements in Heat Exchangers Fitted with Porous Media. Part I: Constant Wall Temperature. Int. J. Therm. Sci. 2003, 42, 385–395. [Google Scholar] [CrossRef]
  4. Mahmoudi, Y.; Maerefat, M. Analytical Investigation of Heat Transfer Enhancement in a Channel Partially Filled with a Porous Material under Local Thermal Non-Equilibrium Condition. Int. J. Therm. Sci. 2011, 50, 2386–2401. [Google Scholar] [CrossRef]
  5. Vafai, K.; Thiyagaraja, R. Analysis of Flow and Heat Transfer at the Interface Region of a Porous Medium. Int. J. Heat Mass Transf. 1987, 30, 1391–1405. [Google Scholar] [CrossRef]
  6. Kiwan, S.; Khodier, M. Natural Convection Heat Transfer in an Open-Ended Inclined Channel-Partially Filled with Porous Media. Heat Transf. Eng. 2008, 29, 67–75. [Google Scholar] [CrossRef]
  7. Yang, K.; Wang, Y.; Mao, Y.; Zhang, W. Heat and Moisture Transfer in a Rectangular Cavity Partially Filled with Hygroscopic Porous Media. Heat Transf. Eng. 2020, 41, 814–824. [Google Scholar] [CrossRef]
  8. Beavers, G.S.; Joseph, D.D. Boundary Conditions at a Natural Permeable Wall. J. Fluid Mech. 1967, 30, 197–207. [Google Scholar] [CrossRef]
  9. Jang, J.Y.; Chen, J.L. Forced Convection in a Parallel Plate Channel Partially with a High Porosity Porous Medium. Int. Comm. Heat Mass Trasnf. 1992, 19, 263–273. [Google Scholar] [CrossRef]
  10. Kuznetsov, A.V. Analytical Investigation of Couette Flow in a Composite Channel Partially Filled with a Porous Medium and Partially with a Clear Fluid. Int. J. Heat Mass Transf. 1998, 41, 2556–2560. [Google Scholar] [CrossRef]
  11. Ochoa-Tapia, J.A.; Whitaker, S. Momentum Transfer at the Boundary between a Porous Medium and a Homogeneous Fluid-I. Theoretical Development. Int. J. Heat Mass Transf. 1995, 38, 2635–2646. [Google Scholar] [CrossRef]
  12. Cekmer, O.; Mobedi, M.; Ozerdem, B.; Pop, I. Fully Developed Forced Convection in a Parallel Plate Channel with a Centered Porous Layer. Transp. Porous Media 2012, 93, 179–201. [Google Scholar] [CrossRef]
  13. Ucar, E.; Mobedi, M.; Pop, I. Effect of an Inserted Porous Layer Located at a Wall of a Parallel Plate Channel on Forced Convection Heat Transfer. Transp. Porous Media 2013, 98, 35–57. [Google Scholar] [CrossRef] [Green Version]
  14. Maerefat, M.; Mahmoudi, S.Y.; Mazaheri, K. Numerical Simulation of Forced Convection Enhancement in a Pipe by Porous Inserts. Heat Transf. Eng. 2011, 32, 45–56. [Google Scholar] [CrossRef]
  15. Guo, Z.; Kim, S.Y.; Sung, H.J. Pulsating Flow and Heat Transfer in a Pipe Partially Filled with a Porous Medium. Int. J. Heat Mass Transf. 1997, 40, 4209–4218. [Google Scholar] [CrossRef]
  16. Amiri, A.; Vafai, K.; Kuzay, T.M. Effects of Boundary Conditions on Non-Darcian Heat Transfer through Porous Media and Experimental Comparisons. Numer. Heat Transf. Part A Appl. 1995, 27, 651–664. [Google Scholar] [CrossRef]
  17. Ochoa-Tapia, J.A.; Whitaker, S. Heat Transfer at the Boundary between a Porous Medium and a Homogeneous Fluid: The One-Equation Model. J. Porous Media 1998, 1, 31–46. [Google Scholar] [CrossRef]
  18. Alazmi, B.; Vafai, K. Analysis of Fluid Flow and Heat Transfer Interfacial Conditions between a Porous Medium and a Fluid Layer. Int. J. Heat Mass Transf. 2001, 44, 1735–1749. [Google Scholar] [CrossRef]
  19. Yang, K.; Vafai, K. Analysis of Temperature Gradient Bifurcation in Porous Media—An Exact Solution. Int. J. Heat Mass Transf. 2010, 53, 4316–4325. [Google Scholar] [CrossRef]
  20. Yang, K.; Vafai, K. Restrictions on the Validity of the Thermal Conditions at the Porous-Fluid Interface—An Exact Solution. J. Heat Transf. 2011, 133, 112601. [Google Scholar] [CrossRef] [Green Version]
  21. Yang, K.; Vafai, K. Analysis of Heat Flux Bifurcation inside Porous Media Incorporating Inertial and Dispersion Effects—An Exact Solution. Int. J. Heat Mass Transf. 2011, 54, 5286–5297. [Google Scholar] [CrossRef]
  22. Xu, H.J.; Qu, Z.G.; Tao, W.Q. Analytical Solution of Forced Convective Heat Transfer in Tubes Partially Filled with Metallic Foam Using the Two-Equation Model. Int. J. Heat Mass Transf. 2011, 54, 3846–3855. [Google Scholar] [CrossRef]
  23. Yang, C.; Nakayama, A.; Liu, W. Heat Transfer Performance Assessment for Forced Convection in a Tube Partially Filled with a Porous Medium. Int. J. Therm. Sci. 2012, 54, 98–108. [Google Scholar] [CrossRef]
  24. Mahmoudi, Y.; Karimi, N.; Mazaheri, K. Analytical Investigation of Heat Transfer Enhancement in a Channel Partially Filled with a Porous Material under Local Thermal Non-Equilibrium Condition: Effects of Different Thermal Boundary Conditions at the Porous-Fluid Interface. Int. J. Heat Mass Transf. 2014, 70, 875–891. [Google Scholar] [CrossRef] [Green Version]
  25. Mahmoudi, Y.; Karimi, N. Numerical Investigation of Heat Transfer Enhancement in a Pipe Partially Filled with a Porous Material under Local Thermal Non-Equilibrium Condition. Int. J. Heat Mass Transf. 2014, 68, 161–173. [Google Scholar] [CrossRef] [Green Version]
  26. Li, Q.; Hu, P. Analytical Solutions of Fluid Flow and Heat Transfer in a Partial Porous Channel with Stress Jump and Continuity Interface Conditions Using LTNE Model. Int. J. Heat Mass Transf. 2019, 128, 1280–1295. [Google Scholar] [CrossRef]
  27. Xu, H.J.; Qu, Z.G.; Lu, T.J.; He, Y.L.; Tao, W.Q. Thermal Modeling of Forced Convection in a Parallel-Plate Channel Partially Filled with Metallic Foams. J. Heat Transfer 2011, 133, 20–22. [Google Scholar] [CrossRef]
  28. Qu, Z.G.; Xu, H.J.; Tao, W.Q. Fully Developed Forced Convective Heat Transfer in an Annulus Partially Filled with Metallic Foams: An Analytical Solution. Int. J. Heat Mass Transf. 2012, 55, 7508–7519. [Google Scholar] [CrossRef]
  29. Torabi, M.; Karimi, N.; Zhang, K. Heat Transfer and Second Law Analyses of Forced Convection in a Channel Partially Filled by Porous Media and Featuring Internal Heat Sources. Energy 2015, 93, 106–127. [Google Scholar] [CrossRef] [Green Version]
  30. Lu, W.; Zhang, T.; Yang, M. Analytical Solution of Forced Convective Heat Transfer in Parallel-Plate Channel Partially Filled with Metallic Foams. Int. J. Heat Mass Transf. 2016, 100, 718–727. [Google Scholar] [CrossRef]
  31. Yang, K.; Chen, H.; Wang, J. Analysis of Heat Transfer and Entropy Generation in a Channel Partially Filled with N-Layer Porous Media. J. Heat Transf. 2018, 140, 082601. [Google Scholar] [CrossRef]
  32. Karimi, N.; Agbo, D.; Khan, A.T.; Younger, P.L. On the Effects of Exothermicity and Endothermicity upon the Temperature Fields in a Partially-Filled Porous Channel. Int. J. Therm. Sci. 2015, 96, 128–148. [Google Scholar] [CrossRef]
  33. Torabi, M.; Zhang, K.; Yang, G.; Wang, J.; Wu, P. Heat Transfer and Entropy Generation Analyses in a Channel Partially Filled with Porous Media Using Local Thermal Non-Equilibrium Model. Energy 2015, 82, 922–938. [Google Scholar] [CrossRef]
  34. Baig, M.F.; Chen, G.M.; Tso, C.P. Viscous Dissipative Forced Convection in a Channel Partially Filled with Porous Medium. J. Thermophys. Heat Transf. 2021, 36, 276–290. [Google Scholar] [CrossRef]
  35. Dehghan, M. Effects of Heat Generations on the Thermal Response of Channels Partially Filled with Non-Darcian Porous Materials. Transp. Porous Media 2015, 110, 461–482. [Google Scholar] [CrossRef]
  36. Hunt, G.; Karimi, N.; Torabi, M. Analytical Investigation of Heat Transfer and Classical Entropy Generation in Microreactors—The Influences of Exothermicity and Asymmetry. Appl. Therm. Eng. 2017, 119, 403–424. [Google Scholar] [CrossRef]
  37. Bhargavi, D.; Reddy, J.S.K. Analytical Investigation of Laminar Forced Convection with Viscous Dissipation in Parallel Plate Channels Partially Filled with Porous Material: Constant Wall Heat Flux. J. Nanofluids 2019, 8, 238–251. [Google Scholar] [CrossRef]
  38. Hu, P.; Li, Q. Effect of Heat Source on Forced Convection in a Partially-Filled Porous Channel under LTNE Condition. Int. Commun. Heat Mass Transf. 2020, 114, 104578. [Google Scholar] [CrossRef]
  39. Al-Hadhrami, A.K.; Elliott, L.; Ingham, D.B. A New Model for Viscous Dissipation in Porous Media across a Range of Permeability Values. Transp. Porous Media 2003, 53, 117–122. [Google Scholar] [CrossRef]
  40. Chen, G.M.; Tso, C.P. Forced Convection with Viscous Dissipation Using a Two-Equation Model in a Channel Filled by a Porous Medium. Int. J. Heat Mass Transf. 2011, 54, 1791–1804. [Google Scholar] [CrossRef]
  41. Chen, G.M.; Tso, C.P. A Two-Equation Model for Thermally Developing Forced Convection in Porous Medium with Viscous Dissipation. Int. J. Heat Mass Transf. 2011, 54, 5406–5414. [Google Scholar] [CrossRef]
  42. Ting, T.W.; Hung, Y.M.; Guo, N. Viscous Dissipative Forced Convection in Thermal Non-Equilibrium Nanofluid-Saturated Porous Media Embedded in Microchannels. Int. Commun. Heat Mass Transf. 2014, 57, 309–318. [Google Scholar] [CrossRef]
  43. Farrukh B, M.; Chen, G.M.; Tso, C.P. Viscous Dissipation Effect on CuO-Water Nanofluid-Cooled Microchannel Heat Sinks. Case Stud. Therm. Eng. 2021, 26, 101159. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the problem.
Figure 1. Schematic diagram of the problem.
Nanomaterials 12 03821 g001
Figure 2. Maximum velocity achievable with porous medium thickness by having D a as a parameter.
Figure 2. Maximum velocity achievable with porous medium thickness by having D a as a parameter.
Nanomaterials 12 03821 g002
Figure 3. Temperature distributions for fluid and solid phases by having B r as a parameter for D a = 10 4 (a) κ = 10 ,   B i = 0.1 ; (b)   κ = 10 ,   B i = 10 ; (c) κ = 0.1 ,   B i = 0.1 ; (d) κ = 0.1 ,   B i = 10 .
Figure 3. Temperature distributions for fluid and solid phases by having B r as a parameter for D a = 10 4 (a) κ = 10 ,   B i = 0.1 ; (b)   κ = 10 ,   B i = 10 ; (c) κ = 0.1 ,   B i = 0.1 ; (d) κ = 0.1 ,   B i = 10 .
Nanomaterials 12 03821 g003aNanomaterials 12 03821 g003b
Figure 4. Nusselt number versus thickness of porous medium [24] (a) κ = 0.1 , B i = 0.1 ; (b)   κ = 0.1 ,   B i = 10 ; (c) κ = 10 ,   B i = 0.1 ; (d) κ = 10 ,   B i = 10 .
Figure 4. Nusselt number versus thickness of porous medium [24] (a) κ = 0.1 , B i = 0.1 ; (b)   κ = 0.1 ,   B i = 10 ; (c) κ = 10 ,   B i = 0.1 ; (d) κ = 10 ,   B i = 10 .
Nanomaterials 12 03821 g004
Figure 5. Nusselt number versus Brinkman number with porous medium thickness η = 0.9 . (a) κ = 0.1 ,   B i = 0.1 (b)   κ = 0.1 ,   B i = 10 (c) κ = 10 ,   B i = 0.1 (d) κ = 10 ,   B i = 10 .
Figure 5. Nusselt number versus Brinkman number with porous medium thickness η = 0.9 . (a) κ = 0.1 ,   B i = 0.1 (b)   κ = 0.1 ,   B i = 10 (c) κ = 10 ,   B i = 0.1 (d) κ = 10 ,   B i = 10 .
Nanomaterials 12 03821 g005
Figure 6. Schematic diagram of the composite microchannel.
Figure 6. Schematic diagram of the composite microchannel.
Nanomaterials 12 03821 g006
Figure 7. Heat transfer coefficient as a function of porous medium thickness with different volume fractions of Al2O3-Water nanofluid (a) R e = 25 (b)   R e = 50 .
Figure 7. Heat transfer coefficient as a function of porous medium thickness with different volume fractions of Al2O3-Water nanofluid (a) R e = 25 (b)   R e = 50 .
Nanomaterials 12 03821 g007
Figure 8. Temperature line plots in Kelvin for Re = 25 in a microchannel. (a) 0% volume fraction Al2O3-Water without viscous dissipation. (b) 0% volume fraction Al2O3-Water with viscous dissipation. (c) 4% volume fraction Al2O3-Water without viscous dissipation. (d) 4% volume fraction Al2O3-Water with viscous dissipation.
Figure 8. Temperature line plots in Kelvin for Re = 25 in a microchannel. (a) 0% volume fraction Al2O3-Water without viscous dissipation. (b) 0% volume fraction Al2O3-Water with viscous dissipation. (c) 4% volume fraction Al2O3-Water without viscous dissipation. (d) 4% volume fraction Al2O3-Water with viscous dissipation.
Nanomaterials 12 03821 g008aNanomaterials 12 03821 g008b
Table 1. Nusselt number comparison with literature for LTNE model.
Table 1. Nusselt number comparison with literature for LTNE model.
Thickness of
Porous   Medium ,   η
Brinkman   Number ,   B r Biot   Number ,   B i Kappa ,   κ Darcy   Number ,   D a N u ,   Present   Study N u [40] N u [41]
100.01100104.1224.117-
1010.011123.9-123.9
101100 10 4 5.901-5.901
10110014.171-4.171
110.01100100.24440.2439-
11100.011132.4-132.4
1110.01149.33-49.33
11110011.622-1.622
1110.010.01143.3-143.3
1100.01100100.02580.0257-
11010.010.01122.6-122.6
11010.010.0485.37-85.37
Table 2. Heat flux bifurcation for fluid and solid phases in porous region at interface η = 0.9 .
Table 2. Heat flux bifurcation for fluid and solid phases in porous region at interface η = 0.9 .
B r D a = 10 3 D a = 10 4
θ f θ s θ f θ s
κ = 0.1 ,   B i = 0.1 101.64500.16855−3.31450.09494
20−3.438280.17674−10.3260.10406
κ = 0.1 ,   B i = 10 10−2.91010.62406−5.97370.36085
20−7.56180.58943−12.9870.37010
Table 3. Specifications of a microchannel [42].
Table 3. Specifications of a microchannel [42].
NanofluidAl2O3-Water
SolidSilicon
Length of the channel, L (m)0.03
Channel height, 2 H   μ m 50
Heat flux, q w (W/m2)1 × 104
Thermal conductivity of solid, k s (W/m·K)148
Porosity, ϵ 0.9
Darcy number, D a 0.1
Biot number, B i 0.01
Table 4. Thermophysical properties of nanofluid [42].
Table 4. Thermophysical properties of nanofluid [42].
PropertiesNanoparticle
(Al2O3)
Base Fluid
(Water)
ρ (kg/m3)3975997
c p (J/kg·K)778.64179
k (W/m·K)360.613
μ (N·s/m2)- 8.55 × 10−4
Table 5. Comparison of heat transfer coefficient for a composite channel.
Table 5. Comparison of heat transfer coefficient for a composite channel.
Working Fluid R e   B r   h ,   W / m 2 · at
η = 0.45  
%   Decrease   in   h
Water250116,1910.10
0.006116,086
4% volume fraction Al2O3-Water nanofluid0132,8002.44
0.064129,550
Water500116,1910.36
0.025115,773
4% volume fraction Al2O3-Water nanofluid0132,8005.62
0.255125,054
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Baig, M.F.; Chen, G.M.; Tso, C.P. The Thermal Performance Analysis of an Al2O3-Water Nanofluid Flow in a Composite Microchannel. Nanomaterials 2022, 12, 3821. https://doi.org/10.3390/nano12213821

AMA Style

Baig MF, Chen GM, Tso CP. The Thermal Performance Analysis of an Al2O3-Water Nanofluid Flow in a Composite Microchannel. Nanomaterials. 2022; 12(21):3821. https://doi.org/10.3390/nano12213821

Chicago/Turabian Style

Baig, Mirza Farrukh, Gooi Mee Chen, and Chih Ping Tso. 2022. "The Thermal Performance Analysis of an Al2O3-Water Nanofluid Flow in a Composite Microchannel" Nanomaterials 12, no. 21: 3821. https://doi.org/10.3390/nano12213821

APA Style

Baig, M. F., Chen, G. M., & Tso, C. P. (2022). The Thermal Performance Analysis of an Al2O3-Water Nanofluid Flow in a Composite Microchannel. Nanomaterials, 12(21), 3821. https://doi.org/10.3390/nano12213821

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