In the present work the classical problem of the kinetic theory of gases (the Smoluchowsky’ problem about temperature jump in rarefied gas) is considered. The rarefied gas fills half-space over a flat firm surface. logarithmic gradient of temperature is set far from surface. The kinetic equation with modelling integral of collisions in the form of BGK-model (Bhatnagar, Gross and Krook) is used.

The general mirror-diffuse boundary conditions of molecules reflexions of gas from a wall on border of half-space (Maxwell conditions) are considered. Expanding distribution function on two orthogonal directions in space of velocities, the Smoluchowsky’ problem to the solution of the homogeneous vector one-dimensional and one-velocity kinetic equation with a matrix kernel is reduced.

Then generalization of source-method is used and boundary conditions include in non-homogeneous vector kinetic equation. The solution in the form of Fourier integral is searched. The problem is reduced to the solution of vector Fredholm integral equation of the second sort with matrix kernel.

The solution of Fredholm equation in the form of Neumann’s polynoms with vector coefficients is searched. The system vector algebraic interengaged equations turns out. The solution of this system is under construction in the form of Neumann’s polynoms. Comparison with well-known Barichello—Siewert’ high-exact results is made. Zero and the first approach of jumps of temperature and numerical density are received. It is shown, that transition from the zero to the first approach raises 10 times accuracy in calculation coefficients of temperature and concentration jump.

Key words: the Smoluchowski’ problem, collisional gas, temperature and concentration jump, vector Fredholm equation of second sort.

PACS numbers: 05.20.Dd Kinetic theory, 47.45.-n Rarefied gas dynamics, 02.30.Rz Integral equations, 51. Physics of gases, 51.10.+y Kinetic and transport theory of gases.

A new method for solving of vector problems for kinetic equations with Maxwell boundary conditions

A. V. Latyshev111

Faculty of Physics and Mathematics,

Moscow State Regional University, 105005,

Moscow, Radio str., 10–A

1 Introduction

The problem about temperature jump is known from the end of XIX century [1]. M.Smoluhovsky has constructed in [1] the theory of temperature jump in the rarefied gas. Since then this problem the invariable attention already draws for a long time to itself (history of this question see in [2] and [3]).

This problem has been solved with use of approximate and numerical methods as for the modelling equations, and for full (non-linear) Boltzmann equation [4] - [18].

In 1972 in work [19] the problem about temperature jump with use of the modelling Boltzmann equation with collisional integral BGK (Bhatnagar, Gross, Krook) and with frequency of collisions of the molecules, proportional to the module of velocities of molecules has analytically been solved.

Attempts of the exact solution of this problem about temperature jump with diffusion boundary conditions and with use of the modelling Boltzmann equation with collisional integral BGK with constant frequency of collisions of molecules [20] - [23] begin with this moment.

The solution of this problem encounters considerable difficulties. This problem is formulated in the form of a vector boundary problem. The solution of last problem meets the solution of a vector boundary value Riemann—Hilbert problem with the matrix coefficient, having points of branchings. These difficulties have been overcome only in 1990 in our work [24] where the analytical solution of the Smoluchowsky problem has been received.

By the method developed in [24], further problems about temperature jump in molecular gases [25], and also in metal [26] have been analytically solved.

Along with the Smoluchowsky’ problem the big interest represents studying of behaviour of gas at weak evaporation (condensation) from a surface. These problems are called as the generalized Smoluchowsky’ problem in view of that boundary conditions in these problems differ slightly.

Let us notice, that in works [17, 18] the various kinetic models were used, in particular, the model of Shakhov (or, S-model, see [27]).

For the solution of boundary half-space problems with accommodation some methods [28] - [30] have been developed, allowing to receive the solution of this problem with any degree of accuracy.

In the present work the generalized source-method from [30] extends on a vector case to which the problem about temperature jump is reduced. Thus the effective method of the solution of boundary problems with mirror - diffusion boundary conditions (Maxwell conditions) develops. We will notice, that the method from [30] has already been applied in problems of electrodynamics of plasma [31] and in condensate problems of Bose—Einstein.

At the heart of an offered method the idea lays to include the boundary condition in the form of a source in the kinetic equation.

The method basis consists in the following. At first in half-space are formulated a vector problem about the temperature jump with boundary Maxwell conditions. Then unknown function continuations in conjugated half-space in the even method on spatial and on velocity variables. In half-space also are formulated the problem about temperature jump.

Now let us expand unknown function (which we will name also distribution function) on two composed: Chapman—Enskog’ distribution function and the second part of function distributions , correspoding to continuous spectrum (see [30] )


Owing to that Chapman—Enskog’ distribution function there is a linear combination of discrete solutions of the initial equation, function also is the solution of the kinetic equations. Function vanishes in zero far from wall. On a wall this function satisfies to boundary Maxwell condition.

Further we will transform the equation for function . We include in this equation boundary condition on wall for function in the form of a member of source-type laying in a plane .

We will underline, that function satisfies to the received equation in both conjugated half-spaces and .

We solve this equation in the second and the fourth quarters of a phase plane as the linear differential equation of the first order, considering known the right part of the equation . From the received solutions we found the boundary values of unknown function at , entering into the equation.

Now we expand by Fourier integrals unknown function , an unknown right part and Dirac delta-function. Boundary values of the unknown functions are thus expressed by the same integral on the spectral density functions .

Substitution of Fourier integrals in the kinetic equation and expression for the right part leads to the vector characteristic system of equations. If to exclude from this system the spectral density of function , we will receive vector Fredholm integral equation of the second sort.

Believing the gradient of the logarithm of temperature is setting, we will expand the unknown quantities of temperature and concentration jumps and also spectral density by polynoms on degrees of coefficient of diffusion (these are Neumann’s polynoms). On this way we receive system of the hooked equations on coefficients of polynoms for spectral density. Thus all equations on coefficients of spectral density have singularity (a pole of the second order in zero). Excepting these singularities consistently, we will construct all members of the polynoms for quantities of temperature and concentration jumps and for spectral density .

2 Statement problem

Let the rarefied one-nuclear gas occupies half-space over the flat firm surface laying in a plane . Far from a wall the logarithmic gradient of temperature is set

We take the stationary kinetic equation of relaxation type with collisional integral BGK (Bhatnagar, Gross and Krook)

where is the time between two consecutive collisions of molecules, is the collisional frequence of gaseous molecules, is the equilibrium distribution function,

where is the mass of molecule, is the Boltzmann constant, is the gas temperature,

is the gas number density (concentration),

Further we will be linearize the kinetic equation and search distribution function in the form

where is the absolute Maxwellian,

where are number density and gas temperature in some point, for example, in origin of coordinates.

Let us be linerize distribution of numerical density and temperature concerning parametres and

According to (2.1) for distribution of numerical density we have


For distribution of temperature we receive

We notice that

Now for temperature we receive

We notice that

Hence, for relative change of temperature it is had

We will be linearize equilibrium function of distribution

We receive the following equation after linearizing the kinetic BGK–equation according to (2.1)

Let us enter dimensionless velocities and parametres — dimensionless velocity of molecules , where , dimensionless time , dimensionless coordinate

where is the mean free path of gaseous molecules, is the thermal velocity of the molecules movements, having an order of velocity of a sound.

Now the kinetic equation will be transformed to the form


Further a variable we will designate again through .

Let us transform the linear kinetic equation to the form

with kernel

It is easy to check up, that the equation (2.2) has the following partial solutions

Let us construct asymptotic Chapman—Enskog distribution in the form of the linear combination of partial solutions of the equation (2.2) with arbitrary constants

where are arbitrary constants.

For finding of these constants we will take advantage of definition of the macroscopical parametres. From definition of numerical density (concentration)

follows, that the extrapolated concentration of gas on the wall is equal

From here we have


Hence, the quantity of jump of concentration is searched under the formula

Substituting expression (2.3) in this equality, we have

Setting gradient of temperature far from a wall means, that temperature distribution in half-space looks like

This distribution we will present in the form



is the required quantity of temperature jump.

From expression (2.4) it is visible, that relative change of temperature far from walls it is described by linear function

Relative change of temperature we will present in the form

Far from a wall relative change of temperature transforms as follows

Substituting in this equality expression (2.3) for , we find, that

Comparing expressions (2.5) and (2.6), we find, that

Thus, asymptotic part of function of distribution (at ) it is constructed and on the basis stated above transforms in the form

Let us formulate down boundary conditions to the equation (2.2). At first let us formulate mirror–diffusion boundary condition on a wall for full function of distribution

Here is the accommodation coefficient, i.e. a part of the molecules flying after reflexion from a wall with Maxwell distribution on velocities, is the part of the molecules reflected from a wall purely mirror.

Using (2.1), from here we receive a boundary condition of a problem onto wall

Let us demand, that far from a wall distribution function passed into Chapman—Enskog distribution with coordinate growth

From here according to (2.1) for function we receive the following boundary conditions

Here is the asymptotic Chapman—Enskog distribution, entered above.

So, the boundary problem about finding of jumps of temperature and concentration of gas (vapor) over a flat surface consists in finding of the such solution of the equation (2.2), which satisfies to boundary conditions (2.7) and (2.8).

3 Reduction to vector boundary problem

If to use substitution

that equation (2.2) is reduced breaks up to two equations


This equation we will present in the vector form

Here is the vector-column

and matrix kernel have the following form

The right part of the equation (3.2)

has clear physical sense. The vector-column looks like

i.e., components of this vector consist of the relative changes of numerical density of gas and relative change temperatures (concerning equilibrium values ). It is possible to present this vector in the form


According to (3.1) from (2.7) and (2.8) for the vector-functions we receive the following vector boundary conditions



Function is the solution of the equation (3.2). Hence, if to search for the solution of the equation (3.2) in the form

then function is still searched from the equation (3.2)

and boundary conditions (3.3) and (3.4) will be transformed thanking (3.5) to the following form


Let us solve further the problem consisting of the solution of the equation (3.6) with boundary conditions (3.7) – (3.10).

4 Kinetic equation with source

For solution of this problem the auxiliary problem is required to us in "negative"  half-space . That it to formulate, we will continue function as follows

Let us notice, that at continuation (4.1) logarithmic temperature gradient , which for "positive"  half-spaces we will designate through , changes the sign

Besides, we will notice, that function automatically satisfies to equality (4.1): . This equality means, that the equality (3.1) is carried out for function : .

Hence, boundary conditions in "negative"  space are formulated as follows


Let us unite both problems — in "positive"  and "negative"  half-spaces — in one, having included boundary conditions in the kinetic equation by means of member of type of the source


These function are finding from equalities

5 Vector Fredholm equation of second sort

The solution of the equations (4.2) and (4.3) we search in the form of Fourier integrals

From equalities (4.3) and (5.1) follows, that

Two following equalities follow from equalities (4.4)

From the kinetic equation (4.2) by means of (5.4) it is found

Substituting (5.5) in (5.3), we come to the vector integral Fredholm equation of the second sort

Here is the dispersion matrix-function

where is the unit matrix of the second order,

the matrix kernel of integral Fredholm equation is defined by integral expression

It is obvious, that

6 Solution of vector Fredholm equation

For the solution of the equation (5.6) we will search in the form of Neumann’s polynoms

Let us substitute (6.1) and (6.2) in the equation (5.6). We receive system of the hooked equations

Let us calculate in an explicit form the matrixes entering into the equation (5.6)

Here integrals are entered

Let us notice, that the dispersion matrix is proportional to . Let us notice, that

Therefore the dispersion matrix–function is equal

Let us write out matrix elements


Further we find a matrix