Study of Particle Size Measurement by the Extinction Method in Flame

: The laser extinction method (LEM) is particularly suitable for measuring particle sizes in ﬂames because this method, which is based on the Beer–Lambert law, is non-intrusive and easy to implement. In the LEM, the interpretation of the extinction data is usually developed under the assumption that light extinction due to scattering is a result of the superposition of single scattering by individual particles; however, this could be violated for ﬂames with dense concentrations of particles in which multiple scattering could occur. Quantifying the effect of multiple scattering under general conditions is still a formidable problem. In this work, we carried out a series of careful measurements of the laser extinction using standard particles of various known sizes, number densities and optical path lengths, all under the condition that the acceptance angle of the detector was limited to nearly zero. Combined with a four-ﬂux model, we quantitatively analyzed the effect of multiple scattering on the size measurement using the LEM. The results show that the effect of multiple scattering could be ignored when the optical thickness is less than two under strict restrictions on the detector acceptance angle. Guided by this, the size distribution of an alumina (Al 2 O 3 ) particle sample was measured by the LEM with dual wavelengths. Parameterized distributions were solved with the help of graph plotting, and the results compared well with the measurement from the Malvern particle size analyzer. The same method was then used to measure the particle size distribution in the plume of a solid rocket motor (SRM). The use of an off-axis parabolic mirror in the experimental setup could suppress the jitter of light passing through the SRM plume, and the particle size in the plume of the measured SRM was in the order of microns.


Introduction
In industrial processes, many combustion processes include particle combustion, and the factors affecting the combustion efficiency are not only related to the characteristics of the fuel itself but also to the flow state and size of the particles.For example, the combustion of fuel in an internal combustion engine is the combustion of fuel particles after atomization, and the size of the fuel particles after atomization affects the combustion efficiency and exhaust emissions [1][2][3].For instance, coal combustion in a boiler is the combustion process of pulverized coal in the form of granules, and the morphology and kinematic state of the particles change dramatically during the combustion process, while the morphology, size and kinematic parameters of the fuel particles have a significant impact on the performance of the burner [4,5].The kinematic state and size of the fuel particles are a prerequisite for the efficient organization of the mass flow and are important conditions for the optimization of the burner.
In the aerospace field, rockets and certain spacecraft mostly use composite-modified double base (CMDB) propellants with specific properties or double base (DB) propellants as fuel, and metal particles such as lithium, magnesium and aluminum are generally added to the propellant to improve the specific impulse of the propellant; non-metallic particles of boron can also be added [6,7].In addition to increasing the specific impulse of a propellant, the particles contain propellant that can also suppress the oscillatory combustion of the propellant [8,9]; though this principle is complex and has not been fully understood, there is a consensus that the effect of metal particles suppressing the oscillatory combustion of the propellant can be related to the size of the metal particles after combustion.Therefore, the size of the metal particles after combustion is a parameter that needs to be considered when designing solid rocket motors (SRMs).However, the addition of metal particles to propellants can lead to some negative effects, which include two-phase flow losses in the SRMs [10][11][12] and nozzle erosion caused by added particles [13][14][15].Some solid rocket motors work in outer space, where the large particles or debris produced by these SRMs are a kind of space debris, and these particles pose hazards to space transportation [16].Large solid particles pose a threat to equipment that operates in space, where particles larger than 100 µm in diameter can penetrate astronauts' space suits and particles larger than 5 mm in diameter can penetrate the weak parts of the space station [17].To date, there are relatively few internal and external flow field parameters that can be measured for SRMs.The curve of the motor thrust and front-end pressure in SRM over time can only almost be measured [18,19].Only a few sets of experimental data in the published literature can be used to verify the theoretical analysis and numerical calculation of the two-phase flow in solid rocket motors, and the measurement of particle velocity and particle size in SRMs' plume flames can provide empirical support for the theory of a new energy loss model [20].For these reasons, there has been wide demand for the online measurement of particle size in the plume flames of SRMs.
To measure the size of the particles in flames, optical measurement technologies should be prioritized due to their non-intrusive and non-contact advantages [21,22]; however, combusting environments bring measurements great challenges [22] that include where the number density of the measured particles is uncontrollable and usually very high; on the other hand, light passing through the flame can be deflected in an uncertain direction due to the variation in the gas density when raised from the temperature gradient of the field, especially for large-scale measurement objects.An example of a high particle number density leading to failure in the particle size measurement is the failure of PDPA in dense diesel fuel [23].As a result of the above challenges, many studies on particle size in flames use sampling methods [16,22,24].However, among the existing online particle size measurement technologies, the light extinction method (LEM) may be the most suitable approach for measuring the particle size of dense particles in combusting environments.The advantages of this technique over other existing online particle size measurement techniques such as angular scattering method, diffraction method, PDPA, etc., include: (1) The technique can be used for dense particle size measurement.( 2) By means of a special experimental setup design, the technique can suppress the jitter of the laser passing through the flame.(3) The technique is simple to set up and easy to measure, requiring only the data of the laser transmittance [25].(4) Compared to particle size measurement techniques that require the scattered signal of the particles to be measured, usually the intensity of the transmitted laser is higher than that of the scattered laser, making it easier to measure and eliminating the need for an expensive detector at the receiving end.The measurement principle of the LEM is the Beer-Lambert law [26], which is based on the hypothesis that the scattered light between the measured particles satisfies the single scattering condition; that is, the detector does not receive any scattered light from particles in the measurement volume.In other words, the multiple scattering effect is completely ignored.In earlier examples in the literature, researchers used the LEM to measure the particle size without discussing the multiple scattering effects [27,28].The fact that many parameters can affect the multiple scattering effect complicates this problem, and the random distribution of particle positions can lead to difficulties in determining not only the intensity of the scattered light but also the phase of the electromagnetic field.These factors make it almost impossible to theoretically determine the light intensity distribution in the case of multiple scattering.There is also no universal conclusion as to when multiple scattering occurs and when the effects of multiple scattering can be ignored for the LEM.
In 1969, Kerker [29] showed that the effect of multiple scattering could be ignored when the transmission of photons through a sample of particles was greater than 90%.In 1981, Bayvel and Jones [30] showed that the effect of backscattering could be ignored if the transmittance was greater than 60%.In 2007, the scattering of eight particles with a wavelength size and randomly assigned positions was calculated by Mishchenko et al. [31] using the T-matrix method, and the results showed that the single-scattering approximation (SSA) was satisfied only when the distance between particles was large and the number density of particles was low.When the volume fraction of particles was less than 1% and x × l pp < 30, where x represents the wave number of the particle size in the surrounding medium and l pp is the average distance between two particles, the effects of multiple scattering were negligible.Su Mingxu et al. [32] used the Monte Carlo method to simulate different light wavelengths, particle sizes and different detection sizes and detector distances and compared the calculation results with those calculated using the Beer-Lambert law.The results showed that the influence of the multiple scattering effects on the particle size measurement using the LEM not only depended on the particle concentration but also on the size of the detector and the distance of the detector's placement, which had a great influence on the measurement results.Reducing the particle concentration and the detector receiving area and increasing the detector distance could diminish the effect of multiple scattering on the measurement results.Swanson et al. [33] showed that by limiting the half-angle of acceptance θ 1/2 in the detector, the effect of multiple scattering could be avoided up to an optical thickness of 10, where θ 1/2 satisfied: where λ represents the laser wavelength in the medium and d p is the particle diameter.However, this conclusion was based on only the limited transmittance data of standard polystyrene particles and did not take into account the sedimentation of polystyrene during the experiment.Determining the conditions under which the particle multiple scattering effects can be ignored is a key point in the particle size experiment using an LEM.In this paper, this condition was investigated by measuring the transmittance of the polystyrene aqueous solution with the half-angle of acceptance θ 1/2 in the detector, which satisfied Equation (1).To prevent the sedimentation of polystyrene particles from interfering with the measurement results, a magnetic stirrer was used to continuously stir the polystyrene solution.The particle size distribution of alumina was then measured using an LEM with two laser wavelengths and was compared with the result measured using the Malvern laser particle size analyzer.Finally, the size distribution of the particles in the SRM plume flame was measured.For this experiment, the problem to be solved was that the laser passing through the flame might not be received by the detector due to refraction caused by the density differences in the flame field.

Materials and Methods
The LEM's measurement principle and two sets of experimental setups were used: one to determine the conditions for neglecting multiple scattering and another to measure the size distribution of alumina powder and particles in the SRM plume flame, which are described in this section.

Measurement Principle of the LEM (Beer-Lambert Law)
When a laser beam with light intensity I 0 irradiates a continuous homogeneous medium, the outgoing light intensity I is attenuated due to the absorption and scattering of the light by the medium.Assuming that the light intensity attenuation of the laser proceeding through a medium with a thickness of dl is −dI, the relationship between the two is thought to be: where τ represents turbidity.When Equation ( 2) is integrated along the direction of light propagation and boundary conditions are used, we obtain: where ξ = τl is the optical thickness and l is the optical length.When the medium is a cloud of particles, and the multiple scattering effects between particles is ignored, τ can be expressed as: where N p is the number density of the sample particles, D p is the particle diameter and Q ext is the extinction coefficient.For polydisperse particles, τ can also be easily derived.
Q ext is the function of particle size D p , laser wavelength λ and complex refractive index of m, which can be calculated according to the Mie theory [34].Generally Q ext and are expressed as: where where m p is the refractive index of the particle, m 0 is the refractive index of the medium around the particle and k is a parameter that is related to the absorption coefficient of light.

Experimental Setup and Experimental Methods
Two experiments are presented here.The first experiment was designed to determine under what conditions the effects of multiple scattering could be ignored when measuring the particle size using the LEM, and the second was designed to measure the size distribution of the particles in flame.

Experimental Setup Used to Study the Multiple Scattering Effect
The transmittance of the laser passing through standard polystyrene particles and suspended in water with different optical thicknesses was measured.A schematic diagram of the experimental setup is shown in Figure 1a, and some photos of the experimental setup are shown in Figure 1b.The laser (generated by SUMLINEG-800), with a wavelength of 532 nm and a laser beam diameter of 1 mm, was split into two beams after passing through a beamsplitter, whose reflectance-to-transmittance ratio, T:R, was 1:9.The transmitted light was received by detector 1 (Thorlabs, PDA36A2), whose function was to monitor the actual output power of the laser in real time, and a bandpass filter with a center wavelength of 532 nm was placed in front of the detector.The reflected laser passed through an acrylic rectangular container filled with purified water, and a triangular sample box to which polystyrene particles were added was placed inside.The triangular box, the main dimensions of which are shown in Figure 1c, could be moved back and forth in the rectangular container so that the optical length l of the laser in the triangular box could be easily changed.A magnetic stirrer was placed in the triangular box, and the sedimentation of polystyrene particles can be prevented by continuous stirring.The laser light, after passing through the polystyrene solution, was received by detector 2; in addition, a bandpass filter with a center wavelength of 532 nm and an aperture of 3 mm diameter was placed in front of detector 2. The diameter of the polystyrene that was used and the optical length l are listed in Table 1.The distance between detector 2 and the triangular box was about 1.5 m; therefore, Equation (1) could be satisfied.The light intensity ratio α l between the two detectors for every optical length l needed to be calibrated in advance before the particles were added in the triangular box, α l , as follows: where V 1 and V 2 are the value of detector 1 and detector 2 when no particles are added in the triangular box, respectively; V d1 and V d1 , close to 10 mV, are the dark voltage of these two detectors; and the symbol '-' denotes the averaging operation.
When polystyrene particles of a known mass were added to the triangular box, the triangular box was moved to the locations of different optical lengths l, and the voltages of the two detectors, V l1 and V l2 , were recorded for different optical lengths l; then, the laser transmittance T tot of the particles was:

Experimental Setup Used to Measure the Size Distribution of Particles in Flame
The experimental setup for measuring the size distribution of the particles in the SRM plume flame using the LEM is shown in Figure 2.This experimental setup mainly included two parts: the transmitter and the receiver.Figure 2a is the schematic diagram of the experimental setup, and Figure 2b,c show the photos of the transmitter and receiver, respectively.In the transmitter, the dual-wavelength laser-emitted laser beams had wavelengths of 532 nm and 1064 nm simultaneously.The two lasers' wavelengths were separated after passing through the beam splitter, which reflected the 532 nm laser and transmitted the 1064 nm laser.The two separated lasers again passed through each beamsplitter, whose ratio of transmission to reflection, T:R, was approximately 1:9.The transmitted laser light was received by detectors 1 and 2, in front of which the corresponding bandpass filters were placed.Both the 532 nm and 1064 nm reflected lasers became collinear lasers again after passing through the beamsplitter.
As mentioned earlier, one of the problems that needed to be solved when measuring the size of the particles in the SRM plume flame was the jittering of the laser after passing through the flame due to the irregular laser deflection caused by the density difference of the flame.Without special treatment, the detectors in the receiver may not have been able to receive the transmitted light.To solve this problem, two off-axis parabolic mirrors were used in the receiver because the jittering laser reflected by an off-axis parabolic mirror always passed through the focal point of the mirror where the detectors were placed.The off-axis parabolic mirror had a cross-sectional diameter of 50.8 mm (2 inches) and a reflective focal length of 50.8 mm (2 inches).The main optical path at the receiver was as follows: the light from the dual-wavelength laser passing through the flame was split into two beams, a 532 nm one and a 1064 nm one, which were reflected by the off-axis parabolic mirror and received by detector 3 and detector 4, respectively.A bandpass filter and an aperture of a 3 mm diameter were placed in front of the detector.In this experiment, although bandpass filters were placed in front of each detector in the receiver to ensure that the detectors were as unaffected as possible by the flame's own radiation, the laser was operated in pulse mode at 10 Hz for both wavelengths with a duty cycle of 50%.In the pulsed laser mode, if the detector value was I on when there was a laser output and I o f f when there was no laser output, the effect of the flame radiation on the measurement results could be eliminated by subtracting I o f f from I on .
In order to verify the accuracy of the particle size measurement taken using a dualwavelength extinction method, the particle size of the alumina powder was measured using this experimental setup, and the measurement results were compared with the results measured by the Malvern laser particle size analyzer (MASTERSIZER 2000).The alumina powder needed to be dispersed when measured by the dual-wavelength extinction method; however, it could not be suspended in water because the density of alumina is about 4 times that of the water.Therefore, the alumina particles were cured in an epoxy resin.When the alumina particles were in the epoxy, the absorption effect of the epoxy on the laser had to be taken into account.To calculate the absorption of the laser light by the epoxy resin, an epoxy resin sample that did not contain particles was also made; the epoxy resin samples are shown in Figure 3.The sample in Figure 3a contained no aluminum oxide particles, whereas the sample in Figure 3b contained aluminum oxide particles.The two samples were identical in size, and the sample without alumina particles, shown in Figure 3a, was used to calibrate the absorption of the laser light by the epoxy resin; i.e., if a laser with an intensity of I 0 changed to I 1 after passing through the sample without alumina and to I 2 after passing through the sample containing alumina, the intensity ratio of the light caused by the alumina particles in the sample would be I 2 I 1 .When measuring the particle size in the SRM plume flame, the measurement volume was 0.9 m from the nozzle outlet of the SRM.The receiver was placed approximately 2.2 m from the flame to ensure the measurement condition satisfied Equation (1).

Results and Discussion
In this section, the conditions under which the multiple scattering effects could be neglected is discussed, and the alumina particle sizes measured by the dual-wavelength extinction method and the Malvern laser particle size analyzer are compared.Finally, the results of the particle size distribution in the plume of a full SRM are presented.

Conditions for Neglecting the Multiple Scattering Effect
A four-flux model [35,36] was used here to study multiple scattering effects.In this model, the light in the sample could be classified into four fluxes: 1.
A collimated beam of intensity I c propagating to negative e z ; 2.
A diffuse radiation of intensity I d propagating to negative e z ; 3.
A collimated beam of intensity J c propagating to positive e z ; 4.
A diffuse radiation of intensity J d propagating to positive e z .
The negative e z is the direction of incidence in the laser, and the relationship between these four fluxes is as follows: where where N is the particle number density, R is the radius of the particle, Q abs is the absorption coefficient, Q sca is the scattering coefficient calculated using Mie theory, ε is the average crossing parameter, and ζ is the forward scattering coefficient, e.g., where p(θ) is the scattered light intensity distribution function of θ, which is also calculated using Mie theory.According to Equation (15) and Mie theory, the values of ζ for the polystyrene particles with diameters of 1, 2 and 3 µm in the water were 0.985, 0.942 and 0.955, respectively.In this paper, ε was taken as a constant value of two, just as in the paper [35].Moreover, k a ≈ 0 because the k of polystyrene in Equation ( 6) is approximately equal to zero [37].According to the four-flux model, the transmittance after the laser passed through the sample could be classified into two parts, collimate-collimate transmittance T cc and collimate-diffuse transmittance T cd due to the absence of diffuse-diffuse incident light.T cc and T cd were calculated as follows, depending on the boundary conditions: In the four-flux model, the area of the detector in the receiver is assumed to be infinitely large; however, in practice, the area of the detector could not be infinitely large.Instead, to avoid the effects of multiple scattering, it needed to be as small as possible.T cc was independent of the detector area but T cd depended on it.For the experimental measurements, the transmittance T tot needed to be: where η is the parameter that should depend on the geometry of the sample box; however, this would be difficult to analyze theoretically for a complicated structure.The value of η was determined by fitting the polystyrene data with a diameter of 1 µm using Equation ( 18  As can be seen in Figure 4, the experimental values of laser transmittance gradually deviated from the Beer-Lambert law as the optical thickness ξ increased, even when the detector reception angle satisfied Equation (1).When the optical thickness was ξ > 2, the difference between the experimental values of the transmittance and theoretical values based on the Beer-Lambert law started to be noticeable.This meant that the effect of multiple scattering needed to be taken into account when using the LEM for particle size measurements.However, in a similar study, Swanson et al. [33] argued that the applicability of the LEM using the Beer-Lambert theorem was ξ < 10.A possible reason for this discrepancy was that Swanson et al. did not take into account the sedimentation of polystyrene particles in the water during the experiment, where sedimentation could result in a higher particle number density at the measurement point compared to a theoretical one.The higher particle number density resulted in a lower measured transmittance compared to the theoretical value, while in the region of greater optical thickness ξ the presence of the multiple scattering effects resulted in a higher transmittance, and the combined effect of these two opposing effects on the transmittance ultimately resulted in the measured transmittance not being significantly different from the Beer-Lambert law.
When the optical thickness was ξ < 4, the four-flux model agreed relatively well with the experimental transmittance values.However, as the optical thickness ξ increased, the model's predictions began to deviate from the experimental values.In order to make the four-flux model fit the experimental values better, a correction to ε could be considered; in addition, one could set η as a function of T cd , as was conducted in the reference [35], but this would complicate the problem further.

Alumina Particle Size Distribution Measurement Results
In this section, the size distribution of alumina particles, which were cured in an epoxy resin because they could not be suspended in water, are measured using the dual wavelengths, 532 nm and 1064 nm, while the extinction method and its results are compared with those measured by the Malvern laser particle size analyzer.
To invert the particle size distribution using the LEM, the extinction coefficient Q ext needed to be calculated.The complex refractive indices of the alumina for the laser lengths of 532 nm and 1064 nm ( m λ532 and m λ1064 , respectively) for this calculation using Mie theory were [38,39]: The extinction coefficient Q ext of alumina particles at two wavelengths of incident light versus the particle size is shown in Figure 5, which indicates that Q ext was different for different incident wavelengths and even for the same particle size.Provided that the measured particle size was polydispersed, the volume distribution functions of the particle sizes are represented as f V (D p ), with the total mass of the sample denoted as m total , the particle density denoted as ρ p , the total volume occupied by the dispersed particles in the sample denoted as V total , and the optical length denoted as l; then, the transmittance of the laser T was given by where m v means the mass of particles per unit of dispersed volume, which could be defined as follows: The value m v of the sample in Figure 3b was 0.016 mg/mL.
When the transmittance of the sample was measured at two wavelengths, the equations for f V (D p ) are as follows: where Q ext1064 and Q ext532 are the extinction coefficients for these two wavelengths, as shown in Figure 5. T λ 1064 and T λ 532 are the transmittances for the 532 nm and 1064 nm lasers, respectively.Here, the values of these two were 0.690 ± 0.0085 and 0.611 ± 0.0056, and the corresponding optical thicknesses ξ were 0.37 and 0.49; as these two optical thickness ξ were both less than two, the multiple scattering effects could be ignored when measuring the particle size using the LEM according to the conclusions in Section 3.1.
In theory, f V (D p ) could be solved by Equation ( 23) if f V (D p ) is a function of two parameters.The extinction coefficient Q ext is expressed in terms of the infinite series of Bessel functions, which corresponds with the integral term in Equation (23); for example, dD p , which cannot be expressed in the primitive function form.This means that Equation (23) could not be solved directly.Here, we used the graph-plotting method to solve the system of Equation ( 23) through the discretization method, which converts the integral terms above into a sum containing the particle size distribution function parameters, i.e., where P 1 and P 2 are the two parameters of the two-parameter particle size distribution; λ is the incident light wavelength; D min and D max are the minimum and maximum particle sizes in the measured particle population, which can be estimated in advance; ∆D i is the discrete step size; and N is the number of discrete steps: After discretization, the left terms of Equation ( 23) became two functions of the particle size distribution parameter (P 1 , P 1 ), i.e., The surfaces of the functions in Equations ( 26) and ( 27) could be plotted with respect to the particle size distribution parameters; then, the intersection of these two surfaces with the plane of "Z = 0" was the solution of Equation (23).
In this paper, we selected the commonly used two-parameter particle size distributions, Rosin-Rammler (R-R) and Logarithmic-Normal (L-N) distribution, as seen in Equations ( 28) and ( 29), to describe the particle size distributions of the measured alumina powder samples.
The R-R distribution was stated as: where the parameter k describes the distribution uniformity and is hence called the uniformity constant and D is the "characteristic particle size", defined as the size at which 63.2% of the particles (by volume) are smaller.The L-N distribution could be stated as: where σ p and ln D are the mean and standard deviation of the logarithm of the particle diameter, respectively.The solutions of f V (D) for the R-R and L-N distributions in Equation ( 23) for the mean transmittance at two wavelengths are shown in Figure 6, where the parameters are k = 1.392 and D = 5.545 µm for the R-R distribution and σ p = 0.869 and D = 4.255 µm for the L-N distribution.Similarly, the particle size distribution parameters that took into account the transmittance errors could also be solved.Therefore, the measurement error could be estimated from this.The comparison between the particle size measurement results of the R-R and L-N distributions obtained by the LEM and those of the Malvern laser particle size analyzer is shown in Figure 7.To view the size distribution of the small particles, the particle number probability density function (PDF) was also calculated, as shown in Figure 8.
From this comparison, it was found that the two hypothetical distributions were similar to the measurements of the Malvern laser particle size analyzer.In some applications, it was not the particle size distribution that was of interest but rather the mean particle size.The mean particle size parameters D[4, 3] and D [3,2] of these particle size distributions and the Malvern measurements are compared in Table 2.For the higher-order mean particle sizes, the LEM particle size measurement results converged with those of the Malvern laser particle size analyzer.The comparisons above show that assuming a two-parameter distribution of the particle size in the sample, these results were in good agreement with those of the Malvern laser analyzer.This indicates that relatively accurate information on the particle size distribution can be obtained using only the two wavelengths of the laser transmittance data.For this particle size measurement technology, the error mainly comes from two aspects.One is the transmittance measurement error [25], which could be suppressed by rigorous calibration and careful measurements.The other error was that the true shape of the particle size distribution of the measured particle cloud could not be known, and the true shape of the particle size distribution might be multi-parametric or even mathematically rigorous [40]; therefore, there seemed to be no good way to eliminate this error.

Particle Size Measurement in SRM Plume Flames
In this section, the particle size in the SRM plume was measured using the LEM.The transmitted laser passing through flame could be received with the help of off-axis parabolic mirrors in the receiver.In this experiment, a laser was operated in pulse mode to eliminate the influence of flame radiation on the measurement results.The detector signal for the 1064 nm laser at the receiver is shown in Figure 9 when the pulsed laser transmitted the SRM plume.In Figure 9, I on was the detector value when there was a laser output, and I o f f was the detector value when there was no laser output, as introduced earlier.From this figure, it can be seen that the detector value I o f f was very small, at about 6 mV, and almost equal to the detector's dark voltage, although the brightness of the SRM plume was very high.The reasons for this included the fact that there was a narrow bandpass filter in front of the detector and that the acceptance angle of the detector was very small.This resulted in a very limited effect of the SRM plume radiation on the detector at the receiver end.The effect of plume radiation could be eliminated by subtracting I o f f from I on .The data of the 532 nm laser was treated in the same way.The transmittances in the SRM flame were measured when the SRM was in steady operation, after the effects of flame radiation were eliminated.The transmittance for the 532 nm and 1064 nm lasers was 0.205 ± 0.022 and 0.263 ± 0.012, respectively, which corresponded with the optical thickness ξ of 1.58 and 1.34, as these two optical thicknesses ξ were both less than two and the multiple scattering effects could be ignored when measuring particle size using the LEM according to the conclusions in Section 3.1.
As in Section 3.2, we assumed that the particle size distribution in the SRM plume flame was either an R-R or L-N distribution, and the main procedure for solving these distribution parameters by plotting was almost the same as in Section 3.2.Somewhat differently, m total in Equation ( 22) was the mass flow rate of the alumina particles of SRM, which were estimated by design parameters, and here, V total was: where R pl is the radius of the flame in the measurement area, estimated from the flame image.V p is the mean axial velocity of the particles in the plume flame estimated by the simulation.Moreover, l in Equation ( 22) was 2R pl .
The solution results for R-R and L-N distributions for the mean transmittance at two wavelengths by plotting are shown in Figure 10.According to Figure 10, the parameter k in the R-R distribution was 1.13 and D was 8.32 µm; the parameter σ p in the L-N distribution was 1.51 and D was 7.45 µm.Based on the particle size distribution parameters solved above, the mean particle size D [3,2] was 1.0 ± 0.3 µm for the R-R distribution and 2.3 ± 0.4 µm for the L-N distribution.The volume PDF of the particles in the SRM plume flame is shown in Figure 11.These measurement results indicate that when the main particle size in the SRM flame is in the order of microns, the results were in very good agreement with those reported in the references [16,[41][42][43].

Figure 1 .
Figure 1.Experimental setup for multiple scattering study.(a) Schematic diagram of the experimental setup.(b) Photographs of part of the experimental setup.(c) Dimensions of triangular sample box.

Figure 2 .
Figure 2. Experimental setup for measuring the size distribution of particles in SRM plume flame.(a) Schematic diagram of the experimental setup.(b) Photographs of the transmitter.(c) Photographs of the receiver.
), where η was the only fitting parameter; thus, η was 9.2 × 10 −4 .This η value was then used to determine the predicted values under the four-flux model of the transmittance of the polystyrene particles for diameters of 2 and 3 µm.The experimental transmittance values for the different particle sizes and optical thicknesses ξ and the theoretical values of the four-flux model are shown in Figure 4.

Figure 4 .
Figure 4. Comparison of transmittance between experimental values and four-flux model prediction.

Figure 5 .
Figure 5. Q ext varies with particle diameter D p for 532 and a 1064 nm wavelength lasers.

Figure 6 .
Figure 6.Solution for f V (D) obtained by plotting.(a) Solution for R-R particle size distribution.(b) Solution for L-N particle size distribution.

Figure 7 .
Figure 7.The volume PDF f V (D) comparison between particle size measurement results of R-R and L-N distribution obtained by the LEM and those of the Malvern laser particle size analyzer.

Figure 8 .
Figure 8.The numerical PDF comparison between particle size measurement results of the R-R and L-N distribution obtained by the LEM and those of the Malvern laser particle size analyzer.

Figure 9 .
Figure 9. Detector signal for 1064 nm laser at the receiver when the pulsed laser transmitted the SRM plume.

Figure 10 .
Figure 10.Solution for f V (D) of the particle in SRM plume flame by plotting.(a) Solution for R-R particle size distribution.(b) Solution for L-N particle size distribution.

Figure 11 .
Figure 11.Volume PDF of particles in SRM plume flame.

Table 2 .
Comparison of measurement results of LEM with that of Malvern laser particle size analyzer.