
Citation: Madhura Kulkarni-Chitnis, Leah Mitchell-Bush, Remmington Belford, Jenaye Robinson, Catherine A. Opere, Sunny E. Ohia, Ya Fatou N. Mbye. Interaction between hydrogen sulfide, nitric oxide, and carbon monoxide pathways in the bovine isolated retina[J]. AIMS Neuroscience, 2019, 6(3): 104-115. doi: 10.3934/Neuroscience.2019.3.104
[1] | Yasong Sun, Jiazi Zhao, Xinyu Li, Sida Li, Jing Ma, Xin Jing . Prediction of coupled radiative and conductive heat transfer in concentric cylinders with nonlinear anisotropic scattering medium by spectral collocation method. AIMS Energy, 2021, 9(3): 581-602. doi: 10.3934/energy.2021028 |
[2] | Hong-Wei Chen, Fu-Qiang Wang, Yang Li, Chang-Hua Lin, Xin-Lin Xia, He-Ping Tan . Numerical design of dual-scale foams to enhance radiation absorption. AIMS Energy, 2021, 9(4): 842-853. doi: 10.3934/energy.2021039 |
[3] | Chao-Jen Li, Peiwen Li, Kai Wang, Edgar Emir Molina . Survey of Properties of Key Single and Mixture Halide Salts for Potential Application as High Temperature Heat Transfer Fluids for Concentrated Solar Thermal Power Systems. AIMS Energy, 2014, 2(2): 133-157. doi: 10.3934/energy.2014.2.133 |
[4] | Hong-Yu Pan, Chuang Sun, Xue Chen . Transient thermal characteristics of infrared window coupled radiative transfer subjected to high heat flux. AIMS Energy, 2021, 9(5): 882-898. doi: 10.3934/energy.2021041 |
[5] | Caliot Cyril, Flamant Gilles . Pressurized Carbon Dioxide as Heat Transfer Fluid: In uence of Radiation on Turbulent Flow Characteristics in Pipe. AIMS Energy, 2014, 1(2): 172-182. doi: 10.3934/energy.2014.2.172 |
[6] | Matteo Moncecchi, Davide Falabretti, Marco Merlo . Regional energy planning based on distribution grid hosting capacity. AIMS Energy, 2019, 7(3): 264-284. doi: 10.3934/energy.2019.3.264 |
[7] | Khadim Ndiaye, Stéphane Ginestet, Martin Cyr . Thermal energy storage based on cementitious materials: A review. AIMS Energy, 2018, 6(1): 97-120. doi: 10.3934/energy.2018.1.97 |
[8] | Kokou Aménuvéla Toka, Yawovi Nougbléga, Komi Apélété Amou . Optimization of hybrid photovoltaic-thermal systems integrated into buildings: Impact of bi-fluid exchangers and filling gases on the thermal and electrical performances of solar cells. AIMS Energy, 2024, 12(5): 1075-1095. doi: 10.3934/energy.2024051 |
[9] | Muluken Biadgelegn Wollele, Abdulkadir Aman Hassen . Design and experimental investigation of solar cooker with thermal energy storage. AIMS Energy, 2019, 7(6): 957-970. doi: 10.3934/energy.2019.6.957 |
[10] | Gemma Graugés Graell, George Xydis . Solar Thermal in the Nordics. A Belated Boom for All or Not?. AIMS Energy, 2022, 10(1): 69-86. doi: 10.3934/energy.2022005 |
Abbreviations: ${A_{ij}}$: area of wall cell ($i, j$), m2; $e$: energy carried by sample Monte Carlo bundle, W; Eb(T): spectral blackbody emissive power at temperature T, W/(m2·m); $H$: total exercise times; $L$: length of the domain, m; M, N: lattice indices in square mesh corresponding to x, y position, respectively; n: refractive index of medium; ${N_{ij}}$: number of sample bundles emitted by cell ($i, j$); ${N_{ij \to k{\kern 1pt} l}}$: number of sample bundles emitted by cell ($i, j$) and absorbed by cell ($k, l$); $N_{ij}^*$: number of sample bundles emitted by cell ($i, j$) in the BDMC method; $N_{ij \to kl}^*$: number of sample bundles emitted by cell ($i, j$) and absorbed by cell ($k, l$)in the BDMC method; ${N_0}$: the number density of sampling bundles for a reference case; ${N_r}$: the number of sample bundles chosen for a reference volume cell ${V_r}$; $P$: performance parameter for Monte Carlo simulation; $R{D_{\lambda, ij \to kl}}$: radiative exchange factor of cell ($i, j$) to cell ($k, l$) for non-gray medium in the TMC method; ${\overline {RD} _{\lambda, ij \to kl}}$: modified form of $R{D_{\lambda, ij \to kl}}$; $R{D_{ij \to kl}}$: radiative exchange factor of cell ($i, j$) to cell ($k, l$) for gray medium in the TMC method; ${\overline {RD} _{ij \to kl}}$: modified form of $R{D_{ij \to kl}}$; $RD_{ij \to kl}^ * $: radiative exchange factor of cell ($i, j$) to cell ($k, l$) in the BDMC method; $\overline {RD} _{ij \to kl}^ * $: modified form of $RD_{ij \to kl}^ * $; ${s_n}$, ${\hat s_n}$: bundles trajectories and coupled bundles trajectories; $t$: time, s; $T$: absolute temperature, K; ${V_{ij}}$: volume of medium cell ($i, j$), m3; ${W_{ij}}$: bundle weight of cell ($i, j$); ${\delta _{ij \to \, kl}}$: reciprocity error of a couple of radiative exchange factor by the TMC method; $\delta _{ij \to \, kl}^ * $: reciprocity error of a couple of radiative exchange factor by the BDMC method; $\Delta P$: variation of performance parameter; ${\Delta _r}$: relative error of reciprocity; ${\varepsilon _\lambda }$: spectral emissivity; ${{\mathit{\Phi }} _{\lambda, ij}}$: spectral radiative energy emitted from cell ${V_{ij}}$ (or ${A_{ij}}$); ${{\mathit{\Phi }} _{\lambda, ij \to \, kl}}$: spectral radiative energy emitted from cell ${V_{ij}}$ (or ${A_{ij}}$), that is absorbed by cell ${V_{kl}}$ (or ${A_{kl}}$); ${\mathit{\Phi }} _{ij \leftrightarrow kl}^ * $: net radiative energy exchange, W; $\gamma $: standard deviation; ${\kappa _{a\lambda }}$: spectral absorption coefficient of medium; ${\kappa _e}$: extinction coefficient, m-1; $\sigma $: Stefan-Boltzmann constant, $\left({\sigma = 5.67 \times {{10}^{ - 8}} {\rm{W/{m^2}{K^4}}}} \right)$; $\tau $: optical thickness $(\tau = {\kappa _e}L)$; $\omega $: scattering albedo of medium. Subscripts, $a$: medium absorption; $h$: hth sampling; $i{\kern 1pt} {\kern 1pt}, \, {\kern 1pt} j$: order number of cell; r: reference cell; $w$: wall
Radiative heat transfer in participating medium is described by radiative transfer equation (RTE), which is an integro-differential equation. Many methods have been developed to solve the RTE, the common used numerical methods such as the discrete ordinate method (DOM) and the finite volume method (FVM), rely on both spatial discretization and angular discretization, which are difficult to solve the RTE with high accuracy. Monte Carlo (MC) method is a stochastic statistical method based on the physical processes, it has been applied to solve radiative heat transfer in various participating media [1,2,3,4,5,6]. Moreover, the results predicted by MC method can often be treated as benchmark solutions due to its high accuracy of solution [7,8].
However, for the statistical nature of the MC method, the high computational cost is still a considerable disadvantage of MC simulation. Many attempts have been made in order to improve the computational efficiency of the method [9]. For example, several sampling approaches were developed to improve the speed of convergence, such as the importance sampling method [10], the rejection sampling method [11], the differential sampling method [12] and the weight-equivalent sampling method [13]. In addition, parallel computing technique was introduced into the MC simulation to improve computation efficiency evidently [1,9]. Howell [9] analyzed the advantages of various programming strategies of the MC method for radiative heat transfer in absorbing and scattering medium. All of these improvements were mainly mathematical efforts and hardware improvement of computer, which are universal to the MC simulations.
It has been noticed that the physical feature of thermal radiation transfer can also provide possibility to improve the MC simulation of radiative heat transfer in participating medium. For example, Walters [14] developed a reverse method based on the reciprocity principle for radiative heat transfer in a generalized enclosure containing an absorbing, emitting and scattering medium, the reverse method was proved to be efficient. Cherkaoui et al. [15,16] developed a net exchange Monte Carlo method based on a net-exchange formulation, provided an efficient way of systematically fulfilling the reciprocity principle, the computing time was proved much smaller than the conventional Monte Carlo approach. Lataillade [17] and his cooperators applied the net exchange Monte Carlo approach for radiative heat transfer in optically thick medium with spectral dependent radiative properties. Eymet et al. [18] extended this method to absorbing, emitting, and scattering media. Tessé et al. [19] improved the forward Monte Carlo (FMC) method based on the reciprocity principle, the method was used for radiative transfer in real gases, and proved to be a better choice for optically thick or nearly isothermal media compared with the forward Monte Carlo method. In addition, the reverse Monte Carlo (RMC) methods, based on the reversibility of radiative transfer trajectory, has been developed to solve the radiative transfer in absorbing and scattering media [20,21,22,23,24]. Kovtanyuk et al. [25] presented a recursive algorithm based on modification of Monte Carlo method, the modified method was used to solve the coupled radiative and conductive heat transfer in an absorbing and scattering medium, and was proved to be more accurate. Soucasse et al. [26] proposed a Monte Carlo formulation for radiative transfer in quasi-isothermal media which consists in directly computing the difference between the actual radiative field and the equilibrium radiative field at the minimum temperature in the medium.
In the present study, a bidirectional Monte Carlo (BDMC) method based on reversibility of bundle trajectory and reciprocity of radiative energy exchange was developed to solve thermal radiation transfer in absorbing and scattering medium. Two types of sampling models for MC simulation were presented, namely the equivalent sampling and the weight sampling. The equivalent sampling was chosen for the uniform mesh while the weight sampling was more suitable to the non-uniform mesh. The bidirectional information of tracing a sampling bundle was utilized by the BDMC method, the solution precision or efficiency can be evidently improved. Radiative heat transfer in a two-dimensional rectangular domain with absorbing and scattering media was solved by the BDMC method and the TMC method, respectively. The radiative exchange factors and the temperature profiles were investigated, in addition, the performance parameter defined by Howell [9] was also calculated to evaluate the two MC methods.
Radiative heat transfer in a two-dimensional (2-D) rectangular domain with absorbing, emitting and/or scattering medium was investigated. Figure 1 shows the 2-D rectangular geometry as well as the coordinate system. The four walls were assumed to be diffuse and gray. Radiative properties, such as the absorption coefficient, the scattering coefficient were assumed to be constant. The rectangular domain was divided into M × N = MN cells, any of which was depicted by Vij (medium cell) or Aij (wall cell), wherein the subscripts I∈[1, M] and j∈[1, N]. Monte Carlo method has robust adaptability and can be extended to more complex cases.
In the present study, the radiative exchange factor $R{D_{\lambda, \; ij \to kl}}$ was introduced to decouple the solution of radiative transfer from that of the temperature profile. It was defined as a fraction of the spectral radiative energy emitted from cell Vij (or Aij), that is absorbed by cell Vkl (or Akl) [27,28].
$R{D_{\lambda , \;ij \to kl}} = {{{{\mathit{\Phi }} _{\lambda , \;ij \to kl}}} / {{{\mathit{\Phi }} _{\lambda , \;ij}}}}$ | (1) |
As Eq 1, where ${{\mathit{\Phi }} _{\lambda, ij}}$ is the spectral radiative energy emitted from cell Vij (or Aij), ${{\mathit{\Phi }} _{\lambda, ij \to kl}}$ is the spectral radiative energy emitted from cell Vij (or Aij), and absorbed by cell Vkl (or Akl), taking into account possible wall reflections. It is obvious that the conservation relation $\sum\limits_{k = 1, l = 1}^{M, N} {R{D_{\lambda, ij \to kl}}} = 1$ is tenable according to the definition of radiative exchange factors. Note that the radiative exchange factor depends only on the system geometry and the radiative properties distribution of the medium [29]. The reciprocity relation between a couple of radiative exchange factors was given by [30].
${({\kappa _{a\lambda }})_{ij}}{V_{ij}}R{D_{\lambda , ij \to kl}} = {({\kappa _{a\lambda }})_{kl}}{V_{kl}}R{D_{\lambda , kl \to ij}}$ | (2) |
or
$4{({\kappa _{a\lambda }})_{ij}}{V_{ij}}R{D_{\lambda , ij \to k{\kern 1pt} l}} = {({\varepsilon _\lambda })_{kl}}{A_{kl}}R{D_{\lambda , kl \to ij}}$ | (3) |
As Eqs 2 and 3, where, ${\kappa _{a\lambda }}$ is the spectral absorption coefficient of the medium, while ${\varepsilon _\lambda }$ is the spectral emissivity of the boundary wall. For convenience, the radiative exchange factors were usually modified as
${\overline {RD} _{\lambda , ij \to kl}} = 4{({\kappa _{a\lambda }})_{ij}}{V_{ij}}R{D_{\lambda , ij \to kl}}$ | (4) |
or
${\overline {RD} _{\lambda , ij \to kl}} = {({\varepsilon _\lambda })_{ij}}{A_{ij}}R{D_{\lambda , ij \to kl}}$ | (5) |
here, Eqs 4 and 5 are used for volume and wall cells, respectively. The reciprocity relation can be written as
${\overline {RD} _{\lambda , ij \to kl}} = {\overline {RD} _{\lambda , kl \to ij}}$ | (6) |
As Eq 6, after solving the radiative exchange factors, the net radiative energy exchange ${\mathit{\Phi }} _{ij \leftrightarrow kl}^ * $ from ${V_{ij}}$ to ${V_{kl}}$ can be calculated from
${\mathit{\Phi }} _{ij \leftrightarrow kl}^ * = \int\limits_0^\infty {{{\overline {RD} }_{\lambda , ij \to kl}}\left[ {{E_{b\lambda }}\left( {{T_{ij}}} \right) - {E_{b\lambda }}\left( {{T_{kl}}} \right)} \right]} \;d\lambda $ | (7) |
As Eq 7, where, is the Stefan-Boltzmann constant ( = 5.67×10-8 W/m2K4), Ebλ(T) is the spectral blackbody emissive power at temperature T.
In the MC simulation, ${\mathit{\Phi } _{\lambda, ij}}$ is represented by a large number of independent sampling bundles, the propagating process of each sampling bundle can be tracked and counted. For example, the number of total sampling bundles for ${{\mathit{\Phi }} _{\lambda, ij}}$ is ${N_{ij}}$, each bundle with the same energy of ${e_{ij}} = {{\mathit{\Phi }} _{\lambda, ij}}/{N_{ij}}$. If ${N_{ij \to kl}}$ bundles among them are finally absorbed by ${V_{kl}}$, then, ${{\mathit{\Phi }} _{\lambda, ij \to kl}} = {e_{ij}} \cdot {N_{ij \to kl}}$, the radiative exchange factor can be calculated from [31].
$R{D_{\lambda , ij \to kl}} = {N_{ij \to kl}}/{N_{ij}}$ | (8) |
As Eq 8, the number ${N_{ij \to kl}}$ can be counted from the MC simulation, however, only an approximate value of ${N_{ij \to kl}}$ can be obtained because of the pseudo randomicity [32] in the sampling process. In fact, the reciprocity correlation shown in Eq 6 cannot be strictly satisfied, and random errors always exist. In order to obtain more accurate results, one need to increase the sampling bundles ${N_{ij}}$, which results in the increasing computing cost.
In the TMC method, the propagation trajectories of ${N_{ij}}$ sampling bundles for cell ${V_{ij}}$ are tracked and counted to get the value of ${N_{ij \to kl}}$ for any cell ${V_{kl}}$. Similarly, the number ${N_{kl \to \, ij}}$ is obtained after tracing another ${N_{kl}}$ propagation trajectories of the sampling bundles for cell ${V_{kl}}$. Wherein, only the forward propagation information of a trajectory is used, the number of sampling bundles for cell ${V_{ij}}$ is ${N_{ij}}$.
According to the reversibility of light propagation, for a tracing trajectory from ${V_{kl}}$ to ${V_{ij}}$, a bundle from ${V_{ij}}$ can also propagate to ${V_{kl}}$ along the reverse direction of the tracing trajectory. The information of the reverse direction can be used in the MC simulation, the bidirectional information can be obtained from the forward tracing, which results in the development of the present BDMC method.
In order to simplify the description of the BDMC method, the formula in the following sections were given for gray medium with gray walls. The basic idea of BDMC is to calculate the radiative exchange factor by the reversibility of light beams. According to the reversibility of light beams, the transmission path of the bundles received by cells can also be used as the transmission path of the bundles transmitted by cells. In this way, the effective information of sampling bundles if fully utilized without increasing the number of sampling bundles and calculation amount. Theoretically, the statistical sampling bundles in calculation can be doubled. And it can better satisfy the relationship of the reciprocal of radiative exchange factor. Figure 1 shows the schematic diagram of radiative transfer in a rectangular domain with participating medium in the BDMC method, it can be noticed that there are ${N_{ij}}$ tracing trajectories departing from cell ${V_{ij}}$, and $\sum\limits_{k = 1, l = 1}^{M, N} {{N_{kl \to \, ij}}} $ tracing trajectories getting into cell ${V_{ij}}$ at the same time. According to the reversibility of bundle trajectories, there must be $\sum\limits_{k = 1, l = 1}^{M, N} {{N_{kl \to \, ij}}} $ tracing trajectories depart from cell ${V_{ij}}$ in the reverse directions. It is believed that the $\sum\limits_{k = 1, l = 1}^{M, N} {{N_{kl \to \, ij}}} $ trajectories are independent of the emitted ${N_{ij}}$ trajectories because of the randomness of the sampling process. Thus, the number of total sampling bundles for ${V_{ij}}$ can be counted anew as
$N_{ij}^ * = {N_{ij}} + \sum\limits_{k = 1, l = 1}^{M, N} {{N_{kl \to \, ij}}} $ | (9) |
while the number of bundles emitted by ${V_{ij}}$ and absorbed by ${V_{kl}}$ is counted anew as
$N_{ij \to kl}^ * = {N_{ij \to kl}} + {N_{kl \to ij}}$ | (10) |
then, the radiative exchange factor $RD_{ij \to kl}^ * $ in the BDMC method can be calculated by
$RD_{ij \to kl}^ * = \frac{{N_{ij \to kl}^ * }}{{N_{ij}^ * }} = \frac{{{N_{ij \to k{\kern 1pt} l}} + {N_{kl \to ij}}}}{{{N_{ij}} + \sum\limits_{k = 1, l = 1}^{M, N} {{N_{kl \to \, ij}}} }}$ | (11) |
where, $RD_{ij \to kl}^ * $ is used to denote the radiative exchange factor in the BDMC method. The information of the bidirectional trajectories has been used, as is shown in Eq 11, it indicates that the effective sampling bundles are increased, while the computing cost remains unchanged. Therefore, more accurate results can be predicted from Eq 11 than that from Eq 8 under the same computing cost.
It should be pointed out that Eqs 9 and 11 are valid only if the sampling bundles for cell Vij and any other cells are equivalent. According to the reversibility of light beams, the transmission path of the bundles received by cells can also be used as the transmission path of the bundles transmitted by cells. In this way, the effective information of sampling bundles if fully utilized without increasing the number of sampling bundles and calculation amount. In other words, the forward and reverse bundles in the BDMC method should have the same contribution for radiative exchange factors calculation. To use Eq 11, the equivalent sampling was introduced, in which the number of sampling bundles for cell (i, j) was determined by
${N_{ij}} = 4{({\kappa _a})_{ij}}{V_{ij}}{N_0} \;{\rm{or}}\; {N_{ij}} = {\varepsilon _{ij}}{A_{ij}}{N_0}$ | (12) |
where ${N_0}$ is the sampling density for a reference case defined as
${N_0} = \frac{{{N_r}}}{{4{{({\kappa _a})}_r}{V_r}}}$ | (13) |
where, ${N_r}$ is the number of sample bundles of the reference cell ${V_r}$, the subscript r refers to the reference cell.
The bidirectional counting of sampling bundles employed in the BDMC method is equivalent to double the sampling number if the random error does not exist. In fact, any of the Monte Carlo methods always accompanied by random error. For the BDMC method, that leads to $N_{ij}^ * \ne 2{N_{ij}}$ and the energy of a bundle is not strictly equal to other cells. For any cell, it can be easily demonstrated by Eq 11, and the radiative exchange factors calculated by the BDMC method satisfied the conservation relation strictly. But the reciprocity relation between a couple of radiative exchange factors cannot be satisfied strictly. According to Eq 11, the following equation can be obtained
$RD_{ij \to kl}^ * \left( {{N_{ij}} + \sum\limits_{k = 1, l = 1}^{M, N} {{N_{kl \to \, ij}}} } \right) = RD_{kl \to ij}^ * \left( {{N_{kl}} + \sum\limits_{i = 1, j = 1}^{M, N} {{N_{ij \to \, kl}}} } \right)$ | (14) |
According to Eq 8, for gray medium, ${N_{ij \to kl}} = {N_{ij}}R{D_{ij \to kl}}$, ${N_{kl \to ij}} = {N_{kl}}R{D_{kl \to ij}}$. Therefore,
$RD_{ij \to kl}^ * \left( {{N_{ij}} + \sum\limits_{k = 1, l = 1}^{M, N} {{N_{kl}}R{D_{kl \to ij}}} } \right) = RD_{kl \to ij}^ * \left( {{N_{kl}} + \sum\limits_{i = 1, j = 1}^{M, N} {{N_{ij}}R{D_{ij \to kl}}} } \right)$ | (15) |
replace Nij and Nkl according to Eq 12, in addition, considering Eqs 4 and 5, then Eq 15 can be transformed into
$\overline {RD} _{ij \to kl}^ * + RD_{ij \to kl}^ * \sum\limits_{k = 1, l = 1}^{M, N} {{{\overline {RD} }_{kl \to \, ij}}} = \overline {RD} _{kl \to ij}^ * + RD_{kl \to ij}^ * \sum\limits_{i = 1, j = 1}^{M, N} {{{\overline {RD} }_{ij \to \, kl}}} $ | (16) |
where, $\overline {RD} _{ij \to kl}^ * $ and $\overline {RD} _{kl \to ij}^ * $ are the modified forms for $RD_{ij \to kl}^ * $ and $RD_{kl \to ij}^ * $, respectively.
Supposing the reciprocity errors for a couple of radiative exchange factors ${\overline {RD} _{ij \to \, kl}}$ and ${\overline {RD} _{kl \to \, ij}}$ predicted in the TMC method are ${\delta _{ij \to \, kl}}$ and ${\delta _{kl \to \, ij}}$, respectively, one can write
${\overline {RD} _{kl \to \, ij}} = {\overline {RD} _{ij \to \, kl}} + {\delta _{ij \to \, kl}} \; {\rm{and}} \; {\overline {RD} _{ij \to \, kl}} = {\overline {RD} _{kl \to \, ij}} + {\delta _{kl \to ij\, }}$ | (17) |
where, it is obvious that ${\delta _{kl \to ij\, }} = - {\delta _{ij \to \, kl}}$.
For radiative exchange factors predicted in the BDMC method, similar equations can be written as follows
$\overline {RD} _{kl \to \, ij}^ * = \overline {RD} _{ij \to \, kl}^ * + \delta _{ij \to \, kl}^ * \; {\rm{and}} \; \overline {RD} _{ij \to \, kl}^ * = \overline {RD} _{kl \to \, ij}^ * + \delta _{kl \to \, ij}^ * $ | (18) |
where, $\delta _{ij \to \, kl}^ * $ and $\delta _{kl \to \, ij}^ * $ are the reciprocity errors for a couple of radiative exchange factors $\overline {RD} _{ij \to \, kl}^ * $ and $\overline {RD} _{kl \to \, ij}^ * $ in the BDMC method, and $\delta {_{kl \to ij}^ * } = - \delta _{ij \to \, kl}^ * $.Combine Eqs 16–18, then
$\delta _{kl \to \, ij}^ * = RD_{kl \to ij}^ * \sum\limits_{i = 1, j = 1}^{M, N} {\left( {{{\overline {RD} }_{kl \to \, ij}} + {\delta _{kl \to ij\, }}} \right)} - RD_{ij \to kl}^ * \sum\limits_{k = 1, l = 1}^{M, N} {\left( {{{\overline {RD} }_{ij \to \, kl}} + {\delta _{ij \to \, kl}}} \right)} $ | (19) |
Considering the conservation relation $\sum\limits_{i = 1, j = 1}^{M, N} {R{D_{\lambda, kl \to ij}}} = 1$, $\sum\limits_{k = 1, l = 1}^{M, N} {R{D_{\lambda, ij \to kl}}} = 1$, ${\delta _{kl \to ij\, }} = - {\delta _{ij \to \, kl}}$, $\delta {_{kl \to ij}^ * } = - \delta _{ij \to \, kl}^ * $ and Eq 4, then, $\sum\limits_{i = 1, j = 1}^{M, N} {{{\overline {RD} }_{kl \to \, ij}} = } 4{({\kappa _a})_{kl}}{V_{kl}}$, and $\sum\limits_{k = 1, l = 1}^{M, N} {{{\overline {RD} }_{ij \to \, kl}}} = 4{({\kappa _a})_{ij}}{V_{ij}}$, moreover, considering Eq 18, the following error relation can be derived from Eq 19
$\delta _{ij \to \, kl}^ * = \frac{1}{2}\left( {RD_{ij \to kl}^ * \sum\limits_{k = 1, l = 1}^{M, N} {{\delta _{ij \to \, kl}}} + RD_{k{\kern 1pt} l \to ij}^ * \sum\limits_{i = 1, j = 1}^{M, N} {{\delta _{ij \to kl\, }}} } \right)$ | (20) |
therefore,
$\left| {\delta _{ij \to \, kl}^ * } \right| \leqslant \frac{1}{2}\left( {RD_{ij \to kl}^ * \cdot \left| {\sum\limits_{k = 1, l = 1}^{M, N} {{\delta _{ij \to \, kl}}} } \right| + RD_{kl \to ij}^ * \cdot \left| {\sum\limits_{i = 1, j = 1}^{M, N} {{\delta _{ij \to kl\, }}} } \right|} \right)$ | (21) |
Because of the randomicity of the reciprocity error ${\delta _{ij \to \, kl}}$, the value of $\left| {\sum\limits_{k = 1, l = 1}^{M, N} {{\delta _{ij \to \, kl}}} } \right|$ would be very close to 0 if M and N tend to infinity. The formulae $\left| {\sum\limits_{k = 1, l = 1}^{M, N} {{\delta _{ij \to \, kl}}} } \right| = O\left({\left| {{\delta _{ij \to \, kl}}} \right|} \right)$ and $\left| {\sum\limits_{k = 1, l = 1}^{M, N} {{\delta _{ij \to \, kl}}} } \right| = O\left({\left| {{\delta _{ij \to \, kl}}} \right|} \right)$ would be valid for most of the discrete elements as $M \times N$ was large in the present study, where $O\left({\left| {{\delta _{ij \to \, kl}}} \right|} \right)$ is less than $\left| {{\delta _{ij \to \, kl}}} \right|$. Note that the values of $RD_{ij \to kl}^ * $ and $RD_{kl \to ij}^ * $ do not greater than unity, therefore
$\left| {\delta _{ij \to \, kl}^ * } \right| \leqslant O\left( {\left| {{\delta _{ij \to \, kl}}} \right|} \right)$ | (22) |
The reciprocity error for a couple of radiative exchange factors predicted by the BDMC method is always smaller than those predicted by the TMC method.
For some cases, the equivalent sampling may be inconvenient or inaccurate. For example, if the differences of apparent radiation characteristics and/or cell volume are very large among different cells, the equivalent sampling will result in a very small sampling number for a cell, and may be a very large sampling number for another. The former is not expected for the statistical analysis, while the latter increase computational cost.
The weight sampling was introduced into the BDMC method to avoid the above shortages of the equivalent sampling. For weight sampling, the number of sampling bundles for a cell was determined only based on the statistical request, but the weight of a sampling bundle was taken into account in the final statistical calculation. If the number of sampling bundles for cell Vij was Nij, then, its bundle weight Wij was defined as
${W_{ij}} = \frac{{{{({\kappa _a})}_{ij}}{V_{ij}}{N_r}}}{{{{({\kappa _a})}_r}{V_r}{N_{ij}}}} , \;{\rm{or}} \; {W_{ij}} = \frac{{{\varepsilon _{ij}}{A_{ij}}{N_r}}}{{4{{({\kappa _a})}_r}{V_r}{N_{ij}}}}$ | (23) |
where Nr is the number of sample bundles for reference cell Vr. Then, the radiative exchange factor $RD_{ij \to kl}^ * $ in the BDMC method should be calculated by
$RD_{ij \to kl}^ * = \frac{{{N_{ij \to kl}} + {N_{kl \to ij}}{W_{kl}}/{W_{ij}}}}{{{N_{ij}} + \sum\limits_{k = 1, l = 1}^{M, N} {\left( {{N_{kl \to \, ij}}{W_{kl}}/{W_{ij}}} \right)} }}$ | (24) |
Eq 11 should be substituted by Eq 24 if the weight sampling is employed.
Radiative transfer in a rectangular domain with absorbing, emitting and isotropic scattering gray medium was solved separately by the BDMC method and the TMC method. The radiative exchange factors, the radiative equilibrium temperature profiles, and the performance parameter defined by Farmer and Howell [1,9], predicted by the BDMC method were compared with those predicted by the TMC method. In addition, the performance of the equivalent sampling and the weight sampling in the MC simulation were also examined.
($i \in [2, M-1]$, $j \in [2, N-1]$) | ||
mesh number | mesh size | |
uniform mesh | $M \times N = 21 \times 21$ | $\Delta {x_{i, j}} = \Delta {y_{i, j}} = const$ |
non-uniform mesh | $M \times N = 21 \times 21$ | $\Delta {x_{i, j}} = \Delta {x_c} \cdot {K^{ - \left| {i - \frac{{M + 1}}{2}} \right|}}$, $\Delta {y_{i, j}} = \Delta {y_c} \cdot {K^{ - \left| {j - \frac{{N + 1}}{2}} \right|}}$ |
The radiative properties of the medium considered were supposed to be uniform and depicted by constant extinction coefficient ${\kappa _e}$, scattering albedo $\omega $, and refractive index n = 1.0. The length and width of the rectangular domain were ${L_x} = {L_y} = L$. Optical thickness was defined as $\tau = {\kappa _e}L$. The walls were assumed to be diffuse and gray, with constant emissivity of $\varepsilon $. The surfaces were imposed to constant temperature of ${T_{WE}}$, ${T_{WW}}$, ${T_{WN}}$, and ${T_{WS}}$, respectively, see Figure 1. The dimensionless coordinates were defined as $X = x/{L_x}$, and $Y = y/{L_y}$. The simulations were conducted using both uniform and non-uniform grids. The discrete parameters of the two types of mesh systems were listed in Table 1, where, xc = yc were the width and length of the largest volume cell in the non-uniform mesh system, and xc = yc, the parameter K = 1.04, M = N = 21, the mesh sizes in the non-uniform mesh system were expressed as $\Delta {x_{ij}} = \frac{{\Delta {x_c}}}{{{K^{\left| {i - 11} \right|}}}}$, and $\Delta {y_{ij}} = \frac{{\Delta {y_c}}}{{{K^{\left| {j - 11} \right|}}}}$. For medium cells, i.e., $i \in [2, M-1]$, and $j \in [2, N-1]$, the maximum width of the cells was $\Delta {x_{i, j}} = \Delta {x_c}$ for i = 11, and the minimum width of the cells was $\Delta {x_{i, j}} = \frac{{\Delta {x_c}}}{{{K^9}}}$ for i = 2 or 20. Similarly, the maximum length of the cells was $\Delta {y_{i, j}} = \Delta {y_c}$ for j = 11, while the minimum length of the cells was $\Delta {y_{i, j}} = \frac{{\Delta {y_c}}}{{{K^9}}}$ for j = 2 or 20. The maximum/minimum ratio of the mesh size is $\left({\Delta {x_c} \cdot \Delta {y_c}} \right)/\left({\frac{{\Delta {x_c}}}{{{K^9}}} \cdot \frac{{\Delta {y_c}}}{{{K^9}}}} \right)$ = 2.026.
The sampling bundles needed in the MC simulation should be estimated as the results were calculated based on statistical computation. Figures 2 and 3 show the radiative equilibrium temperature profiles at dimensionless locations of Y = 0.3 and Y = 0.5 predicted by the TMC method with different sampling bundles of Nb = 104, 105, and 106, respectively. Figure 2 shows the results for absorbing medium, while the results shown in Figure 3 were for absorbing and isotropic scattering medium with scattering albedo of = 0.5. For both the two cases, the wall temperatures were imposed separately as TWE = TWW = 1000 K and TWN = TWS = 1500 K, the emissivity of the walls were assumed to be constant and equal to 0.5, the optical thickness = keL = 15. For both the absorbing and absorbing-scattering medium, the predicted temperature profiles were found to be close to each other for different sampling bundles of Nb = 104, 105, and 106, moreover, the temperature difference for Nb = 105 and Nb = 106 were less than 0.5%, therefore, the results were considered to reach convergence for sampling bundles Nb = 105 in the present study, and the sampling bundles of Nb = 105 was also chosen as the reference sampling bundles Nr in the BDMC method.
In this section, the radiative exchange factors predicted by the BDMC method were compared with those predicted by the TMC method, the numerical simulations were implemented using both uniform and non-uniform grids, meanwhile, the equivalent sampling and the weight sampling were also employed separately. The predicted radiative exchange factors always satisfy the conservation relation strictly, therefore, only the reciprocity for radiative exchange factors was examined in the present study. The relative error of the reciprocity for a couple of radiative exchange factors was defined as
${\Delta _r} = \frac{{2\left| {{{\overline {RD} }_{ij \to kl}} - {{\overline {RD} }_{kl \to ij}}} \right|}}{{{{\overline {RD} }_{ij \to kl}} + {{\overline {RD} }_{kl \to ij}}}} \times 100\% $ | (25) |
where, ${\overline {RD} _{ij \to kl}}$ refers not only to the modified form of the radiative exchange factors given by Eq 8 for the TMC method, but also to those given by Eq 11 (BDMC method, equivalent sampling) and by Eq 24 (BDMC method, weight sampling).
Table 2 shows the modified radiative exchange factors predicted by the BDMC method and the TMC method employing the equivalent sampling and the uniform grids as well as the corresponding relative error Dr defined by Eq 25. The emissivity of the walls were taken as = 0.5, the optical thicknesses of the rectangle medium along the length and width direction were $\tau = 1.0$, the reference sampling bundles was taken as ${N_r} = {10^5}$. It indicates that the reciprocity error for a couple of radiative exchange factors predicted by the BDMC method was at least an order of magnitude less than those predicted by the TMC method with the same sampling bundles. The maximum relative error for all the elements considered did not exceed 0.58% in the BDMC method, while it reached 30.1% in the TMC method.
($\varepsilon = 0.5$, $\tau = 1.0$, $\omega = 0.0$, ${N_r} = {10^5}$, $\sum {{N_{i, j}}} = {\rm{7}}{\rm{.22}} \times {\rm{1}}{{\rm{0}}^{\rm{7}}}$) | |||||||
Cell 1 | Cell 2 | TMC method | BDMC method | ||||
(i, j) | (k, l ) | ${\overline {RD} _{i, j \to k, {\kern 1pt} l}}$ | ${\overline {RD} _{k, l \to i, j}}$ | ${\Delta _r}\left(\% \right)$ | ${\overline {RD} _{i, j \to k, {\kern 1pt} l}}$ | ${\overline {RD} _{k, l \to i, j}}$ | ${\Delta _r}\left(\% \right)$ |
(11, 11) | (1, 11) | 1.115 × 10-4 | 1.065 × 10-4 | 4.59 | 1.091 × 10-4 | 1.091 × 10-4 | 0.00 |
(11, 11) | (11, 1) | 1.044 × 10-4 | 1.098 × 10-4 | 5.04 | 1.071 × 10-4 | 1.066 × 10-4 | 0.47 |
(11, 11) | (11, 21) | 1.014 × 10-4 | 1.065 × 10-4 | 4.90 | 1.040 × 10-4 | 1.041 × 10-4 | 0.10 |
(11, 11) | (21, 11) | 1.021 × 10-4 | 1.098 × 10-4 | 7.27 | 1.060 × 10-4 | 1.061 × 10-4 | 0.09 |
(11, 11) | (2, 2) | 1.410 × 10-5 | 1.747 × 10-5 | 21.3 | 1.579 × 10-5 | 1.586 × 10-5 | 0.44 |
(11, 11) | (2, 20) | 1.936 × 10-5 | 1.684 × 10-5 | 13.9 | 1.811 × 10-5 | 1.813 × 10-5 | 0.11 |
(11, 11) | (20, 2) | 1.305 × 10-5 | 1.768 × 10-5 | 30.1 | 1.537 × 10-5 | 1.546 × 10-5 | 0.58 |
(11, 11) | (20, 20) | 1.452 × 10-5 | 1.873 × 10-5 | 25.3 | 1.663 × 10-5 | 1.670 × 10-5 | 0.42 |
Table 3 shows the radiative exchange factors predicted by the BDMC method using the equivalent sampling and the weight sampling in a non-uniform mesh system. The computing parameters were $\varepsilon = 0.5$, $\tau = 1.0$, $\omega = 0.0$, and ${N_r} = {10^5}$. The reciprocity error for a couple of radiative exchange factors predicted by the BDMC method with the equivalent sampling was unacceptable for the non-uniform mesh, it indicates that the equivalent sampling was not recommended to solve radiative transfer in the BDMC method with non-uniform grids. If there is obvious difference of heat capacity between cells, the equivalent sampling will lead to big difference of sampling number for the same cells. The calculation accuracy will come under influence if the equivalent sampling is applied. However, the BDMC method with the weight sampling can predict quite satisfactory results for all the cells considered.
($\varepsilon = 0.5$, $\tau = 1.0$, $\omega = 0.0$, ${N_r} = {10^5}$, $\sum {{N_{i, j}}} = {\rm{8}}{\rm{.12}} \times {\rm{1}}{{\rm{0}}^{\rm{7}}}$) | |||||||
Cell 1 | Cell 2 | BDMC (equivalent sampling) | BDMC (weight sampling) | ||||
(i, j) | (k, l ) | ${\overline {RD} _{i, j \to k, {\kern 1pt} l}}$ | ${\overline {RD} _{k, l \to i, j}}$ | ${\Delta _r}\left(\% \right)$ | ${\overline {RD} _{i, j \to k, {\kern 1pt} l}}$ | ${\overline {RD} _{k, l \to i, j}}$ | ${\Delta _r}\left(\% \right)$ |
(11, 11) | (1, 11) | 6.730 × 10-4 | 5.000 × 10-4 | 29.5 | 9.413 × 10-4 | 9.414 × 10-4 | 0.01 |
(11, 11) | (11, 1) | 6.697 × 10-4 | 4.950 × 10-4 | 30.0 | 9.364 × 10-4 | 9.364 × 10-4 | 0.00 |
(11, 11) | (11, 21) | 6.541 × 10-4 | 4.855 × 10-4 | 29.6 | 9.257 × 10-4 | 9.239 × 10-4 | 0.19 |
(11, 11) | (21, 11) | 6.654 × 10-4 | 4.940 × 10-4 | 29.6 | 9.335 × 10-4 | 9.332 × 10-4 | 0.03 |
(11, 11) | (2, 2) | 6.993 × 10-5 | 2.837 × 10-5 | 84.6 | 4.374 × 10-5 | 4.384 × 10-5 | 0.23 |
(11, 11) | (2, 20) | 7.047 × 10-5 | 2.855 × 10-5 | 84.7 | 5.052 × 10-5 | 5.049 × 10-5 | 0.06 |
(11, 11) | (20, 2) | 6.849 × 10-5 | 2.783 × 10-5 | 84.4 | 4.442 × 10-5 | 4.452 × 10-5 | 0.22 |
(11, 11) | (20, 20) | 7.317 × 10-5 | 2.971 × 10-5 | 84.5 | 4.378 × 10-5 | 4.386 × 10-5 | 0.18 |
Table 4 shows the radiative exchange factors predicted by the BDMC method and the TMC method employing the weight sampling and the non-uniform grids. It can be seen that the reciprocity error in the TMC method was larger than those in the BDMC method for all the cells considered. The BDMC method showed greatly superior to the TMC method when the weight sampling and the non-uniform grids was used.
($\varepsilon = 0.5$, $\tau = 1.0$, $\omega = 0.0$, ${N_r} = {10^5}$, $\sum {{N_{i, j}}} = {\rm{8}}{\rm{.12}} \times {\rm{1}}{{\rm{0}}^{\rm{7}}}$) | |||||||
Cell 1 | Cell 2 | TMC method | BDMC method | ||||
(i, j) | (k, l) | ${\overline {RD} _{i, j \to k, {\kern 1pt} l}}$ | ${\overline {RD} _{k, l \to i, j}}$ | (i, j) | (k, l) | ${\overline {RD} _{i, j \to k, {\kern 1pt} l}}$ | ${\overline {RD} _{k, l \to i, j}}$ |
(11, 11) | (1, 11) | 9.492 × 10-4 | 9.344 × 10-4 | (11, 11) | (1, 11) | 9.492 × 10-4 | 9.344 × 10-4 |
(11, 11) | (11, 1) | 9.485 × 10-4 | 9.252 × 10-4 | (11, 11) | (11, 1) | 9.485 × 10-4 | 9.252 × 10-4 |
(11, 11) | (11, 21) | 9.180 × 10-4 | 9.344 × 10-4 | (11, 11) | (11, 21) | 9.180 × 10-4 | 9.344 × 10-4 |
(11, 11) | (21, 11) | 9.428 × 10-4 | 9.252 × 10-4 | (11, 11) | (21, 11) | 9.428 × 10-4 | 9.252 × 10-4 |
(11, 11) | (2, 2) | 4.156 × 10-5 | 4.596 × 10-5 | (11, 11) | (2, 2) | 4.156 × 10-5 | 4.596 × 10-5 |
(11, 11) | (2, 20) | 5.664 × 10-5 | 4.446 × 10-5 | (11, 11) | (2, 20) | 5.664 × 10-5 | 4.446 × 10-5 |
(11, 11) | (20, 2) | 4.223 × 10-5 | 4.667 × 10-5 | (11, 11) | (20, 2) | 4.223 × 10-5 | 4.667 × 10-5 |
(11, 11) | (20, 20) | 3.988 × 10-5 | 4.772 × 10-5 | (11, 11) | (20, 20) | 3.988 × 10-5 | 4.772 × 10-5 |
The advantage of the BDMC method was further verified by comparing the predicted temperature profiles. For gray medium, the radiative equilibrium temperature ${T_{ij}}$ of cell (i, j) satisfies the energy equation
$4{({\kappa _a})_{ij}}{V_{ij}}T_{ij}^4 = \sum\limits_{k = 1, l = 1}^{M, N} {{{\overline {RD} }_{kl \to \, ij}}T_{kl}^4} $ | (26) |
or
${\varepsilon _{ij}}{A_{ij}}T_{ij}^4 = \sum\limits_{k = 1, l = 1}^{M, N} {{{\overline {RD} }_{kl \to \, ij}}T_{kl}^4} $ | (27) |
Utilizing the predictions of radiative exchange factors, the radiative equilibrium temperature profile can be solved iteratively from Eq 26 or Eq 27. In the iteration, the iteration is stopped by setting the residuals. The number of iterations is not certain, which is related to the residuals.
Figure 4 shows the radiative equilibrium temperature profiles for an absorbing, emitting, and non-scattering medium predicted by the TMC method and the BDMC method employing the equivalent sampling and the uniform grids. The wall temperatures were imposed separately as TWE = TWW =1000 K and TWN = TWS =1500 K, the emissivity of the walls were assumed to be constant and equal to 0.5, the optical thickness = keLx = keLy = 15. First, the BDMC predictions were very close to the TMC predictions, it indicates that the development of the present BDMC for predicting temperature profiles is correct, in addition, the BDMC predictions were shown smoother than those predicted by the TMC method, the BDMC method converged faster than the TMC method. Similar superior of the BDMC method to the TMC method can be seen from Figure 5, where the weight sampling and the non-uniform grids were employed and the medium was absorbing and isotropic scattering with scattering albedo of = 0.5.
Figure 6 shows the comparison of temperature profiles predicted by the BDMC method with non-uniform grids employing different sampling models of the equivalent sampling and the weight sampling. The temperature profiles predicted by using the weight sampling was smooth, while the predictions using equivalent sampling model contains some small rectangles, this may due to the fact that the number of the sampling bundles for different cells employing equivalent sampling were significantly different, which lead to large random error and affected the temperature profiles.
Figure 7 shows the radiative equilibrium temperature profiles predicted by the TMC method and the BDMC method with a relative small number of weight sampling bundles using the non-uniform grids. The temperature profiles predicted by the TMC method became unacceptable as the reference sampling bundles decreased to Nr = 103, while the BDMC method can still predict better results, the BDMC method converged faster than the TMC method. Therefore, the BDMC method would be a better choice if the sampling bundles for radiation computation is limited or efficient computation is required.
Farmer and Howell introduced the performance parameter P to evaluate various MC methods or strategies [1,9], the performance parameter was defined as
${\rm{P}} = {\gamma ^2} \cdot t$ | (28) |
where, t is the CPU time spent by the concerned MC simulation and ${\gamma ^2}$ is the variance of the results. A good method or strategy for the MC simulation tends to minimize the performance parameter P. The variance ${\gamma ^2}$ is given by
${\gamma ^2} = \frac{1}{H}\sum\limits_{h = 1}^H {\gamma _h^2} $ | (29) |
where, $\gamma _h^2$ is the variance of the hth exercise, $h \in [1, H]$, and H is the total exercise times, $\gamma _h^2$ can be calculated from
$\gamma _h^2 = \frac{1}{{{M^2} \cdot {N^2}}}{\sum\limits_{i = 1, j = 1}^{M, N} {\sum\limits_{k = 1, l = 1}^{M, N} {\left[ {{{\left( {{{\overline {RD} }_{ij \to \, kl}}} \right)}_h} - {{\left( {{{\overline {RD} }_{ij \to \, kl}}} \right)}_a}} \right]} } ^2}$ | (30) |
where, ${\left({{{\overline {RD} }_{ij \to \, kl}}} \right)_h}$ is the predicted value of ${\overline {RD} _{ij \to \, kl}}$ in the hth exercise, and ${\left({{{\overline {RD} }_{ij \to \, kl}}} \right)_a}$ is the "exact value" of ${\overline {RD} _{ij \to \, kl}}$.
It has been discussed that the acceptable numerical results can be obtained if the sampling bundles Nb in the MC simulation reached 105. In order to compute the performance parameters for the TMC method and the BDMC method, the "exact solution" of the radiative exchange factors should be firstly obtained. The approximate "exact solution" can be obtained by increasing the sampling bundles of the MC simulation due to the high precision of the MC method. To do so, the reference sampling bundles were taken as Nr = 106 and Nr = 108, respectively, the predicted results for the two cases were found to be nearly unchanged, and the maximum relative error for any of the radiative exchange factors correspond to different cells was less than 0.1% as the reference sampling bundles Nr increased from 106 to 108. Therefore, the results employing Nr = 108 were treated as the "exact solution" in the present study.
For each exercise, only the initial random number and sampling order were changed, therefore, the CPU time spent by each exercise should be nearly unchanged. The total exercise times was set as H = 10. The performance parameter of the BDMC method and the TMC method were written as ${{\rm{P}}_{BDM}}$ and ${{\rm{P}}_{TM}}$, respectively, while the performance increment $\Delta {\rm{P}}$ was defined as
$\Delta {\rm{P}} = \frac{{{{\rm{P}}_{TM}} - {{\rm{P}}_{BDM}}}}{{{{\rm{P}}_{TM}}}} \times 100\% $ | (31) |
Table 5 reports the performance parameter for the TMC method and the BDMC method using the equivalent sampling and the uniform grids. It indicates that the performance parameter for the TMC method was always larger than that for the BDMC method, in addition, the performance increment decrease with the increasing sampling bundles. Table 6 shows the performance parameter for the TMC method and the BDMC method employing the weight sampling and the non-uniform grids. Similar conclusions can be drawn as those employing the equivalent sampling and the uniform grids.
($\varepsilon = 0.5$, $\omega = 0.5$, ${N_r} = {10^3}$) | |||
PBDM | PTM | $\Delta {\rm{P}}$ | |
$\tau = 1.0$ $\sum {{N_{i, j}}} = {\rm{7}}{\rm{.22}} \times {\rm{1}}{{\rm{0}}^{\rm{5}}}$ |
8.24 × 10-10 | 1.38 × 10-9 | 40.2 % |
$\tau = 0.01$ $\sum {{N_{i, j}}} = {\rm{3}}{\rm{.64}} \times {\rm{1}}{{\rm{0}}^{\rm{7}}}$ |
6.13 × 10-10 | 7.91 × 10-10 | 22.4% |
($\varepsilon = 0.5$, $\omega = 0.5$, ${N_r} = {10^3}$) | |||
PBDM | PTM | $\Delta {\rm{P}}$ | |
$\tau = 1.0$ $\sum {{N_{i, j}}} = 8.{\rm{12}} \times {\rm{1}}{{\rm{0}}^{\rm{5}}}$ |
9.98 × 10-10 | 1.57 × 10-9 | 36.56 % |
$\tau = 0.01$ $\sum {{N_{i, j}}} = {\rm{4}}{\rm{.54}} \times {\rm{1}}{{\rm{0}}^{\rm{7}}}$ |
1.78 × 10-10 | 2.11 × 10-10 | 15.68% |
Table 7 shows the performance parameter for the BDMC method with equivalent sampling and the weight sampling using the non-uniform grids. The performance parameter for BDMC with the equivalent sampling was found to be larger than that with the weight sampling, this indicate that the weight sampling was more suitable for non-uniform grids in the BDMC simulations.
($\varepsilon = 0.5$, $\omega = 0.5$, ${N_r} = {10^3}$) | |||
weight sampling | equivalent sampling | $\Delta {\rm{P}}$ | |
$\tau = 1.0$ $\sum {{N_{i, j}}} = 8.{\rm{12}} \times {\rm{1}}{{\rm{0}}^{\rm{5}}}$ |
9.98 × 10-10 | 2.55 × 10-9 | 60.84% |
$\tau = 0.01$ $\sum {{N_{i, j}}} = {\rm{4}}{\rm{.54}} \times {\rm{1}}{{\rm{0}}^{\rm{7}}}$ |
1.78 × 10-10 | 6.91 × 10-10 | 74.24% |
The BDMC method along with the equivalent sampling and the weight sampling were developed to solve radiative transfer in absorbing and scattering medium. The BDMC method approximately doubles the information from the bundle tracing to possess a high efficiency by making the best use of the reversibility of the bundle trajectory. The formula for the BDMC method and the corresponding error analysis were derived and presented. Radiative heat transfer in a two-dimensional rectangular domain of absorbing and/or scattering medium were solved by the TMC method and the BDMC method, respectively. The reciprocity of the radiative exchange factors, the radiative equilibrium temperature profiles, and the performance parameter predicted by the two MC methods were examined and compared. The results showed that (1) the BDMC method can greatly improve the reciprocity satisfaction of radiative exchange factors, which is helpful for temperature profile solution, (2) the BDMC method was more accurate or efficient than the TMC method, (3) in the BDMC simulations, the weight sampling was found to be more flexible than the equivalent sampling.
The authors acknowledge the financial support by National Science Foundation of China (No. 51776052) and Aeronautical Science Foundation of China (No.201927077002).
The authors declare no conflict of interest.
Xiaofeng Zhang: Software, Formal analysis Data Curation, Writing - Original Draft, Visualization. Qing Ai: Conceptualization, Methodology, Validation, Resources, Writing - Review & Editing, Supervision, Funding acquisition. Kuilong Song: Validation, Formal analysis, Writing - Review & Editing, Heping Tan: Writing - Review & Editing, Supervision, Project administration.
[1] |
Wang R (2010) Hydrogen sulfide: the third gasotransmitter in biology and medicine. Antioxid Redox Sign 12: 1061–1064. doi: 10.1089/ars.2009.2938
![]() |
[2] |
Ritter JM (2010) Human pharmacology of hydrogen sulfide, putative gaseous mediator. Br J Clin Pharmacol 69: 573–575. doi: 10.1111/j.1365-2125.2010.03690.x
![]() |
[3] |
Li L, Moore PK (2007) An overview of the biological significance of endogenous gases: new roles for old molecules. Biochem Soc Trans 35: 1138–1141. doi: 10.1042/BST0351138
![]() |
[4] |
Szabo C (2016) Gasotransmitters in cancer: from pathophysiology to experimental therapy. Nat Rev Drug Discov 15: 185–203. doi: 10.1038/nrd.2015.1
![]() |
[5] |
van den Born JC, Hammes HP, Greffrath W, et al. (2016) Gasotransmitters in vascular complications of diabetes. Diabetes 65: 331–345. doi: 10.2337/db15-1003
![]() |
[6] | Lowicka E, Beltowski J (2007) Hydrogen sulfide (H2S)-the third gas of interest for pharmacologists. Pharmacol Rep 59: 4–24. |
[7] |
Kolluru GK, Shen X, Bir SC, et al. (2013) Hydrogen sulfide chemical biology: pathophysiological roles and detection. Nitric Oxide 35: 5–20. doi: 10.1016/j.niox.2013.07.002
![]() |
[8] |
Kram L, Grambow E, Mueller-Graf F, et al. (2013) The anti-thrombotic effect of hydrogen sulfide is partly mediated by an upregulation of nitric oxide synthases. Thromb Res 132: e112–117. doi: 10.1016/j.thromres.2013.07.010
![]() |
[9] |
Wang R (2003) The gasotransmitter role of hydrogen sulfide. Antioxid Redox Signal 5: 493–501. doi: 10.1089/152308603768295249
![]() |
[10] |
Pong WW, Stouracova R, Frank N, et al. (2007) Comparative localization of cystathionine beta-synthase and cystathionine gamma-lyase in retina: differences between amphibians and mammals. J Comp Neurol 505: 158–165. doi: 10.1002/cne.21468
![]() |
[11] |
Kimura H, Shibuya N, Kimura Y (2012) Hydrogen sulfide is a signaling molecule and a cytoprotectant. Antioxid Redox Signal 17: 45–57. doi: 10.1089/ars.2011.4345
![]() |
[12] |
Kimura H, Nagai Y, Umemura K, et al. (2005) Physiological roles of hydrogen sulfide: synaptic modulation, neuroprotection, and smooth muscle relaxation. Antioxid Redox Signal 7: 795–803. doi: 10.1089/ars.2005.7.795
![]() |
[13] |
Kimura H (2002) Hydrogen sulfide as a neuromodulator. Mol Neurobiol 26: 13–19. doi: 10.1385/MN:26:1:013
![]() |
[14] |
Kimura H (2011) Hydrogen sulfide: its production, release and functions. Amino Acids 41: 113–121. doi: 10.1007/s00726-010-0510-x
![]() |
[15] |
Persa C, Osmotherly K, Chen KC-W, et al. (2006) The distribution of cystathionine beta-synthase (CBS) in the eye: implication of the presence of a trans-sulfuration pathway for oxidative stress defense. Exp Eye Res 83: 817–823. doi: 10.1016/j.exer.2006.04.001
![]() |
[16] | Shibuya N , Kimura H (2013) Production of hydrogen sulfide from d-cysteine and its therapeutic potential. Front Endocrinol (Lausanne) 4: 87. |
[17] |
Shibuya N, Koike S, Tanaka M, et al. (2013) A novel pathway for the production of hydrogen sulfide from D-cysteine in mammalian cells. Nat Commun 4: 1366. doi: 10.1038/ncomms2371
![]() |
[18] |
Shibuya N, Mikami Y, Kimura Y, et al. (2009) Vascular endothelium expresses 3-mercaptopyruvate sulfurtransferase and produces hydrogen sulfide. J Biochem 146: 623–626. doi: 10.1093/jb/mvp111
![]() |
[19] |
Shibuya N, Tanaka M, Yoshida M, et al. (2009) 3-Mercaptopyruvate sulfurtransferase produces hydrogen sulfide and bound sulfane sulfur in the brain. Antioxid Redox Signal 11: 703–714. doi: 10.1089/ars.2008.2253
![]() |
[20] |
Tanizawa K (2011) Production of H2S by 3-mercaptopyruvate sulphurtransferase. J Biochem 149: 357–359. doi: 10.1093/jb/mvr018
![]() |
[21] |
Ali MY, Ping CY, Mok YY, et al. (2006) Regulation of vascular nitric oxide in vitro and in vivo; a new role for endogenous hydrogen sulphide? Br J Pharmacol 149: 625–634. doi: 10.1038/sj.bjp.0706906
![]() |
[22] |
Li L, Salto-Tellez M, Tan CH, et al. (2009) GYY4137, a novel hydrogen sulfide-releasing molecule, protects against endotoxic shock in the rat. Free Radic Biol Med 47: 103–113. doi: 10.1016/j.freeradbiomed.2009.04.014
![]() |
[23] |
Wang R (2002) Two's company, three's a crowd: can H2S be the third endogenous gaseous transmitter? FASEB J 16: 1792–1798. doi: 10.1096/fj.02-0211hyp
![]() |
[24] |
Whiteman M, Le Trionnaire S, Chopra M, et al. (2011) Emerging role of hydrogen sulfide in health and disease: critical appraisal of biomarkers and pharmacological tools. Clin Sci (Lond) 121: 459–488. doi: 10.1042/CS20110267
![]() |
[25] |
Huang X, Meng XM, Liu DH, et al. (2013) Different regulatory effects of hydrogen sulfide and nitric oxide on gastric motility in mice. Eur J Pharmacol 720: 276–285. doi: 10.1016/j.ejphar.2013.10.017
![]() |
[26] |
Yong QC, Cheong JL, Hua F, et al. (2011) Regulation of heart function by endogenous gaseous mediators-crosstalk between nitric oxide and hydrogen sulfide. Antioxid Redox Signal 14: 2081–2091. doi: 10.1089/ars.2010.3572
![]() |
[27] |
Whiteman M, Li L, Kostetski I, et al. (2006) Evidence for the formation of a novel nitrosothiol from the gaseous mediators nitric oxide and hydrogen sulphide. Biochem Biophys Res Commun 343: 303–310. doi: 10.1016/j.bbrc.2006.02.154
![]() |
[28] |
Zhang QY, Du JB, Zhou WJ, et al. (2004) Impact of hydrogen sulfide on carbon monoxide/heme oxygenase pathway in the pathogenesis of hypoxic pulmonary hypertension. Biochem Biophys Res Commun 317: 30–37. doi: 10.1016/j.bbrc.2004.02.176
![]() |
[29] |
Jin HF, Du JB, Li XH, et al. (2006) Interaction between hydrogen sulfide/cystathionine gamma-lyase and carbon monoxide/heme oxygenase pathways in aortic smooth muscle cells. Acta Pharmacol Sin 27: 1561–1566. doi: 10.1111/j.1745-7254.2006.00425.x
![]() |
[30] | Kulkarni Chitnis M, Belford R, Robinson J, et al. (2014) Interaction between hydrogen sulfide and nitric oxide in isolated bovine retina (1060.2). FASEB J 28, 1 Supplement: 1060–1062. |
[31] |
Yetik-Anacak G, Dereli MV, Sevin G, et al. (2015) Resveratrol stimulates hydrogen sulfide (H2S) formation to relax murine corpus cavernosum. J Sex Med 12: 2004–2012. doi: 10.1111/jsm.12993
![]() |
[32] |
Dufton N, Natividad J, Verdu EF, et al. (2012) Hydrogen sulfide and resolution of acute inflammation: A comparative study utilizing a novel fluorescent probe. Sci Rep 2: 499. doi: 10.1038/srep00499
![]() |
[33] |
Fitzgerald R, DeSantiago B, Lee DY, et al. (2014) H2S relaxes isolated human airway smooth muscle cells via the sarcolemmal K(ATP) channel. Biochem Biophys Res Commun 446: 393–398. doi: 10.1016/j.bbrc.2014.02.129
![]() |
[34] | Zheng Y, Liao F, Du JB, et al. (2012) Modified methylene blue method for measurement of hydrogen sulfide level in plasma. Sheng Li Xue Bao: [Acta Physiologica Sinica] 64: 681–686. |
[35] |
Whiteman M, Moore PK (2009) Hydrogen sulfide and the vasculature: a novel vasculoprotective entity and regulator of nitric oxide bioavailability? J Cell Mol Med 13: 488–507. doi: 10.1111/j.1582-4934.2009.00645.x
![]() |
[36] |
Xia M, Chen L, Muh RW, et al. (2009) Production and actions of hydrogen sulfide, a novel gaseous bioactive substance, in the kidneys. J Pharmacol Exp Ther 329: 1056–1062. doi: 10.1124/jpet.108.149963
![]() |
[37] |
Kulkarni M, Njie-Mbye YF, Okpobiri I, et al. (2011) Endogenous production of hydrogen sulfide in isolated bovine eye. Neurochem Res 36: 1540–1545. doi: 10.1007/s11064-011-0482-6
![]() |
[38] |
Wu D, Hu Q, Zhu Y (2016) Therapeutic application of hydrogen sulfide donors: the potential and challenges. Front Med 10: 18–27. doi: 10.1007/s11684-015-0427-6
![]() |
[39] | Kolesnikov SI, Vlasov BY, Kolesnikova LI (2015) Hydrogen sulfide as a third essential gas molecule in living tissues. Vestn Ross Akad Med Nauk 2: 237–241. |
[40] |
van Goor H, van den Born JC, Hillebrands JL, et al. (2016) Hydrogen sulfide in hypertension. Curr Opin Nephrol Hypertens 25: 107–113. doi: 10.1097/MNH.0000000000000206
![]() |
[41] |
Ahmad A, Sattar MA, Rathore HA, et al. (2015) A critical review of pharmacological significance of hydrogen sulfide in hypertension. Indian J Pharmacol 47: 243–247. doi: 10.4103/0253-7613.157106
![]() |
[42] |
Kida K, Ichinose F (2015) Hydrogen sulfide and neuroinflammation. Handb Exp Pharmacol 230: 181–189. doi: 10.1007/978-3-319-18144-8_9
![]() |
[43] | Cui Y, Duan X, Li H, et al. (2015) Hydrogen sulfide ameliorates early brain injury following subarachnoid hemorrhage in rats. Mol Neurobiol 53: 3646–3657. |
[44] |
Wang YF, Mainali P, Tang CS, et al. (2008) Effects of nitric oxide and hydrogen sulfide on the relaxation of pulmonary arteries in rats. Chinese Med J 121: 420–423. doi: 10.1097/00029330-200803010-00010
![]() |
[45] |
Altaany Z, Yang G, Wang R (2013) Crosstalk between hydrogen sulfide and nitric oxide in endothelial cells. J Cell Mol Med 17: 879–888. doi: 10.1111/jcmm.12077
![]() |
[46] |
Pong WW, Eldred WD (2009) Interactions of the gaseous neuromodulators nitric oxide, carbon monoxide, and hydrogen sulfide in the salamander retina. J Neurosci Res 87: 2356–2364. doi: 10.1002/jnr.22042
![]() |
[47] | Salomone S, Foresti R, Villari A, et al. (2014) Regulation of vascular tone in rabbit ophthalmic artery: cross talk of endogenous and exogenous gas mediators. Biochem Pharmacol 92: 4661–4668. |
[48] |
Zhao W, Wang R (2002) H(2)S-induced vasorelaxation and underlying cellular and molecular mechanisms. Am J Physiol-Heart C 283: H474–480. doi: 10.1152/ajpheart.00013.2002
![]() |
[49] |
Zhao W, Zhang J, Lu Y, et al. (2001) The vasorelaxant effect of H(2)S as a novel endogenous gaseous K(ATP) channel opener. EMBO J 20: 6008–6016. doi: 10.1093/emboj/20.21.6008
![]() |
[50] |
Zhao W, Ndisang JF, Wang R (2003) Modulation of endogenous production of H2S in rat tissues. Can J Physiol Pharmacol 81: 848–853. doi: 10.1139/y03-077
![]() |
[51] | Guo W, Kan JT, Cheng ZY, et al. (2012) Hydrogen sulfide as an endogenous modulator in mitochondria and mitochondria dysfunction. Oxid Med Cell Longev 2012. |
[52] |
Yanfei W, Lin S, Junbao D, et al. (2006) Impact of L-arginine on hydrogen sulfide/cystathionine-gamma-lyase pathway in rats with high blood flow-induced pulmonary hypertension. Biochem Biophys Res Commun 345: 851–857. doi: 10.1016/j.bbrc.2006.04.162
![]() |
[53] | Brancaleone V, Roviezzo F, Vellecco V, et al. (2008) Biosynthesis of H2S is impaired in non-obese diabetic (NOD) mice. Br J Pharmacol 155: 673–680. |
[54] |
Dyson RM, Palliser HK, Latter JL, et al. (2015) Interactions of the gasotransmitters contribute to microvascular tone (dys)regulation in the preterm neonate. PLoS One 10: e0121621. doi: 10.1371/journal.pone.0121621
![]() |
[55] |
Holwerda KM, Faas MM, van Goor H, et al. (2013) Gasotransmitters: a solution for the therapeutic dilemma in preeclampsia? Hypertension 62: 653–659. doi: 10.1161/HYPERTENSIONAHA.113.01625
![]() |
[56] |
Hosoki R, Matsuki N, Kimura H (1997) The possible role of hydrogen sulfide as an endogenous smooth muscle relaxant in synergy with nitric oxide. Biochem Biophys Res Commun 237: 527–531. doi: 10.1006/bbrc.1997.6878
![]() |
[57] |
Privitera MG, Potenza M, Bucolo C, et al. (2007) Hemin, an inducer of heme oxygenase-1, lowers intraocular pressure in rabbits. J Ocul Pharmacol Ther 23: 232–239. doi: 10.1089/jop.2006.101
![]() |
[58] |
Stagni E, Bucolo C, Motterlini R, et al. (2010) Morphine-induced ocular hypotension is modulated by nitric oxide and carbon monoxide: role of μ3 receptors. J Ocul Pharmacol Ther 26: 31–36. doi: 10.1089/jop.2009.0081
![]() |
[59] |
Peng YJ, Nanduri J, Raghuraman G, et al. (2010) H2S mediates O2 sensing in the carotid body. Proc Natl Acad Sci USA 107: 10719–10724. doi: 10.1073/pnas.1005866107
![]() |
[60] |
Shintani T, Iwabuchi T, Soga T, et al. (2009) Cystathionine beta-synthase as a carbon monoxide-sensitive regulator of bile excretion. Hepatology 49: 141–150. doi: 10.1002/hep.22604
![]() |
[61] |
Morikawa T, Kajimura M, Nakamura T, et al. (2012) Hypoxic regulation of the cerebral microcirculation is mediated by a carbon monoxide-sensitive hydrogen sulfide pathway. Proc Natl Acad Sci USA 109: 1293–1298. doi: 10.1073/pnas.1119658109
![]() |
[62] |
Monjok EM, Kulkarni KH, Kouamou G, et al (2008) Inhibitory action of hydrogen sulfide on muscarinic receptor-induced contraction of isolated porcine irides. Exp Eye Res 87: 612–616. doi: 10.1016/j.exer.2008.09.011
![]() |
[63] |
Kulkarni-Chitnis M, Njie-Mbye YF, Mitchell L, et al. (2015) Inhibitory action of novel hydrogen sulfide donors on bovine isolated posterior ciliary arteries. Exp Eye Res 134: 73–79. doi: 10.1016/j.exer.2015.04.001
![]() |
[64] |
Chitnis MK, Njie-Mbye YF, Opere CA, et al. (2013) Pharmacological actions of the slow release hydrogen sulfide donor GYY4137 on phenylephrine-induced tone in isolated bovine ciliary artery. Exp Eye Res 116: 350–354. doi: 10.1016/j.exer.2013.10.004
![]() |
1. | V. A. Kuznetsov, Modern Methods for Numerical Simulation of Radiation Heat Transfer in Selective Gases (Review), 2022, 69, 0040-6015, 702, 10.1134/S0040601522080043 |
($i \in [2, M-1]$, $j \in [2, N-1]$) | ||
mesh number | mesh size | |
uniform mesh | $M \times N = 21 \times 21$ | $\Delta {x_{i, j}} = \Delta {y_{i, j}} = const$ |
non-uniform mesh | $M \times N = 21 \times 21$ | $\Delta {x_{i, j}} = \Delta {x_c} \cdot {K^{ - \left| {i - \frac{{M + 1}}{2}} \right|}}$, $\Delta {y_{i, j}} = \Delta {y_c} \cdot {K^{ - \left| {j - \frac{{N + 1}}{2}} \right|}}$ |
($\varepsilon = 0.5$, $\tau = 1.0$, $\omega = 0.0$, ${N_r} = {10^5}$, $\sum {{N_{i, j}}} = {\rm{7}}{\rm{.22}} \times {\rm{1}}{{\rm{0}}^{\rm{7}}}$) | |||||||
Cell 1 | Cell 2 | TMC method | BDMC method | ||||
(i, j) | (k, l ) | ${\overline {RD} _{i, j \to k, {\kern 1pt} l}}$ | ${\overline {RD} _{k, l \to i, j}}$ | ${\Delta _r}\left(\% \right)$ | ${\overline {RD} _{i, j \to k, {\kern 1pt} l}}$ | ${\overline {RD} _{k, l \to i, j}}$ | ${\Delta _r}\left(\% \right)$ |
(11, 11) | (1, 11) | 1.115 × 10-4 | 1.065 × 10-4 | 4.59 | 1.091 × 10-4 | 1.091 × 10-4 | 0.00 |
(11, 11) | (11, 1) | 1.044 × 10-4 | 1.098 × 10-4 | 5.04 | 1.071 × 10-4 | 1.066 × 10-4 | 0.47 |
(11, 11) | (11, 21) | 1.014 × 10-4 | 1.065 × 10-4 | 4.90 | 1.040 × 10-4 | 1.041 × 10-4 | 0.10 |
(11, 11) | (21, 11) | 1.021 × 10-4 | 1.098 × 10-4 | 7.27 | 1.060 × 10-4 | 1.061 × 10-4 | 0.09 |
(11, 11) | (2, 2) | 1.410 × 10-5 | 1.747 × 10-5 | 21.3 | 1.579 × 10-5 | 1.586 × 10-5 | 0.44 |
(11, 11) | (2, 20) | 1.936 × 10-5 | 1.684 × 10-5 | 13.9 | 1.811 × 10-5 | 1.813 × 10-5 | 0.11 |
(11, 11) | (20, 2) | 1.305 × 10-5 | 1.768 × 10-5 | 30.1 | 1.537 × 10-5 | 1.546 × 10-5 | 0.58 |
(11, 11) | (20, 20) | 1.452 × 10-5 | 1.873 × 10-5 | 25.3 | 1.663 × 10-5 | 1.670 × 10-5 | 0.42 |
($\varepsilon = 0.5$, $\tau = 1.0$, $\omega = 0.0$, ${N_r} = {10^5}$, $\sum {{N_{i, j}}} = {\rm{8}}{\rm{.12}} \times {\rm{1}}{{\rm{0}}^{\rm{7}}}$) | |||||||
Cell 1 | Cell 2 | BDMC (equivalent sampling) | BDMC (weight sampling) | ||||
(i, j) | (k, l ) | ${\overline {RD} _{i, j \to k, {\kern 1pt} l}}$ | ${\overline {RD} _{k, l \to i, j}}$ | ${\Delta _r}\left(\% \right)$ | ${\overline {RD} _{i, j \to k, {\kern 1pt} l}}$ | ${\overline {RD} _{k, l \to i, j}}$ | ${\Delta _r}\left(\% \right)$ |
(11, 11) | (1, 11) | 6.730 × 10-4 | 5.000 × 10-4 | 29.5 | 9.413 × 10-4 | 9.414 × 10-4 | 0.01 |
(11, 11) | (11, 1) | 6.697 × 10-4 | 4.950 × 10-4 | 30.0 | 9.364 × 10-4 | 9.364 × 10-4 | 0.00 |
(11, 11) | (11, 21) | 6.541 × 10-4 | 4.855 × 10-4 | 29.6 | 9.257 × 10-4 | 9.239 × 10-4 | 0.19 |
(11, 11) | (21, 11) | 6.654 × 10-4 | 4.940 × 10-4 | 29.6 | 9.335 × 10-4 | 9.332 × 10-4 | 0.03 |
(11, 11) | (2, 2) | 6.993 × 10-5 | 2.837 × 10-5 | 84.6 | 4.374 × 10-5 | 4.384 × 10-5 | 0.23 |
(11, 11) | (2, 20) | 7.047 × 10-5 | 2.855 × 10-5 | 84.7 | 5.052 × 10-5 | 5.049 × 10-5 | 0.06 |
(11, 11) | (20, 2) | 6.849 × 10-5 | 2.783 × 10-5 | 84.4 | 4.442 × 10-5 | 4.452 × 10-5 | 0.22 |
(11, 11) | (20, 20) | 7.317 × 10-5 | 2.971 × 10-5 | 84.5 | 4.378 × 10-5 | 4.386 × 10-5 | 0.18 |
($\varepsilon = 0.5$, $\tau = 1.0$, $\omega = 0.0$, ${N_r} = {10^5}$, $\sum {{N_{i, j}}} = {\rm{8}}{\rm{.12}} \times {\rm{1}}{{\rm{0}}^{\rm{7}}}$) | |||||||
Cell 1 | Cell 2 | TMC method | BDMC method | ||||
(i, j) | (k, l) | ${\overline {RD} _{i, j \to k, {\kern 1pt} l}}$ | ${\overline {RD} _{k, l \to i, j}}$ | (i, j) | (k, l) | ${\overline {RD} _{i, j \to k, {\kern 1pt} l}}$ | ${\overline {RD} _{k, l \to i, j}}$ |
(11, 11) | (1, 11) | 9.492 × 10-4 | 9.344 × 10-4 | (11, 11) | (1, 11) | 9.492 × 10-4 | 9.344 × 10-4 |
(11, 11) | (11, 1) | 9.485 × 10-4 | 9.252 × 10-4 | (11, 11) | (11, 1) | 9.485 × 10-4 | 9.252 × 10-4 |
(11, 11) | (11, 21) | 9.180 × 10-4 | 9.344 × 10-4 | (11, 11) | (11, 21) | 9.180 × 10-4 | 9.344 × 10-4 |
(11, 11) | (21, 11) | 9.428 × 10-4 | 9.252 × 10-4 | (11, 11) | (21, 11) | 9.428 × 10-4 | 9.252 × 10-4 |
(11, 11) | (2, 2) | 4.156 × 10-5 | 4.596 × 10-5 | (11, 11) | (2, 2) | 4.156 × 10-5 | 4.596 × 10-5 |
(11, 11) | (2, 20) | 5.664 × 10-5 | 4.446 × 10-5 | (11, 11) | (2, 20) | 5.664 × 10-5 | 4.446 × 10-5 |
(11, 11) | (20, 2) | 4.223 × 10-5 | 4.667 × 10-5 | (11, 11) | (20, 2) | 4.223 × 10-5 | 4.667 × 10-5 |
(11, 11) | (20, 20) | 3.988 × 10-5 | 4.772 × 10-5 | (11, 11) | (20, 20) | 3.988 × 10-5 | 4.772 × 10-5 |
($\varepsilon = 0.5$, $\omega = 0.5$, ${N_r} = {10^3}$) | |||
PBDM | PTM | $\Delta {\rm{P}}$ | |
$\tau = 1.0$ $\sum {{N_{i, j}}} = {\rm{7}}{\rm{.22}} \times {\rm{1}}{{\rm{0}}^{\rm{5}}}$ |
8.24 × 10-10 | 1.38 × 10-9 | 40.2 % |
$\tau = 0.01$ $\sum {{N_{i, j}}} = {\rm{3}}{\rm{.64}} \times {\rm{1}}{{\rm{0}}^{\rm{7}}}$ |
6.13 × 10-10 | 7.91 × 10-10 | 22.4% |
($\varepsilon = 0.5$, $\omega = 0.5$, ${N_r} = {10^3}$) | |||
PBDM | PTM | $\Delta {\rm{P}}$ | |
$\tau = 1.0$ $\sum {{N_{i, j}}} = 8.{\rm{12}} \times {\rm{1}}{{\rm{0}}^{\rm{5}}}$ |
9.98 × 10-10 | 1.57 × 10-9 | 36.56 % |
$\tau = 0.01$ $\sum {{N_{i, j}}} = {\rm{4}}{\rm{.54}} \times {\rm{1}}{{\rm{0}}^{\rm{7}}}$ |
1.78 × 10-10 | 2.11 × 10-10 | 15.68% |
($\varepsilon = 0.5$, $\omega = 0.5$, ${N_r} = {10^3}$) | |||
weight sampling | equivalent sampling | $\Delta {\rm{P}}$ | |
$\tau = 1.0$ $\sum {{N_{i, j}}} = 8.{\rm{12}} \times {\rm{1}}{{\rm{0}}^{\rm{5}}}$ |
9.98 × 10-10 | 2.55 × 10-9 | 60.84% |
$\tau = 0.01$ $\sum {{N_{i, j}}} = {\rm{4}}{\rm{.54}} \times {\rm{1}}{{\rm{0}}^{\rm{7}}}$ |
1.78 × 10-10 | 6.91 × 10-10 | 74.24% |
($i \in [2, M-1]$, $j \in [2, N-1]$) | ||
mesh number | mesh size | |
uniform mesh | $M \times N = 21 \times 21$ | $\Delta {x_{i, j}} = \Delta {y_{i, j}} = const$ |
non-uniform mesh | $M \times N = 21 \times 21$ | $\Delta {x_{i, j}} = \Delta {x_c} \cdot {K^{ - \left| {i - \frac{{M + 1}}{2}} \right|}}$, $\Delta {y_{i, j}} = \Delta {y_c} \cdot {K^{ - \left| {j - \frac{{N + 1}}{2}} \right|}}$ |
($\varepsilon = 0.5$, $\tau = 1.0$, $\omega = 0.0$, ${N_r} = {10^5}$, $\sum {{N_{i, j}}} = {\rm{7}}{\rm{.22}} \times {\rm{1}}{{\rm{0}}^{\rm{7}}}$) | |||||||
Cell 1 | Cell 2 | TMC method | BDMC method | ||||
(i, j) | (k, l ) | ${\overline {RD} _{i, j \to k, {\kern 1pt} l}}$ | ${\overline {RD} _{k, l \to i, j}}$ | ${\Delta _r}\left(\% \right)$ | ${\overline {RD} _{i, j \to k, {\kern 1pt} l}}$ | ${\overline {RD} _{k, l \to i, j}}$ | ${\Delta _r}\left(\% \right)$ |
(11, 11) | (1, 11) | 1.115 × 10-4 | 1.065 × 10-4 | 4.59 | 1.091 × 10-4 | 1.091 × 10-4 | 0.00 |
(11, 11) | (11, 1) | 1.044 × 10-4 | 1.098 × 10-4 | 5.04 | 1.071 × 10-4 | 1.066 × 10-4 | 0.47 |
(11, 11) | (11, 21) | 1.014 × 10-4 | 1.065 × 10-4 | 4.90 | 1.040 × 10-4 | 1.041 × 10-4 | 0.10 |
(11, 11) | (21, 11) | 1.021 × 10-4 | 1.098 × 10-4 | 7.27 | 1.060 × 10-4 | 1.061 × 10-4 | 0.09 |
(11, 11) | (2, 2) | 1.410 × 10-5 | 1.747 × 10-5 | 21.3 | 1.579 × 10-5 | 1.586 × 10-5 | 0.44 |
(11, 11) | (2, 20) | 1.936 × 10-5 | 1.684 × 10-5 | 13.9 | 1.811 × 10-5 | 1.813 × 10-5 | 0.11 |
(11, 11) | (20, 2) | 1.305 × 10-5 | 1.768 × 10-5 | 30.1 | 1.537 × 10-5 | 1.546 × 10-5 | 0.58 |
(11, 11) | (20, 20) | 1.452 × 10-5 | 1.873 × 10-5 | 25.3 | 1.663 × 10-5 | 1.670 × 10-5 | 0.42 |
($\varepsilon = 0.5$, $\tau = 1.0$, $\omega = 0.0$, ${N_r} = {10^5}$, $\sum {{N_{i, j}}} = {\rm{8}}{\rm{.12}} \times {\rm{1}}{{\rm{0}}^{\rm{7}}}$) | |||||||
Cell 1 | Cell 2 | BDMC (equivalent sampling) | BDMC (weight sampling) | ||||
(i, j) | (k, l ) | ${\overline {RD} _{i, j \to k, {\kern 1pt} l}}$ | ${\overline {RD} _{k, l \to i, j}}$ | ${\Delta _r}\left(\% \right)$ | ${\overline {RD} _{i, j \to k, {\kern 1pt} l}}$ | ${\overline {RD} _{k, l \to i, j}}$ | ${\Delta _r}\left(\% \right)$ |
(11, 11) | (1, 11) | 6.730 × 10-4 | 5.000 × 10-4 | 29.5 | 9.413 × 10-4 | 9.414 × 10-4 | 0.01 |
(11, 11) | (11, 1) | 6.697 × 10-4 | 4.950 × 10-4 | 30.0 | 9.364 × 10-4 | 9.364 × 10-4 | 0.00 |
(11, 11) | (11, 21) | 6.541 × 10-4 | 4.855 × 10-4 | 29.6 | 9.257 × 10-4 | 9.239 × 10-4 | 0.19 |
(11, 11) | (21, 11) | 6.654 × 10-4 | 4.940 × 10-4 | 29.6 | 9.335 × 10-4 | 9.332 × 10-4 | 0.03 |
(11, 11) | (2, 2) | 6.993 × 10-5 | 2.837 × 10-5 | 84.6 | 4.374 × 10-5 | 4.384 × 10-5 | 0.23 |
(11, 11) | (2, 20) | 7.047 × 10-5 | 2.855 × 10-5 | 84.7 | 5.052 × 10-5 | 5.049 × 10-5 | 0.06 |
(11, 11) | (20, 2) | 6.849 × 10-5 | 2.783 × 10-5 | 84.4 | 4.442 × 10-5 | 4.452 × 10-5 | 0.22 |
(11, 11) | (20, 20) | 7.317 × 10-5 | 2.971 × 10-5 | 84.5 | 4.378 × 10-5 | 4.386 × 10-5 | 0.18 |
($\varepsilon = 0.5$, $\tau = 1.0$, $\omega = 0.0$, ${N_r} = {10^5}$, $\sum {{N_{i, j}}} = {\rm{8}}{\rm{.12}} \times {\rm{1}}{{\rm{0}}^{\rm{7}}}$) | |||||||
Cell 1 | Cell 2 | TMC method | BDMC method | ||||
(i, j) | (k, l) | ${\overline {RD} _{i, j \to k, {\kern 1pt} l}}$ | ${\overline {RD} _{k, l \to i, j}}$ | (i, j) | (k, l) | ${\overline {RD} _{i, j \to k, {\kern 1pt} l}}$ | ${\overline {RD} _{k, l \to i, j}}$ |
(11, 11) | (1, 11) | 9.492 × 10-4 | 9.344 × 10-4 | (11, 11) | (1, 11) | 9.492 × 10-4 | 9.344 × 10-4 |
(11, 11) | (11, 1) | 9.485 × 10-4 | 9.252 × 10-4 | (11, 11) | (11, 1) | 9.485 × 10-4 | 9.252 × 10-4 |
(11, 11) | (11, 21) | 9.180 × 10-4 | 9.344 × 10-4 | (11, 11) | (11, 21) | 9.180 × 10-4 | 9.344 × 10-4 |
(11, 11) | (21, 11) | 9.428 × 10-4 | 9.252 × 10-4 | (11, 11) | (21, 11) | 9.428 × 10-4 | 9.252 × 10-4 |
(11, 11) | (2, 2) | 4.156 × 10-5 | 4.596 × 10-5 | (11, 11) | (2, 2) | 4.156 × 10-5 | 4.596 × 10-5 |
(11, 11) | (2, 20) | 5.664 × 10-5 | 4.446 × 10-5 | (11, 11) | (2, 20) | 5.664 × 10-5 | 4.446 × 10-5 |
(11, 11) | (20, 2) | 4.223 × 10-5 | 4.667 × 10-5 | (11, 11) | (20, 2) | 4.223 × 10-5 | 4.667 × 10-5 |
(11, 11) | (20, 20) | 3.988 × 10-5 | 4.772 × 10-5 | (11, 11) | (20, 20) | 3.988 × 10-5 | 4.772 × 10-5 |
($\varepsilon = 0.5$, $\omega = 0.5$, ${N_r} = {10^3}$) | |||
PBDM | PTM | $\Delta {\rm{P}}$ | |
$\tau = 1.0$ $\sum {{N_{i, j}}} = {\rm{7}}{\rm{.22}} \times {\rm{1}}{{\rm{0}}^{\rm{5}}}$ |
8.24 × 10-10 | 1.38 × 10-9 | 40.2 % |
$\tau = 0.01$ $\sum {{N_{i, j}}} = {\rm{3}}{\rm{.64}} \times {\rm{1}}{{\rm{0}}^{\rm{7}}}$ |
6.13 × 10-10 | 7.91 × 10-10 | 22.4% |
($\varepsilon = 0.5$, $\omega = 0.5$, ${N_r} = {10^3}$) | |||
PBDM | PTM | $\Delta {\rm{P}}$ | |
$\tau = 1.0$ $\sum {{N_{i, j}}} = 8.{\rm{12}} \times {\rm{1}}{{\rm{0}}^{\rm{5}}}$ |
9.98 × 10-10 | 1.57 × 10-9 | 36.56 % |
$\tau = 0.01$ $\sum {{N_{i, j}}} = {\rm{4}}{\rm{.54}} \times {\rm{1}}{{\rm{0}}^{\rm{7}}}$ |
1.78 × 10-10 | 2.11 × 10-10 | 15.68% |
($\varepsilon = 0.5$, $\omega = 0.5$, ${N_r} = {10^3}$) | |||
weight sampling | equivalent sampling | $\Delta {\rm{P}}$ | |
$\tau = 1.0$ $\sum {{N_{i, j}}} = 8.{\rm{12}} \times {\rm{1}}{{\rm{0}}^{\rm{5}}}$ |
9.98 × 10-10 | 2.55 × 10-9 | 60.84% |
$\tau = 0.01$ $\sum {{N_{i, j}}} = {\rm{4}}{\rm{.54}} \times {\rm{1}}{{\rm{0}}^{\rm{7}}}$ |
1.78 × 10-10 | 6.91 × 10-10 | 74.24% |