The extensive engineering applications of Micro-Electro-Mechanical systems have promoted abundant studies of its fluid flow and heat transfer characteristics. This work presented numerical simulation of fully developed flow for characteristic laminar slip flow and heat transfer in rhombus microchannels with the presence of radiation. The slip velocity and temperature jump boundary conditions at the wall are performed. The effects of Reynolds number, velocity slip and temperature jump on Poiseuille number and impact of the presence of radiation on Nusselt number for different aspect ratio is reported.
Analysis of fluid transport at microscale is of a great importance not only for playing a key role in the biological systems, but also in a wide variety of contemporary engineering applications, involving micro-scale devices such as microsensors, micropumps, and cooling electronic equipments. Microchannels however are the basic structure of these systems. People realize that there are many differences between heat and flow transfer in Microchannels and channels with conventional size.
Get your grade
or your money back
using our Essay Writing Service!
A number of practical situations involve heat transfer between a low-density gas and a solid surface. In employing the term low density, we shall mean those circumstances where the mean free path of the gas molecules is no longer small in comparison with a characteristic dimension of the heat transfer surface. The mean free path is the distance a molecule travels, on the average, between collisions. The larger this distance becomes, the greater the distance required to communicate the temperature of a hot surface to a gas in contact with it. Evidently, the parameter which is of principal interest is a ratio of the mean free path to a characteristic body dimension. This grouping is called the Knudson number which is defined as:
Beskok and Karniadakis  classified the gas flow in microchannels into four flow regimes. They are continuum flow regime (Kn≤0.001), slip flow regime (0.001<Kn≤0.1), transition flow regime (0.1<Kn≤10) and free molecular flow regime (kn>10).
In slip flow regime, the Navier-Stokes and energy equations remain applicable, provided a velocity slip and temperature jump are taken into account at the walls. Researchers have investigated theoretical, numerical and experimental studies in microchannels along with velocity slip and temperature jump boundary condition at the walls [2-8]. Numerical analysis of fully developed laminar slip flow and heat transfer in trapezoidal microchannels had been studied by Bin et al.  with uniform wall heat flux boundary condition. In their investigation, the influence of velocity slip and temperature jump on friction factor and Nusselt number were investigated in detail. M. Shams and C. Aghanajafi has studied numerical simulation of slip flow through rhombus microchannels without the impact of heat transfer through radiation . Wei et al.  investigated the steady state convective heat transfer for laminar, two dimensional, incompressible rarefied gas flow by the finite volume finite difference scheme with slip flow and temperature jump boundary conditions.
In this paper, a fully developed laminar flow in a rhombus microchannel with the presence of radiation is investigated. Numerical results are obtained using a continuum based three dimensional, incompressible, steady model which is solved by a finite-volume method with slip velocity and temperature jump boundary conditions applied to the momentum and energy equations, respectively. The effects of rarefaction and aspect ratio on Poiseuille number and Nusselt number are studied. The influence of the presence of radiation on Nusselt number is obtained.
2. Model setup
In this paper, microchannels with rhombus cross section are analyzed. The simulations are performed based on the following assumptions:
(1) The governing equations based on Navier-Stokes equations associated with slip boundary conditions that can describe slip flow regimes in microchannels.
(2) The flow is laminar.
(3) The process is three-dimensional incompressible steady flow.
(4) The body forces are neglected.
(5) The energy equation with the presence of radiation is applied.
The resulting governing equations are following:
where, u is the velocity component, p is pressure, ρ is the fluid density, µ is the dynamic viscosity.
where, T is temperature, and cp is the specific heat at constant pressure.
where is the Stefan-Boltzman constant and k is the Rosseland mean absorption coefficient. Assuming that the optical depth of the gas is sufficiently large and the temperature gradients are sufficiently small so that the local intensity results from local emissions and also can be expressed as a linear function of temperature
Always on Time
Marked to Standard
where the higher order terms of the expansion are neglected.
The most common means of analytically or numerically modeling a rarified flow within the slip flow regime, 0.01≤Kn≤0.1, is through the use of slip velocity and temperature jump boundary conditions applied to the conventional continuum momentum and energy equations. The original slip velocity boundary condition and temperature jump boundary condition were derived by Maxwell  and Smoluchowski , respectively.
where, u is the tangential velocity, λ is the mean-free path of the molecules, σν is the tangential momentum accommodation coefficient, σT is the thermal accommodation coefficient, γ is the specific heat ratio of the gas, Pr is the Prandtl number, µ is the viscosity, ρ is the density and T is the temperature of the gas at the wall.
The first term in eq.(7) is the velocity slip due to the shear stress at the solid surface, and the second term is the thermal creep velocity due to a temperature gradient tangential to the wall. These equations are presented in a format assuming a Cartesian coordinate system, a wall normal direction (y), and a streamwise direction (x).
It should be noted that the streamwise velocity and the temperature at the inlet are uniform and other component of velocity are zero. The wall temperature is constant. The velocity gradients along the axial flow direction are zero which means the boundary condition at the outlet boundary is outflow. Obviously, the channel length is chosen so that the developing length can be neglected in comparison with the developed flow length and there is thermally and hydrodynamically developed flow at the outlet boundary.
The tangential momentum accommodation and the thermal accommodation coefficients vary between zero and unity. The tangential momentum accommodation coefficient is zero for specular molecular reflection at the wall and it is unity for diffuse reflection . The thermal accommodation coefficient is zero if molecules hold their original temperature upon collision and it is unity if they gain the wall temperature. In this paper these coefficients are assumed as unity. The specific heat ratio of the gas γ is equal to 1.4.
3. Numerical model
In this paper the rhombus microchannel is studied. The aspect ratio AR is defined as:
where, H and b are horizontal and vertical diameters of the rhombus channel, respectively. In order to reach fully developed condition the length is sufficiently long.
The 3D, incompressible, laminar, steady state momentum and energy equations are used. The commercial CFD software FLUENTTM  version 6.0.12 is employed to solve the governing equations. The SIMPLEC is used to link the pressure and velocities. In order to ensure the fineness of the computational grid spacing a grid independency test is performed. To investigate the influences of the slip boundary conditions a code is included in FLUENTTM v.6.0.12. The residuals are considered 10-6 to reach accurate results and convergence of the solution. A grid independency check is included to ensure that the solution is independent of grid size.
To establish that the numerical procedure accurately models the slip flow regime, some comparisons are performed. Poiseuille numbers have been compared with analytical values according to Po=24/(1+12Kn) in Table 1 that shows good agreement. In Table 2, the Nusselt number have been compared with the results of Reinksizbulut  for Kn=0 through a range of Reynolds number. The comparisons have been performed for parallel plates. The aspect ratio has been chosen sufficiently large that the channel geometry tends to parallel plates. Another case of validation is presented in Table 3 in which contains the results of this research and the results of Morini. The results are completely in accord with the Morini's results for different value of Φ (Φ is the ratio of Poiseuille number at Kn=0 to Poiseuille number at non-zero Kn).
Fully developed Poiseuille number for parallel plates.
Fully developed Nusselt number for parallel plates.
Results of Reinksizbulut 
This Essay is
a Student's Work
This essay has been submitted by a student. This is not an example of the work written by our professional essay writers.Examples of our work
Fully developed Poiseuille number for rectangular microchannel.
Results of 
Results of 
Present results Φ
Results of 
Present results Φ
Results of 
Present results Φ
5. Results and discussion
5.1. Fully developed non-dimensional velocity
The non-dimensional maximum velocity and velocity slip for AR=1 and varying Knudson number are shown in Fig. 1. As shown in this figure the velocity at the wall is zero for Kn=0. This figure shows that the velocity slip at the wall increases with increasing Kn, which means that the slip velocity increases as the effect of rarefaction becomes more significant. Also it is obvious in this figure that the maximum velocity in the center of the channel decreases, when Knudson number increases. The variations of the fully developed non-dimensional maximum velocity and velocity slip with different aspect ratio and Kn=0.01 are shown in Fig. 2. As shown in this figure, the maximum velocity in the center of the channel decreases when the aspect ratio of the channel increases. It is shown that, increasing the value of aspect ratio augments the slip velocity.
5.2. Hydrodynamic characteristics
Poiseuille number is one of the most important parameters in fluid flow which is the product of friction factor and local Reynolds number. The effects of rarefaction, aspect ratio and Reynolds number on Poiseuille number on rhombus microchannels are shown below. Fig. 3 shows the influence of Reynolds number on Poiseuille number for different rarefactions and unity aspect ratio. As seen in Fig.
3, the variations of Reynolds number have little effect on Poiseuille number for a fixed rarefaction. It can be concluded from this figure that Poiseuille number decreases with increasing value of Knudson number. Fig. 4 shows the Poiseuille number as a function of aspect ratio for different Reynolds number and at Kn=0.01. As seen in Fig. 4, Poiseuille number increases with increasing value of aspect ratio, and as shown in this figure, the variations of Reynolds number have little influences on Poiseuille number. Variation of fully developed Poiseuille number with aspect ratio for different Knudsen number at Re=10 is the purpose of Fig. 5. It is obvious in Fig. 5 that Poiseuille increases with increasing value of aspect ratio for a fixed Knudsen number. As seen in this figure the aspect ratio is more effective when it is less than 0.7 and at greater values than this aspect ratio, the Poiseuille number is constant for a fixed rarefaction. This figure also shows that Knudsen number can considerably affect on Poiseuille number, in which Poiseuille number reaches its maximum when Knudsen number is equal to zero, i.e. no slip condition, and then increasing the value of rarefaction reduces the value of Poiseuille number.
5.3. Heat transfer
The most important parameter in heat transfer is the ratio of convective to conductive heat transfer across the boundary which is called Nusselt number. The average local heat transfer or the local Nusselt number is defined as:
where, h is the heat transfer coefficient, is the circumferentially average temperature gradient to the wall and is the bulk mean temperature.
In order to study the influence of the presence of Radiation the heat transfer coefficient is defined as:
where, is the convective heat transfer coefficient and is the heat transfer coefficient of radiation. Obviously, when the radiation is neglected, the total heat transfer coefficient is equal to convective heat transfer coefficient.
The figures following are the results of investigations which show the effects of Reynolds number, Prandtl number, Knudsen number and channel aspect ratio on fully developed Nusselt number for microchannel with rhombus cross section. The effect of the presence of radiation on Nusselt number is shown and discussed in details. Fig. 6 shows the effect of Prandtl number on Nusselt number for different aspect ratio at Re=10 and Kn=0.01. It is obvious that increasing the value of Prandtl number causes the reduction in temperature jump and consequently reduces the heat transfer rate and the Nusselt number. It is shown in this figure that the Nusselt number increases by increasing the value of aspect ratio of the channel. Also it can be seen in this figure that the presence of the radiation augments the value of Nusselt number. Fig. 7 shows the Nusselt numbers as a function of aspect ratio for different Reynolds
numbers for Kn=0.01 and Pr=0.74.it is clear in this figure that the heat transfer rate decreases with increasing the value of Reynolds number and this reduction is greater between lower value of Reynolds number. Also it is shown that the Nusselt number increases when the aspect ratio of the channel increases for a fixed Reynolds number. The influence of the presence and absence of radiative heat transfer is also depicted in this figure. Fig. 8 shows the Nusselt number as a function of rarefaction for different Reynolds number at Pr=0.74 and aspect ratio of unity. The effect of radiation is clearly shown in this
figure. It is shown in this figure that the value of Nusselt number is greater for lower Reynolds number, and also the Nusselt number decreases with increasing the value of rarefaction. The variation of fully developed Nusselt number as a function of aspect ratio and Knudsen number is investigated and shown in Fig. 9 at Re=10 and Pr=0.74. As seen in Fig. 9, the Nusselt number has its maximum value when there are no rarefaction and it is obvious in this figure that the rarefaction causes reduction in heat transfer rate and Nusselt number for a fixed aspect ratio. It is clearly shown in this figure that decreasing the value of aspect ratio for a fixed Knudsen number affects the Nusselt number in which it reduces and lower aspect ratios have more considerable effects on this reduction. We can also see in this figure that the presence of radiation has a magnificent effect on the heat transfer rate i.e. Nusselt number.
A three dimensional incompressible steady model has been developed to investigate the effects of Reynolds number, aspect ratio, rarefaction, Prandtl number and radiation heat transfer on Poiseuille and Nusselt number for fully developed laminar flow over a range of slip flow regime in a microchannel with rhombus cross section. Constant wall temperature, slip velocity and temperature jump boundary conditions were included. The Nusselt and Poiseuille number both decrease with decreasing the value of aspect ratio, but rarefaction has inverse effect on Poiseuille and Nusselt number in which Poiseuille and Nusselt number decreases with increasing in Kn. It is observed from the figures that Reynolds number has more effect on Nusselt number then Poiseuille number. It is concluded from the figures and results that the presence of the radiation augments the heat transfer rate and Nusselt number to about 3% in comparison by the absence of radiation.