求解任意方向辐射强度的广义多流法

本文提出一种基于有限体积法的精确求解任意方向辐射强度的广义多流法(Multi-Flux Method,MFM).以圆柱形吸收、散射、发射性介质为例,采用MFM模拟其沿任意方向的出射辐射强度,并与反向蒙特卡洛法(Backward MonteCarlo,BMC)的计算结果进行了比较分析.结果表明,MFM法与BMC法的计算结果吻合较好,但MFM法的计算效率明显优于BMC法.因此,MFM模型是一种适用于计算任意方向辐射强度问题的高效数值模型.     

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法处理瞬态入射辐射问题的能力．     

任意入射辐射条件下介质特性对瞬态辐射传递的影响

用基于Monte Carlo法的DRESOR法在平行平板系统内具有吸收、无发射介质中研究不同波形入射、壁面反射、介质散射率、光学厚度、各向异性散射等条件对瞬态辐射传递的影响.任意连续波形入射辐射是目前大多数数值方法很难处理的瞬态辐射问题,而DRESOR法通过在系统内计算一单位入射辐射能对介质的DRESOR数分布,就能计算任意连续波形入射辐射条件下高方向分辨率的瞬态辐射强度结果.DRESOR法和Monte Carlo法计算的结果进行了比较验证,两者吻合较好,证明了DRESOR法处理瞬态入射辐射问题的正确性和有效性.     

具有漫反射边界一维灰性平行平板介质中辐射传递方程的求解

本文提出一种基于蒙特卡洛法(Monte Carlo Method, MCM)的新方法在具有漫反射边界一维灰性平行平板介质中求解辐射传递方程(Radiative Transfer Equation,RTE)。该方法能以较高的方向分辨率精确地计算任意点的辐射强度而不需要辐射能量平衡。验证结果显示了该方法的有效性和正确性。并且这种方法具有的一些特点,例如不同辐射源对辐射强度分布贡献的可加性和自动满足边界条件等,使得用该方法更容易处理复杂的辐射换热并行计算问题和复杂的边界条件问题。     

三维燃烧介质和壁面温度的非接触联合重建研究

温度分布在线实时测量对于燃烧过程优化和污染物控制具有重要意义, 针对以往非接触三维温度分布重建过程的耗时性问题和忽略壁面辐射的不足, 本文提出了一种新的离散重建模型, 用于三维吸收、 发射和散射性高温燃烧介质以及壁面温度的快速联合非接触测量. 该模型以四个CCD(Charge Coupled Device) 为测量传感器, 通过构建辐射逆问题求解方程, 从CCD输出的辐射投影图像重建温度分布. 介质中不同投影方向内的辐射传递过程通过离散传递法来描述, 介质的散射和壁面反射则通过离散坐标法来近似. 离散后计算局部辐射强度的病态方程通过最小二乘余量法来求解, 论文对其计算速度进行了优化. 通过非对称温度分布测量算例分析了该模型的有效性, 讨论了测量噪音、 介质和壁面辐射特性对重建精度的影响, 并与其他方法对比分析了模型的重建速度. 计算结果表明本文提出的离散模型可以有效地用于大型高温燃烧介质和壁面温度分布的联合非接触测量. 即使在有噪声的情况下, 该模型也能获得准确的测量结果, 与其他计算方法相比, 采用改进的最小二乘余量法, 能有效地提高温度分布的重建计算速度.     

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

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

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

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

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

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

测量层流扩散火焰温度与烟黑浓度分布的实验研究

本文提出了一种通过线性规划中的内点法(仿射变换法)来计算烟黑浓度和温度分布的模型.根据烟黑辐射特性,利用火焰单色辐射强度图像信息采用此模型同时重建轴对称含烟黑火焰的温度与烟黑浓度分布,对层流乙烯扩散火焰的温度与烟黑浓度进行测量,得到了较好的结果.     

求解辐射传递的非结构混合有限体积／有限元法

本文给了一种适用于任意非结构网格的有限体积/有限元法的混合算法用于求解多维半透明吸收、发射、散射性灰矩形介质内的辐射传递.该方法使用有限元法进行角度离散,有限体积法进行空间离散.与基于辐射传递离散坐标方程的方法不同的是,该方法在迭代求解的过程中,针对每一个空间体元,所有角度方向的辐射强度同时耦合求出.通过两个算例验证了该解法的正确性.     

Finite element approach for radiative transfer in multi-layer graded index cylindrical medium with Fresnel surfaces

Because the optical plane defined by the incidence and reflection direction at a cylindrical surface has a complicated relation with the local azimuthal angle and zenith angle in the traditional cylindrical coordinate system, it is difficult to deal with the specular reflective boundary condition in the solution of the traditional radiative transfer equation for cylindrical system. In this paper, a new radiative transfer equation for graded index medium in cylindrical system (RTEGCN) is derived based on a newly defined cylindrical coordinate system. In this new cylindrical coordinate system, the optical plane defined by the incidence and reflection direction is just the isometric plane of the local azimuthal angle, which facilitates the RTEGCN in dealing with cylindrical specular reflective boundaries. A least squares finite element method (LSFEM) is developed for solving radiative transfer in single and multi-layer cylindrical medium based on the discrete ordinates form of the RTEGCN. For multi-layer cylindrical medium, a radial basis function interpolation method is proposed to couple the radiative intensity at the interface between two adjacent layers. Various radiative transfer problems in both single and multi-layer cylindrical medium are tested. The results show that the present finite element approach has good accuracy to predict the radiative heat transfer in multi-layer cylindrical medium with Fresnel surfaces.     

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 radiative transfer in cylindrical medium and partition allocation method by overlap regions

In this paper the radiative transfer in one-dimensional and two-dimensional cylindrical medium is simulated by the Monte-Carlo (M-C) method. Our results agree with the previous ones very well. It indicates that our cylindrical M-C model is creditable. In this paper, the partition allocation method of radiative heat transfer in participating cylindrical media is presented, in which every sub-domain is isolated by imaginary black wall at certain equivalent temperature and overlaps each other. The stitched results of all sub-domains can predict the results of a whole zone accurately. The partitioned treatment by overlap regions can achieve reasonable result, save memory efficiently, and compute parallel.     

辐射换热的分区计算研究

本文利用蒙特卡洛法模拟圆柱形参与性介质内的辐射换热过程,分区计算了温度仅沿轴向(一维)变化和同时沿轴向和径向(二维)变化的吸收、发射、各向异性散射介质的边界热流分布,将各子区域的边界面假想为等效黑体面,采用由子区间的投射热流计算假想等效黑体面温度的方法,通过分别对介质温度分布一维和二维两种变化情况的研究表明,分区方法可以大大节省计算机内存,同时子区间介质的温度水平,几何尺寸比例以及各区间重合的光学厚度是影响热流分布精度的主要因素。     

The Iterative-DRESOR method to solve radiative transfer in a plane-parallel, anisotropic scattering medium with specular-diffuse boundaries

Using the intensity with high directional resolution obtained by the Basic-DRESOR method as an initial guess, which is substituted into the integrated radiative transfer equation (IRTE), an iterative algorithm is proposed, called the Iterative-DRESOR method. This method can reduce the error levels of the intensity from several percent using the Basic-DRESOR method to a level of less than 1.0×10⁻⁶ with acceptable computation costs. The method is also validated against the exact heat flux in literature in some cases. It further clarifies some uncertain results for the reflectance in a pure, linearly anisotropic scattering medium with specular-diffuse boundaries. The directional distributions of intensity are obviously influenced by the reflecting modes of the boundary, especially in the zone near the boundary. The reflecting mode of an emitting boundary has little effect on the transmittance or reflectance. The reflecting mode of a non-emitting boundary also has little effect on the transmittance, but it obviously influences the reflectance. The difference between the reflectance for specular and diffuse boundaries increases at first, and then decreases, as the optical thickness of the medium increases. The difference will decrease as the scattering albedo of the medium increases, and it is negligible when the medium is pure scattering. The effect of the scattering phase function of the medium on the difference can also not be ignored. The Iterative-DRESOR method is expected to strengthen the capability of the Monte Carlo method to produce accurate results and to validate the results of other methods to solve RTE.     

POSITIVE DEFINITE PROBLEM OF ENERGY DENSITY AND RADIATIVE ENERGY FLUX FOR PULSE CYLINDRICAL GRAVITATIONAL WAVE

By using the general expressions of energy momentum pseudo-tensor of the cylin-drical gravitational waves (GW) given by Rosen and Virbhadra in Cartesian coor-dinates, the concrete forms of energy density and radiative energy flux of the pulse cylindrical GW are obtained. Their physical properties, suitable range and asymp-totic behaviour are considered. It is found that: For the region in which space radial coordinates to origin are greater than the pulse width of the pulse cylindrical GW, the energy density and radiative energy flux of the outward travelling pulse cylindrical GW propagating along at light-cone are positive definite. However, for the region in which the space radial coordinates are less than the pulse width, there is no guarantee for the positive definite property of the radiative energy flux of the outward travelling wave. Moreover, we show that the asymptotic behaviour of the energy and energy flux densities of the pulse cylindrical GW and that of the Riemann curvature tensor have good self-consistancy in space like, time like and null infinity regions.