Effect of the mode of gas-surface interaction on intensive monatomic gas evaporation

The direct Monte Carlo simulation method is used for investigating the effect of the thermal accommodation coefficient α _E on the relations on the Knudsen layer boundary in the presence of intensive evaporation. The model of mirror reflection of molecules from the surface is considered. It is shown that diffuse reflection with α _E = 0 leads to almost the same relations on the Knudsen layer boundary as mirror reflection. The accuracy of the moment method is estimated in application to the problems of intensive evaporation with diffuse and mirror reflection from the surface.     

Solution of the linearized couette-flow problem in a rarefied gas by the integral-diffusion method

In ordinary diffusion theory the transfer of properties is determined by the local gradients of the corresponding fields. As the mean free path increases, the flux density becomes an integral quantity and is determined by a neighborhood of the point under consideration of the order of a few mean free paths. In a previous article [1], the author proposed a model for a one-dimensional transfer process in linear rarefield-gas problems based on the analogy with radiative transfer. The same approach, though without directional averaging, is used in the present paper to analyze the linearized Couette flow problem. The solution obtained here has the properties of the solution obtained by more exact methods based on the solution of the Boltzmann equation [3-4].Nomenclature p_xy shear stress - c mean thermal velocity of molecules - 2/3 A mean free path - d half-width of channel - ±w₀ plate velocity - c nonequilibriumvalue of momentum flux density - y transverse coordinate - ratio of specific heats - W dimensionless velocity - P_xy shear stress scaled with respect to the shear stress in free-molecule flow - Y dimensionless coordinate - W₁(y) velocity distribution according to Millikan's solution - coefficient of viscosity - R Reynolds number - K Knudsen number     

Simulations for gas flows in microgeometries using the direct simulation Monte Carlo method

Micro gas flows are often encountered in MEMS devices and classical CFD could not accurately predict the flow and thermal behavior due to the high Knudsen number. Therefore, the gas flow in microgeometries was investigated using the direct simulation Monte Carlo (DSMC) method. New treatments for boundary conditions are verified by simulations of micro-Poiseuille flow, compared with the previous boundary treatments, and slip analytical solutions of the continuum theory. The orifice flow and the corner flow in microchannels are simulated using the modified DSMC codes. The predictions were compared with existing experimental phenomena as well as predictions using continuum theory. The results showed that the channel geometry significantly affects the microgas flow. In the orifice flow, the flow separation occurred at very small Reynolds numbers. In the corner flow, no flow separation occurred even with a high driving pressure. The DSMC results agreed well with existing experimental information.     

Supersonic condensation of a monatomic gas

The problem of steady supersonic condensation of a monatomic gas on a plane evaporating surface is solved in the Knudsen layer by the direct statistical modeling method. The domain of existence of the solution of the problem is determined. The results of calculating the structure of the Knudsen layer near the surface are presented. A topological picture of the solutions of the strong evaporation and subsonic and supersonic strong condensation problems is given as a function of the Mach number, determined from the normal velocity component, and the other governing parameters.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 171–175, May–June, 1990.     

Analytical solution of the strong evaporation (condensation) problem

An exact solution of the model Boltzmann equation with BGK (Bhatnagar-Gross-Krook) collision operator is obtained in the problem of strong evaporation (condensation) from a plane surface. The generalized eigenvectors of the corresponding characteristic equation are found. The existence and uniqueness of the expansion of the solution in eigenvectors of the continuous and discrete spectra are demonstrated. This expansion reduces to a vector Riemann-Hilbert boundary-value problem with matrix coefficient. An apparatus for the diagonalization and factorization of the boundary-value problem coefficient is developed. The matrix diagonalizing the problem coefficient has branch points in the complex plane which depend parametrically on the evaporation (condensation) rateu. The solution of the problem is investigated in terms ofu and the physical characteristics of the evaporation process are described.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.6, pp. 143–155, November–December, 1993.     

Use of the Monte Carlo method to calculate the flow of a highly rarefied gas in systems with arbitrary wall configuration

A method is described for modeling collisions of gas molecules with the walls of a system whose geometry is given numerically in the computer memory from a graphical representation of the walls (for example, from their working drawings). When changing from the calculation of one system to another only the information on the wall changes; the computational program remains the same. The method is applicable to problems of rarefied gasdynamics which are solvable by the Monte Carlo method; in the following it is used to calculate the conductance of elements of high vacuum lines and the compression ratio created by molecular vacuum pumps, and also to calculate the forces acting on a body which rotates in a cavity filled with a highly rarefied gas.In conclusion the authors wish to thank Yu. I. Neimark for formulating the problem and discussions of the results obtained.     

Intensive subsonic condensation of a monatomic gas

The best developed method for solving the complete Boltzmann equation without the a priori assumptions concerning the form of the distribution functions inherent in moment methods is direct statistical modeling. Solutions of the Boltzmann equation for intensive subsonic condensation are obtained by this method over a fairly broad interval of the parameters. The accuracy and region of applicability of the moment methods are estimated on the basis of the exact solutions obtained.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 165–169, January–February, 1989.     

Monte Carlo simulation of the heat flow in a dense hard sphere gas

The one-dimensional steady heat flow in a dense hard sphere gas is studied solving the Enskog equation numerically by a recently proposed DSMC-like particle scheme. The accuracy of the solutions is assessed through a comparison with solutions obtained from a semi-regular method which combines finite difference discretization with Monte Carlo quadrature techniques. It is shown that excellent agreement is found between the two numerical methods. The solutions obtained from the Enskog equation have also been found in good agreement with the results of molecular dynamics simulations.     

Molecular models for simulation of rarefied gas flows using direct simulation Monte Carlo method

The direct simulation Monte Carlo (DSMC) method is a technique for the numerical simulation of the rarefied gas flows by employing simulated molecules in simulated physical spaces. In the procedures involved in DSMC, the accuracy of the simulation of intermolecular collisions depends on the collision model adopted in the collision routine. The simplest molecular model is the hard-sphere model. In order to improve the accuracy of the simulations, more and more refined collision models were introduced for the use in DSMC. Thus, the variable hard-sphere, the variable soft-sphere, the generalised hard-sphere, the generalised soft-sphere and the variable sphere models were put forward by various researchers. And, all these models have met with varying degrees of success. Meanwhile, the Borgnakke-Larsen model, statistical inelastic cross-section models for both continuous and discrete internal energy and the dynamic molecular collision model were proposed for the treatment of polyatomic molecules in which transfer of energy among translational, rotational and vibrational degrees of freedom is possible. This paper gives a brief introduction to the intermolecular potentials based on which the molecular models have been constructed. Then the various models are introduced in the chronological sequence; finally concluding with a brief summary of the progress that has been made so far in this area.     

Monte Carlo method of radiative transfer applied to a turbulent flame modeling with LES

Radiative transfer plays an important role in the numerical simulation of turbulent combustion. However, for the reason that combustion and radiation are characterized by different time scales and different spatial and chemical treatments, the radiation effect is often neglected or roughly modelled. The coupling of a large eddy simulation combustion solver and a radiation solver through a dedicated language, CORBA, is investigated. Two formulations of Monte Carlo method (Forward Method and Emission Reciprocity Method) employed to resolve RTE have been compared in a one-dimensional flame test case using three-dimensional calculation grids with absorbing and emitting media in order to validate the Monte Carlo radiative solver and to choose the most efficient model for coupling. Then the results obtained using two different RTE solvers (Reciprocity Monte Carlo method and Discrete Ordinate Method) applied on a three-dimensional flame holder set-up with a correlated-k distribution model describing the real gas medium spectral radiative properties are compared not only in terms of the physical behavior of the flame, but also in computational performance (storage requirement, CPU time and parallelization efficiency). To cite this article: J. Zhang et al., C. R. Mecanique 337 (2009).     

Solution of the problem of a flow of rarefied gas of constant density around a semiinfinite plate by the integral diffusion method

Several papers [1–4] have proposed approximate diffusion models which can be used to examine the transport process in a rarefied gas where the mean free path is large and transport is not determined by the local gradient of the particular quantity.In this paper the integral diffusion model [2] is used to solve the problem of determination of the friction stress and velocity of a flow of an incompressible gas around a plane semi-infinite plate in the whole range of Knudsen numbers. The obtained solution is compared with published solutions and experimental data [9].     

Solution of the laminar duct flow problem by a boundary integral equation method

Summary A boundary integral equation method is proposed for approximate numerical and exact analytical solutions to fully developed incompressible laminar flow in straight ducts of multiply or simply connected cross-section. It is based on a direct reduction of the problem to the solution of a singular integral equation for the vorticity field in the cross section of the duct. For the numerical solution of the singular integral equation, a simple discretization of it along the cross-section boundary is used. It leads to satisfactory rapid convergency and to accurate results. The concept of hydrodynamic moment of inertia is introduced in order to easily calculate the flow rate, the main velocity, and the fRe-factor. As an example, the exact analytical and, comparatively, the approximate numerical solution of the problem of a circular pipe with two circular rods are presented. In the literature, this is the first non-trivial exact analytical solution of the problem for triply connected cross section domains. The solution to the Saint-Venant torsion problem, as a special case of the laminar duct-flow problem, is herein entirely incorporated.     

Monte Carlo analysis of dissociation and recombination behind strong shock waves in nitrogen

Computations are presented for the relaxation zone behind strong, one-dimensional shock waves in nitrogen. The analysis is performed with the direct simulation Monte Carlo method (DSMC). The DSMC code is vectorized for efficient use on a supercomputer. The code simulates translational, rotational and vibrational energy exchange and dissociative and recombinative chemical reactions. A new model is proposed for the treatment of three body recombination collisions in the DSMC technique which usually simulates binary collision events. The new model represents improvement over previous models in that it can be employed with a large range of chemical rate data, does not introduce into the flow field troublesome pairs of atoms which may recombine upon further collision (pseudo-particles) and is compatible with the vectorized code. The computational results are compared with existing experimental data. It is shown that the derivation of chemical rate coefficients must account for the degree of vibrational nonequilibrium in the flow. A nonequilibrium chemistry model is employed together with equilibrium rate data to compute successfully the flow in several different nitrogen shock waves.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1990.     

Solution of nonaxisymmetric three-dimensional thermoplasticity problem by the secondary-stress method

The nonaxisymmetric thermoplasticity problem for a laminar solid of revolution is solved by successive approximation. The theory of deformation along slightly curved trajectories linearized by the secondary-stress method is employed. Numerical examples show that the proposed procedure reduces the number of successive approximations by 20% relative to the traditional approach. S. P. Timoshenko Institute of Mechanics, National Academy of Sciences of Ukraine, Kiev. Translated from Prikladnaya Mekhanika, Vol. 35, No. 12, pp. 19–25, December, 1999.