# Solution of Boltzmann Equation

### From CFD-Wiki

(New page: == Definition == Consider the classical Boltzmann equation for a simple, dilute gas of particles <math>f_t + (v,{\rm grad}_x f) = Q(f,f)</math> which describes the time evolution of th...) |
|||

(3 intermediate revisions not shown) | |||

Line 44: | Line 44: | ||

One of the first discrete versions of the Boltzmann equation was published by D. Goldstein, B. Sturtevant and | One of the first discrete versions of the Boltzmann equation was published by D. Goldstein, B. Sturtevant and | ||

- | J.E. Broadwell. Many authors then published different ideas to lead to a discrete version of the Boltzmann collision operator Rogier, Schneider in 1994; Ohwada in 1993, Wagner in 1995, Platkowski, Illner in 1988, Palczewski, Schneider in 1998, Panferov in 1997, Panferov, Heintz in 2002 | + | J.E. Broadwell. Many authors then published different ideas to lead to a discrete version of the Boltzmann collision operator Rogier, Schneider in 1994; Ohwada in 1993, Wagner in 1995, Platkowski, Illner in 1988, Palczewski, Schneider in 1998, Panferov in 1997, Panferov, Heintz in 2002. |

- | + | ||

- | + | ||

The main difficulty with the deterministic approximation of the Boltzmann collision integral besides its | The main difficulty with the deterministic approximation of the Boltzmann collision integral besides its | ||

- | high dimensionality is the fact that any grid for the integration over the whole space <math>{\rm R}^3</math> will not fit for the integration over the unit sphere <math>S^2</math>. Thus only <math>O(n)</math> irregularly distributed integration points belong to the unit sphere if <math>n</math> regular points in one direction are used for the approximation of the <math>{\rm R}^3</math> integral. | + | high dimensionality is the fact that any grid for the integration over the whole space <math>{\rm R}^3</math> will not fit for the integration over the unit sphere <math>S^2</math>. Thus only <math>O(n)</math> irregularly distributed integration points belong to the unit sphere if <math>n</math> regular points in one direction are used for the approximation of the <math>{\rm R}^3</math> integral. |

A. Bobylev, A. Palczewski and J. Schneider in 1995 considered this direct approximation of the Boltzmann collision integral and showed that the corresponding numerical method is consistent. The arithmetical work is <math>O(n^7)</math> per time step and the formal accuracy is <math>O(n^{-1/2})</math>. | A. Bobylev, A. Palczewski and J. Schneider in 1995 considered this direct approximation of the Boltzmann collision integral and showed that the corresponding numerical method is consistent. The arithmetical work is <math>O(n^7)</math> per time step and the formal accuracy is <math>O(n^{-1/2})</math>. | ||

- | + | L. Pareschi and B. Perthame in 1996 considered deterministic spectral methods for the Boltzmann equation and in 2000 [http://www.jstor.org/stable/2587302 L. Pareschi and G. Russo] show that these method are consistent and are able to provide spectrally accurate solutions at <math>O(n^6)</math> computational cost. Moreover the method are able to treat also the difficult case of non cut-off collision kernel as shown by L.Pareschi, G.Toscani and C.Villani in 2003. | |

- | Recent results that | + | Later, Ibragimov and Rjasanow in 2002, F.Filbet, L. Pareschi and G.Toscani in 2004 studied spectral methods for numerical solution of Boltzmann equations and Ibragimov, Rjasanow in 2002 proved an <math>O(n^{-2})</math> numerical accuracy with <math>O(n^6)</math> arithmetic complexity. |

+ | |||

+ | Recent results that improve the numerical solution of Boltzmann equations are done in 2007 by [http://www.elegant-mathematics.com/gpu-boltzmann.pdf Ibragimov and Ibragimova] who suggests to use a multilinear approximation to approximate a velocity subspace. Hence, the function state <math>f</math> is approximated as one of the following methods: | ||

<math>f(t,x,v) = \sum_{i_1, i_2, i_3}^{r_1,r_2,r_3} f_{i_1}^{(1)}(t,x,v_1) f_{i_2}^{(2)}(t,x,v_2) f_{i_3}^{(3)}(t,x,v_3) g_{i_1, i_2, i_3}, | <math>f(t,x,v) = \sum_{i_1, i_2, i_3}^{r_1,r_2,r_3} f_{i_1}^{(1)}(t,x,v_1) f_{i_2}^{(2)}(t,x,v_2) f_{i_3}^{(3)}(t,x,v_3) g_{i_1, i_2, i_3}, | ||

Line 68: | Line 68: | ||

Several real time industrial applications solved with the help of this method are mentioned on the web site of [http://www.elegant-mathematics.com/ Elegant Mathematics] who develops such solvers and provide test calculations for free. | Several real time industrial applications solved with the help of this method are mentioned on the web site of [http://www.elegant-mathematics.com/ Elegant Mathematics] who develops such solvers and provide test calculations for free. | ||

+ | |||

+ | Finally let us mention that an algorithm which strongly reduces the computational cost of spectral methods to <math>O(n^3logn)</math> in the 3D case for hard spheres without loosing the spectral accuracy has been obtained recently by F.Filbet, C.Mouhot and L.Pareschi in 2006 and several applications of this method has been presented to show its potential in computing solution of CFD problems involving the Boltzmann equation. |

## Latest revision as of 03:08, 13 May 2010

## Definition

Consider the classical Boltzmann equation for a simple, dilute gas of particles

which describes the time evolution of the particle density

Here denotes the set of non-negative real numbers and is a domain in physical space. The right-hand side of the Boltzmann equation, known as the collision integral or the collision term, is of the form

where

are the pre-collision velocities,

is a unit vector,

are the post-collision velocities and

is the collision kernel. The operator represents the change of the distribution function due to the binary collisions between particles. A single collision results in a change of the velocities of the colliding partners with

where denotes the relative velocity. The Boltzmann equation is subjected to an initial condition

and to the boundary conditions on .

The kernel can be written as

The function is the differential cross-section and is the scattering angle.

One of the first discrete versions of the Boltzmann equation was published by D. Goldstein, B. Sturtevant and J.E. Broadwell. Many authors then published different ideas to lead to a discrete version of the Boltzmann collision operator Rogier, Schneider in 1994; Ohwada in 1993, Wagner in 1995, Platkowski, Illner in 1988, Palczewski, Schneider in 1998, Panferov in 1997, Panferov, Heintz in 2002.

The main difficulty with the deterministic approximation of the Boltzmann collision integral besides its high dimensionality is the fact that any grid for the integration over the whole space will not fit for the integration over the unit sphere . Thus only irregularly distributed integration points belong to the unit sphere if regular points in one direction are used for the approximation of the integral.

A. Bobylev, A. Palczewski and J. Schneider in 1995 considered this direct approximation of the Boltzmann collision integral and showed that the corresponding numerical method is consistent. The arithmetical work is per time step and the formal accuracy is .

L. Pareschi and B. Perthame in 1996 considered deterministic spectral methods for the Boltzmann equation and in 2000 L. Pareschi and G. Russo show that these method are consistent and are able to provide spectrally accurate solutions at computational cost. Moreover the method are able to treat also the difficult case of non cut-off collision kernel as shown by L.Pareschi, G.Toscani and C.Villani in 2003.

Later, Ibragimov and Rjasanow in 2002, F.Filbet, L. Pareschi and G.Toscani in 2004 studied spectral methods for numerical solution of Boltzmann equations and Ibragimov, Rjasanow in 2002 proved an numerical accuracy with arithmetic complexity.

Recent results that improve the numerical solution of Boltzmann equations are done in 2007 by Ibragimov and Ibragimova who suggests to use a multilinear approximation to approximate a velocity subspace. Hence, the function state is approximated as one of the following methods:

where are almost constant.

Indeed, when are equal to 1, we can easily get Navies-Stocks equation, when ranks are equal to 2, we have very robust model for the boundary layer, with small enough rank (3-5) we can definitely explain a shock waves, etc.

Several real time industrial applications solved with the help of this method are mentioned on the web site of Elegant Mathematics who develops such solvers and provide test calculations for free.

Finally let us mention that an algorithm which strongly reduces the computational cost of spectral methods to in the 3D case for hard spheres without loosing the spectral accuracy has been obtained recently by F.Filbet, C.Mouhot and L.Pareschi in 2006 and several applications of this method has been presented to show its potential in computing solution of CFD problems involving the Boltzmann equation.