Heat source reconstruction from noisy temperature fields using a gradient anisotropic diffusion filter

This paper presents a post-processing technique for noisy temperature maps based on a gradient anisotropic diffusion (GAD) filter in the context of heat source reconstruction. The aim is to reconstruct heat source maps from temperature maps measured using infrared (IR) thermography. Synthetic temperature fields corrupted by added noise are first considered. The GAD filter, which relies on a diffusion process, is optimized to retrieve as well as possible a heat source concentration in a two-dimensional plate. The influence of the dimensions and the intensity of the heat source concentration are discussed. The results obtained are also compared with two other types of filters: averaging filter and Gaussian derivative filter. The second part of this study presents an application for experimental temperature maps measured with an IR camera. The results demonstrate the relevancy of the GAD filter in extracting heat sources from noisy temperature fields.     

红外无损探测中多宗量多热源反演问题的研究

为研究红外无损探测稳态多热源反演逆问题, 建立不同形状的均质与非均质稳态热传导模型, 其中内热源个数、位置、强度、面积均为未知项. 基于数值算法中有限元算法对模型进行离散分析, 化简有限元矩阵方程, 最终转化为对Ax = b高度欠定方程的求解. 首次利用分段多项式谱截断奇异值分解法处理内热源逆问题, 并对算法进行改进, 有效改善了该算法在处理多热源反演时存在的严重的热源叠加效应. 根据反演出的内热源信息, 利用有限元算法计算重构出整个模型内所有节点的温度分布. 运用数值仿真Comsol软件和具体实物实验对算法进行有效性评估, 并验证算法在不同热传导模型中的表现. 结果表明, 算法能够准确反演出多热源各参量信息, 在非均质材料模型中仍能准确地反演出热源项, 并有效重构出模型内温度场. 该算法可应用于材料无损检测及人体红外医学成像等领域.     

Independent component analysis in the presence of noise in fMRI

A noisy version of independent component analysis (noisy ICA) is applied to simulated and real functional magnetic resonance imaging (fMRI) data. The noise covariance is explicitly modeled by an autoregressive (AR) model of order 1. The unmixing matrix of the data is determined using a variant of the FastICA algorithm based on Gaussian moments. The sources are estimated using the principle of maximum likelihood by modeling the source densities as asymmetric exponential functions. Effect of dimensionality reduction on the effective noise covariance used, accuracy of the obtained mixing matrix and degree of improvement in estimating fMRI sources are investigated. The primary conclusions after using this method of evaluation are as follows: (a) weighting matrix estimates are similar for noisy and conventional ICA in the realm of typical fMRI data, and (b) source estimates are improved by 5% (as measured by the correlation coefficient) in realistic simulated data by explicitly modeling the source densities and the noise, even when just a simple white noise model is used.     

Dimensional scaling of flame propagation in discrete particulate clouds

The critical dimension necessary for a flame to propagate in suspensions of fuel particles in oxidiser is studied analytically and numerically. Two types of models are considered: First, a continuum model, wherein the individual particulate sources are not resolved and the heat release is assumed spatially uniform, is solved via conventional finite difference techniques. Second, a discrete source model, wherein the heat diffusion from individual sources is modelled via superposition of the Green's function of each source, is employed to examine the influence of the random, discrete nature of the media. Heat transfer to cold, isothermal walls and to a layer of inert gas surrounding the reactive medium are considered as the loss mechanisms. Both cylindrical and rectangular (slab) geometries of the reactive medium are considered, and the flame speed is measured as a function of the diameter and thickness of the domains, respectively. In the continuum model with inert gas confinement, a universal scaling of critical diameter to critical thickness near 2:1 is found. In the discrete source model, as the time scale of heat release of the sources is made small compared to the interparticle diffusion time, the geometric scaling between cylinders and slabs exhibits values greater than 2:1. The ability of the flame in the discrete regime to propagate in thinner slabs than predicted by continuum scaling is attributed to the flame being able to exploit local fluctuations in concentration across the slab to sustain propagation. As the heat release time of the sources is increased, the discrete source model reverts back to results consistent with the continuum model. Implications of these results for experiments are discussed.     

Modeling Chemically Reactive Flow of Sutterby Nanofluid by a Rotating Disk in Presence of Heat Generation/Absorption

In this article we investigate the flow of Sutterby liquid due to rotating stretchable disk. Mass and heat transport are analyzed through Brownian diffusion and thermophoresis. Further the effects of magnetic field, chemical reaction and heat source are also accounted. We employ transformation procedure to obtain a system of nonlinear ODE's. This system is numerically solved by Built-in-Shooting method. Impacts of different involved parameter on velocity, temperature and concentration are described. Velocity, concentration and temperature gradients are numerically computed. Obtained results show that velocity is reduced through material parameter. Temperature and concentration are enhanced with thermophoresis parameter.     

Heat flow in interference filters

Presented here is a linear theory of heat conduction in an interference filter with an absorbing spacer. The source of heat is a Gaussian beam. The effects of thermal properties and construction of the system on the temperature are studied. A comparison with and conditions of validity of the single-film theory are also discussed.     

Ultrasonic imaging in air using fan-beam tomography and electrostatic transducers

Air-coupled ultrasonic capacitance transducers operating at frequencies of up to 1 MHz have been employed in a fan-beam configuration for the cross-sectional tomographic imaging of temperature fields and flow fields in air, and the location of solid objects. Separate transmitter and receiver transducers were manufactured using thin polymer dielectric membranes and polished metal backplates, and used to acquire through-transmission data. The fan-beam reconstruction was developed in LabVIEW using a re-bin routine combined with a filtered backprojection algorithm and a difference technique to generate the cross-sectional images. The system was first used to reconstruct images showing the locations of solid objects positioned within the scanned region through interpretation of the arrival time of the transmitted ultrasound. The technique was then extended to image the temperature fields produced in air above a small heat source and the flow field produced by a nozzle connected to a regulated compressed air source. Reconstructed temperatures were within 4% of the measured background air temperature and 9% of the air temperature measured above the heat source. Reconstructed images of the flow field above a small nozzle were also presented, showing that the horizontal component of the flow velocity could be resolved using this method.     

Magnetoelectric Coupling Induced Electric Dipole Glass State in Heisenberg Spin Glass

Multiferroic behavior in an isotropic Heisenberg spin glass with Gaussian random fields, incorporated by magneto-electric coupling derived from the Landau symmetry argument, are investigated. Electric dipole glass transitions at finite temperature, due to coupling, are demonstrated by Monte Carlo simulation. This electric dipole glass state is solely ascribed to the coupling term with chiral symmetry of the magnetization, while the term associated with the spatial derivative of the squared magnetization has no contribution.     

局域热源对磁岛热输运的影响

研究托卡马克等离子体磁岛内的热输运行为.应用局域高斯热源对磁岛加热来模拟电子回旋共振加热.对同时存在背景热源与局域高斯热源的情况,观察局域高斯热源对径向电子温度分布及热输运产生的变化,分析局域高斯热源对温度扰动及磁岛约束能量的影响.     

一种新的叠栅层析迭代算法

利用叠栅层析具有的测量动态范围大,对震动环境不敏感的特点,从其偏折原理的本质出发,提出一种利用代数迭代法来重建三维流场的叠栅层析重建算法。在此基础上进行了双峰函数的模拟实验,分别进行了6方向,12方向以及叠加高斯噪声的12方向的数值重建,并在相同条件下与滤波反投影方法和对偏折投影数据积分的代数迭代算法这两种已有的叠栅层析重建算法进行了比较,同时用该新算法结合属性矩阵对包含遮挡物的模拟场进行重建。对比重建结果,本算法具有较强的抗噪声能力,并且对非完全数据下的层析重建也有较好的处理效果。     

Infrared thermography applied to the resolution of inverse heat conduction problems: recovery of heat line sources and boundary conditions

In this paper we present an application of infrared thermography for inverse heat conduction problems resolution. The approach described in the paper is based on a Boundary Element Method formulation of the transient heat diffusion equation. The inverse problems under investigation concern the time and space reconstruction of unknown boundary conditions or heat line source strength. As there is a lack of information in the system, some additional measurements are necessary to solve the problem. In the examples proposed in this paper the extra information is provided by an infrared scanner. The measurements contained in the infrared pictures are used in the model as a Dirichlet boundary condition or as a special boundary condition prescribing both temperature and heat flux density on the scanned boundary. We present some experimental results concerning line source strength identification and the reconstruction of unknown heat fluxes applied on an out of reach boundary. All the examples presented in this paper are related to 2D transient diffusion. As the inverse problem is ill-posed, time and space regularization techniques are used to stabilize the solution and reduce the sensitivity of the latter to measurement errors.     

Experimental and Numerical Investigation of Mixed-Convection Heat Transfer from Protruding Heat Sources in an Enclosure

Heat transfer from a 5-by-4 array of protruding heat sources in the bottom wall of an adiabatic enclosure is investigated experimentally and numerically. The enclosure has a fixed inlet and different outlets. Flow velocity and applied heat fluxes are varied. Steady-state temperatures of the heat sources, air inlet and outlet temperatures, and enclosure wall temperature are measured. Results indicate that the first row is the coolest, and the Nusselt numbers of the heat sources, which are near the enclosure wall, are higher than those in the inner side. There is a good agreement between the experimental and numerical results.     

轴对称介质内辐射源项反演算法

结合Abel变换和离散坐标法,提出了一种基于CCD相机采集的单幅辐射图像重建轴对称发射-吸收介质内辐射源项分布的反演算法。通过在求解辐射正问题得到的准确值的基础上,添加随机噪声模拟试验测量数据,分析了网格数目、辐射源项分布形式、吸收系数和测量误差对算法反演精度的影响。测试结果表明:该算法对测量误差不敏感,在有测量误差的情况下也能够准确的重建介质内的辐射源项分布。     

Inverse radiation problem of axisymmetric turbulent sooting free flame

An iteration method is extended to reconstruct the time-averaged temperature distribution in turbulent axisymmetric sooting free flame by the multi-wavelength measured data of low time-resolution outgoing emission and transmission radiation intensities. A Gaussian probability density function is used to simulate the turbulent fluctuation of temperature. The reconstruction of time-averaged temperature profile consists of three steps. First, the time-averaged spectral absorption coefficient is retrieved from the time-averaged transmissivity data by iteration method. Then the time-averaged blackbody spectral radiation intensity is estimated from the low time-resolution outgoing spectral emission radiation intensities. Finally, the time-averaged temperature and its standard deviation are approximately reconstructed from the multi-wavelength time-averaged spectral emission radiation data by the least-square method. Both exact and noisy input data have been used to test the performance of the proposed inversion method. The results show that the time-averaged temperature profiles can be estimated with good accuracy by the presented inversion method, even with noisy input data, and the standard deviation of temperature is more sensitive to the measurement errors. In the case of large temperature fluctuation, the errors of estimation for time-averaged temperature profile are large if the turbulent fluctuation is not taken into account.     

Fundamental solutions to time-fractional heat conduction equations in two joint half-lines

Heat conduction in two joint half-lines is considered under the condition of perfect contact, i.e. when the temperatures at the contact point and the heat fluxes through the contact point are the same for both regions. The heat conduction in one half-line is described by the equation with the Caputo time-fractional derivative of order α, whereas heat conduction in another half-line is described by the equation with the time derivative of order β. The fundamental solutions to the first and second Cauchy problems as well as to the source problem are obtained using the Laplace transform with respect to time and the cos-Fourier transform with respect to the spatial coordinate. The fundamental solutions are expressed in terms of the Mittag-Leffler function and the Mainardi function.     

On unsteady-state heat-and-mass transfer in polyatomic gases

The problem of a nonstationary source of heat or particles in a polyatomic gas is solved analytically. The temperature and concentration distributions vs. distance to the source are constructed. It is shown that a nonstationary heat source may generate acoustic vibration. The solution obtained can be helpful in studying acoustic wave propagation, temperature and concentration fields produced by heat or particle sources of various configuration (in particular, by a laser beam in an absorbing media), etc.     

基于二阶导数谱与特征吸收窗的红外光谱定量检测方法

红外光谱分析技术已被液体、固体和气体混合物各组分的定量检测广泛应用。针对红外光谱定量检测过程中存在的计算过程复杂、受背景及噪声干扰、多组分谱线混叠的问题,一方面,通过获取谱线二阶导数去除背景及噪声干扰,并一定程度上区分混叠谱峰。在获取二阶导数谱线过程中运用Savitzky-Golay滤波器进行平滑滤波处理,并根据谱线频率特征选取Savitzky-Golay滤波器最佳滤波参数,解决了滤波参数缺乏标准化选取方法的问题。另一方面,根据被测组分的吸收分布情况生成特征吸收窗分别对待测原始吸收谱线以及二阶导数谱线加以处理,提取出在浓度计算过程中占据重要性更大的特征吸收的区域,从而排除非吸收区域中背景、噪声及其他组分对计算结果的干扰。以一组丙烷、丙烯和甲苯混合气体为例,运用二阶导数方法以及特征吸收窗方法进行定量检测,并与直接利用原始谱线的检测结果进行比较分析。实验结果表明采用二阶导数谱与特征吸收窗的定量检测方法精度更高,相对误差基本在5%以内。     

采用振速测量的改进统计最优近场声全息方法

针对常规统计最优近场声全息在空间多源声场重建过程中所需波函数项数多、重建精度不理想的问题,本文提出了一种基于振速测量的改进统计最优近场声全息算法。与常规算法不同,改进算法主要根据声源特点选取合适的波函数组合来计算声场传递矩阵。通过数值仿真初步验证了该方法的准确性和有效性,并与常规算法进行了详细的对比分析。仿真结果表明,改进算法重建精度高,随频率的变化相对误差波动较小,且随着频率的增大相对误差有逐渐减小的趋势;此外,不同的波函数组合,重建效果差异很大,当选取的波函数与声源共形且数量一致时重建效果最好。     

基于曲率滤波和梯度变换的图像增强

针对基于梯度变换的图像增强算法抗噪声干扰能力差的问题,引入曲率滤波理论,提出了基于高斯曲率滤波和梯度变换的图像增强算法.该算法通过对图像梯度场进行非线性变换来增强图像对比度,通过构造能量泛函,采用梯度下降法从变换后的梯度场重构出增强后的图像,并利用高斯曲率滤波对梯度下降法迭代过程中的重构图像及其各阶偏微分进行平滑,有效解决了图像重构过程中的噪声非线性放大和扩散问题,同时保留了丰富的细节信息.采用多组边缘模糊图像进行仿真实验,实验结果表明该算法在增强图像边缘对比度的同时,能够有效抑制噪声.     

Photon path length analysis of radiative heat transfer in planar layers with arbitrary temperature distributions

The internal source analytical technique is extended to predict the radiative heat transfer for a layer having an arbitrary temperature distribution. By combining a number of internal sources distributed at various optical depths in the layer and weighting them appropriately, a nonisothermal layer is modeled. Heat flux and intensity distributions within layers having a single internal source are presented. The distributions are found to present trends unique to the internal source problem. Isothermal layers are modeled and compare very well with published results. Increased accuracy is attained for all cases and particularly for larger optical depths and smaller albedos by increasing the number of internal sources. The technique is applied to a nonisothermal layer having a temperature distribution similar to that for a hot medium with a cold boundary region. The effect of the boundary region on the normalized heat flux leaving the layer is seen to collapse to a single line for small layer optical thicknesses and large albedos, the slope of which is governed by the temperature ratio T_max/T_min.