Inverse Transient Radiative Analysis in Two-Dimensional Turbid Media by Particle Swarm Optimizations

Three intelligent optimization algorithms, namely, the standard Particle SwarmOptimization (PSO), the Stochastic Particle Swarm Optimization (SPSO), and the hybrid Differential Evolution-Particle Swarm Optimization (DE-PSO), were applied to solve the inverse transient radiation problem in two-dimensional (2D) turbid media irradiated by the short pulse laser. The time-resolved radiative intensity signals simulated by finite volume method (FVM) were served as input for the inverse analysis. The sensitivities of the time-resolved radiation signals to the geometric parameters of the circular inclusions were also investigated. To illustrate the performance of these PSO algorithms, the optical properties, the size, and location of the circular inclusion were retrieved, respectively. The results showed that all these radiative parameters could be estimated accurately, even with noisy data. Compared with the PSO algorithm with inertia weights, the SPSO and DE-PSO algorithm were demonstrated to be more effective and robust, which had the potential to be implemented in 2D transient radiative transfer inverse problems.


Introduction
The problem of transient radiative transfer in participating media has attracted increasing interest over the last two decades due to the emergence of ultrashort pulse laser and its application to the picosecond level or higher accuracy level of time-resolved techniques [1].Lasers with pulse durations of pico-to femtoseconds have been widely utilized in the areas of laser light-scattering flame diagnose, the prediction of the temperature distribution in combustion chambers, atmospheric science, remote sensing, nondestructive testing, laser machining, and medical diagnostics, to name a few [2][3][4].In recent years, the medical applications of short pulse laser, such as the laser tissue welding [5,6], the laser tissue ablation [7,8], the optical tomography in medical imaging [9], and the photodynamic therapy [10], have drawn more and more attentions all over the world.In particular, the short pulse laser of near-infrared spectrum is extensively used in noninvasive techniques to retrieve the characteristics of the semitransparent media such as the biological tissues and turbid media [11][12][13][14].Theoretically speaking, the research emphases of noninvasive reconstruction technology are focused on using the model-based iterative imaging reconstruction techniques, which employ a precise and efficient numerical forward model based on solving complete transient radiative transfer equation (TRTE) and then establish a proper inverse algorithm to retrieve the optical property parameters or internal geometry of the medium by using the measured time-resolved transmittance and reflectance signals.Compared with the traditional measurement methods, its biggest advantages are noninvasive, in situ measurement, containing a lot of time-resolved information about internal composition and characteristic of the participating media.
For solving the inverse transient radiation problems in the noninvasive reconstruction, the first step is to solve the direct TRTE model accurately and efficiently.Nowadays, several numerical strategies have been developed to solve the forward model, including the discrete ordinate method (DOM) [15][16][17][18], the finite volume method (FVM) [18][19][20][21], the integral equation (IE) method [22], the finite-element method FEM [23,24], and the Monte Carlo method [25].Each of these methods has its own relative advantages and disadvantages and none of them is superior to others in all aspects.However, many numerical complications may 2 Mathematical Problems in Engineering arise from incompatibilities between the radiative transport model and the discrete formulation employed for the fluid dynamics and other heat transfer modes.The FVM for RTE avoids many of these complications and has the advantage of greater compatibility with existing finite volume based heat transfer software.Besides, the FVM also could simplify the problems of complex geometries using the unstructured meshes.Through decades of development, FVM has become a sophisticated technology in mechanics and thermal analysis.There have been several investigations focused on the interaction of ultrafast radiation with participating media using FVM [19,20].In these studies, both the homogeneous and inhomogeneous media have been considered.However, most of the studies of inverse radiation problems concentrate on one-dimensional (1D) cases and there are very few researches about multidimensional problems, especially the inverse transient radiative problems.
To solve the inverse problems, a wide variety of inverse techniques have been successfully employed in the inverse radiation analyses which can be grouped into two categories, that is, the traditional algorithm based on gradient and the intelligent optimizations [26].To date, many gradient-based techniques have been employed in the inverse radiation analysis such as the constrained least-squares method [27], the conjugate gradient (CG) method [28], and Levenberg-Marquardt method [29].However, all these traditional methods depend on the initial value or the derivatives and gradients which are difficult to be solved accurately by numerical simulation in some cases.Furthermore, if the initial value is not chosen properly, the solution may be infeasible in the original domain, or in a worse case, not be convergent.In addition, these conventional methods using derivatives are generally restricted to the conditions of continuity, sensitivity, convexity, monotonicity, and nonlinearity of both objective functions and constraints [30].Compared with traditional gradient-based methods, the intelligent optimization algorithms have a couple of outstanding characteristics; that is, the ill-posed inverse problem could also be solved, the derivative of the objective function is not necessary, and the a priori information is not needed.A typical characteristic of the intelligent optimization methods is that they can solve the global optimal problem reliably and obtain high quality global solutions with enough computational time, especially for higher problem dimensions [31].This is primarily due to the fact that the random search technique is used to modify locally examined regions, which gives the algorithm a global optimization capability.The Particle Swarm Optimization (PSO), as an intelligent swarm optimization, first introduced in 1995 by Eberhart and Kennedy [32], has been studied extensively by many researchers in recent years.It is a potential heuristic bionic evolutionary algorithm originally inspired by the locking behavior of birds.Generally speaking, PSO is characterized to be simple in concept, easy to implement, and computationally efficient.Early studies have found the implementation of PSO to be effective and robust in solving problems featuring nonlinearly, nondifferentiability, and high dimensionality.Many modifications have been made to improve the convergence rate of the original PSO algorithm.More recently, our research group has done several studies on the PSO algorithm to solve the one-dimensional inverse radiation problems [33][34][35][36][37].However, to the best of our knowledge, few researches are presented to solve the multidimensional inverse transient radiation problems by the intelligent optimization methods.
The objective of present work is to apply the PSO-based algorithms to solve the inverse transient radiation problem in two-dimensional (2D) participating medium.The optical properties, the size, and location of the circular inclusion were retrieved, respectively.The remainder of the paper is organized as follows.The detailed mathematical formulation and computational steps of the direct model are described in Section 2. The theoretical models of the standard PSO, the Stochastic PSO, and the hybrid Differential Evolution and Particle Swarm Optimization (DE-PSO) algorithms are described in Section 3. The solving model for direct problem and the inverse transient radiation analysis by the PSO-based algorithms are examined in Section 4. In Section 5, the main conclusions and perspectives are provided.

The Direct Model
In the present work, the thermal effect caused by the incident laser was ignored and only the transient radiative transfer was considered in the numerical model.Thus, the transient radiative transfer equation can be expressed as [38] 1 where  is the length along the direction Ω, radiative intensity  is the function of position , direction Ω, and time , ,   , and   denotes extinction, absorption, and scattering coefficients, respectively, and  =   +   .Φ(Ω  , Ω) represents the scattering phase function between incoming direction Ω  and scattering direction Ω. Ω  is the solid angle in the direction Ω  .In the 2D transient radiative transfer problem caused by collimated irradiation, the radiative intensity within the medium is separated into two parts [35]: (a)   , the remnant of the collimated beam after partial extinction, by absorption and scattering, along its path; (b)   , the diffuse part, which is the result of emission from the boundaries, emission from the medium inside, and the radiation scattered away from the collimated irradiation.Thus, the radiative intensity in the media can be written as For the square pulse laser, the intensity within the participating media can be expressed as [35] where  is Heaviside function [35];  denotes Dirac function [35];  in is the heat flux of incident laser at boundary  = 0; Ω 0 is the incident direction; and  0 represents the geometric distance along the direction Ω 0 .
For the 2D participating media, the discrete equation could be obtained using the FVM model: where In ( 4),   =   +   denotes the total source term.The source terms   and   result from the collimated radiation   and the diffused radiation   , respectively.For the diffuse reflection boundary, the boundary condition is given as [35]  (  , Ω, ) =     (  , ) where (  , Ω, ) and (  , Ω  , ) are the outgoing and the incoming intensities at the boundary, respectively; n is the outward unit normal vector at the boundary;   is the emissivity of the boundary.Substituting the length-dimensional time term  * =  into (4), we can get The total source term   can be depicted as where Ω 0 denotes the direction of incident radiation.Using the fully implicit scheme, the diffusion intensity   for the control volume  in the direction s  at time  * can be written as where  represents the velocity of light; Δ * is the time step;   represents the direction weight;  and  are set as Δ * /(1+ Δ * ) and 1/(1+Δ * ), respectively.  , and   , denote the -direction radiative intensities of node  at the upstream boundaries of control volume which are along axis  and axis , respectively.The time-domain thermal signals, that is, the transmittance signal   and reflectance signal   , can be expressed as where  is the direction cosine and   denotes the direct transmission flux of the collimated light.The FVM was chosen to solve the equation of TRT.For the sake of simplicity, the details of FVM are available in [20] and are not repeated here.

The Inverse Model
Commonly, there are three ways to obtain the internal information of media, that is, the continuous wave method, the time-domain method, and the frequency-domain method [39].All these methods have the similar procedure in the experimental studies and the numerical simulations.Taking the time-domain method as an example, the reconstruction scheme consists of three major parts calculation is performed to get the corresponding detector predictions.The reconstruction process is completed when the value of the objective function is less than a preset value or the number of iterations exceeds the maximum number of iterations.Foremost, the forward model must be solved precisely enough so that the measurement data obtained by detectors could be simulated correctly.Consequently, the forward model based on the complete transport-theory TRTE should be utilized [40].The details of the inverse optimization algorithms are shown as follows.

Theoretical Model of Basic PSO Algorithm.
The PSO algorithm was first introduced by Eberhart and Kennedy [32] inspired by social behavior simulations of flocking birds.The flocking population is called a swarm, and the individuals are called the particles.The basic idea of PSO can be described simply as follows.Each particle in a swarm represents a potential solution which changes its position and updates its velocity based on its own and its neighbors' experience.The velocity and position of the th particle at generation  + 1 could be expressed as [32] V where  1 and  2 represent two positive constants called acceleration coefficients;  1 and  2 are uniformly distributed random numbers in the interval [0, 1];  is the inertia weight factor which is used to control the impact of the previous velocities on the current velocity.In the present work,  is defined as  = 1.0 − 0.6 × /, where  represents the maximum generation number; X  () = [ 1 (),  2 () ⋅ ⋅ ⋅  in ()] T denotes the present location of particle , which represents a potential solution; T is the present velocity of particle , which is based on its own and its neighbors' flying experience; P  () = [ 1 (),  2 (), . . .,  in ()] T is the "local best" in the th generation of the swarm.The global best position of the swarm is expressed as P  .

Theoretical Model of DE-PSO Algorithm.
In order to avoid premature convergence of the standard PSO and ensure the overall convergence, a considerable number of modified PSO algorithms were proposed in recent years, like the Stochastic PSO (SPSO), the Multiphase PSO (MPPSO), the Quantum-Behaved PSO (QPSO), the hybrid Ant Colony Optimization and PSO (ACO-PSO), and so forth [33,36,41].All these improved PSO-based algorithms are proved to be adaptive and robust parameter-searching techniques.In the present work, the SPSO and the hybrid Differential Evolution and PSO (DE-PSO) are utilized in the inverse procedure [42].
The DE-PSO is the combination of two stochastic algorithms and is expected to accelerate the convergence process by introducing the DEA which has strong global searching ability to the PSO.The detailed model description of SPSO is The object function Newly generated vector available in [41] and will not be repeated here.The theoretical model of DE-PSO is described as follows.
The differential evolution algorithm (DEA) is a novel parallel, direct, and stochastic searching method which was proposed to solve the continuous domain problems in 1995 [43].The basic idea of DEA is to use the differential of two individuals in the current generation, namely, vectors  2 and  3 , to combine with a third vector  1 to generate a new particle, which is called the trial vector.Then, the objective function was calculated to decide whether the trial vector is to replace vector  or not.The trial vector T is generated according to [43] where  is a real and constant factor which controls the amplification of the differential variation [X  2 ()−X  3 ()], and in this paper it was set as 0.5.The integers  1 ,  2 , and  3 are chosen randomly from the interval [1, ] and are different from running index .Other parameters are defined as same as those of the PSO. Figure 1 illustrates the vector-generation process defined by (13).
In the process of DE-PSO, the differential evolution operator is introduced to the position updating equation, which can be expressed as [42] where  = [X  2 () − X  3 ()], which is named as differential evolution operator;  min represents the minimum value of objective function in the present generation;  DE denotes a preset constant value which satisfies  DE ≤  min .
The main function of the third term on the right side of ( 14) is to avoid local convergence.However, the differential evolution operation may lead to instability of the algorithm, so only when the objective function of the new particle is smaller than that of the last generation, ( 14) is applied as the position update equation.Otherwise, the position updates according to (12).The implementation of the approach for solving the inverse transient radiation problem by DE-PSO algorithm can be carried out according to the following routine.
Step 1. Initialize the position X  () and the speed V  () of each particle and calculate the objective function.Afterwards, the global best position P  and every local best position P  () can be determined.
Step 2. Check if the value of the objective function is smaller than the preset parameter  or the maximum number of iterations is reached.If so, go to Step 6; otherwise, go to the next step.
Step 3. Update the velocity and position of each particle according to (11) and ( 14), respectively.Each velocity and position component is bounded by the maximum and minimum values to which each particle in the space is constrained.The range of the particle position is decided by the physical situation; for example, the single scattering albedo can only be set in the range of [0, 1].Then, calculate the objective function of each particle and update the global beat position and local best positions.
Step 4. Calculate the objective function of each newly generated particle.If the value of the objective function is smaller than that of the corresponding particle in the last generation, then the new particle can be maintained.Otherwise, regenerate the particle according to (12).
Step 5. Update the generation from  to  + 1 and go to Step 2.
Step 6.Output the global best position and its corresponding value of objective function.
To illustrate the difference between PSO, SPSO, and DE-PSO algorithms more clearly, the flowcharts of the three algorithms are shown in Figure 2.

Inverse Analysis
Procedure.The transient radiative inverse problems for estimating the internal radiative properties of the 2D participating media are solved by minimizing the objective function which can be defined as where  represents the different wall of the 2D media,  ,mea (a) denotes the measured time-resolved transmittance or reflectance signals, which will be simulated by the forward model using the FVM, and  ,est (a) is the estimated timeresolved transmittance and reflectance for an estimated vector a = ( 0 ,  1 , . . .,   ) T .Moreover,  1 and  2 are start and end point of the sampling span, respectively.
To demonstrate the effects of measurement errors on the predicted terms, the random errors were considered.In the present paper, the simulated measured signals with random errors were obtained by adding normally distributed errors to exact reflectance and transmittance as follows: where  mea and  exa represent the measured reflectance or transmittance signals with and without noisy data, respectively, and rand( ) denotes a normal distributed random variable with zero mean and unit standard deviation.The standard deviation of the measured values is denoted by , which is defined as where % is the measured error and 2.576 arises from the fact that 99% of a normally distributed population is contained within ±2.567 standard deviation of the mean value.The absorption and scattering coefficients or the inclusion locations were estimated by minimizing the objective function using SPSO and DE-PSO algorithms.The inverse algorithms were adopted for minimization of the objective function, which is also called the fitness function.The stop criterion of the inverse algorithms is that the number of iterations exceeds the maximum number of iterations or the best fitness function is less than a specified small value.

Results and Discussion
4.1.Homogeneous Media.Consider a 2D homogeneous semitransparent medium with length 1.0 m × 1.0 m whose left side exposes to a collimated square pulse laser beam in the -direction with pulse width   = 0.1 m.The absorption coefficient   and the scattering coefficient   are set as 0.02 m −1 and 9.98 m −1 , respectively.The emissions of the medium and its boundaries are neglected.The media is assumed to be anisotropic scattering with black boundaries and the phase function is defined as where  0 = 1,  0 = 0, and  0 = −1 represent the backward scattering, isotropic scattering, and forward scattering, respectively., , and  and   ,   , and   denote the direction cosines of the incoming direction Ω and scattering direction Ω  , respectively.The control volumes were set as   ×   = 500×500 and solid angles were divided into   ×  = 20×10 in the FVM.The total observed time span  was 10 m.As shown in Figure 3, the results of FVM model are compared with the solutions of Discrete Ordinate Method (DOM) [16] and Least Square Finite Element Method (LSFEM) [44].It can be seen that the results calculated by FVM approximation in the present study are in good agreement with those in [16,44] which means the FVM model in this paper is proved to be valid.

Nonhomogeneous Media with Circular Inclusions.
The estimation in this paper is expected to be applied in the nondestructive detection of human or small animals.Figure 4 shows the model of a rat liver with tumor which is simplified to a rectangle medium with a circular inclusion.To illustrate the influence of the inclusions to the radiative signals, four cases are investigated in this section that is, two circular inclusions with different radiuses, location and optical properties in a rectangular medium.The media size is 1.0 m × 1.0 m which exposed to a pulse laser with pulse width   = 0.1 m.The absorption coefficient and scattering coefficient of the 2D medium are set as  0 = 2.0 m −1 and  0 = 3.0 m −1 .The geometric locations of the two inclusions are identified by three parameters, namely,  1 ,  2 , and , among which  1 and  2 represent the coordinate values of the inclusion center and  denotes the radius of the circular inclusions.Table 1 shows the control parameters of the four cases.It is worthy noticing that the four cases are designed on the principle that, compared to Case 1, the other three cases only change one set of parameters, that is, the size, location, or optical properties, for both inclusions.For the case of clarity, the diagrams of the four cases are shown in Figure 5.  were considered.For the laser incidence, two different conditions were included (see Figure 7): (1) point-irradiated (PI) condition which means one single laser beam irradiates one spot of the media wall; (2) face-irradiated (FI) condition which means the laser beams incident on the whole boundary wall.In this section, the two cases were proceeded under both two laser incident conditions.In Case 1, the inclusion was in the center of the medium which means  1 =  2 = 0.5 m and  = 0.25 m.The responses of detectors D1, D2, D3, and D4, which are shown in Figure 8(a), were simulated by the forward model.The properties of the inclusion were estimated and the inverse analysis was taken by using the PSO, SPSO, and DE-PSO algorithms.Meanwhile, the optical properties of the 2D medium were taken as  0 = 0.20 m −1 and  0 = 0.30 m −1 , respectively.The true values of the inclusion properties were  1 = 1.0 m −1 and  1 = 9.0 m −1 , respectively.The searching intervals were set as [0, 10 m −1 ].The estimated results are listed in Table 2.In Case 2, two inclusions inside the medium were considered as shown in Figure 8(b).The optical properties of both inclusions were estimated by PSO, SPSO, and DE-PSO.In this case, the optical properties of the medium were taken as  0 = 0.20 m −1 and  0 = 0.30 m −1 , respectively.The true value of the inclusions which needed to be estimated was ( 1 ,  1 ) = (1.0 m −1 , 9.0 m −1 ) and ( 2 ,  2 ) = (4.5 m −1 , 5.5 m −1 ), respectively.The size and location of the two inclusions were ( 1,1 ,  2,1 ,  1 ) = (0.3, 0.3, 0.2) m and ( 1,2 ,  2,2 ,  2 ) = (0.7, 0.6, 0.15) m.The retrieving results of optical parameters are shown in Tables 3 and 4.
It can be seen from Table 2 that  1 and  1 can be retrieved accurately under both PI and FI conditions by PSO, SPSO, and DE-PSO when there is only one inclusion.Furthermore, it also demonstrates that the retrieving results under FI condition are more accurate than those under PI condition, the reason of which is that, under FI condition, the reflectance and transmittance signals are stronger and are carrying more information available.Meanwhile, the DE-PSO algorithm shows its priority in reconstructing the inclusion properties.Under both conditions, the relative errors of retrieval results are smaller than those of the SPSO and PSO.more accurate than those of the PSO.Furthermore, under FI condition, for Case 2, the optical properties of both inclusions can be estimated accurately without measurement errors and with the fact that the estimated relative errors of results for inclusion B (see Figure 8) are larger than those of inclusion A. However, the calculation time required of the DE-PSO is more than that of the SPSO and PSO.Thus, it can be concluded that SPSO and DE-PSO can retrieved the optical properties in a 2D medium accurately but with longer calculation time.The value of the objective function of PSO, DE-PSO, and SPSO with two inclusions is shown in Figure 9.
It can be seen that the PSO fall into local convergence quickly and stop searching.In the earlier stage of the inverse process the objective function of SPSO declines faster than that of DE-PSO.However, in the later period of the inverse process, the SPSO almost stops searching and the objective function of DE-PSO continues declining which leads to the result that the retrieving results of DE-PSO are more accurate than those of the SPSO.

Estimation of Geometrical Parameters.
In the inverse procedure, the long sampling span will reduce the computational efficiency.Hence, it is of great significance to select the proper sampling span.Thus, the sensitivity of radiative signal to the inverse parameters needs to be analyzed.The sensitivity coefficient is one of the most important characteristic parameters in the sensitivity analysis, which is the first derivative of the radiative signals to a certain inverse parameter.The sensitivity coefficient is defined as where   denotes the independent variable which stands for  1 ,  2 , and ; Δ represents a tiny change;   is the radiative signals of each boundary.The sensitivities of the inclusion location ( 1 ,  2 ) and size  for the signals D1-D4 are shown in Figure 10.The parameters are set the same as Case 1 described in Section 4.3.1.It demonstrates that the sensitivity coefficient is relatively larger in the span of  ∈ [0, 3 m].So this span can be chosen as the sampling span for inverse analysis.The same results are expected in other cases.
In many conditions, the shapes and optical properties of inclusions in the media are known in the actual application cases of reconstructing the medium's inside information.For example, a tumor or cyst in the biological tissue is spherical or spherical-like, and the defects in the specimen, such as the air pore, are mainly the shape of a cube, sphere, or ellipsoid.In addition, the optical properties of human and animals tissues, in vitro or in vivo, are studied thoroughly [45].Consequently, the geometrical feature and optical properties of inclusions could be used as a priori information in these cases, which will reduce the complexity of the reconstruction technique and improve the quality and efficiency of the reconstruction results.Therefore, on the condition that the inclusion shape and the optical properties are known, the inverse problem turns into the retrieval of the size and location of the inclusions.
In this section, the inclusion location ( 1 ,  2 ) and size  of cases ( 1 =  2 = 0.25 m,  = 0.25 m), ( 1 =  2 = 0.5 m,  = 0.25 m), and ( 1 =  2 = 0.75 m,  = 0.25 m) were estimated simultaneously.The size of the 2D medium was set as 1.0 m × 1.0 m.The absorption coefficient  0 and scattering coefficient  0 of the medium were 0.2 m −1 and 0.3 m −1 , respectively.There was one circular inclusion in the 2D medium with absorption coefficient  1 of 1.0 m −1 and scattering coefficient  1 of 9.0 m −1 , but the inclusion size and location were unknown, which needed to be retrieved.The FI incidence type and the data of the detectors 1, 2, 3, and 4 were used for inverse analysis by using the PSO, SPSO, and DE-PSO algorithms with measurement errors.The search span of  1 ,  2 , and  was set as [0, 1] m.The results are shown in Tables 5 and 6. Figure 11 shows the comparison of the exact and reconstructed reflectance signal when  1 =  2 =  = 0.25 m with 10% measurement error.
It can be seen that the size and location can be estimated accurately using PSO, SPSO, and DE-PSO without measurement errors under FI condition.Furthermore, the results of DE-PSO are more accurate than those of the SPSO and PSO.The relative errors are reasonable even with 10% measurement error.However, it is worth noting that the relative errors of DE-PSO with 10% measurement error are bigger than those of the SPSO, which means the SPSO algorithm is more robust than the DE-PSO.

Conclusions
In the present study, the standard PSO, SPSO, and DE-PSO algorithms were applied to estimate the size, location, and optical parameters of the circular inclusions in a 2D rectangular medium.The conclusion was obtained that the retrieving    results of absorption coefficient and scattering coefficient under FI condition are more accurate than those under PI condition, the reason of which is that, under FI condition, the reflectance and transmittance signals are stronger and are carrying more information available.In addition, the DE-PSO algorithm shows its priority in reconstructing the inclusion properties.When only one inclusion is considered, under both laser incident conditions, the results' relative errors of DE-PSO are smaller than those of the SPSO.However, when estimating the size and location of inclusion media, the SPSO is more robust than DE-PSO.Furthermore, compared with the performance of standard PSO, both the SPSO and DE-PSO are proved to be more accurate and robust, which have the potential to be applied in the field of 2D inverse transient radiative problems.

Nomenclature a:
The vector of estimated properties : The speed of light, m/s  1 ,  2 : The acceleration coefficients of PSO : The control variable of DEA : The objective function : The radiation intensity, W/(m ⋅ sr) : Thenumberofthepopulation   ,   : The number of grids   ,   : The number of polar angles and azimuthal angles P  (): The global best position discovered by all particles at generation  P  (): The local best position of particle  discovered at generation  or earlier : The radius of the inclusions  1 ,  2 : The uniformly distributed random numbers : The sensitivity coefficient : Time or generation in PSO algorithm   : The incident pulse width, s  1 ,  2 : The start and end point of sampling span Δ: Thetimestep,s   (): The trial vector of DEA V  (): Thevelocityarrayoftheth particle at the th generation in PSO : The inertia weight coefficient X  (): Thepositionarrayoftheth particle at the th generation in PSO  1 ,  2 : The coordinate values of the inclusions' center.

1 𝑐𝜕𝐼
(, Ω, )  +   (, Ω, )  = −  (, Ω, ) +   (, Ω, ) +   (, Ω, ) = −  (, Ω, ) +   (, Ω, ) , [40]:(1) a forward model that predicts the detector signals based on the solution of the transient radiative transfer equation; (2) an objective function that provides criterion of the differences between the detected and the predicted data; (3) the reconstruction technique which can minimize the objective function to get new guesses of the estimated parameters such as the optical properties or the geometrical parameters.Based on the new guesses of the optical properties, a new forward

Figure 1 :
Figure 1: The schematic diagram of vector-generation process for DEA.

Figure 8 :Figure 9 :
Figure 8: The schematic of pulse laser irradiation and the measurement position of time-domain signals of the two test cases.

Figure 10 :
Figure 10: The sensitivity of the radiation signals for the geometric parameters of the inclusion.

Figure 11 :
Figure 11: The comparison of the exact and reconstructed reflectance signals.

Table 1 :
The control parameters of the four cases.

Table 2 :
The retrieving results of the optical properties for Case 1 using PSO, SPSO, and DE-PSO without measurement error.
Table 3 demonstrates the retrieving results of the optical properties for Case 2 of different generations.As shown in Table 4, the results of SPSO and DE-PSO are Mathematical Problems in Engineering

Table 5 :
The inverse results of size and location of circular inclusion using PSO, SPSO, and DE-PSO without measurement error under FI condition.

Table 6 :
The inverse results of size and location of circular inclusion using PSO, SPSO, and DE-PSO with 10% measurement error under FI condition.