DRESOR法对瞬态辐射传递问题的研究

用基于蒙特卡洛法(Monte Carlo Method,MCM)的DRESOR法(Distributions of Ratios of Energy Scattered by the medium Or Reflected by the boundary surface)求解入射辐射经过介质散射、壁面反射传递后辐射强度随时间变化的瞬态辐射传递方程(Transient RadiativeTransfer Equation,TRTE)问题。通过在系统内计算一单位瞬态入射辐射对介质的DRESOR数分布,就能计算任意时间内入射辐射在系统内时间响应特性,这样有效提高数值方法处理瞬态辐射问题的通用性。并且能够获得高方向分辨率的辐射强度随时间变化的结果,这是目前大多数数值处理方法比较难做到的,显示出了DRESOR法处理瞬态入射辐射问题的能力．     

DRESOR法对二维辐射传递问题研究

高方向辐射强度信息对一些逆问题非常有帮助,例如,在工业炉中用辐射成像技术重建二维/三维温度场.DRESOR法被发展用来在二维矩形各向同性/异性散射介质中求解辐射传递方程.用DRESOR法可以在矩形区域边界处可以计算得到辐射强度在半球立体角空间内6658个方向上的分布.用DRESOR法计算得到的无维辐射热流和有关文献用中离散坐标法、近似法及有限体积法计算得到的结果吻合得很好.从计算结果中可以观察边界上,特别是靠近发射源区域的辐射强度随方位角的变化情况.     

DRESOR法对平行入射辐射问题的研究

本采用一种基于蒙特卡洛法(Monte Carlo Method,MCM)求解辐射传递方程(Radiative Transfer Equation, RTE)的快捷、有效的方法-DRESOR法(Distributions of Ratios of Energy Scattered Or Reflected)在一维充满吸收、各向同性散射介质平行平板中,外部有平行入射条件下,求解计算空间点的辐射强度沿空间方向角的分布,而不需要辐射平衡和在空间位置坐标和方向角度坐标上同时离散辐射传递方程进行迭代求解。     

DRESOR法求解超短脉冲在一维非均匀介质中的瞬态辐射传输

本文中DRESOR法(Distributions of Ratios of Energy Scattered Or Reflected)应用于求解一维非均匀介质内的瞬态辐射传输问题.结果表明,反射热流的双峰分布现象能够体现介质内部的非均匀分布特征,并且能够通过分析其特点推断非均匀介质的中间界面位置.通过DRSESOR法计算得到的高方向分辨率的辐射强度体现了其在决定非均匀介质表面位置的独特优势.最后,对双峰分布现象的存在条件进行了探讨.     

矩形介质内辐射换热的有限元法

利用有限单元法离散求解辐射传递方程和能量控制方程.分别计算了边界为黑体和灰体条件下矩形吸收、发射、各向同性散射介质内的平均入射强度和温度分布,并同蒙特卡罗法(M-C法)计算结果进行了比较.     

圆柱坐标系下任意方向辐射强度的源项六流法模拟

基于传统热流法,提出一种圆柱坐标系下的源项六流模型(Source Six Flux,SSF),可快速准确地计算参与性介质内任意方向的出射辐射强度.详细介绍SSF模型的基本原理和求解步骤,以圆柱形吸收、散射、发射性介质为例,模拟其沿任意方向的出射辐射强度,并与反向蒙特卡罗法(Backward Monte Carlo,BMC)和二流法(Two Flux Method,TFM)的计算结果进行比较.结果表明,SSF法与BMC法的计算结果吻合较好,计算精度均高于TFM法,但SSF法的计算效率明显优于BMC法.因此,SSF模型是一种适用于计算任意方向辐射强度问题的高效数值模型.     

一维瞬态辐射传输的有限元法模拟

随着超短脉冲激光的快速发展,吸收散射性介质内的瞬态辐射传输引起了人们的广泛关注.本文基于离散坐标法和最小二乘有限元法(LSFEM),提出了模拟多维吸收散射性介质内瞬态辐射传输的数值模型.该模型有效地克服了在标准Galerkin有限元法(GFEM)中发生的伪振荡现象,在时间步长较大的情况下仍然可以得到光滑无振荡的解.而且,最小二乘法产生的求解系数矩阵是对称正定的,与GFEM中的系数矩阵相比,仅需要存储一半的非零系数,可以应用许多高效的迭代求解方法进行求解.为了检验模型,本文研究了一维吸收散射性介质内瞬态辐射传输问题,其结果与蒙特卡洛法(MCM)和积分模型法(IE)的结果进行了比较,结果证实:本文的方法可以精确、高效地模拟参与性介质内的瞬态辐射传输.     

散射性介质瞬态矢量辐射传输的蒙特卡洛模拟

瞬态效应和偏振特性对于短脉冲激光在散射性介质中的传输有重要影响。本文采用蒙特卡洛法来求解一维散射性介质内的瞬态偏振辐射传输问题。采用拒绝法确定光束的散射方向。定义了瞬态矢量辐射传输矩阵(TVRTM)来描述瞬态偏振辐射的传输特性,并以此得到了Stokes矢量的角度与时间分布。在蒙特卡洛模型中采用时间平移和叠加原理,可大幅度提高计算精度和计算效率。     

超短脉冲激光在二维非均匀介质内的瞬态辐射传输

随着超短脉冲激光的快速发展,吸收散射性介质内的瞬态辐射传输引起了广泛的关注.本文采用最小二乘有限元法模拟了超短脉冲激光局部入射条件下,具有高散射核的二维非均匀介质内的瞬态辐射传输.研究了不同边界位置上反射和透射信号随时间的变化情况.结果表明,对于具有高散射核的非均匀介质,能够揭示散射核位置的双峰现象可能在早期的反射信号中发生,早期的反射信号比后期的、经过介质衰减后信号更重要.因此,在利用短脉冲激光进行光谱分析和成像等技术中,人们应该重视早期反射信号的测量.     

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

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

The influence of anisotropic scattering on the radiative intensity in a gray, plane-parallel medium calculated by the DRESOR method

Even though there have been many ways to treat complex anisotropic scattering problems, in most of the cases only the radiation flux or its dimensionless data were provided, and radiative intensity with high directional resolution could merely be seen. In this paper, a comprehensive formulation for the DRESOR method was proposed to deal with the anisotropic scattering, emitting, absorbing, plane-parallel media with different boundary conditions. The method was validated by the data from literature and the integral formulation of RTE. The DRESOR value plays an important role in the DRESOR method, and how it is determined by the anisotropic scattering was demonstrated by some typical results. The intensities with high directional resolution at any point can be given by the present method. It was found that the scattering phase function has little effect on the intensity for thin optical thickness, for example, 0.1. And there is the largest boundary intensity for the medium with the largest forward scattering capability, and the smallest one with the largest backward scattering capability. An attractive phenomenon was observed that the scattering of the medium makes the intensity at boundary can not reach the blackbody emission capability with the same temperature, even if the optical thickness tends to very large. It was also revealed that the scattering of the medium does not mean it cannot alter the magnitude of the energy; actually, stronger scattering causes the energy to have more chance to be absorbed by the medium, and indirectly changes the energy magnitude in the medium. Finally, it is easy to deduce all the associated quantities such as the radiation flux, the incident radiation and the heat source from the intensity, just as done in literature.     

The DRESOR method for radiative heat transfer in a one-dimensional medium with variable refractive index

This paper extends the DRESOR (Distribution of Ratios of Energy Scattered by the medium Or Reflected by the boundary surface) method to radiative transfer in a variable refractive index medium. In this method, the intensity is obtained from the source term along the curved integration paths determined only by the variable refractive index, and the DRESOR values are calculated by the Monte Carlo method in which the propagation of the energy bundles are affected by Snell's law. With given temperatures on the black boundaries of a one-dimensional medium, the temperature distribution inside the medium with a variable scattering property is calculated under the condition of radiative equilibrium. It is shown that the DRESOR method has a good accuracy in the cases studied. For an isotropic-scattering medium with the same optical thickness, the scattering albedo has no effect on the temperature distribution, which can be obtained from the general equations and can be seen as an extension of what exists for a constant refractive index; however, the different refractive index causes obvious changes in the temperatures inside the medium. The effect of anisotropic scattering on the temperature distribution cannot be ignored, although it is still weaker than the effect caused by variation in the refractive index.     

Simulation of the infrared radiation characteristics of high-temperature exhaust plume including particles using the backward Monte Carlo method

The analysis of the infrared radiation characteristics of high-temperature free-stream flow including particles is very significant for the field of target detection, combustion diagnosis and temperature measurement of flame. In this paper, the infrared radiation characteristics of high-temperature free-stream flow are calculated and analyzed using the backward Monte Carlo method, considering the effect of the directional radiation heat flux due to the particle scattering and the different boundary conditions. The calculation results of emitting, absorbing and anisotropically scattering media are compared with the forward Monte Carlo and finite-volume methods results, which shows the superiority on computational efficiency with the backward Monte Carlo method.     

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

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

Discrete ordinate and Monte Carlo simulations for polarized radiative transfer in a coupled system consisting of two media with different refractive indices

We present an algorithm for polarized radiative transfer in a vertically stratified system consisting of two plane-parallel media with different refractive indices. It is based on the discrete ordinate method and includes multiple elastic scattering, thermal radiation, Fresnel reflection and transmission, incident parallel-beam or isotropic radiation at the top of the upper medium and bidirectional reflection at the bottom of the lower medium. Comparisons with results from Monte Carlo simulations show that the discrete-ordinate code provides accurate results for all four elements of the Stokes vector (I, Q, U, and V) at a speed that is orders of magnitude faster.     

Backward Monte Carlo simulation for apparent directional emissivity of non-isothermal semitransparent slab

A backward Monte Carlo method based on radiation distribution factor is extended to compute the apparent directional emissivity of one-dimensional absorbing-emitting-scattering semitransparent slab with specular semitransparent surface and opaque diffuse substrate. A non-isothermal semitransparent slab is taken as an example to examine the efficiency of backward Monte Carlo method. The apparent directional emissivity of the semitransparent slab is determined by backward Monte Carlo simulation and compared with that determined by forward Monte Carlo simulation. The results show that the backward Monte Carlo method is very efficient in solving the apparent directional emissivity of the semitransparent slab. The backward Monte Carlo method converges more quickly than the forward Monte Carlo method.