Mathematical Modeling and Computer Simulations of Nanofluid Flow with Applications to Cooling and Lubrication

There is a growing range of applications of nanoparticle-suspension flows with or without heat transfer. Examples include enhanced cooling of microsystems with low volume-fractions of nanoparticles in liquids, improved tribological performance with lubricants seeded with nanoparticles, optimal nanodrug delivery in the pulmonary as well as the vascular systems to combat cancer, and spray-coating using plasma-jets with seeded nanoparticles. In order to implement theories that explain experimental evidence of nanoparticle-fluid dynamics and predict numerically optimum system performance, a description of the basic math modeling and computer simulation aspects is necessary. Thus, in this review article, the focus is on the fundamental understanding of the physics of nanofluid flow and heat transfer with summaries of microchannel-flow applications related to cooling and lubrication.


Introduction
Over the last several decades, there has been an ever-increasing interest in nanofluids, i.e., liquids with all sorts of solid nanoparticles (NPs) well dispersed at low concentrations. Application areas for nanofluids range from engineering to medicine, taking advantage of their unique properties in heat transfer [1][2][3], drug delivery [4][5][6][7], mass transport [8,9], boiling phenomena [10,11], absorption and radiation [12][13][14][15], optics [16,17], reacting surfaces and catalysts [18], spray-coating [19,20], and lubrication [21][22][23][24]. Though initially investigated in the heat transfer community (see [25]), the concept of "nanofluids" is being continuously expanded in developing more and more useful applications. Three application-schematics are given in Figure 1. Clearly, in order to achieve optimal system performance in present and future applications, the underlying physics has to be clearly understood via benchmark experiments as well as sound theories. Focusing here on various aspects of mathematical modeling and computer simulations, transport processes of nanofluids are analyzed and system performances numerically predicted.

Experimental Evidence
Experimental evidence indicating the advantages of using nanofluids for cooling and lubrication is abundant. In this section, we discuss the measured performance enhancement when using nanofluids in cooling and lubrication in contrast to employing base fluids alone.

Cooling with Nanofluids
One of the main potential applications of nanofluids is enhanced heat transfer, especially the cooling of microsystems with very large heat fluxes. This has been the focus of a number of studies, where microchannel heat sinks (MCHS's) were chosen as heat removal devices. Convection heat transfer in microchannels using conventional fluids such as water and ethylene glycol have been recently reviewed in [27,28]. Nanofluids promise better performance for heat transfer in microchannels than the base liquids, provided that all system components are carefully designed. Largely supportive experimental evidence of improved heat transfer performance when using nanofluids in microsystems is summarized in chronological order as follows.
Chein and Chuang [29] measured the heat transfer and fluid flow characteristics of 0.2-0.4 vol% CuO-water nanofluids in a trapezoidal silicon MCHS. They observed a large enhancement of heat transfer at low flow rates, but little benefit at high flow rates when compared to the base fluid. They attributed this result to the fact that at high flow rates heat is transferred rapidly, whereas nanoparticles (NPs) contribute little to the thermal exchange. While only a slight pressure drop increase was observed in the presence of NPs, nanoparticle deposition was recorded. Jung et al. [30] measured the heat transfer and pressure drop of Al2O3-water nanofluid flow in rectangular microchannels in the laminar flow regime with NP-volume fractions up to 1.8%. They found that the Nusselt number of the nanofluid flow can be correlated by Nu 0.014 Re Pr  . The measured increase in convection heat transfer coefficient at 1.8% NP volume fraction was 32% over that of the distilled water, while the measured friction losses were similar for both the nanofluid and the base fluid. Ho et al. [31] experimentally assessed the thermal and hydraulic performance of Al2O3-water nanofluid flow in a copper MCHS with rectangular microchannels in the laminar flow regime. It was found that the friction factor in the microchannel heat sink cooled by the nanofluid, containing up to 2 vol% NPs, slightly increased relative to the pure water. On the other hand, a 70% increase in average heat transfer coefficient and 25% decrease in MCHS thermal resistance were found by adding 1 vol% NPs to the base fluid at Re = 1676. Using nanofluids was found to be more beneficial at higher flow rates. In contrast, the experimental study of Anoop et al. [32] showed that, when compared to water, nanofluid heat transfer enhancements occurred at lower flow rates in a microchannel, and heat transfer degradation occurred at higher flow rates due to NP deposition on the microchannel wall. Byrne et al. [33] measured the fluid flow and heat transfer characteristics of CuO-water nanofluids in microchannels with NP volume fractions from 0.005% to 0.1%. While a modest improvement in heat transfer was found, the use of a surfactant was essential in maintaining a proper suspension of NPs in the liquid. Rimbault et al. [34] evaluated transitional and turbulent CuO-water flow in a MCHS

Experimental Evidence
Experimental evidence indicating the advantages of using nanofluids for cooling and lubrication is abundant. In this section, we discuss the measured performance enhancement when using nanofluids in cooling and lubrication in contrast to employing base fluids alone.

Cooling with Nanofluids
One of the main potential applications of nanofluids is enhanced heat transfer, especially the cooling of microsystems with very large heat fluxes. This has been the focus of a number of studies, where microchannel heat sinks (MCHS's) were chosen as heat removal devices. Convection heat transfer in microchannels using conventional fluids such as water and ethylene glycol have been recently reviewed in [27,28]. Nanofluids promise better performance for heat transfer in microchannels than the base liquids, provided that all system components are carefully designed. Largely supportive experimental evidence of improved heat transfer performance when using nanofluids in microsystems is summarized in chronological order as follows.
Chein and Chuang [29] measured the heat transfer and fluid flow characteristics of 0.2-0.4 vol % CuO-water nanofluids in a trapezoidal silicon MCHS. They observed a large enhancement of heat transfer at low flow rates, but little benefit at high flow rates when compared to the base fluid. They attributed this result to the fact that at high flow rates heat is transferred rapidly, whereas nanoparticles (NPs) contribute little to the thermal exchange. While only a slight pressure drop increase was observed in the presence of NPs, nanoparticle deposition was recorded. Jung et al. [30] measured the heat transfer and pressure drop of Al 2 O 3 -water nanofluid flow in rectangular microchannels in the laminar flow regime with NP-volume fractions up to 1.8%. They found that the Nusselt number of the nanofluid flow can be correlated by Nu " 0.014ϕ 0.095 Re 0.4 Pr 0.6 . The measured increase in convection heat transfer coefficient at 1.8% NP volume fraction was 32% over that of the distilled water, while the measured friction losses were similar for both the nanofluid and the base fluid. Ho et al. [31] experimentally assessed the thermal and hydraulic performance of Al 2 O 3 -water nanofluid flow in a copper MCHS with rectangular microchannels in the laminar flow regime. It was found that the friction factor in the microchannel heat sink cooled by the nanofluid, containing up to 2 vol % NPs, slightly increased relative to the pure water. On the other hand, a 70% increase in average heat transfer coefficient and 25% decrease in MCHS thermal resistance were found by adding 1 vol % NPs to the base fluid at Re = 1676. Using nanofluids was found to be more beneficial at higher flow rates. In contrast, the experimental study of Anoop et al. [32] showed that, when compared to water, nanofluid heat transfer enhancements occurred at lower flow rates in a microchannel, and heat transfer degradation occurred at higher flow rates due to NP deposition on the microchannel wall. Byrne et al. [33] measured the fluid flow and heat transfer characteristics of CuO-water nanofluids in microchannels with NP volume fractions from 0.005% to 0.1%. While a modest improvement in heat transfer was found, the use of a  [34] evaluated transitional and turbulent CuO-water flow in a MCHS with NP-diameters of 29 nm and volume fractions ranging from 0.24% to 4.5%. Little heat transfer enhancement was obtained by adding NPs at low concentrations, while a clear decrease was found at high concentrations. On the other hand, a large pressure drop increase (up to 70% at 4.5 vol %) with respect to water was found for all NP-loading ranges. This is the only study presenting results of transitional and turbulent nanofluid flow in microsystems, which is opposite to the beneficial findings for thermal nanofluid flow in the laminar regime. Indeed, the observed deterioration in heat transfer was confirmed by the authors of [35], who analyzed Al 2 O 3 -water flow in a mini-channel. Rimbault et al. [34] suggested that particle migration towards the wall and enhanced local particle deposition may be the reasons. Manay and Sahin [36] experimentally determined the effect of the aspect ratio of microchannels on nanofluid heat transfer. While nanofluids provided better heat transfer performance without causing excessive pressure drop, reduced heat transfer and increased pressure drop were generated by nanofluid flow in microchannels with high aspect ratios. Azizi et al. [37] studied the thermal performance and friction factor of a cylindrical MCHS cooled by Cu-water nanofluids. The results showed that the presence of NPs enhanced the entrance Nusselt number up to 43%; however, the friction factor also increased up to 45.5% when compared to the base fluid.
Some of the experimental evidence of the heat transfer performance of nanofluids in microsystems is shown in Figures 2-4. It can be seen that, based on available data, nanofluids significantly enhance the heat transfer performance in microsystems as compared with that of the base fluids at elevated Reynolds numbers, i.e., Re > 100. Yet in the low Re-number region, both enhanced and reduced heat transfer were measured. Severe fouling in microchannels was suggested to be the reason for heat transfer deterioration. Little data exist for turbulent nanofluid flow in microsystems. Concerning friction characteristics of nanofluid flow in microsystems, the data are more scattered. Most of the studies concluded that a minor pressure drop increase was detected when adding NPs to the base fluid. However, Chevalier et al. [38] measured the rheological properties of nanofluids flowing through microchannels and obtained much higher effective viscosities than those predicted by the Einstein equation (see Equation (26) in Section 3). This suggests that the thermal performance of nanofluids may be compromised. with NP-diameters of 29 nm and volume fractions ranging from 0.24% to 4.5%. Little heat transfer enhancement was obtained by adding NPs at low concentrations, while a clear decrease was found at high concentrations. On the other hand, a large pressure drop increase (up to 70% at 4.5 vol %) with respect to water was found for all NP-loading ranges. This is the only study presenting results of transitional and turbulent nanofluid flow in microsystems, which is opposite to the beneficial findings for thermal nanofluid flow in the laminar regime. Indeed, the observed deterioration in heat transfer was confirmed by the authors of [35], who analyzed Al2O3-water flow in a mini-channel.
Rimbault et al. [34] suggested that particle migration towards the wall and enhanced local particle deposition may be the reasons. Manay and Sahin [36] experimentally determined the effect of the aspect ratio of microchannels on nanofluid heat transfer. While nanofluids provided better heat transfer performance without causing excessive pressure drop, reduced heat transfer and increased pressure drop were generated by nanofluid flow in microchannels with high aspect ratios. Azizi et al. [37] studied the thermal performance and friction factor of a cylindrical MCHS cooled by Cu-water nanofluids. The results showed that the presence of NPs enhanced the entrance Nusselt number up to 43%; however, the friction factor also increased up to 45.5% when compared to the base fluid. Some of the experimental evidence of the heat transfer performance of nanofluids in microsystems is shown in Figures 2-4. It can be seen that, based on available data, nanofluids significantly enhance the heat transfer performance in microsystems as compared with that of the base fluids at elevated Reynolds numbers, i.e., Re > 100. Yet in the low Re-number region, both enhanced and reduced heat transfer were measured. Severe fouling in microchannels was suggested to be the reason for heat transfer deterioration. Little data exist for turbulent nanofluid flow in microsystems. Concerning friction characteristics of nanofluid flow in microsystems, the data are more scattered. Most of the studies concluded that a minor pressure drop increase was detected when adding NPs to the base fluid. However, Chevalier et al. [38] measured the rheological properties of nanofluids flowing through microchannels and obtained much higher effective viscosities than those predicted by the Einstein equation (see Equation (26) in Section 3). This suggests that the thermal performance of nanofluids may be compromised.

Nanoparticles in Lubrication
Friction is one of the main causes of energy loss in mechanical systems. In fact, machine elements usually fail due to excessive wear. Hence, it is essential to provide superior lubrication to improve the efficiency of energy utilization and the reliability of mechanical systems. Adding nanoparticles to the lubricants can enhance their tribological properties, especially in the boundary and mixed lubrication regimes. The anti-friction and anti-wear behaviors of a variety of nanoparticles used as lubricant additives have been reported [39], including metals (e.g., Zn, Cu, Ni, Fe, Co, Pd, and Au), metal oxides (e.g., TiO2, Al2O3, CuO, ZnO, ZrO2, and ZnAl2O4) and their composites, sulfides (e.g., WS2, MoS2, CuS, and ZnS), nonmetals (e.g., diamond, SiO2, carbon nanotubes, and fullerenes), inorganic oxythiomolybdates (e.g., Cs2MoO2S2, ZnMoO2S2), and baron-based nanoparticles (e.g., calcium borate, zinc borate, boron nitride, and alkali borates). Most of the evidence indicates that adding NPs to the lubricants can significantly improve lubrication performance by reducing both friction and wear, thereby allowing for an increase in load-carrying capacity [21,[40][41][42][43][44][45][46][47][48][49][50][51]. A number of mechanisms have been proposed to explain the superior performance of lubricants with nanoparticle additives, including surface adsorption, penetration into asperities, and tribo-chemical reaction to

Nanoparticles in Lubrication
Friction is one of the main causes of energy loss in mechanical systems. In fact, machine elements usually fail due to excessive wear. Hence, it is essential to provide superior lubrication to improve the efficiency of energy utilization and the reliability of mechanical systems. Adding nanoparticles to the lubricants can enhance their tribological properties, especially in the boundary and mixed lubrication regimes. The anti-friction and anti-wear behaviors of a variety of nanoparticles used as lubricant additives have been reported [39], including metals (e.g., Zn, Cu, Ni, Fe, Co, Pd, and Au), metal oxides (e.g., TiO2, Al2O3, CuO, ZnO, ZrO2, and ZnAl2O4) and their composites, sulfides (e.g., WS2, MoS2, CuS, and ZnS), nonmetals (e.g., diamond, SiO2, carbon nanotubes, and fullerenes), inorganic oxythiomolybdates (e.g., Cs2MoO2S2, ZnMoO2S2), and baron-based nanoparticles (e.g., calcium borate, zinc borate, boron nitride, and alkali borates). Most of the evidence indicates that adding NPs to the lubricants can significantly improve lubrication performance by reducing both friction and wear, thereby allowing for an increase in load-carrying capacity [21,[40][41][42][43][44][45][46][47][48][49][50][51]. A number of mechanisms have been proposed to explain the superior performance of lubricants with nanoparticle additives, including surface adsorption, penetration into asperities, and tribo-chemical reaction to

Nanoparticles in Lubrication
Friction is one of the main causes of energy loss in mechanical systems. In fact, machine elements usually fail due to excessive wear. Hence, it is essential to provide superior lubrication to improve the efficiency of energy utilization and the reliability of mechanical systems. Adding nanoparticles to the lubricants can enhance their tribological properties, especially in the boundary and mixed lubrication regimes. The anti-friction and anti-wear behaviors of a variety of nanoparticles used as lubricant additives have been reported [39], including metals (e.g., Zn, Cu, Ni, Fe, Co, Pd, and Au), metal oxides (e.g., TiO 2 , Al 2 O 3 , CuO, ZnO, ZrO 2 , and ZnAl 2 O 4 ) and their composites, sulfides (e.g., WS 2 , MoS 2 , CuS, and ZnS), nonmetals (e.g., diamond, SiO 2 , carbon nanotubes, and fullerenes), inorganic oxythiomolybdates (e.g., Cs 2 MoO 2 S 2 , ZnMoO 2 S 2 ), and baron-based nanoparticles (e.g., calcium borate, zinc borate, boron nitride, and alkali borates). Most of the evidence indicates that adding NPs to the lubricants can significantly improve lubrication performance by reducing both friction and wear, thereby allowing for an increase in load-carrying capacity [21,[40][41][42][43][44][45][46][47][48][49][50][51]. A number of mechanisms have been proposed to explain the superior performance of lubricants with nanoparticle additives, including surface adsorption, penetration into asperities, and tribo-chemical reaction to reduce wear, as well  reduce wear, as well as size effect, colloidal effect, exfoliation, protective film, and third-body effect to reduce friction (see Figures 5 and 6).  Exfoliation can also happen on nano-sheets from aggregates, from smearing the nanotubes, and by breaking large NPs into smaller pieces (see [52]). The arrows indicate the sliding direction.
The focus of this review is on the enhanced lubrication performance when using nanoparticles and the underlying mechanisms. The different factors that affect the tribological properties as well as different applications are discussed.
Early attempts in using small particles to enhance lubricants were made decades ago. For example, Bartz [53] investigated the influence of the particle-diameter range of 0.7-7 µ m on the lubricating effectiveness of MoS2-suspensions using a four-ball tester. The author found that finer particles only improved wear performance at high loads with smooth surfaces. He suggested that temperature, material, and surface roughness must be taken into account in order to transfer these results to conditions different from the four-ball tester. Cusano and Sliney [54,55] added sub-micron graphite and MoS2 particles to a mineral oil and found deleterious effects under boundary lubrication conditions of a Hertzian contact.
Tao et al. [21] added 1% diamond nanoparticles with a 5-nm average diameter to paraffin oil. The results indicated that, under boundary-lubrication conditions, the nanofluid lubricant exhibited excellent load-carrying capacity, anti-wear, and friction-reduction properties. The authors postulated that the nanoparticles lodge into the rubbing surfaces due to their small size and in effect polish the surface. At the same time, some nanoparticles may form an inlaid pattern into the surface, which raises the hardness of the rubbing surface. Under extended rubbing time, the nanofluid forms a boundary lubricating film between the gliding surfaces. The NPs in the film not only bear the loading but also separate the solid parts and hence prevent their direct contact. They can also roll between the polished surfaces. This ball-bearing effect of nanoparticles also increases the load-bearing capacity while reducing wear and friction.
Subsequent studies [40][41][42][43] confirmed that adding nanoparticles to paraffin and water at low concentrations greatly improved the load-carrying capacity as well as the anti-wear and frictionreduction properties. A boundary protective film, consisting of nanoparticles, was found to form between the gliding surfaces [56,57]. Qiu et al. [58] suggested that the forming of a third body due to the deposition of nanoparticles on the contacting surfaces led to the observed lubrication behavior, which was supported in [44,45]. Some studies supported the mending effect, i.e., the healing of the surface during the tribological test, which may be due to exfoliation of NPs [46]. Nanoparticles fill the dimples and somewhat smooth the asperities, with a tribo-chemical reaction on the friction  reduce wear, as well as size effect, colloidal effect, exfoliation, protective film, and third-body effect to reduce friction (see Figures 5 and 6).  Exfoliation can also happen on nano-sheets from aggregates, from smearing the nanotubes, and by breaking large NPs into smaller pieces (see [52]). The arrows indicate the sliding direction.
The focus of this review is on the enhanced lubrication performance when using nanoparticles and the underlying mechanisms. The different factors that affect the tribological properties as well as different applications are discussed.
Early attempts in using small particles to enhance lubricants were made decades ago. For example, Bartz [53] investigated the influence of the particle-diameter range of 0.7-7 µ m on the lubricating effectiveness of MoS2-suspensions using a four-ball tester. The author found that finer particles only improved wear performance at high loads with smooth surfaces. He suggested that temperature, material, and surface roughness must be taken into account in order to transfer these results to conditions different from the four-ball tester. Cusano and Sliney [54,55] added sub-micron graphite and MoS2 particles to a mineral oil and found deleterious effects under boundary lubrication conditions of a Hertzian contact.
Tao et al. [21] added 1% diamond nanoparticles with a 5-nm average diameter to paraffin oil. The results indicated that, under boundary-lubrication conditions, the nanofluid lubricant exhibited excellent load-carrying capacity, anti-wear, and friction-reduction properties. The authors postulated that the nanoparticles lodge into the rubbing surfaces due to their small size and in effect polish the surface. At the same time, some nanoparticles may form an inlaid pattern into the surface, which raises the hardness of the rubbing surface. Under extended rubbing time, the nanofluid forms a boundary lubricating film between the gliding surfaces. The NPs in the film not only bear the loading but also separate the solid parts and hence prevent their direct contact. They can also roll between the polished surfaces. This ball-bearing effect of nanoparticles also increases the load-bearing capacity while reducing wear and friction.
Subsequent studies [40][41][42][43] confirmed that adding nanoparticles to paraffin and water at low concentrations greatly improved the load-carrying capacity as well as the anti-wear and frictionreduction properties. A boundary protective film, consisting of nanoparticles, was found to form between the gliding surfaces [56,57]. Qiu et al. [58] suggested that the forming of a third body due to the deposition of nanoparticles on the contacting surfaces led to the observed lubrication behavior, which was supported in [44,45]. Some studies supported the mending effect, i.e., the healing of the surface during the tribological test, which may be due to exfoliation of NPs [46]. Nanoparticles fill the dimples and somewhat smooth the asperities, with a tribo-chemical reaction on the friction Figure 6. Schematic presentation of the exfoliation process: individual nanotube exfoliation. Exfoliation can also happen on nano-sheets from aggregates, from smearing the nanotubes, and by breaking large NPs into smaller pieces (see [52]). The arrows indicate the sliding direction.
The focus of this review is on the enhanced lubrication performance when using nanoparticles and the underlying mechanisms. The different factors that affect the tribological properties as well as different applications are discussed.
Early attempts in using small particles to enhance lubricants were made decades ago. For example, Bartz [53] investigated the influence of the particle-diameter range of 0.7-7 µm on the lubricating effectiveness of MoS 2 -suspensions using a four-ball tester. The author found that finer particles only improved wear performance at high loads with smooth surfaces. He suggested that temperature, material, and surface roughness must be taken into account in order to transfer these results to conditions different from the four-ball tester. Cusano and Sliney [54,55] added sub-micron graphite and MoS 2 particles to a mineral oil and found deleterious effects under boundary lubrication conditions of a Hertzian contact.
Tao et al. [21] added 1% diamond nanoparticles with a 5-nm average diameter to paraffin oil. The results indicated that, under boundary-lubrication conditions, the nanofluid lubricant exhibited excellent load-carrying capacity, anti-wear, and friction-reduction properties. The authors postulated that the nanoparticles lodge into the rubbing surfaces due to their small size and in effect polish the surface. At the same time, some nanoparticles may form an inlaid pattern into the surface, which raises the hardness of the rubbing surface. Under extended rubbing time, the nanofluid forms a boundary lubricating film between the gliding surfaces. The NPs in the film not only bear the loading but also separate the solid parts and hence prevent their direct contact. They can also roll between the polished surfaces. This ball-bearing effect of nanoparticles also increases the load-bearing capacity while reducing wear and friction.
Subsequent studies [40][41][42][43] confirmed that adding nanoparticles to paraffin and water at low concentrations greatly improved the load-carrying capacity as well as the anti-wear and friction-reduction properties. A boundary protective film, consisting of nanoparticles, was found to form between the gliding surfaces [56,57]. Qiu et al. [58] suggested that the forming of a third body due to the deposition of nanoparticles on the contacting surfaces led to the observed lubrication behavior, which was supported in [44,45]. Some studies supported the mending effect, i.e., the healing of the surface during the tribological test, which may be due to exfoliation of NPs [46]. Nanoparticles fill the dimples and somewhat smooth the asperities, with a tribo-chemical reaction on the friction surfaces [47][48][49]59]. This is similar to the surface polishing effect proposed in [21]. Mosleh et al. [60] suggested two mechanisms of reducing asperity contact using nanoparticles: (1) filling the micro-valleys of contacting surfaces [51]; and (2) shearing NPs at the interface [61].
The mechanisms may be different, however, if nanoparticles with high aspect ratios or special material structures are used, e.g., inorganic nanotubes and fullerene-like nanoparticles [46,62]. In such cases, exfoliation of the platelet-shaped NPs may play a more significant role in friction and wear reduction [51,52,[63][64][65]. Figure 6 illustrates the mechanism of friction and wear reduction due to NP exfoliation. Specifically, NPs exfoliate under load and form an adherent film of thin flakes on the surfaces in contact. At high pressure, when greases or oils are squeezed out, this film can still act as a lubricant [66,67].
The lubrication effectiveness of the nanofluid matrix is strongly affected by a number of factors, including the physical and chemical properties of the nanoparticles, particle size, temperature in the contact region, particle concentration, particle aggregation, the presence of extreme pressure agents, type of base fluid, and the test conditions. The effects of some of these factors on lubrication properties are shown in Figures 7-9. Finer particles tend to provide better tribological properties, while an optimal particle concentration seems to exist for each NP-lubricant pair. At elevated temperatures, the viscosity of the base fluid is greatly reduced, leading to direct contact of the nanoparticles with the sliding surfaces, hence improving the lubrication effectiveness. Mechanically hard materials exhibit better load-carrying capacities and friction-reducing properties, due to their higher yield strength. However, they may lead to more surface wear. Nanofluid lubricants may be more effective in the boundary-lubrication regime, due to the high contact pressure that drains the fluids. In elasto-hydrodynamic lubrication and film-lubrication regimes, the liquid film plays a more important role [68][69][70]. surfaces [47][48][49]59]. This is similar to the surface polishing effect proposed in [21]. Mosleh et al. [60] suggested two mechanisms of reducing asperity contact using nanoparticles: (1) filling the microvalleys of contacting surfaces [51]; and (2) shearing NPs at the interface [61]. The mechanisms may be different, however, if nanoparticles with high aspect ratios or special material structures are used, e.g., inorganic nanotubes and fullerene-like nanoparticles [46,62]. In such cases, exfoliation of the platelet-shaped NPs may play a more significant role in friction and wear reduction [51,52,[63][64][65]. Figure 6 illustrates the mechanism of friction and wear reduction due to NP exfoliation. Specifically, NPs exfoliate under load and form an adherent film of thin flakes on the surfaces in contact. At high pressure, when greases or oils are squeezed out, this film can still act as a lubricant [66,67].
The lubrication effectiveness of the nanofluid matrix is strongly affected by a number of factors, including the physical and chemical properties of the nanoparticles, particle size, temperature in the contact region, particle concentration, particle aggregation, the presence of extreme pressure agents, type of base fluid, and the test conditions. The effects of some of these factors on lubrication properties are shown in Figures 7-9. Finer particles tend to provide better tribological properties, while an optimal particle concentration seems to exist for each NP-lubricant pair. At elevated temperatures, the viscosity of the base fluid is greatly reduced, leading to direct contact of the nanoparticles with the sliding surfaces, hence improving the lubrication effectiveness. Mechanically hard materials exhibit better load-carrying capacities and friction-reducing properties, due to their higher yield strength. However, they may lead to more surface wear. Nanofluid lubricants may be more effective in the boundary-lubrication regime, due to the high contact pressure that drains the fluids. In elastohydrodynamic lubrication and film-lubrication regimes, the liquid film plays a more important role [68][69][70]. The effective use of nanofluid lubricants have been also examined for other applications. For example, the lubrication performance of nanofluids for minimal quantity lubrication (MQL) in grinding processes was evaluated [71][72][73][74]. It was shown that nanoparticle-jet MQL achieved cooling and lubrication effects close to that of flood lubrication. Hybrid nanoparticles have better lubrication effects caused by the "physical synergistic effect," which integrates the lubrication advantages of  different NPs for better overall performance. Yu et al. [75,76] studied the role of the oxide-scale on the tribological behavior in the nanoparticle lubrication of hot rolling. Again, nanoparticles were found to be trapped in the cracks, which changed the wear rates. Zareh-Desari et al. [77] studied the effect of NP additives on lubrication performance in a deep drawing process and found that adding NPs to the base lubricant reduces the required forming load and improved the surface quality of the formed cups.  Stability of nanofluids is critical to their consistent performance in enhancing the tribological properties of the lubricants. A typical method to stabilize the nanofluid is to modify the NP surface tribological behavior in the nanoparticle lubrication of hot rolling. Again, nanoparticles were found to be trapped in the cracks, which changed the wear rates. Zareh-Desari et al. [77] studied the effect of NP additives on lubrication performance in a deep drawing process and found that adding NPs to the base lubricant reduces the required forming load and improved the surface quality of the formed cups.  Stability of nanofluids is critical to their consistent performance in enhancing the tribological properties of the lubricants. A typical method to stabilize the nanofluid is to modify the NP surface The effective use of nanofluid lubricants have been also examined for other applications. For example, the lubrication performance of nanofluids for minimal quantity lubrication (MQL) in grinding processes was evaluated [71][72][73][74]. It was shown that nanoparticle-jet MQL achieved cooling and lubrication effects close to that of flood lubrication. Hybrid nanoparticles have better lubrication effects caused by the "physical synergistic effect," which integrates the lubrication advantages of different NPs for better overall performance. Yu et al. [75,76] studied the role of the oxide-scale on the tribological behavior in the nanoparticle lubrication of hot rolling. Again, nanoparticles were found to be trapped in the cracks, which changed the wear rates. Zareh-Desari et al. [77] studied the effect of NP additives on lubrication performance in a deep drawing process and found that adding NPs to the base lubricant reduces the required forming load and improved the surface quality of the formed cups.
Stability of nanofluids is critical to their consistent performance in enhancing the tribological properties of the lubricants. A typical method to stabilize the nanofluid is to modify the NP surface [78][79][80][81]. To deal with the adverse effect of nanoparticle aggregation on the lubrication properties, Verma et al. [82] designed a low-cost approach to allow the intercalation of MoS 2 nanoparticles with organic agents, activating the inorganic nanoparticles. They proposed that the nanoparticles could bond to specific active agents and act as nano-reservoirs to store these additives. Once the NPs steer through the asperities of the contacting surfaces, they can deliver the organic agents and thus achieve specific functions. This idea is in the same vein as the methodology of targeted nanodrug delivery [5]. Recent experimental studies, focusing on reduced friction, anti-wear, and enhanced load-carrying capacity, have explored different NP materials and various surfactants to stabilize the mixture [52,70,[83][84][85][86][87][88][89][90].

Theoretical Background
The laboratory examples presented in Section 2.1 indicate that the use of nanoparticles in coolants and lubricants can be quite advantages, i.e., leading to improved efficiency and reliability, respectively. In this section, the underlying physics in terms of forces acting on nanoparticles as well as nanofluid properties are discussed in preparation for the subsequent sections on math modeling and computer simulations [91].

Forces Acting on Nanoparticles
In nanofluids, Brownian agitation of nanoparticles typically overcomes particle-sedimentation, resulting in mixture properties that are determined by both the nanoparticles and the base fluid. The very large surface-to-volume ratio and the very small inter-particle distance endow nanofluids with many unique properties. Fundamental are the forces experienced by each nanoparticle in the base fluid. Here, we give a brief discussion of the various NP-displacement velocities due to forces acting on the nanoparticles, explaining the microscopic mechanisms that give rise to the unique properties of nanofluids.
Random force. The random force ÝÑ F B ptq that derives from the collisions between particles and the randomly moving molecules of the base fluid gives rise to the Brownian motion of the nanoparticles. Based on the Langevin equation and the equi-partition theorem, the stochastic process causes an average instantaneous nanoparticle velocity of Ñ v p , which can be expressed as where κ B = 1.38ˆ10´2 3 m 2¨k g¨s´2¨K´1 is the Boltzmann constant, T is the local temperature, and Ñ v BM is the Brownian motion velocity of a typical nanoparticle with mass m p in a time scale of t ! τ, where τ is the particle momentum relaxation time. The collective effect of the random motion of nanoparticles is the dispersion of particles, which can be expressed by a dispersion coefficient, D 0 : where d p is the particle diameter. The effects of the Brownian motion can be incorporated by a random force, Ñ F Br , whose effect is to cause a dispersion equivalent to that predicted by Equation (3): where Ñ R is a random vector whose components are Gaussian random numbers with zero mean and unit variance, and ∆t is the integration time which is much higher than the characteristic time of molecular collisions.
Thermophoretic force. Particles can diffuse under the effect of a temperature gradient, where particles migrate from hotter to colder regions and cause a higher concentration of nanoparticles in the colder regions. This phenomenon is called thermophoresis and is equivalent to the Soret effect for gaseous or liquid mixtures. The thermophoresis effect can be expressed in terms of a thermophoretic velocity: or, equivalently, a thermophoretic force where the proportionality factor K tp is a function of the thermal conductivities of both the fluid and the solid particles. The expression of K tp for aerosols is given as where the parameters C s , C m , and C t can be determined empirically with numerical values of 1.17, 1.14, and 2.18, respectively, for Kn ă 0.1 [92]. Here, Kn " {d p is the Knudsen number, and is the mean free path of molecules. For thermophoresis of relatively large particles (1 µm) in liquids, the following expression can be used [93]: It should be noted that thermophoresis is a direct consequence of the Brownian motion of particles. It has the same cause as the Brownian diffusion of particles as discussed above.
Body force. Nanoparticles may experience different field forces depending on the types of applications. Typically, the gravitational force applies due to the density difference between the nanoparticles and the base fluids: which gives rise to a settling velocity, where µ b f is the base fluid viscosity, ∆ρ is the density difference between the nanoparticle and the base fluid, d p is the particle diameter, and g is the gravitational constant.
In an applied electrical field, particles will undergo electrophoretic motion with a velocity of v el , i.e., v el " where σ is the surface charge density of the particle, and E is the applied electric field. This force is due to an externally applied field. Forces due to the charges of NPs will be discussed later. Stokes force. When a nanoparticle moves relative to the base fluid, it experiences a drag force. At particle Reynolds numbers of Re ď 1, the drag force can be expressed according to Stokes' law as where Ñ v 1 p is the instantaneous particle velocity relative to the base-fluid motion. An important parameter in the analysis of particle motion in a fluid is the relaxation time defined as This characterizes the time required for a particle to adjust to or relax its velocity to a new condition of forces.
Inter-particle forces. When two particles approach each other due to Brownian motion, the attractive van der Waals forces and the repulsive electrostatic force become important.
The van der Waals forces cause attraction between particles. The attraction potential is proportional to the particle diameter, the material constant (the Hamaker constant A H ), and is inversely proportional to the distance of separation: where x " h{d p . For small particle separations (h ! d p ), Equation (13) is reduced to The average (minimum) separation distance between two spheres in a well-dispersed suspension can be obtained via where ϕ is the nanoparticle volume fraction. According to Equation (15), the average particle distance scales with the particle diameter and hence becomes very small for nanoparticles. When two particles approach each other, the ionic atmospheres overlap and create an ion concentration difference between the overlapped region and the bulk. This difference gives rise to an osmotic pressure acting to force the particles apart. The potential of this repulsive force can be obtained by integrating the force over distance [94]: where κ B is the Boltzmann constant, σ is the valence of the ions, κ is the Debye constant which is the inverse of the Debye length, ρ el is the concentration of electrons, e is the electronic charge, and φ 0 is the electric double layer (EDL) potential at the surface of the nanoparticles. The expression is valid when the value of the product, κd p , is greater than 20, and becomes invalid when the double-layer thickness is about the same as or larger than the particle radius. For the latter situation, the potential equations for the interaction between two double layers must be solved specifically for spherical particles. This is especially important for nanofluids because, for metal-oxide nanofluids and metal nanofluids, κd p is usually less than 4 [95].
Relative importance of these forces. An order-of-magnitude analysis of the Stokes force and the inter-particle forces in nanofluids has been provided in [96]. In dilute nanofluids, the surface forces, i.e., the van der Waals forces and the electrostatic force, are of the same order as the Stokes force if multi-particle interactions are considered. Meanwhile, Buongiorno [97] has shown that body forces are negligible compared to the random force that causes Brownian diffusion. However, the surface forces are only important when nanoparticles are well dispersed in the base fluid. When nanoparticle aggregation occurs, the surface forces become less important, while wake interactions due to Brownian motion dominate.

Nanoparticle Aggregation
Nanoparticle aggregation is one of the most important factors that may affect any consistent performance of the majority of nanofluids [98,99]. Due to the large surface area of NPs and the small average particle distance, NPs are prone to aggregation due to the inter-particle forces. Aside from potentially enhancing the effective thermal conductivity of a nanofluid, NP-aggregation may also lead to heterogeneous mixture properties and particle settling. However, it is very difficult to qualify, let alone quantify, the aggregation process due to the inevitable heterogeneous nature of nanofluids. It is a phenomenon caused by the inter-particle forces; however, it cannot be categorized as a nanofluid property. Statistical methods can be used only to a certain extent because of the high dependence on NP-material type and the NP-preparation method. One way to quantify NP-aggregation is to use the method developed by Nan et al. [100] for colloidal systems, as has been adopted in [14] among others. NP-aggregation affects nearly every property of the nanofluids by modifying the micro-structures of the mixture. It can be mediated by adding surfactants to the nanofluids or by altering the pH value of the mixture, both of which enhance the stability of nanoparticles in aqueous suspensions by introducing repulsive forces between NPs.

Thermo-Physical Properties of Nanofluids
Considering the forces experienced by NPs provides a microscopic view that serves as the basis for the discussion of macroscopic properties exhibited by nanofluids. In this section, we discuss the bulk properties of nanofluids and their mathematical description, which can then be used in the governing equations to predict nanofluid performances for various applications.
The basic particle-fluid mixture properties are greatly a function of the particle volume fraction. If particles are assumed to be spherical, mono-disperse, and forming a homogeneous dilute suspension, the nanofluid density and heat capacity in the case of low NP-volume fractions can be expressed as and where the subscripts bf, p indicate base fluid and particle, respectively. Thermal conductivity of nanofluids. Certain nanofluids have shown significantly enhanced thermal conductivities over their base fluids. It has been shown that the effective thermal conductivity of the mixture k n f increases with the nanoparticle volume fraction and temperature, and with the decrease of particle size (see Figure 5). In addition, particle shape, pH, and aggregation also have significant effects on k n f .
Assuming a monodisperse suspension spherical particles, Maxwell [101] derived the effective thermal conductivity of fluid-particle mixture as where k b f and k p are the thermal conductivities of the base fluid and the nanoparticle material, respectively. This model predicts the k e f f of mixtures with particle sizes greater than 1 µm. However, it significantly underestimates the effective thermal conductivity of nanofluids. Brownian motion and particle aggregation are arguably the most fundamental effects that characterize nanofluids. The Brownian motion of nanoparticles results in micro-scale mixing effects and causes an enhanced effective thermal conductivity, while locally small cluster-like structures due to aggregation can form thermal pathways and augment heat transfer, resulting in enhanced k n f -values.
For example, Kleinstreuer and Feng [96] proposed a model based on Brownian motion-induced micro-convection. Xu and Kleinstreuer [15] enhanced the F-K model to include the local clustering effect as well as the interfacial thermal resistance effect. The model divides the thermal conductivity k n f of the nanofluid into a static part and a micro-mixing part, i.e., k n f " k static`kmm . In this equation, k static is similar to Maxwell's model, except that the nanoparticle thermal conductivity k p in the original Maxwell model is replaced by the effective thermal conductivity of particle aggregate unit k a : while Here, The effective thermal conductivity of the nanoparticles including the thermal contact resistance, k pe f f , is given as In these equations, k de is the effective thermal conductivity of the mixture of base fluid and dead-end particles, ϕ de is the volume fraction of dead-end particles in an aggregate unit, ϕ b is the volume fraction of backbone particles in an aggregate unit, and L ii and β ii are parameters appearing in the derivation of the effective thermal conductivity of arbitrary isotropic particulate composites with interfacial thermal resistance. For a detailed explanation of these parameters, see Section 2.5 of Xu [91]. In addition, R K is the Kapitza resistance, ρ is the density, c p is the specific heat capacity, and T is the time-averaged temperature. The subscripts p and a denote nanoparticle and aggregate unit, respectively. C c («1 ) is a correction factor that depends on nanoparticle-based fluid parings. For example, C c = 1.1 for Al 2 O 3 -water nanofluids. The characteristic time interval τ a is expressed as where R g is the average radius of the aggregate units: Here, d p is the nanoparticle diameter, and T 0 = 273 K is the reference temperature. Comparison with measurements indicated that this model is applicable for several types of metal oxide nanoparticles in water for the range of volume fractions up to 10% and mixture temperatures below 350 K. A few different mechanisms have been proposed to explain the anomalous enhancement of the mixture thermal conductivity. For example, the base-fluid molecules form layered structures at the particle surface, acting as thermal bridge between the liquid and the solid particles, and thus enhance the k n f [102]. However, there are strong debates about whether they are really important to the enhanced thermal conductivity.
Viscosity of nanofluids. Measurements of the rheological properties of nanofluids indicate that most nanofluids behave as Newtonian-fluid mixtures at low particle volume concentrations [103,104]. Particle concentration, size [105], shape [106], temperature, shear rate [107], surfactants, and pH have direct impacts on the viscosity of nanofluids. For example, viscosity increases with higher particle volume fractions, while it decreases with elevated temperatures. However, due to the aggregation effect of nanoparticles, it is difficult to determine how the particle size affects the effective viscosity. Surfactants and pH affect µ n f mainly by enhancing the dispersion of nanoparticles.
Einstein [108] developed a formula for the effective viscosity of fluids containing spherical particles at low volume fractions (<2%). It indicates a linear increase in viscosity with particle volume fraction: where µ e f f and µ b f are the viscosities of the mixture and the base fluid, respectively. Brinkman [109] extended this model to consider moderate particle concentrations up to 4% : Comparison with experimentally measured data indicated that the Einstein model and the Brinkman model tend to underestimate the viscosity of nanofluids. In addition, these two models only account for the effect of particle loading on the mixture viscosity but neglect other important influencing factors such as temperature, particle size, and pH. In addition, newly proposed theoretical formulas are not able to accurately predict the effective viscosity for a wide range of nanofluids [104]. Since theoretical models for the viscosity of the mixture usually rely on the homogeneous-dispersion assumption, it is perhaps impossible to derive a general correlation that is able to accurately predict µ n f , given that the aggregation of nanoparticles almost always lead to a heterogeneous mixture. Still, due to the substantial difference among the available experimental data of µ n f , specific empirical correlations that are applicable to only one or a few types of nanofluids and within specific ranges of the pertinent parameters seem to be a more acceptable option at this moment. For example, Corcione [110] obtained an empirical correlation for µ n f with nanoparticles ranging in diameter from 25 nm to 200 nm, volume fractions ranging from 0.01% to 7.1%, and temperatures ranging from 293 to 333 K. The correlation reads where d bf is the equivalent diameter of a base fluid molecule, given by (29) in which M is the molecular weight of the base fluid, N is the Avogadro number, and ρ b f 0 is the mass density of the base fluid calculated at temperature T 0 = 293 K.
Convection heat transfer coefficient. Convection heat transfer performance is the most important measure of thermal nanofluid flow. Heat transfer performance for Fourier-type processes is characterized by the convective heat transfer coefficient or the Nusselt number, i.e., while where q w is the local wall heat flux, and Nu x is the Nusselt number at location x. The Nusselt number excludes the effect of increased pumping cost due to the dispersion of nanoparticles in the base fluid. Hence, the thermal performance factor may be used to evaluate the performance of thermal nanofluid flow: When corresponding nanofluid properties are used, the friction factor correlations developed for pure fluids may still be valid for nanofluids in both laminar and turbulent flow regimes. However, it is under debate whether the conventional correlations for the Nusselt number can remain accurate in predicting the nanofluid heat transfer characteristics [111,112]. While experimentally determined Nu of laminar nanofluid flow scatters, the Dittus-Boelter equation can still predict quite accurately the enhancement turbulent heat transfer coefficient of nanofluids [113], i.e., where n " 0.4 for heating (wall hotter than the bulk fluid), and n " 0.33 for cooling. Again, in Equation (33), nanofluid properties should be used in calculating Re and Pr.

Transport Equations and Modeling Approaches
The dominant forces that determine the nanofluid properties have been given in Section 2.1, while the bulk properties of nanofluids have been discussed in Section 2.2. In this section, the transport equations and modeling approaches for nanofluids are discussed. While there are a number of numerical techniques that can be used to simulate nanofluid flow, the underlying governing equations are all based on force balances and conservation laws. Focusing here on basic engineering applications, numerical techniques that could provide high-fidelity results and probe into the microscopic phenomena through extensive computations may be found elsewhere [114][115][116][117].

Single-Phase Approach
In the single-phase approach, the suspended nanoparticles are assumed to be in thermal equilibrium with the liquid phase, and the relative velocity between the two phases can be neglected. The reasoning is that nanoparticles are so small, say, d p < 100 nm, that they closely follow the fluid streamlines, causing the nanofluid to behave like a homogeneous mixture [103]. Thus, the single-phase Navier-Stokes equations can be used with the effective thermo-physical properties of nanofluids, as follows: Continuity equation: Momentum equation: Energy equation: (Laminar flow) where Φ is the energy dissipation term, Clearly, the single-phase approach has advantages in ease of use, computational cost, while achieving good comparison with selected experimental data sets. However, the results depend strongly on the selected thermo-physical property models, especially for the thermal conductivity and viscosity. As mentioned, particle migration can significantly affect heat transfer results by creating non-uniform concentration fields, especially in entrance regions and boundary layers, particularly at large Peclet numbers, Pe = Re¨Pr.
Incorporating the slip between the particles and the fluids can improve the single-phase model by adding a virtual term in the thermal conductivity expression [118]. Furthermore, solving the NP-mass transfer equation together with the momentum and energy equations can also increase the accuracy [119].

Two-Phase Approach
Even if the ultra-fine particles can be uniformly suspended in the base fluid, several factors such as Brownian motion and thermophoresis can significantly affect the relative particle motion, leading to particle slip. Hence, two-phase models, which solve one set of governing equations for each of the phases, may be more suitable. The most used two-phase model for nanofluid flow and heat transfer is the so-called Eulerian-Eulerian model, where both phases are considered as "fluids" with volume fractions summing up to unity [120]. Specifically, the governing equations for phase α are as follows: Continuity equation: B pϕ α ρ α q Bt`∇ pϕ α ρ α u α q " 0; Momentum equation: F α is the surface force due to the existence of another phase, and Q α is the interface heat transfer term. Aside from the inclusion of volume fraction ϕ, the difference in Equations (40)-(42) from Equations (34)-(39) is due to interfacial momentum and energy change. The interfacial forces that may apply to nanofluid flow are discussed in the above sections. The rate of interface heat transfer between the nanoparticles (α) and the fluid (β) is where A α,β is the interfacial area per unit volume. Some studies employed the Eulerian-Lagrangian framework, which tracks the particle trajectories in the fluid phase. Still, a mixture model that introduces a term due to the drift of nanoparticles into the momentum equation of the base fluid is feasible [121]. However, the comparisons between different two-phase models suggest no obvious winner in terms of prediction accuracy. Hence, the Eulerian two-phase model is favored due to ease of use and low computational cost. It should be noted that, although the two-phase models are capable of capturing the non-uniform concentration fields, the turbulence intensity of the fluid phase were often assumed unaffected by the presence of the nanoparticle phase.
In some cases, it is more desirable to use a convection-diffusion equation to represent the dispersion process of NPs in the liquid medium. For example, the nanomedicine distribution in the systemic system can be described by where c is the NP concentration, and D NP is the NP diffusion coefficient. This Euler-Euler method, which assumes small NP loading, one-way fluid-particle coupling, and negligible NP coagulation, is often preferred over the Euler-Lagrange approach due to its minimized computational expense.

Nanofluid Transport in Porous Media
Fluid flow in porous medium is governed by the Darcy-Forchheimer law [122,123]: The Darcy-Forchheimer law states that the momentum loss of flow through a porous medium is the result of two effects, i.e., the loss due to inertia and viscous effects. The first term on the right-hand side of the equation is the Darcy term due to the viscous effect, and the second term on the right-hand side is the Forchheimer term due to the inertia effect. Here, µ and ρ are the fluid viscosity and density, respectively, u j is the superficial velocity vector, K is the permeability, and C f is the inertia resistance factor. In the limit of negligible inertia effect (Re < 10, where Re " u ? K{ν), Equation (45) is reduced to Darcy's law, i.e., ∇p "´pµ{Kq Ñ v . The permeability α and the inertia coefficient C f are both functions of the porosity of the medium, which is defined as the ratio of the total voids volume and the total volume of the porous media. For example, Ergun [124] provided where is the size of the pore, and C = 150-180. The mathematical formulation of nanofluid flow in porous medium needs to incorporate the effect of the porosity and permeability. For example, assuming homogeneity and local thermal equilibrium in the porous medium, using Boussinesq approximation and Darcy's law, the governing equations for two-phase nanofluid natural convection can be written as [122,125] B Bx i pρ n f u i q " 0, 1 ε where the buoyancy term can be approximated by Here, the subscript m means the effective property of the porous medium. T is the nanofluid temperature, T 0 is a reference temperature, and β is the volume expansion coefficient. The superficial velocity u i is used in the equations.

Magnetic Nanofluids
Magnetic nanofluids behave like a paramagnetic material. In the absence of a magnetic field, the magnetic moments of particles are randomly orientated due to Brownian motion, and the fluid has no net magnetization. When an external magnetic field is applied, the magnetic moments align along the magnetic field direction, generating net magnetization to the fluid. In the presence of a magnetic field, the governing equations of the fluid flow and heat transfer need to be modified to include the Lorentz force and the corresponding energy. The continuity equation remains unchanged, while the momentum equation becomes where Ñ v is the velocity field, p is the pressure, D{Dt is the total derivative, The energy equation has the form where µ 0 " 4πˆ10´7 H¨m´1 is the magnetic constant, Ñ H is the magnetic field intensity, M is the magnetization of the fluid, and µ n f Φ is the viscous dissipation term. As such, the fluid flow and heat transfer of magnetic nanofluids can be manipulated by applying different magnetic fields.

Nanofluid Transport in Thin Films
The behavior of a nanoparticle embedded in a moving film can be modeled using the Reynolds lubrication equation, which governs the pressure distribution of thin viscous fluid films [126]: where p is the fluid-film pressure, x and y are the bearing width and length coordinates, and z is the vertical coordinate. In addition, h is the variable fluid-film thickness; u, v, and w are the film velocities in x-, y-, and z-direction, respectively; and subscripts a, b denote the top and bottom bounding wall-surfaces, respectively. This equation can be derived from the Navier-Stokes equation, which assumes laminar incompressible flow of Newtonian fluids. In addition, the inertial force is neglected, pressure is assumed to be constant across the film, and there is no velocity-slip at the boundaries. For a journal bearing with variable fluid-film thickness (or spacing h(θ, z)), the Reynolds equation can be simplified and normalized via [127] h " h{C, θ " x{R, µ " µ n f {µ b f , z " z{R, p " pC 2 {µ b f ωR 2 . (56) Then, the modified non-dimensional Reynolds equation can be expressed as where the normalized film thickness can be modified to incorporate the effect of adjustment on the convectional film thickness: h " 1`εcos`θ 1˘`R adj`Tadj .
Here, C is the radial clearance, R is the radius of the journal, L is the length of the bearing, R adj and T adj are the normalized radial and tilt adjustment, and ε " e{C is the normalized eccentricity.

Nanofluid Flow Applications
A major engineering goal in thermal nanofluid flow applications is to achieve high system performance, i.e., optimal efficiency via a reduction in entropy increase. The methodology is briefly discussed in Section 3.1 with an application to nanofluid flow in a microchannel. Nanoparticles often appear in nature or as drugs in liquid form. Thus, under thermal conditions, droplet-size changes due to evaporation or condensation must be accounted for. Section 3.2 provides a modern application to multi-component droplet-vapor interactions.
Finally, in Sections 3.3 and 3.4 microsystem cooling with nanofluids and improved lubrication with metallic nanoparticles are reviewed, respectively.

Entropy Generation
Most studies considering thermal nanofluid flow have focused on the enhanced heat transfer. However, to optimize the thermal systems for the efficient removal of high heat fluxes, entropy generation minimization must be considered along with the heat transfer maximization. Entropy generation is the measure of process irreversibility caused by both heat transfer and friction effects. It can be employed to assess the performance of thermal devices. The local entropy-generation rate (S gen in W/K) for convection heat transfer can be expressed as [103,128] S gen " S gen pthermalq`S gen pfrictionq.
Clearly, Equation (59) encapsulates the irreversibility due to heat transfer and frictional effects, where and The dissipation due to fluctuations in turbulent flow is roughly the same as that due to the mean flow [128]. In fact, the frictional entropy generation rate is only comparable to the thermal entropy generation rate when the flow Reynolds number is very large. It should be noted that, when using the above equations for nanofluids, the corresponding nanofluid properties should be used.
In the laminar regime of nanofluid flow, adding nanoparticles to the base fluid reduces entropy generation. An optimal NP concentration exists that corresponds to minimized entropy generation [15,129]. Moreover, the total entropy generation of the system can be reduced without significantly increasing the required pumping power [130]. Thus, adding nanoparticles to the base fluid can be very beneficial from a perspective of the second law of thermodynamics. In the turbulent flow regime, especially when the Reynolds number is very high, frictional entropy generation exceeds its thermal counterpart. In this case, adding NPs to the base fluid may push up the total entropy generation rate due to the augmented nanofluid viscosity [131]. Clearly, the thermo-physical properties of nanofluids play an important role in determining the entropy generation rate.

Nanodroplet-Vapor-Air Mixture Dynamics
A more specific set of governing equations are for the transport and phase change of inhaled nanodroplets in lung airways, which has rich roots in biomedical and engineering problems including administered drug-aerosols using inhalers, fume-or mist-nanodroplet inhalation from fuel handling, welding operation, spray processes, as well as smoke-aerosol inhalation from conventional cigarettes and electronic cigarettes. Their transport and conversion phenomena involve complex physical and thermodynamic interactions between air, vapors, and nanodroplets with their therapeutic (or toxic) compounds and the respiratory tract walls. This is an area of interdisciplinary applications of modern thermal internal nanofluid-flow engineering. A basic mathematical model can be outlined as follows, which is based on the work of Feng et al. [132].
Any given air-vapor mixture is described as a single continuous phase, i.e., the conservation laws for mass, momentum, and energy describing the air-vapor mixture transport and the advection-diffusion equation describing vapor-species transport. Droplet-vapor interaction, i.e., liquid-vapor mass change due to evaporation or condensation, is implemented by: (i) incorporating the energy generation/destruction due to phase change via source terms in the energy equation and vapor-species transport equations (see Equations (62)-(65)); and (ii) calculating the local vapor-mass fraction in the droplet mass conservation equation (see Equation (66)). Specifically, the energy equation of the air-vapor mixture reads where the energy source term S pEq v´d is due to the latent heat of evaporation or condensation which is released or absorbed by the droplets per local mesh cell: , . - where N d,cell is the total droplet number in a specified mesh cell. The governing equation for advection and diffusion of the sth vapor species is expressed as where Sc T " 0.9 is the turbulence Schmidt number for Y v,s , v T is the turbulent viscosity, and r D a´v,s is the molecular diffusivity of the sth vapor species in the air-vapor mixture. The local vaporized/condensed vapor-mass flow rate of the aerosol components are added to its advection-diffusion equation as a source term S pYq v´d,s (kg¨m´3¨s´1), i.e., S pYq v´d,s " where j s is the average evaporation/condensation mass flux normal to the droplet surface of the sth component (i.e., j s ą 0 for evaporation and j s ă 0 for condensation); and A d is the droplet surface area. dt d represents the droplet phase time differential, and ∆t f is the flow time-step. For the multicomponent droplet size-change and transport, the governing equations for discrete droplets are the translational equation as well as the mass and energy conservation laws.
Newton's second law of motion states where F D i , F L i , F BM i , and F G i are the drag force, lift force, Brownian motion-induced force, and gravity, respectively. Droplet mass changes due to condensation/evaporation can be described as where the average mass flux j s is given by [133]: Here, Sh is the Sherwood number (see the Nomenclature for definition), while Y s,vsur f and Y s,v,cell are the mass fractions of the sth vapor phase at the droplet surface and at the center of the cell, where the droplet currently resides. The correction factor C m for submicron droplets, considering non-continuum effects, can be expressed as where Kn is the Knudsen number, Kn " 2λ a´v {d d , in which λ a´v is the mean free path of the air-vapor mixture surrounding the droplet; and α m is the mass accommodation coefficient, where α m " 1 was used in the present study [134]. The derivation of Equation (68) is based on the assumption that the distance between the droplet mass center and the mesh cell center is much larger than the droplet radius. Specifically, Y s,v,cell is determined by the advection-diffusion equation (see Equation (64)), while Y s,vsur f can be calculated by Here, ρ s,vsur f is the equivalent density of vapor species s in the air-vapor mixture at the droplet surface; γ s is the activity coefficient of species s, which is a correction of the evaporation/condensation characteristics of a certain liquid component due to the molecular bonding changes in the pure droplets and in the multicomponent droplets [135]; X s,d is the mole fraction of the sth component in the droplet; s is the species gas constant; T d is the droplet temperature; and P s,vsat is the temperature-dependent saturation pressure of the pure sth species. Clearly, parameters γ s and X s,d serve as corrections to the vapor pressure at the surface of multicomponent droplets (Raoult's law), while K s is the correction factor for the Kelvin effect [134]: where σ s is the surface tension of component s at the droplet surface. The droplet heat transfer equation takes on the form where L s is the latent heat of liquid-vapor phase transition of the sth species (L s is always larger than zero). Moreover, C T is the correction factor for submicron droplets: where α T is the thermal accommodation coefficient with the value of α T " 1 [134]. The governing equations of the transport and phase change of the multicomponent droplet-vapor mixtures can be solved using computational fluid-particle dynamics (CF-PD) methods [103].

Nanofluids for Microsystem Cooling
Experimental studies have shown somewhat controversial results regarding the heat transfer performance of nanofluids in microsystems. Whether or not nanofluids outperform their base fluids is unclear. On the other hand, numerical studies tend to show that nanofluids have better heat transfer performance in microsystems as compared to their base fluids. This part reviews the numerical studies on microsystem cooling using nanofluids. As the numerical study of nanofluid flow and heat transfer in microsystems are abundant, only selective studies are reviewed here, which provide a good cover of the history and the state of the art.
By assuming that nanofluid has higher thermal conductivity and similar viscosity as base fluid, Lee and Choi [136] demonstrated that nanofluids offer a dramatic enhancement of cooling rates by reducing the thermal resistance of the microchannel heat sink significantly. Using a similar approach with empirical correlations of nanofluid properties, Chein and Huang [137] showed that the performance of a MCHS is greatly improved by using nanofluids in the laminar flow regime, without causing extra pressure loss. Koo and Kleinstreuer [138] conducted a numerical study of nanofluid flow in MCHS's using theoretical models of effective thermal conductivity and dynamic viscosity of nanofluids. They demonstrated that, at low volume fractions (<4%), nanofluid significantly increases the heat transfer performance of microheat-sinks, given that particle accumulation is avoided. High aspect-ration microchannel is recommended for better performance. Jang and Choi [139] numerically studied the effectiveness of using water-diamond nanofluid in a MCHS and found that, at a fixed pumping power, nanofluid with 1 vol % enhanced the cooling performance by 10% as compared to water. A theoretical model for nanofluid thermal conductivity was developed and used in this study, while the Einstein equation was adopted for nanofluid viscosity. Li and Kleinstreuer [140] studied heat transfer in a trapezoidal microchannel using a theoretical k nf model, whose predictions agree well with available measured data in the literature. They suggested that nanofluids measurably enhance the thermal performance of microchannel mixture flow with a small penalty in pumping power. Lelea [141] and Lelea and Nisulescu [142] numerically investigated nanofluid flow in microchannels. They found that the enhancement of heat transfer due to the presence of NPs rises along the channels. The heat transfer enhancement for cooling and heating are different for the same NP concentration. Additionally, neglecting viscous heating effects would give incorrect predictions. Ting et al. [143] investigated the effects of streamwise conduction on the thermal performance of nanofluid flow in MCHS's under exponentially decaying wall heat flux. The Peclet number was found to strongly affect the temperature distribution when the streamwise conduction is incorporated. The effect of streamwise conduction became important when Pe < 10.
MCHS's can be treated as porous media so that, using an averaging approach, the heat transfer and friction characteristics can be obtained without considering each individual channels. Tasi and Chein [144] employed the two-equation averaging approach to study nanofluid cooled MCHS's, where the MCHS's were treated as a porous medium. They concluded that adding nanoparticles to the base fluid reduced the temperature difference between the heated wall and the fluid due to a reduction in the conductive thermal resistance. However, adding NPs increases the convective thermal resistance of the MCHS due to increased viscosity and decreased thermal capacity. Optimum values of channel aspect ratio and porosity exist, above which the thermal resistance of the MCHS is increased by adding NPs to the base fluid. Chen and Ding [145] employed the two-equation porous media approach to evaluate the thermal performance of a nanofluid-cooled MCHS. The model was tuned based on experimental measurements. Specifically, the effect of the inertia force term of the Forchheimer-Brinkman-extended Darcy equation on the thermal performance was assessed. Hatami and Ganji [146] treated MCHS's cooled by nanofluids as a porous medium and suggested high-channel aspect ratios, large NP size, and wide fins for better thermal performance.
Many studies have used the two-phase modeling approach to study nanofluid flow in microsystems. Kalteh et al. [147] used a Eulerian two-phase model to numerically study nanofluid laminar forced convection in a microchannel. The heat transfer enhancement was found to increase with the Re number and NP concentration, and with the decrease of the NP diameter. The model also predicted that the velocity and temperature difference between the fluid and the NP phase was negligible and that the NP concentration was uniform. Kalteh et al. [148] used experimental measurements to validate single-phase and two-phase nanofluid flow in MCHS's. The two-phase model showed better agreement with the measurements, though the velocity and temperature difference between the two phases was shown to be negligible. Fani et al. [149] used the two-phase model to study the size effect of the nanoparticle on the thermal performance of nanofluids in MCHS's. The increment of the NP size was shown to cause an increase in pressure drop and a decrease in heat transfer. Ghale et al. [150] evaluated single-phase and two-phase mixture models for the heat transfer of nanofluids in microchannels and concluded that the two-phase model achieved better agreement with experimental measurements. They also demonstrated that ribs at microchannels increased both the Nusselt number and the friction coefficient of nanofluids.
Alternative modeling methods have also been used to numerically study microsystem cooling using nanofluids. Sun et al. [115] reported an equilibrium molecular dynamics simulation of nanofluid thermal conductivity in a nanochannel. The thermal conductivity of the mixture was found to be anisotropic due to the confinement of the wall. Its value decreased with the increase of the channel width owing to the wakened NP-wall interaction. Ahmed and Eslamian [151] used the multi-phase lattice Boltzmann method to evaluate the importance of particle-fluid slip on the laminar forced convection heat transfer in a microchannel. The results revealed the significance of Brownian force at Re < 10, where the nanofluid flow was heterogeneous and provided enhanced heat transfer as compared to the base fluid. At a higher Re number, the nanofluid flow was homogeneous, and no heat transfer augmentation was achieved. Thermophoresis effect was found to be insignificant at all Re number ranges. Karimipoura et al. [114] numerically studied nanofluid flow in a microchannel in slip flow regime using the lattice Boltzmann method. The results showed that decreasing the values of the slip coefficient enhances the convective heat transfer coefficient.
Entropy generation is a measure of the thermal performance of a system. It is an important factor to consider especially when designing a microsystem for enhanced heat transfer, where system efficiency is paramount [152]. Li and Kleinstreuer [129] conducted entropy generation analysis for nanofluid flow in trapezoidal microchannels with different aspect ratios. Optimal Reynolds number ranges were found to exist to operate the system based on minimization of total entropy generation. A suitable operational Reynolds-number range was found to be lower for microchannels with high aspect ratios. Ting et al. [153] analytically evaluated the effect of stream-wise conduction on the entropy generation of nanofluid flow in circular MCHS's under an exponentially decaying wall heat flux and constant pumping power conditions. The effect of stream-wise conduction was found to be important when the Peclet number was small, i.e., irreversibility due to heat transfer dominants. Ting et al. [154] evaluated the entropy generation of viscous dissipative nanofluid flow in thermal non-equilibrium porous media embedded in microchannels. Results showed that, in a low-aspect-ratio microchannel, adding NPs to the base fluid decreased the thermodynamic efficiency, while using nanofluids in a high-aspect-ratio microchannel enhanced exergetic effectiveness in a low-Re flow regime. Cruz-Duarte et al. [155] demonstrated the design of a MCHS with nanofluids through the entropy generation minimization criterion.
In addition to microchannels and micro-heat sinks, alternative microsystems were considered for enhanced heat transfer using nanofluids [156]. For example, Mohammed et al. [157,158] numerically studied the thermal performance of microchannel heat exchangers using nanofluids for both parallel-flow and counter-flow scenarios. The results indicate that increasing NP loading and Reynolds number cause an increase in heat transfer performance of counter-flow heat exchangers, but deteriorations for parallel-flow scenarios. A slight increase in pressure drop was found due to the presence of NPs. An optimization study of microchannel heat exchangers was performed in [159]. Islami et al. [160,161] studied heat transfer and fluid flow in microchannels containing micromixers. Newtonian and non-Newtonian base fluids were considered. Heat transfer enhancement was shown to be more significant with Newtonian base fluids. Lelea and Laza [162] considered a micro-heat sink with straight microtubes and multiple tangential inlet jets. Employing a nanofluid showed little advantage as compared to the base fluid when a fixed pumping power was used.
Geometric optimization of nanofluid-cooled MCHS's has been conducted in a number of studies using different techniques [163][164][165][166]. For example, Mital [167] analyzed laminar, developing nanofluid flow in microchannels using empirical correlations. The results indicated that for an optimal volume fraction, heat transfer can be maximized at a given Re number.
Nanoparticle distribution in the flow field, especially near the wall, may have a significant effect on the fluid flow and heat transfer characteristics of nanofluids in microchannels. Alvariño et al. [168] assessed the self-diffusion of fully developed nanofluid flow in microchannels using a new CFD solver. They showed that, while thermophoresis may be negligible, self-diffusion can be significant. Furthermore, contrary to what has usually been assumed, self-diffusion, Brownian diffusion, and thermophoresis cannot be superimposed when determining a NP-concentration field. Hedayati et al. [169] investigated the effect of NP migration and asymmetric heating on the convection heat transfer of nanofluid flow in microchannels considering Brownian motion and thermophoretic diffusion. A particle depletion layer was detected near the wall, which contributed to an enhanced heat transfer rate. This augmentation was found to be enhanced by a slip-boundary due to microscopic roughness. Nikkhaha et al. [170] numerically studied nanofluid flow in a microchannel with oscillating heat flux and slip boundary condition. Higher slip-coefficients and NP-volume fractions were found to increase the Nusselt number.
Convection heat transfer characteristics of nanofluids in microsystems in the presence of a magnetic field have also been the subject of a number of studies [171,172]. Malvandi and Ganji [173] and Malvandi et al. [174] found that a strong magnetic field had a negative effect on nanofluid performance. Hajialigol et al. [175] investigated nanofluid magneto-hydrodynamic mixed convection and entropy generation in a microchannel and revealed that the heat transfer increases with the Hartmann number. The total entropy generation was shown to decrease with increasing magnetic strength and volume fraction.
In conclusion, numerical studies indicate a better heat transfer performance of nanofluids in microsystems as compared to that of the base fluids; clearly, it is largely due to the models used for the effective thermal conductivities. Meanwhile, significant penalty in pressure drop were recorded by adding NPs to the base fluids.

Nanoparticles for Enhanced Lubrication
There are a very few numerical studies of enhanced lubrication using nanoparticles. This can be attributed to the fact that the lubrication process often involves the deformation of NPs, which is difficult to model. Nevertheless, some studies have revealed lubrication characteristics using numerical methods. This section reviews these studies on enhanced lubrication using nanoparticles.
The friction-related mechanisms of lubricants with nanoparticles can be deduced by analyzing scanning electron microscopy images and energy-dispersive spectrometry patterns of the wear surface. In addition, it can be probed by molecular level simulations. For example, Lv et al. [176] first used molecular dynamics (MD) simulation to study the friction mechanisms using nanofluid lubricants. They showed the rotations of NPs under low pressure as well as trapping and aggregations of NPs under higher pressures. These effects separated the friction surfaces and led to enhanced lubrication properties. The molecular dynamic simulation of Hu et al. [22] demonstrated that phase transition took place under high load, while the transition pressure of an NP-added lubricant was higher than that of the base lubricant. Nanoparticles could effectively increase the load-carrying capacity, especially for NPs with smaller sizes, due to stronger micro-motions and volume effects. Hu et al. [39] compared the differences in the mechanical properties of contact layers between friction surfaces with and without Cu-nanoparticles. The MD simulation results demonstrated better improvement in friction properties with Cu-nanoparticles at a low velocity compared with that at a high velocity. At low velocities, a Cu-nanofilm that was formed on the friction surface beard the velocity gradient and plastic deformation. Under high velocity conditions, a transfer layer that appeared adjacent to the interface reduced the friction forces. Hu et al. [177] numerically studied the tribological properties of two hard NPs confined by two iron blocks using an MD simulation. They showed that, at a low velocity and low load, the NPs separated the two blocks from each other and acted as ball bearings. At high velocity and high load, the support effect of these nanoparticles was lost in a short sliding time due to the crushing of the NPs. In another study using an MD simulation, Hu et al. [178] demonstrated that the improved load-carrying capacity is due to the organized molecular structures of the lubricant near the NPs, and to the deformation of nanoparticles.
Shenoy et al. [127] used a modified Reynolds equation to numerically investigate the effect of a nanoparticle additive on the fluid film of a journal bearing. The results indicated that dispersing TiO 2 nanoparticles in engine oil can reduce end-leakage as well as friction, thus improving the load capacity of the bearing featuring negative radial and tilt adjustments. Binu et al. [179] used the modified Krieger-Dougherty viscosity model to incorporate the nanoparticle effect on the fluid rheology. They introduced coupled stresses to the Reynolds equation to account for the impact of nanoparticles. The model predicted a significant increase in load-carrying capacity of the journal bearing by using a TiO 2 -based lubricant. Nicoletti [180] found that adding nanoparticles to the base oil modifies the volumetric heat capacity of the lubricant, resulting in lower temperature and thus higher viscosity for an improved load-carrying capacity.
Österle et al. [181] investigated the effect of the presence of soft and hard nanoparticles on the friction performance during dry sliding, employing the method of movable cellular automata. The results revealed that agglomerates of soft NPs decomposed and mixed with the oxide layer and reduced the friction coefficient, while that of hard NPs mixed with the oxide in the third body layer without decomposing. Increasing the amount of hard NPs would increase the friction coefficient. In another study, Österle et al. [182] numerically investigated the effect of the presence of SiC-NPs on the lubrication between tribofilms. They found that, though increasing the thickness of the tribofilm significantly reduced friction, the volume fraction of SiC-NPs only marginally affected the friction. Dmitriev et al. [183] showed that embedding silica-NPs in a polymer matrix tribofilm significantly affected the interface structure and smoothness of the sliding mechanism. Österle et al. [126] demonstrated that reduced friction and wear can be expected if the silica tribofilms are mixed with at least 10 vol % graphite.

Conclusions and Future Work
Adding nanoparticles to base liquids at low concentrations creates a complex mixture, where the small size of the particles and their large surface-to-volume ratio cause new physical phenomena that have not been observed in other suspensions. Hence, as part of the rich applications of nanofluids, the underlying physics of the nanoparticle-fluid suspension and their mathematical description need to be properly addressed. Nanoparticle aggregation poses a great challenge, due to the difficulties in characterizing the process and the potential deterioration of nanofluid performance. Clearly, more refined mathematical models are needed to accurately account for the aggregation process.
Nanofluids can provide superior cooling performance in microsystems for laminar intermediate-to-high Reynolds number flow conditions, without causing an excessive increase in pressure drop. However, possible nanoparticle deposition, i.e., microchannel fouling, need to be properly addressed to ensure good performance of nanofluid flow and heat transfer. Numerical studies of nanofluid flow and heat transfer in microsystems are providing insight to improved system performance.
Using nanoparticles as additives to lubricants can greatly enhance their tribological properties by reducing friction and wear, as well as increase the load-carrying capacity. Low weight percentages of nanoparticles seem to provide the best performance. The mechanisms for enhanced lubrication performance are manifold and depend on the material structure of the nanoparticles. Under film-lubrication conditions, where the load is supported by the lubricant and any solid-solid contact is avoided, the potential increase of the lubricant's viscosity due to nanoparticle addition needs to be considered. Again, a homogeneous suspension of nanoparticles in lubricants is vital for proper applications, as large aggregates of nanoparticles may be detrimental to the lubrication process.
Acknowledgments: Zelin Xu acknowledges partial financial support via a CSC (Chinese Scholarship Council) award.
Author Contributions: Zelin Xu conducted the literature review and drafted and formatted the manuscript; Clement Kleinstreuer designed and revised the manuscript and finalized the submitted version.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: MCHS microchannel heat sink MD molecular dynamics NP nanoparticle