Lattice Boltzmann in micro- and nano- flow simulations

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.One of the fundamental difficulties in micro- and nano-flow simulations is that the validity’s of the continuum assumption and the hydro-dynamic equations start to become questionable in this flow regime. The lower-level kinetic/molecular alternatives are often either prohibitively expensive for practical purposes or poorly justified from a fundamental perspective. The lattice Boltzmann (LB) method, which originated from a simplistic Boolean kinetic model, is recently shown to converge asymptotically to the continuum Boltzmann-BGK equation and therefore offers a theoretically sound and computationally effective approach for micro- and nano-flow simulations. In addition, its kinetic nature allows certain microscopic physics to be modeled at the macroscopic level, leading to a highly efficient model for multiphase flows with phase transitions. With the inherent computational advantages of a lattice model, e.g., the algorithm simplicity and parallelizability, the ease of handling complex geometry and so on, the LB method has found many applications in various areas of Computational Fluid Dynamics (CFD) and matured to the extend of commercial applications. In this talk, I shall give an introduction to the LB method with the emphasis given to the theoretical justifications for its applications in micro- and nano-flow simulations. Some recent examples will also be reported

Mechanics of blood flow in capillaries

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.Blood is a concentrated suspension of red blood cells (RBCs). Motion and deformation of RBCs can be analyzed based on knowledge of their mechanical characteristics. Models for single-file motion of RBCs in capillaries yield predictions of apparent viscosity in good agreement with experimental results for diameters up to about 8 μm. In living microvessels, flow resistance is also strongly influenced by the presence of a ~ 1-micron layer of macromolecules bound to the inner lining of vessel walls, the endothelial surface layer. Two-dimensional simulations, in which each RBC is represented as a set of interconnected viscoelastic elements, predict that off-center RBCs take asymmetric shapes and drift toward the center-line. Predicted trajectories agree closely with observations in microvessels of the rat mesentery. Realistic simulation of multiple interacting RBCs in microvessels remains as a major challenge for future work.This work was supported by NIH Grant HL034555

Similarities and differences between flow boiling in microchannels and pool boiling

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.Recent literature indicates that under certain conditions the heat transfer coefficient during flow boiling in microchannels is quite similar to that under pool boiling conditions. This is rather unexpected as microchannels are believed to provide significant heat transfer enhancement under single-phase as well as flow boiling conditions. This paper explores the underlying heat transfer mechanisms and illustrates the similarities and differences between the two processes. Formation of elongated bubbles and their passage over the microchannel walls have similarities to the bubble ebullition cycle in pool boiling. During the passage of elongated bubbles, the longer duration between two successive liquid slugs leads to wall dryout and a critical heat flux that may be lower than that under pool boiling conditions. A clear understanding of the similarities and differences will help in overcoming some of these limiting factors and in developing strategies for enhancing heat transfer during flow boiling in microchannels

Multifactoring concept – A key to investigation of forced-eoiling in microsystems

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.In the present paper forced-boiling in microsystems is considered in the light of fundamentals of boiling heat transfer such as local temperature pulsations of heating surface (Moore and Mesler, 1961), pumping effect of growing bubble (PEGB)(Shekriladze, 1966), a model of “the theatre of director” (MTD) (Shekriladze and Ratiani, 1966) and multifactoring concept (MFC) (Shekriladze, 2006). An attempt is made to resolve a contradiction between accordance of heat transfer process to developed boiling heat transfer law in the major part of experiments and qualitatively differing trends in the other part of processes. The problem of interpretation of generation of strong reverse vapor flows, related cyclical oscillations and flow instabilities also is touched. According to presented analysis leading role in specific thermo-hydrodynamic characteristics of boiling microsystems is played by so-called duration-dependent multifactoring which, by its part, is linked to transition to prolonged action of microlayer evaporation (MLE) and PEGB. As a result drastically increases a number of influencing heat transfer factors extremely complicating description of the process. At the same time prolongation of intensive stage of acting of MLE and PEGB creates prerequisites for specific thermo-hydrodynamic appearances

On emerging micro- and nanoscale thermofluidic technologies

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.This paper highlights examples of my current research in heat transfer and fluidics at the interface of energy applications and micro- and nanoscale technologies. It is not the scope of this paper to present an exhaustive account of all current and past activities related to its title. It is rather an account of current research in my laboratory in this area, containing both the underlying scientific challenges as well as the hoped final outcome in terms of applications. To this end, examples from the areas of energy conversion, as well as energy transport will be discussed. In the area of energy conversion an original, deformable, direct methanol microfuel cell will be presented made of lightweight, flexible, polymer-based materials. A basic understanding and control of two-phase flows (in this case methanol and carbon dioxide) in microchannels as well as novel materials processing and microfabrication methods are directly related to the performance of such energy conversion devices. In the area of energy conservation and reuse, examples from the information technology are employed. Specifically, new concepts of liquid (water) cooling of chips reaching heat removal rates in excess of 700 W/cm2 in domains with restricted heights of the order of one mm will be presented. One additional advantage of using water to cool high density electronics is energy reuse, due to the potentially much higher exergy content of the coolant compared to air cooled technologies. The last part of the paper focuses on the employment of functional nanostructures such as carbon nanotubes and nanowires of conductive and semiconductive materials for the efficient transport of electricity and heat and the need for the development of novel technologies for the manufacturing, characterization as well as handling of such nanostructures

A way to visualise heat transfer in 3D unsteady flows

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.Heat transfer in fluid flows traditionally is examined in terms of temperature field and heat-transfer coefficients. However, heat transfer may alternatively be considered as the transport of thermal energy by the total convective-conductive heat flux in a way analogous to the transport of fluid by the flow field. The paths followed by the total heat flux are the thermal counterpart to fluid trajectories and facilitate heat-transfer visualisation in a similar manner as flow visualisation. This has great potential for applications in which insight into the heat fluxes throughout the entire configuration is essential (e.g. cooling systems, heat exchangers). To date this concept has been restricted to 2D steady flows. The present study proposes its generalisation to 3D unsteady flows by representing heat transfer as the 3D unsteady motion of a virtual fluid subject to continuity. The heat transfer visualisation is provided with a physical framework and demonstrated by way of representative examples. Furthermore, a fundamental analogy between fluid motion and heat transfer is addressed that may pave the way to future heat-transfer studies by well-established geometrical methods from laminar-mixing studies

Instabilities in free-surface electroosmotic flows

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.With the recent development of novel microfluidic devices electroosmotic flows with fluid/fluid interfaces have emerged as very important subjects of investigation. Two immiscible fluids may need to be transported in a microchannel, or one side of a channel may be open to air for various purposes, including adsorption of airborne molecules to liquid for high-sensitivity substance detection. The liquid/liquid or liquid/gas interface in these cases can deform, resulting in significant corrugations followed sometimes by incipient rupture of liquid layers. For electroosmotic flow the rupture, leading to shortcircuit, can cause overall failure of the device. It is thus imperative to know the conditions for the rupture as well as the initial interfacial instability. Studies based on the Debye-Huckle approximation reveal that all free-surface electroosmotic flows of thickness larger than the Debye screening length are unstable and selectively lead to rupture. Layers of the order of Debye screening length, however, are not properly described by the Debye-Huckle approximation. Even for micro-scale layers, the rupture phenomenon can make local layer thickness to be nanoscale. A fully coupled system of hydrodynamics, electric field, and ionic distribution need to be analyzed. In this paper linear instability and subsequent nonlinear developments of a nanoscale free-surface electroosmotic flow are reported.This study is sponsored by the Ministry of Education, Science and Technology of Korea through the World Class University Grant

Variational approach to gas flows in microchannels on the basis of the Boltzmann equation for hard-sphere molecules

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.The objective of the present paper is to provide an analytic expression for the first- and second-order velocity slip coefficients. Therefore, gas flow rates in microchannels have been rigorously evaluated in the near-continuum limit by means of a variational technique which applies to the integrodifferential form of the Boltzmann equation based on the true linearized collision operator. The diffuse-specular reflection condition of Maxwell’s type has been considered in order to take into account the influence of the accommodation coefficient on the slip parameters. The polynomial form of Knudsen number obtained for the Poiseuille mass flow rate and the values of the second order velocity slip coefficients found on the basis of our variational solution of the linearized Boltzmann equation for hardsphere molecules are analyzed in the frame of potential applications of classical continuum numerical tools (as lattice Boltzmann methods) in simulations of microscale flows

Local heating effects on flow and heat transfer in microchannels

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.A series of numerical investigations was conducted to explore the effects of temperature-dependent viscosity and thermal conductivity on two-dimensional low Reynolds number convection of water in microchannels with locally heating. An emphasis was addressed on the fundamental characteristics of flow and thermal re-development at different localized heat fluxes and different inlet temperatures. The velocity field is highly coupled with temperature distribution and distorted through the variations of viscosity and thermal conductivity. The induced cross-flow velocity has a marked contribution to the convection. The heat transfer enhancement due to viscosity-variation is pronounced, though the axial convection introduced by thermal-conductivity-variation is insignificant unless for the cases of very low Reynolds numbers. The heat transfer enhancement is described by defining the peak value and location of relative Nusselt number distribution as ΔNu%max and Xmax. Strong nonlinear interaction mechanism prevails in the correlation of ΔNu%max and Xmax due to high heat flux condition and dramatic rise of liquid temperature.This study is supported by the National Natural Science Foundation of China (Grant No. 50636030)

Numerical modelling of heat generated by electroosmotic flows in micro-channels

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.In this paper, numerical modeling of Joule heating in electroosmotic flows is described in some detail. The finite element method is used for the spatial discretization along with the characteristic based split (CBS) time discretization. A new non-dimensional scaling is also introduced. In addition to standard problems of micro channel flows, flow and heat generation in a T-mixer are also discussed in this paper