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

A direct method for the Boltzmann equation based on a pseudo-spectral velocity space discretization

A deterministic method is proposed for solving the Boltzmann equation. The method employs a Galerkin discretization of the velocity space and adopts, as trial and test functions, the collocation basis functions based on weights and roots of a Gauss-Hermite quadrature. This is defined by means of half- and/or full-range Hermite polynomials depending whether or not the distribution function presents a discontinuity in the velocity space. The resulting semi-discrete Boltzmann equation is in the form of a system of hyperbolic partial differential equations whose solution can be obtained by standard numerical approaches. The spectral rate of convergence of the results in the velocity space is shown by solving the spatially uniform homogeneous relaxation to equilibrium of Maxwell molecules. As an application, the two-dimensional cavity flow of a gas composed by hard-sphere molecules is studied for different Knudsen and Mach numbers. Although computationally demanding, the proposed method turns out to be an effective tool for studying low-speed slightly rarefied gas flows

Solving the Boltzmann Equation on GPU

We show how to accelerate the direct solution of the Boltzmann equation using Graphics Processing Units (GPUs). In order to fully exploit the computational power of the GPU, we choose a method of solution which combines a finite difference discretization of the free-streaming term with a Monte Carlo evaluation of the collision integral. The efficiency of the code is demonstrated by solving the two-dimensional driven cavity flow. Computational results show that it is possible to cut down the computing time of the sequential code of two order of magnitudes. This makes the proposed method of solution a viable alternative to particle simulations for studying unsteady low Mach number flows.Comment: 18 pages, 3 pseudo-codes, 6 figures, 1 tabl

Rarefied flow between plates of finite length via the coupling approach

This paper was presented at the 3rd Micro and Nano Flows Conference (MNF2011), which was held at the Makedonia Palace Hotel, Thessaloniki in Greece. The conference was organised by Brunel University and supported by the Italian Union of Thermofluiddynamics, Aristotle University of Thessaloniki, University of Thessaly, IPEM, the Process Intensification Network, the Institution of Mechanical Engineers, the Heat Transfer Society, HEXAG - the Heat Exchange Action Group, and the Energy Institute.The coexistence of rarefied continuum flow regime areas and relatively small elements in which rarefaction effects become important is a typical feature of many complex gas flows micro systems. In rarefied domains, the mean free path of gas molecules is comparable or larger than a characteristic scale of the system. These domains are naturally described by kinetic equation for the velocity distribution function, which involve a considerable effort in terms of CPU time and memory requirements, due to the discretization in both physical and velocity space. The continuum domains are best described by the fluid Navier Stokes (NS) equations in terms of average flow velocity, gas density and temperature. These equations are more efficient, but less accurate in critical rarefied areas. Thus, the development of hybrid solver combining kinetic and continuum models is of great interest especially for applications range from gas flows in micro systems to the aerospace applications, such as high altitude flights. The pressure–driven gas flow of rarified monatomic gas through a two-dimensional short microchannel is considered using hybrid solver. The calculations have been carried out for pressure ratios 0.1, 0.5 and 0.9 and fixed relatively large Knudsen number. The applicability of the solver is discussed via comparison with the kinetic and NS solutions.The European Community's Seventh Framework Programme FP7/2007-2013 under grant agreement ITN GASMEMS no 215504