间断有限元法求解一维矢量辐射传输

采用间断有限元法(Discontinuous finite element method,DFEM)求解一维散射性介质内矢量辐射传输问题。推导了散射性介质内矢量辐射传输方程的间断有限元离散格式,空间离散采用间断有限元离散将求解域划分为相互独立的单元,角度离散在传统的均匀分段光滑近似(Piecewise Constant Approximation,PCA)角度离散基础上进行局部加密以得到关键方向上的数值解。采用两个稳态矢量辐射传输算例对间断有限元法求解矢量辐射传输方程的正确性和计算效率进行了验证,在此基础上拓展间断有限元法应用于求解散射性介质内瞬态矢量辐射传输方程,将随着时间推移过程中的瞬态辐射信息和最终的稳态结果相比较验证了本文方法的正确性。     

梯度折射率介质瞬态辐射-导热耦合传热

多维辐射-导热耦合传热问题有广泛的应用背景。本文采用有限元法对辐射传递方程的离散坐标形式进行离散,并采用有限元法求解瞬态辐射-导热耦合能量方程,研究二维方腔半透明参与性梯度折射率介质内瞬态耦合传热,考虑第一类热边界条件,分析了折射率分布、衰减系数、散射反照率及散射相函数对对瞬态耦合传热的影响。结果表明,梯度折射率对瞬态耦合传热影响显著。     

基于微分单元的圆柱内辐射导热耦合传热分析

本文充分利用微分单元法不涉及积分、较为稳定、易于实施的特点,进行圆柱坐标系下辐射与导热耦合传热分析.在求解过程中,将热辐射的空间和角度耦合成三维计算域求解,实现空间和角度的高阶离散.针对具有强对流特性的辐射传递方程,采用一种迎风格式可有效抑制数值振荡,实现稳定计算.由于辐射边界条件与导热和对流边界条件不同,它是具有方向性的单向边界条件,且在边界处会存在强烈的间断性和奇异点.为了有效抑制间断并处理奇异点,本文提出双层间断节点方案.通过与文献中圆柱坐标系下辐射与导热耦合传热算例进行对比分析,发现微分单元模型在空间和角度方向上均能实现高分辨率刻画,获得较好的计算精度.     

有限体积法求解圆柱形散射介质内辐射与导热耦合换热

将谱带模型与有限体积解法相结合;求吸收、发射、散射性非灰介质圆柱体内辐射传递方程。考虑辐射强度场与热扩散温度场的耦合,将控制容积法与有限体积法结合,求解辐射与导热耦合换热。经与光线踪迹法、离散传递法的计算结果比较表明,谱带模型与有限体积解法相结合能处理多场耦合下非灰介质内的辐射换热。     

氧化铝纤维传热方式及反问题研究

本文分析了典型热防护绝热材料氧化铝纤维内部的传热方式及传热反问题。氧化铝纤维内部的传热主要以固体纤维导热、纤维孔隙气体导热和辐射传热为主。本文以气-固耦合导热模型和辐射模型为依据,利用数值方法求解了氧化铝纤维内部的温度场,并拟合得出耦合气体导热,固体导热和辐射三种方式的有效热导率。本文利用遗传算法求解氧化铝纤维的稳态传热反问题,根据氧化铝纤维的温度数据,优化得到有效热导率的预测值,并与参考值进行对比,结果表明在一定进化代数内优化误差随进化代数的增加而减小。     

三维立方体内辐射磁流体流动与传热研究

采用配置点谱方法求解了在存在壁面温差的立体方腔内有外加磁场作用下热辐射对磁流体流动与传热的影响。采用谱投影算法处理速度和压力的耦合,辐射传递方程的角向离散采用球带等差数列微元等分离散坐标(SRAPN)法进行离散,最后在相同的节点下采用配置点谱方法进行求解。在验证数值方法之后,分析了Ha数、介质光学厚度、以及导热辐射参数对等温面、流线投影等的详细影响规律。     

隐-显积分因子间断Galerkin方法求解二维辐射扩散方程

提出一种求解二维非平衡辐射扩散方程的数值方法.空间离散上采用加权间断Galerkin有限元方法,其中数值流量的构造采用一种新的加权平均;时间离散上采用隐-显积分因子方法,将扩散系数线性化,然后用积分因子方法求解间断Galerkin方法离散后的非线性常微分方程组.数值试验中在非结构网格上求解了多介质的辐射扩散方程.结果表明:对于强非线性和强耦合的非线性扩散方程组,该方法是一种非常有效的数值算法.     

离散统一气体动理学格式稳态辐射输运应用

工程应用中的介质热辐射问题是典型的多尺度问题. 基于Boltzmann输运方程建立的各类气体动理学格式, 在多尺度瞬态问题中得到了广泛应用. 为了克服显式求解方案中CFL条件等的限制, 文章通过气体动理学格式实现稳态辐射输运方程的直接求解. SDUGKS格式由离散统一气体动理学格式(discrete unified gas kinetic scheme, DUGKS)的核心思想发展而来, 应用于稳态问题计算. 将SDUGKS格式进一步拓展到多尺度的稳态热辐射输运计算. SDUGKS格式继承了DUGKS格式沿特征线离散实现的界面重构, 并通过隐式增量格式的单元更新实现对辐射强度的较正, 采用逐次迭代法将辐射强度渐近收敛到稳定值. 选用多组一维和二维不同尺度的辐射传热算例, 通过与特定的解析解以及其他数值方法结果对比, 检验了SDUGKS的计算精度和计算效率, 并论证了它在多尺度问题中的渐进保持性质.      

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

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

高温空气燃烧炉内耦合传热的数值模拟

用自行开发的三维湍流流动、燃烧、辐射传热和NO_x湍流生成的计算程序,用离散坐标方法(DOM)模拟炉内的辐射传热并与其它传热方式相耦合,与各个输运方程共同求解,对燃烧室内的温度场进行了数值模拟.通过改变高温空气的预热温度,数值分析了高温空气燃烧炉内采用燃气直接喷射技术(FDI)的温度分布特性,检验了所采用的耦合数值模拟技术的有效性.模拟计算结果表明,提高空气预热温度,炉内温度的峰值相应增大,温度梯度降低,温度分布更均匀,火焰更长.相关的实验结果与数值预报的结果对比表明了相一致的规律.     

Numerical analyses of radiative heat transfer in any arbitrarily-shaped axisymmetric enclosures

A numerical approach for the treatment of radiative heat transfer in any irregularly-shaped axisymmetric enclosure filled with absorbing, emitting and scattering gray media is developed. Radiative transfer equation (RTE) is formulated for a general axisymmetric geometrical configurations, and the discretized equation is conducted using an unstructured meshes, generated by an appropriate computer algorithm, and the control volume finite element method which frequently adopted in CFD problems. A computer procedure has been done to solve the discretized RTE and to examine the accuracy and the computational efficiency of the proposed numerical approach. By using this computer algorithm, five test cases, a cylindrical enclosure with absorbing and emitting medium, a diffuser shaped axisymmetric enclosure, a finite axisymmetric cylindrical enclosure with a curved wall, a furnace with axially varying medium temperature and a rocket nozzle, are treated and the obtained results agree very well with other published works. Furthermore, the developed computer procedure has an accurate CPU time and it can be coupled easily with CFD codes.     

多孔介质快速干燥过程热质耦合方程的代数显式解析解

对多孔介质快速干燥过程的传热与传质耦合方程组导出了两套代数显式解析特解。这些解首先可以作为计算传热传质学的标准解,用以检验数值计算的准确性、收敛性与稳定性等,还可以启发数值工作者改进计算技巧例如差分格式与网格生成技术等。当然,解析解还会有其相应的理论价值。     

Prediction of radiative heat transfer in 3D complex geometries using the unstructured control volume finite element method

In this paper, a 3D algorithm for the treatment of radiative heat transfer in emitting, absorbing, and scattering media is developed. The numerical approach is based on the utilization of the unstructured control volume finite element method (CVFEM) which, to the knowledge of the authors, is applied for the first time to simulate radiative heat transfer in participated media confined in 3D complex geometries. This simulation makes simultaneously the use of the merits of both the finite element method and the control volume method. Unstructured 3D triangular element grids are employed in the spatial discretization and azimuthal discretization strategy is employed in the angular discretization. The general discretization equation is presented and solved by the conditioned conjugate gradient squared method (CCGS). In order to test the efficiency of the developed method, several 3D complex geometries including a hexahedral enclosure, a 3D equilateral triangular enclosure, a 3D L-shaped enclosure and 3D elliptical enclosure are examined. The results are compared with the exact solutions or published references and the accuracy obtained in each case is shown to be highly satisfactory. Moreover, this approach required a less CPU time and iterations compared with those of even parity formulation of the discrete ordinates method.     

Generalized Modal Scattering Matrix of Bandpass Filters Using Ridged Waveguides

Simple and accurate integral equation formulation using a three dimensional finite element (3DFEM) to determine the propagation of microwave passive structures is presented. This technique is applied to solve the scattering problems due to various forms of irises in lossless rectangular waveguide.The analysis of these components is highly dependent of the numbers of modes considered in the neighbourhood of the discontinuity. The theoretical model is compared with measurements and Mode Matching Method (MMM). Results are found to be good agreement.     

Least-squares finite element formulations for one-dimensional radiative transfer

We present least-squares-based finite element formulations for the numerical solution of the radiative transfer equation in its first-order primitive variable form. The use of least-squares principles leads to a variational unconstrained minimization problem in a setting of residual minimization. In addition, the resulting linear algebraic problem will always have a symmetric positive definite coefficient matrix, allowing the use of robust and fast iterative methods for its solution. We consider space-angle coupled and decoupled formulations. In the coupled formulation, the space-angle dependency is represented by two-dimensional finite element expansions and the least-squares functional minimized in the continuous space-angle domain. In the decoupled formulation the angular domain is represented by discrete ordinates, the spatial dependence represented by one-dimensional finite element expansions, and the least-squares functional minimized continuously in space domain and at discrete locations in the angle domain. Numerical examples are presented to demonstrate the merits of the formulations in slab geometry, for absorbing, emitting, anisotropically scattering mediums, allowing for spatially varying absorption and scattering coefficients. For smooth solutions in space-angle domain, exponentially fast decay of error measures is demonstrated as the p-level of the finite element expansions is increased. The formulations represent attractive alternatives to weak form Galerkin finite element formulations, typically applied to the more complicated second-order even- and odd-parity forms of the radiative transfer equation.     

Finite element formulation of the first- and second-order discrete ordinates equations for radiative heat transfer calculation in three-dimensional participating media

Two finite element methods (FEMs), FEDOM1 and FEDOM2 (standing for the first and the second finite element discrete ordinates methods, respectively), are formulated and numerically tested. The reference second-order discrete equation is modified in its scattering terms and is applied to the problems of absorbing/emitting and anisotropically scattering media by using the FEM. Numerical features of the developed FEMs are compared with one of the discrete ordinates interpolation method (DOIM), which uses a finite difference scheme. Prediction results of radiative heat transfer by these two FEMs are compared with reference solutions and verified in three-dimensional enclosures containing participating media. The results of FEDOM1 and FEDOM2 agree well with exact solutions for the problem of absorbing/emitting medium with various range of optical thickness. Generally, the two FEMs show more accurate results than DOIM. And FEDOM1 shows more accurate results than FEDOM2 in most of the test problems. Both of the developed FEMs show reasonable results compared with published Monte Carlo solutions for the tested absorbing/emitting and anisotropically scattering media. Although the FEDOM2 is not as accurate as the FEDOM1, it shows its own advantages that it reduces CPU time and memory space of dependent variable to half.     

多孔复合介质周期结构热传导和质扩散问题的多尺度数值方法

本文针对一类复杂的多孔复合介质的热传导和质扩散问题,给出具体的多尺度渐近展开公式,并在此基础上设计了有限元算法格式,它是宏观和细观相结合的数值方法。理论分析和数值实验均表明：多尺度数值方法对求解多孔复合介质周期结构的热传导和质扩散问题是可行的和有效的。     

Three-dimensional radiative transfer modeling using the control volume finite element method

In the present study, a three-dimensional algorithm for the treatment of radiative heat transfer in emitting, absorbing and scattering media is developed. The approach is based on the utilization of control volume finite element method (CVFEM) which, to the knowledge of the authors, is applied at the first time to 3D radiative heat transfer in participating media. The accuracy of the present algorithm is tested by comparing its predictions to other published works. Comparisons show that CVFEM produces good results. Moreover, this approach permits compatibility with other numerical methods used for computational fluids mechanics problems.     

Étude numérique de l'influence du transfert radiatif sur la convection naturelle laminaire,bidimensionnelle, permanente,dans un espace annulaire d'axe horizontal,délimité par deux cylindres circulaires isothermes

Combined radiation and natural convection in a participating medium between concentric or vertically eccentric horizontal cylinders is investigated numerically. The annular medium is considered as a gray, emitting, absorbing, and isotropically scattering gas. The equations of steady, laminar, two-dimensional, thermal, natural convection are written by using a two-cylindrical coordinate system, the stream function, and the vorticity. The finite volumes method is used to discretize the coupled equations of momentum, energy, and radiative transfer. To solve the global nonlinear algebraic equations the successive-over-relaxation iterating scheme is applied. Numerical solutions are obtained for a Rayleigh number in the range 10³–10⁵ and radiation-conduction parameter ranging from 0 to ∞. The influences of radiation-conduction parameter, Rayleigh number and other parameters on flow and temperature distributions and heat transfer are discussed.