Heat Transfer of Oil / MWCNT Nanoﬂuid Jet Injection Inside a Rectangular Microchannel

: In the current study, laminar heat transfer and direct ﬂuid jet injection of oil / MWCNT nanoﬂuid were numerically investigated with a ﬁnite volume method. Both slip and no-slip boundary conditions on solid walls were used. The objective of this study was to increase the cooling performance of heated walls inside a rectangular microchannel. Reynolds numbers ranged from 10 to 50; slip coe ﬃ cients were 0.0, 0.04, and 0.08; and nanoparticle volume fractions were 0–4%. The results showed that using techniques for improving heat transfer, such as ﬂuid jet injection with low temperature and adding nanoparticles to the base ﬂuid, allowed for good results to be obtained. By increasing jet injection, areas with eliminated boundary layers along the ﬂuid direction spread in the domain. Dispersing solid nanoparticles in the base ﬂuid with higher volume fractions resulted in better temperature distribution and Nusselt number. By increasing the nanoparticle volume fraction, the temperature of the heated surface penetrated to the ﬂow centerline and the ﬂuid temperature increased. Jet injection with higher velocity, due to its higher ﬂuid momentum, resulted in higher Nusselt number and a ﬀ ected lateral areas. Fluid velocity was higher in jet areas, which diminished the e ﬀ ect of the boundary layer.


Introduction
In recent years, industrial developments have led scientists to search for methods to improve heat transfer in heat exchangers and industrial equipment. Therefore, a new generation of cooling fluids, called nanofluids, is used in industrial and commercial applications. Cooling fluid jet is used in turbine blade cooling and indirect access surfaces. Metal and non-metal nanoparticles have a higher thermal conductivity coefficient than water and lead to higher conductive heat transfer coefficients of fluid, as well as improve the temperature distribution of the nanofluid. Experimental results show that adding nanoparticles to the base fluid increases the heat transfer coefficient of nanofluids. Using nanofluids is one of the novel heat transfer improvement methods with high efficiency [1][2][3][4].
Numerous researchers have investigated the thermal or hydrodynamic performance of nanofluids in different microchannel heat sinks [5][6][7]. Hang et al. [8] numerically investigated the heat transfer performance of a microchannel heat sink with different nanofluids. Akbari et al. [9] studied laminar flow and heat transfer parameters of water/Al 2 O 3 nanofluid with different nanoparticle volume fractions inside a rectangular microchannel and found that using rough surfaces in microchannel leads to higher heat transfer. Behnampour et al. [10] numerically investigated laminar flow and heat transfer parameters of water/AgO nanofluid with different nanoparticle volume fractions in a rectangular microchannel and showed that by increasing fluid velocity, an optimized trade-off can be obtained between heat transfer, hydrodynamic behavior of nanofluid, and the performance evaluation criteria (PEC) variations. Geravandian et al. [11] numerically simulated the laminar heat transfer of nanofluid flow in a rectangular microchannel and revealed that by increasing TiO 2 nanoparticles, heat transfer, friction coefficient, PEC, and pressure drop increase. Studies on the effect of using cooling fluid jet injection on heated surfaces [12] and other methods to increase heat transfer, such as using dimples, rough surfaces [13,14], and twisted tapes [15,16], have been conducted for different industrial and experimental geometries. These studies show that by creating vortexes, uniform temperature distribution can be obtained. Fluid jet plays an important role in cooling technologies and by creating better mixtures of cooling fluid flow, thermal performance can be enhanced [17][18][19]. Chen et al. [20] numerically and experimentally investigated the forced convection heat transfer inside a rectangular channel for determining fluid flow and heat transfer properties. They also compared the performance of heat sinks with solid and perforated pins and showed that by increasing the number of perforations and their diameter, pressure drop decreases and the Nusselt number increases. In their experiments, thermal performance of the heat sink with perforated pins was better than with solid pins. Nafon et al. [21] experimentally studied the effects of inlet temperature, Reynolds number, and heat flux on heat transfer properties of a water/TiO 2 nanofluid jet in a semi-rectangular heat sink, and showed that the average heat transfer coefficient of nanofluid is higher than base fluid, and the pressure drop increases by increasing the nanoparticle volume fraction. Jasperson et al. [22] studied the thermal and hydrodynamic performance of a copper microchannel and a pin fin microchannel and showed that by increasing the volume rate of flow, the thermal resistance of a pin fin heat sink decreases. Zhuwang et al. [23] studied heat transfer inside a microchannel with fluid jet and different coolants and showed that using fluid jet results in higher heat transfer compared with ordinary parallel flows.
Few studies have investigated the thermal and hydrodynamic performance of laminar nanofluid flow in a rectangular microchannel with cooling fluid jet injection. Therefore, the main purpose of this study is to investigate heat transfer improvement methods, such as using nanofluid and direct fluid jet injection. We also investigated the effect of the hydrodynamic velocity boundary condition on the cooling performance of the microchannel. In the present numerical study, laminar flow and heat transfer of oil/MWCNT nanofluid inside a two-dimensional microchannel with nanofluid jet injection were simulated using a finite volume method. Cooling nanofluid jet injection in a microchannel disturbs the thermal boundary layer and increases the heat transfer rate. The use of fluid jet injection in various sections of the microchannel is one of the novelties of this study. Fluid flow and heat transfer behavior are separately simulated in cases with no jets and with one, two, and three jets. The effect of applying slip and no-slip boundary conditions on solid walls of microchannel on the flow and heat transfer were studied.

Problem Definition
In the present numerical study, oil/MWCNT nanofluid flow was simulated in 0-4% nanoparticle volume fractions and Reynolds numbers of 10 and 50. We considered different numbers of fluid jets on the insulated bottom wall of the microchannel. The top heated wall of the microchannel had a constant temperature of T h = 303 K. The inlet cold fluid entered from the left side of the microchannel with the temperature of T c = 293 K. Figure 1 shows the microchannel in the present study. The length of studied microchannel was L = 2500 µm, the width of inlet jet entrance was 10 µm, the height of microchannel was H = 50 µm, and the jet pitch was L1 = 500 µm. The velocity of the inlet nanofluid jet at all studied Reynolds numbers was constant (Re = 10). Fluid flow and heat transfer inside the microchannel were separately simulated in cases with no jet and with 1, 2, and 3 jets. Nanofluids with high velocity through the lower wall were injected into the micro-channel. As the cooling fluid jet flowed on the surface, the temperature of this surface decreased. Also, in this investigation, the effect of applying slip and no-slip boundary conditions on solid walls of the microchannel in slip coefficients of B = 0.0, 0.04, and 0.08 were investigated. The flow was laminar, forced, Newtonian, single-phase, and incompressible; the nanofluid was homogeneous and uniform; and also, the radiation effects were neglected. The properties of the base fluid and nanofluid in different nanoparticle volume fractions are presented in Table 1 [24].  The length of studied microchannel was L = 2500 µm, the width of inlet jet entrance was 10 µm, the height of microchannel was H = 50 µm, and the jet pitch was L 1 = 500 µm. The velocity of the inlet nanofluid jet at all studied Reynolds numbers was constant (Re = 10). Fluid flow and heat transfer inside the microchannel were separately simulated in cases with no jet and with 1, 2, and 3 jets. Nanofluids with high velocity through the lower wall were injected into the micro-channel. As the cooling fluid jet flowed on the surface, the temperature of this surface decreased. Also, in this investigation, the effect of applying slip and no-slip boundary conditions on solid walls of the microchannel in slip coefficients of B = 0.0, 0.04, and 0.08 were investigated. The flow was laminar, forced, Newtonian, single-phase, and incompressible; the nanofluid was homogeneous and uniform; and also, the radiation effects were neglected. The properties of the base fluid and nanofluid in different nanoparticle volume fractions are presented in Table 1 [24].

Numerical Details
Finite volume method [27,28] and the second-order upwind discretization method with maximum residual of 10 −6 [29,30] was used. For coupling velocity-pressure equations, the SIMPLEC [31,32] algorithm was used. Also, the effect of nanoparticle volume fraction, Reynolds number, and jet number on flow and heat transfer parameters were investigated; and flow parameters, temperature, and streamlines contours are the presented.
The meshes in this study were rectangular and regular grids ranging from 25,000 to 100,000 grids were used for Reynolds numbers of 10 and 50, 2% volume fraction, and no-slip boundary conditions (Table 2). In order to ensure the independence of flow and heat transfer parameters from mesh number, Nusselt number and average velocity on flow centerline were studied for different mesh numbers. The 850 × 50 mesh number has better performance than 900 × 90, therefore, the former was used in the present simulations.  Figures 2 and 3 show the validation of velocity and dimensionless flow temperature at the outlet section of microchannel against the results of Raisi et al. [33], who investigated laminar and forced flow of water/CuO nanofluid with different nanoparticle volume fractions in a two-dimensional rectangular microchannel using a finite volume method. The results of Raisi et al. [33] for 3% nanoparticle volume fraction, Re = 100, and different slip velocity coefficients were used for the validation. There is good agreement between the results of the present study and Raisi et al. [33]. Therefore, the boundary conditions and assumptions of this study are correct. mesh number, Nusselt number and average velocity on flow centerline were studied for different mesh numbers. The 850 × 50 mesh number has better performance than 900 × 90, therefore, the former was used in the present simulations.  3 show the validation of velocity and dimensionless flow temperature at the outlet section of microchannel against the results of Raisi et al. [33], who investigated laminar and forced flow of water/CuO nanofluid with different nanoparticle volume fractions in a two-dimensional rectangular microchannel using a finite volume method. The results of Raisi et al. [33] for 3% nanoparticle volume fraction, Re = 100, and different slip velocity coefficients were used for the validation. There is good agreement between the results of the present study and Raisi et al. [33]. Therefore, the boundary conditions and assumptions of this study are correct.

Results and Discussion
In the present numerical study, the effect of nanofluid jet injection on laminar flow and heat transfer of oil/MWCNT nanofluid for 0-4% nanoparticle volume fractions, at Reynolds numbers of 10 and 50, and different slip velocity coefficients was numerically studied. The effects of the variations of jet number, nanoparticle volume fractions, and velocity boundary layer on solid walls of microchannel were also studied. Reynolds number, nanoparticle volume fraction, and slip velocity coefficient were input variables and their effect on Nusselt number, temperature domain, velocity domain, flow parameters, and temperature contours were investigated and compared based on their local and average values.

Streamlines and Isothermal Contours
In Figure 4, streamline contours (right side) and constant temperature contours (left side) are shown for Re = 10, B = 0.04, and φ = 2% for cases 1, 2, 3, and 4. In these contours, the effect of fluid motion with minimum velocity along the microchannel at different nanofluid jet injection on the heated surface was investigated. The temperature difference between inlet cold fluid and the heated wall improves heat transfer and creates temperature gradients. Fluid motion and its impact with the solid walls of the microchannel creates velocity gradients in areas close to the wall. In general, fluid motion along the microchannel leads to the creation of velocity and thermal boundary layers. Nanofluid jet injection on fluid direction results in the elimination of velocity and thermal boundary layers. In all contours, by increasing the number of fluid jets, areas with eliminated boundary layers become larger. Areas with a high number of jet injections have smaller velocity boundary layers. In the cases with more jets, due to the high volume of crossing fluid, axial velocity parameters of fluid are higher than the cases with fewer jets.

Results and Discussion
In the present numerical study, the effect of nanofluid jet injection on laminar flow and heat transfer of oil/MWCNT nanofluid for 0-4% nanoparticle volume fractions, at Reynolds numbers of 10 and 50, and different slip velocity coefficients was numerically studied. The effects of the variations of jet number, nanoparticle volume fractions, and velocity boundary layer on solid walls of microchannel were also studied. Reynolds number, nanoparticle volume fraction, and slip velocity coefficient were input variables and their effect on Nusselt number, temperature domain, velocity domain, flow parameters, and temperature contours were investigated and compared based on their local and average values.

Streamlines and Isothermal Contours
In Figure 4, streamline contours (right side) and constant temperature contours (left side) are shown for Re = 10, B = 0.04, and φ = 2% for cases 1, 2, 3, and 4. In these contours, the effect of fluid motion with minimum velocity along the microchannel at different nanofluid jet injection on the heated surface was investigated. The temperature difference between inlet cold fluid and the heated wall improves heat transfer and creates temperature gradients. Fluid motion and its impact with the solid walls of the microchannel creates velocity gradients in areas close to the wall. In general, fluid motion along the microchannel leads to the creation of velocity and thermal boundary layers. Nanofluid jet injection on fluid direction results in the elimination of velocity and thermal boundary layers. In all contours, by increasing the number of fluid jets, areas with eliminated boundary layers become larger. Areas with a high number of jet injections have smaller velocity boundary layers. In the cases with more jets, due to the high volume of crossing fluid, axial velocity parameters of fluid are higher than the cases with fewer jets. Streamline contours (right side) and constant temperature contours (left side) at Re = 50, B = 0.04, and φ = 2% are presented in Figure 5. The effect of increasing jet number on the heated surface on these contours is shown for cases 1, 2, 3, and 4. In the inlet areas of the microchannel and jet injection areas is 5 times higher, contrary to Figure 4. Increasing fluid velocity decreases the thermal boundary layer thickness. Increasing the fluid momentum results in higher crossing fluid velocity. At Re = 50, by increasing the number of fluid jets, due to the better mixture of fluid and reduction of temperature curve slope, a significant reduction of the thermal boundary layer thickness is observed. According to the streamline contours, nanofluid jet injection with higher velocity results in better mixture of flow and creation of local vortices. By increasing Reynolds number, especially in the last jet, the effect of vertical jet injection decreases due to the higher of volumetric rate of crossing fluid and creation of a smaller vortex behind each jet.

Isothermal line Stream function
Case (1) Streamline contours (right side) and constant temperature contours (left side) at Re = 50, B = 0.04, and φ = 2% are presented in Figure 5. The effect of increasing jet number on the heated surface on these contours is shown for cases 1, 2, 3, and 4. In the inlet areas of the microchannel and jet injection areas is 5 times higher, contrary to Figure 4. Increasing fluid velocity decreases the thermal boundary layer thickness. Increasing the fluid momentum results in higher crossing fluid velocity. At Re = 50, by increasing the number of fluid jets, due to the better mixture of fluid and reduction of temperature curve slope, a significant reduction of the thermal boundary layer thickness is observed. According to the streamline contours, nanofluid jet injection with higher velocity results in better mixture of flow and creation of local vortices. By increasing Reynolds number, especially in the last jet, the effect of vertical jet injection decreases due to the higher of volumetric rate of crossing fluid and creation of a smaller vortex behind each jet. Streamline contours (right side) and constant temperature contours (left side) at Re = 50, B = 0.04, and φ = 2% are presented in Figure 5. The effect of increasing jet number on the heated surface on these contours is shown for cases 1, 2, 3, and 4. In the inlet areas of the microchannel and jet injection areas is 5 times higher, contrary to Figure 4. Increasing fluid velocity decreases the thermal boundary layer thickness. Increasing the fluid momentum results in higher crossing fluid velocity. At Re = 50, by increasing the number of fluid jets, due to the better mixture of fluid and reduction of temperature curve slope, a significant reduction of the thermal boundary layer thickness is observed. According to the streamline contours, nanofluid jet injection with higher velocity results in better mixture of flow and creation of local vortices. By increasing Reynolds number, especially in the last jet, the effect of vertical jet injection decreases due to the higher of volumetric rate of crossing fluid and creation of a smaller vortex behind each jet.

Local Nusselt Number
In Figure 6, the local Nusselt number on the heated wall of the microchannel is shown for Re = 10, B = 0.04, and different nanoparticles volume fractions. In this figure, we have different jet numbers with minimum fluid velocity. Nanoparticle volume fractions were 0-4%. Heat transfer occurs because of temperature differences between fluid and the top heated wall. The value of heat transfer has a direct relationship with Nusselt number and temperature line slope. Uniform temperature distribution and removal of temperature gradients are desirable and lead to higher Nusselt numbers. The local jump in the Nusselt number in Figure 6 was because of the low temperature of the injected fluid jet on the heated surface. In these areas, the heat transfer distribution was more uniform and the Nusselt number increased. In all curves for Figure 6, by adding nanoparticle volume fraction, the Nusselt number and heat transfer were significantly increased.
In Figure 7, the local Nusselt numbers at Re = 50, B = 0.04, and different nanoparticle volume fractions are shown. In this figure, the effect of the increasing fluid velocity on Nusselt number is compared with Figure 6. Increasing ϕ and the number of jets resulted in the higher Nusselt number graphs in Figure 7. Figures 6 and 7 show the elimination of thermal and velocity boundary layers due to the fluid jet injection, which is more significant at higher Reynolds numbers. Jet injection with higher velocity leads to the increase of the Nusselt number due to the higher fluid momentum which affects the lateral regions. Increasing the nanoparticle volume fraction at Re = 50, compared to Re = 10, resulted in a significant rise of the Nusselt number. In general, the local Nusselt number on heated walls showed that using fluid jet injection with low temperature can increase the heat transfer rate. Compared to the primary jets, the heat transfer performance of the last jets was low, because by adding more jets, the effects of fluid momentum and crossing fluid become more significant, which diminishes the lateral flow by vertical jets as they need higher momentum. Therefore, in cases 2 and 3, the lateral cooling areas were narrower.

Local Nusselt Number
In Figure 6, the local Nusselt number on the heated wall of the microchannel is shown for Re = 10, B = 0.04, and different nanoparticles volume fractions. In this figure, we have different jet numbers with minimum fluid velocity. Nanoparticle volume fractions were 0-4%. Heat transfer occurs because of temperature differences between fluid and the top heated wall. The value of heat transfer has a direct relationship with Nusselt number and temperature line slope. Uniform temperature distribution and removal of temperature gradients are desirable and lead to higher Nusselt numbers. The local jump in the Nusselt number in Figure 6 was because of the low temperature of the injected fluid jet on the heated surface. In these areas, the heat transfer distribution was more uniform and the Nusselt number increased. In all curves for Figure 6, by adding nanoparticle volume fraction, the Nusselt number and heat transfer were significantly increased.
In Figure 7, the local Nusselt numbers at Re = 50, B = 0.04, and different nanoparticle volume fractions are shown. In this figure, the effect of the increasing fluid velocity on Nusselt number is compared with Figure 6. Increasing φ and the number of jets resulted in the higher Nusselt number graphs in Figure 7. Figures 6 and 7 show the elimination of thermal and velocity boundary layers due to the fluid jet injection, which is more significant at higher Reynolds numbers. Jet injection with higher velocity leads to the increase of the Nusselt number due to the higher fluid momentum which affects the lateral regions. Increasing the nanoparticle volume fraction at Re = 50, compared to Re = 10, resulted in a significant rise of the Nusselt number. In general, the local Nusselt number on heated walls showed that using fluid jet injection with low temperature can increase the heat transfer rate. Compared to the primary jets, the heat transfer performance of the last jets was low, because by adding more jets, the effects of fluid momentum and crossing fluid become more significant, which diminishes the lateral flow by vertical jets as they need higher momentum. Therefore, in cases 2 and 3, the lateral cooling areas were narrower.

Temperature along the Symmetry Plane
In Figure 8, local dimensionless temperature distribution on flow centerline for Re = 10, B = 0.04, and volume fractions of 0-4% is presented. The dimensionless temperature distribution in the central section of the microchannel was influenced by the reduction of fluid velocity and jet number for cases 1, 2, 3, and 4, and different nanoparticle volume fractions. Because of the temperature difference between the fluid and heated surface, the penetration of the thermal boundary layer affected all areas of the microchannel, especially the areas in the central line of flow. Figure 8 shows that increasing ϕ resulted in higher thermal conductivity of the microchannel and uniform temperature distribution in all areas of the microchannel, especially in the heated areas. Due to the higher thermal conductivity of cooling fluid in higher nanoparticle volume fractions, temperature gradients in heated areas were diminished, and also, the conductivity of the fluid layers was higher. Therefore, by increasing ϕ, heat transfer penetrated the flow centerline, and because of the higher conductivity of the fluid layers, the fluid temperature rose.

Temperature along the Symmetry Plane
In Figure 8, local dimensionless temperature distribution on flow centerline for Re = 10, B = 0.04, and volume fractions of 0-4% is presented. The dimensionless temperature distribution in the central section of the microchannel was influenced by the reduction of fluid velocity and jet number for cases 1, 2, 3, and 4, and different nanoparticle volume fractions. Because of the temperature difference between the fluid and heated surface, the penetration of the thermal boundary layer affected all areas of the microchannel, especially the areas in the central line of flow. Figure 8 shows that increasing φ resulted in higher thermal conductivity of the microchannel and uniform temperature distribution in all areas of the microchannel, especially in the heated areas. Due to the higher thermal conductivity of cooling fluid in higher nanoparticle volume fractions, temperature gradients in heated areas were diminished, and also, the conductivity of the fluid layers was higher. Therefore, by increasing φ, heat transfer penetrated the flow centerline, and because of the higher conductivity of the fluid layers, the fluid temperature rose.

Temperature along the Symmetry Plane
In Figure 8, local dimensionless temperature distribution on flow centerline for Re = 10, B = 0.04, and volume fractions of 0-4% is presented. The dimensionless temperature distribution in the central section of the microchannel was influenced by the reduction of fluid velocity and jet number for cases 1, 2, 3, and 4, and different nanoparticle volume fractions. Because of the temperature difference between the fluid and heated surface, the penetration of the thermal boundary layer affected all areas of the microchannel, especially the areas in the central line of flow. Figure 8 shows that increasing ϕ resulted in higher thermal conductivity of the microchannel and uniform temperature distribution in all areas of the microchannel, especially in the heated areas. Due to the higher thermal conductivity of cooling fluid in higher nanoparticle volume fractions, temperature gradients in heated areas were diminished, and also, the conductivity of the fluid layers was higher. Therefore, by increasing ϕ, heat transfer penetrated the flow centerline, and because of the higher conductivity of the fluid layers, the fluid temperature rose. In Figure 9, local dimensionless temperature distribution at Re = 50, B = 0.04, and 0-4% nanoparticle volume fractions along the flow centerline are shown for cases 1, 2, 3, and 4. Increasing the number of jets resulted in lower temperature of the fluid, and thus the dimensionless temperature graphs in the microchannel centerlines were decreased, which was due to the reduction of temperature associated with increased fluid velocity and transferred temperature to the flow centerline. By increasing the velocity of injected cold fluid jet, the cold and hot fluids were better mixed. Consequently, the fluid temperature decreased significantly. According to Figures 8 and 9, by decreasing the fluid velocity, the fluid becomes thermally better developed, which is due to the penetration of heat to the inlet section of the microchannel. By increasing the fluid velocity, the undeveloped length of flow at the inlet section of microchannel increases. In Figure 9, local dimensionless temperature distribution at Re = 50, B = 0.04, and 0-4% nanoparticle volume fractions along the flow centerline are shown for cases 1, 2, 3, and 4. Increasing the number of jets resulted in lower temperature of the fluid, and thus the dimensionless temperature graphs in the microchannel centerlines were decreased, which was due to the reduction of temperature associated with increased fluid velocity and transferred temperature to the flow centerline. By increasing the velocity of injected cold fluid jet, the cold and hot fluids were better mixed. Consequently, the fluid temperature decreased significantly. According to Figures 8 and 9, by decreasing the fluid velocity, the fluid becomes thermally better developed, which is due to the penetration of heat to the inlet section of the microchannel. By increasing the fluid velocity, the undeveloped length of flow at the inlet section of microchannel increases. In Figure 9, local dimensionless temperature distribution at Re = 50, B = 0.04, and 0-4% nanoparticle volume fractions along the flow centerline are shown for cases 1, 2, 3, and 4. Increasing the number of jets resulted in lower temperature of the fluid, and thus the dimensionless temperature graphs in the microchannel centerlines were decreased, which was due to the reduction of temperature associated with increased fluid velocity and transferred temperature to the flow centerline. By increasing the velocity of injected cold fluid jet, the cold and hot fluids were better mixed. Consequently, the fluid temperature decreased significantly. According to Figures 8 and 9, by decreasing the fluid velocity, the fluid becomes thermally better developed, which is due to the penetration of heat to the inlet section of the microchannel. By increasing the fluid velocity, the undeveloped length of flow at the inlet section of microchannel increases.

Axial Velocity along the Symmetric Plate
In Figure 10, dimensionless axial velocity in the flow centerline at Reynolds numbers of 10 and 50, B = 0.0, and φ=2% were investigated. Fluid motion along the microchannel was affected by velocity boundary layer caused by the solid walls. Contrary to the inlet velocity, fluid velocity was higher in case 1, which was due to the effects of the velocity boundary layer along the microchannel walls. In case 1, fluid velocity in developed areas at Re = 50 was lower than Re = 10. Fluid jet increased the volume rate of the crossing fluid downstream of the jets compared to areas upstream of the jets. Therefore, in cases 1, 2, and 3, fluid velocity increased in areas after the jets. In fact, high fluid velocity decreased the size of the velocity boundary layer. The velocity curves for fluid with Re = 50 were higher than fluid with Re = 10. In general, fluid jet eliminates the velocity boundary layer and prevents local flow development.

Axial Velocity along the Symmetric Plate
In Figure 10, dimensionless axial velocity in the flow centerline at Reynolds numbers of 10 and 50, B = 0.0, and φ=2% were investigated. Fluid motion along the microchannel was affected by velocity boundary layer caused by the solid walls. Contrary to the inlet velocity, fluid velocity was higher in case 1, which was due to the effects of the velocity boundary layer along the microchannel walls. In case 1, fluid velocity in developed areas at Re = 50 was lower than Re = 10. Fluid jet increased the volume rate of the crossing fluid downstream of the jets compared to areas upstream of the jets. Therefore, in cases 1, 2, and 3, fluid velocity increased in areas after the jets. In fact, high fluid velocity decreased the size of the velocity boundary layer. The velocity curves for fluid with Re = 50 were higher than fluid with Re = 10. In general, fluid jet eliminates the velocity boundary layer and prevents local flow development.

Axial Velocity along the Symmetric Plate
In Figure 10, dimensionless axial velocity in the flow centerline at Reynolds numbers of 10 and 50, B = 0.0, and φ=2% were investigated. Fluid motion along the microchannel was affected by velocity boundary layer caused by the solid walls. Contrary to the inlet velocity, fluid velocity was higher in case 1, which was due to the effects of the velocity boundary layer along the microchannel walls. In case 1, fluid velocity in developed areas at Re = 50 was lower than Re = 10. Fluid jet increased the volume rate of the crossing fluid downstream of the jets compared to areas upstream of the jets. Therefore, in cases 1, 2, and 3, fluid velocity increased in areas after the jets. In fact, high fluid velocity decreased the size of the velocity boundary layer. The velocity curves for fluid with Re = 50 were higher than fluid with Re = 10. In general, fluid jet eliminates the velocity boundary layer and prevents local flow development.
Case(1) In Figure 11, dimensionless axial velocity along the flow centerline at Reynolds numbers of 10 and 50, B = 0.08, and φ = 2% are shown. Unlike Figure 10, the slip boundary condition was applied to the solid walls which facilitated fluid motion and reduced the velocity boundary layer on the solid walls of the microchannel. Therefore, the velocity profile in the flow centerline for all cases of 1, 2, 3, and 4 were lower than Figure 10. According to Figures 10 and 11, increasing the slip velocity coefficient on the solid walls of the microchannel at Re = 10 had less effect than Re = 50. At Re = 50, the lateral mixture of flow covered more areas in the case with slip boundary conditions (BCs) compared with the no-slip BC case. According to Figure 11, this factor resulted in the microchannel contraction which increased the effect of jet injection areas and local axial velocity. In Figure 11, dimensionless axial velocity along the flow centerline at Reynolds numbers of 10 and 50, B = 0.08, and φ = 2% are shown. Unlike Figure 10, the slip boundary condition was applied to the solid walls which facilitated fluid motion and reduced the velocity boundary layer on the solid walls of the microchannel. Therefore, the velocity profile in the flow centerline for all cases of 1, 2, 3, and 4 were lower than Figure 10. According to Figures 10 and 11, increasing the slip velocity coefficient on the solid walls of the microchannel at Re = 10 had less effect than Re = 50. At Re = 50, the lateral mixture of flow covered more areas in the case with slip boundary conditions (BCs) compared with the no-slip BC case. According to Figure 11, this factor resulted in the microchannel contraction which increased the effect of jet injection areas and local axial velocity. In Figure 11, dimensionless axial velocity along the flow centerline at Reynolds numbers of 10 and 50, B = 0.08, and φ = 2% are shown. Unlike Figure 10, the slip boundary condition was applied to the solid walls which facilitated fluid motion and reduced the velocity boundary layer on the solid walls of the microchannel. Therefore, the velocity profile in the flow centerline for all cases of 1, 2, 3, and 4 were lower than Figure 10. According to Figures 10 and 11, increasing the slip velocity coefficient on the solid walls of the microchannel at Re = 10 had less effect than Re = 50. At Re = 50, the lateral mixture of flow covered more areas in the case with slip boundary conditions (BCs) compared with the no-slip BC case. According to Figure 11, this factor resulted in the microchannel contraction which increased the effect of jet injection areas and local axial velocity.

Average Nusselt Number on the Heated Surface
In Figure 12, the average Nusselt number at Reynolds numbers of 10 and 50, volume fractions of 0, 2, and 4%, and B = 0.04 are shown. Using fluid jet injection with low temperature and adding solid nanoparticles to the base fluid, heat transfer was improved. Increasing the number of jets resulted in a better mixture of cold and hot flows in the microchannel, reduction of total temperature of the microchannel, and better cooling performance of the nanofluid. Adding solid nanoparticles to the base fluid with higher volume fractions resulted in better temperature distribution, and the Nusselt number increased as well. Also, increasing the Reynolds number resulted in higher fluid velocity and better temperature distribution, especially in heated areas, and also, the average Nusselt number increased. Among investigated cases, case 4 had the highest value of Nusselt number due to the maximum number of cold fluid jet injection. In case 1, due to the lack of nanofluid jet injection, the effects of thermal and velocity boundary layers became significant, and this case had the lowest heat transfer and Nusselt number.

Average Nusselt Number on the Heated Surface
In Figure 12, the average Nusselt number at Reynolds numbers of 10 and 50, volume fractions of 0, 2, and 4%, and B = 0.04 are shown. Using fluid jet injection with low temperature and adding solid nanoparticles to the base fluid, heat transfer was improved. Increasing the number of jets resulted in a better mixture of cold and hot flows in the microchannel, reduction of total temperature of the microchannel, and better cooling performance of the nanofluid. Adding solid nanoparticles to the base fluid with higher volume fractions resulted in better temperature distribution, and the Nusselt number increased as well. Also, increasing the Reynolds number resulted in higher fluid velocity and better temperature distribution, especially in heated areas, and also, the average Nusselt number increased. Among investigated cases, case 4 had the highest value of Nusselt number due to the maximum number of cold fluid jet injection. In case 1, due to the lack of nanofluid jet injection, the effects of thermal and velocity boundary layers became significant, and this case had the lowest heat transfer and Nusselt number.

Average Nusselt Number on the Heated Surface
In Figure 12, the average Nusselt number at Reynolds numbers of 10 and 50, volume fractions of 0, 2, and 4%, and B = 0.04 are shown. Using fluid jet injection with low temperature and adding solid nanoparticles to the base fluid, heat transfer was improved. Increasing the number of jets resulted in a better mixture of cold and hot flows in the microchannel, reduction of total temperature of the microchannel, and better cooling performance of the nanofluid. Adding solid nanoparticles to the base fluid with higher volume fractions resulted in better temperature distribution, and the Nusselt number increased as well. Also, increasing the Reynolds number resulted in higher fluid velocity and better temperature distribution, especially in heated areas, and also, the average Nusselt number increased. Among investigated cases, case 4 had the highest value of Nusselt number due to the maximum number of cold fluid jet injection. In case 1, due to the lack of nanofluid jet injection, the effects of thermal and velocity boundary layers became significant, and this case had the lowest heat transfer and Nusselt number.

Effect of Slip Coefficient on Axial Velocity
Dimensionless axial velocity in the central line of the microchannel at Re = φ = 2% and different slip coefficients was demonstrated in Figure 13. As it is seen, nanofluid jet injection at Re = 10 leads to some changes in jet areas. At this Reynolds number, increasing slip velocity coefficient results in lower fluid momentum dissipation and higher axial velocity. By increasing the slip velocity coefficient, because of the smaller velocity boundary layer, fluid moves on the microchannel walls with less friction. By increasing the jet number with a higher fluid velocity and volumetric rate of crossing fluid, the heat transfer was increased.

Effect of Slip Coefficient on Axial Velocity
Dimensionless axial velocity in the central line of the microchannel at Re = φ = 2% and different slip coefficients was demonstrated in Figure 13. As it is seen, nanofluid jet injection at Re = 10 leads to some changes in jet areas. At this Reynolds number, increasing slip velocity coefficient results in lower fluid momentum dissipation and higher axial velocity. By increasing the slip velocity coefficient, because of the smaller velocity boundary layer, fluid moves on the microchannel walls with less friction. By increasing the jet number with a higher fluid velocity and volumetric rate of crossing fluid, the heat transfer was increased. In Figure 14, dimensionless axial velocity in the flow centerline at Re = 50 and φ = 2% for cases 1, 2, 3, and 4 are shown. Increasing the slip coefficient on the solid surfaces parallel with the fluid direction, especially at higher Reynolds numbers, leads to the higher of fluid momentum. By increasing the slip velocity coefficient, due to the local flow contraction and significant changes in axial fluid velocity, fluid jet injection creates larger areas in the longitudinal direction. In Figure 14, because of the higher value of Reynolds number, the level of graphs are higher than those in Figure 13.  In Figure 14, dimensionless axial velocity in the flow centerline at Re = 50 and φ = 2% for cases 1, 2, 3, and 4 are shown. Increasing the slip coefficient on the solid surfaces parallel with the fluid direction, especially at higher Reynolds numbers, leads to the higher of fluid momentum. By increasing the slip velocity coefficient, due to the local flow contraction and significant changes in axial fluid velocity, fluid jet injection creates larger areas in the longitudinal direction. In Figure 14, because of the higher value of Reynolds number, the level of graphs are higher than those in Figure 13. In Figure 15, dimensionless axial velocity profiles at the outlet section of the microchannel at Re = 10, φ = 4%, and slip coefficients of 0.0 and 0.08 are shown. Applying the slip coefficient on the solid walls resulted in a slight reduction of the axial velocity. By increasing the slip velocity coefficient, the velocity increased and the size of the velocity boundary layer on the heated surface decreased. According to the axial velocity profiles with the no-slip boundary condition, the effects of the velocity boundary layer was significant for the layers further away from the solid walls, and the velocity profile was stronger in the flow centerline. Therefore, from solid walls to flow centerline, the slope of the axial velocity variations was higher than the velocity profile with B = 0.08. Applying the slip velocity coefficient on the solid walls resulted in lower axial velocity gradients. The slope of the velocity curve from the walls to the flow centerline indicated lower gradients than the case with the no-slip boundary condition. The velocity profile in the central areas of the flow confirmed this issue and the corresponding curve had a milder slope. In Figure 15, dimensionless axial velocity profiles at the outlet section of the microchannel at Re = 10, φ = 4%, and slip coefficients of 0.0 and 0.08 are shown. Applying the slip coefficient on the solid walls resulted in a slight reduction of the axial velocity. By increasing the slip velocity coefficient, the velocity increased and the size of the velocity boundary layer on the heated surface decreased. According to the axial velocity profiles with the no-slip boundary condition, the effects of the velocity boundary layer was significant for the layers further away from the solid walls, and the velocity profile was stronger in the flow centerline. Therefore, from solid walls to flow centerline, the slope of the axial velocity variations was higher than the velocity profile with B = 0.08. Applying the slip velocity coefficient on the solid walls resulted in lower axial velocity gradients. The slope of the velocity curve from the walls to the flow centerline indicated lower gradients than the case with the no-slip boundary condition. The velocity profile in the central areas of the flow confirmed this issue and the corresponding curve had a milder slope.

Conclusions
In the present numerical study, laminar flow of oil/MWCNT nanofluid inside a two-dimensional microchannel with nanofluid jet injection was investigated. This simulation was carried out in a two-dimensional domain for laminar and Newtonian flow and heat transfer of oil/MWCNT nanofluid, with volume fractions of 0-4%, at Reynolds numbers of 10 and 50. Various jet numbers were used on the insulated bottom wall of the microchannel. Results showed that nanofluid jet injection with higher velocity leads to better flow mixture through the formation of local vortexes. In Re = 50, by increasing the number of fluid jets, the thermal boundary layer thickness becomes smaller because of the better mixture of fluid. Also, the temperature line slope was milder. Heat transfer occurs due to the temperature difference between fluid and the top heated wall. Heat transfer had a direct relationship with Nusselt number and the slope of the temperature line. Fluid jets eliminated the thermal and velocity boundary layers, which was more significant in higher Reynolds numbers. Because of the temperature differences between fluid and heated surfaces, penetration of the thermal boundary layer affected all sections of the microchannel, especially areas on the flow centerline. Because of higher conductivity of the fluid layer, by increasing the nanoparticle volume fraction, the heat transfer reached the flow centerline and the fluid temperature increased in that area. By decreasing fluid velocity, due to the higher temperature

Conclusions
In the present numerical study, laminar flow of oil/MWCNT nanofluid inside a two-dimensional microchannel with nanofluid jet injection was investigated. This simulation was carried out in a two-dimensional domain for laminar and Newtonian flow and heat transfer of oil/MWCNT nanofluid, with volume fractions of 0-4%, at Reynolds numbers of 10 and 50. Various jet numbers were used on the insulated bottom wall of the microchannel. Results showed that nanofluid jet injection with higher velocity leads to better flow mixture through the formation of local vortexes. In Re = 50, by increasing the number of fluid jets, the thermal boundary layer thickness becomes smaller because of the better mixture of fluid. Also, the temperature line slope was milder. Heat transfer occurs due to the temperature difference between fluid and the top heated wall. Heat transfer had a direct relationship with Nusselt number and the slope of the temperature line. Fluid jets eliminated the thermal and velocity boundary layers, which was more significant in higher Reynolds numbers. Because of the temperature differences between fluid and heated surfaces, penetration of the thermal boundary layer affected all sections of the microchannel, especially areas on the flow centerline. Because of higher conductivity of the fluid layer, by increasing the nanoparticle volume fraction, the heat transfer reached the flow centerline and the fluid temperature increased in that area. By decreasing fluid velocity, due to the higher temperature penetration, the developed length of flow became shorter and the fluid developed faster. Higher fluid velocity in areas with jets eliminated the velocity boundary layer.
Unlike the case with the no-slip boundary condition, by applying the slip coefficient on the solid walls of the microchannel, the lateral mixture of flow at Re = 50 covered larger areas. Increasing the Reynolds number resulted in higher fluid velocity and better temperature distribution, especially in heated areas, and this led to higher average Nusselt number. In fluid with higher velocity, applying the slip boundary condition on solid walls significantly increases the crossing fluid momentum.
Author Contributions: All authors contributed equally.