Application of the lattice Boltzmann method to transient conduction and radiation heat transfer in cylindrical media

In this paper, the lattice Boltzmann method (LBM) is applied to solve the energy equation of a transient conduction-radiation heat transfer problem in a two-dimensional cylindrical enclosure filled with an emitting, absorbing and scattering media. The control volume finite element method (CVFEM) is used to obtain the radiative information. To demonstrate the workability of the LBM in conjunction with the CVFEM to conduction-radiation problems in cylindrical media, the energy equation of the same problem is also solved using the finite difference method (FDM). The effects of different parameters, such as the grid size, the scattering albedo, the extinction coefficient and the conduction-radiation parameter on temperature distribution within the medium are studied. Results of the present work are compared with those available in the literature. LBM-CVFEM results are also compared with those given by the FDM-CVFEM. In all cases, good agreement has been obtained.     

An inverse analysis to determine conductive and radiative properties of a fibrous medium

In present paper, a modified factor of extinction coefficient and an equivalent albedo of scattering were defined taking into account anisotropic scattering in fibrous insulation. An inverse conduction-radiation analysis in an absorbing, emitting and scattering medium was conducted for the simultaneous estimation of the conductive and radiative properties using the experimentally measured temperature responses for external temperatures up to 980 K. The estimated properties were validated by comparing the predicted and measured results under transient and steady-state condition. It was found that the calculated results corresponded well with the experimental data within an average of 3.1% under transient condition and 9.8% under steady-state condition. This confirms the good behavior of the model and the validity of results.     

Identifying time-dependent damping and stiffness functions by a simple and yet accurate method

The inverse vibration problem is a mathematical process to determine unknown mechanical parameters from measured vibration data. In this study the data of displacement are chosen in order to identify a time-dependent function of damping or stiffness. However, when both functions are to be identified we require both the data of displacement and velocity. This is the first time that a closed-form estimation method for the inverse vibration problems of estimating time-dependent parameters has been constructed. We are able to transform the inverse vibration problem into an identification problem governed by a parabolic-type partial differential equation (PDE). Then, a one-step group-preserving scheme (GPS) for the semi-discretization of PDE is established, which can be used to derive a closed-form solution for estimating parameters. The new Lie-group estimation method has three further advantages: it does not require any prior information on the functional forms of unknown functions; no initial guesses are required; and no iterations are required. Numerical examples were examined to show that the present approach is highly accurate and efficient even for identifying discontinuous and oscillatory parameters. Against the noise is good when only one function is estimated; however, the present approach is slightly weak against the noise when both functions are identified.     

An inverse radiative transfer problem of simultaneously estimating profiles of temperature and radiative parameters from boundary intensity and temperature measurements

A parallel-plane space filled with absorbing, emitting, isotropically scattering, gray medium is studied in this paper. The boundary intensity and boundary temperature profiles are calculated for the inverse analysis. For the simultaneous estimation of temperature, absorption and scattering coefficient profiles in the medium, the sum of residuals of boundary intensity and temperature after being weighted by a balance factor is minimized through using a Newton-type iteration algorithm and the least-squares method. To avoid over-updating for the parameters, the relative updating magnitude during the iteration process is constrained not to be >0.5. It is shown that the boundary intensity measurement alone is not enough to estimate simultaneously the temperature (source) and the radiative properties (both absorption and scattering coefficients) when the measurement data contain sensitive random errors. The boundary temperature measurement can serve as a necessary supplementation to the boundary intensity to make this kind of inverse radiative transfer problem resolvable. It was shown that a compensation relationship between absorption and scattering coefficients makes it difficult to fix them accurately. Parabolic profiles for the three parameters are used to validate the estimation method. When the optical thickness approaches 4.0, the results for the radiative properties are not acceptable, although the result for temperature profile is reasonable. This means the method needs further improvements.     

A sensitivity function-based conjugate gradient method for optical tomography with the frequency-domain equation of radiative transfer

The Sensitivity Function-based Conjugate Gradient Method (SFCGM) is described. This method is used to solve the inverse problems of function estimation, such as the local maps of absorption and scattering coefficients, as applied to optical tomography for biomedical imaging. A highly scattering, absorbing, non-reflecting, non-emitting medium is considered here and simultaneous reconstructions of absorption and scattering coefficients inside the test medium are achieved with the proposed optimization technique, by using the exit intensity measured at boundary surfaces. The forward problem is solved with a discrete-ordinates finite-difference method on the framework of the frequency-domain full equation of radiative transfer. The modulation frequency is set to 600 MHz and the frequency data, obtained with the source modulation, is used as the input data. The inversion results demonstrate that the SFCGM can retrieve simultaneously the spatial distributions of optical properties inside the medium within a reasonable accuracy, by significantly reducing a cross-talk between inter-parameters. It is also observed that the closer-to-detector objects are better retrieved.     

Accurate reconstruction of the optical parameter distribution in participating medium based on the frequency-domain radiative transfer equation

Reconstructing the distribution of optical parameters in the participating medium based on the frequency-domain radiative transfer equation(FD-RTE) to probe the internal structure of the medium is investigated in the present work.The forward model of FD-RTE is solved via the finite volume method(FVM). The regularization term formatted by the generalized Gaussian Markov random field model is used in the objective function to overcome the ill-posed nature of the inverse problem. The multi-start conjugate gradient(MCG) method is employed to search the minimum of the objective function and increase the efficiency of convergence. A modified adjoint differentiation technique using the collimated radiative intensity is developed to calculate the gradient of the objective function with respect to the optical parameters. All simulation results show that the proposed reconstruction algorithm based on FD-RTE can obtain the accurate distributions of absorption and scattering coefficients. The reconstructed images of the scattering coefficient have less errors than those of the absorption coefficient, which indicates the former are more suitable to probing the inner structure.     

Synthesis of long-period fiber gratings with a Lagrange multiplier optimization method

Based on the Lagrange multiplier optimization (LMO) method, a new synthesis approach for designing complex long-period fiber grating (LPG) filters is developed and demonstrated. The proposed synthesis method is a simple and direct approach. It was used to efficiently search for optimal solutions and constrain various parameters of the designed LPG filters according to practical requirements. The inverse scattering algorithm is a good tool for designing FBG filters; as well as for designing transmission-type fiber grating filters like LPGs. Compared to the results of the discrete layer-peeling (DLP) inverse scattering algorithm for LPGs, the synthesized LPGs, using the LMO approach, are more flexible and workable for the constraint conditions can be easily set in the user-defined cost functional, such as the limitation on the maximum value of the refractive index modulation and the parameters of the initial guess in LMO algorithm. Linear transmission LPG filters and EDFA gain flattening filters are synthesized and analyzed systematically by using the proposed method with different parameters. Moreover, as a variation-based method, we find that the convergence and the synthesized results of the LMO method are strongly dependent on the initial guess parameters.     

抛物型方程的演化参数识别方法

给出了一种利用演化计算方法求解微分方程中的参数识别类型反问题的方法。该方法把参数识别问题转化为泛函的优化问题用演化算法来求解,指定待定参数的函数类形式,用遗传算法(Genetic Algorithms)来演化待求参数的最优估计值,并将该方法运用于线性扩散方程和拟线性对流扩散反方程反问题的数值模拟中。     

Efficient inverse radiation analysis in a cylindrical geometry using a combined method of hybrid genetic algorithm and finite-difference Newton method

An inverse radiation problem was considered to estimate boundary conditions such as temperature distribution and emissivity in axisymmetric absorbing, emitting, and scattering medium, given the measured incident radiative heat fluxes. The finite-volume method was employed to solve a direct radiative transfer equation for a two-dimensional axisymmetric geometry. Various parameter estimators, such as conjugate-gradient method, hybrid genetic algorithm, and finite-difference Newton method, were employed to solve the inverse problems, while discussing their performances in terms of estimation accuracy and computational efficiency. Based on this, we proposed, as a best inverse analysis tool, a new combined method that adopted the hybrid genetic algorithm as an initial value selector and used the finite-difference Newton method as a parameter estimator.     

二维平板通道中流动与传热的格子-Boltzmann模拟

格子-Boltzmann数值模拟方法(LBM),在最近十几年来得到迅速发展。本文发展了LBM的流动与传热模型,并对二维平板通道中的流动与传热进行了模拟,采用密度分布函数得到速度场,用单独的内能分布函数得到温度场,并与传统FVM方法所得到的多个特征量结果进行了比较,模拟结果与FVM解均吻合很好。鉴于LBM边界条件处理简单和易于实施等特点,该方法可望成为求解流动与传热的一种有效数值模拟手段。     

On the finite volume method and the discrete ordinates method regarding radiative heat transfer in acute forward anisotropic scattering media

The ability of the finite volume method (FVM) and the discrete ordinates method (DOM) to model radiative heat transfer in acute forward anisotropic scattering media has been investigated. The test case involves a purely scattering medium in a cubic enclosure, irradiated by one boundary with diffuse emission. Four phase functions have been considered: three of the Henyey-Greenstein type with respective asymmetry factors of 0.2, 0.8 and 0.93, and a Mie phase function with a strong forward scattering peak (computed for a size parameter of 245 and corresponding to an asymmetry factor of 0.93). Results obtained with the FVM are in good agreement with Monte Carlo reference solutions, whatever the level of acute anisotropic scattering (for asymmetry factors up to 0.93). The DOM combined with the renormalization procedures of the phase function proposed by Kim and Lee (Effect of anisotropic scattering on radiative heat transfer in two-dimensional rectangular enclosures. Int J Heat Mass Transfer 1988;31:1711-21. [1]) and Wiscombe (On initialization error and flux conservation in the doubling method. JQSRT 1976;18:637-58. [2]) provides accurate results only for the smallest asymmetry factor. As the asymmetry factor increases, the renormalization procedures induce strong modifications in the values of the discretized phase function resulting in an underestimation of the effective attenuation by scattering. This error has been found to increase with optical thickness. In fact, when using the DOM, results would be more accurate combining this method with a Delta-Eddington approximation of the phase function, instead of using the actual phase function which is altered too much by renormalization.     

Estimation of the optical properties of seawater from measurements of exit radiance

An inverse analysis for simultaneous estimation of the radiation phase function, single scattering albedo and optical thickness in natural waters, from the knowledge of the exit radiance measurements, is presented. A forward and an inverse model are utilized in our analysis. The forward model uses an analytical discrete-ordinates method for solving the radiative transfer equation and the inverse model contains an algorithm for least-squares estimation that is iteratively solved for retrieving the desired optical properties. The experimental data are simulated with synthetic data corrupted with noise. The results show that the optical properties, with the exception of the optical thickness, can be recovered with high accuracy, even for data with up to 10% noise.     

用改进的遗传算法和高斯牛顿法联合反演三维地下水流模型参数

遗传算法在处理非线性优化问题时具有较好的全局搜索性能,但在局部搜索时搜索效率不高,解的精度亦不高,高斯牛顿法在处理非线性优化问题时的性质正好和遗传算法相反,利用遗传算法和高斯牛顿法的优点,用改进的遗传算法和高斯牛顿法联合反演地下水数值模型参数.首先用遗传算法求出地下水模型参数的初值,然后利用这组初值用高斯牛顿法进行数值模型参数的反演,并以一非均质各向同性三维承压非稳定流理想模型为例,结合有限元法讨论了用遗传算法和高斯牛顿法联合反演地下水数值模型参数的过程.计算结果表明,联合参数反演方法,具有收敛速度快、解的精度高的特点,在地下水渗流和水资源评价等领域可广泛应用.     

Combined spectroscopic method for determining the fluorophore concentration in highly scattering media

The spectroscopic analysis of composition of highly scattering media, in particular, biological tissues, is associated with difficult interpretation of the spectrum of radiation subjected to both absorption and multiple scattering in a medium. The proposed technique of the combined analysis of spectra of laser radiation fluorescence and diffuse reflection by highly scattering biological media allows quantitative determination of concentrations of basic chromophores and fluorophores in tissues with simultaneous estimation of structural changes in objects under study. To determine the effect of structural parameters on the intensity of backscattered laser and fluorescent radiation, numerical simulation and physical modeling of the interaction of exciting and fluorescent radiation with highly scattering media are performed. The dependences of the measured signal on the concentrations of scatterers and fluorophore (protoporphyrin IX) under study are measured.     

基于瞬态辐射信号的各向同性散射均匀介质光学特性PSO反演

利用激光入射半透明介质时有限体积法（FVM）得到的时域反射和透射信号,结合标准微粒群（PSO）、随机微粒群（SPSO）和多相微粒群（MPPSO）三种不同的智能算法反演短脉冲激光辐照下各项同性散射均匀参与性介质的吸收系数和散射系数,结果表明三种微粒群算法都能够高效准确地反演均匀介质内的光学参数。     

Analysis of radiative transport in a cylindrical enclosure—An application of the modified discrete ordinate method

Application of the modified discrete ordinate method (MDOM) proposed by Mishra et al. [Mishra SC, Roy HK, Misra N. Discrete ordinate method with a new and simple quadrature scheme. J Quant Spectrosc Radiat Transfer 2006;101:249-262.] has been extended for calculation of volumetric radiative information in a cylindrical enclosure. Radiatively, the medium inside a diffuse gray 1-D concentric cylinder is absorbing, emitting and scattering. Three types of problems, viz., an isothermal medium representing non-radiative equilibrium case, a non-isothermal medium representing radiative equilibrium situation and the case of a combined mode conduction and radiation heat transfer have been used to test the robustness of the MDOM. Temperature/emissive power and heat flux/energy flow rate distributions in the medium have been found for the effects of various parameters like the extinction coefficient, the scattering albedo, the boundary emissivity and the conduction-radiation parameter. To check the accuracy of the results of the MDOM, results have been compared with those available in the literature and also by obtaining the radiative information using the finite volume method. MDOM has been found to provide accurate results.     

Inverse radiation problem of temperature distribution in one-dimensional isotropically scattering participating slab with variable refractive index

In this paper, an inverse analysis is performed for estimation of source term distribution from the measured exit radiation intensities at the boundary surfaces in a one-dimensional absorbing, emitting and isotropically scattering medium between two parallel plates with variable refractive index. The variation of refractive index is assumed to be linear. The radiative transfer equation is solved by the constant quadrature discrete ordinate method. The inverse problem is formulated as an optimization problem for minimizing an objective function which is expressed as the sum of square deviations between measured and estimated exit radiation intensities at boundary surfaces. The conjugate gradient method is used to solve the inverse problem through an iterative procedure. The effects of various variables on source estimation are investigated such as type of source function, errors in the measured data and system parameters, gradient of refractive index across the medium, optical thickness, single scattering albedo and boundary emissivities. The results show that in the case of noisy input data, variation of system parameters may affect the inverse solution, especially at high error values in the measured data. The error in measured data plays more important role than the error in radiative system parameters except the refractive index distribution; however the accuracy of source estimation is very sensitive toward error in refractive index distribution. Therefore, refractive index distribution and measured exit intensities should be measured accurately with a limited error bound, in order to have an accurate estimation of source term in a graded index medium.     

Analysis of transport of collimated radiation in a participating media using the lattice Boltzmann method

Application of the lattice Boltzmann method (LBM) recently proposed by Asinari et al. [Asinari P, Mishra SC, Borchiellini R. A lattice Boltzmann formulation to the analysis of radiative heat transfer problems in a participating medium. Numer Heat Transfer B 2010; 57:126–146] is extended to the analysis of transport of collimated radiation in a planar participating medium. To deal with azimuthally symmetric radiation in planar medium, a new lattice structure for the LBM is used. The transport of the collimated component in the medium is analysed by two different, viz., flux splitting and direct approaches. For different angles of incidence of the collimated radiation, the LBM formulation is tested for the effects of the extinction coefficient, the anisotropy factor, and the boundary emissivities on heat flux and emissive power distributions. Results are compared with the benchmark results obtained using the finite volume method. Both the approaches in LBM provide accurate results.     

On a method for solving the inverse problem in potential scattering

A method for solving the inverse problem in the non-relativistic elastic scattering theory, using the analytic and asymptotic properties of the scattering amplitude is proposed and the influence of the discontinuity parameters of the scattering amplitude on the properties of the resulting potentials is discussed. The case with spherically symmetric forces and without bound states is considered. The possibility for solving the inverse problem by this method, leading to the singular repulsive potentials is mentioned.     

凝固前沿和气泡相互作用的大密度比格子玻尔兹曼方法模拟

建立了二维双组分两相流的大密度比格子玻尔兹曼方法 (lattice Boltzmann method, LBM)模型. 该模型基于改进的Shan-Chen伪势多相流LBM模型, 结合采用不同时间步长的方法, 实现密度比达800以上的气液两相流模拟. 为了对模型进行验证, 模拟了在不同气液相互作用系数和密度比条件下气泡内外压力差与其半径之间的关系, 其结果满足Laplace定律. 将所建立的大密度比LBM与介观尺度的元胞自动机(cellular automaton, CA)和有限差分法(FDM)相耦合, 用LBM模拟气液两相流, 用CA方法模拟固相生长, 用有限差分法模拟温度场, 采用LBM-CA-FDM耦合模型对定向凝固过程中凝固前沿的气泡与液-固界面之间的相互作用进行模拟研究. 结果表明, 绝热气泡的存在影响了温度场分布, 使得凝固前沿接近气泡时, 液-固界面凸起, 在不同的固相生长速度条件下, 出现凝固前沿淹没气泡或气泡脱离凝固前沿的不同情况, 模拟结果与实验结果符合良好. 关键词： 格子玻尔兹曼方法 元胞自动机 凝固 气泡