-
Compton and inverse Compton scattering from relativistic Maxwellian electrons both have an important feature, i.e. calculating the radiation transport in high-temperature and full-ionized plasma. Description and evaluation of relativistic photon-Maxwellian electron scattering are numerically complex and computationally time consuming. A Monte Carlo method is proposed to simulate photon scattering with relativistic Maxwellian electron and compute the scattering cross-sections. To compute the total cross-section of a photon of energy hν interacting with electrons at temperature Te in the laboratory coordinate, the calculation steps of Monte Carlo scheme are described as follows. The first step is to sample the velocity of an electron, the directions are isotropically sampled, and the speed is sampled from the relativistic Maxwellian distribution at temperature Te. The second step is to transform the photon energy hν into the photon energy in the coordinate in which the electron is at rest. The third step is to use the exact Klein-Nishina formula to compute the cross-sections. The fourth step is go back to the first step, and cycle this many times. The last step is to summarize all computed cross-sections and averaged them, and the average value is what we need. The operation and corresponding formula for each step are described in this paper. A better method of sampling the speed of a relativistic electron is expected to be found.A Monte Carlo processor is developed to compute the scattering cross-section of a photon of any energy, interacting with electrons at any temperature. To check this method, scatterings of the photons of various energies with electrons with various temperatures are simulated, and the results are compared with those from the numerical integration method. The comparison indicates that the simulated cross-sections are in pretty good agreement with those from the multiple integration method for the cases of electron temperature less than 25 keV. But unfortunately, the difference is obvious for the case of temperature more than 25 keV, and the error increases with temperature increasing. Why so? When the temperature is more than 25 keV, the sampling of electron speed is inaccurate when using the present method, which maybe results in this difference. So, we need to find a more accurate method of sampling relativistic electron speed to solve this problem in the future.
-
Keywords:
- Compton scattering /
- relativistic Maxwellian electrons /
- radiation transport calculations /
- Monte Carlo method
[1] Yu M 1996 Selected Papers of Yu Min (Beijing: Institute of Applied Physics and Computational Mathematics) p102 (in Chinese)[于敏 1996 于敏论文集(北京: 北京应用物理与计算数学研究所) 第102页]
[2] Pomraning G C 1973 The Equations of Radiation Hydrodynamics (Oxford: Pergamon Press) p183
[3] Evans R D 1955 The Atomic Nucleus (New York: McGraw-Hill Press) p677
[4] Lux I, Koblinger L 1991 Monte Carlo Particle Transport Methods: Neutron and Photon Calculations (Boston: CRC Press) p45
[5] Hastings C J 1955 Approximations for Digital Computers (Princeton: Princeton University Press) p154
[6] Wienke B R 1973 Nuclear Science and Engineering 52 247
[7] Cooper G E 1974 J. Quant. Spectr. Rad. Transfer 14 887
[8] Wienke B R 1975 J. Quant. Spectr. Rad. Transfer 15 151
[9] Wienke B R, Lathrop B L 1984 J. Comp. Phys. 53 331
[10] Brinkmann W 1984 J. Quant. Spectrosc. Radiat. Transfer 31 417
[11] Wienke B R, Hendricks J S, Booth T E 1985 J. Quant. Spectr. Rad. Transfer 33 555
[12] Wienke B R, Lathrop B L, Devaney J J 1986 Radiation Effects 94 303
[13] Prasad M K, Kershaw D S, Beason J D 1986 Applied Physics Letters 48 1193
[14] Kershaw D S 1987 J. Quant. Spectr. Rad. Transfer 38 347
[15] Shestakov A I, Kershaw D S, Prasad M K 1988 J. Quant. Spectr. Rad. Transfer 40 577
[16] Webster J B, Stephan B G, Bridgman C J 1973 Trans. Amer. Nucl. Soc. 17 574
[17] Wienke B R, Lathrop B L, Devaney J J 1984 Nuclear Sci. Eng. 88 71
[18] Booth T E, Hendricks J S 1985 Nuclear Sci. Eng. 90 248
[19] Pomraning G C 1972 J. Quant. Spectr. Rad. Transfer 12 1047
[20] Pomraning G C 1973 The Equations of Radiation Hydrodynamics (Oxford: Pergamon Press) p185
[21] Mohamed N M A 2014 Theory Probab. Appl. 58 698
[22] Wienke B R 1975 Am. J. Phys. 43 317
-
[1] Yu M 1996 Selected Papers of Yu Min (Beijing: Institute of Applied Physics and Computational Mathematics) p102 (in Chinese)[于敏 1996 于敏论文集(北京: 北京应用物理与计算数学研究所) 第102页]
[2] Pomraning G C 1973 The Equations of Radiation Hydrodynamics (Oxford: Pergamon Press) p183
[3] Evans R D 1955 The Atomic Nucleus (New York: McGraw-Hill Press) p677
[4] Lux I, Koblinger L 1991 Monte Carlo Particle Transport Methods: Neutron and Photon Calculations (Boston: CRC Press) p45
[5] Hastings C J 1955 Approximations for Digital Computers (Princeton: Princeton University Press) p154
[6] Wienke B R 1973 Nuclear Science and Engineering 52 247
[7] Cooper G E 1974 J. Quant. Spectr. Rad. Transfer 14 887
[8] Wienke B R 1975 J. Quant. Spectr. Rad. Transfer 15 151
[9] Wienke B R, Lathrop B L 1984 J. Comp. Phys. 53 331
[10] Brinkmann W 1984 J. Quant. Spectrosc. Radiat. Transfer 31 417
[11] Wienke B R, Hendricks J S, Booth T E 1985 J. Quant. Spectr. Rad. Transfer 33 555
[12] Wienke B R, Lathrop B L, Devaney J J 1986 Radiation Effects 94 303
[13] Prasad M K, Kershaw D S, Beason J D 1986 Applied Physics Letters 48 1193
[14] Kershaw D S 1987 J. Quant. Spectr. Rad. Transfer 38 347
[15] Shestakov A I, Kershaw D S, Prasad M K 1988 J. Quant. Spectr. Rad. Transfer 40 577
[16] Webster J B, Stephan B G, Bridgman C J 1973 Trans. Amer. Nucl. Soc. 17 574
[17] Wienke B R, Lathrop B L, Devaney J J 1984 Nuclear Sci. Eng. 88 71
[18] Booth T E, Hendricks J S 1985 Nuclear Sci. Eng. 90 248
[19] Pomraning G C 1972 J. Quant. Spectr. Rad. Transfer 12 1047
[20] Pomraning G C 1973 The Equations of Radiation Hydrodynamics (Oxford: Pergamon Press) p185
[21] Mohamed N M A 2014 Theory Probab. Appl. 58 698
[22] Wienke B R 1975 Am. J. Phys. 43 317
Catalog
Metrics
- Abstract views: 5084
- PDF Downloads: 81
- Cited By: 0
下载:
