Microchannel fluid flow and heat transfer by lattice boltzmann method

This paper was presented at the 4th Micro and Nano Flows Conference (MNF2014), which was held at University College, London, UK. The conference was organised by Brunel University and supported by the Italian Union of Thermofluiddynamics, IPEM, the Process Intensification Network, the Institution of Mechanical Engineers, the Heat Transfer Society, HEXAG - the Heat Exchange Action Group, and the Energy Institute, ASME Press, LCN London Centre for Nanotechnology, UCL University College London, UCL Engineering, the International NanoScience Community, www.nanopaprika.eu.Micro flow has become a popular field of interest due to the advent of micro electromechanical systems (MEMS). In this work, the lattice Boltzmann method, a particle-based approach, is applied to simulate the two-dimensional micro channel fluid flow. We simulated fluid flow and heat transfer inside microchannel, the prototype application of this study is micro-heat exchangers. The main incentive to look at fluidic behaviour at micron scale is that micro devices tend to behave much differently from the objects we are used to handling in daily life. The choice of using LBM for micro flow simulation is a good one owing to the fact that it is based on the Boltzmann equation which is valid for the whole range of the Knudsen number. Slip velocity and temperature jump boundary conditions are used for the microchannel simulations with Knudsen number values covering the slip flow. The lattice Bhatnagar-Gross-Krook single relaxation time approximation was used. The results found are compared with the Navier-Stokes analytical and numerical results available in the literature and good matches are observed

Simulating fluid flows in micro and nano devices : the challenge of non-equilibrium behaviour

We review some recent developments in the modelling of non-equilibrium (rarefied) gas flows at the micro- and nano-scale, concentrating on two different but promising approaches: extended hydrodynamic models, and lattice Boltzmann methods. Following a brief exposition of the challenges that non-equilibrium poses in micro- and nano-scale gas flows, we turn first to extended hydrodynamics, outlining the effective abandonment of Burnett-type models in favour of high-order regularised moment equations. We show that the latter models, with properly-constituted boundary conditions, can capture critical non-equilibrium flow phenomena quite well. We then review the boundary conditions required if the conventional Navier-Stokes-Fourier (NSF) fluid dynamic model is applied at the micro scale, describing how 2nd-order Maxwell-type conditions can be used to compensate for some of the non-equilibrium flow behaviour near solid surfaces. While extended hydrodynamics is not yet widely-used for real flow problems because of its inherent complexity, we finish this section with an outline of recent 'phenomenological extended hydrodynamics' (PEH) techniques-essentially the NSF equations scaled to incorporate non-equilibrium behaviour close to solid surfaces-which offer promise as engineering models. Understanding non-equilibrium within lattice Boltzmann (LB) framework is not as advanced as in the hydrodynamic framework, although LB can borrow some of the techniques which are being developed in the latter-in particular, the near-wall scaling of certain fluid properties that has proven effective in PEH. We describe how, with this modification, the standard 2nd-order LB method is showing promise in predicting some rarefaction phenomena, indicating that instead of developing higher-order off-lattice LB methods with a large number of discrete velocities, a simplified high-order LB method with near-wall scaling may prove to be just as effective as a simulation tool

A thermal lattice Boltzmann model for micro/nano-flows

The dynamic behavior of charged micro and nanofluids plays a crucial role in a large variety of industrial and biological processes. Such dynamic behavior is characterized by the simultaneous occurrence of several competing mechanisms, such as electrostatic interactions, viscous dissipation and hydrodynamic effects, often taking place in complex geometries. This paper focuses on a thermal lattice Boltzmann model for micro/nano-flows

Lattice ellipsoidal statistical BGK model for thermal non-equilibrium flows

A thermal lattice Boltzmann model is constructed on the basis of the ellipsoidal statistical Bhatnagar-Gross-Krook (ES-BGK) collision operator via the Hermite moment representation. The resulting lattice ES-BGK model uses a single distribution function and features an adjustable Prandtl number. Numerical simulations show that using a moderate discrete velocity set, this model can accurately recover steady and transient solutions of the ES-BGK equation in the slip-flow and early transition regimes in the small Mach number limit that is typical of microscale problems of practical interest. In the transition regime in particular, comparisons with numerical solutions of the ES-BGK model, direct Monte Carlo and low-variance deviational Monte Carlo simulations show good accuracy for values of the Knudsen number up to approximately 0:5. On the other hand, highly non-equilibrium phenomena characterized by high Mach numbers, such as viscous heating and force-driven Poiseuille flow for large values of the driving force, are more difficult to capture quantitatively in the transition regime using discretizations that have been chosen with computational efficiency in mind such as the one used here, although improved accuracy is observed as the number of discrete velocities is increased

Visual Simulation of Flow

We have adopted a numerical method from computational fluid dynamics, the Lattice Boltzmann Method (LBM), for real-time simulation and visualization of flow and amorphous phenomena, such as clouds, smoke, fire, haze, dust, radioactive plumes, and air-borne biological or chemical agents. Unlike other approaches, LBM discretizes the micro-physics of local interactions and can handle very complex boundary conditions, such as deep urban canyons, curved walls, indoors, and dynamic boundaries of moving objects. Due to its discrete nature, LBM lends itself to multi-resolution approaches, and its computational pattern, which is similar to cellular automata, is easily parallelizable. We have accelerated LBM on commodity graphics processing units (GPUs), achieving real-time or even accelerated real-time on a single GPU or on a GPU cluster. We have implemented a 3D urban navigation system and applied it in New York City with real-time live sensor data. In addition to a pivotal application in simulation of airborne contaminants in urban environments, this approach will enable the development of other superior prediction simulation capabilities, computer graphics and games, and a novel technology for computational science and engineering

Accuracy analysis of high-order lattice Boltzmann models for rarefied gas flows

In this work, we have theoretically analyzed and numerically evaluated the accuracy of high-order lattice Boltzmann (LB) models for capturing non-equilibrium effects in rarefied gas flows. In the incompressible limit, the LB equation is shown to be able to reduce to the linearized Bhatnagar–Gross–Krook (BGK) equation. Therefore, when the same Gauss–Hermite quadrature is used, LB method closely resembles the discrete velocity method (DVM). In addition, the order of Hermite expansion for the equilibrium distribution function is found not to be directly correlated with the approximation order in terms of the Knudsen number to the BGK equation for incompressible flows. Meanwhile, we have numerically evaluated the LB models for a standing-shear-wave problem, which is designed specifically for assessing model accuracy by excluding the influence of gas molecule/surface interactions at wall boundaries. The numerical simulation results confirm that the high-order terms in the discrete equilibrium distribution function play a negligible role in capturing non-equilibrium effect for low-speed flows. By contrast, appropriate Gauss–Hermite quadrature has the most significant effect on whether LB models can describe the essential flow physics of rarefied gas accurately. Our simulation results, where the effect of wall/gas interactions is excluded, can lead to conclusion on the LB modeling capability that the models with higher-order quadratures provide more accurate results. For the same order Gauss–Hermite quadrature, the exact abscissae will also modestly influence numerical accuracy. Using the same Gauss–Hermite quadrature, the numerical results of both LB and DVM methods are in excellent agreement for flows across a broad range of the Knudsen numbers, which confirms that the LB simulation is similar to the DVM process. Therefore, LB method can offer flexible models suitable for simulating continuum flows at the Navier–Stokes level and rarefied gas flows at the linearized Boltzmann model equation level

DSMC-LBM mapping scheme for rarefied and non-rarefied gas flows

We present the formulation of a kinetic mapping scheme between the Direct Simulation Monte Carlo (DSMC) and the Lattice Boltzmann Method (LBM) which is at the basis of the hybrid model used to couple the two methods in view of efficiently and accurately simulate isothermal flows characterized by variable rarefaction effects. Owing to the kinetic nature of the LBM, the procedure we propose ensures to accurately couple DSMC and LBM at a larger Kn number than usually done in traditional hybrid DSMC-Navier-Stokes equation models. We show the main steps of the mapping algorithm and illustrate details of the implementation. Good agreement is found between the moments of the single particle distribution function as obtained from the mapping scheme and from independent LBM or DSMC simulations at the grid nodes where the coupling is imposed. We also show results on the application of the hybrid scheme based on a simpler mapping scheme for plane Poiseuille flow at finite Kn number. Potential gains in the computational efficiency assured by the application of the coupling scheme are estimated for the same flow.Comment: Submitted to Journal of Computational Scienc