Optimization of three-dimensional branching networks of microchannels for thermal management of microelectronics

The aim of this work is to present a methodology to develop cost-effective thermal management solutions for microelectronic devices, capable of removing maximum amount of heat and delivering maximally uniform temperature distributions. The topological and geometrical characteristics of multiple-story three-dimensional branching networks of microchannels were developed using multi-objective optimization. A conjugate heat transfer analysis software package and an automatic 3D microchannel network generator were developed and coupled with a modified version of a particle-swarm optimization algorithm with a goal of creating a design tool for 3D networks of optimized coolant flow passages. Numerical algorithms in the conjugate heat transfer solution package include a quasi-ID thermo-fluid solver and a steady heat diffusion solver, which were validated against results from high-fidelity Navier-Stokes equations solver and analytical solutions for basic fluid dynamics test cases. Pareto-optimal solutions demonstrate that thermal loads of up to 500 W/cm2 can be managed with 3D microchannel networks, with pumping power requirements up to 50% lower with respect to currently used high-performance cooling technologies.

In vitro assessment of migratory behavior of two cell populations in a simple multichannel microdevice

Recent literature suggests that mesenchymal stem/stromal cells (MSC) could be used as Trojan Horses to deliver “death-signals” to cancer cells. Herein, we describe the development of a novel multichannel cell migration device, and use it to investigate the relative migration rates of bone marrow-derived MSC and breast cancer cells (MCF-7) towards each other. Confluent monolayers of MSC and MCF-7 were established in adjacent chambers separated by an array of 14 microchannels. Initially, culture chambers were isolated by air bubbles (air-valves) contained within each microchannel, and then bubbles were displaced to initiate the assay. The MCF-7 cells migrated preferentially towards MSC, whilst the MSC did not migrate preferentially towards the MCF-7 cells. Our results corroborate previous literature that suggests MSC migration towards cancer cells in vivo is in response to the associated inflammation rather than directly to signals secreted by the cancer cells themselves.

Optimisation methods for coupled thermodynamic and 1D design of radial-inflow turbines

Optimisation is a fundamental step in the turbine design process, especially in the development of non-classical designs of radial-inflow turbines working with high-density fluids in low-temperature Organic Rankine Cycles (ORCs). The present work discusses the simultaneous optimisation of the thermodynamic cycle and the one-dimensional design of radial-inflow turbines. In particular, the work describes the integration between a 1D meanline preliminary design code adapted to real gases and the performance estimation approach for radial-inflow turbines in an established ORC cycle analysis procedure. The optimisation approach is split in two distinct loops; the inner operates on the 1D design based on the parameters received from the outer loop, which optimises the thermodynamic cycle. The method uses parameters including brine flow rate, temperature and working fluid, shifting assumptions such as head and flow coefficients into the optimisation routine. The discussed design and optimisation method is then validated against published benchmark cases. Finally, using the same conditions, the coupled optimisation procedure is extended to the preliminary design of a radial-inflow turbine with R143a as working fluid in realistic geothermal conditions and compared against results from commercially-available software RITAL from Concepts-NREC.

Time dependence of effective slip on textured hydrophobic surfaces

In this paper, we present results on water flow past randomly textured hydrophobic surfaces with relatively large surface features of the order of 50 µm. Direct shear stress measurements are made on these surfaces in a channel configuration. The measurements indicate that the flow rates required to maintain a shear stress value vary substantially with water immersion time. At small times after filling the channel with water, the flow rates are up to 30% higher compared with the reference hydrophilic surface. With time, the flow rate gradually decreases and in a few hours reaches a value that is nearly the same as the hydrophilic case. Calculations of the effective slip lengths indicate that it varies from about 50 µm at small times to nearly zero or “no slip” after a few hours. Large effective slip lengths on such hydrophobic surfaces are known to be caused by trapped air pockets in the crevices of the surface. In order to understand the time dependent effective slip length, direct visualization of trapped air pockets is made in stationary water using the principle of total internal reflection of light at the water-air interface of the air pockets. These visualizations indicate that the number of bright spots corresponding to the air pockets decreases with time. This type of gradual disappearance of the trapped air pockets is possibly the reason for the decrease in effective slip length with time in the flow experiments. From the practical point of usage of such surfaces to reduce pressure drop, say, in microchannels, this time scale of the order of 1 h for the reduction in slip length would be very crucial. It would ultimately decide the time over which the surface can usefully provide pressure drop reductions. ©2009 American Institute of Physics

Underwater Sustainability of the ``Cassie'' State of Wetting

A rough hydrophobic surface when immersed in water can result in a ``Cassie'' state of wetting in which the water is in contact with both the solid surface and the entrapped air. The sustainability of the entrapped air on such surfaces is important for underwater applications such as reduction of flow resistance in microchannels and drag reduction of submerged bodies such as hydrofoils. We utilize an optical technique based oil total internal reflection of light at the water-air interface to quantify the spatial distribution of trapped air oil such a surface and its variation with immersion time. With this technique, we evaluate the sustainability of the Cassie state on hydrophobic surfaces with four different kinds of textures. The textures studied are regular arrays of pillars, ridges, and holes that were created in silicon by a wet etching technique, and also a texture of random craters that was obtained through electrodischarge machining of aluminum. These surfaces were rendered hydrophobic with a self-assembled layer Of fluorooctyl trichlorosilane. Depending on the texture, the size and shape of the trapped air pockets were found to vary. However, irrespective of the texture, both the size and the number of air pockets were found to decrease with time gradually and eventually disappear, suggesting that the sustainability of the ``Cassie'' state is finite for all the microstructures Studied. This is possibly due to diffusion of air from the trapped air pockets into the water. The time scale for disappearance of air pockets was found to depend on the kind of microstructure and the hydrostatic pressure at the water-air interface. For the surface with a regular array of pillars, the air pockets were found to be in the form of a thin layer perched on top of the pillars with a large lateral extent compared to the spacing between pillars. For other surfaces studied, the air pockets are smaller and are of the same order as the characteristic length scale of the texture. Measurements for the surface with holes indicate that the time for air-pocket disappearance reduces as the hydrostatic pressure is increased.

A thermogravimetric study of the oxidative growth of Al2O3/Al alloy composites

The oxidation of liquid Al–Mg–Si alloys at 900–1400 °C was studied by thermogravimetric analysis (TGA). The development of a semi-protective surface layer of MgO/MgAl2O4 allows the continuous formation of an Al2O3-matrix composite containing an interpenetrating network of metal microchannels at 1000–1350 °C. An initial incubation period precedes bulk oxidation, wherein Al2O3 grows from a near-surface alloy layer by reaction of oxygen supplied by the dissolution of the surface oxides and Al supplied from a bulk alloy reservoir through the microchannel network. The typical oxidation rate during bulk growth displays an initial acceleration followed by a parabolic deceleration in a regime apparently limited by Al transport to the near-surface layer. Both regimes may be influenced by the Si content in this layer, which rises due to preferential Al and Mg oxidation. The growth rates increase with temperature to a maximum at ~1300 °C, with a nominal activation energy of 270 kJ/mole for an Al-2.85 wt. % Mg-5.4 wt. % Si alloy in O2 at furnace temperatures of 1000–1300 °C. An oscillatory rate regime observed at 1000–1075 °C resulted in a banded structure of varying Al2O3-to-metal volume fraction.

Nucleation and growth of Al2O3/metal composites by oxidation of aluminum alloys

The nucleation and growth mechanisms during high temperature oxidation of liquid Al-3% Mg and Al-3% Mg-3% Si alloys were studied with the aim of enhancing our understanding of a new composite fabrication process. The typical oxidation sequence consists of an initial event of rapid but brief oxidation, followed by an incubation period of limited oxide growth after which bulk Al2O3/Al composite forms. A duplex oxide layer, MgO (upper) and MgAl2O4 (lower), forms on the alloy surface during initial oxidation and incubation. The spinel layer remains next to the liquid alloy during bulk oxide growth and is the eventual repository for most of the magnesium in the original alloy. Metal microchannels developed during incubation continuously supply alloy through the composite to the reaction interface. During the growth process, a layered structure exists at the upper extremity of the composite, consisting of MgO at the top surface, MgAl2O4 (probably discontinuous), Al alloy, and finally the bulk Al2O3 composite containing microchannels of the alloy. The bulk oxide growth mechanism appears to involve continuous formation and dissolution of the Mg-rich oxides at the surface, diffusion of oxygen through the underlying liquid metal, and epitaxial growth of Al2O3 on the existing composite body. The roles of Mg and Si in the composite growth process are discussed.

Kinetic Assessment of Measured Mass Flow Rates and Streamwise Pressure Distributions in Microchannel Gas Flows

Measured mass flow rates and streamwise pressure distributions of gas flowing through microchannels were reported by many researchers. Assessment of these data is crucial before they are used in the examination of slip models and numerical schemes, and in the design of microchannel elements in various MEMS devices. On the basis of kinetic solutions of the mass flow rates and pressure distributions in microchannel gas flows, the measured data available are properly normalized and then are compared with each other. The 69 normalized data of measured pressure distributions are in excellent agreement, and 67 of them are within 1 +/- 0.05. The normalized data of mass flow-rates ranging between 0.95 and 1 agree well with each other as the inlet Knudsen number Kn (i) < 0.02, but they scatter between 0.85 and 1.15 as Kn (i) > 0.02 with, to some extent, a very interesting bifurcation trend.

磁性纳米颗粒对微槽道液体流动的混合增强作用

为了改进微流控芯片中液体混合效率,尝试在一种液体中加入磁性纳米颗粒,增强与另外一种液体混合.在外加磁场作用下,对Y型微槽道中磁性与非磁性液体流动混合进行了实验观测,所需混合长度较无磁场时缩短了约2个数量级.利用CFD软件进行数值模拟,验证了外磁场作用下纳米磁性液体对流动混合有着明显的强化作用.

Experiments about Simple Liquid Flows in Microtubes

The behavior of micro-scale flow is significant for the performance of Micro-Electro-Mechanical- Systems (MEMS) devices. Some experiments about liquid flow through microtubes with diameters about 3similar to20mum are presented here. The liquids used in our experiments include some simple liquids with small molecules, such as non-ion water and several kinds of organic liquids (CCL4, C6H5C2H5 and Isopropanol etc.). The flow rate and the normalized friction cocfficients were measured in micro-flow experimental apparatus. The results show that when the driven pressure varies from 0 to 1Mpa, the flow behaviors in 20mum microtube for both polar and non-polar liquids are in agreement with Hagen-Poiseuille law of the classical theory. It means that N-S equation based on continuous medium still acts well in this case. For higher pressure drop from 1 to 30Mpa, in the microtubes with diameter of 3similar to10mum, the normalized friction coefficients of organic liquids can't keep constant with pressure increases. However the non-ion water reveals different trends.

Microfluidic Rheometry

The development and growth of microfluidics has stimulated interest in the behaviour of complex liquids in micro-scale geometries and provided a rich platform for rheometric investigations of non-Newtonian phenomena at small scales. Microfluidic techniques present the rheologist with new opportunities for material property measurement and this review discusses the use of microfluidic devices to measure bulk rheology in both shear and extensional flows. Capillary, stagnation and contraction flows are presented in this context and developments, limitations and future perspectives are examined. (C) 2008 Elsevier Ltd. All rights reserved.

An Atomistic-Continuum Hybrid Simulation Of Fluid Flows Over Superhydrophobic Surfaces

Recent experiments have found that slip length could be as large as on the order of 1 mu m for fluid flows over superhydrophobic surfaces. Superhydrophobic surfaces can be achieved by patterning roughness on hydrophobic surfaces. In the present paper, an atomistic-continuum hybrid approach is developed to simulate the Couette flows over superhydrophobic surfaces, in which a molecular dynamics simulation is used in a small region near the superhydrophobic surface where the continuum assumption is not valid and the Navier-Stokes equations are used in a large region for bulk flows where the continuum assumption does hold. These two descriptions are coupled using the dynamic coupling model in the overlap region to ensure momentum continuity. The hybrid simulation predicts a superhydrophobic state with large slip lengths, which cannot be obtained by molecular dynamics simulation alone.

An atomistic-continuum hybrid simulation of fluid flows over superhydrophobic surfaces

Recent experiments have found that slip length could be as large as on the order of 1 mu m for fluid flows over superhydrophobic surfaces. Superhydrophobic surfaces can be achieved by patterning roughness on hydrophobic surfaces. In the present paper an atomistic-continuum hybrid approach is developed to simulate the Couette flows over superhydrophobic surfaces in which a molecular dynamics simulation is used in a small region near the superhydrophobic surface where the continuum assumption is not valid and the Navier-Stokes equations are used in a large region for bulk flows where the continuum assumption does hold. These two descriptions are coupled using the dynamic coupling model in the overlap region to ensure momentum continuity. The hybrid simulation predicts a superhydrophobic state with large slip lengths which cannot be obtained by molecular dynamics simulation alone.

Solute transport in nanochannels with roughness-like structures

Newfound attention has been given to solute transport in nanochannels. Because the electric double layer (EDL) thickness is comparable to characteristic channel dimensions, nanochannels have been used to separate ionic species with a constant charge-to-size ratio (i.e., electrophoretic mobility) that otherwise cannot be separated in electroosmotic or pressure- driven flow along microchannels. In nanochannels, the electrical fields within the EDL cause transverse ion distributions and thus yield charge-dependent mean ion speeds in the flow. Surface roughness is usually inevitable during microfabrication of microchannels or nanochannels. Surface roughness is usually inevitable during the fabrication of nanochannels. In the present study, we develop a numerical model to investigate the transport of charged solutes in nanochannels with hundreds of roughness-like structures. The model is based on continuum theory that couples Navier-Stokes equations for flows, Poisson-Boltzmann equation for electrical fields, and Nernst-Planck equation for solute transports. Different operating conditions are considered and the solute transport patterns in rough channels are compared with those in smooth channels. Results indicate that solutes move slower in rough nanochannels than in smooth ones for both pressure- driven and electroosmotic flows. Moreover, solute separation can be significantly improved by surface roughness under certain circumstances.

The ionic exclusion-enrichment in a hybrid micro-/nano-channel

An ionic exclusion-enrichment phenomenon has been found at the ends of a nano-channel when electric-driven fluid passes through a micro-/nano-hybrid channel [1-3]. In our experiments, the hybrid channels are fabricated with two poly-dimethysiloxane (PDMS) monoliths microchannels (100um X20um X 9mm) and a nanoporous polycarbonate nuclear track-etched (PCTE) membrane (with 50nm pores). The flows are driven under different electrical potential and the test liquids with different PH values are used. The ion depletion in the source channel is observed by the MicroPIV system. In addition, the numerical simulations about ionic exclusion-enrichment in the hybrid channel are carried out. Some results are as followed: