Forced convection in micro-channels filled with porous media in local thermal non-equilibrium conditions
Introduction
In recent years, research activity in heat transfer at micro and nanoscale geometries has been strongly developing due to the incredible growth of micro-electro-mechanical systems (MEMS). Several engineering and biomedical applications have determined an increasing research interest in micro and nano flow, as recently reviewed in Refs. [1], [2], [3], [4], [5], [6], [7], [8], [9], [10], [11]. A general introduction on the importance and developments of microfluidics is reported in Ref. [1]. A complete description of the fundamentals in microchannels is provided in Ref. [2]. A bio-MEMS application in terms of nanofluid flow in microchannels is presented in Ref. [3]. A review focused on molecular momentum transport at fluid–solid interfaces mainly related to microfluidics and nanofluidics in micro-/nano-electromechanical systems (MEMS/NEMS) is performed in Ref. [4]. A survey on existing methods for the characterization of mixing and flow in microchannels, micromixers and microreactors is accomplished in Ref. [5]. A review which emphasizes gas–liquid two-phase flow in minichannels or microchannels related to PEM fuel cell applications is reported in Ref. [6]. An analysis of published data on a gas–liquid two-phase flow in capillaries of various shapes to systematize the collected body of information is provided in Ref. [7]. A review on one and two phases flows and methods to improve heat transfer in microchannels is given in Ref. [8]. A comprehensive review on the work done regarding to heat transfer and fluid flow behaviors in microchannels heat exchanger using conventional fluids as well as nanofluids is described in Ref. [9]. The microscale transport processes that arise in the fabrication of advanced materials are reviewed in Ref. [10]. The basic considerations on transport phenomena at micro- or nanoscale levels are outlined. Different approaches found in literature for measuring microflow characteristics and the geometry of the microchannels are explored and categorized in Ref. [11]. It is pointed out the advantages and disadvantages inherent to each experimental technique.
At micro and nanoscales the phenomenon of rarefaction or Knudsen flow have been taken into account and the Navier–Stokes but energy equations in the continuum flow model could not be completely appropriated to describe the fluid flow and heat transfer [12], [13], [14], [15], [16], [17]. A critical review on the status of the understanding of fluid flow phenomena in microdevices is given in Ref. [12]. A new general boundary condition that accounts for the reduced momentum and heat exchange with wall surfaces in a wide range of Knudsen number at low Mach number is proposed and its validity is investigated in Ref. [13]. The importance of different approaches in terms of length scales for macro and micro scales is underlined in Refs. [14], [15]. A review of the available experimental works on the convective heat transfer through microchannels is presented in Ref. [16]. A survey on the main theoretical and experimental results about steady pressure-driven gas microflows is reported in Ref. [17].
The modeling in microflow depends on the Knudsen number, Kn, defined as the ratio of the fluid mean free path to the macroscopic length scale of the physical system, λ/Lref. It allows to have a measure of the degree of rarefaction of gases which flow through very small channels and, consequently, of the degree of the validity of the continuum model. A classification of gas flow regimes depends on the Knudsen number values and the following expressions are usually accepted, as indicated in Refs. [2], [17], [18], [19]:
Kn < 10−3 for the continuum no-slip flow,
10−3 < Kn < 10−1 for the continuum slip flow,
10−1 < Kn < 101 for the transition flow and
Kn > 10 for the molecular regime.
A greater part of the investigations, both experimental and numerical, have been accomplished on forced convection gas flow in order to evaluate the pressure drop and the convective heat transfer in microchannels and microtubes [2], [5], [6], [9], [12], [13], [14], [15], [16], [17], [18], [19], [20], [21], [22], [23], [24], [25], [26], [27], [28], [29].
Several studies and analysis on heat transfer in porous media are based on local thermal equilibrium (LTE) condition, assuming that fluid phase and solid phase local temperatures are equal and one energy equation model is considered. In many real cases temperature differences between the two phases, fluid and solid, are not negligible and the one equation model is not valid. In such cases, the solid phase and the fluid phase have different local temperatures and it is assumed a local thermal non-equilibrium (LTNE) condition in the porous medium which involves a two energy equations model, one for the fluid phase and one for the solid phase. The temperature difference determines a convective heat transfer between solid and fluid phase inside the porous medium, i.e. an inter-phase heat transfer or interstitial heat transfer coefficient.
Research interest has been recently developing in microchannels and microtubes filled with porous medium due to their applications in micro filtration, fractionation, catalysis and microbiology, as underlined in Ref. [30]. However, forced convection in porous channels and tubes, in rarefied condition or in porous microchannels and microtubes, have not been widely studied. Analytical and/or numerical solutions on different geometrical and thermal conditions should be examined. Moreover, some phenomenological aspects related to LTNE in porous media inside microchannels and microtubes have received only small attention.
One of the first papers on porous microchannel was accomplished by Haddad et al. [31]. Fluid dynamics and thermal behaviors of developing natural convection in a vertical open-ended parallel-plate microchannel, filled with porous media, were numerically investigated considering the slip flow regime. Local thermal non-equilibrium (LTNE) was assumed and the flow in porous medium was modeled employing the Darcy–Brinkman–Forchheimer model which accounts for viscous and inertia effects. An analytical solutions for forced convection with slip-flow in parallel plate channel and circular duct saturated by a rarefied gas was carried out by Kuznetsov and Nield [32]. They considered constant heat flux at the walls and studied the fully developed situation in local thermal equilibrium (LTE). Forced convection slip flow in parallel-plate channel, at uniform wall temperature, filled with porous media was studied numerically by Haddad et al. [33]. The analysis was accomplished in slip flow regime for Kn values from 10−3 to 10−1, assuming LTE condition and the Darcy–Brinkman–Forchheimer model for momentum equation. Transient Couette flow, pulsating Poiseuille flow, Stokes second problem flow and transient natural convection flow in a porous media, in the slip flow regime, were investigated by Haddad et al. [34]. In the study the Darcy–Brinkman model, which accounts for viscous effects, and LTE condition were assumed. Results were carried out in the Knudsen number range 10−3–10−1. The analysis of the problem investigated in Ref. [33] was extended to LTNE by Haddad et al. [35]. Forced convection slip flow in a circular micro-channel filled with porous media was accomplished by Haddad et al. [36]. Wall uniform temperature, Darcy–Brinkman–Forchheimer model and LTE condition were assumed. The investigation was performed for Kn values between 10−3 and 10−1.
Stability analysis on a slip flow in a channel filled with a hyperporous medium, saturated by a rarefied gas, was studied by Avramenko et al. [37]. Numerical analysis for two and three-dimensional disturbances were carried out using the collocation method and the LTE was assumed. Fully developed forced convection in a rectangular microchannel, filled with or without a porous medium in LTE, was investigated analytically by Hooman [38], [39]. The Darcy–Brinkman flow model, slip flow and prescribed local heat flux, with H2 boundary condition [38] or H1 boundary condition [39], were assumed. The analysis was accomplished in the range: 10−3 ≤ Kn ≤ 10−1. Forced convection in a parallel plate channel or a circular tube filled by a hyperporous medium saturated by a rarefied gas, with walls held at constant heat flux was investigated by Kuznetsov and Nield [40]. They found analytic solutions for thermally developing flow in LTE condition. Fully developed forced convection in a circular channel filled with a highly porous medium, saturated with a rarefied gas and in the slip-flow regime, was studied by Chauhan and Kumar [41]. Uniform heat flux at tube wall, Darcy–Brinkman–Forchheimer momentum equation and local thermal equilibrium were assumed. The entropy generation due to heat transfer and fluid friction was also evaluated. A numerical simulation of hydrodynamically developed and thermally developing forced convection in a circular micro/nano channel filled with porous media was performed by Shokouhmand et al. [30]. The Darcy–Brinkman–Forchheimer flow model, slip boundary condition and LTE condition were assumed. A unified flow model was considered for all flow regimes in a large Knudsen numbers range. Hydrodynamically and thermally fully developed flow of a dilute gas in forced convection in a porous annular microduct was investigated analytically by Hashemi et al. [42]. Darcy–Brinkman equation and local thermal equilibrium were used. Uniform heat flux at the outer cylinder and insulated inner cylinder and vice versa were considered as thermal boundary conditions.
The main contributions on the microchannels filled with porous media are summarized in Table 1. Except in Refs. [35], it seems that there are no studies on microchannels or channels in rarefied gas filled with the assumption of LTNE on porous media, mainly for assigned heat flux at the microchannel or channel walls.
As indicated by Hashemi et al. [42], several investigations on microchannel heat sinks analyze these devices by means of porous media models in LTNE [43], [44], [45], [46], [47], [48], [49], [50]. The no-slip condition is, however, assumed in all studied cases. Moreover, very few works in microchannels filled with porous media have been accomplished with the entropy generation analysis [41] which was carried out in LTE.
In the present study a hydrodynamically and thermally fully developed flow of a dilute gas in a parallel-plate microchannel, assuming the slip-flow and jump temperature on the microchannel walls, filled with a porous medium, is carried out adopting the LTNE condition on the porous medium. An analytical solution for uniform wall heat flux and Darcy–Brinkman model is evaluated and the entropy generation analysis is performed. Results, for Knudsen number, provided that 10−3 ≤ Knp ≤ 10−1 (Knp = λ/K1/2), in terms of dimensionless wall temperature profiles, Nusselt number, dimensionless fluid and solid phases and entropy generation number are given for various parameters values such as Darcy number, Biot number and thermal conductivity ratio. It should be underlined that the Biot number is referred to the convective heat transfer between the solid matrix and fluid at solid–fluid interface in the porous medium i.e. the interstitial heat transfer coefficient [35]. Moreover, the interstitial heat transfer coefficient is an assigned value.
Section snippets
Governing equations and mathematical model
The considered physical problem is a forced convection steady laminar flow through a parallel plates micro-channel filled with a porous medium as shown in Fig. 1. The distance between the two parallel walls is 2H and a uniform wall heat flux is applied on both the channel walls. Fluid flow is assumed incompressible and hydrodynamically and thermally fully developed with thermo-physical properties constant with temperature. Moreover, buoyancy force, and consequently natural convection,
Solution procedure
In dimensionless form, the momentum Eq. (1) and related boundary conditions (4), (5) are as follow:by employing the following dimensionless variables:
The solution of the problem (8), (9) is:with
The friction factor is defined as follows [65]:by definition of G and Eq. (10), it is obtained
The average velocity, the bulk temperature and
Entropy generation analysis
In irreversible process, such as heat convection in a porous medium, entropy generation caused by fluid motion due to the viscous effect and heat transfer under finite temperature gradients is always present. Following Bejan [68] and taking into account the analysis proposed in Refs. [41], [69], [70] for porous media in LTE [41], [70] and in LTNE [69], the volumetric rate of entropy generation for the present problem of forced convection in parallel plates microchannels, filled with a porous
Results and discussion
In the present investigation, results are given for 10−2 ≤ Bi ≤ 106, over the following ranges: 0 ≤ ω ≤ 103 and 10−4 ≤ κ ≤ 102 with α and β in the [0,1] range. As previously emphasized in Eqs. (10), (20), α and β are a combination of accommodation coefficients, Knudsen and Prandtl numbers and coefficient γ. In the following analysis the considered Knp range is from Knp < 10−3 (no slip) to Knp = 0.1 (slip continuum).
A comparison with the two following asymptotic cases is reported in Fig. 2: 1)
Conclusions
Fully developed, steady-state forced convection, in parallel plates micro-channels, filled with a porous medium saturated with rarefied gas, in local thermal non-equilibrium condition was investigated in the slip flow regime (10−3 ≤ Knp ≤ 0.1). An analytic solution was proposed for the Darcy–Brinkman model with assigned uniform heat flux at micro-channel walls. Moreover, the entropy generation analysis was provided both in the rarefied gas and in solid matrix. The effects of the main
Acknowledgments
This work was funded by MIUR with Articolo 12 D. M. 593/2000 Grandi Laboratori ‘‘EliosLab” and by Italian Government, MIUR grant PRIN-2009KSSKL3. The authors would also like to thank very much the Reviewers; their very interesting and useful comments and suggestions have significantly allowed to improve the article. Oronzio Manca thanks to UPEMLV for the support during his stay at Laboratoire Modélisation et Simulation Multi-Echelle as Visiting Professor.
References (73)
- et al.
Microfluidics of nano-drug delivery
Heat Mass Transfer
(2008) - et al.
Current methods for characterizing mixing and flow in microchannels
Chem. Eng. Sci.
(2010) - et al.
A critical review of two-phase flow in gas flow channels of proton exchange membrane fuel cells
J. Power Sources
(2010) - et al.
Heat transfer and fluid flow characteristics in microchannels heat exchanger using nanofluids: a review
Renew. Sustain. Energy Rev.
(2011) - et al.
A critical review of the measurement techniques for the analysis of gas microflows through microchannels
Exp. Therm. Fluid Sci.
(2011) Comments on “Critical view on new results in micro-fluid mechanics”
Int. J. Heat Mass Transfer
(2003)- et al.
Size effect on microscale single-phase flow and heat transfer
Int. J. Heat Mass Transfer
(2003) Single-phase convective heat transfer in microchannels: a review of experimental results
Int. J. Therm. Sci.
(2004)- et al.
Analytical and numerical studies on microscale heat sinks for electronic applications
Appl. Therm. Eng.
(2005) - et al.
Size effect on single-phase channel flow and heat transfer at microscale
Int. J. Heat Fluid Flow
(2003)