A FORMAC program is described which calculates the derivatives of appropriate generating functions to obtain the matrix elements with respect to the Burnett functions of the linearized Boltzmann cross-collision operators of a gaseous mixture. The general procedure discussed, however, is applicable t
Computer programs for the Boltzmann collision matrix elements
β Scribed by P. Das
- Publisher
- Elsevier Science
- Year
- 1989
- Tongue
- English
- Weight
- 589 KB
- Volume
- 55
- Category
- Article
- ISSN
- 0010-4655
No coin nor oath required. For personal study only.
β¦ Synopsis
When the distribution function in the kinetic theory of gases is expanded in a basis of orthogonal functions, the Boltzmann collision operators can be evaluated in terms of appropriate matrix elements. These matrix elements are usually given in terms of highly complex algebraic expressions. When Burnett functions, which Consist of Sonine polynomials and spherical harmonics, are used as the basis, the irreducible tensor formalism provides expressions for the matrix elements that are algebraically simple, possess high symmetry, and are computationally more economical than in any other basis.
The package reported here consists of routines to compute such matrix elements in a Burnett function basis for a mixture of hard sphere gases, as also the loss integral of a Burnett mode and the functions themselves. The matrix elements involve the Clebsch-Gordan and Brody-Moshinsky coefficients, both of which are used here for unusually high values of their arguments. For the purpose of validation both coefficients are computed using two different methods.
Though written for hard sphere molecules the package can, with only slight modification, be adapted to more general molecular models as well.
π SIMILAR VOLUMES
A. Salin / Wave functions and collision matrix elements for one-electron diatomic molecules Overlay structure: Optionally overlaid Nature of the physical problem: Determination of the norm of the wave-functions and of radial and rotational coupling No. of magnetic tapes required: none matrix element
## Abstract Steep high angular momentum Gaussian basis functions in the vicinity of a nucleus whose inner electrons are replaced by an effective core potential may lead to numerical instabilities when calculating matrix elements of the core potential. Numerical roundoff errors may be amplified to a