管内高温介质层流入口段中的热辐射作用

数值研究了高温介质密度随温度变化时,管内层流入口段耦合换热中的热辐射作用。采用离散坐标法、控制容积法耦合求解辐射传递方程、能量方程及N-S方程。考察了中等大小光学厚度下,热辐射作用对介质内速度分布、温度分布以及换热的影响。结果表明,即使在不大的光学厚度下,热辐射作用对管内高温介质层流入口段耦合换热的速度场与换热强度都有明显影响。     

Non-Fourier heat conduction with radiation in an absorbing, emitting, and isotropically scattering medium

This article numerically analyses the combined conductive and radiative heat transfer in an absorbing, emitting, and isotropically scattering medium. The non-Fourier heat conduction equation, which includes the time lag between heat flux and the temperature gradient, is used to model the conductive heat transfer in the medium. It predicts that a temperature disturbance will propagate as a wave at finite speed. The radiative heat transfer is solved using the P₃ approximation method. In addition, the MacCormack's explicit predictor-corrector scheme is used to solve the non-Fourier problem. The effects of radiation including single scattering albedo, conduction-to-radiation parameter, and optical thickness of the medium on the transient and steady state temperature distributions are investigated in detail. Analysis results indicate that the internal radiation in the medium significantly influences the wave nature. The thermal wave nature in the combined non-Fourier heat conduction with radiation is more obvious for large values of conduction-to-radiation parameter, small values of optical thickness and higher scattering medium. The results from non-Fourier-effect equation are also compared to those obtained from the Fourier equation. Non-Fourier effect becomes insignificant as either time increases or the effect of radiation increases.     

A Monte Carlo implementation to solve radiation heat transfer in non-uniform media with spectrally dependent properties

This paper presents the application of the Monte Carlo method to solve the radiative heat exchange in non-homogeneous, non-isothermal gases with spectrally dependent properties. Among others models, the absorption-line blackbody (ALB) distribution function, originally defined and derived for the spectral line-based weighted-sum-of-gray-gases (SLW) model, allows an immediate, simple implementation of the Monte Carlo method to account the spectral dependence of the radiative properties. This work shows how the Monte Carlo method can be combined to the ALB distribution function, and provides results for heat transfer in a mixture of water vapor, carbon dioxide and nitrogen that have satisfactory agreement with the SLW method and with line-by-line integration. Finally, the solution technique is employed to solve two examples aiming at demonstrating the effect of the absorbing species concentration on the thermal radiative exchanges. The method is of great interest for the computation of radiative transfer in combustion systems where the chemical species concentration and the temperature are not uniform.     

Radiation heat transfer in axisymmetric quartz glass tubes

A mathematical model is proposed to describe the heat transfer in quartz glass axisymmetric tubes. Heat is transferred inside the glass by radiation and conduction. Scattering of thermal radiation inside the glass is ignored. At the boundaries of the tube the radiative intensity is specularly reflected. The spectral dependent radiative intensity and the temperature distribution inside the tube are determined. The model is applied to simulate the cooling process of a quartz glass tube. The calculated temperature is in agreement with that obtained from an experiment. Furthermore, steady-state temperature distributions in quartz glass tubes of different lengths have been determined.     

Temperature measurements of high-temperature semi-transparent infrared material using multi-wavelength pyrometry

Temperature measurements inside semi-transparent materials are important in many fields. This study investigates the measurements of interior temperature distributions in a one-dimensional semi-transparent material using multi-wavelength pyrometry based on the Levenberg–Marquardt method (LMM). The investigated material is semi-transparent Zinc Sulfide (ZnS), an infrared-transmitting optical material operating at long wavelengths. The radiation properties of the one-dimensional semi-transparent ZnS plate, including the effective spectral–directional radiation intensity and the proportion of emitted radiation, are numerically discussed at different wavelengths (8.0–14.0 μm) and temperature distributions (400–800 K) to provide the basic data for the temperature inversion problem. Multi-wavelength pyrometry was combined with the Levenberg–Marquardt method to resolve the temperature distribution along the radiative transfer direction based on the line-of-sight spectral radiation intensities at multiple wavelengths in the optimized spectral range of (11.0–14.0 μm) for the semi-transparent ZnS plate. The analyses of the non-linear inverse problem show that with less than 5.0% noise, the inversion temperature results using the Levenberg–Marquardt method are satisfactory for linear or Gaussian temperature distributions in actual applications. The analysis provides valuable guidelines for applications using multi-wavelength pyrometry for temperature measurements of semi-transparent materials.     

辐射(火积)耗散与空间辐射器温度场均匀化的关系

(火积)理论是针对传热优化发展起来的,并获得了越来越多的应用。本文基于辐射(?)原理,对空间辐射器散热过程中的散热量分布、发射率分布和散热面积分布问题进行了分析。对于以上三类优化分布问题,理论分析和数值计算均表明,辐射器最小的辐射(?)耗散和辐射热阻均对应于散热表面均匀的温度场。因此,辐射(?)原理可用于空间辐射器的温度场均匀化设计。     

Theoretical study of combined radiative conductive and convective heat transfer in semi-transparent porous medium in a spherical enclosure

A new method for the solution of the radiative transfer equation in spherical media based on a modified discrete ordinates method is extended to study radiative, conductive and convective heat transfer in a semi-transparent scattering porous medium. The set of differential equations is solved using the fourth-order Runge-Kutta method. Various results are obtained for the case of combined radiative and conductive heat transfer, as well as for the interaction of those modes with convection. The effects of some radiative properties of the medium on the heat transfer rate are examined.     

Simulation of radiative heating of snow and ice coating

Numerical simulation of thermal state was performed for semitransparent medium; the thermal state is formed under impact of incident radiation and convective heat transfer to the ambient medium. The model takes into account the heat transfer to the semi-infinite opaque substrate below the layer of semitransparent material. Computation was performed for governing parameters typical of snow and ice coating during winter period. Solving the radiative part of problem uses the modified average flux method. This method takes into account the dependence of optical properties on the wavelength of incident radiation, scattering, and reflective properties of the layer boundaries.     

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.     

Radiative gas dynamics of the Fire-II superorbital space vehicle

The rates of convective and radiative heating of the Fire-II reentry vehicle are calculated, and the results are compared with experimental flight data. The computational model is based on solving a complete set of equations for (i) the radiative gas dynamics of a physically and chemically nonequilibrium viscous heatconducting gas and (ii) radiative transfer in 2D axisymmetric statement. The spectral optical parameters of high-temperature gases are calculated using ab initio quasi-classical and quantum-mechanical methods. The transfer of selective thermal radiation in terms of atomic lines is calculated using the line-by-line method on a specially generated computational grid that is nonuniform in radiation wavelength.     

Determination of true temperature of opaque materials via spectral distribution of thermal radiation intensity: application of Wien’s displacement law

The equation for the derivative connecting surface spectral emissivity, wavelength, and thermodynamic (true) temperature of an opaque heated body at the point of spectral maximum of thermal radiation was obtained. It is suggested to solve the problem of determining the true temperature of an opaque surface in two stages. At the first stage, the spectral range, most comfortable for approximation of body emissivity, is distinguished using a special function (relative emissivity), and the true temperature is determined. At the second stage, the true temperature is determined again using the resulting equation for the derivative. The dimensionless parameter that connects the radiative properties of material with the peak wavelength and characterizes deviation from Wien’s displacement law was found. If the absolute value of this parameter is low, the value of true temperature obtained at the first step can be specified at the second stage. This approach is illustrated by experimental data obtained at comparison of spectral radiance of the temperature lamps.     

气动加热下高温陶瓷材料内的瞬态耦合换热

采用控制容积法结合蒙特卡罗法与谱带模型,数值模拟了气动加热下高温陶瓷内的辐射与导热瞬态耦合换热.分析了冷却边界条件、各向同性散射和表面特性对高温陶瓷内瞬态温度分布的影响.结果表明,高温下陶瓷内的辐射换热作用十分明显,陶瓷材料的表面特性对瞬态温度分布影响很大,各向同性散射也有明显影响.     

DIFFERENTIALLY MOVING MEDIA WITH MANY SPECTRAL LINES: STOCHASTIC APPROACH

Based upon the analytical solution of the radiative transfer equation for a given source function and a new approach to account for very many spectral lines contributing to the extinction, the connection between line properties and the emergent intensity is derived under the assumption that the wavelengths of the line centers follow a Poisson point process, whereas the other line parameters may have arbitrary distribution functions.A comparison with the widely used list of Kurucz shows that the Poisson distribution well describes deterministic “real” lines. The presentation by a Poisson point process requires only a modest number of parameters and is very flexible. It allows most operations to be carried out analytically and hence is very suitable to study the intricate influence of many lines on radiation fields in differentially moving media.We consider a simplified case of the solution of the radiative transfer equation in order to demonstrate the basic effects of the velocity field upon the emerging radiation field. Expressions for the expectation value of the intensity are derived, and examples are given for Lorentz line profiles and infinitely sharp lines, in particular as functions of the velocity gradient and the mean line density.     

Thermally driven sensitive surface phonon polariton modes in thin-film SiC

The results of this work show that significant changes to the spectral properties of thermal emission may be accomplished by taking advantage of surface phonon polariton (SPP) modes that are very sensitive to small refractive index changes. Surface phonon polaritons are interface restricted wave modes that do not usually participate in thermal radiative energy exchange. In this work, the possibility of varying the spectral properties by using a thin-film surface active material, layered between dielectrics, is explored. Unlike surface gratings, the multilayered structure discussed in this work couples the surface modes into radiative modes using subtle changes in refractive index. The proposed multilayer effect is studied by first developing the SPP mode dispersion relation to identify the spectral location of the surface modes. The relationship between SPP coupling and refractive index change is then developed. The fluctuation-dissipation theorem and an impedance boundary condition is used to theoretically show how thermal radiation can excite SPP modes. Results are then presented that show how film thickness and substrate properties act to support SPP modes that are particularly sensitive to changes in refractive index and therefore offers opportunities for dramatic changes to the spectral properties.     

Inversion of combustion gas temperature/concentration profile with radiation/turbulence interaction using SRS

The spectral remote sensing (SRS) method is applied to the combustion gas with radiation/turbulence interaction to invert the temperature and concentration profiles. The flame is made and controlled to be optically thin per each fluctuation length so that spatially fluctuating gas layer can be treated equivalently as timewisely fluctuating one sharing identical expression in the radiative transfer equation. Then, the spectral intensity, temperature and concentration distributions are measured for the inversion and as the reference solutions. From the inversion results, we find that SRS can successfully invert the coupled temperature/concentration fluctuation amplitudes and mean values. For the two cases of experiment, inverted values are in good agreement with measured ones. However, SRS cannot find the detailed local fluctuation parameter such as pattern or phase, etc. as far as they do not affect the resulting radiative intensity. So, it is deduced that the pattern or phase is not so influential, whereas the mean value and the fluctuation amplitude of temperature/concentration profile are important parameters for the radiative intensity in the case of the optically thin condition. Further, the radiation/turbulence interaction is verified to play an important role in the radiation.     

热辐射输运问题的隐式蒙特卡罗方法求解

热辐射与物质相互作用及辐射光子在物质中的传输是惯性约束聚变研究中的重要课题. 介绍了基于隐式蒙特卡罗方法的辐射输运方程,在该方程的积分-微分形式基础上,推导了利于蒙特卡罗方法模拟的等价的积分输运方程;基于积分方程设计数值模拟流程,编写三维蒙特卡罗数值模拟程序;针对热辐射输运典型问题及benchmark问题开展了数值实验,计算结果验证了方法的适应性及程序的正确性. 关键词： 热辐射 惯性约束聚变 输运方程 隐式蒙特卡罗     

金属冲击温度测量中金属/窗口界面表观光谱辐亮度的再研究

 在冲击压缩下，理想接触的金属/窗口界面温度历史是时间无关的，因而一般认为界面的表观光谱辐亮度也是时间无关的。研究表明：（1）在冲击压缩下，将伴随金属对光的吸收系数减小和光学厚度增大。因此，在考虑辐射输运效应后，理想接触的金属/窗口界面处的表现光谱辐亮度是时间相关的。非常明显的时间相关过程自冲击波到达界面后持续约10 ns。（2）若把存在空间温度梯度的金属界面辐射看作是透光厚度内平均温度的等效辐射效应，其结果与辐射输运效应的计算结果是很相近的。     

碳氢扩散火焰实验研究：燃烧气氛对火焰光谱与温度的影响

辐射是各种燃烧过程中热传递的主要方式。在不同的火焰中,辐射光谱分布十分复杂。在这项工作中,利用光谱仪测量了可见光(200～900 nm),近红外(900～1 700 nm)和中红外(2 500～5 000 nm)波段火焰的光谱强度,分析了空气和富氧气氛下扩散火焰的光谱特征。并基于光谱分析,定量得到了火焰中碳烟以及气体发射的辐射力,计算了火焰的温度分布。结果表明,空气燃烧中的火焰温度低于富氧燃烧中的火焰温度。在空气气氛下,火焰中的碳烟和气体均对中的热辐射起着重要作用。而在富氧气氛下,气体对于火焰热辐射更为重要。在可见光和近红外波段,由于在空气气氛下火焰中碳烟的大量形成,光谱曲线显示出了良好连续性。而富氧气氛下火焰的辐射光谱降低。在中红外波段,空气气氛下火焰的气体辐射明显弱于富氧气氛下火焰的气体辐射。     

Numerical study of heat transfer in an optically thick semi-transparent spherical porous medium

In this paper, heat transfer by simultaneous convection, conduction and radiation in a semi-transparent spherical porous medium is investigated. The ROSSELAND approximation is adopted to take account of radiation in the heat transfer rate. The routine used here to solve the set of differential equations is taken from the IMSL MATH/LIBRARY. Various results are obtained for the dimensionless temperature profiles in the solid and fluid phases, the radiative, conductive, convective and total heat fluxes. The effects of some radiative properties of the medium on the heat transfer rate are examined.     

热辐射输运问题的高效蒙特卡罗模拟方法

惯性约束聚变研究中,热辐射光子在介质中的输运以及热辐射光子与介质的相互作用是重要研究课题,蒙特卡罗方法是该类问题的重要研究手段之一.隐式蒙特卡罗方法虽然能正确地模拟热辐射在介质中的输运过程,但当模拟重介质(材料的吸收系数大)问题时,该方法花费的计算时间将变得很长,导致模拟效率很低.本文以离散扩散蒙特卡罗方法为基础,开发了"离散扩散蒙特卡罗方法辐射输运模拟程序",可以较好地解决重介质区的计算效率问题,但是离散扩散蒙卡罗方法在模拟轻介质区时精度不够高.辐射输运问题中通常既有轻介质也有重介质,为了能同时解决蒙特卡罗方法模拟的效率和精度问题,本文研究了离散扩散蒙特卡罗方法与隐式蒙特卡罗方法相结合的模拟方法,并提出了新的扩散区与输运区界面处理方法,研制了混合蒙特卡罗方法的辐射输运模拟程序.典型辐射输运问题模拟显示:在模拟重介质问题时,该程序能大幅缩短模拟时间,且能取得与隐式蒙特卡罗方法一致的结果;在模拟轻重介质均存在的问题时,与隐式蒙特卡罗方法相比,混合蒙特卡罗方法的模拟精度与其相当且计算效率同样能够得到显著提升.