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
The problem about temperature jump is known from the end of XIX century . M.Smoluhovsky has constructed in  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  and ).
In 1972 in work  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  -  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  where the analytical solution of the Smoluchowsky problem has been received.
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.
In the present work the generalized source-method from  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  has already been applied in problems of electrodynamics of plasma  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  )
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
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