Lattice Boltzmann methods

Lattice Boltzmann methods

Lattice Boltzmann methods (LBM) is a class of computational fluid dynamics (CFD) methods for fluid simulation. Instead of solving the Navier–Stokes equations, the discrete Boltzmann equation is solved to simulate the flow of a Newtonian fluid with collision models such as Bhatnagar-Gross-Krook (BGK). By simulating the interaction of a limited number of particles the viscous flow behavior emerges automatically from the intrinsic particle streaming and collision processes.

Algorithm

LBM is a relatively new simulation technique for complex fluid systems and has attracted interest from researchers in computational physics. Unlike the traditional CFD methods, which solve the conservation equations of macroscopic properties (i.e., mass, momentum, and energy) numerically, LBM models the fluid consisting of fictive particles, and such particles perform consecutive propagation and collision processes over a discrete lattice mesh. Due to its particulate nature and local dynamics, LBM has several advantages over other conventional CFD methods, especially in dealing with complex boundaries, incorporating of microscopic interactions, and parallelization of the algorithm.

Development from the LGA method

LBM originated from the lattice gas automata (LGA) method, which can be considered as a simplified fictitious molecular dynamics model in which space, time, and particle velocities are all discrete. Each lattice node is connected to its neighbors by 6 lattice velocities, using for example the hexagonal FHP model. There can be either 0 or 1 particles at a lattice node moving in a lattice direction. After a time interval, each particle will move to the neighboring node in its direction, this process is called the propagation or streaming step. When there are more than one particles arriving at the same node from different directions, they collide and change their directions according to a set of collision rules. Suitable collision rules should conserve the particle number (mass), momentum, and energy before and after the collision. However, it was found that LGA suffers from several native defects: lack of Galilean invariance, statistical noise, exponential complexity for three-dimensional lattices, etc.

The main motivation for the transition from LGA to LBM was the desire to remove the statistical noise by replacing the Boolean particle number in a lattice direction with its ensemble average, the so-called density distribution function. Accompanying this replacement, the discrete collision rule is also replaced by a continuous function known as the collision operator. In the LBM development, an important simplification is to approximate the collision operator with the Bhatnagar-Gross-Krook (BGK) relaxation term. This lattice BGK (LBGK) model makes simulations more efficient and allows flexibility of the transport coefficients. On the other hand, it has been shown that the LBM scheme can also be considered as a special discretized form of the continuous Boltzmann equation. Through a Chapman-Enskog analysis, one can recover the governing continuity and Navier-Stokes equations from the LBM algorithm. In addition, the pressure field is also directly available from the density distributions and hence there is no extra Poisson equation to be solved as in traditional CFD methods.

Lattices and the DnQm classification

Lattice Boltzmann models can be operated on a number of different lattices, both cubic and triangular, and with or without rest particles in the discrete distribution function.

A popular way of classifying the different methods by lattice is the DnQm scheme. Here "Dn" stands for "n dimensions" while "Qm" stands for "m speeds". For example, D3Q15 is a three-dimensional Lattice Boltzmann model on a cubic grid, with rest particles present. Each node has a crystal shape, and can deliver particles to each of the six neighboring nodes which share a surface, the eight neighboring nodes sharing a corner, and itself. [Succi, p 68] (The D3Q15 model does not contain particles moving to the twelve neighboring nodes which share an edge; adding those would create a "D3Q27" model.)

imulation of mixtures

Simulating multiphase/multicomponent flows has always been a challenge to conventional CFD because of the moving and deformable interfaces. More fundamentally, the interfaces between different phases (liquid and vapor) or components (e.g, oil and water) originate from the specific interactions among fluid molecules. Therefore it is difficult to implement such microscopic interactions into the macroscopic Navier–Stokes equation. However, in LBM, the particulate kinetics provides a relatively easy and consistent way to incorporate the underlying microscopic interactions by modifying the collision operator. Several LBM multiphase/multicomponent models have been developed. Here phase separations are generated automatically from the particle dynamics and no special treatment is needed to manipulate the interfaces as in traditional CFD methods. Successful applications of multiphase/multicomponent LBM models can be found in various complex fluid systems, including interface instability, bubble/droplet dynamics, wetting on solid surfaces, interfacial slip, and droplet electrohydrodynamic deformations.

Limitations

Despite the increasing popularities of LBM in simulating complex fluid systems, this novel approach has some limitations. At present, high-Mach number flows in aerodynamics are still difficult for LBM, and a consistent thermo-hydrodynamic scheme is absent. However, as with Navier–Stokes based CFD, LBM methods have been successfully coupled to thermal-specific solutions to enable heat transfer (solids-based conduction, convection and radiation) simulation capability. For multiphase/multicomponent models, the interface thickness is usually large and the density ratio across the interface is small when compared with real fluids. Nevertheless, the wide applications and fast advancements of this method during the past twenty years have proven its potential in computational physics, including microfluidics: LBM demonstrates promising results in the area of high Knudsen number (defined by the ratio of the mean free path between molecules to a geometric length scale) flows.

Mathematical details

The Boltzmann equation is an evolution equation for a single particle probability distribution function f(x,v,t):

:partial_t f + vpartial_x f + Fpartial_v f=Omega

where F is an external force and Omega is a collision integral. The lattice Boltzmann method discretizes this equation by limiting space to a lattice and the velocity space to a discrete set of velocities v_i. The discretized Boltmann equation, which is the Lattice Boltzmann equation then reads:

:f_i(x+v_i,t+1)-f(x,t) + F_i=Omega

The collision operator is often approximated by a BGK collision operator:

:Omega = frac{1}{ au} (f_i^0-f_i)

where f_i^0 is the local equilibrium distribution.

The moments of the f_i give the local conserved quantities. The density is given by

: ho=sum_i f_i

and the local momentum is given by

: ho u = sum_i f_i v_i.

For the popular isothermal lattice Boltmann methods these are the only conserved quantities. Thermal models also conserve energy and therefore have an additional conserved quantity:

: ho heta + ho u u =sum_i f_i v_i v_i.

The collision operator has to respect the conservation laws. Therefore the equilibirum distribution f_i^0 must have the same conserved moments as the f_i.

oftware

* [http://www.exa.com/pages/pflow/pflow_physics.html PowerFLOW] : commercial CFD code which uses LBM
* [http://elbeem.sourceforge.net/ El'Beem] : free CFD code (GPL) which uses LBM
* [http://thomas-pohl.info/work/lba.html J-Lattice-Boltzmann] : interactive Java applet for experimenting with LBM
* [http://physics.ndsu.edu/people/faculty/wagner/lattice_boltzmann_codes/ C examples] : Some simple example LBM code in C.
* [http://www.lbmethod.org/openlb/ OpenLB] : Open source lattice Boltzmann code

External links

* [http://www.lbmethod.org/ LBMethod.org] : A site with various resources related to LBM, including a forum.
* [http://www.science.uva.nl/research/scs/projects/lbm_web/lbm.html LBM Method]
* [http://physics.ndsu.edu/fileadmin/physics.ndsu.edu/Wagner/LBbook.pdf Lattice Boltzmann summary]
* [http://www.lstm.uni-erlangen.de/lbm/lb-methods/ Lattice Boltzmann mailing list]
* [http://www.icmmes.org/ Website of the annual conference on Lattice Boltzmann methods and their applications ]
* [http://dsfd.lmpt.ufsc.br/ DSFD series of international conferences on the Discrete Simulation of Fluid Dynamics ]

Further reading

*
*
*

Notes


Wikimedia Foundation. 2010.

Игры ⚽ Нужно решить контрольную?

Look at other dictionaries:

  • Lattice gas automaton — Lattice gas automata (LGA) or lattice gas cellular automata (LGCA) methods are a series of cellular automata methods used to simulate fluid flows. It was the precursor to the lattice Boltzmann methods. From the LGCA, it is possible to derive the… …   Wikipedia

  • Boltzmann equation — For other uses, see Boltzmann s entropy formula, Stefan–Boltzmann law and Maxwell–Boltzmann distribution The Boltzmann equation, also often known as the Boltzmann transport equation, devised by Ludwig Boltzmann, describes the statistical… …   Wikipedia

  • Ludwig Boltzmann — Infobox Scientist name = Ludwig Boltzmann |225px image width = 225px caption = Ludwig Eduard Boltzmann (1844 1906) birth date = birth date|1844|2|20|mf=y birth place= Vienna, Austrian Empire death date = death date and age|1906|9|5|1844|2|20|mf=y …   Wikipedia

  • Monte Carlo methods for electron transport — The Monte Carlo method for electron transport is a semiclassical Monte Carlo(MC) approach of modeling semiconductor transport. Assuming the carrier motion consists of free flights interrupted by scattering mechanisms, a computer is utilized to… …   Wikipedia

  • Prabhu Lal Bhatnagar — Born 18 August 1912(1912 08 18) Kota, India Died 5 October 1976(1976 10 05) …   Wikipedia

  • List of numerical analysis topics — This is a list of numerical analysis topics, by Wikipedia page. Contents 1 General 2 Error 3 Elementary and special functions 4 Numerical linear algebra …   Wikipedia

  • Computational fluid dynamics — Computational physics Numerical analysis  …   Wikipedia

  • Метод решёточных уравнений Больцмана — Методы решёточных уравнений Больцмана (англ. Lattice Boltzmann methods, LBM) класс методов вычислительной гидродинамики для моделирования жидкостей. В отличие от многих других методов, метод LBM не решает уравнения Навье Стокса, а моделирует …   Википедия

  • Parallel computing — Programming paradigms Agent oriented Automata based Component based Flow based Pipelined Concatenative Concurrent computing …   Wikipedia

  • GPGPU — General purpose computing on graphics processing units (GPGPU, also referred to as GPGP and to a lesser extent GP²) is the technique of using a GPU, which typically handles computation only for computer graphics, to perform computation in… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”