逆正向蒙特卡罗方法及应用

大气与环境光学学报 ›› 2022, Vol. 17 ›› Issue (5): 521-532.

逆正向蒙特卡罗方法及应用

施国栋

上海财经大学金融学院, 上海 200433

收稿日期:2021-04-28 修回日期:2022-06-25 出版日期:2022-09-28 发布日期:2022-10-17
通讯作者: E-mail: ivysgd@sandau.edu.cn E-mail:ivysgd@126.com
作者简介:施国栋 (1980 - ), 江苏东台人, 博士, 讲师, 主要从事算法方面的研究。 E-mail: ivysgd@sandau.edu.cn
基金资助:
上海财经大学中央高校基本科研业务费专项资金

Reverse forward Monte Carlo method and its application

SHI Guodong

School of Finance, Shanghai University of Finance and Economics, Shanghai 200433, China

Received:2021-04-28 Revised:2022-06-25 Published:2022-09-28 Online:2022-10-17
Supported by:
Supported by the General Program of China Postdoctoral Science Foundation (中国博士后科学基金面上项目, 2020M681253)

摘要/Abstract

摘要： 在打破传统蒙特卡罗方法在光线跟踪计算时各个子过程连贯性的基础上, 发展了灵活的逆正向蒙特卡罗法, 进而利用该方法研究了介质的矢量辐射传输问题。通过建立一维梯度折射率半透明介质模型, 模拟了大气层折射率分布, 得到了考虑偏振时介质自身的 Stokes 矢量。研究结果表明, 梯度折射率介质不同表面的发射率存在一定的区别, 散射性介质在考虑偏振时表面发射率的结果与标量方法存在一定的差异。逆正向蒙特卡罗方法可用以解决矢量辐射传输这类复杂问题。

关键词: 蒙特卡罗法, 梯度折射率, 斯托克斯矢量

Abstract: On the basis of breaking the coherence of each subprocess in traditional Monte Carlo method in ray tracing, a flexible reverse forward Monte Carlo method is developed to study vector radiative transfer in media. The distribution of atmospheric refractive index is simulated by establishing a one-dimensional gradient refractive index translucent medium model, and the Stokes vector of the medium considering polarization is studied. The results show that the emissivity of different surfaces of graded-index medium is different, and the emissivity of scattering medium when polarization is considered is different from that of scalar method. It is shown that the reverse forward Monte Carlo method can be used to solve the complex problems such as vector radiative transfer.

Key words: Monte Carlo method, gradient index, Stokes vector

中图分类号:

O435.1

施国栋. 逆正向蒙特卡罗方法及应用[J]. 大气与环境光学学报, 2022, 17(5): 521-532.

SHI Guodong. Reverse forward Monte Carlo method and its application[J]. Journal of Atmospheric and Environmental Optics, 2022, 17(5): 521-532.

参考文献 28

[1]	Liu L H. Benchmark numerical solutions for radiative heat transfer in two-dimensional medium with graded index distribution
[J]	Journal of Quantitative Spectroscopy and Radiative Transfer, 2006, 102(2): 293-303.
[2]	Huang Y, Wang K, Wang J. Two methods for solving radiative transfer in a two-dimensional graded-index medium [J]. Journal
	of Thermophysics and Heat Transfer, 2009, 23(4): 703-710.
[3]	Modest M F. Backward Monte Carlo simulations in radiative heat transfer [J]. Journal of Heat Transfer, 2003, 125(1): 57-62.
[4]	Lau C W, Watson K M. Radiation transport along curved ray paths [J]. Journal of Mathematical Physics, 1970, 11(11): 3125-
	3137.
[5]	Zhao J M, Tan J Y, Liu L H. On the derivation of vector radiative transfer equation for polarized radiative transport in graded
	index media [J]. Journal of Quantitative Spectroscopy and Radiative Transfer, 2012, 113(3): 239-250.
[6]	Zhao J M, Tan J Y, Liu L H. Monte Carlo method for polarized radiative transfer in gradient-index media [J]. Journal of
	Quantitative Spectroscopy and Radiative Transfer, 2015, 152: 114-126.
[7]	Ben X, Yi H L, Tan H P. Approximate solution to vector radiative transfer in gradient-index medium [J]. Applied Optics, 2014,
53	(3): 388-401.
[8]	Kravtsov Y A, Orlov Y I. Geometrical Optics of Inhomogeneous Media [M]. Berlin, Heidelberg: Springer Berlin Heidelberg,
	1990.
[9]	Ambirajan A, Look D C. A backward Monte Carlo study of the multiple scattering of a polarized laser beam [J]. Journal of
	Quantitative Spectroscopy and Radiative Transfer, 1997, 58(2): 171-192.
[10]	Collins D G, Blattner W G, Wells M B, ¨ et al. Backward Monte Carlo calculations of the polarization characteristics of the
	radiation emerging from spherical-shell atmospheres [J]. Applied Optics, 1972, 11(11): 2684-2696.
[11]	Gratiy S L, Walker A C, Levin D A, et al. Multi-wavelength simulations of atmospheric radiation from Io with a 3-D sphericalshell backward Monte Carlo radiative transfer model [J]. Icarus, 2010, 207(1): 394-408.
[12]	Gay B, Vaillon R, Menguc¸ M P. Polarization imaging of multiply-scattered radiation based on integral-vector Monte Carlo ¨
	method [J]. Journal of Quantitative Spectroscopy and Radiative Transfer, 2010, 111(2): 287-294.
[13]	Emde C, Buras R, Mayer B, et al. The impact of aerosols on polarized sky radiance: Model development, validation, and
	applications [J]. Atmospheric Chemistry and Physics, 2010, 10(2): 383-396.
[14]	Oikarinen L. Polarization of light in UV-visible limb radiance measurements [J]. Journal of Geophysical Research: Atmospheres, 2001, 106(D2): 1533-1544.
[15]	Garc´ıa Munoz A, Mills F P. Pre-conditioned backward Monte Carlo solutions to radiative transport in planetary atmospheres ˜
[J]	Astronomy & Astrophysics, 2015, 573: A72.
[16]	Zhu K Y, Huang Y, Wang J. Curved ray tracing method for one-dimensional radiative transfer in the linear-anisotropic scattering
	medium with graded index [J]. Journal of Quantitative Spectroscopy and Radiative Transfer, 2011, 112(3): 377-383.
[17]	Huang Y, Liang X G, Xia X L. Monte Carlo simulation of radiative transfer in scattering, emitting, absorbing slab with gradient
	index [J]. Journal of Quantitative Spectroscopy and Radiative Transfer, 2005, 92(1): 111-120.
[18]	Chandrasekhar S. Radiative Transfer [M]. NewYork: Oxford University Press, 1950.
[19]	Kattawar G W, Adams C N. Stokes vector calculations of the submarine light field in an atmosphere-ocean with scattering
	according to a Rayleigh phase matrix: Effect of interface refractive index on radiance and polarization [J]. Limnology and
	Oceanography, 1989, 34(8): 1453-1472.
[20]	Zhai P W, Hu Y X, Trepte C R, et al. A vector radiative transfer model for coupled atmosphere and ocean systems based on
	successive order of scattering method [J]. Optics Express, 2009, 17(4): 2057.
[21]	Zhai P W, Hu Y X, Chowdhary J, et al. A vector radiative transfer model for coupled atmosphere and ocean systems with a
	rough interface [J]. Journal of Quantitative Spectroscopy and Radiative Transfer, 2010, 111(7/8): 1025-1040.
[22]	Hulst H C, Van De Hulst H C. Light Scattering by Small Particles [M]. NewYork: Dover publications, 1981.
[23]	Tynes H H, Kattawar G W, Zege E P, et al. Monte Carlo and multicomponent approximation methods for vector radiative
	transfer by use of effective mueller matrix calculations [J]. Applied Optics, 2001, 40(3): 400-412.
[24]	Barkstrom B R. An efficient algorithm for choosing scattering directions in Monte Carlo work with arbitrary phase functions
[J]	Journal of Quantitative Spectroscopy and Radiative Transfer, 1995, 53(1): 23-38.
[25]	Whitney B A . Monte Carlo radiative transfer [J]. Bulletin of the Astronomical Society of India, 2011, 39(1): 101-127.
[26]	Ben X, Yi H L, Tan H P. Polarized radiative transfer in an arbitrary multilayer semitransparent medium [J]. Applied Optics,
20	14, 53(7): 1427-1441.
[27]	Green R M. Spherical Astronomy [M]. New York: Cambridge University Press, 1985.
[28]	Abdallah P B, Dez V L. Thermal emission of a semi-transparent slab with variable spatial refractive index [J]. Journal of
	Quantitative Spectroscopy and Radiative Transfer, 2000, 67(3): 185-198.