Numerical investigation on the heat transfer enhancement using micro/nano phase-change particulate flow

The introduction of phase change material fluid and nanofluid in micro-channel heat sink design can significantly increase the cooling capacity of the heat sink because of the unique features of these two kinds of fluids. To better assist the design of a high performance micro-channel heat sink using phase change fluid and nanofluid, the heat transfer enhancement mechanism behind the flow with such fluids must be completely understood. ^ A detailed parametric study is conducted to further investigate the heat transfer enhancement of the phase change material particle suspension flow, by using the two-phase non-thermal-equilibrium model developed by Hao and Tao (2004). The parametric study is conducted under normal conditions with Reynolds numbers of Re = 90–600 and phase change material particle concentrations of ϵp ≤ 0.25, as well as extreme conditions of very low Reynolds numbers (Re < 50) and high phase change material particle concentration (ϵp = 50%–70%) slurry flow. By using the two newly-defined parameters, named effectiveness factor ϵeff and performance index PI, respectively, it is found that there exists an optimal relation between the channel design parameters L and D, particle volume fraction ϵp, Reynolds number Re, and the wall heat flux qw. The influence of the particle volume fraction ϵp, particle size dp, and the particle viscosity μ p, to the phase change material suspension flow, are investigated and discussed. The model was validated by available experimental data. The conclusions will assist designers in making their decisions that relate to the design or selection of a micro-pump suitable for micro or mini scale heat transfer devices. ^ To understand the heat transfer enhancement mechanism of the nanofluid flow from the particle level, the lattice Boltzmann method is used because of its mesoscopic feature and its many numerical advantages. By using a two-component lattice Boltzmann model, the heat transfer enhancement of the nanofluid is analyzed, through incorporating the different forces acting on the nanoparticles to the two-component lattice Boltzmann model. It is found that the nanofluid has better heat transfer enhancement at low Reynolds numbers, and the Brownian motion effect of the nanoparticles will be weakened by the increase of flow speed. ^

Lattice Boltzmann method for flow and heat transfer in microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

A study of branching fluid networks for enhancing the performance of thermal -fluid devices

Compact thermal-fluid systems are found in many industries from aerospace to microelectronics where a combination of small size, light weight, and high surface area to volume ratio fluid networks are necessary. These devices are typically designed with fluid networks consisting of many small parallel channels that effectively pack a large amount of heat transfer surface area in a very small volume but do so at the cost of increased pumping power requirements. ^ To offset this cost the use of a branching fluid network for the distribution of coolant within a heat sink is investigated. The goal of the branch design technique is to minimize the entropy generation associated with the combination of viscous dissipation and convection heat transfer experienced by the coolant in the heat sink while maintaining compact high heat transfer surface area to volume ratios. ^ The derivation of Murray's Law, originally developed to predict the geometry of physiological transport systems, is extended to heat sink designs which minimze entropy generation. Two heat sink designs at different scales are built, and tested experimentally and analytically. The first uses this new derivation of Murray's Law. The second uses a combination of Murray's Law and Constructal Theory. The results of the experiments were used to verify the analytical and numerical models. These models were then used to compare the performance of the heat sink with other compact high performance heat sink designs. The results showed that the techniques used to design branching fluid networks significantly improves the performance of active heat sinks. The design experience gained was then used to develop a set of geometric relations which optimize the heat transfer to pumping power ratio of a single cooling channel element. Each element can be connected together using a set of derived geometric guidelines which govern branch diameters and angles. The methodology can be used to design branching fluid networks which can fit any geometry. ^

Lattice Boltzmann Method for Flow and Heat Transfer in Microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

Design construction and testing of a heat loop for study of fouling and its interrelationship to corrosion in heat exchanger tubes

A heat loop suitable for the study of thermal fouling and its relationship to corrosion processes was designed, constructed and tested. The design adopted was an improvement over those used by such investigators as Hopkins and the Heat Transfer Research Institute in that very low levels of fouling could be detected accurately, the heat transfer surface could be readily removed for examination and the chemistry of the environment could be carefully monitored and controlled. In addition, an indirect method of electrical heating of the heat transfer surface was employed to eliminate magnetic and electric effects which result when direct resistance heating is employed to a test section. The testing of the loop was done using a 316 stainless steel test section and a suspension of ferric oxide and water in an attempt to duplicate the results obtained by Hopkins. Two types of thermal ·fouling resistance versus time curves were obtained . (i) Asymptotic type fouling curve, similar to the fouling behaviour described by Kern and Seaton and other investigators, was the most frequent type of fouling curve obtained. Thermal fouling occurred at a steadily decreasing rate before reaching a final asymptotic value. (ii) If an asymptotically fouled tube was cooled with rapid cir- ·culation for periods up to eight hours at zero heat flux, and heating restarted, fouling recommenced at a high linear rate. The fouling results obtained were observed to be similar and 1n agreement with the fouling behaviour reported previously by Hopkins and it was possible to duplicate quite closely the previous results . This supports the contention of Hopkins that the fouling results obtained were due to a crevice corrosion process and not an artifact of that heat loop which might have caused electrical and magnetic effects influencing the fouling. The effects of Reynolds number and heat flux on the asymptotic fouling resistance have been determined. A single experiment to study the effect of oxygen concentration has been carried out. The ferric oxide concentration for most of the fouling trials was standardized at 2400 ppM and the range of Reynolds number and heat flux for the study was 11000-29500 and 89-121 KW/M², respectively.

Diatom abundance in surface sediments of the Southern Ocean

Based on the quantitative analysis of diatom assemblages preserved in 274 surface sediment samples recovered in the Pacific, Atlantic and western Indian sectors of the Southern Ocean we have defined a new reference database for quantitative estimation of late-middle Pleistocene Antarctic sea ice fields using the transfer function technique. The Detrended Canonical Analysis (DCA) of the diatom data set points to a unimodal distribution of the diatom assemblages. Canonical Correspondence Analysis (CCA) indicates that winter sea ice (WSI) but also summer sea surface temperature (SSST) represent the most prominent environmental variables that control the spatial species distribution. To test the applicability of transfer functions for sea ice reconstruction in terms of concentration and occurrence probability we applied four different methods, the Imbrie and Kipp Method (IKM), the Modern Analog Technique (MAT), Weighted Averaging (WA), and Weighted Averaging Partial Least Squares (WAPLS), using logarithm-transformed diatom data and satellite-derived (1981-2010) sea ice data as a reference. The best performance for IKM results was obtained using a subset of 172 samples with 28 diatom taxa/taxa groups, quadratic regression and a three-factor model (IKM-D172/28/3q) resulting in root mean square errors of prediction (RMSEP) of 7.27% and 11.4% for WSI and summer sea ice (SSI) concentration, respectively. MAT estimates were calculated with different numbers of analogs (4, 6) using a 274-sample/28-taxa reference data set (MAT-D274/28/4an, -6an) resulting in RMSEP's ranging from 5.52% (4an) to 5.91% (6an) for WSI as well as 8.93% (4an) to 9.05% (6an) for SSI. WA and WAPLS performed less well with the D274 data set, compared to MAT, achieving WSI concentration RMSEP's of 9.91% with WA and 11.29% with WAPLS, recommending the use of IKM and MAT. The application of IKM and MAT to surface sediment data revealed strong relations to the satellite-derived winter and summer sea ice field. Sea ice reconstructions performed on an Atlantic- and a Pacific Southern Ocean sediment core, both documenting sea ice variability over the past 150,000 years (MIS 1 - MIS 6), resulted in similar glacial/interglacial trends of IKM and MAT-based sea-ice estimates. On the average, however, IKM estimates display smaller WSI and slightly higher SSI concentration and probability at lower variability in comparison with MAT. This pattern is a result of different estimation techniques with integration of WSI and SSI signals in one single factor assemblage by applying IKM and selecting specific single samples, thus keeping close to the original diatom database and included variability, by MAT. In contrast to the estimation of WSI, reconstructions of past SSI variability remains weaker. Combined with diatom-based estimates, the abundance and flux pattern of biogenic opal represents an additional indication for the WSI and SSI extent.

Heat transport in low dimensions : introduction and phenomenology

Review paper, to appear in the Springer Lecture Notes in Physics volume "Thermal transport in low dimensions: from statistical physics to nanoscale heat transfer" (S. Lepri ed.)

Investigação experimental das propriedades termofísicas e da convecção forçada de nanofluido de grafeno

This work presents an experimental investigation of thermal hydraulic performance of the nanofluid composed by graphene nanoparticles dispersed in a mixture of water and ethylene glycol at a ratio of 70:30% by volume. The tests were carried out under forced convection inside a circular tube with uniform heat flux on the wall for the laminar-turbulent transition regime. The mass flow rate ranged from 40 to 70 g/s corresponding to Reynolds numbers between 3000 and 7500. The heat flux was maintained constant at values of 11, 16 and 21 kW/m², as well as the inlet temperature of 15, 20 and 25°C. Three samples were produced with the nanofluid volumetric concentration of 0.05%, 0.10% and 0.15%. Thermophysical properties were experimentaly measured for all samples that were critically compared and discussed with theoretical models most commonly used in the literature. Initially, experiments with distilled water confirmed the validity of the experimental equipment for the thermo-hydraulic tests. Therefore, nanofluid samples that showed the highest thermal conductivity, corresponding to the volumetric concentrations of 0.15% and 0.10%, were subjected to the tests. The thermal-hydraulic performance for both samples was unsatisfactory. The heat transfer coefficients for convection of nanofluids reduced 21% in average, for the sample with = 0.15% and 26% and for =0.10%. The pressure drop of the samples was higher than the base fluid. Finally, the pressure drop and heat transfer coefficient by convection of both samples were also compared to theoretical models. The models used for pressure drop showed an excellent agreement with experimental results, which is remarkable considering the transitional flow.

Constructal Design of Energy Systems

This dissertation shows the use of Constructal law to find the relation between the morphing of the system configuration and the improvements in the global performance of the complex flow system. It shows that the better features of both flow and heat transfer architecture can be found and predicted by using the constructal law in energy systems. Chapter 2 shows the effect of flow configuration on the heat transfer performance of a spiral shaped pipe embedded in a cylindrical conducting volume. Several configurations were considered. The optimal spacings between the spiral turns and spire planes exist, such that the volumetric heat transfer rate is maximal. The optimized features of the heat transfer architecture are robust. Chapter 3 shows the heat transfer performance of a helically shaped pipe embedded in a cylindrical conducting volume. It shows that the optimized features of the heat transfer architecture are robust with respect to changes in several physical parameters. Chapter 4 reports analytically the formulas for effective permeability in several configurations of fissured systems, using the closed-form description of tree networks designed to provide flow access. The permeability formulas do not vary much from one tree design to the next, suggesting that similar formulas may apply to naturally fissured porous media with unknown precise details, which occur in natural reservoirs. Chapter 5 illustrates a counterflow heat exchanger consists of two plenums with a core. The results show that the overall flow and thermal resistance are lowest when the core is absent. Overall, the constructal design governs the evolution of flow configuration in nature and energy systems.

Design of Latent Thermal Energy Storage Systems

This dissertation documents the results of a theoretical and numerical study of time dependent storage of energy by melting a phase change material. The heating is provided along invading lines, which change from single-line invasion to tree-shaped invasion. Chapter 2 identifies the special design feature of distributing energy storage in time-dependent fashion on a territory, when the energy flows by fluid flow from a concentrated source to points (users) distributed equidistantly on the area. The challenge in this chapter is to determine the architecture of distributed energy storage. The chief conclusion is that the finite amount of storage material should be distributed proportionally with the distribution of the flow rate of heating agent arriving on the area. The total time needed by the source stream to ‘invade’ the area is cumulative (the sum of the storage times required at each storage site), and depends on the energy distribution paths and the sequence in which the users are served by the source stream. Chapter 3 shows theoretically that the melting process consists of two phases: “invasion” thermal diffusion along the invading line, which is followed by “consolidation” as heat diffuses perpendicularly to the invading line. This chapter also reports the duration of both phases and the evolution of the melt layer around the invading line during the two-dimensional and three-dimensional invasion. It also shows that the amount of melted material increases in time according to a curve shaped as an S. These theoretical predictions are validated by means of numerical simulations in chapter 4. This chapter also shows that the heat transfer rate density increases (i.e., the S curve becomes steeper) as the complexity and number of degrees of freedom of the structure are increased, in accord with the constructal law. The optimal geometric features of the tree structure are detailed in this chapter. Chapter 5 documents a numerical study of time-dependent melting where the heat transfer is convection dominated, unlike in chapter 3 and 4 where the melting is ruled by pure conduction. In accord with constructal design, the search is for effective heat-flow architectures. The volume-constrained improvement of the designs for heat flow begins with assuming the simplest structure, where a single line serves as heat source. Next, the heat source is endowed with freedom to change its shape as it grows. The objective of the numerical simulations is to discover the geometric features that lead to the fastest melting process. The results show that the heat transfer rate density increases as the complexity and number of degrees of freedom of the structure are increased. Furthermore, the angles between heat invasion lines have a minor effect on the global performance compared to other degrees of freedom: number of branching levels, stem length, and branch lengths. The effect of natural convection in the melt zone is documented.

Numerical Investigation on the Heat Transfer Enhancement Using Micro/Nano Phase-Change Particulate Flow

The introduction of phase change material fluid and nanofluid in micro-channel heat sink design can significantly increase the cooling capacity of the heat sink because of the unique features of these two kinds of fluids. To better assist the design of a high performance micro-channel heat sink using phase change fluid and nanofluid, the heat transfer enhancement mechanism behind the flow with such fluids must be completely understood. A detailed parametric study is conducted to further investigate the heat transfer enhancement of the phase change material particle suspension flow, by using the two-phase non-thermal-equilibrium model developed by Hao and Tao (2004). The parametric study is conducted under normal conditions with Reynolds numbers of Re=600-900 and phase change material particle concentrations ¡Ü0.25 , as well as extreme conditions of very low Reynolds numbers (Re < 50) and high phase change material particle concentration (0.5-0.7) slurry flow. By using the two newly-defined parameters, named effectiveness factor and performance index, respectively, it is found that there exists an optimal relation between the channel design parameters, particle volume fraction, Reynolds number, and the wall heat flux. The influence of the particle volume fraction, particle size, and the particle viscosity, to the phase change material suspension flow, are investigated and discussed. The model was validated by available experimental data. The conclusions will assist designers in making their decisions that relate to the design or selection of a micro-pump suitable for micro or mini scale heat transfer devices. To understand the heat transfer enhancement mechanism of the nanofluid flow from the particle level, the lattice Boltzmann method is used because of its mesoscopic feature and its many numerical advantages. By using a two-component lattice Boltzmann model, the heat transfer enhancement of the nanofluid is analyzed, through incorporating the different forces acting on the nanoparticles to the two-component lattice Boltzmann model. It is found that the nanofluid has better heat transfer enhancement at low Reynolds numbers, and the Brownian motion effect of the nanoparticles will be weakened by the increase of flow speed.

Air-sea transfer of gas phase controlled compounds

Gases in the atmosphere/ocean have solubility that spans several orders of magnitude. Resistance in the molecular sublayer on the waterside limits the air-sea exchange of sparingly soluble gases such as SF6 and CO2. In contrast, both aerodynamic and molecular diffusive resistances on the airside limit the exchange of highly soluble gases (as well as heat). Here we present direct measurements of air-sea methanol and acetone transfer from two open cruises: the Atlantic Meridional Transect in 2012 and the High Wind Gas Exchange Study in 2013. The transfer of the highly soluble methanol is essentially completely airside controlled, while the less soluble acetone is subject to both airside and waterside resistances. Both compounds were measured concurrently using a proton-transfer-reaction mass spectrometer, with their fluxes quantified by the eddy covariance method. Up to a wind speed of 15 m s-1, observed air-sea transfer velocities of these two gases are largely consistent with the expected near linear wind speed dependence. Measured acetone transfer velocity is ~30% lower than that of methanol, which is primarily due to the lower solubility of acetone. From this difference we estimate the "zero bubble" waterside transfer velocity, which agrees fairly well with interfacial gas transfer velocities predicted by the COARE model. At wind speeds above 15 m s-1, the transfer velocities of both compounds are lower than expected in the mean. Air-sea transfer of sensible heat (also airside controlled) also appears to be reduced at wind speeds over 20 m s-1. During these conditions, large waves and abundant whitecaps generate large amounts of sea spray, which is predicted to alter heat transfer and could also affect the air-sea exchange of soluble trace gases. We make an order of magnitude estimate for the impacts of sea spray on air-sea methanol transfer.

Air-sea transfer of gas phase controlled compounds

Gases in the atmosphere/ocean have solubility that spans several orders of magnitude. Resistance in the molecular sublayer on the waterside limits the air-sea exchange of sparingly soluble gases such as SF6 and CO2. In contrast, both aerodynamic and molecular diffusive resistances on the airside limit the exchange of highly soluble gases (as well as heat). Here we present direct measurements of air-sea methanol and acetone transfer from two open cruises: the Atlantic Meridional Transect in 2012 and the High Wind Gas Exchange Study in 2013. The transfer of the highly soluble methanol is essentially completely airside controlled, while the less soluble acetone is subject to both airside and waterside resistances. Both compounds were measured concurrently using a proton-transfer-reaction mass spectrometer, with their fluxes quantified by the eddy covariance method. Up to a wind speed of 15 m s-1, observed air-sea transfer velocities of these two gases are largely consistent with the expected near linear wind speed dependence. Measured acetone transfer velocity is ~30% lower than that of methanol, which is primarily due to the lower solubility of acetone. From this difference we estimate the "zero bubble" waterside transfer velocity, which agrees fairly well with interfacial gas transfer velocities predicted by the COARE model. At wind speeds above 15 m s-1, the transfer velocities of both compounds are lower than expected in the mean. Air-sea transfer of sensible heat (also airside controlled) also appears to be reduced at wind speeds over 20 m s-1. During these conditions, large waves and abundant whitecaps generate large amounts of sea spray, which is predicted to alter heat transfer and could also affect the air-sea exchange of soluble trace gases. We make an order of magnitude estimate for the impacts of sea spray on air-sea methanol transfer.

The development of dynamic CFD model of the magnetic semi-levitation of liquid metals

In the casting of reactive metals, such as titanium alloys, contamination can be prevented if there is no contact between the hot liquid metal and solid crucible. This can be achieved by containing the liquid metal by means of high frequency AC magnetic field. A water cooled current-carrying coil, surrounding the metal can then provide the required Lorentz forces, and at the same time the current induced in the metal can provide the heating required to melt it. This ‘attractive’ processing solution has however many problems, the most serious being that of the control and containment of the liquid metal envelope, which requires a balance of the gravity and induced inertia forces on the one side, and the containing Lorentz and surface tension forces on the other. To model this process requires a fully coupled dyna ic solution of the flow fields, magnetic field and heat transfer/melding process to account for. A simplified solution has been published previously providing quasi-static solutions only, by taking the irrotational ‘magnetic pressure’ term of the Lorentz force into account. The authors remedy this deficiency by modelling the full problem using CFD techniques. The salient features of these techniques are included in this paper, as space allows.

Moisture Desorption Studies on Polymer Hydrated and Vacuum Extruded Bentonite Clay Mat

Moisture desorption observations from two bentonite clay mats subjected to ten environmental zones with individually different combinations of laboratory-controlled constant temperatures (between 20 °C and 40 °C) and relative humidity (between 15% and 70%) are presented. These laboratory observations are compared with predictions from mathematical models, such as thin-layer drying equations and kinetic drying models proposed by Page, Wang and Singh, and Henderson and Pabis. The quality of fit of these models is assessed using standard error (SE) of estimate, relative percent of error, and coefficient of correlation. The Page model was found to better predict the drying kinetics of the bentonite clay mats for the simulated tropical climates. Critical study on the drying constant and moisture diffusion coefficient helps to assess the efficacy of a polymer to retain moisture and control desorption through water molecule bonding. This is further substantiated with the Guggenheim–Anderson–De Boer (GAB) desorption isotherm model which is presented.