Next Article in Journal
GCM Solver (Ver. 3.0): A Mathematica Notebook for Diagonalization of the Geometric Collective Model (Bohr Hamiltonian) with Generalized Gneuss–Greiner Potential
Next Article in Special Issue
Numerical Investigation of Film Cooling Enhancement Using an Upstream Sand-Dune-Shaped Ramp
Previous Article in Journal
Advanced Bio-Inspired System for Noninvasive Cuff-Less Blood Pressure Estimation from Physiological Signal Analysis
Previous Article in Special Issue
The Hydraulic Cavitation Affected by Nanoparticles in Nanofluids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Determination of Aerosol Size Distribution from Angular Light-Scattering Signals by Using a SPSO-DE Hybrid Algorithm

Aero-engine Thermal Environment and Structure Key Laboratory of Ministry of Industry and Information Technology, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
*
Author to whom correspondence should be addressed.
Submission received: 6 July 2018 / Revised: 19 August 2018 / Accepted: 27 August 2018 / Published: 29 August 2018
(This article belongs to the Special Issue Computational Heat, Mass and Momentum Transfer)

Abstract

:
The comparison of the angular light-scattering method (ALSM) and the spectral extinction method (SEM) in solving the inverse problem of aerosol size distribution (ASD) are studied. The inverse problem is solved by a SPSO-DE hybrid algorithm, which is based on the stochastic particle swarm optimization (SPSO) algorithm and differential evolution (DE) algorithm. To improve the retrieval accuracy, the sensitivity analysis of measurement signals to characteristic parameters in ASDs is studied; and the corresponding optimal measurement angle selection region for ALSM and optimal measurement wavelength selection region for SEM are proposed, respectively. Results show that more satisfactory convergence properties can be obtained by using the SPSO-DE hybrid algorithm. Moreover, short measurement wavelengths and forward measurement angles are beneficial to obtaining more accurate results. Then, common monomodal and bimodal ASDs are estimated under different random measurement errors by using ALSM and SEM, respectively. Numerical tests show that retrieval results by using ALSM show better convergence accuracy and robustness than those by using SEM, which is attributed to the distribution of the objective function value. As a whole, considering the convergence properties and the independence on prior optical information, the ALSM combined with SPSO-DE hybrid algorithm provides a more effective and reliable technique to obtain the ASDs.

1. Introduction

Atmospheric aerosols are important drivers of air quality, atmospheric heat radiative balance, human health and climate, etc. Generally, there are three vitally important properties of aerosols, i.e., aerosol size distributions (ASDs), aerosol optical depths (AODs) and Ångström exponents [1,2]. The knowledge of ASDs is essential to a better understanding of their potential contributions to climate change, source identification and environmental quality, necessary to control pollution and preserve citizen health. Therefore, in situ measuring the ASDs presented in the atmosphere is of great necessity [3,4]. Nowadays, although there are several global ground-based aerosol observation networks and satellite instrument networks that have been established to study the properties of the atmospheric aerosols, e.g., AERONET, MODIS [5,6,7,8], determining the ASDs accurately is still regarded as an unsolved problem and needs further research.
During the past several years, various types of methods have been proposed to study the ASDs, e.g., the aerodynamic method, the optical measurement method, electrical mobility and condensation method, electrical sensing zone method and the electron microscopy method [9]. Owing to offering a useful and effective approach, reliable measurement results and broad size ranges, the optical measurement method coupled with inverse methods has drawn much attention in the field of determining the ASDs [10]. Usually the optical measurement method could be divided into several categories, such as spectral extinction method (SEM), angular light-scattering method (ALSM) and optical depth measurement method, etc. During our previous work, the ASDs were reconstructed effectively by using SEM, and the retrieval accuracy was improved by increasing the number of measurement wavelengths [11]. However, all the satisfactory results obtained by using SEM were based on having sufficient spectral measurement information about the aerosol dispersion system beforehand, e.g., the multi-spectral transmittance signals. To obtain these information measurements, the multi-spectral light source and detector should be needed, and the multi-spectral optical constants of aerosol should be known, which may increase the difficulty and the cost of experiment. Unlike the SEM, only single spectral information of aerosol is required by using ALSM, and the ASDs can be recovered by measuring multi-angle light-scattering signals [12,13].
In the present study, the comparison of ALSM and SEM in retrieval of ASDs are studied. The inverse problem is solved by using a SPSO-DE hybrid algorithm, which is based on the stochastic particle swarm optimization (SPSO) algorithm and differential evolution (DE) algorithm. The research is organized as follows. First, the principles of ALSM and SPSO-DE hybrid algorithm are introduced, and that of SEM, available in Reference [14], is not repeated here. Then, the sensitivity analysis of optical measurement signals to characteristic parameters in ASDs is studied, and the corresponding optimal measurement angle selection region for ALSM and optimal measurement wavelength selection region for SEM are proposed. Finally, common monomodal and bimodal ASDs are retrieved by SEM and ALSM, respectively, and the main conclusions and perspectives are provided.

2. Methodology

2.1. The principle of the ALSM

According to the principle of the ALSM, when a suspension particle system is illuminated by a parallel light beam with intensity I 0 and the refractive index of the particle is different from those of the surrounding medium, the angular light-scattering intensity I ( θ ) can be written as:
I ( θ ) = π N 0   λ D min D max i ( θ , m λ , D , λ ) f ( D ) d D
i ( θ , m λ , D , λ ) = i 1 ( θ , m λ , D , λ ) + i 2 ( θ , m λ , D , λ )   8 π 2 λ 2 I 0
where i 1 ( θ , m λ , D , λ ) and i 2 ( θ , m λ , D , λ ) can be derived in terms of the Mie scattering functions [12,15]. Detailed information about the ALSM can be found in References [12,13,16].

2.2. The Principle of the SPSO

The standard PSO algorithm, introduced by Eberhart and Kennedy [17] in 1995, firstly, can be easily implemented by adjusting a small number of parameters. However, in the standard PSO algorithm, there is only one term which regulates historic velocity during the evolution process, and therefore the particle can only fly along one direction until it finds a better position, which will result in premature convergence. To overcome the difficulty, the SPSO algorithm was introduced, and showed better retrieval accuracy in solving the problem of estimating the ASD and other inverse radiation problems [18,19,20].
In the SPSO algorithm, the velocity-update formula and the position-update formula can be described as [18]:
V i ( t + 1 ) = C 1 r a n d 1 [ P i ( t ) X i ( t )   ] + C 2 r a n d 2 [ P g ( t ) X i ( t ) ]
X i ( t + 1 ) = X i ( t ) + V i ( t + 1 )  
where C 1 and C 2 denote the positive acceleration coefficients; r a n d 1 and r a n d 2 are the random numbers selected within [0, 1]; P i ( t ) and P g ( t ) denote the personal best position and global best position at the t th iteration, respectively; X i ( t ) denotes the t th iteration for old position, while X i ( t + 1 ) is the ( t + 1 ) th iteration for the new one. Detailed information of the SPSO is available in Reference [18].

2.3. The Principle of the SPSO-DE Hybrid Algorithm

In this study, to avoid the swarm diversity declining rapidly and the local optima, the Differential Evolution (DE) algorithm is proposed to modify the personal best food source positions. The detail of the DE is available in Reference [21,22,23]. During the DE algorithm, the mutant vector can be derived as [23]:
P i , mut   ( t ) = P k ( t ) + F · [ P j ( t ) P m ( t ) ]
where k , l , m are random integers uniformly selected from [ 1 , N s ] , which are different from each other; F [ 0 , 2 ] , a real and constant value, denotes the mutant factor. After the crossover process, the trial vector P i , tri ( t ) = ( p i 1 , tri , p i 2 , tri , p i 3 , tri , p i N p , tri ) can be expressed as [23]:
p i j , tri   ( t ) = { p i j , mut ( t )       if r a n d 3 ( j ) C R or j = r n b r ( i ) p i j ( t )              if r a n d 3 ( j ) > C R or j r n b r ( i )
where r a n d 3 , a uniform random number, is selected within [0, 1]; C R [ 0 , 1 ] is the user-defined crossover constant; r n b r [ 1 , N s ] is a random integer. The detailed search procedures of the SPSO-DE hybrid algorithm are shown below.
Step 1.
Initialize the system control variables of SPSO-DE hybrid algorithm, such as N s , N c , N p , ε , C 1 , C 2 , F and C R .
Step 2.
Obtain the initial particles’ position and velocity for all the swarm randomly; calculate and evaluate the corresponding objective function value F obj , i of each particle; record the initial personal best position P i and objective function value p b e s t i of each particle; record the global best position P g and objective function value g b e s t of the swarm.
Step 3.
Finish the iteration, if the calculation and iteration satisfy the stop criteria:
(i)
g b e s t is less than the tolerance ε;
(ii)
iteration number is not less than the user-defined iteration limit N c .
Step 4.
Update and obtain the new particle velocity and position according to Equations (3) and (4).
Step 5.
Calculate new F obj , i and evaluate the new positions for each particle.
Step 6.
Compare the old F obj , i with the new one for each particle. Update the corresponding particle information if the new F obj , i is superior
Step 7.
Use the DE algorithm to modify the personal best food source positions according to Equations (5) and (6).
Step 8.
Determine if the current personal best position and global best position are superior than the old ones, if so, update them. Loop to Step 3.

3. The Aerosol Size Distribution and Optical Constants

The mathematical representations of particle size distribution functions studied here are shown below [11]:
f L - N   mono ( D ) = 1 2 π D ln σ × exp [ 1 2 ( ln D ln D ¯ ln σ ) 2 ]
f Gamma   mono ( D ) = D α × exp ( β D γ )
f L - N   bi ( D ) = n × 1 2 π D ln σ 1 × exp [ 1 2 ( ln D ln D ¯ 1 ln σ 1 ) 2 ] + ( 1 n ) × 1 2 π D ln σ 2 × exp [ 1 2 ( ln D ln D ¯ 2 ln σ 2 ) 2 ]
f Gamma   bi ( D ) = n × D α 1 × exp ( β 1 D γ 1 ) + ( 1 n ) × D α 2 × exp ( β 2 D γ 2 )
where the superscript ‘mono’ and ‘bi’ denote the monomodal and bimodal distributions; the superscript ‘L-N’ and ‘Gamma’ denote the Logarithmic Normal distribution and the Gamma distribution. The detailed description of these functions is available in our previous work [11]. To meet the need of optical measurement methods, the particle size range is set within [0.001, 10.0] μm [24]. Table 1 lists the true values of the ASDs mentioned above.
According to the standard radiation atmospheric model (SRA) of atmospheric aerosols, which was recommended by the radiation committee of the International Association of Meteorology and Atmospheric Physics (IAMAP) in 1983, the atmospheric aerosols with different composition were divided into six categories, i.e., dust-like aerosol, water-soluble aerosol, oceanic aerosol, soot aerosol, volcanic ash aerosol and sulfuric acid aerosol [25]. Later, D’Almeida and coworkers [26] proposed a new kind of global aerosol model and added the mineral aerosol into the model, which represents the aerosol of one-third of global terrestrial deserts and grassland (arid or semi-arid) areas. Figure 1 depicts the spectral optical constants of four kinds of common aerosols, i.e., dust-like aerosol, water-soluble aerosol, oceanic aerosol and soot aerosol, which are available in Reference [25,26]. In the present work, the spectral optical constants are selected according to soot aerosol.

4. Sensitivity Analysis of Optical Measurement Signals to Characteristic Parameters in ASDs

Because the inverse accuracy of ASDs is relevant to the information obtained from the optical measurement signals, i.e., spectral transmittance signals or angular light-scattering signals, it is necessary to optimize the measurement signals.
Generally, to improve retrieval accuracy, the optical measurement signals are usually selected irrelevant to each other and sensitive to changes of inverse parameters. With this idea, the sensitivity analysis of optical measurement signals to unsolved characteristic parameters in ASDs at a different measurement wavelength or angle is employed to find the optimal measurement wavelength selection region for SEM and optimal measurement angle selection region for ALSM. The sensitivity coefficient is one of the most important characteristic parameters in the sensitivity analysis, which is the first derivative of the optical measurement signals, i.e., τ ( λ ) or I ( θ ) , to a certain inverse characteristic parameter. The sensitivity coefficient at every wavelength or angle χ a ( η ) is defined as [27]:
χ a ( η ) = ϕ ( a , η )   a = ϕ ( a + Δ a , η ) ϕ ( a Δ a , η ) 2 Δ a ,    ϕ = τ or I ,    η = λ or θ
where a denotes the inverse characteristic parameters in ASDs; Δ represents a tiny change and is set at 0.5% in this study, Δ λ , Δ θ > 0 ; I ( θ ) is the angular light-scattering measurement signals; τ ( λ ) is the spectral transmittance obtained in the SEM. It is obvious that to improve the retrieval accuracy, the optical measurement signals with larger absolute values of sensitivity coefficient should better be selected.
Figure 2 and Figure 3 depict the sensitivity coefficient for different monomodal ASDs with respect to spectral transmittance signals and angular light-scattering signals, respectively. From Figure 2, it can be found that increasing the measurement wavelength will reduce the sensitivity coefficient in general, so to ensure the retrieval accuracy, the measurement wavelengths in this study are selected from [ 0.2 , 1.4 ] μ m . From Figure 3, it is obvious that the value of the sensitivity coefficient will decline with an increase in the measurement angles. According to References [28,29], for a small angle, the scattered intensities are weakly dependent upon the particle’s composition but strongly dependent upon the particle’s roughness and size, so the measurement angles in this study are selected from [ 0 , 30 ] deg to obtain more accurate retrieval results. Similar conclusions can be obtained in studying bimodal ASDs, and the measurement wavelengths and angles are also selected from [ 0.2 , 1.4 ] μ m and [ 0 , 30 ] deg , respectively. The details in studying the bimodal ASDs are not shown here and interested readers can contact us.

5. Numerical Simulation

The true values of characteristic parameters of monomodal and bimodal ASDs in this study are listed in Table 1. According to our former studies [11], the number of measurement wavelengths and angles are set as 4 for monomodal ASDs and 7 for bimodal ASDs to make sure the retrieval accuracy, and the corresponding measurement wavelengths and angles are listed in Table 2.
The retrieval of the ASDs is solved by minimizing the objective function value, which is the sum of the square residuals between the estimated and measured transmittance:
F o b j   = i = 1 M { [ ϕ ( η ) ] est [ ϕ ( η ) ] mea [ ϕ ( η ) ] mea } 2 ,    ϕ = τ or I , η = λ or θ .
where [ ϕ ( a , η ) ] est and [ ϕ ( a , η ) ] mea denote the estimated and measured optical measurement signals, respectively. For the algorithm is a stochastic optimization method and all optimizations have certain randomness, all the calculations are repeated 30 times. Moreover, the relative deviation of ASDs δ is used to evaluate the quality of estimated results:
δ = { i = 1   N [ f est ( D ˜ i ) f true ( D ˜ i ) ] 2 } 1 / 2 { i = 1 N [ f true ( D ˜ i ) ] 2 } 1 / 2
where N denotes the number of subintervals which the particle size range [ D min , D max ] is divided into; D ˜ i is the midpoint of the ith subinterval [ D i , D i + 1 ]; f true ( D ˜ i ) is the true volume frequency distribution in the ith subinterval; and f est ( D ˜ i ) is the estimated volume frequency distribution in the ith subinterval.

5.1. Comparison of SPSO and Hybrid SPSO-DE Algorithms

At first, the performance of the hybrid SPSO-DE algorithm is investigated through comparison with the SPSO algorithm in retrieving the monomodal Gamma distribution by using an angular light-scattering method (See Figure 4). Table 3 lists the system control parameters. When the objective function value is less than 10−16 or the iteration number is larger than 1000, finish the calculation and output the results. The investigation shows that the convergence speed and the value of the objective function of the SPSO-DE hybrid algorithm are superior to the SPSO algorithm, which means the local optima and low convergence accuracy exiting in the SPSO can be avoided in the hybrid SPSO-DE algorithm.

5.2. Retrieval of the Aerosol Size Distribution

With the help of hybrid SPSO-DE algorithm, the ASDs mentioned above are reconstructed under different random measurement errors by using the SEM and ALSM, respectively, which are showed in Table 4 and Figure 5. Satisfactory retrieval accuracy can be obtained by both SEM and ASLM without measurement errors. The retrieval accuracy deteriorates with more random measurement errors, but are still reasonable even with 10% errors. Comparing the convergence accuracy of retrieval results by using the SEM with that by using the ALSM, it can be found that the robustness of retrieval results by using the ASLM is better than those by using the SEM. Moreover, in studying the bimodal ASDs, there is more satisfactory agreement between the estimated results and the true values of the major peak than those of the other peak.
Figure 6 and Figure 7 depict the distributions of the objective function values for monomodal ASDs by using two different type of optical measurement methods. Figure 8 shows the corresponding distributions of estimated results. From Figure 6 and Figure 7, it is obvious that the minimum value region of the objective function is a point, which means the optimal retrieval result is unique, especially for the monomodal Gamma distribution. Moreover, comparing the distributions of objective function by using the SEM and that by using the ALSM, it is obvious that the objective function value by using the SEM is smaller than that by using the ALSM. So, to obtain more accurate retrieval results by using the SEM, the lower objective function value should be obtained during inverse calculation, which will increase the convergence difficulty and time. This phenomenon can explain the reason for the better retrieval accuracies by using the ALSM than those by using the SEM. From Figure 8, it can be found that the distributions of the results retrieved by using the SEM show more dispersion than those retrieved by using ALSM, which also confirms the results obtained above.

6. Conclusions

Combined with a SPSO-DE hybrid algorithm, the convergence properties of angular light-scattering method (ALSM) and spectral extinction method (SEM) in retrieving ASDs are studied. Moreover, the sensitivity analysis of measurement signals to characteristic parameters in ASDs is also studied. Compared with the SEM, there are several advantages of using ALSM. The following conclusions can be drawn:
(1)
To obtain more accurate results, the measurement wavelengths should be better selected from the short wavelength region for SEM (e.g., λ [ 0.2 , 1.4 ] μ m in this study), and the measurement angles should better be selected from the forward angle region for ALSM (e.g., θ [ 0 , 30 ] deg in this study).
(2)
The SPSO-DE hybrid algorithm can converge much faster and obtain a lower fitness function value within a smaller number of generations than the SPSO algorithm in retrieving the ASDs.
(3)
The retrieval results by using ALSM show better convergence accuracy and robustness than those by using SEM, which attributes to the different distributions of objective function value for these methods. Moreover, considering only single spectral information of aerosol is required, the ALSM provides a more effective and reliable technique to obtain spherical ASDs. Future research will focus on the application of ALSM in studying spheroidal ASDs.

Author Contributions

All authors have contributed equally to this work. Z.-Z.H., J.-K.M. and X.-S.H. conceived the structure of paper; Z.-Z.H. and J.-K.M. wrote the simulation program; Z.-Z.H. and X.-S.H. analyzed the data; Z.-Z.H., J.-K.M. and X.-S.H. wrote the paper.

Funding

The supports of this work by National Natural Science Foundation of China (No: 51806103 and No: 51606095) and Jiangsu Provincial Natural Science Foundation (No: BK20170800, BK20160794) are gratefully acknowledged.

Acknowledgments

A very special acknowledgement is made to the editors and referees who make important comments to improve this paper. Moreover, the authors are grateful to Deepak Adarsh and D’Almeida Guillaume A. for supplying the spectral optical constants of aerosols.

Conflicts of Interest

The authors declare no conflict of interest.

Notation

The following notations are used in this manuscript:
C R crossover constant
D particle diameter, μm
f ( D ) unknown volume frequency distribution
F obj objective function value
F mutant factor in the DE algorithm
I ( θ ) angular light-scattering intensity at the angle θ , W / ( m 2 sr )
I 0 total intensity of the laser, W / ( m 2 sr )
i , i 1 , i 2 Mie scattering functions
k absorption index of optical constant
m λ optical constant of particle at wavelength λ
M number of the measurement wavelengths or angles
n refractive index of optical constant
N 0 total number concentration of the suspended particle system
N c maximum number of iterations
N s total number of the particles in the swarm
N p number of the inversion parameters (dimensions)
P i , P g , P m the individual, global and mean best positions
Q ext extinction efficiency
V i velocity of the i th particle
X i position of the i th particle

Greeks Symbols

λ incident wavelength of the laser, μm
χ sensitivity coefficient
θ scattering light measurement angle
δ relative deviation of the aerosol size distribution
τ transmittance of the particle system

Abbreviations

The following abbreviations are used in this manuscript:
estestimated value
L-NLog-Normal distribution
GammaGamma distribution
maxmaximum value
minminimum value
truetrue value

References

  1. Chen, Q.X.; Yuan, Y.; Shuai, Y.; Tan, H.P. Graphical aerosol classification method using aerosol relative optical depth. Atmos. Environ. 2016, 135, 84–91. [Google Scholar] [CrossRef]
  2. Gong, W.; Zhang, S.; Ma, Y. Aerosol Optical Properties and Determination of Aerosol Size Distribution in Wuhan, China. Atmosphere 2014, 5, 81–91. [Google Scholar] [CrossRef] [Green Version]
  3. Esparza, A.E.; Fitzgerald, R.M.; Gill, T.E.; Polanco, J. Use of light-extinction method and inverse modeling to study aerosols in the Paso del Norte Airshed. Atmos. Environ. 2011, 45, 7360–7369. [Google Scholar] [CrossRef]
  4. Ladji, R.; Yassaa, N.; Balducci, C.; Cecinato, A. Particle size distribution of n-alkanes and polycyclic aromatic hydrocarbons (PAHS) in urban and industrial aerosol of Algiers, Algeria. Environ. Sci. Pollut. Res. 2014, 21, 1819–1832. [Google Scholar] [CrossRef] [PubMed]
  5. Kokhanovsky, A.A.; Leeuw, G.H. Satellite Aerosol Remote Sensing over Land; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  6. Chen, Q.X.; Yuan, Y.; Huang, X.; Jiang, Y.Q.; Tan, H.P. Estimation of surface-level PM 2.5 concentration using aerosol optical thickness through aerosol type analysis method. Atmos. Environ. 2017, 159, 26–33. [Google Scholar] [CrossRef]
  7. Aerosol Robotic Network (AERONET). Available online: http://aeronet.gsfc.nasa.gov/ (accessed on 6 July 2018).
  8. Moderate Resolution Imaging Spectroradiometer (MODIS). Available online: https://modis.gsfc.nasa.gov/ (accessed on 6 July 2018).
  9. Wang-Li, L.; Cao, Z.; Buser, M.; Whitelock, D.; Parnell, C.B.; Zhang, Y. Techniques for measuring particle size distribution of particulate matter emitted from animal feeding operations. Atmos. Environ. 2013, 66, 25–32. [Google Scholar] [CrossRef]
  10. He, Z.Z.; Mao, J.K.; Han, X.S. Non-parametric estimation of particle size distribution from spectral extinction data with PCA approach. Powder Technol. 2018, 325, 510–518. [Google Scholar] [CrossRef]
  11. He, Z.Z.; Qi, H.; Yao, Y.C.; Ruan, L.M. Inverse estimation of the particle size distribution using the Fruit Fly Optimization Algorithm. Appl. Therm. Eng. 2015, 88, 306–314. [Google Scholar] [CrossRef]
  12. Ye, M.; Wang, S.; Lu, Y.; Hu, T.; Zhu, Z.; Xu, Y. Inversion of particle-size distribution from angular light-scattering data with genetic algorithms. Appl. Opt. 1999, 38, 2677–2685. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Wang, Y.; Liang, G.; Pan, Z. Inversion of particle size distribution from light-scattering data using a modified regularization algorithm. Particuology 2010, 8, 365–371. [Google Scholar] [CrossRef]
  14. He, Z.Z.; Hong, Q.; Chen, Q.; Ruan, L.M. Retrieval of aerosol size distribution using Improved Quantum-behaved Particle Swarm Optimization from spectral extinction measurement. Particuology 2014, 133, 245–252. [Google Scholar]
  15. Dombrovsky, L.A.; Baillis, D. Thermal Radiation in Disperse Systems: An Engineering Approach; Begell House: New York, NY, USA, 2010. [Google Scholar]
  16. Vargas-Ubera, J.; Aguilar, J.F.; Gale, D.M. Reconstruction of particle-size distributions from light-scattering patterns using three inversion methods. Appl. Opt. 2007, 46, 124–132. [Google Scholar] [CrossRef] [PubMed]
  17. Kennedy, J.; Eberhart, R. Particle swarm optimization. In Proceedings of the IEEE International Conference on Neural Networks, Perth, Australia, 27 November–1 December 1995; pp. 1942–1948. [Google Scholar]
  18. Qi, H.; Ruan, L.M.; Zhang, H.; Wang, Y.; Tan, H.P. Inverse radiation analysis of a one-dimensional participating slab by stochastic particle swarm optimizer algorithm. Int. J. Therm. Sci. 2007, 46, 649–661. [Google Scholar] [CrossRef]
  19. Yuan, Y.; Yi, H.L.; Shuai, Y.; Liu, B.; Tan, H.P. Inverse problem for aerosol particle size distribution using SPSO associated with multi-lognormal distribution model. Atmos. Environ. 2011, 45, 4892–4897. [Google Scholar] [CrossRef]
  20. Wei, L.Y.; Qi, H.; Ren, Y.T.; Ruan, L.M. Application of stochastic particle swarm optimization algorithm to determine the graded refractive index distribution in participating media. Infrared Phys. Technol. 2016, 79, 74–84. [Google Scholar] [CrossRef]
  21. Wu, Y.; Lu, J.; Sun, Y. An improved differential evolution for optimization of chemical process. Chin. J. Chem. Eng. 2008, 16, 228–234. [Google Scholar] [CrossRef]
  22. Lu, S.; Sun, C.; Lu, Z. An improved quantum-behaved particle swarm optimization method for short-term combined economic emission hydrothermal scheduling. Energy Convers. Manag. 2010, 51, 561–571. [Google Scholar] [CrossRef]
  23. Storn, R.; Price, K. Differential evolution-a simple and efficient heuristic for global optimization over continuous spaces. J. Glob. Optim. 1997, 11, 341–359. [Google Scholar] [CrossRef]
  24. Yuan, Y.; Yi, H.L.; Shuai, Y.; Wang, F.Q.; Tan, H.P. Inverse problem for particle size distributions of atmospheric aerosols using stochastic particle swarm optimization. J. Quant. Spectrosc. Radiat. Transf. 2010, 111, 2106–2114. [Google Scholar] [CrossRef]
  25. Deepak, A.; Gerber, H.E. Report of the Experts Meeting on Aerosols and Their Climatic Effects, Williamsburg, VA, 28–30 March 1983; World Meteorological Organization: Geneva, Switzerland, 1983. [Google Scholar]
  26. D’Almeida, G.A.; Koepke, P.; Shettle, E.P. Atmospheric Aerosols: Global Climatology and Radiative Characteristics; A Deepak Pub.: Hampton, VA, USA, 1991. [Google Scholar]
  27. Kamrunnahar, M.; Braatz, R.D.; Alkire, R.C. Parameter Sensitivity Analysis of Pit Initiation at Single Sulfide Inclusions in Stainless Steel. J. Electrochem. Soc. 2004, 151, 90–97. [Google Scholar] [CrossRef]
  28. Renard, J.B.; Dulac, F.; Berthet, G.; Lurton, T.; Vignelles, D.; Jégou, F.; Tonnelier, T.; Jeannot, M.; Couté, B.; Akiki, R. LOAC: A small aerosol optical counter/sizer for ground-based and balloon measurements of the size distribution and nature of atmospheric particles—Part 1: Principle of measurements and instrument evaluation. Atmos. Meas. Tech. 2016, 9, 1721–1742. [Google Scholar] [CrossRef] [Green Version]
  29. Lurton, T.; Renard, J.B.; Vignelles, D.; Jeannot, M.; Akiki, R.; Mineau, J.L.; Tonnelier, T. Light scattering at small angles by atmospheric irregular particles: Modelling and laboratory measurements. Atmos. Meas. Tech. 2014, 7, 931–939. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Spectral optical constants of four different kind of aerosols [25,26].
Figure 1. Spectral optical constants of four different kind of aerosols [25,26].
Computation 06 00047 g001
Figure 2. Sensitivity analysis of spectral transmittance signals to characteristic parameters in ASDs: (a,c) for L-N distribution, (b,d) for Gamma distribution.
Figure 2. Sensitivity analysis of spectral transmittance signals to characteristic parameters in ASDs: (a,c) for L-N distribution, (b,d) for Gamma distribution.
Computation 06 00047 g002aComputation 06 00047 g002b
Figure 3. Sensitivity analysis of angular light-scattering signals to characteristic parameters in ASDs: (a,c) for L-N distribution, (b,d) for Gamma distribution.
Figure 3. Sensitivity analysis of angular light-scattering signals to characteristic parameters in ASDs: (a,c) for L-N distribution, (b,d) for Gamma distribution.
Computation 06 00047 g003
Figure 4. Comparison of objective function values of SPSO and hybrid SPSO-DE algorithms in retrieving monomodal Gamma distribution.
Figure 4. Comparison of objective function values of SPSO and hybrid SPSO-DE algorithms in retrieving monomodal Gamma distribution.
Computation 06 00047 g004
Figure 5. Reproducibility of different ASDs by using different optical measurement methods: (a) monomodal L-N distribution; (b) monomodal Gamma distribution; (c) bimodal L-N distribution; (d) bimodal Gamma distribution.
Figure 5. Reproducibility of different ASDs by using different optical measurement methods: (a) monomodal L-N distribution; (b) monomodal Gamma distribution; (c) bimodal L-N distribution; (d) bimodal Gamma distribution.
Computation 06 00047 g005
Figure 6. Distribution of the objective function values for different ASDs by using the SEM: (a) L-N distribution; (b) Gamma distribution.
Figure 6. Distribution of the objective function values for different ASDs by using the SEM: (a) L-N distribution; (b) Gamma distribution.
Computation 06 00047 g006
Figure 7. Distribution of the objective function values for different ASDs by using the ALSM: (a) L-N distribution; (b) Gamma distribution.
Figure 7. Distribution of the objective function values for different ASDs by using the ALSM: (a) L-N distribution; (b) Gamma distribution.
Computation 06 00047 g007
Figure 8. Distribution of retrieval results of monomodal ASDs for 30 calculations by using the SEM and ALSM: (a) and (b) for SEM; (c) and (d) for ALSM.
Figure 8. Distribution of retrieval results of monomodal ASDs for 30 calculations by using the SEM and ALSM: (a) and (b) for SEM; (c) and (d) for ALSM.
Computation 06 00047 g008aComputation 06 00047 g008b
Table 1. The true values of characteristic parameters of ASDs.
Table 1. The true values of characteristic parameters of ASDs.
Function TypesParametersTrue Values
MonomodalL-N ( D ¯ , σ ) ( 1.0 , 2.99 )
Gamma ( α , β ) ( 6.0 , 3.0 )
BimodalL-N ( D ¯ 1 , σ 1 , D ¯ 2 , σ 2 , n ) ( 2.5 , 1.1 , 6.0 , 1.2 , 0.2 )
Gamma ( α 1 , β 1 , α 2 , β 2 , n ) ( 6.0 , 2.2 , 0.8 , 2.0 , 0.15 )
Table 2. The measurement wavelengths and angles for retrieval of ASDs.
Table 2. The measurement wavelengths and angles for retrieval of ASDs.
Function TypesSpectral Extinction Method, λ (μm)Angular Light-Scattering Method, θ (deg)
MonomodalL-N ( 0.2 , 0.4 , 0.6 , 0.8 ) ( 5 , 10 , 15 , 20 )
Gamma ( 0.2 , 0.4 , 0.6 , 0.8 ) ( 5 , 10 , 15 , 20 )
BimodalL-N ( 0.2 , 0.4 , 0.6 , 0.8 , 1.0 , 1.2 , 1.4 ) ( 4 , 8 , 12 , 16 , 20 , 24 , 28 )
Gamma ( 0.2 , 0.4 , 0.6 , 0.8 , 1.0 , 1.2 , 1.4 ) ( 4 , 8 , 12 , 16 , 20 , 24 , 28 )
Table 3. System control parameters of the SPSO and hybrid SPSO-DE algorithms in retrieving monomodal Gamma distribution.
Table 3. System control parameters of the SPSO and hybrid SPSO-DE algorithms in retrieving monomodal Gamma distribution.
ParametersNsNpNc ε C1C2FCR
SPSO502100010−161.01.0
SPSO-DE502100010−161.01.00.50.4
Table 4. Retrieval results of ASDs under different optical measurement methods.
Table 4. Retrieval results of ASDs under different optical measurement methods.
Different FunctionsRandom ErrorSEMASLM
Retrieval Error, δRetrieval Error, δ
MonomodalL-N0%0.0001250.000312
5%0.1304010.064109
10%0.2151680.162801
20%0.3876210.265412
Gamma0%0.0003220.000545
5%0.1627210.144108
10%0.1875420.169336
20%0.3657290.287532
BimodalL-N0%0.0001560.003148
2%0.0965420.085159
5%0.2444940.232532
10%0.3562790.328765
Gamma0%0.0252420.090019
2%0.1456230.110330
5%0.1815290.144053
10%0.3197320.298621

Share and Cite

MDPI and ACS Style

He, Z.-Z.; Mao, J.-K.; Han, X.-S. Determination of Aerosol Size Distribution from Angular Light-Scattering Signals by Using a SPSO-DE Hybrid Algorithm. Computation 2018, 6, 47. https://0-doi-org.brum.beds.ac.uk/10.3390/computation6030047

AMA Style

He Z-Z, Mao J-K, Han X-S. Determination of Aerosol Size Distribution from Angular Light-Scattering Signals by Using a SPSO-DE Hybrid Algorithm. Computation. 2018; 6(3):47. https://0-doi-org.brum.beds.ac.uk/10.3390/computation6030047

Chicago/Turabian Style

He, Zhen-Zong, Jun-Kui Mao, and Xing-Si Han. 2018. "Determination of Aerosol Size Distribution from Angular Light-Scattering Signals by Using a SPSO-DE Hybrid Algorithm" Computation 6, no. 3: 47. https://0-doi-org.brum.beds.ac.uk/10.3390/computation6030047

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop