Simultaneous estimation of the temperature and heat rate distributions within the combustion region by a new inverse radiation analysis

A new inverse radiation analysis is presented for estimating the heat rate and temperature distributions in the combustion region from the information of the temperature and heat flux profiles of wall elements in the system. The Monte Carlo method is employed to solve the radiative heat transfer equation. The inverse radiation problem is posed as a minimization problem of the least squares criterion, which is solved by the conjugate gradient method. The performance of the present technique of inverse analysis is evaluated and the effects of the errors of the absorption coefficient, emissivity and convective heat transfer coefficient on the inverse analysis are investigated. The results show that the present technique is robust and yields accurate estimation even with noisy measurement.     

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

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

Inverse estimation of the pulse parameters of a time-varying laser pulse to obtain desired temperature at the material surface

A proper selection of the pulse parameters is essential to achieve desired temperature at the material surface. This leads to obtain the required metallurgical changes in the surface vicinity when a time-varying laser pulse is used in a heating process such as surface modification. In this paper, the conjugate gradient method (CGM) for parameter estimation is successfully applied to estimate the unknown laser pulse parameters for those purposes during laser heating process. The determination of the pulse parameters is treated as a one-dimensional, transient, non-linear inverse heat conduction problem (IHCP). Based on a sensitivity analysis, the inverse problem is solved as an optimization problem comparing a desired temperature at the surface and a calculated one where the objective function is minimized by CGM. The method has been applied to a test case of a heating process on steel, appropriate pulse parameters and desired temperature distribution can also be returned.     

共轭梯度法求解双曲传热多宗量反问题

提出双曲传热反问题热物性参数和边界条件多宗量联合反演的一般数值求解模式,考虑了非均质和分布参数的影响,时域上采用时域精细算法进行离散,建立了便于敏度分析的有限元正演模型.由最小二乘原理建立反演模型,应用共轭梯度法进行求解.探讨了时间步长和测量误差对反演结果的影响,并进行了数值验证.     

Regularized variable metric method versus the conjugate gradient method in solution of radiative boundary design problem

In this paper, an appropriate distribution of the heating elements’ strengths in a radiation furnace is estimated using inverse methods so that a pre-specified temperature and heat flux distribution is attained on the design surface. Minimization of the sum of the squares of the error function is performed using the variable metric method (VMM), and the results are compared with those obtained by the conjugate gradient method (CGM) established previously in the literature. It is shown via test cases and a well-founded validation procedure that the VMM, when using a “regularized” estimator, is more accurate and is able to reach at a higher quality final solution as compared to the CGM. The test cases used in this study were two-dimensional furnaces filled with an absorbing, emitting, and scattering gas.     

一维辐射系统吸收系数的反问题

辐射反问题是当前的研究热点之一。本文采用离散坐标法求解正问题,用共轭梯度法求解一维系统吸收系数的反问题。用正问题的计算结果加上误差作为反问题的求解条件。数值计算的结果表明,共轭梯度法可以比较好的反演辐射系统的吸收系数。在考虑散射时,反演效果比不考虑散射的反演效果好。     

一维梯度折射率介质光学特性的反演

在单个方波脉冲入射情况下,利用共轭梯度法对一维梯度折射率介质的折射率、吸收系数以及散射系数进行了反演.正问题采用间断有限元法求解,反问题的解则在正问题的基础上通过共轭梯度法得到.研究结果表明,利用单个方波脉冲入射情况下的时域半球反射率以及时域半球透射率作为测试值能够有效地反演一维梯度折射率介质的光学特性及其分布情况.     

炉膛三维温度场重建中Tikhonov正则化和截断奇异值分解算法比较

在煤粉锅炉诊断中火焰辐射能图像扮演着越来越重要的角色, 通过电荷耦合器件(CCD)获得的辐射能图像可以重建出炉内火焰三维温度场, CCD 用于获取视场角内的辐射能图像. 温度场重建的矩阵方程是一个严重病态的方程, 本文使用两种算法(Tikhonov正则化算法和截断奇异值分解(TSVD)算法)来重建温度场. 应用广义交叉检验算法来选取正确的正则化参数. 数值模拟的环境为一个10 m×10 m×10 m的三维炉膛, 系统被划分为10×10×10的1000个网格, 每个网格单元都是边长为1 m的立方体. 在正问题求解所得到的CCD接受信号基础上加上不同随机误差以模拟测量时的CCD接受信号. 研究两种算法重建后的温度重建误差、两者的重建时间, 以及最高温度的重建效果. 初步的研究结果显示, 一般情况下基于Tikhonov算法重建的温度场比基于TSVD算法重建的温度场误差要小, 计算所需时间短, 最高温度重建更准确.     

Estimating thermal transport in deep X-ray lithography with an inversion method

This study proposes a general methodology for estimating the depth profile of the heat source of the thermal transport system during deep X-ray lithography. The exposure process in a lithography system is considered as an inverse heat conduction problem with an unknown heat source. The conjugate gradient method is used to solve the inverse problem. Numerical results confirm that the method proposed herein can accurately estimate the heat source even involving the inevitable measurement errors. Furthermore, this methodology can also be applied to estimate the local distribution of temperatures when using scanning thermal microscopy (SThM) to microthermally machine materials and will contribute to increase the quality of microthermally machined products. In addition, a thermomechanical data-storage system, which utilizes a resistively heated atomic-force-microscopy (AFM) cantilever tip to read and write data bits, can also adopt this inverse methodology to control the temperature of a polymer substrate.     

Estimation of thermal contact resistance and thermally induced optical effects in single-coated optical fibers

In this study, a conjugate gradient method based on an inverse algorithm is applied to estimate the unknown time-dependent thermal contact resistance in a single-coated optical fiber, which is subjected to transient thermal loading. While knowing the temperature history at the measuring position, no prior information is needed on the functional form of the unknown contact resistance. The temperature data obtained from the direct problem are used to simulate the temperature measurement. The influence of measurement errors, initial guess values, and measurement locations upon the precision of the estimated results is also investigated. Results show that an excellent estimation on the time-dependent thermal contact resistance, temperature distributions, thermally induced microbending loss, and refractive index changes can be obtained for the case considered in this study.     

管道内壁侵蚀形状识别的无网格法研究

本文使用直接配点无网格法结合共轭梯度法对管道内壁面侵蚀状况进行了反演识别。直接配点无网格法使用节点离散求解区域,采用移动最小二乘近似构造试函数,直接配点法构造线性方程组进行导热正问题求解;反演过程采用共轭梯度法使目标函数最小化,Akima三次样条插值将连续的几何边界反演问题转化为离散点几何位置的反演,并最终将这些离散点拟合成为光滑曲线。文中选择两个典型算例对数值方法进行验证,模拟结果表明使用直接配点无网格法结合共轭梯度法进行管道内壁几何边界识别具有较高精度。     

基于红外测温的试件内部缺陷的识别算法研究

内部缺陷形状位置的识别是传热逆问题的重要研究课题.本文对带有内部缺陷的试件建立了二维传热模型,通过有限元法求解了在一定加热条件下试件表面的温度分布;并通过共轭梯度法,提出了根据表面所测温度分布识别试件内部缺陷形状、位置的计算方法.数值实验证明了算法的有效性.     

A numerical procedure for estimation of the melt depth in laser material processing

In this study, an algorithm based on conjugate gradient method (CGM) is applied to estimate the unknown time dependent melt depth during laser material processing in liquid phase. The determination of the melt depth is treated as a one-dimensional, transient, inverse heat conduction problem (IHCP). It is assumed that no prior information is available for the functional form of the unknown melt depth, but it can be estimated by an inverse analysis with temperature measurements near the heated surface. The algorithm has been applied to aluminum, titanium and fused quartz and accurate melting depth and temperature distributions can also be returned. In addition, this methodology can also be applied to solve other problems such as calculating the cutting forces in nanomachining by atomic force microscopy (AFM), and estimating the heat sources in a X-ray lithographic process.     

ESTIMATED TEMPERATURE ON A MACHINED SURFACE USING AN INVERSE APPROACH

An experiment devoted to the heat flux estimation in a workpiece during a machining process by turning is presented. The method is based on temperature measurement from thermocouples embedded in the workpiece, close to the heated surface. A model that expresses the heat flux according to the temperature at the sensors is developed. The stationary and linearity assumptions are used in order to decompose the three-dimensional original problem into two bi-dimensional problems. This decomposition can be realized given the difference between the cutting speed and feed velocity in two orthogonal directions. The temperature on the machined surface is calculated from the estimated heat flux and the heat transfer model in the workpiece. The application concerns hard steel machining, using a CBN insert tool. Three parameters are placed into evidence from this application: the temperature magnitude on the machined surface, the thermal gradient in the workpiece, and the `thermal persistence' that represent the heating time of the machined surface. This study leads to a better understanding of the influence of temperature during a hard steel turning process.     

半透明介质中辐射传递方程的反演计算及数值模拟

本文介绍了一侧为半透明、另一侧为非透明界面时一维半透明介质的辐射强度计算式。采用辐射与导热复合换热模型计算半透明介质内温度场。利用已知的温度场求半透明介质的辐射强度一正问题计算。将此辐射强度代入辐射反问题计算模型，引入测量误差，采用Ｃｈａｈｉｎｅ方法及演半透明介质内温度场一反问题计算。数值模拟表明，本文所采用的辐射反演法具有较高的精度及稳定性。     

Inverse radiation problem in one-dimensional slab by time-resolved reflected and transmitted signals

A time-domain inverse approach is proposed for estimating the distribution of absorbing and scattering coefficients in one-dimensional inhomogeneous media. The temporal reflected and transmitted signals are detected when an ultra-short pulse irradiates on the boundary of semi-transparent scattering media. Forward computation and inverse algorithm employ the least-squares finite element method and conjugate gradient method, respectively. As the prevalent diffusion approximation is not employed in our model, the present approach can be extended to more comprehensive application. The investigation about detected signals indicates that the reflected signals play a significant role in reconstructing optical properties; the signals in early sampling time are more important than those at long-time logarithm slope, and so, more attention should be paid to the early signals in the solution of inverse radiation problem. Three different inverse radiation problems are investigated to show the ability of the present approach to deal with the two-layer, three-layer and continuous inhomogeneous media. The effect of measured errors on the accuracy of reconstruction is investigated by adding artificial random errors. The results indicate that accurate reconstruction depends on not only precise numerical simulation but also quality of detected data.     

Acoustical inverse problems regularization: Direct definition of filter factors using Signal-to-Noise Ratio

Acoustic imaging aims at localization and characterization of sound sources using microphone arrays. In this paper a new regularization method for acoustic imaging by inverse approach is proposed. The method first relies on the singular value decomposition of the plant matrix and on the projection of the measured data on the corresponding singular vectors. In place of regularization using classical methods such as truncated singular value decomposition and Tikhonov regularization, the proposed method involves the direct definition of the filter factors on the basis of a thresholding operation, defined from the estimated measurement noise. The thresholding operation is achieved using modified filter functions. The originality of the approach is to propose the definition of a filter factor which provides more damping to the singular components dominated by noise than that given by the Tikhonov filter. This has the advantage of potentially simplifying the selection of the best regularization amount in inverse problems. Theoretical results show that this method is comparatively more accurate than Tikhonov regularization and truncated singular value decomposition.     

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.     

Modelling of the solid–liquid interface during laser processing using an inverse methodology

This study presents a methodology for estimating the melt depth during laser processing of solid materials. The determination of the melt depth is treated as an inverse heat conduction problem, which includes the solid and liquid phases. The conjugate gradient method is applied to treat the inverse problem using the available temperature measurements. Without the inverse methodology the melt depth is very difficult to obtain with precision. The proposed method can also be applied during microthermal machining to determine the location of the solid–liquid interface and the temperature distributions of the two phases by using scanning thermal microscopy.     

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.