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

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

关于具有定向自旋粒子的契连科夫辐射

量子电动力学方法研究在介质中运动的带电粒子自旋状态对契连科夫辐射强度的影响,指出:在横向自旋状态下,辐射强度具有新的量子补充值。和纵向粒子一样,辐射由两部分组成:极化部分(在辐射阈上cosθ=1等于零)和非极化部分(在阈上区别于零)。     

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

半透明介质中辐射传递方程的反演计算及数值模拟

本文介绍了一侧为半透明、另一侧为非透明界面时一维半透明介质的辐射强度计算式。采用辐射与导热复合换热模型计算半透明介质内温度场。利用已知的温度场求半透明介质的辐射强度一正问题计算。将此辐射强度代入辐射反问题计算模型，引入测量误差，采用Ｃｈａｈｉｎｅ方法及演半透明介质内温度场一反问题计算。数值模拟表明，本文所采用的辐射反演法具有较高的精度及稳定性。     

基于辐射光谱法的高温颗粒辐射传热参数在线测量方法研究

针对在高温燃烧环境中的颗粒辐射传热问题,基于普朗克辐射定律,提出了用于高温颗粒辐射传热参数在线测量的辐射光谱法,根据高温颗粒可见波段辐射光谱随波长变化情况,通过参数拟合方法直接获得颗粒温度及辐射强度等辐射传热参数。为验证该方法测量准确性,搭建了高温黑体炉辐射测量系统,实验测量结果显示：温度测量值与设定温度相对偏差小于3%;辐射强度测量值与理论计算值相对偏差小于5%。以此为基础,设计了应用于高温燃烧环境下的颗粒辐射传热参数测量的水冷结构探针,并利用该探针开展了高温燃烧环境气固两相流200～1 100 nm波段辐射光谱测量,基于上述方法,直接获得了高温颗粒温度、辐射强度等辐射传热参数沿截面分布情况,有效剥离了高温气体对流传热的影响,为高温颗粒辐射传热研究提供数据支撑。     

基于辐射光谱法的高温颗粒辐射传热参数在线测量方法研究（英文）

针对在高温燃烧环境中的颗粒辐射传热问题,基于普朗克辐射定律,提出了用于高温颗粒辐射传热参数在线测量的辐射光谱法,根据高温颗粒可见波段辐射光谱随波长变化情况,通过参数拟合方法直接获得颗粒温度及辐射强度等辐射传热参数。为验证该方法测量准确性,搭建了高温黑体炉辐射测量系统,实验测量结果显示:温度测量值与设定温度相对偏差小于3%;辐射强度测量值与理论计算值相对偏差小于5%。以此为基础,设计了应用于高温燃烧环境下的颗粒辐射传热参数测量的水冷结构探针,并利用该探针开展了高温燃烧环境气固两相流200～1 100nm波段辐射光谱测量,基于上述方法,直接获得了高温颗粒温度、辐射强度等辐射传热参数沿截面分布情况,有效剥离了高温气体对流传热的影响,为高温颗粒辐射传热研究提供数据支撑。     

各向异性介质中契连科夫辐射的量子理论

本文在各向异性介质中电磁场二次量子化的基础上,探讨做高匀速运动(超过该介质中的光速)的带电粒子辐射状况,得到了寻常波和非常波辐射强度的表达式,解释了粒子的自旋状态对辐射强度的影响。在粒子的运动方向沿晶体光轴和垂直于晶体光轴的两种情况下,对结果进行了分析。结果表明,量子方法的计算有别于经典结果,当考虑到粒子自旋时,在沿光轴这动的情况下,既辐射寻常波又辐射非常波。     

复燃对液体火箭尾焰红外辐射特性的影响

为了定量研究复燃对液体火箭尾焰红外辐射特性的影响,建立了一个可以计算液体火箭尾焰复燃流场和红外辐射特性的模型.首先,使用FLUENT软件计算液体火箭尾焰复燃流场,其中尾焰中的复燃反应使用有限速率化学反应模型计算;然后,使用基于HITEMP数据库的窄带模型计算尾焰内气体的辐射参量;最后,使用有限体积法求解尾焰中的辐射传输方程.通过比较该模型计算的Titan ⅢB尾焰光谱辐射强度与(美国)国家航空航天局公布结果的一致性,证明了该模型的正确性.最后,利用该模型计算了复燃对某液体火箭尾焰光谱和波段红外辐射强度的影响,结果表明,复燃反应可以显著增加尾焰红外光谱辐射强度,在2.5～3.0 μm和4.2～4.7 μm两个主要辐射波段平均辐射强度的增加比例分别达到了30.8％和28.3％,所以,在计算液体火箭尾焰准确的红外辐射特性时,需要考虑复燃的影响.     

用MEI方法研究脉冲线源辐射问题

 脉冲线电流源的辐射问题是研究高功率微波辐射、传播和散射问题的基础。用时域有限差分方法结合时域MEI吸收边界条件来对该问题进行求解，通过与不同种类边界条件配合得出的数值解的比较可以证明，利用时域MEI方法的一阶吸收边界条件对线源辐射问题进行处理时，可以使截断边界离源更近，从而减少了计算量，并达到提高计算精度的目的。     

高温热辐射输运模拟的蒙特卡罗全局降方差方法

针对隐式蒙特卡罗方法模拟高温热辐射输运问题时存在的“辐射强度计算误差时间空间分布严重不均匀”现象,通过理论分析和数值模拟手段,找到决定误差大小的主要因素为“网格内的辐射径迹长度记录数”.据此提出“隐式蒙特卡罗全局降方差方法”并推导相应的计算公式.该方法主要包含如下3个关键技术:1)针对辐射输运蒙特卡罗模拟的自适应动态分配源粒子方法; 2)与自适应动态分配源粒子相匹配的动态权窗设计技术及粒子权无偏估计算法; 3)辐射强度的解析估计降方差方法.针对这3个关键技术,设计蒙特卡罗数值模拟方案,编写相应的数值模拟程序模块.典型辐射输运问题模拟结果显示:隐式蒙特卡罗全局降方差方法能够使网格辐射强度计算误差在整个时空范围内分布相对比较均匀,最大误差可控,计算效率提升10倍左右.新方法在激光惯性约束聚变的黑腔辐射输运模拟应用中取得了显著效果.     

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.     

An efficient diffusion approximation for 3D radiative transfer parameterization: application to cloudy atmospheres

The three-dimensional (3D) diffusion radiative transfer equation, which utilizes a four-term spherical harmonics expansion for the scattering phase function and intensity, has been efficiently solved by using the full multigrid numerical method. This approach can simulate the transfer of solar and thermal infrared radiation in inhomogeneous cloudy conditions with different boundary conditions and sharp boundary discontinuity. The correlated k-distribution method is used in this model for incorporation of the gaseous absorption in multiple-scattering atmospheres for the calculation of broadband fluxes and heating rates in the solar and infrared spectra. Comparison of the results computed from this approach with those computed from plane-parallel and 3D Monte Carlo models shows excellent agreement. This 3D radiative transfer approach is well suited for radiation parameterization involving 3D and inhomogeneous clouds in climate models.     

Two analytic methods applied to radiative transfer in the linear-anisotropic scattering medium with graded index and diffuse gray boundaries

The curved ray-tracing method is extended to radiative transfer in the graded index medium with diffuse gray boundary conditions instead of black boundary conditions and the pseudo-source adding method is extended to the case of the linear-anisotropic scattering medium with graded index from non-scattering medium. Furthermore, the equivalence of the two methods is verified by formulation derivation. As exact analytical solutions, both the methods have high accuracy and fast computational speed. The predicted temperature distributions and dimensionless radiative heat flux at radiative equilibrium are determined by the proposed methods, and the numerical results are compared with the data in references. The results show that the present methods have a good accuracy. Influences of various combinations of refractive index and boundary emissivities on the temperature distributions and dimensionless radiative heat flux are also investigated.     

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.     

Radiative intensity solution and thermal emission analysis of a semitransparent medium layer with a sinusoidal refractive index

The radiative intensity in a sinusoidal refractive index semitransparent medium layer is solved by the curved ray-tracing method in combination with the pseudo-source adding method. One boundary of the medium layer is an opaque diffuse substrate wall. The other boundary is a semitransparent specular or diffuse surface, from which the medium thermal emission emerges. With considering a linear temperature distribution, the radiative intensity formulae are, respectively, deduced under the two boundary conditions. On the basis of the radiative intensity solutions, the directional and hemispherical emission of the medium layer with a specular surface as well as the hemispherical emission of that with a diffuse surface are calculated. The influences of the optical thickness, sinusoidal refractive index distribution and linear temperature distribution on the thermal emission are investigated. The results show that the effects of refractive index and temperature distribution are significant and are different under the two reflecting modes of the surface.     

An inverse radiative transfer problem of simultaneously estimating profiles of temperature and radiative parameters from boundary intensity and temperature measurements

A parallel-plane space filled with absorbing, emitting, isotropically scattering, gray medium is studied in this paper. The boundary intensity and boundary temperature profiles are calculated for the inverse analysis. For the simultaneous estimation of temperature, absorption and scattering coefficient profiles in the medium, the sum of residuals of boundary intensity and temperature after being weighted by a balance factor is minimized through using a Newton-type iteration algorithm and the least-squares method. To avoid over-updating for the parameters, the relative updating magnitude during the iteration process is constrained not to be >0.5. It is shown that the boundary intensity measurement alone is not enough to estimate simultaneously the temperature (source) and the radiative properties (both absorption and scattering coefficients) when the measurement data contain sensitive random errors. The boundary temperature measurement can serve as a necessary supplementation to the boundary intensity to make this kind of inverse radiative transfer problem resolvable. It was shown that a compensation relationship between absorption and scattering coefficients makes it difficult to fix them accurately. Parabolic profiles for the three parameters are used to validate the estimation method. When the optical thickness approaches 4.0, the results for the radiative properties are not acceptable, although the result for temperature profile is reasonable. This means the method needs further improvements.     

Discontinuous finite element method for radiative heat transfer in semitransparent graded index medium

To avoid the complicated and time-consuming computation of curved ray trajectories, a discontinuous finite element method based on discrete ordinate equation is extended to solve the radiative transfer problem in a multi-dimensional semitransparent graded index medium. Two cases of radiative heat transfer in two-dimensional rectangular gray semitransparent graded index medium enclosed by opaque boundary are examined to verify this discontinuous finite element method. Special layered and radial graded index distributions are considered. The predicted dimensionless net radiative heat fluxes and dimensionless temperature distributions are determined by the discontinuous finite element method and compared with the results obtained by the curved Monte Carlo method in references. The results show that the discontinuous finite element method has a good accuracy in solving the multi-dimensional radiative transfer problem in a semitransparent graded index medium.     

Transient coupled radiative and conductive heat transfer in a semitransparent layer of ceramic

This work considers transient radiative and conductive heat transfer in a semitransparent layer of ceramic, submitted to several thermal and radiative boundary conditions. Each side of the layer is exposed to hot or cold radiative surroundings, while each boundary is heated or cooled by convection. The solution procedure must provide accurate temperature distribution in the layer, so a nodal analysis based on Hottel's zonal method extended by ray tracing method is carried out. A finite difference method with non-uniform space and time increments is used to solve the transient energy equation, including a radiative heat source, coupled to a equation of radiative transfer. Variable spacing was used to concentrate grid points in regions with large temperature gradients. The influence of refractive index, optical thicknesses and conduction-radiation parameters is investigated.     

Attractive and repulsive Casimir vacuum energy with general boundary conditions

The infrared behaviour of quantum field theories confined in bounded domains is strongly dependent on the shape and structure of space boundaries. The most significant physical effect arises in the behaviour of the vacuum energy. The Casimir energy can be attractive or repulsive depending on the nature of the boundary. We calculate the vacuum energy for a massless scalar field confined between two homogeneous parallel plates with the most general type of boundary conditions depending on four parameters. The analysis provides a powerful method to identify which boundary conditions generate attractive or repulsive Casimir forces between the plates. In the interface between both regimes we find a very interesting family of boundary conditions which do not induce any type of Casimir force. We also show that the attractive regime holds far beyond identical boundary conditions for the two plates required by the Kenneth–Klich theorem and that the strongest attractive Casimir force appears for periodic boundary conditions whereas the strongest repulsive Casimir force corresponds to anti-periodic boundary conditions. Most of the analysed boundary conditions are new and some of them can be physically implemented with metamaterials.