The boundary element method for light scattering by ice crystals and its implementation in BEM++

Groth, Samuel P.; Baran, Anthony J.; Betcke, Timo; Havemann, Stephan; Śmigaj, Wojciech

doi:10.1016/j.jqsrt.2015.08.001

Cited by 30 publications

(28 citation statements)

References 30 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Macke et al, 1996a, b;McFarquhar and Heymsfield, 1996;Yang et al, 2000Yang et al, , 2005Yang et al, , 2013Um and McFarquhar, 2007Nousiainen et al, 2011;Baum et al, 2005Baum et al, , 2011Baran and C.-Labonnote, 2007;Ishimoto et al, 2012b;Liu et al, 2014a). Numerous light-scattering computation methods have been employed to calculate the single-scattering properties of the various ice particles, including the finitedifference time-domain (FDTD) method (Yee, 1966;Yang and Liou, 1998a;Sun et al, 1999, Ishimoto et al, 2012a, the T-matrix (Baran et al, 2001;Bi and Yang, 2014a, b), the discrete dipole approximation method (Purcell and Pennypacker, 1973;Draine and Flatau, 1994;Yurkin et al, 2007), the boundary element method (Mano, 2000;Groth et al, 2015), the pseudo-spectral time-domain method (Liu, 1997(Liu, , 1998Chen et al, 2008;Liu et al, 2012), the surface-integral equation method (Nakajima et al, 2009), the improved geometric optics method (IGOM) (Yang and Liou, 1996), geometric optics integral equation (GOIE) (Yang and Liou, 1996;Ishimoto et al, 2012a), and the ray-tracing geometric optics method (GOM) Liou, 1989, 1993;Macke, 1993;Macke et al, 1996a;Yang and Liou, 1998b;Masuda et al, 2012).…”

Section: Introductionmentioning

confidence: 99%

Investigation of ice particle habits to be used for ice cloud remote sensing for the GCOM-C satellite mission

et al. 2016

View full text Add to dashboard Cite

Abstract. In this study, various ice particle habits are investigated in conjunction with inferring the optical properties of ice clouds for use in the Global Change Observation Mission-Climate (GCOM-C) satellite programme. We develop a database of the single-scattering properties of five ice habit models: plates, columns, droxtals, bullet rosettes, and Voronoi. The database is based on the specification of the Second Generation Global Imager (SGLI) sensor on board the GCOM-C satellite, which is scheduled to be launched in 2017 by the Japan Aerospace Exploration Agency. A combination of the finite-difference time-domain method, the geometric optics integral equation technique, and the geometric optics method is applied to compute the single-scattering properties of the selected ice particle habits at 36 wavelengths, from the visible to the infrared spectral regions. This covers the SGLI channels for the size parameter, which is defined as a single-particle radius of an equivalent volume sphere, ranging between 6 and 9000 µm. The database includes the extinction efficiency, absorption efficiency, average geometrical cross section, single-scattering albedo, asymmetry factor, size parameter of a volume-equivalent sphere, maximum distance from the centre of mass, particle volume, and six nonzero elements of the scattering phase matrix. The characteristics of calculated extinction efficiency, single-scattering albedo, and asymmetry factor of the five ice particle habits are compared. Furthermore, sizeintegrated bulk scattering properties for the five ice particle habit models are calculated from the single-scattering database and microphysical data. Using the five ice particle habit models, the optical thickness and spherical albedo of ice clouds are retrieved from the Polarization and Directionality of the Earth's Reflectances-3 (POLDER-3) measurements, recorded on board the Polarization and Anisotropy of Reflectances for Atmospheric Sciences coupled with Observations from a Lidar (PARASOL) satellite. The optimal ice particle habit for retrieving the SGLI ice cloud properties is investigated by adopting the spherical albedo difference (SAD) method. It is found that the SAD is distributed stably due to the scattering angle increases for bullet rosettes with an effective diameter (D eff ) of 10 µm and Voronoi particles with D eff values of 10, 60, and 100 µm. It is confirmed that the SAD of small bullet-rosette particles and all sizes of Voronoi particles has a low angular dependence, indicating that a combination of the bullet-rosette and Voronoi models is sufficient for retrieval of the ice cloud's spherical albedo and optical thickness as effective habit models for the SGLI Published by Copernicus Publications on behalf of the European Geosciences Union. 12288H. Letu et al.: Investigation of ice particle habits to be used for ice cloud remote sensing sensor. Finally, SAD analysis based on the Voronoi habit model with moderate particle size (D eff = 60 µm) is compared with the conventional general habit mixture mode...

show abstract

Section: Introductionmentioning

confidence: 99%

Investigation of ice particle habits to be used for ice cloud remote sensing for the GCOM-C satellite mission

et al. 2016

View full text Add to dashboard Cite

show abstract

“…Early applications of BEM to the simulation of light scattering by simple ice crystals include [37] and [38]. More comprehensive studies of complex crystal shapes (including hexagonal columns with conventional and stepped cavities, bullet rosettes and Chebyshev ice particles) have been given recently by Groth et al [34] and Baran and Groth [35], using the open-source BEM software library Bempp [39], available at www.bempp.com. In [34,35] it was shown (by comparison with a T-matrix method) that, using a discretization with at least 10 boundary elements per wavelength, Bempp can compute far-field quantities with 1% relative error, with reciprocity also satisfied to a similar accuracy.…”

Section: Introductionmentioning

confidence: 99%

“…More comprehensive studies of complex crystal shapes (including hexagonal columns with conventional and stepped cavities, bullet rosettes and Chebyshev ice particles) have been given recently by Groth et al [34] and Baran and Groth [35], using the open-source BEM software library Bempp [39], available at www.bempp.com. In [34,35] it was shown (by comparison with a T-matrix method) that, using a discretization with at least 10 boundary elements per wavelength, Bempp can compute far-field quantities with 1% relative error, with reciprocity also satisfied to a similar accuracy. The results in [34,35] were obtained using a standard desktop machine and were limited to size parameters up to 15, but Bempp supports parallelization and hence with HPC architectures much larger problems can be tackled.…”

Section: Introductionmentioning

confidence: 99%

“…In [34,35] it was shown (by comparison with a T-matrix method) that, using a discretization with at least 10 boundary elements per wavelength, Bempp can compute far-field quantities with 1% relative error, with reciprocity also satisfied to a similar accuracy. The results in [34,35] were obtained using a standard desktop machine and were limited to size parameters up to 15, but Bempp supports parallelization and hence with HPC architectures much larger problems can be tackled. We end this brief literature review of BEM in atmospheric physics applications by mentioning recent work by Yu et al [36], in which boundary integral equations were applied to scattering by multiple dielectric particles illuminated by unpolarized high-order Bessel vortex beams, and also by Groth et al [40], where a "hybrid numerical-asymptotic" BEM was presented for 2D high frequency scalar dielectric scattering problems which achieves fixed accuracy with frequency-independent computational cost.…”

Section: Introductionmentioning

confidence: 99%

See 1 more Smart Citation

Calderón preconditioning of PMCHWT boundary integral equations for scattering by multiple absorbing dielectric particles

Kleanthous

Betcke

Hewett

et al. 2019

Journal of Quantitative Spectroscopy and Radiative Transfer

View full text Add to dashboard Cite

We consider the simulation of electromagnetic scattering by single and multiple isotropic homogeneous dielectric particles using boundary integral equations. Galerkin discretizations of the classical Poggio-Miller-Chang-Harrington-Wu-Tsai (PMCHWT) boundary integral equation formulation provide accurate solutions for complex particle geometries, but are well-known to lead to ill-conditioned linear systems. In this paper we carry out an experimental investigation into the performance of Calderón preconditioning techniques for single and multiple absorbing obstacles, which involve a squaring of the PMCHWT operator to produce a well-conditioned second-kind formulation. For single-particle scattering configurations we find that Calderón preconditioning is actually often outperformed by simple "mass-matrix" preconditioning, i.e. working with the strong form of the discretized PMCHWT operator. In the case of scattering by multiple particles we find that a significant saving in computational cost can be obtained by performing block-diagonal Calderón preconditioning in which only the self-interaction blocks are preconditioned. Using the boundary element software library Bempp (www.bempp.com) the numerical performance of the different methods is compared for a range of wavenumbers, particle geometries and complex refractive indices relevant to the scattering of light by atmospheric ice crystals.

show abstract

“…So, it can be found that light propagation in biotissues is a complicated process including photon absorption and scattering. Because of complexity of light propagation, the numerical or statistical methods such as Monte Carlo (László & Manuela, ), finite element analysis (Lin & Jao, ), boundary element method (Groth, Baran, Betcke, Havemann, & Śmigaj, ) and transmission diffusion approximation equation (Zhan, Tang, Li, Wang, & Li, ) are frequently used to modeling and analyze the transmission characteristic of light in bio‐tissues. However, expensive instrument and complex calculation limit application of these methods.…”

Section: Introductionmentioning

confidence: 99%

Near‐infrared light penetration depth analysis inside melon with thick peel by a novel strategy of slicing combining with least square fitting method

Zhang

2018

J Food Process Engineering

View full text Add to dashboard Cite

The light penetration depth in melon was explored by near‐infrared (NIR) diffuse reflectance spectroscopy (500–1,100 nm) combining a novel strategy of slicing and least square fitting. The depths of peak wavelength and effective waveband (720–880 nm) based on three source‐detector distances (d = 10, 20, and 30 mm) were calculated and compared by least square fitting. The results indicated that the obtained depth were similar both slicing method and least square fitting. The depth values of melon fruits were 12.3, 14.9, and 19.0 mm for peak wavelength, and the depth value ranges were 11.7–12.4, 14.2–15.8, and 16.8–21.0 mm for the effective waveband of 720–880 nm when d was 10, 20, and 30 mm, respectively. The overall results demonstrated that least square fitting method could be accurately applied to calculate the light penetration depth inside melon and the depth value increased with increasement of d. And 20 mm source‐detector distance could be the optimal selection. This study can provide the valuable reference for development of spectral instrument used to effectively detect the internal quality of fruits, especially for fruits with thick peel. Practical application The light penetration depth in fruit tissues was very important for obtaining the more effective spectral information used to accurately and robustly assess the internal quality of fruits. Currently, there was no effective and rapid technology to measure the light penetration depth of fruits, especially fruits with thick peel. Hence, in this research, we provided a novel method of slicing and least square fitting to calculate the light penetration depth in melon based on different source‐detector distance. The result indicated that the proposed method could be applied to calculate the penetration depth in melon and provided a valuable reference to develop the spectral equipment.

show abstract

The boundary element method for light scattering by ice crystals and its implementation in BEM++

Cited by 30 publications

References 30 publications

Investigation of ice particle habits to be used for ice cloud remote sensing for the GCOM-C satellite mission

Investigation of ice particle habits to be used for ice cloud remote sensing for the GCOM-C satellite mission

Calderón preconditioning of PMCHWT boundary integral equations for scattering by multiple absorbing dielectric particles

Near‐infrared light penetration depth analysis inside melon with thick peel by a novel strategy of slicing combining with least square fitting method

Contact Info

Product

Resources

About