Statistic modeling and Monte Carlo methods
Non-stationary process of penetration of neutral interstellar atoms of Hydrogenium (H-atoms) into the solar system is considered. The model of the media is defined by two-shock axially symmetric gas dynamical model of interaction of the solar wind plasma with local interstellar medium (LISM) plasma. The given parameters are: the number density, the temperature and the vector of mean velocity of the plasma protons. The velocity distribution of protons is shifted Maxwellian. The plasma parameters are assumed to be stationary and independent of the H-atoms. There are two types of interactions of H-atoms with the medium: a resonance charge-exchange process being equivalent to scattering process and a photoionization process being equivalent to absorption process. The phase density of H-atoms satisfies to the non-stationary kinetic transport equation.
The Monte Carlo algorithms and computer program for calculation of time-spatial-velocity distribution of H-atoms in plasma medium were developed on the basis of a technique of "fiction collisions". For computation of the appropriate functionals of the phase density some new Monte Carlo estimators were constructed. For functionals concentrated on special surfaces (disks, flat rings) the randomized estimators of the generalized method of "mathematical expectations" were used. The auxiliary random samplings of different types were used: points on a surface, vectors in a cone and in all these cases auxiliary random modules of the velocity vector were sampled. The process of absorption was taken into account by the weight with the technique of modelling "on scatterings".
Note. Abstracts are published in author's edition
Mail to Webmaster
|||Home Page| |English Part|||
Go to Home