吸收散射性三维矩形介质内辐射源项的反问题

提出了一种由边界出射辐射强度反演吸收散射性三维矩形介质内辐射源项分布的方法。该方法是在辐射传递方程离散坐标近似的基础上,用求目标函数极小值的共轭梯度法进行反演计算。通过对介质辐射特性、光学厚度等参数对反演精度影响的分析,结果表明,即使存在测量误差,本文所提出的方法可较精确地反演辐射源项。     

半透明介质场参数的非接触测量算法研究

本文在考虑半透明介质吸收、发射以及散射衰减作用的情况下采用视在光线法计算某探测方向上介质边界出射辐射强度分布作为逆问题的输入数据,采用LSQR方法在介质辐射特性参数已知情况下重建介质三维温度分布,重建结果显示,无论有无测量误差存在,利用LSQR方法都可以在已知辐射特性参数的条件下很好地重建出介质的三维温度场。在此基础上提出LSQR-SPSO混合算法,并应用于同时重建三维温度场及辐射特性参数(吸收系数、散射系数)。计算结果显示,无论有无测量误差存在,利用LSQR-SPSO混合算法都可以很好地同时重建出介质的吸收系数、散射系数及三维温度场,相比较而言,介质的温度场更容易被重建出来。     

基于频域辐射传输模型的光学成像研究

本文以超短脉冲激光照射参与性介质的光学成像为研究背景,分别构建了短脉冲激光在参与性介质内的频域辐射传输正问题模型和根据边界探测所得频域信号重建介质内部光学参数的逆问题模型。在瞬态辐射传输方程的基础上,利用傅里叶变换得到频域辐射传输方程,采用有限体积法求解频域传输方程,模拟超短脉冲激光在二维参与性介质内传输的过程,得到介质边界的出射频域辐射信号。选取共轭梯度法作为反演算法,采用伴随差分模型求解目标函数梯度,重建了二维非均匀参与性介质内不同位置内含物的光学参数分布。结果表明,基于频域辐射传输方程的伴随差分模型能够较为准确地反演多维参与性介质内的光学参数。     

用射线踪迹法解黑表面矩形介质耦合换热

本文用射线踪迹-节点分析法研究了二维黑体表面矩形、各向同性散射半透明介质内辐射与导热瞬态耦合换热。采用全隐格式的有限差分法离散二维瞬态微分能量方程,用辐射传递系数来表示辐射源项,结合谱带模型并采用射线踪迹法求解辐射传递系数。采用Patankar线性化方法将辐射源项及不透明边界条件线性化,并采用附加源项法处理边界条件,运用ADI方法求解名以上的线性化方程组,从而解得二维矩形介质内的瞬态温度分布。     

激光脉冲在各向异性散射介质内的瞬态热效应

本文考察了激光脉冲在吸收、发射、各向异性散射介质内引起的瞬态热效应。将激光脉冲辐射在介质内的传递过程分为两个子过程：发射－衰减－反射过程和吸收－散射过程。用光线踪迹法结合节点分析导出辐射传递系数和辐射源项,用控制容积法解瞬态能量方程。检验结果表明,本文的计算方法准确。在此基础上,考察了散射特性、初始温度对激光脉冲响应的影响。     

辐射模型和弥散效应对多孔介质内燃烧影响的数值模拟

本文根据气固两相局部非热平衡假设,建立了甲烷/空气预混气在惰性多孔介质内的一维层流燃烧数学模型。分别采用附加导热、Rossland模型和二通量法模型求解固体能量方程中的辐射源项,研究了热辐射模型和弥散效应对多孔介质内燃烧火焰结构的影响。结果表明,多孔固体表面辐射的影响不可忽略,辐射模型对火焰温度的预测影响显著,而弥散效应对气体温度的分布及反应热影响则较小。     

基于时域光信息的近红外光学成像

利用时域辐射传输方程,模拟了脉冲激光在弥散介质内的传输。利用广义Guess-Markov随机场模型,构建了正则化项,克服了反问题的病态特性。将步长加速法与序列二次规划算法相结合来求解反问题,构建了一种不依赖初值的光学成像算法,重建了二维弥散介质内的吸收系数和散射系数的分布,并模拟了癌变生物组织的近红外光学成像。研究结果表明,所提算法能够较好地解决初值依赖性问题且具备高精度,重建得到的图像能够清晰地分辨出介质内部的结构。     

表面不透明三层漫反射散射性介质耦合换热

本文采用射线踪迹、节点分析法研究了三层吸收、各向同性散射性介质层内的一维辐射和导热瞬态耦合换热,复合层表面不透明漫反射,介质层交界面半透明漫反射,且半透明漫反射交界面的反射率采用Fresnel反射定律确定。采用一层和二层辐射能量传递模型跟踪辐射能量在三层介质内的传递,从而推导出辐射传递系数。运用辐射传递系数求解辐射源项,在辐射对流边界条件下、采用全隐格式求解瞬态能量方程,并从机理上研究了辐射和导热耦合换热过程。     

三维燃烧介质和壁面温度的非接触联合重建研究

温度分布在线实时测量对于燃烧过程优化和污染物控制具有重要意义, 针对以往非接触三维温度分布重建过程的耗时性问题和忽略壁面辐射的不足, 本文提出了一种新的离散重建模型, 用于三维吸收、 发射和散射性高温燃烧介质以及壁面温度的快速联合非接触测量. 该模型以四个CCD(Charge Coupled Device) 为测量传感器, 通过构建辐射逆问题求解方程, 从CCD输出的辐射投影图像重建温度分布. 介质中不同投影方向内的辐射传递过程通过离散传递法来描述, 介质的散射和壁面反射则通过离散坐标法来近似. 离散后计算局部辐射强度的病态方程通过最小二乘余量法来求解, 论文对其计算速度进行了优化. 通过非对称温度分布测量算例分析了该模型的有效性, 讨论了测量噪音、 介质和壁面辐射特性对重建精度的影响, 并与其他方法对比分析了模型的重建速度. 计算结果表明本文提出的离散模型可以有效地用于大型高温燃烧介质和壁面温度分布的联合非接触测量. 即使在有噪声的情况下, 该模型也能获得准确的测量结果, 与其他计算方法相比, 采用改进的最小二乘余量法, 能有效地提高温度分布的重建计算速度.     

各向异性介质辐射特性参数联合反演

本文基于BP神经网络方法结合蒙特卡洛和BEER定律辐射传输模拟方法建立了联合反演各向异性散射介质的辐射特性参数模型。首先采用半球透射率结合半球反射率反演模型反演了各向同性介质的吸收系数和散射系数,在此基础上增加准直透射率,建立了联合反演各向异性介质的吸收系数、散射系数和散射不对称因子三参数联合反演模型。反演结果表明该模型能准确反演出介质辐射特性参数,具有实用意义。此外,为了检验测量误差对模型的反演准确性的影响,分别在不同程度测量误差情况下进行反演,结果显示测量误差对散射不对称因子反演值影响较大。     

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.     

An inverse source problem in radiative transfer

The spherical-harmonics method is used to develop a solution to an inverse source problem in radiative transfer. It is assumed that, with the exception of the inhomogeneous source term, all aspects of the radiation-transport problem are known, and we seek to determine the inhomogeneous source term from specified angular distributions of radiation exiting the two surfaces of a homogeneous plane-parallel medium. Anisotropic scattering is included in the monochromatic radiative-transfer model and general reflecting boundary conditions are considered.     

Application of multi-phase particle swarm optimization technique to inverse radiation problem

The multi-phase particle swarm optimization (MPPSO) technique is applied to the inverse radiation problem in the present paper. The directional radiative intensities are served as the measurement data to estimate the radiative source term, optical thickness, scattering albedo, and phase function in one-dimensional semitransparent plane-parallel media by the inverse simulation. To check the performance and accuracy in retrieval, a comparison is presented between three PSO methods, i.e. the MPPSO, the standard PSO, and the Stochastic PSO. The results confirm the potential of the proposed approach MPPSO and show its effectiveness and superiority over the other two PSO algorithms. Furthermore, the effects of swarm size, searching space, phase change frequency, and velocity-reinitializing frequency on the convergence velocity and computational accuracy of MPPSO are also investigated.     

气粒混合物非灰辐射特性合并宽窄谱带K分布模型

气粒混合物辐射问题具有全场性、非灰性、耦合性等特点,准确预估高温燃气/粒子非灰辐射特性是非常重要的。本文将合并宽窄谱带K分布模础(CWNBCK)与离散坐标法(DOM)结合,开展了非灰气粒混合物辐射换热问题的模拟工作,分别验证了一维和三维情况下应用该模型的准确性,给出不同工况下的热流源项、壁面热流或辐射热流等。结果表明:该模型能够给出与SNB模型精度基本相同的结果,考虑其计算效率的提高,可以在工程实际中应用该模型计算非灰气粒混合物辐射换热。     

An efficient and robust reconstruction method for optical tomography with the time-domain radiative transfer equation

An efficient and robust method based on the complex-variable-differentiation method (CVDM) is proposed to reconstruct the distribution of optical parameters in two-dimensional participating media. An upwind-difference discrete-ordinate formulation of the time-domain radiative transfer equation is well established and used as forward model. The regularization term using generalized Gaussian Markov random field model is added in the objective function to overcome the ill-posed nature of the radiative inverse problem. The multi-start conjugate gradient method was utilized to accelerate the convergence speed of the inverse procedure. To obtain an accurate result and avoid the cumbersome formula of adjoint differentiation model, the CVDM was employed to calculate the gradient of objective function with respect to the optical parameters. All the simulation results show that the CVDM is efficient and robust for the reconstruction of optical parameters.     

An Inverse Source Problem with Sparsity Constraint for the Time-Fractional Diffusion Equation

In this paper, an inverse source problem for the time-fractional diffusion equation is investigated. The observational data is on the final time and the source term is assumed to be temporally independent and with a sparse structure. Here the sparsity is understood with respect to the pixel basis, i.e., the source has a small support. By an elastic-net regularization method, this inverse source problem is formulated into an optimization problem and a semismooth Newton (SSN) algorithm is developed to solve it. A discretization strategy is applied in the numerical realization. Several one- and two- dimensional numerical examples illustrate the efficiency of the proposed method.     

Transient coupled heat transfer in a rectangular medium with black surfaces

The main objective of this paper is to extend to two-dimensional (2-D) medium the ray tracing-node analyzing method, which has already been successfully used to solve one-dimensional (1-D) problem of coupled heat transfer in a semitransparent medium. For simplicity, an infinitely long rectangular semitransparent medium with four black opaque surfaces is chosen as our studying object. A control volume method in the implicit scheme is adopted for discretizing the partial transient energy equation. In combination with spectral band model, the radiative heat source term is calculated using the radiative transfer coefficients (RTCs), which are deduced by the ray tracing method. The Partankar's linearization method is used to linearize the radiative source term and the opaque boundary condition, and the linearized equations are solved by the ADI method. Effects of absorption coefficient, refractive index and conductivity on transient cooling process in the 2-D gray rectangular medium are investigated under the condition that the radiation and convection processes cool one side of the rectangular medium while heat the remaining three sides.     

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.