Design and Verification of Geometric Optical Approximation Model Based on Ray Tracing Technique and Electromagnetic Equivalence Principle
Zhao Jiaqi,Hu Shuai,Dang Ruijun,Liu Lei,Liu Xichuan,Xiao Yao
DOI: https://doi.org/10.3788/aos231247
2024-01-01
Acta Optica Sinica
Abstract:Objective Cirrus clouds and some top clouds typically contain a large number of ice crystal particles, which have a strong scattering and absorption effect for visible and infrared radiation and play an important role in the balance of the atmospheric energy budget. Understanding the radiation characteristics of cirrus clouds must start with understanding the single scattering characteristics of non-spherical ice crystals. However, due to the irregular shape of ice crystal particles such as hexagonal columns and bullet flowers, the calculation of their scattering characteristics has significant uncertainty, making them one of the most uncertain factors in radiation transfer simulation. Therefore, accurately simulating the light scattering process of ice crystals is a current research focus. At present, for ice crystal particles with small and medium-sized parameters, scattering calculation models for non- spherical particles are gradually developing, such as discrete dipole approximation (DDA), invariant imbedding T-matrix (IIMT matrix), and finite-difference time- domain (FDTD) method. The DDA and FDTD models are mainly employed for particles with size parameters less than 40, and the IIMT matrix model is mainly for particles with size parameters less than 100. For large- size parameter particles, most of the traditional geometric optical approximation models are adopted in China to calculate their scattering parameters. However, this method is based on the fact that the particle size is much larger than the incident light wavelength, and the effect becomes worse during calculating particles with size parameters less than 300. Therefore, it is extremely important to independently develop the geometrical optics approximation model and expand its applicable size parameter range. In this regard, we develop an improved geometric-optical approximation model by combining the electromagnetic equivalence principle with ray tracing techniques and considering the effects of diffraction and particle absorption. Methods Firstly, a ray tracing algorithm considering polarization is constructed using Monte Carlo technology to simulate the beam reflection and refraction processes and to track the propagation direction and electric field vector of the beam. Secondly, based on ray tracing, a calculation scheme for far-field scattering electric field is studied. For particles with size parameters greater than 300, a direct ray statistics scheme is adopted, and for small particles, a far- field electric field calculation scheme based on the electromagnetic equivalence principle is designed. Furthermore, the study of diffraction calculation techniques for irregular cross-sections involves projecting particles in the direction of incident light and numerically solving the diffraction equation to obtain the diffraction electric field. Finally, based on the ray tracing electric field, the far field electric field of scattered light, and the diffraction electric field, the calculation scheme of particle scattering characteristics is designed to realize the calculation of particle extinction, absorption cross section, S-matrix, and other scattering parameters. The calculation model of the edge diffraction effect is built to realize the compensation and correction of extinction and absorption cross section. Results and Discussions We compare the calculation results of IGOA with the physical geometrical optics approximation model (PGOM) developed by Yang et al. and the independently developed IIMT model, with the phase matrix elements computed by IGOA and PGOM compared. The particle is a hexagonal column with an aspect ratio of 1.0, and the size parameter is 100 and 300. The refractive index is 1.308+ i1.43x10- 9, and the incident light wavelength is 0.65 mu m. The results of the IGOA model and the PGOM model are basically consistent, and their scattering phase function (P 11) varying with scattering angle is basically in sound agreement, especially in the forward scattering angle (0 degrees. 90 degrees), which is more significant. Generally, P 12, P 22, P33, P34, and P 44 have a high degree of agreement. A hexagonal prism particle with a size parameter of 300 is calculated, with a bottom length of 572.95779 mu m, a height of 286.47889 mu m, an incident light wavelength of 12 mu m, and a particle complex refractive index of 1.2762+i0.4133. The calculated results are compared with those of PGOM, as shown in Fig. 7. For strongly absorbing particles, the calculation results of IGOA and PGOM maintain high consistency at various scattering angles, especially for polarization characteristics, where the calculation curves of the two models basically coincide. This means that the IGOA model can also achieve a high level of computational accuracy for strongly absorbing particles. To analyze the calculation results of smaller parameter particles with different shapes, we calculate the scattering phase matrices of hexagonal and dodecagonal prism particles by IGOA and IIMT models, with sound consistency between the two models. Conclusions The geometrical optics approximation model is an important tool for calculating the light scattering characteristics of large size parameter particles, but it is based on the fact that the incident light wavelength is longer than the particle size. For particles with size parameters of 100. 300, the simulation accuracy is relatively poor because they are just in the transition scale range from physical optics to geometrical optics. To solve this problem, we combine the ray tracing technology with the electromagnetic equivalence principle, consider the strong absorption and diffraction of particles and other factors, and independently compile the improved IGOA model. The calculation results are compared with those of the IIMT and PGOM models, and the following conclusions are drawn: 1) For particles with size parameters ranging from 80 to 300, the calculation results of the IGOA model are highly consistent with those of the PGOM and IIMT models, which indicates that the calculation results of the model have high accuracy. Generally, the calculation accuracy of the model is maintained at a high level for particles with different size parameters and complex refractive indexes. 2) For scattering characteristic calculation of strongly absorbing particles, the IGOA model shows sound simulation performance, and its calculation results basically coincide with the calculation curve of PGOM.