## 1. Introduction

Radiative heating devices are encountered in various industrial fields, such as industrial boilers, spacecraft, infrared reflecting ovens, metallurgical equipment, and so on and the design of radiative enclosure has a direct impact on security issues [1]. Inverse design technique is a new method in recent years, whose solving process is to establish an objective function according to the design requirements at first, and then optimize the objective function by some optimization methods, and achieve the pre-specified purpose finally. Inverse design technique has the advantages of simple process, short design circle, good optimization results, etc., and it has got more and more attentions and applications.

Inverse design problems can be divided into two categories according to the design prerequisites [2]. One is inverse boundary design problems, in which the geometry shape of the radiative enclosures is fixed and the boundary conditions are need to be deigned [3–6]. The other is inverse geometry design problems, where the boundary conditions are predetermined and the geometry shape of the design surface needs to be designed [7, 8]. For the reason that the geometry shapes of the boundaries of radiative enclosures are different at each iteration and the grids in the computational domain must be re-meshed in every iterative calculation, the inverse geometry design of radiative enclosures is the most complex inverse radiative problem [2, 9].

During the last few decades, some inverse design techniques have been successfully used for solving inverse geometry design problems. Howell et al. [10, 11] proposed inverse design ideas and applied inverse Monte Carlo techniques, Tikhonov method, truncated singular value decomposition (TSVD), Modified TSVD (MTSVD), artificial neural networks (ANN) and conjugate gradient method (CGM) to solve the inverse design problem of a three-dimensional industrial furnace, which greatly improved the practical design of thermal and environmental systems. Franca and Howell [12] studied a transient inverse design problem that finding the optimal location of a heater on the top surface of a three-dimensional enclosure to produce a prescribed time-dependent temperature distribution on the bottom surface of the enclosure, the TSVD method is used to regularize the ill-conditioned system of linear equations the pre-specified temperature curve is obtained with an error of less than 1.0%. Tan et al. [13] applied meshless method to solve the coupled conductive and radiative heat transfer problem in heating devices, in which a series of nodes are used for discretizing the computational domain to overcome the tedious re-mesh works, and CGM is adopted to optimize the height of the adiabatic diffuse reflection surfaces and the geometry shape of the heating surface to satisfy the required total heat flux on the pre-appointed region of the low temperature heated surface. Sarvari and Mansouri [14] used CGM to minimize the objective function which is expressed as the sum of square residuals between estimated and desired heat fluxes on the design surface to satisfy the specified temperature and heat flux distributions. The radiative heat transfer problem in the two-dimensional irregular enclosure filled with gray participating media with uniform absorption coefficient is solved by discrete transfer method, the effect of optical depth and angular refinement on the inverse design results are also investigated.

However, it can be found that most of the above researches about inverse geometry design problems are solved based on the gradient-based methods. All these methods have the common disadvantages that the computational process of getting the gradient is complicated and the retrieval results are strongly depended on the initial guessed values. On the contrast, intelligent algorithms can abolish the complex calculations about gradient and randomly generate the potential solutions in the search space to overcome the drawback of depending on the initial value. In recent years, some intelligent algorithms have been successfully applied for solving inverse radiative problems, including PSO, genetic algorithm (GA), ant colony optimization (ACO), difference evolution (DE), fruit fly optimization algorithm (FOA), to name a few [15–22]. Compared with conventional techniques, intelligent algorithms can obtain much more potential solutions at each iteration and all searches are executed in parallel, which greatly improved the computational efficiency, especially for solving some high dimensional problems. Intelligent algorithms also have been successfully used to solve inverse design problems. For example, Moparthi et al. [23] solved the coupled radiative and conductive heat transfer problem in one-dimensional planar system based on finite volume method (FVM) and lattice Boltzmann method (LBM) and applied GA to optimize the heater temperature to produce the desired heat flux and temperature distributions on the design surface. The retrieval results show that the temperature or heat flux of the heater surface has a significant effect on the design surface condition, the medium properties and the distance between the two surfaces. Amiri et al. [24] adopted modified discrete ordinate method to solve the radiative transfer equation (RTE), the micro genetic algorithm (MGA) is employed to optimize the objective function which is defined as the sum of the square of the difference between estimated and desired heat fluxes on the design surfaces. The design purpose is finding the best number and locations of the discretized heaters to meet the desired temperature and heat flux distributions on the design surface. Sarvari et al. [25] discretized the design surface into a series of control points and used *B*-spline to approximate the geometry shape of the boundary, the MGA is employed to optimize the locations of the control points to produce a desired heat flux distribution on the temperature-specified surface. The effects of corresponding parameters on the inverse design results are also investigated and the angular meshes are recommended as

The PSO algorithm is a kind of biologically inspired algorithm whose search process is similar to foraging of birds and it was proposed in 1995 by Eberhart and Kennedy [26]. The physical model of PSO is very simple and the computational program is easy to be implemented, it also has strong robustness and achieves good performance on computational efficiency and accuracy. In addition, PSO algorithm can well balance the global and local search of particles, which enhance the global convergence of the algorithm. Farahmand et al. [27] investigated the inverse geometry design of two-dimensional radiative enclosures with diffuse gray surfaces based on the PSO and the retrieval results show that PSO algorithm obtains better performance in satisfying the design goal based in terms of computational accuracy and CPU time compared with MGA. However, the standard PSO algorithm also suffers from easily trapping into local optima in solving high dimensional problems. In order to strengthen the applicability of PSO, some improvements have been proposed and widely applied, including Stochastic PSO (SPSO), Differential Evolution PSO (DEPSO), Multi Phase PSO (MPPSO), and so on. However, to the authors’ best knowledge, there are few reports concerning about the applications of improved PSO algorithms for solving inverse geometry design problems of radiative enclosures.

In this chapter, the application of PSO algorithms in solving inverse geometry design problems of two-dimensional radiative enclosures filled with participating media is investigated. The design goal is to satisfy a uniform distribution of radiative heat flux on the designed surface. The discrete ordinate method (DOM) with a body-fitted coordinate system is used to solve the RTE. The standard PSO, SPSO and DEPSO algorithms are applied to optimize the locations of the control points, and Akima cubic interpolation is adopted to obtain the boundary geometry shape through these points. A typical inverse geometry design test is studied to demonstrate the good performance of PSO algorithms and the effects of corresponding parameters are also discussed.

The remainder of this chapter is organized as follows: the theoretical principles of PSO algorithms are introduced in Section 2. The feasibility of PSO algorithms by four famous benchmark functions is verified in Section 3. The inverse geometry design of two-dimensional radiative enclosures and the influences of the number of control points and the radiative properties of media on the inverse design results are investigated in Section 4. The main conclusions of the researches in this chapter are summarized in Section 5.

## 2. Theoretical overview of PSO algorithms

Bird individuals will communicate with each other to share their information about food when they are foraging, which can help birds find food faster. The advantages of cooperation of bird swarm are much greater than the disadvantages of competition among bird individuals. Based on the features of bird foraging behavior, Kennedy and Eberhart proposed PSO algorithm in 1995 [26]. The solving process of PSO algorithm is similar to the foraging behavior of birds, and the corresponding relationships are shown in **Table 1**.

There are two dominant parameters in the PSO algorithm, namely, the speed and the location of particles. The moving speed decides the direction and distance of the particles, and every location of particles can be considered as the potential a solution of optimization problems. PSO adopts a combination of local and global searches and shares the evolutionary information among particle individuals to find the optimal solution.

### 2.1. Basic PSO algorithm

At the beginning of the optimization of PSO, every location and velocity of particles are randomly generated. During each iteration, there will be two extreme values. One is best location that an individual particle found so far, which is called local best location. Another is the best location that the whole particle swarm found so far, which is called global best location. The velocity and location of each particle are stochastically accelerated according to these two extremes and the evolutionary formula can be expressed as follows [26]

where **V**_{i}(*t*) and *i*th particle at iterations *t* and *t*+1, respectively, and **V**_{i} ∈ [−*V*_{max}, *V*_{max}]. **P**_{i}(*t*) and **P**_{g}(*t*) indicate the local best location and global best location, respectively. **X**_{i}(*t*) denotes the location of *i*th particle, which is depending on the search experience of *i*th particle and surrounding particles. The evolutionary formula of *i*th particle’s location is defined as [26]

According to Eq. (1), we can find that the velocity of *i*th particle consists of three parts: the first part on the right side is the current velocity of *i*th particle, which can counterpoise the local search and the global search; the second part on the right side denotes the influence of the search memory of *i*th particle, which makes individual particle has the ability of global search; the third part on the right side indicates the influence of the cooperation among particles. The updating process of the *i*th particle’ location is shown in **Figure 1**.

The main procedure of PSO algorithm for solving optimization problems can be carried out according to the following steps:

**Step 1**: Initialization: Randomly initialize the location and the velocities of particle of every particle in the searching space, input the number of particle swarm, the maximum iteration numbers*t*_{max}and the stop criterion*ε*, etc. Set the current iteration as*t*=1.**Step 2**: Fitness evaluation: Evaluate each particle’s fitness according to its location and determine the local best location**P**_{i}(*t*) and global best location**P**_{g}(*t*).**Step 3**: Updating: Update the velocity and the location of each particle according to Eqs. (1) and (2), respectively. Calculate the new objective function value of each particle and update the local and global best locations**P**_{i}(*t*) and**P**_{g}(*t*).**Step 4**: Comparison: Compare the objective function value of each newly obtained particle with the corresponding values at the last iteration. If the new objective function value is better than the one in the last generation, then the new location and velocity of this particle is updated. Otherwise, the new location is abandoned.**Step 5**: Repeating: Check whether one of the following two stop criterion is reached: (1) the objective function value is less than the value of*ε*, and (2) the iteration number reaches the maximum iteration number. If so, go to the next step; otherwise, go to**Step 3**.**Step 6**: Update of iteration number: Update the number of iteration from*t*to*t*+1.**Step 7**: Termination of iteration: Output the global best optima and its corresponding results of optimization problems and then stop the calculation.

The flowchart of the basic PSO algorithm is shown in **Figure 2**.

However, there are some obvious shortcomings in the basic PSO algorithm, such as slow convergence, easy to fall into local optimum, and even the velocities tend to be infinity in some occasions. Thus, many modifications have been proposed to overcome these drawbacks.

### 2.2. Standard PSO algorithm

There are two important capabilities in PSO algorithm, namely exploration and exploitation of particles. Exploration is the phenomenon that particles leave the original orbit and search for new space. Exploitation is the phenomenon that particles look for better locations along the original track. In order to better take advantage of these two search way, Shi and Eberhart put forward the standard PSO algorithm on the basis of basic PSO in 1998 [28], in which an inertia weight coefficient *w* is introduced to control the impact of the current velocity on the next velocity. The velocity formula of *i*th particle can be expressed as [28]

where *w* is the inertia weight coefficient, which can directly affect the balance between the global and local exploration abilities. At the initial stage of the search process, a big inertia weight coefficient is recommended to improve the global exploration ability in the relatively large space; whereas the inertia weight should be reduced with the iteration number increases to strengthen the local exploitation ability. It is worth pointing that a linearly decreasing inertia weight coefficient can successfully prevent particles from oscillating near the global best location [29]. Therefore, the inertia weight coefficient can be defined as

Comparing Eq. (1) with Eq. (3) we can find that basic PSO algorithm is a special circumstance of standard PSO algorithm that inertia weight is set as *w*=1. The searching efficiency is significantly improved by introducing inertia weight *w*. However, the proportional relation of particle velocities at every generation is not all the same, and the standard PSO algorithm can’t be successfully applied for solving some complicated optimization problems. For breaking through the limitation of PSO, many modified techniques have been developed out and widely applied in engineering fields.

### 2.3. Stochastic PSO algorithm

In order to overcome the drawback that PSO algorithm converges too early and make sure to reach the goal of global convergence, Zeng and Cui proposed SPSO algorithm in 2004 [30, 31], in which a stopped changing particle is utilized to improve the global searching ability of particle swarms.

In SPSO algorithm, the inertia weight coefficient is set as *w*=0. Hence, the velocity of *i*th particle at *t*+1 iteration is determined by three parameters of *t* iteration, namely **X**_{i}(*t*), **P**_{i}(*t*) and **P**_{g}(*t*). The new velocity of *i*th particle can be expressed as [30]

According to Eq. (5), it can be found that the local searching ability of SPSO is increased compared with standard PSO. However, the global searching ability is reduced significantly. In order to further strengthen the global searches of SPSO, the algorithm randomly generates a particle in the searching space whose location is **X**_{j}(*t* + 1), and other particles’ locations are updated based on the Eq. (5). The whole amending process can be expressed by the following equations

After above updates, the following criterions are executed:

(1) If **P**_{g} = **P**_{j}, which demonstrates that the random location **X**_{j} is the global best location. In this situation the *j*th particle will not be updated on the basis of Eq. (5) and the algorithm will randomly generate a location **X**_{j} in the searching space at the next iteration. At the same time, the velocities and locations of other particles are updated according to Eqs. (5) and (2) after **P**_{g} and **P**_{j} are updated, respectively.

(2) If **P**_{g} ≠ **P**_{j} and **P**_{g} has not been updated, which indicates that there is no better global best location found, compared with the last iteration, then all the particles’ the velocities and locations are updated based on the Eqs. (5) and (2), respectively.

(3) If **P**_{g} ≠ **P**_{j} and **P**_{g} has been updated, which demonstrates that there is a location of *k*th particle (*k* ≠ *j*) meets the requirement *k*th particle is the global best location and it is better than the global best location at the last iteration. In this situation, the *k*th particle stops evolving and it is used for storing global optima. Other particles’ velocities and locations of other particles are updated according to Eqs. (5) and (2) after **P**_{g} and **P**_{j} are updated, respectively.

Through the above analysis we can find that there is at least one particle’s location reaches the global best location at a particular iteration, which indicates that at least one particle is randomly generated at each iteration. Therefore, SPSO algorithm has been proved with strong global search ability.

### 2.4. Differential evolution PSO algorithm

Differential Evolution (DE) algorithm adopts simple differential operation among potential solutions to produce new candidate solution, which is a parallel, direct and stochastic searching technique. It was first proposed for solving Chebyshev polynomials and global optimization problems over continuous spaces by Storn and Price in 1995 [32]. Taking a cue from DE, the mutation operation is introduced into PSO algorithm to overcome the drawback of trapping in local optima, which is called DEPSO.

In DEPSO algorithm, the differential evolution operator is introduced to increase the diversity of particle swarms which is defined as

where *χ* is the differential mutation operator which controls the magnification of differential variation *M*] and they are not equal to the index *i*. Hence, the location evolution of *i*th particle can be expressed as

where *C* represents a preset constant value which satisfies *C* ≤ *F*_{min}, *F*_{min} indicates the minimum objective function value at the current iteration. The differential evolution term can force the particle to change if only *C* ≠ *F*_{min}, which can effectively prevent PSO algorithm from falling into local optima.

However, the location obtained by mutation operation maybe a worse result which will cause a bad influence on the search of other particles. In order to make sure the rapidity and stability of DEPSO algorithm, the following judgment should be executed before the location is updated:

(1) If the fitness value of the new location is better than the fitness value of the earlier location, which demonstrates the mutation is successful. Then the location of *i*th particle is updated according to Eq. (8).

(2) If the new fitness value is worse than before, which indicates this mutation is failed. Then the location of *i*th particle is updated according to Eq. (2) and the mutation operation of *i*th will continue at the next iteration until the mutation is successful.

In addition, there is a significant difference between DEPSO and basic PSO that the velocities of particles are not limited in the searching process, which can increase the convergence rate.

## 3. Simulation test of PSO algorithms

In order to test the performance of the above PSO algorithms, four benchmark optimization functions are used for verification whose details are shown in **Table 2**. The parameters in PSO algorithms are set as follows: the number of particle swarm population is set as *t*_{max}=1000 or the objective function value is less than 10^{-30}. The test results of the four functions with dimension *n*=2 are shown in **Figure 3**. **Table 3** lists the average objective function values and the corresponding mean squared errors for these four benchmark functions with 1000 independent runs. As shown, both the SPSO and DEPSO algorithms achieve better performance than the standard PSO algorithm in terms of computational accuracy.

## 4. Inverse geometry design of two-dimensional radiative enclosures

### 4.1. Description of the inverse geometry design problem

Considering a radiative equilibrium problem in two-dimensional irregular enclosures filled with participating media whose schematic diagram is shown in **Figure 4**. The curve EF represents the design surface and the design purpose is to produce a uniform distribution of radiative heat flux on the design surface. The bottom surface AD is the heating surface which can considered as the radiative heat source and its temperature is fixed as *T*_{S}. The two side surfaces AB and CD are cold with temperatures of 0K.

In order to optimize the geometry shape of the design surface to meet the specified requirement, the objective function (being equal to the fitness function in PSO algorithms) is defined as the square residuals between the estimated and average dimensionless radiative heat flux values which can be expressed as

where *N* is the number of the computational node on the design surface, and *Q*_{wi} and *Q*_{av} are the dimensionless radiative heat flux of the *i*th node and average value on the design surface, respectively. The iteration stop criterion of PSO is defined as

where *ε* is a small positive value. The smaller the value of *ε* is, the better the homogenization degree will be.

To evaluate the optimization results, the relative error is defined as

### 4.2. Akima cubic interpolation

The design surface is discretized into a series of control points, and Akima cubic interpolation is used to approximate the geometry shape of the surface in optimization process. Akima interpolation is formulated by a cubic polynomial between two control points. First, the prerequisite is introduced as [33]

where *g*_{k} is the slope of the curve at the position *x*_{k}, which can be defined as [33]

and the function *l*_{k} can be expressed as [33]

At the endpoints, the value of the function *l*_{k} can be defined as [33]

If the above equations are satisfied, then the cubic polynomial in the subinterval [

where *C*_{1}, *C*_{2}, *C*_{3}, and *C*_{4} are polynomial coefficients, which can be calculated as [33]

In order to test the performance of the Akima cubic interpolation, three interpolation cases are applied, in which six specified points are used as control points whose coordinates are shown in **Table 4**. The curves in **Figure 5** indicate the Akima cubic interpolation can successfully applied for geometry shape fitting.

### 4.3. Discrete ordinate method with a body-fitted coordinate

The boundary shape changes with the optimization of the design surface and the computational domain must be re-meshed, which greatly increase the solving difficulty of the inverse geometry design problems. In addition, the forward radiative heat transfer problem cannot be precisely solved by the normal numerical method. For the purpose of fitting the irregular boundary shape, the DOM with a body-fitted coordinate system is adopted to solve the RTE. For the participating media, the forward can be written as [34]

which is an integro-differential type, where *I* is the function of position *s* and direction **ŝ**. *β*_{e}, *κ*_{a}, and *κ*_{s} are the extinction, absorption, and scattering coefficients of media, respectively. The *Φ*(**ŝ**_{i}, **ŝ**) is the scattering phase function between incoming direction **ŝ**_{i} and scattering direction **ŝ**, which can be defined as *a*=0, *a*=1 and *a*=-1 for the isotropic scattering, the forward scattering and the backward scattering of media. The forward radiative heat transfer problem in the irregular enclosures is solved by using the computationally feasible DOM with a body-fitted coordinate system in this research. The 2D RTE discretized by the DOM can be expressed as [34, 35]

where *α*^{m} and *β*^{m} are the direction cosines of the discrete direction *m*. *s* refers to the spatial position, *w* is the quadrature weight. The radiative boundary condition can be directly imposed as follows [34]

where *ε*_{w} is the emissivity of boundaries, *I*_{b,w} is radiative intensity of blackbody boundaries, **n**_{w} represents the unit normal vector on the boundary, and **s** denotes the direction of radiative transfer.

The Jacobian matrix is used for coordinate transformation to solve the radiative heat transfer in irregular enclosures

The spatially discretized RTE with a body-fitted coordinate system can be expressed as

(22) |

where *P* represents the central node of the control volume. The subscripts *e*, *w*, *n*, and *s* represent the eastern, western, northern and southern boundaries around *P*, respectively.

The step scheme is applied to solve the above equations, and Eq. (22) can be expressed as

where the subscripts *E*, *W*, *N*, and *S* represent the central node of eastern, western, northern, and southern control volumes around control volume *P*, and

(24) |

Eq. (24) can be expressed in matrix form

where **A** represents the five-diagonal non-symmetric coefficient matrix, **Ψ** represents the vector that consists of the variables *I*^{m} at grid nodes, and **B** represents the vector that consists of the variables *b*^{m} on the right side of Eq. (23). The conjugate gradients stabilized (CGSTAB) method is adopted to solve the final discretized RTE because of its stability and fast convergence rate [35].

Consider a non-radiative equilibrium problem in the two-dimensional irregular enclosure filled with participating media. The four boundaries are cold surfaces whose temperature is 0K and all the boundaries are assumed as blockbody. The temperature of media is set as *T*_{g}, which can be considered as the heat source. The computational grids are meshed as 10 × 10 and the schematic in the computational domain is shown in **Figure 6**. The radiative heat transfer problem for different absorption coefficients in the radiative enclosure is solved by the body-fitted DOM and the retrieval results are compared with the results obtained by FVM in Ref. [36] which is shown in **Figure 7**. The curves show that the retrieval results achieve good agreements with FVM, which verified the accuracy and reliability of the computational program for solving radiative problems in irregular enclosures.

### 4.4. Inverse design results and discussions

The inverse geometry design model in Section 4.1 is considered and the standard PSO algorithm is abbreviated as PSO if there is no special instruction. The initial shape of the enclosure is rectangular, whose size is set as ^{-1} and ^{-1}, respectively. The parameters in the PSO algorithms are set as the same as the ones in test cases of Section 3. The stopping criteria are set as follows: (1) the objective function value is less than 10^{-7}, (2) the iteration times reach a maximum of 1000 and (3) the fitness value remains unchanged in consequent 100 iterations. The initial and final geometry shape of the design surfaces and their dimensionless radiative heat flux distributions are shown in **Figure 8(a)** and **8(b)**, respectively. The curves demonstrate that a uniform distribution of radiative heat flux on the design surface is obtained by means of PSO algorithms. **Figure 8(c)** shows the relative error distributions of the dimensionless radiative heat flux on the design surface. The average relative errors are 0.0310%, 0.0234% and 0.0279%, and the maximum relative errors are 0.0728%, 0.0540% and 0.0807% for PSO, SPSO and DEPSO algorithms, respectively. Overall, the retrieved results clearly reveal that the specified requirement of producing a uniform radiative heat flux distribution on the design surface can be obtained by PSO algorithm.

Algorithms | CPU time (s) | Fitness values | Average relative error (%) | Maximum relative error (%) |
---|---|---|---|---|

PSO | 40037 | 7.28×10^{-8} | 0.0313 | 0.0737 |

SPSO | 31572 | 4.05×10^{-8} | 0.0244 | 0.0518 |

DEPSO | 36944 | 5.62×10^{-8} | 0.0275 | 0.0669 |

In view of the random characteristic of intelligent algorithms, all the tests are repeated 50 trials to compare the performance of PSO algorithms. **Table 5** shows the comparison of the results obtained by PSO, SPSO and DEPSO algorithms. It can be found that all the PSO algorithms have reached the special design requirement and both SPSO and DEPSO achieve better performance than the initial PSO in terms of computational accuracy and efficiency.

In order to enhance the computational efficiency of inverse geometry design problems, the effects of corresponding parameters are investigated in this study. For the fact that the design surface is discretized into a series of control points in the inverse design process, the number of control point has a direct impact on the inverse geometry design results. The radiative physical parameters of media are kept as the same as the above typical example and the numbers of control point are set as *N*_{d}=1, 3, 5, and 7, respectively. The SPSO algorithm is adopted as the inverse design method and the initial and optimized geometry shape and dimensionless radiative heat flux on the design surface are shown in **Figure 9,** respectively. The curves show that the special design requirement is satisfied under the conditions that *N*_{d}=3, 5 and 7, whereas the homogenization degree is relatively poor in the case that *N*_{d}=1. **Table 6** lists the iteration numbers, fitness values and relatives errors for different numbers of control point. It can be found that both the iteration number and the relative error are smallest when *N*_{d}=3, which is because too many control points will decrease the sensitivity of radiative heat flux to the shape changing of design surface and one control point cannot provide enough necessary information for the geometry shape of the boundary [2, 9].

The physical properties of the media have an important influence on the energy transfer and then affect the optimization results of radiative enclosures. The effects of the extinction coefficient and the scattering albedo on the inverse design results are studied here. The scattering albedo of media is fixed as **Figure 10**. It can be seen that three PSO algorithms achieve similar boundaries and the geometry shapes are significantly different for different extinction coefficients. The extinction coefficient of media is fixed as **Figure 11**. It can be seen that the designed boundaries are very close to each other for different scattering albedos. **Tables 7** and **8** list the detailed inverse design results for different extinction coefficients and scattering albedos, respectively. As shown, both the optimized geometry the retrieval results for different extinction coefficients are significantly different, and both the height of the final designed surface and the dimensionless radiative heat flux decrease with the extinction coefficient increases. However, the retrieval results for different scattering albedo are close. Therefore, the scattering albedo of the media has little influence on the inverse geometry design results when the extinction coefficient is determined.

The scattering of media will affect the original transfer direction of radiative heat or energy, so the scattering characteristic of media also has influence on the inverse geometry design of radiative enclosures. The extinction coefficient and scattering albedo of the media are set as **Figure 12**. It can be found that the optimal geometry shapes of the design surface are close. **Table 9** plots the inverse geometry design results for different scattering characteristics. As shown, the values of dimensionless radiative heat flux on the design surface and relative errors are significantly different. Therefore, the scattering characteristics of media mainly affect the radiative heat flux on the boundaries and have no obvious effect on the final geometry shape of inverse radiative design problems.

## 5. Conclusions

In this chapter, the basic theoretical principles of PSO algorithm is introduced in detail and three kinds of PSO algorithms—standard PSO, SPSO and DEPSO—are applied for solving the inverse geometry design problem of a two-dimensional radiative enclosure filled with participating media. The design purpose is to produce a uniform distribution of radiative heat flux on the designed surface. The design surface is discretized into a series of control points and the Akima cubic interpolation is used to approximate the geometry shape of the boundary. The radiative heat transfer problem in the irregular enclosures is solved by the DOM with a body-fitted coordinate system. The pre-required radiative heat flux distribution is satisfied by optimizing the positions of control points based on the PSO algorithms. The retrieval results show that PSO algorithms can be successfully applied to solve inverse geometry design problems and SPSO achieves the best performance on computational time. Meanwhile, the scattering albedo and scattering characteristics of media have little effect on the geometry shape of the design surface. To improve the computational efficiency, the number of control points is recommended as *N*_{d}=3.