辐射换热的输运模拟研究

本文从辐时与传输的基本原理出发，引进辐射推动力及流热参数的概念，建立了一般情况下辐射传热的二阶微分方程的数学模型。该方程具有传输方程的特点，使得辐射传热的处理在方法上与导热、流动问题相一致。     

再生核粒子法在辐射传输无网格计算中的应用

基于光滑粒子动力学(SPH)方法发展起来的再生核粒子法(RKPM)是一种拉格朗日型无网格方法。该方法在处理复杂形体内流动、导热、相变、应力应变等多物理场耦合问题时获得了较广泛的应用,但在辐射传输中的适用性尚无研究。而高温半透明材料的相变过程,除以上多物理过程外还需要考虑辐射传输过程的影响。本文基于建立高温半透明颗粒相变过程光热力耦合过程统一的无网格模拟方法的目的,重点研究了RKPM模拟辐射传输的可行性问题。文中首先建立了辐射强度的RKPM拟合公式,然后建立了辐射传输方程的直接配点离散格式。计算了一维参与性灰介质的无量纲热流和无因次温度分布。通过和文献值对比,初步验证了RKPM能够模拟辐射传输方程,且具有较好的计算精度。     

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

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

地表双向反射对天基矢量辐射探测的影响分析

矢量辐射传输方程定量描述了辐射在地表-大气耦合介质中的传输过程,是定量遥感的基础.在处理辐射和离散介质相互作用时,如何处理多次散射、辐射偏振效应和耦合地表模型是研究的重点,直接影响定量化遥感反演的精度.文中基于逐次散射近似法求解了矢量辐射传输方程,求解过程中耦合典型地表的非偏双向反射（BRDF）模型和偏振双面反射（BPDF）模型.采用相对误差因子定量分析了地表双向反射效应和大气偏振效应对天基矢量辐射的影响.为进一步研究地表-大气耦合介质系统的偏振特性及地表大气参数的定量遥感反演提供理论支持. 关键词： 矢量辐射传输方程 逐次散射近似 定量遥感 偏振遥感     

高超声速飞行器气体辐射噪声计算方法

高超声速飞行器周围的激波层内高温气体会发生剧烈的物理化学变化,伴随强烈的光辐射过程,直接影响红外导引头的光学成像效果。采用流体力学Navier-Stokes方程和热化学非平衡模型模拟高温非平衡流动,考虑电子跃迁和振转跃迁以窄带法求解气体辐射特性参数,基于有限体积法离散辐射传输方程,在分别验证流场解算、辐射参数求解和辐射传输计算的基础上,进行了飞行器高速飞行的流动和辐射模拟。数值求解得到了飞行器流场特征和粒子数空间分布。计算的选定波长范围内的气体辐射发射系数空间分布显示其与激波形状和波后气体温度分布相似。通过传输得到的飞行器光学窗口视线路径上的气体辐射噪声成轴对称分布,发现辐射噪声和飞行速度、气体成分等密切相关,马赫数增加时气体辐射噪声显著增强。     

高超声速飞行器气体辐射噪声计算方法

高超声速飞行器周围的激波层内高温气体会发生剧烈的物理化学变化,伴随强烈的光辐射过程,直接影响红外导引头的光学成像效果。采用流体力学Navier-Stokes方程和热化学非平衡模型模拟高温非平衡流动,考虑电子跃迁和振转跃迁以窄带法求解气体辐射特性参数,基于有限体积法离散辐射传输方程,在分别验证流场解算、辐射参数求解和辐射传输计算的基础上,进行了飞行器高速飞行的流动和辐射模拟。数值求解得到了飞行器流场特征和粒子数空间分布。计算的选定波长范围内的气体辐射发射系数空间分布显示其与激波形状和波后气体温度分布相似。通过传输得到的飞行器光学窗口视线路径上的气体辐射噪声成轴对称分布,发现辐射噪声和飞行速度、气体成分等密切相关,马赫数增加时气体辐射噪声显著增强。     

任意阶球谐函数法求解辐射传输问题的误差分析

针对各向异性散射介质内的红外辐射传输开展数值方法研究.应用球谐函数展开推导了一维吸收、发射、散射性灰介质的辐射传输近似方程组,采用差分方法结合三对角矩阵解法建立了任意阶辐射传输近似方程数值解法.并通过计算算例验证了高阶PN方法的数值无关性问题.计算结果表明,高阶球谐函数数值解与理论解吻合得很好,同时计算验证了不同展开阶...     

随机介质中的矢量辐射传输理论

本文讨论了极化电磁波在随机介质中多次散射，传输和热辐射的斯托克斯矢量的辐射传输理论。其中包括随机分布离散的球形和非球形粒子的矢量辐射传输方程，离散坐标——特征分析法，付利叶变换，迭代法等数值解。讨论了非球形粒子的穆勒矩阵。并研究了密集分布的散射粒子介质的辐射传输理论，考虑了密集粒子散射的相干性，计算了有效传播常数。理论及数值结果与实验用了很好的比较。     

二维稳态辐射传输方程的有限差分求解法

针对扩散光学层析在小动物成像中的应用问题并基于混浊介质空间光子三维散射的实际物理效应,提出的二维稳态辐射传输方程的有限差分数值求解新方法.在此基础上,研究了不同的空间剖分网格和角度离散密度对求解准确度的影响,并通过将所提方法与蒙特卡洛模拟进行比对,验证方法的正确性.研究表明:在均匀组织体内,当离散角度达到一定数量时,由辐射传输方程的有限差分解获得的透射面和侧面的光子密度对空间网格大小并不敏感,而在反射面上光子密度计算则需要较密的空间网格才能够达到一定准确度.本研究为发展基于辐射传输方程的扩散光学层析理论奠定了基础.     

二维稳态辐射传输方程的有限差分求解法

针对扩散光学层析在小动物成像中的应用问题并基于混浊介质空间光子三维散射的实际物理效应,提出的二维稳态辐射传输方程的有限差分数值求解新方法.在此基础上,研究了不同的空间剖分网格和角度离散密度对求解准确度的影响,并通过将所提方法与蒙特卡洛模拟进行比对,验证方法的正确性.研究表明：在均匀组织体内,当离散角度达到一定数量时,由辐射传输方程的有限差分解获得的透射面和侧面的光子密度对空间网格大小并不敏感,而在反射面上光子密度计算则需要较密的空间网格才能够达到一定准确度.本研究为发展基于辐射传输方程的扩散光学层析理论奠定了基础.     

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.     

Numerical prediction of heat transfer by natural convection and radiation in an enclosure filled with an isotropic scattering medium

This paper deals with the numerical solution for natural convection and volumetric radiation in an isotropic scattering medium within a heated square cavity using a hybrid thermal lattice Boltzmann method (HTLBM). The multiple relaxation time lattice Boltzmann method (MRT-LBM) has been coupled to the finite difference method (FDM) to solve momentum and energy equations, while the discrete ordinates method (DOM) has been adopted to solve the radiative transfer equation (RTE) using the S8 quadrature. Based on these approaches, the effects of various influencing parameters such as the Rayleigh number (Ra), the wall emissivity (ε_ι), the Planck number (Pl), and the scattering albedo (ω), have been considered. The results presented in terms of isotherms, streamlines and averaged Nusselt number, show that in absence of radiation, the temperature and the flow fields are centro-symmetrics and the cavity core is thermally stratified. However, radiation causes an overall increase in the temperature and velocity gradients along both thermally active walls. The maximum heat transfer rate is obtained when the surfaces of the enclosure walls are regarded as blackbodies. It is also seen that the scattering medium can generate a multicellular flow.     

Solidification of a semitransparent planar layer subjected to radiative and convective cooling

This work deals with the analysis of solidification of a semi-transparent planar layer subjected to radiative and convective cooling. Isothermal mushy-zone model is considered. Enthalpy formulation of the energy equation is solved using the lattice Boltzmann method. To compare the results, the same equation is also solved using the fully implicit finite volume method. Discrete ordinate method is used to compute the radiative information in both the approaches. Effects of radiative properties such as the extinction coefficient, the scattering albedo and refractive index on the solid fraction and temperature are analyzed. Results are validated with those available in the literature. Lattice Boltzmann method and the discrete ordinate method were found to work satisfactorily.     

An inverse analysis of a transient 2-D conduction-radiation problem using the lattice Boltzmann method and the finite volume method coupled with the genetic algorithm

This article deals with the simultaneous estimation of parameters in a 2-D transient conduction-radiation heat transfer problem. The homogeneous medium is assumed to be absorbing, emitting and scattering. The boundaries of the enclosure are diffuse gray. Three parameters, viz. the scattering albedo, the conduction-radiation parameter and the boundary emissivity, are simultaneously estimated by the inverse method involving the lattice Boltzmann method (LBM) and the finite volume method (FVM) in conjunction with the genetic algorithm (GA). In the direct method, the FVM is used for computing the radiative information while the LBM is used to solve the energy equation. The temperature field obtained in the direct method is used in the inverse method for simultaneous estimation of unknown parameters using the LBM-FVM and the GA. The LBM-FVM-GA combination has been found to accurately predict the unknown parameters.     

Monte Carlo ray-tracing simulation for radiative heat transfer in turbulent fluctuating media under the optically thin fluctuation approximation

The Monte Carlo ray-tracing method (MCRT) based on the concept of radiation distribution factor is extended to solve radiative heat transfer problem in turbulent fluctuating media under the optically thin fluctuation approximation. A one-dimensional non-scattering turbulent fluctuating media is considered, in which the mean temperature and absorption coefficient distribution are assumed and the shape of probability density function is given. The distribution of the time-averaged volume radiation heat source is solved by MCRT and direct integration method. It is shown that the results of MCRT based on the concept of radiation distribution factor agree with these of integration solution very well, but results of MCRT based on the concept of radiative transfer coefficient do not agree with these of integration solution. The solution of time-averaged radiative transfer equation by the concept of radiative transfer coefficient should be treated with caution.     

格子Boltzmann方法模拟融化相变过程

建立融化和自然对流耦合的格子Boltzmann双分布函数模型,采用焓法迭代求解相变非线性源项,并通过变松弛时间方法处理固液两相变热物性间题.热传导融化过程的计算值与分析解的对比分析说明该模型能够准确地模拟融化过程.自然对流条件下融化过程的模拟结果表明自然对流在一定程度上影响了融化传热及融化速率等,体现出与热传导融化不同...     

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.     

半透明梯度折射率介质内辐射熵传递方程及其数值模拟

在非相干辐射条件下,基于Planck光谱辐射熵强度定义,导出半透明梯度折射率介质内光谱辐射熵传递方程,以及局部辐射熵产率理论表达式.基于离散坐标法对辐射熵传递方程进行数值求解.以一维半透明梯度介质平板为例,对辐射熵方程及其算法进行验证.平板整体无因次辐射熵产的计算结果与宏观热力学定律的结果一致.     

A Lattice Boltzmann Model for Fluid-Solid Coupling Heat Transfer in Fractal Porous Media

We report a lattice Boltzmann model that can be used to simulate fluid-solid coupling heat transfer in fractal porous media. A numerical simulation is conducted to investigate the temperature evolution under different ratios of thermal conductivity of solid matrix of porous media to that of fluid. The accordance of our simulation results with the solutions from the conventional CFD method indicates the feasibility and the reliability for the developed lattice Boltzmann model to reveal the phenomena and rules of fluid-solid coupling heat transfer in complex porous structures.     

Laser short-pulse heating of silicon-aluminum thin films

Energy transport in silicon-aluminum thin films is examined during the laser short-pulse irradiation subjected to the silicon film. The silicon film is considered to be at the top of the aluminum film. Thermal boundary resistance at the interface of the films is incorporated in the analysis. The absorption of laser radiation in the silicon and aluminum films is modeled using the transfer matrix method. Since the silicon film is dielectric, the phonon radiative transport basing the Boltzmann transport equation is incorporated to determine equivalent equilibrium temperature in the film while modified two-equation model is used to account for the non-equilibrium energy transport due to thermal separation of electron and phonon sub-systems in the aluminum film during the laser short-pulse heating process.