Next Article in Journal
Real-Time Detection and Correction of Abnormal Errors in GNSS Observations on Smartphones
Previous Article in Journal
Multi-Stage Feature Fusion of Multispectral and SAR Satellite Images for Seasonal Crop-Type Mapping at Regional Scale Using an Adapted 3D U-Net Model
Previous Article in Special Issue
Landsat-8/9 Atmospheric Correction Reliability Using Scene Statistics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Non-Dominated Sorting Genetic Algorithm II (NSGA2)-Based Parameter Optimization of the MSMGWB Model Used in Remote Infrared Sensing Prediction for Hot Combustion Gas Plume

School of Energy and Power Engineering, Beihang University, Beijing 102206, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(17), 3116; https://doi.org/10.3390/rs16173116
Submission received: 27 June 2024 / Revised: 12 August 2024 / Accepted: 20 August 2024 / Published: 23 August 2024

Abstract

:
The Multi-Scale Multi-Group Wide-Band (MSMGWB) model was used to calculate radiative transfer in strongly non-isothermal and inhomogeneous media such as the remote infrared sensing of aircraft exhaust system and jet plume scenario. In this work, the reference temperature was introduced into the model as an independent variable for each spectral subinterval group. Then, to deal with the exceedingly vast parameter sample space (i.e., the combination of spectral subinterval grouping results, reference temperatures and Gaussian quadrature schemes), an MSMGWB model’s parameter optimization process superior to the exhaustive approach employed in previous studies was established, which was consisted of the Non-dominated Sorting Genetic Algorithm II method (NSGA2) and an iterative scan method. Through a series of 0-D test cases and two real 3-D remote infrared imaging results of an aircraft exhaust system, it was observed that the MSMGWB model established and optimiazed in current work demonstrated notable improvements in both accuracy and computational efficiency.

Graphical Abstract

1. Introduction

The remote infrared detection of aircraft is extensively applied in domains such as military and aerospace [1]. Among them, the solid walls and jets of the aircraft exhaust system, which have relatively high temperatures, are one of the main sources of infrared signals for detection [2]. The calculation and prediction of their remote infrared radiation possess significant engineering significance. The calculation of remote infrared radiation mainly encounters two issues: computational efficiency and accuracy [3]. Regarding the computational efficiency, as the gas radiation absorption spectrum differs from that of solids and is discontinuous, the calculation of the precise radiation transfer equation necessitates the adoption of a line-by-line (LBL) calculation approach [4]. However, this approach demands a considerable amount of computing capacity, severely restricting its application in engineering issues. Hence, several spectral band models have been developed and employed in recent years, such as k-distribution models [5,6,7], l-distribution models [8], statistical spectral band models [9,10,11], etc. From the aspect of accuracy, as the gas absorption coefficient oscillates intensely along the spectrum and changes unpredictably with the variation of thermodynamic states, these models are challenged by the radiation transfer calculation in strongly non-isothermal and inhomogeneous media.
For the given remote infrared radiation detection within the atmospheric radiation window band that does not require detailed spectral radiation characteristics, the wide-band model possesses the advantage in terms of computational effort in contrast to the narrow-band model [3]. The early wide-band models, such as the box model [12] and the exponential wide-band model [13,14], directly obtain the average parameters within the spectral band. They are unable to capture and utilize the detailed spectral line information, frequently resulting in an average error of over 20% thus mainly concentrating on qualitative research. Yin employed the k-distribution model based on the principle of rearrangement of absorption coefficients to deal with the wide-band issue [15]. They made use of the spectral information to obtain the cumulative k-distribution function and decreased the average error to below 20%. The k-distribution model possesses two major advantages [16,17,18]: Firstly, it exhibits remarkable computational accuracy for radiation transfer in homogeneous gases and inhomogeneous gases that fulfill the correlated-k (CK) assumption, and the computational cost is merely a small portion of that of the LBL approach; Secondly, it is compatible with the methods of discretizing and solving the radiation transfer equation (RTE) in both differential and integral forms. However, the principal challenge in employing the aforementioned k-distribution model to solve the radiance issue such as the remote infrared imaging of the hot jet from an aircraft nozzle using hydrocarbon fuel is that the substantial disparity in the molar ratio and temperature of the radiative participating components between the combustion gas and the atmosphere can severely diminish the CK characteristics [19,20,21]. Studies have shown that the presence of “hot lines” in the absorption spectra of radiatively participating components exacerbates the disruption of CK characteristics due to temperature non-uniformity [22]. Current solutions to this issue can be categorized into the Multiple Line Group (MLG) [23,24,25] method and the Spectral Mapping Method (SMM) [26,27,28]. Both approaches involve dividing the concerned absorption spectra (SMM) or absorption lines (MLG) into several subsets to ensure that the absorption spectra of each subset maintain high CK characteristics across various thermodynamic states. The CK characteristic disruption caused by non-uniformity in component molar ratios is due to differences in the absorption spectra of the radiative participating components. Existing solutions include methods based on joint distribution functions [29], multiple integration [30] and various convolution [31]. Pal [32] incorporated both the convolution method [33] and the SMM method in the Full-Spectrum k-distribution (FSK) model, and the established Multi-Scale Multi-Group Full-Spectrum k-distribution (MSMGFSK) model can cope with the two CK characteristic failure mechanisms at the same time. Hu refined the same method to develop the MSMGWB model for wide-band applications, which was successfully applied to remote infrared sensing calculations [21]. However, any of the aforementioned methods for addressing the mechanisms that degrade CK characteristic come with increased computational demands, when using both the MS and MG methods to address two CK characteristic degradation issues. Therefore, both the MSMGFSK model and the MSMGWB model must be optimized to achieve the best balance between computational accuracy and cost.
Previous works have indicated that the random grouping initialization has a significant impact on the division results of the wavenumber subintervals, thus further affects the performance of the corresponding MSMGWB model [21]. Hu established an automatic optimization scheme to identify the optimal combination within the grouping results of the wavenumber subintervals and the Gaussian quadrature scheme [34]. Wang enhanced the scheme by augmenting the quantity of cases for the optimization objective function and the selection scope of Gaussian quadrature schemes, thereby obtaining model parameters with superior performance [35]. However, all the schemes of the aforementioned optimization platforms employ the exhaustive approach. As the optimization sample space increases, this approach has been incapable of obtaining optimization results within a time period measured in years.
The research contents of this study are as follows: In contrast to the previous work where the same reference temperature was employed among RTEs corresponding to each spectral subinterval group, the reference temperature within each group of the MSMGWB model was regarded as an independent variable for the first time, giving rise to three factors that require matching and optimization: the Gaussian quadrature points quantity, the reference temperature, and the spectral subinterval grouping results. The NSGA2 genetic algorithm was employed to obtain the optimal Gaussian quadrature points quantity and reference temperature at given grouping results [36,37]. An iterative scan method was proposed to determine the most suitable grouping combination result, then the optimization results of three factors were obtained. Finally, the calculation results and costs of the MSMGWB model optimized in this work were compared with other models in a series of 0-D cases, as well as two 3-D cases where the parameters of the combustion gas and ambient air are completely different.

2. MSMGWB Model

The MSMGWB model was proposed to predict the remote infrared radiation emitted by hot combustion gas jets. The RTE without scattering gain can be written as
d I η d s = ( k η ( ϕ ) k s η ( ϕ ) ) I b η k η ( ϕ ) I η I Δ η = η L η U I η d η
where I η is the spectral radiance, I b η the spectral radiance of blackbody, and s the location through the radiation transfer path. k η and k s η are the spectral extinction and scattering coefficient at wavenumber η , respectively. In this work, k η = k s η + n = 1 4 k n η , and k s η = k 4 s η , where n ranges from 1 to 4, corresponding to water vapor, carbon dioxide, carbon monoxide, and aerosol particles, respectively. ϕ = ϕ ( T , p , x , N ) is the gas and particle thermodynamic state at a certain space location, including temperature, pressure, component concentration, and particle number density, which determines the local k η and k s η . η L , η U , and Δ η are the lower wavenumber, upper wavenumber, wavenumber range of current infrared bands user concerned, respectively. In this study, k 1 ~ 3 η is calculated by LBL based on HITEMP2010 database for gas components [38], while k 4 η and k s η by Mie scattering theory for aerosol [39,40].
In hot combustion gas remote detection scenario, water vapor and carbon dioxide are treated as participation in both emission and absorption process, while carbon monoxide and aerosols only contribute to the attenuation part, because they exist only in the atmosphere rather than in turbo gas. Besides, the emission and scattering gain from atmosphere are not considered here since they are commonly treated as background radiation and calculated separately. Therefore, Equation (1) can be simplified into [19,41]:
d I n η d s = k n η I b η ( T ) k η I n η I η = n = 1 2 I n η k 3 η = k 4 η = k 4 s η = 0 , in combustion gas I b η = 0 , in atmosphere
where I η is calculated by emission of water vapor and carbon dioxide while absorption and scattering of all four components.
Next, to establish and later enhance the relation between the cumulative distribution function of n-thcomponents g n ( k n , T 0 ) and the total extinction coefficient k of the mixed medium, the equation is given as
0 g n k n * g n , T 0 , T k n g n , T d g n = 1 I b , Δ η ( T ) η L η U k n η H k k η I b η ( T ) d η g n k n , T = 1 I b , Δ η ( T ) η L η U H k n k n η I b η ( T ) d η
where I b , Δ η ( T ) = η L η U I b η ( T ) d η , k n is the absorption coefficient of n-th components, H the Heaviside step function, and T 0 the reference temperature. Both k η and k are monotonically increasing with g n ( T 0 ) by Equation (3). Consequently, as derived in [34] Equation (2) is transformed based on the k-distribution theory to
d I n g d s = k n * ( g n , T 0 ) a n ( k n * , T , T 0 ) I b , Δ η ( T ) k ( g n , T 0 ) I n g I Δ η = n = 1 2 I n = n = 1 2 0 1 I n g d g n ( T 0 )
where k n * is a pseudo absorption coefficient of n-th components, and a n ( k n * , T , T 0 ) = d g n ( k n * , T ) d g n ( k n * , T 0 ) is the non-gray stretch factor.
When the inhomogeneous medium contains only one participating component, it is obvious that Equation (4) is quivalent to Equation (1) if the extinction spectra is correlated, i.e., g n ( ϕ , k n η , T 0 ) = F ( k n η ( ϕ r e f ) ) | η , ϕ n 1 where F is a monotonically increasing function, ϕ r e f is a reference thermodynamic state. Therefore, Equations (1) and (4) are equivalent only when equation
g n ( ϕ , k η , T 0 ) = F k η ( ϕ r e f ) | η , ϕ
holds for cases involving two or more participating components in inhomogeneous media [34,35]. Generally, g n ( ϕ , k n , T 0 ) is not a constant for each η under various ϕ , introducing variations between the calculation results derived from Equations (1) and (4). This discrepancy is especially noticeable in scenarios with pronounced ϕ inhomogeneities.
The critical aspect of the MSMGWB model involves partitioning the wavenumber subintervals that constitute the absorption spectrum of the n-th component into M distinct subsets. The primary objective is to ensure that Equation (5) is satisfied within each subset of wavenumber subintervals to the greatest extent possible under the thermodynamic states of interest. Consequently, Equations (2)–(4) can be reformulated as
d I n m η d s = k n m η I b η ( T ) k η I n m η I η = n = 1 2 m = 1 M n I n m η k 3 η = k 4 η = k 4 s η = 0 , in combustion gas I b η = 0 , in atmosphere
0 g n m k n m * g n m , T n m 0 , T k n m g n m , T d g n m = 1 I b , Δ η ( T ) η L η U k n m η H k k η I b η ( T ) d η g n m k n m , T = 1 I b , Δ η ( T ) η L η U H k n m k n m η I b η ( T ) d η
d I n m g d s = k n m * ( g n m , T n m 0 ) a n m ( k n m * , T , T n m 0 ) k ( g n m , T n m 0 ) I n m g I Δ η = n = 1 2 m = 1 M n I n m = n = 1 2 m = 1 M n 0 1 I n m g d g n m ( T n m 0 ) a n m ( k n m * , T , T n m 0 ) = d g n m ( k n m * , T ) d g n m ( k n m * , T n m 0 )
It is different from the previous version MSMGWB model that an independent reference temperature was employed here, thus each group of each component has its own reference temperature T n m 0 and Gaussian quadrature scheme. M 1 and M 2 are not constrained to be identical; they can take different values since the grouping of k 1 η and k 2 η is independent. The width of the wavenumber subintervals must be smaller than that in LBL, which was set at 0.005 cm 1 in this study. The grouping method is described as follows:
  • Step i. Select the representative thermodynamic states encountered during the computation of atmospheric transmission characteristics for radiation emitted by hydrocarbon fuel combustion gases, as listed in Table 1.
  • Step ii. Stochastically assign an initial group membership for each wavenumber subinterval within the absorption spectrum of the water vapor.
  • Step iii. Compute M 1 group tags, which are given by
    ξ n m η = i = 1 5 i = 6 9 l n 2 k η ( ϕ i ) k ( ϕ i , g n m ( ϕ i , k η , T 0 ) , T 0 )
    for each wavenumber subinterval, followed by their redistribution to the group corresponding to the minimum ξ n m η .
  • Step iv. Iterate Step iii until fewer than 0.1% of the total k 1 η numbers change their group assignment.
  • Step v. Repeat Steps ii~iv for each k 2 η .
It is obvious that ξ n m η = 0 | η  m- t h g r o u p   o f  n- t h   c o m p o n e n t means Equation (5) is true within m-th group of n-th component, i.e., g n m ( ϕ , k η , T n m 0 ) = F k η ( ϕ r e f ) | η n m ϕ , where the η n m is the wavenumber subinterval of the m-th group of the n-th component. In practice, the integration in the middle formula of Equation (8) is calculated by Gauss-Legendre quadrature, to give
d I n m g ¯ , q ¯ d s = A ¯ n m , q ¯ ( T ) I b , Δ η ( T ) k ( g ¯ n m , q ¯ , T n m 0 ) I n m g ¯ , q ¯ I Δ η = n = 1 2 m = 1 M n I n m = n = 1 2 m = 1 M n q ¯ = 1 N n m g I n m g ¯ , q ¯ I n m g ¯ , q ¯ = g ˜ , m , q ¯ 1 g ˜ n m , q ¯ I n m g d g n m ( T n m 0 ) A ¯ n m , q ¯ ( T ) = 1 I b , Δ η ( T ) η U η L k n m η I b η ( T ) H k n m ( g ˜ n m , q ¯ , T n m 0 ) k n m η H k n m η k n m ( g ˜ n m , q ¯ 1 , T n m 0 ) 0 = g ˜ n m , 0 < g ¯ n m , 1 < g ˜ n m , 1 < g ¯ n m , 2 < g ˜ n m , q ¯ 1 < g ¯ n m , q ¯ < g ˜ n m , q ¯ < < g ˜ n m , N n m g 1 < g ¯ n m , N n m g < g ˜ n m , N n m g = g n m ( , T n m 0 )
where the N n m g is the total number of g ¯ n m , q ¯ ( T n m 0 ) . Both g ¯ n m , q ¯ g n m ( , T n m 0 ) and g ˜ n m , q ¯ g n m ( , T n m 0 ) are constants that solely denpend on the chosen Gauss quadrature schemes. It is worth noting that in contrast to the equation presented in [35] for the MSMGWB model, the sole disparity in the current study is that the reference temperature employed in each group is mutually independent.
The boundary solution of gray walls involving diffuse emitted and reflected radiation is transformed from Equation (8) as
I n m g ¯ , q ¯ , w = ε k ϕ w , g ¯ n m , q ¯ , T w A ¯ n m , q ¯ ϕ w , T w I b , Δ η T w + 1 ε π 0 4 π I n m g ¯ , q ¯ ( s ) max n w · s , 0 d Ω
where ϕ w = ( T w , p w , x w , N w ) represents the temperature and pressure at wall, gas component mole fractions and particle number densities of the medium adjacent to the wall, respectively. ε is the gray wall emissivity, s the (unit) direction vector of the incident radiation ( I n m g ( s ) ), and Ω the solid angle corresponding to the s , n w the (unit) normal vector of the wall.

3. Parameter Optimization Method

The iteration process for group division of wavenumber subintervals, as mentioned earlier, resembles the standard k-means clustering algorithm. Like other methods of the k-means family, a distance needs to be defined as the basis for measuring the relationship between sample point and cluster center, and the distance in this work is defined by Equation (9). Consequently, akin to the standard k-means clustering, the final results of group division which directly determine the performance of MSMGWB model are significantly impacted by its initialization (Step ii in Section 2). Besides, the quadrature accuracy of Equation (8) may not exhibit a monotonic increase with increasing N n m g , especially when the k - g n m function contain stair-like zones, as depicted in Figure 1.
Furthermore, the reference temperature value can also substantially impact correlated k-distribution models as reported in [34] and Figure 1. To identify the optimal combination of these three factors, we formulated an error function to gauge the performance of the MSMGWB model as
f e r r = j = 1 56 e r r o r j , max 0.08 + 10 ( e r r o r j , max 0.08 ) 0.48 + 100 ( e r r o r j , max 0.12 ) e r r o r j , max 0.08 0.08 < e r r o r j , max < 0.12 e r r o r j , max 0.12 e r r o r j , max = max I Δ η , j , MSMGWB ( L a ) I Δ η , j , LBL ( L a ) I Δ η , j , LBL ( L a )
where j is the 56 0-D cases number (see in Appendix A), and it is obvious that the larger the e r r o r j , m a x is, the more punishment (represented as slope) of f e r r gains. These 0-D cases is selected by typical scenarios of combustion gas and ambient atmosphere to comprehensively evaluate the performance of the MSMGWB model [34]. Figure 2 shows the effects of the wavenumber subinterval grouping results, the Gaussian quadrature points quantity, and the reference temperature on the f e r r , respectively, when the other two factors are fixed at their optimal values. It is elucidated that all three factors significantly impact f e r r . Besides, since we are pursuing the minimum value of f e r r , the left part of the figure is more noteworthy. The red line shows lower probability density but big slope at small f e r r area than the other two lines, i.e., the optimal Gaussian quadrature scheme is more difficult to be discovered than the other two factors.
In this section, the NSGA2 method will be employed to optimize the Gaussian quadrature points quantity and reference temperature for each group within an arbitrarily given water vapor and carbon dioxide grouping results combination. The optimization is based on the dual objective of quantitative accuracy and efficiency of the MSMGWB model. Concurrently, an iterative scheme driven by the NSGA2 optimization results will be established to select the optimal combination of water vapor and carbon dioxide grouping schemes.

3.1. Genotype Modeling

The genes chain consists of quadrature points quantity and reference temperature from each group. In 3~5 μ m band, for instance, the spectra of H 2 O and CO 2 are divided into 5 groups and 10 groups, respectively. Each group has its own Gaussian quadrature points quantity and reference temperature, therefore the length of the genes chain is 30. The quadrature points quantity range of each group is prescribed from 3 to 9 for feasible computational cost, and the reference temperature range is constrained from 100 K to 2500 K represented by 9 typical points. The variables in dark background color of the top chain in Figure 3 stand for quadrature points quantity of each group while the variables in light background color were the number of the corresponding reference temperature. Here we defined
f N = n = 1 2 m = 1 M n N n m g
which is directly proportional to the total Gaussian quadrature point quantity while solving RTE. Now, a certain gene chain is corresponding to certain parameters of the MSMGWB model, which also yields certain f e r r and f N .

3.2. Non-Dominated Genetic Algorithm Process

The non-dominated genetic algorithm process is depicted in Figure 4 and basically consists of three steps as initializing, generating offspring, and ending iteration.

3.2.1. Initializing

Randomly generate an initial value for each gene in the genotype of every individual in the population.

3.2.2. Generating Offspring

As Figure 3 shows, generating offspring contains selection, crossover and mutation process. Firstly, conduct non-dominated sorting for all individuals in the parent population. An individual A is considered to dominate individual B only when both objective function values f e r r and f N of individual A are smaller than those of individual B. It is evident that each individual may dominate others and, in turn, be dominated by different individuals. The number of times an individual is dominated is referred to as its rank. A rank of 0 indicates that the individual is not dominated by any other individuals, and all individuals with a rank of 0 form the Pareto front of the current solution set.
Secondly, select parent individuals with the current rank of 0, duplicate them directly into the new offspring, and decrease the rank values of the individuals they dominate by 1. After this process a new generation of the current Pareto front is identified and duplicated into the new offspring. Repeat this operation until the number of parent individuals directly copied into the offspring reaches 50% of the entire population. Individuals within the same rank are first reordered then prioritized based on their crowding distance, given to
d i = , i { 1 , N r } d i = ( f e r r , i + 1 f e r r , i 1 ) ( f e r r , m a x f e r r , m i n ) + ( f N , i + 1 f N , i 1 ) ( f N , m a x f N , m i n ) , i { 1 , N r }
where 1 and N r represent the first and the last individual after ascending reordering of f e r r and f N , respectively. A larger crowding distance implies higher genetic diversity in the solutions after the individual is added to the offspring, facilitating the expansion of the search space and reducing the risk of the iterative process getting trapped in local optima.
Finally, Another half of offspring are generated by crossover and mutation. Two parent individuals are randomly selected from the copied 50% advantageous population obtained in the second step, then they crossover by the probability P c . During the crossover process, two gene location points are randomly chosen, as depicted in Figure 3, to execute a two-point crossover and generate a new individual. In the absence of crossover by the probability P c , one of the parent individuals is directly duplicated to create a new individual. Due to the stochastic crossover point selection, overlapping points can lead to single-point crossover. Moreover when both points are situated at the head or tail location of the genes chain, it reverts to no crossover, involving a direct duplication of the parent individual. After crossover, the individual randomly mutate by the number of mutation points and the probability P m to finally generated a new individual into the offspring, as shown in Figure 5.
In this study, P c is fixed at 0.8, while the P m varies with the number of stagnant iterations of updating best solutions ( N S I ). When the N S I is less than 20, P m is set to 0.05, and the number of mutation points is set to 1, while 0.2 and 2 as N S I no less than 20. Note that when the number of mutation points is set to 2, the mutation process will occur with an equal probability (50%–50%) between single-point mutation and two-points mutation. Upon each update of the best solutions, the N S I is reset to zero, which ensures the P m stay at low value to avoid impeding convergence. When the N S I becomes high, it often indicates that the iteration process is nearing its conclusion or has reached a local optimum. In such scenarios, implementing a higher number of mutation points and an increased mutation probability aids in breaking out of local optima, thereby expanding the search space.

3.2.3. Iteration Ending

The iteration process is finished when the N S I reaches to 50 for it was proved by several repeated tests that 50 is large enough for the scale of this study to retain a convergence result.

3.2.4. Iteration Results

Figure 6a shows the fronts of a certain H 2 O - CO 2 grouping results combination. It is observed that the spacing between fronts is exceedingly thin, corresponding to a convergent and dense result by genetic algorithm. Therefore, the smallest error function value of the first front, defined as f e r r 0 , was served as the benchmark for assessing the iteration progress. Convergence process of f e r r 0 for ten random grouping results combinations in 20,000 population size is depicted in Figure 6b. All of them began with a large f e r r 0 for randomly initializing, descend rapidly, and eventually finish at a small stable f e r r 0 .
To assess the accuracy of the genetic algorithm, we compared the results obtained by the genetic algorithm and exhaustive search in the 7.7~9.7 μ m band which contains 11 H 2 O groups and 2 CO 2 groups. Each group has 3 options for quadrature points quantity (3/5/7), and 2 options for reference temperatures (300 K/1900 K), resulting in a total of 13,060,694,016 possible combinations. As shown in Figure 7, the genetic algorithm produced results identical to those of the exhaustive search with a population size of 1000, while it costs five orders of magnitude less computational time than the exhaustive search.

3.3. Grouping Result Selection Based on an Iterative Scan Method

Due to the stochastic process in Step ii of wavenumber subintervals group division, it is nearly impossible to achieve good enough grouping results combinations between H 2 O and CO 2 in few grouping experiments. Therefore, we produced about 10 to 10,000 grouping results for both H 2 O and CO 2 , then filter them for good combination by iteration method, as shown in Table 2 (S for the number of grouping results, subscript 1 and 2 for the H 2 O and CO 2 respectively). Besides, for each group, the quadrature points quantity option is set to 7 (3 to 9), and the reference temperature is set to 7 (100 K to 2500 K). As a result, the combination quantity of quadrature points quantity, reference temperature, and grouping results is exceedingly large so that exhaustive search is totally incompetent. Two methods were proposed below to address the excessive computational complexity issue, reducing the computational burden from quadratic growth with respect to the number of grouping results to linear growth, and finally collaborated with the non-dominated genetic algorithm to solve the optimization.

3.3.1. Iterative Scan Method

To illustrate the problem in hand, 100 grouping results of H 2 O and 400 grouping results for CO 2 were selected in 3~5 μ m band, resulting in a total of 40,000 grouping combinations. Each combination’s f e r r 0 were computed using the genetic algorithm with a population size of 40,000 individuals per generation. The results, as shown in Figure 8, depict smaller f e r r 0 in blue and larger f e r r 0 in red.
Four profile lines are randomly chosen for both H 2 O and CO 2 . These lines depict the f e r r 0 of the local grouping result combination traversing all CO 2 ( H 2 O ) grouping results while keeping the H 2 O ( CO 2 ) grouping result fixed. In most cases, if CO 2 grouping result serial number ( i C O 2 ) at a certain position, for example i C O 2 = 10 as indicated by the purple dashed line in the figure, the inequality relation f e r r 0 i H 2 O f e r r 0 i H 2 O holds, then for any value of i C O 2 , f e r r 0 i H 2 O f e r r 0 i H 2 O holds. The same is true for lines with i H 2 O = Const . Therefore, once the profile lines with low f e r r 0 values for both H 2 O and CO 2 are identified, the corresponding grouping results combinations are highly likely to contain the combination with the globally lowest f e r r 0 . This probability increases as the number of selected low f e r r 0 value profile lines grows. The scanning approach is as follows:
  • Step i. Randomly select an i x , H 2 O 1 , combine it with all i C O 2 and conduct the genetic algorithm iteration, arrange them in ascending order based on the f e r r 0 , and identify the i 1 , C O 2 1 that results in the minimum f e r r 0 .
  • Step ii. Select the i 1 , C O 2 1 from step i, combine it with all i H 2 O , and conduct the genetic algorithm iteration, arrange them in ascending order based on the f e r r 0 , and identify the i 1 , H 2 O 2 that results in the minimum f e r r 0 .
  • Step iii. Select the i 1 , H 2 O 2 from step ii, combine it with all i C O 2 and conduct the genetic algorithm iteration, and record the first p serial numbers with the smallest f e r r 0 as a set I C O 2 = { i 1 , C O 2 3 , i 2 , C O 2 3 , , i p , C O 2 3 } .
  • Step iv. Select the i 1 , C O 2 3 in I C O 2 from step iii, combine it with all i H 2 O and conduct the genetic algorithm iteration, and record the first q serial numbers with the smallest f e r r 0 as a set I H 2 O = { i 1 , H 2 O 4 , i 2 , H 2 O 4 , , i q , H 2 O 4 } .
  • Step v. Combine the set I C O 2 and I H 2 O , then find out the optimal grouping result combination from the total p × q combinations.
where the superscript of i (e.g., i 1 , H 2 O 2 ) represents the step number in which the value is generated, and the first item of the subscript represents the order sorted by f e r r 0 .
Considering the stochastic selection of i x , H 2 O 1 in step i, the scanning approach introduces an additional iteration process to ensure the stability of the results. Clearly, with an increaseing number of iterations, both the results and stability improve, but the computational cost increases simultaneously. Therefore, we compared 4 different iteration approaches, evaluating the results when the total groups results for H 2 O and CO 2 are taken in sets of 100 and 400, respectively. Each plan differs in the way I C O 2 and I H 2 O are obtained, while the remaining steps are the same. Plan A directly obtains I C O 2 and I H 2 O through steps i and ii only; Plan B follows the whole steps as described above; Plan C involves one more steps iii and iv than Plan B; Plan D repeats steps iii and iv multiple times, and filters out the grouping results corresponding to the higher 50% of f e r r 0 at each step, until the remaining number of grouping results decreases to predefined numbers of p and q. As shown in Figure 9, blue cells represent scanned sample points, green cells represent the currently scanning sample point, and red cells represent the sample point with the lowest f e r r 0 in the current scanning row (column), where the smaller f e r r 0 the higher saturation of red color. Orange cells represent the finally selected sample points, i.e., the global optimal sample points obtained by the algorithm. Black cells represent excluded sample points by Plan D.
The performance of the 4 plans is demonstrated in Figure 10. The f e r r results was obtained by averaging 100,000 repeated experiments for stability. The number of sample points is determined by the iterative scan method itself and the sizes of the sets I C O 2 and I H 2 O . Figure 9 illustrates the execution process of the four plans. It can be observed that Plan A, due to the randomness of Step i, cannot consistently obtain results with lower f e r r 0 , and it descends slowly as the number of sample points increases. Plan C, involving multiple executions of steps iii and iv, leads to an increase in computational complexity, but quickly approaches the global optimum as the number of sample points increases. Plan D exhibits a good descending trend with the early increase in sample points, but the descent slows in the later stages, as some good grouping results are excluded prematurely. Plan B demonstrates results close to Plan C with a relatively small number of computation sample points. Therefore, it is adopted as a compromise between performance and efficiency.

3.3.2. Population Size Selection

Due to the enormous genotype sample space in this study, the genetic algorithm requires a large population size to achieve satisfactory results. However, the computational complexity of the NSGA2 increases quadratically with the population size. Therefore, the computational workload remains significant even with the iterative scanning approach adopted. As a result, some improvements need to be made at the population size level. In this section, 20 grouping results of both H 2 O and CO 2 were randomly selected in 3~5 μ m band, forming 400 grouping combination sample points. The genetic algorithm was then employed with 9 different population sizes (100, 200, 500, 1000, 2000, 5000, 10,000, 20,000, 40,000), and the f e r r 0 obtained for each population was taken as the representative result.
As shown in Figure 11, the angular axis is divided into 400 parts representing 400 grouping results combinations, and the radial axis represents the population size. The results for a population size of 40,000 is taken as the baseline, and the contour values represent the ratio between the f e r r 0 of the model at current population size and the f e r r 0 at the baseline population size. It can be observed that as the population size increases, the ratio gradually decreases towards 1, indicating an improvement in algorithm performance. When the population size reaches 20,000, the improvement is not significant. Additionally, Figure 12 analyzes the results for each grouping results combination at population size of 5000 and 40,000. The Spearman correlation analysis shows a high correlation coefficient of 0.920 and a low p-value of 0.000, indicating a strong correlation between them. It is also evident that for the same grouping results combination, the calculation result for a population size of 40,000 is significantly better than that for a population size of 5000 in the majority of grouping results combinations (corresponding to the points under the dash line), especially at those with smaller f e r r 0 , which are exactly the grouping results combinations we are searching for. There are still few combination points above the dash line, however the most of them hold large f e r r 0 and will not be selected in the iterative scan method. Since the steps i~iv of Plan B place a higher emphasis on the order of the results rather than their specific values, it is reasonable and accurate enough to utilize small population size in these steps. Note that there is approximate 8 2 times of workload between 5000 and 40,000 methods. Therefore, the completed approach is: steps i~iv employ a population size of 5000, while step v employs a population size of 40,000. This approach significantly improves computational efficiency without deteriorating the results.

4. Results

4.1. 0-D Cases

Adopting the methods established in last section, we optimized the MSMGWB model in grouping results combination, Gaussian quadrature points quantity, and reference temperature, at 5 typical atmospheric infrared window bands (2~2.5 μ m, 3~5 μ m, 7.7~9.7 μ m, 8~14 μ m), and Figure 13a, Figure 14a, Figure 15a, Figure 16a and Figure 17a show the f N versus f e r r in each band. It is demonstrated that with increasing f N , the f e r r significantly reduces when f N is at low level, while the reduction of f e r r gradually becomes negligible at high f N level and even showing an upward trend.
When compared to the current MSMGWB model [35], the MSMGWB model optimized in this work achieved significantly better computational accuracy in the vast majority of the test cases as shown in Figure 13b, Figure 14b, Figure 15b, Figure 16b and Figure 17b and Table 3. The theoretical computational cost, i.e., the total number of solved RTE f N of the 2~2.5 μ m band, 3.7~4.8 μ m band, 3~5 μ m band, 7.7~9.7 μ m band, and 8~14 μ m band were decreased by 11.4%, 5.7%, 8.6%, 3.3%, and 11.6%, respectively. Meanwhile, the computational error ( f e r r ) were decreased by 52.8%, 61.4%, 62.4%, 62.3%, 57.3% at the 5 bands, respectively. The finally selected points at the Pareto fronts are a trade-off consideration of f N and f e r r , as shown by the red points in the Figure 13a, Figure 14a, Figure 15a, Figure 16a and Figure 17a.

4.2. Two High-Temperature Exhaust System 3-D Cases

Two 3-D cases of the high-temperature exhaust system cases were calculated. The MSMGWB model optimized in this work was used for remote infrared imaging and the error contours was obtained. These two 3-D cases have different exhaust system sizes, structures, working conditions, atmosphere parameters (see Table 4 and Table 5), and aerosol spectral extinction characteristics (see Figure 18).
For the large-sized case, as shown in Figure 19, the equivalent outlet diameter of the main nozzle of the supersonic exhaust system and the maximum diameter of the outer sleeve are 791 mm and 1220 mm, respectively. The inner wall of the main nozzle contains a heat shield and cooling structure, with a convergence angle of 35°. The working altitude of the nozzle is 7 km, the flight Mach number 0.7, the nozzle pressure ratio 2.38, the exhaust total temperature 900 K, and the surface emissivity of the material 0.8. The infrared detector is located on the ground and is 70 km away from the nozzle. The environmental atmosphere parameters are selected as the equatorial mode and the aerosol mode is selected as the 23-km visibility aerosol mode over the sea surface in the MODTRAN5 software. The infrared imaging resolution is 720 × 260 pixels, and the imaging direction forms an angle of 30° with the axis of the nozzle.
The number of computational grids for the internal and external flow fields of this exhaust system is approximately 920,000. The calculation results show that there is a series of shock waves and expansion waves in internal fluid field and jet plume, as well as a complex-shaped combustion gas and environment air mixing layer, which causes the non-uniformity of the gas temperature, pressure, and mole fraction of components. However, although the mole fractions of water vapor and carbon dioxide change in the fluid field, their ratio is close to 1 everywhere (Figure 20). Due to the inclusion of the gas film cooling structure, solid wall temperature is relatively low (Figure 21).
The calculation results of long-distance infrared imaging at 2~2.5 μ m, 3.7~4.8 μ m, 3~5 μ m, 7.7~9.7 μ m, and 8~14 μ m of this exhaust system and the calculation error distribution of the MSMGWB model are shown in Figure 22, with the latter based on the calculation results of the LBL model. Considering the computational cost of the LBL method, the radiative reflection of the solid wall is ignored.

4.3. Small-Sized High-Temperature Exhaust System 3-D Cases

For the small-sized case, as shown in Figure 23, the diameter of the throat and the maximum diameter of the outer sleeve are 128.6 mm and 152.4 mm, respectively. The high-temperature inner wall of the main nozzle has no cooling structure, with the convergence and divergence angles are 15.05° and 2.12° respectively. The working altitude of the nozzle is 0 km, the nozzle pressure ratio 10.1, the exhaust total temperature 1500 K, and the surface emissivity of the material 0.8. The infrared detector is located at an altitude of 7 km and is 50 km (for 2~2.5 μ m band), 70 km (for 7.7~9.7 μ m band), and 140 km (for other bands) away from the nozzle. The environmental atmosphere parameters are selected as the mid-latitudes winter mode and the aerosol mode is selected as the marine aerosol mode over the sea surface in the MODTRAN5.2.1.0 software.The infrared imaging resolution is 1000 × 400 pixels, and the imaging direction forms an angle of 30° with the axis of the nozzle.
The number of computational grids for the internal and external flow fields of this exhaust system is approximately 890,000. The calculation results show that compared with the large-sized exhaust system, although the size is much smaller, the temperature and pressure of the exhaust jet are higher. Since it has the same V-shaped trailing edge as the large-sized exhaust system, the shape of the mixing layer between the jet plume and the ambient air is complex (Figure 24). Due to the higher exhaust temperature and the absence of a gas film cooling structure, the temperature of its solid structure is much higher than that of the large-sized case, and the temperature distribution of the main nozzle shows a distinct 3-D heat conduction effect. Moreover, the outer sleeve is affected by the radiative heat transfer from the outer wall of the main nozzle, and its temperature is significantly higher than the surrounding ambient temperature (Figure 25).
The same as in the large-sized case, Figure 26 shows the calculation results of long-distance infrared imaging and error distribution.
As shown in Table 6, the max relative calculation error of both the large-sized case and small-sized case are no more than ± 13 % compared to the LBL results, which shows great accuracy with less computational cost than the previous MSMGWB model.

5. Discussion

The Pareto front results of the NSGA2 algorithm show that as the number of RTE equations solved increases, the calculation accuracy is further improved, but there is a marginal diminishing effect (Figure 13a, Figure 14a, Figure 15a, Figure 16a and Figure 17a). Note that we might not have deep enough search for the samples with big f N and overrate the f e r r 0 of them, because the central limit theorem leads few individuals with big f N at the initialization step of NSGA2, meanwhile a part of the driving force for the NSGA2 iterations is directed towards generating individuals with smaller f N . However, due to the compromise considerations of f N and f e r r 0 , the samples with big f N were not considered in practice for their disadvantage in computational efficiency.
Even though this study made progress in the modeling and parameter optimization of the MSMGWB model, several issues warrant further exploration. Firstly, In addition to the four main substances adopted in this work, there are some other radiative absorption components in the atmosphere (such as methane, nitrous oxide, etc.). More accurate calculations need to take all these components into account. Secondly, the number of groups for water vapor and carbon dioxide will also affect the calculation accuracy and efficiency of the model. After adding the remaining atmospheric components in the second point, this parameter needs to be updated. In the future work, we will commence from the aforementioned contents and keep enhancing the performance and universality of the MSMGWB model.

6. Conclusions

In this study, compared with the MSMGWB models where each spectral subinterval group shares the same reference temperature [34,35,42], we proposed to optimize the reference temperature as an independent parameter within each group. To deal with the increase in the number of model parameter combinations, the NSGA2 genetic algorithm and an iterative scan method were introduced, and it was demonstrated that the above algorithms can obtain the optimal model parameters combination. The calculation results of 0-D and 3-D cases show that the MSMGWB model optimized based on the above algorithms has significantly improved computational efficiency and accuracy in 5 atmospheric infrared window wavebands compared to the present MSMGWB model.
Other conclusions were drawn as follows:
The Gaussian quadrature scheme has a greater influence on the computational accuracy of the MSMGWB model than the spectrum grouping results and reference temperature; A population size of 20,000 is sufficient for parameter optimization of the MSMGWB model based on the NSGA2 algorithm; At the scanning and screening stage of spectrum grouping results combinations, 5000 of population is enough and can greatly reduce the calculation time; The spectrum grouping result of water vapor and carbon dioxide has an approximately independent influence on the accuracy of the MSMGWB model. Therefore, the iterative scan method can be used to quickly obtain the optimal grouping results combination.

Author Contributions

Conceptualization, H.H. and Y.L.; methodology, Y.L.; software, Y.L.; validation, Y.L., H.H. and Q.W.; formal analysis, Y.L.; investigation, Y.L.; resources, Y.L.; data curation, Y.L.; writing—original draft preparation, Y.L.; writing—review and editing, Y.L., H.H. and Q.W.; visualization, Y.L.; supervision, Y.L.; project administration, H.H.; funding acquisition, H.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by J2019-III-0009-0053 National Science and Technology Major Project, and the APC was funded by National Science and Technology Major Project.

Data Availability Statement

The source code of this work can be found at this website, https://github.com/1oyue/mdpi_data, accessed on 19 August 2024.

Acknowledgments

We thank the School of Energy and Power Engineering, Beihang University for its help in this project.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A

As shown in Figure A1, there are two types of radiative transfer paths in this study among 56 0-D cases. In the first type, radiance emits through a hot combustion gas with a transmission path length L h (which holds constant in the 0-D case calculation), then attenuates though the atmosphere under various transmission distances L c . In the second type, radiance emits from a black wall and is attenuated by the atmosphere after passing through a hot gas layer and a cold gas layer with constant path lengths of L h and L c , respectively. The detailed parameters of the cases are shown in Table A1.
Figure A1. Two types of radiative transfer paths, diagram of 56 0-D cases.
Figure A1. Two types of radiative transfer paths, diagram of 56 0-D cases.
Remotesensing 16 03116 g0a1
Table A1. Thermal dynamic states table of 56 0-D cases.
Table A1. Thermal dynamic states table of 56 0-D cases.
Hot Gas ParametersCold Gas ParametersEnvironmental Atmospheric Parameters
NumberLh [cm]p [atm]T [K] x H 2 O x C O 2 T w [K] L c [cm] p [atm] T [K] x H 2 O x C O 2 x C O L a [cm] p [atm] T [K] x H 2 O x C O 2 x C O
15018000.10.1 401308.150.05724 3.23 × 10 4 1.47 × 10 7
25016000.10.1 401288.150.00184 3.23 × 10 4 1.47 × 10 7
35018000.10.1 401288.150.00184 3.23 × 10 4 1.47 × 10 7
4801.513000.10.1 1000.9298.150.03226 3.23 × 10 4 1.47 × 10 7
580113000.10.1 1001298.150.03226 3.23 × 10 4 1.47 × 10 7
680213000.10.1 1000.8285.20.00959 3.23 × 10 4 1.39 × 10 7
780113000.10.1 1001298.150.03226 3.23 × 10 4 1.47 × 10 7
88014000.110.11 4013000.0323 3.23 × 10 4 1.47 × 10 7
980116000.110.11 10013000.0323 3.23 × 10 4 1.47 × 10 7
1080216000.110.11 1000.93000.0323 3.23 × 10 4 1.47 × 10 7
1180119000.120.12 10013000.0323 3.23 × 10 4 1.47 × 10 7
12801.519000.120.12 1000.93000.0323 3.23 × 10 4 1.47 × 10 7
13801.519000.120.12 1000.9294.20.0184 3.23 × 10 4 1.47 × 10 7
14701.618000.140.12 1000.83000.03 3.23 × 10 4 1.47 × 10 7
1560110500.10.1 801294.20.0184 3.23 × 10 4 1.47 × 10 7
1660210500.10.1 800.9294.20.0184 3.23 × 10 4 1.47 × 10 7
1760210500.10.1 801294.20.0184 3.23 × 10 4 1.47 × 10 7
1860210500.10.1 800.7279.20.00595 3.28 × 10 4 1.34 × 10 7
191000.4215000.080.08 2000.42254.70.00102 3.30 × 10 4 1.25 × 10 7
201000.17718000.10.1 2000.177215.8 8.01 × 10 6 3.30 × 10 4 6.38 × 10 8
2150115000.10.190015015500.050.05 1000.93000.0323 3.23 × 10 4 1.47 × 10 7
2250115000.10.145015015500.050.05 1000.93000.0323 3.23 × 10 4 1.47 × 10 7
2350215000.10.145015015500.050.05 1000.93000.0323 3.23 × 10 4 1.47 × 10 7
24502.517000.10.165030015500.050.05 2000.62273.20.0038 3.29 × 10 4 1.31 × 10 7
25501.517000.10.165030015500.050.05 2000.62273.20.0038 3.29 × 10 4 1.31 × 10 7
26502.517000.10.165010015500.050.05 2000.62273.20.0038 3.29 × 10 4 1.31 × 10 7
27502.514000.10.145030015500.050.05 2000.62273.20.0038 3.29 × 10 4 1.31 × 10 7
28500.718000.120.125001500.485500.10.1 2000.32241.7 4.13 × 10 4 3.30 × 10 4 1.09 × 10 7
295018000.10.1 401308.150.05724 3.23 × 10 4 1.47 × 10 7
305016000.10.1 401288.150.00184 3.23 × 10 4 1.47 × 10 7
315018000.10.1 401288.150.00184 3.23 × 10 4 1.47 × 10 7
3280113000.10.1 1001298.150.03226 3.23 × 10 4 1.47 × 10 7
3380113000.10.1 1001298.150.03226 3.23 × 10 4 1.47 × 10 7
3480113000.10.1 301298.150.03226 3.23 × 10 4 1.47 × 10 7
358014000.110.11 4013000.0323 3.23 × 10 4 1.47 × 10 7
3680216000.110.14 600.93000.0323 3.23 × 10 4 1.47 × 10 7
3721.618000.140.12 200.83000.024 3.23 × 10 4 1.47 × 10 7
3860110500.10.1 801294.20.0184 3.23 × 10 4 1.47 × 10 7
3960210500.10.1 800.9294.20.0184 3.23 × 10 4 1.47 × 10 7
401000.4215000.080.08 2000.42254.70.00102 3.30 × 10 4 1.25 × 10 7
411000.17718000.10.1 2000.177215.8 8.01 × 10 6 3.30 × 10 4 6.38 × 10 8
4250115000.10.1 15015500.050.05 1000.93000.0323 3.23 × 10 4 1.47 × 10 7
43150.86500.10.1 1000.8288.150.005 3.23 × 10 4 1.47 × 10 7
44517500.10.1 1000.8288.150.01 3.23 × 10 4 1.47 × 10 7
4550.59000.10.1 1200.5263.150.002 3.23 × 10 4 1.47 × 10 7
46100.55000.10.1 401293.150.015 3.23 × 10 4 1.47 × 10 7
471015500.120.1 1000.6273.150.004 3.23 × 10 4 1.47 × 10 7
481015000.10.12 8013000.012 3.23 × 10 4 1.47 × 10 7
491502.516000.10.1 1000.8288.150.015 3.23 × 10 4 1.47 × 10 7
50100.515000.130.1 1200.9293.150.02 3.23 × 10 4 1.47 × 10 7
5170114000.120.124003013000.03 3.4 × 10 4 1.47 × 10 7 1200.62600.0015 3.30 × 10 4 1.31 × 10 7
521500.617000.10.12600100.62600.0015 3.3 × 10 4 1.31 × 10 7 4013000.03 3.30 × 10 4 1.47 × 10 7
5380216000.110.14 1000.93000.0323 3.23 × 10 4 1.47 × 10 7
5420318000.120.121000401.514000.10.1 4013000.03 3.23 × 10 4 1.47 × 10 7
5530.86000.050.0545050.85000.020.02 3013000.02 3.23 × 10 4 1.47 × 10 7
5630.83500.050.054500.84004000.02 4012880.01 3.23 × 10 4 1.47 × 10 7

References

  1. Jiang, W.; Pan, H.; Wang, Y.; Li, Y.; Lin, Y.; Bi, F. A Multi-Level Cross-Attention Image Registration Method for Visible and Infrared Small Unmanned Aerial Vehicle Targets via Image Style Transfer. Remote Sens. 2024, 16, 2880. [Google Scholar] [CrossRef]
  2. Zhou, Y.; Wang, Q.; Li, T. A new model to simulate infrared radiation from an aircraft exhaust system. Chin. J. Aeronaut. 2017, 30, 651–662. [Google Scholar] [CrossRef]
  3. Modest, M.F. The Treatment of Nongray Properties in Radiative Heat Transfer: From Past to Present. J. Heat Transf.-Trans. ASME 2013, 135, 061801. [Google Scholar] [CrossRef]
  4. Hartmann, J.M.; Levi Di Leon, R.; Taine, J. Line-by-line and narrow-band statistical model calculations for H2O. J. Quant. Spectrosc. Radiat. Transf. 1984, 32, 119–127. [Google Scholar] [CrossRef]
  5. Soufiani, A.; Taine, J. High temperature gas radiative property parameters of statistical narrow-band model for H2O, CO2 and CO, and correlated-k model for H2O and CO2. Int. J. Heat Mass Transf. 1997, 40, 987–991. [Google Scholar] [CrossRef]
  6. Modest, M.F. Narrow-band and full-spectrum K-distributions for radiative heat transfer–correlated-k vs. scaling approximation. J. Quant. Spectrosc. Radiat. Transf. 2003, 76, 69–83. [Google Scholar] [CrossRef]
  7. Marin, O.; Buckius, R.O. Wide Band Correlated-k Approach to Thermal Radiative Transport in Nonhomogeneous Media. J. Heat Transf. 1997, 119, 719–729. [Google Scholar] [CrossRef]
  8. André, F. The L-distribution method for modeling non-gray absorption in uniform and non-uniform gaseous media. J. Quant. Spectrosc. Radiat. Transf. 2016, 179, 19–32. [Google Scholar] [CrossRef]
  9. Curtis, A. A statistical model for watervapour absorption. Q. J. R. Met. Soc. 1952, 78, 639–640. [Google Scholar]
  10. Godson, W.L. The computation of infrared transmission by atmospheric water vapor. J. Atmos. Sci. 1955, 12, 272–284. [Google Scholar] [CrossRef]
  11. Malkmus, W. Random Lorentz band model with exponential-tailed S-1 line-intensity distribution function. JOSA 1967, 57, 323–329. [Google Scholar] [CrossRef]
  12. Penner, S.; Landshoff, R. Quantitative molecular spectroscopy and gas emissivities. Phys. Today 1960, 13, 38–40. [Google Scholar] [CrossRef]
  13. Edwards, D.; Menard, W. Comparison of models for correlation of total band absorption. Appl. Opt. 1964, 3, 621–625. [Google Scholar] [CrossRef]
  14. Edwards, D.; Balakrishnan, A. Thermal radiation by combustion gases. Int. J. Heat Mass Transf. 1973, 16, 25–40. [Google Scholar] [CrossRef]
  15. Yin, X.M.; Liu, L.H.; Li, B.X. Numerical simulation of infrared signature emitted by liquid rocket plume using wide band K-distribution model. Adv. Mater. Res. 2012, 516, 41–53. [Google Scholar] [CrossRef]
  16. Modest, M.F.; Mazumder, S. Radiative Heat Transfer; Academic Press: Cambridge, MA, USA, 2021. [Google Scholar]
  17. Zhang, F.; Zhu, M.; Li, J.; Li, W.; Di, D.; Shi, Y.N.; Wu, K. Alternate Mapping Correlated k-Distribution Method for Infrared Radiative Transfer Forward Simulation. Remote Sens. 2019, 11, 994. [Google Scholar] [CrossRef]
  18. Wang, Q.; Li, Y.; Hu, H. Numerical simulation of remote infrared imaging for marine exhaust system using wide-band k-distribution model combined with a new grouping method for spectral sub-intervals. Infrared Phys. Technol. 2019, 96, 276–290. [Google Scholar] [CrossRef]
  19. Zhang, H.; Modest, M.F. A multi-scale full-spectrum correlated-k distribution for radiative heat transfer in inhomogeneous gas mixtures. J. Quant. Spectrosc. Radiat. Transf. 2002, 73, 349–360. [Google Scholar] [CrossRef]
  20. Rivière, P.; Soufiani, A.; Taine, J. Correlated-k fictitious gas model for H2O infrared radiation in the Voigt regime. J. Quant. Spectrosc. Radiat. Transf. 1995, 53, 335–346. [Google Scholar] [CrossRef]
  21. Hu, H.; Wang, Q. Improved spectral absorption coefficient grouping strategy of wide band k-distribution model used for calculation of infrared remote sensing signal of hot exhaust systems. J. Quant. Spectrosc. Radiat. Transf. 2018, 213, 17–34. [Google Scholar] [CrossRef]
  22. Zhang, H.; Modest, M.F. Scalable Multi-Group Full-Spectrum Correlated-k Distributions for Radiative Transfer Calculations. J. Heat Transf. 2003, 125, 454–461. [Google Scholar] [CrossRef]
  23. Khodyko, Y.V.; Kurskov, A.A.; Antiporovich, N.V. Multigroup method for the calculation of selective IR radiation transfer in nonhomogeneous media. J. Appl. Spectrosc. 1986, 44, 284–288. [Google Scholar] [CrossRef]
  24. Pierrot, L.; Rivière, P.; Soufiani, A.; Taine, J. A fictitious-gas-based absorption distribution function global model for radiative transfer in hot gases. J. Quant. Spectrosc. Radiat. Transf. 1999, 62, 609–624. [Google Scholar] [CrossRef]
  25. Soufiani, A.; André, F.; Taine, J. A fictitious-gas based statistical narrow-band model for IR long-range sensing of H2O at high temperature. J. Quant. Spectrosc. Radiat. Transf. 2002, 73, 339–347. [Google Scholar] [CrossRef]
  26. West, R.; Crisp, D.; Chen, L. Mapping transformations for broadband atmospheric radiation calculations. J. Quant. Spectrosc. Radiat. Transf. 1990, 43, 191–199. [Google Scholar] [CrossRef]
  27. André, F.; Hou, L.; Roger, M.; Vaillon, R. The multispectral gas radiation modeling: A new theoretical framework based on a multidimensional approach to K-Distrib. Methods. J. Quant. Spectrosc. Radiat. Transf. 2014, 147, 178–195. [Google Scholar] [CrossRef]
  28. André, F.; Solovjov, V.P.; Webb, B.W.; Galtier, M.; Dubuisson, P. Spectral mapping method based on intervals of comonotonicity for modelling of radiative transfer in non-uniform gaseous media. J. Quant. Spectrosc. Radiat. Transf. 2019, 229, 33–39. [Google Scholar] [CrossRef]
  29. Pierrot, L.; Soufiani, A.; Taine, J. Accuracy of narrow-band and global models for radiative transfer in H2O, CO2, and H2O-CO2 mixtures at high temperature. J. Quant. Spectrosc. Radiat. Transf. 1999, 62, 523–548. [Google Scholar] [CrossRef]
  30. Shi, G.; Xu, N.; Wang, B.; Dai, T.; Zhao, J. An improved treatment of overlapping absorption bands based on the correlated k distribution model for thermal infrared radiative transfer calculations. J. Quant. Spectrosc. Radiat. Transf. 2009, 110, 435–451. [Google Scholar] [CrossRef]
  31. Zhang, H.; Nakajima, T.; Shi, G.; Suzuki, T.; Imasu, R. An optimal approach to overlapping bands with correlated k distribution method and its application to radiative calculations. J. Geophys. Res. Atmos. 2003, 108. [Google Scholar] [CrossRef]
  32. Pal, G.; Modest, M.F. A New Hybrid Full-Spectrum Correlated k-Distribution Method for Radiative Transfer Calculations in Nonhomogeneous Gas Mixtures. In Proceedings of the ASME/JSME 2007 Thermal Engineering Heat Transfer Summer Conference Collocated with the ASME 2007 InterPACK Conference, Vancouver, BC, Canada, 8–12 July 2007. [Google Scholar] [CrossRef]
  33. Wang, L.; Modest, M.F. A Hybrid Multi-Scale Full-Spectrum k-Distribution Method for Radiative Transfer in Inhomogeneous Gas Mixtures. In Proceedings of the ASME International Mechanical Engineering Congress and Exposition, Orlando, FL, USA, 5–11 November 2005; Volume 42223, pp. 175–182. [Google Scholar]
  34. Hu, H.; Li, Y.; Wei, Z.; Wang, Q. Parameter optimization for MSMGWB model used to calculate infrared remote sensing signals emitted by hot combustion gases of hydrocarbon fuel. J. Quant. Spectrosc. Radiat. Transf. 2020, 249, 107003. [Google Scholar] [CrossRef]
  35. Wang, Q.; Hao, J.; Hu, H.; Li, Y.; Wu, Y.; Xie, Y.; Deng, H. Optimization of the MSMGWB models used to predict remote infrared signals of jet engine in various spectral intervals. Infrared Phys. Technol. 2024, 140, 105403. [Google Scholar] [CrossRef]
  36. Deb, K.; Pratap, A.; Agarwal, S.; Meyarivan, T. A fast and elitist multiobjective genetic algorithm: NSGA-II. IEEE Trans. Evol. Comput. 2002, 6, 182–197. [Google Scholar] [CrossRef]
  37. Washington, R.; Garmatyuk, D.; Mudaliar, S.; Narayanan, R.M. Many-Objective RadarCom Signal Design via NSGA-II Genetic Algorithm Implementation and Simulation Analysis. Remote Sens. 2022, 14, 3787. [Google Scholar] [CrossRef]
  38. Rothman, L.S.; Gordon, I.E.; Barber, R.J.; Dothe, H.; Gamache, R.R.; Goldman, A.; Perevalov, V.I.; Tashkun, S.A.; Tennyson, J. HITEMP, the high-temperature molecular spectroscopic database. J. Quant. Spectrosc. Radiat. Transf. 2010, 111, 2139–2150. [Google Scholar] [CrossRef]
  39. Volz, F.E. Infrared Refractive Index of Atmospheric Aerosol Substances. Appl. Opt. 1972, 11, 755–759. [Google Scholar] [CrossRef] [PubMed]
  40. Wiscombe, W.J. Improved Mie scattering algorithms. Appl. Opt. 1980, 19, 1505–1509. [Google Scholar] [CrossRef] [PubMed]
  41. Pal, G.; Modest, M.F. A Narrow Band-Based Multiscale Multigroup Full-Spectrum k-Distribution Method for Radiative Transfer in Nonhomogeneous Gas-Soot Mixtures. J. Heat Transf. 2009, 132, 023307. [Google Scholar] [CrossRef]
  42. Hu, H.; Wang, Q. Improved spectral absorption-coefficient grouping strategies in radiation heat transfer calculations for combustion gases with pressure and temperature inhomogeneity. Numer. Heat Transf. Part B Fundam. 2019, 75, 178–197. [Google Scholar] [CrossRef]
Figure 1. The relationship between k and g n m of a group in different reference temperatures and thermodynamic states at 8~14 μ m band.
Figure 1. The relationship between k and g n m of a group in different reference temperatures and thermodynamic states at 8~14 μ m band.
Remotesensing 16 03116 g001
Figure 2. Relationship between probability density of objective function value and three critical factors (Gaussian quadrature point quantity, reference temperature, and wavenumber subinterval grouping).
Figure 2. Relationship between probability density of objective function value and three critical factors (Gaussian quadrature point quantity, reference temperature, and wavenumber subinterval grouping).
Remotesensing 16 03116 g002
Figure 3. Genotype and crossover process diagram.
Figure 3. Genotype and crossover process diagram.
Remotesensing 16 03116 g003
Figure 4. NSGA2 algorithm workflow diagram.
Figure 4. NSGA2 algorithm workflow diagram.
Remotesensing 16 03116 g004
Figure 5. Offspring generation workflow diagram.
Figure 5. Offspring generation workflow diagram.
Remotesensing 16 03116 g005
Figure 6. Convergence results of the NSGA2 method: (a) the foremost 10 Pareto front results, (b) convergence iteration process of 10 random grouping strategy combinations.
Figure 6. Convergence results of the NSGA2 method: (a) the foremost 10 Pareto front results, (b) convergence iteration process of 10 random grouping strategy combinations.
Remotesensing 16 03116 g006
Figure 7. f e r r 0 results between exhaustive search method and NSGA2 method.
Figure 7. f e r r 0 results between exhaustive search method and NSGA2 method.
Remotesensing 16 03116 g007
Figure 8. The f e r r 0 results among 100 H 2 O and 400 CO 2 grouping strategy combinations.
Figure 8. The f e r r 0 results among 100 H 2 O and 400 CO 2 grouping strategy combinations.
Remotesensing 16 03116 g008
Figure 9. Diagram of 4 iterative scan method process plans.
Figure 9. Diagram of 4 iterative scan method process plans.
Remotesensing 16 03116 g009
Figure 10. Convergence perfomance of 4 plans for scan iteration process.
Figure 10. Convergence perfomance of 4 plans for scan iteration process.
Remotesensing 16 03116 g010
Figure 11. Ratio of the f e r r 0 at the current sample population size to its corresponding baseline value.
Figure 11. Ratio of the f e r r 0 at the current sample population size to its corresponding baseline value.
Remotesensing 16 03116 g011
Figure 12. f e r r 0 results between the same grouping result combination in the NSGA2 model population sizes of 5000 and 40,000.
Figure 12. f e r r 0 results between the same grouping result combination in the NSGA2 model population sizes of 5000 and 40,000.
Remotesensing 16 03116 g012
Figure 13. Optimization results at 2~2.5 μ m band: (a) Pareto front, (b) e r r o r j , m a x in 56 0-D cases.
Figure 13. Optimization results at 2~2.5 μ m band: (a) Pareto front, (b) e r r o r j , m a x in 56 0-D cases.
Remotesensing 16 03116 g013
Figure 14. Optimization results at 3.7~4.8 μ m band: (a) Pareto front, (b) e r r o r j , m a x in 56 0-D cases.
Figure 14. Optimization results at 3.7~4.8 μ m band: (a) Pareto front, (b) e r r o r j , m a x in 56 0-D cases.
Remotesensing 16 03116 g014
Figure 15. Optimization results at 3~5 μ m band: (a) Pareto front, (b) e r r o r j , m a x in 56 0-D cases.
Figure 15. Optimization results at 3~5 μ m band: (a) Pareto front, (b) e r r o r j , m a x in 56 0-D cases.
Remotesensing 16 03116 g015
Figure 16. Optimization results at 7.7~9.7 μ m band: (a) Pareto front, (b) e r r o r j , m a x in 56 0-D cases.
Figure 16. Optimization results at 7.7~9.7 μ m band: (a) Pareto front, (b) e r r o r j , m a x in 56 0-D cases.
Remotesensing 16 03116 g016
Figure 17. Optimization results at 8~14 μ m band: (a) Pareto front, (b) e r r o r j , m a x in 56 0-D cases.
Figure 17. Optimization results at 8~14 μ m band: (a) Pareto front, (b) e r r o r j , m a x in 56 0-D cases.
Remotesensing 16 03116 g017
Figure 18. Aerosol spectral extinction coefficient at 0~7 km altitude and 2~14 μ m: (a) large-sized case, (b) small-sized case.
Figure 18. Aerosol spectral extinction coefficient at 0~7 km altitude and 2~14 μ m: (a) large-sized case, (b) small-sized case.
Remotesensing 16 03116 g018
Figure 19. Diagram of the Large-sized exhaust system with a cooling structure.
Figure 19. Diagram of the Large-sized exhaust system with a cooling structure.
Remotesensing 16 03116 g019
Figure 20. Distribution of temperature (T), pressure (p), carbon dioxide mass fraction ( y C O 2 ), and Mach number ( M a ) in the meridional and axial sections of the fluid field of the large-sized exhaust system with a cooling structure.
Figure 20. Distribution of temperature (T), pressure (p), carbon dioxide mass fraction ( y C O 2 ), and Mach number ( M a ) in the meridional and axial sections of the fluid field of the large-sized exhaust system with a cooling structure.
Remotesensing 16 03116 g020
Figure 21. Temperature (T) distribution of the solid part of the large-sized exhaust system with a cooling structure.
Figure 21. Temperature (T) distribution of the solid part of the large-sized exhaust system with a cooling structure.
Remotesensing 16 03116 g021
Figure 22. Remote infrared imaging of the large-sized exhaust system with a cooling structure in different atmospheric window bands (left), and the distribution of calculation errors of the optimized MSMGWB model (right), (a,b) 2~2.5 μ m band, (c,d) 3.7~4.8 μ m band, (e,f) 3~5 μ m band, (g,h) 7.7~9.7 μ m band, (i,j) 8~14 μ m band.
Figure 22. Remote infrared imaging of the large-sized exhaust system with a cooling structure in different atmospheric window bands (left), and the distribution of calculation errors of the optimized MSMGWB model (right), (a,b) 2~2.5 μ m band, (c,d) 3.7~4.8 μ m band, (e,f) 3~5 μ m band, (g,h) 7.7~9.7 μ m band, (i,j) 8~14 μ m band.
Remotesensing 16 03116 g022
Figure 23. Diagram of the small-sized exhaust system without a cooling structure.
Figure 23. Diagram of the small-sized exhaust system without a cooling structure.
Remotesensing 16 03116 g023
Figure 24. Distribution of temperature (T), pressure (p), carbon dioxide mass fraction ( y C O 2 ), and Mach number ( M a ) in the meridional and axial sections of the fluid field of the small-sized exhaust system without a cooling structure.
Figure 24. Distribution of temperature (T), pressure (p), carbon dioxide mass fraction ( y C O 2 ), and Mach number ( M a ) in the meridional and axial sections of the fluid field of the small-sized exhaust system without a cooling structure.
Remotesensing 16 03116 g024
Figure 25. Temperature (T) distribution of the major components of the small-sized exhaust system without a cooling structure.
Figure 25. Temperature (T) distribution of the major components of the small-sized exhaust system without a cooling structure.
Remotesensing 16 03116 g025
Figure 26. Remote infrared imaging of the small-sized exhaust system without a cooling structure in different atmospheric window bands (left) and the distribution of calculation errors of the optimized MSMGWB model (right), (a,b) 2~2.5 μ m band, (c,d) 3.7~4.8 μ m band, (e,f) 3~5 μ m band, (g,h) 7.7~9.7 μ m band, (i,j) 8~14 μ m band.
Figure 26. Remote infrared imaging of the small-sized exhaust system without a cooling structure in different atmospheric window bands (left) and the distribution of calculation errors of the optimized MSMGWB model (right), (a,b) 2~2.5 μ m band, (c,d) 3.7~4.8 μ m band, (e,f) 3~5 μ m band, (g,h) 7.7~9.7 μ m band, (i,j) 8~14 μ m band.
Remotesensing 16 03116 g026
Table 1. Representative thermodynamic states for aeroengine jet plume and atmosphere.
Table 1. Representative thermodynamic states for aeroengine jet plume and atmosphere.
ScenarioThermodynamic State T ( K ) x H 2 O x C O 2 x C O p ( a t m )
Aeroengine jet plume ϕ 1 19000.120.1202
ϕ 2 19000.120.1201
ϕ 3 15000.10.100.5
ϕ 4 9000.080.0801
ϕ 5 9000.080.0800.5
Atmosphere ϕ 6 3000.034 3.4 × 10 4 1.42 × 10 7 1
ϕ 7 3000.0068 3.4 × 10 4 1.42 × 10 7 1
ϕ 8 2930.02 3.4 × 10 4 1.42 × 10 7 0.9
ϕ 9 2630.002 3.4 × 10 4 1.42 × 10 7 0.5
Table 2. Grouping result quantity and wavenumber group quantity of five typical infrared remote sensing bands for H 2 O and CO 2 .
Table 2. Grouping result quantity and wavenumber group quantity of five typical infrared remote sensing bands for H 2 O and CO 2 .
Parameter2~2.5 μ m3.7~4.8 μ m3~5 μ m7.7~9.7 μ m8~14 μ m
S 1 319624711510,00010,000
S 2 72810,00010,000102172
M 1 15551110
M 2 51010210
Table 3. Comparison between the optimized MSMGWB model and other models.
Table 3. Comparison between the optimized MSMGWB model and other models.
f e r r Number of Solved RTEs/Transmissivities
Wave-bandMSMGWB-newMSMGWB in [35]SNBFGNBKDMSMGWB-newMSMGWB in [35]SNBFGNBKD
2~2.5 μ m8.1917.3554.211212.41091232723280
3.7~4.8 μ m5.1513.33130.951.0828721610,300
3~5 μ m2.105.59216.1111.9647033611,840
7.7~9.7 μ m6.3016.7224.231097.2596165550
8~14 μ m3.47.0112.861111.272951371730
Table 4. Atmosphere thermodynamic state parameters at 0~7 km altitude in the large-sized case.
Table 4. Atmosphere thermodynamic state parameters at 0~7 km altitude in the large-sized case.
Altitude [km]p [atm]T [K] x H 2 O x C O 2 x C O
0~10.947296.7 2.073 × 10 2 3.301 × 10 4 1.476 × 10 7
1~20.843290.7 1.925 × 10 2 3.300 × 10 4 1.426 × 10 7
2~30.750285.7 1.488 × 10 2 3.298 × 10 4 1.378 × 10 7
3~40.665280.4 6.619 × 10 3 3.294 × 10 4 1.331 × 10 7
4~50.588273.6 3.926 × 10 3 3.296 × 10 4 1.306 × 10 7
5~60.519267.0 2.761 × 10 3 3.298 × 10 4 1.295 × 10 7
6~70.456260.3 1.718 × 10 3 3.302 × 10 4 1.268 × 10 7
Table 5. Atmosphere thermodynamic state parameters at 0~7 km altitude in the small-sized case.
Table 5. Atmosphere thermodynamic state parameters at 0~7 km altitude in the small-sized case.
Altitude [km]p [atm]T [K] x H 2 O x C O 2 x C O
0~10.938258.1 1.350 × 10 3 3.300 × 10 4 1.481 × 10 7
1~20.822257.5 1.778 × 10 3 3.298 × 10 4 1.434 × 10 7
2~30.719254.3 1.629 × 10 3 3.295 × 10 4 1.377 × 10 7
3~40.628250.2 9.886 × 10 4 3.295 × 10 4 1.331 × 10 7
4~50.547244.3 6.202 × 10 4 3.293 × 10 4 1.307 × 10 7
5~60.475237.5 3.396 × 10 4 3.295 × 10 4 1.295 × 10 7
6~70.411230.7 1.946 × 10 4 3.302 × 10 4 1.268 × 10 7
Table 6. Max relative error of the optimized MSMGWB model in two 3-D cases at five wave-bands.
Table 6. Max relative error of the optimized MSMGWB model in two 3-D cases at five wave-bands.
Wave-BandMax Relative Error (Large-Sized Case)Max Relative Error (Small-Sized Case)
2~2.5 μ m−8.35/+9.95%−3.24/+10.41%
3.7~4.8 μ m−6.19/+10.19%−5.48/+12.17%
3~5 μ m−4.06/+3.78%−4.43/+7.79%
7.7~9.7 μ m−9.84/+4.86%−6.48/+0.04%
8~14 μ m−6.65/+5.56%−8.49/+2.32%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, Y.; Hu, H.; Wang, Q. Non-Dominated Sorting Genetic Algorithm II (NSGA2)-Based Parameter Optimization of the MSMGWB Model Used in Remote Infrared Sensing Prediction for Hot Combustion Gas Plume. Remote Sens. 2024, 16, 3116. https://doi.org/10.3390/rs16173116

AMA Style

Li Y, Hu H, Wang Q. Non-Dominated Sorting Genetic Algorithm II (NSGA2)-Based Parameter Optimization of the MSMGWB Model Used in Remote Infrared Sensing Prediction for Hot Combustion Gas Plume. Remote Sensing. 2024; 16(17):3116. https://doi.org/10.3390/rs16173116

Chicago/Turabian Style

Li, Yihan, Haiyang Hu, and Qiang Wang. 2024. "Non-Dominated Sorting Genetic Algorithm II (NSGA2)-Based Parameter Optimization of the MSMGWB Model Used in Remote Infrared Sensing Prediction for Hot Combustion Gas Plume" Remote Sensing 16, no. 17: 3116. https://doi.org/10.3390/rs16173116

APA Style

Li, Y., Hu, H., & Wang, Q. (2024). Non-Dominated Sorting Genetic Algorithm II (NSGA2)-Based Parameter Optimization of the MSMGWB Model Used in Remote Infrared Sensing Prediction for Hot Combustion Gas Plume. Remote Sensing, 16(17), 3116. https://doi.org/10.3390/rs16173116

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

Article Metrics

Back to TopTop