Solving the kinetic equation for all Knudsen numbers

被引:51
作者
Beylich, AE [1 ]
机构
[1] Rhein Westfal TH Aachen, D-52062 Aachen, Germany
关键词
D O I
10.1063/1.870322
中图分类号
O3 [力学];
学科分类号
08 ; 0801 ;
摘要
Starting from the Boltzmann kinetic equation, a system is constructed which allows one to calculate flow fields with strong deviation from equilibrium and for Knudsen numbers 0 < Kn(0)<infinity. This is achieved by interlacing the path-integral form of the kinetic equation with the equations of conservation (the first five moments of the equation of transfer). Thus, at the lower level, we maintain the ability of the drifting part of the kinetic equation to drive the state toward nonequilibrium and, at the upper level, we keep the properties of the hydrodynamic (collisional) system to transport information on the grid system in physical (r) space in a nonmaterial way by waves. The distribution function f(r,c,t) is defined on an equally spaced Cartesian grid in r space and in velocity (c) space. An important step in the construction of the interlaced system is to interpolate the distribution function f(r,c,t) at interstitial points in r,c space. Approximations of the path-integral equation and different levels of approximation of the gain function (collision term) allow an adequate and efficient description of flows at any degree of deviation from equilibrium. Introducing a phase function in c space allows one to select those collisions that contribute to the gain function. The determination of collision terms then becomes very efficient, and calculating a time step is reduced by orders of magnitude. The problem of maintaining conservation of the lower moments, generally arising when the gain collision integral is to be determined at interstitial points, is solved by imposing on f a set of integral constraints. The system was developed and tested using the Riemann problem and steady planar Couette flow. Results demonstrate that the system is able to maintain conservation, describe nonequilibrium over a large range of Knudsen numbers and over orders of magnitude of density change. In the example of Couette flow, the system was tested for 0.0001 < Kn(0)< 50 and wall Mach numbers Ma(w)less than or equal to 3. (C) 2000 American Institute of Physics. [S1070-6631(00)00802-3].
引用
收藏
页码:444 / 465
页数:22
相关论文
共 32 条
[1]  
ABE T, 1997, 553 ISAS U TOK, V42, P221
[2]   DENSITY PROFILES IN ARGON AND NITROGEN SHOCK-WAVES MEASURED BY ABSORPTION OF AN ELECTRON-BEAM [J].
ALSMEYER, H .
JOURNAL OF FLUID MECHANICS, 1976, 74 (APR6) :497-513
[3]  
[Anonymous], 1991, ADV KINETIC THEORY C
[4]   A MODEL FOR COLLISION PROCESSES IN GASES .1. SMALL AMPLITUDE PROCESSES IN CHARGED AND NEUTRAL ONE-COMPONENT SYSTEMS [J].
BHATNAGAR, PL ;
GROSS, EP ;
KROOK, M .
PHYSICAL REVIEW, 1954, 94 (03) :511-525
[5]  
Bird G. A., 1994, MOL GAS DYNAMICS DIR
[6]  
Bird G.A., 1976, MOL GAS DYNAMICS
[7]   A GENERAL KINETIC THEORY OF LIQUIDS .1. THE MOLECULAR DISTRIBUTION FUNCTIONS [J].
BORN, M ;
GREEN, HS .
PROCEEDINGS OF THE ROYAL SOCIETY OF LONDON SERIES A-MATHEMATICAL AND PHYSICAL SCIENCES, 1946, 188 (1012) :10-18
[8]   SHOCK STRUCTURE IN A SIMPLE DISCRETE VELOCITY GAS [J].
BROADWELL, JE .
PHYSICS OF FLUIDS, 1964, 7 (08) :1243-1247
[9]  
Broadwell JE, 1964, J FLUID MECH, V19, P401, DOI 10.1017/S0022112064000817
[10]  
CERCIGNANI C, 1988, APPL MATH SERIES, V67