半透明介质内辐射导热耦合换热的加速迭代模型 安巍.docVIP

中国工程热物理学会 传热传质学 编号：143105 半透明介质内辐射导热耦合换热的加速迭代模型 安巍朱彤(同济大学，　) (Tel: 021 anwei@tongji.edu.cn) 摘 要 半透明介质中辐射导热耦合换热问题的数值求解往往需要大量的迭代计算，本文通过对此类问题中迭代低效的原因进行分析，提出了一个基于温度场修正和投射辐射力修正的加速迭代计算模型，用于提高迭代求解效率，并通过两个算例检验了该算法的性能，结果表明：当前的迭代加速计算模型可以在不改变计算精度的前提下，大幅减少迭代次数和计算时间。特别是当壁面发射率较小，对流换热系数较小的情况下，其加速迭代效果更明显。 关键词 半透明介质；辐射导热耦合； (1) 对于灰体介质，则上述方程中将出现温度的四次方项，在数值求解中必须对非线性项进行线性化处理。一般地，可采用Taylor公式对其进行一阶近似。能量方程(1)可变为如下的形式： (2) 其中，下标n代表迭代计算次数，投射辐射力定义如下： (3) 式(2)的边界条件为： 当 (4-a) 当 (4-b) 当 (4-c) 其中，为边界上的辐射热流，可通过(3)式求得， (5) 为边界的内法向量，表示方向与的内积。为介质的吸收系数，当介质为非灰体时是波长的函数。为定温边界条件，为定热流边界条件，为对流换热边界条件，为边界上的对流换热系数，为外部流体的温度。由式(2)—(5)可见，能量方程和其边界条件中包含有未知量辐射强度，其是位置和方向的函数，可以通过式(6)所示的辐射传输方程(Radiative Transfer Equation, RTE)来求解， (6) 如果壁面为漫射边界，则边界条件如式7所示， (7) 从(6)—(7)式可以发现：式6的右端包含着与温度场有关的发射源项，而式7的边界条件中也包含发射项。因此，能量方程与RTE相互耦合，必须采用迭代方法进行求解。 求解导热辐射耦合传热问题最通常的迭代方法称为源项迭代(Source Iteration, SI)。即先假设温度场，代入RTE中求出辐射强度，计算能量方程中的投射辐射和边界上的辐射热流，再代入能量方程求出新的温度场，如此反复直到温度场的误差满足设定的收敛判据。然而，当采用基于DOM的各类数值方法求解RTE时，由于温度场未知，式(6)和(7)中与温度有关的发射项并不能在每次的迭代中被精确的计算，由此导致能量方程(2)中的投射辐射及边界条件(4)中的辐射热流也不能精确获得，这样的误差在能量方程和RTE中反复传递，使得源项迭代收敛缓慢，特别是对光学厚度较大或散射占优介质，其迭代效率进一步恶化。 基于以上的分析可知，如果介质温度和投射辐射力能够被适当地朝准确值修正，减小温度场和辐射场中误差的累积和传递效应，则迭代效率将获得提升。因此，温度场与投射辐射力的修正方法至关重要。 1.1 温度场的修正 为了获得介质温度的修正值，可以假设为介质温度场的准确值，为迭代中的一个近似值，它们两者的差为温度场的修正值，即。事实上，温度的不准确是由迭代中和的不准确导致的。和分别满足方程2，则有： (8-a) (8-b) 上式的两边分别相减，可以得到： (9) 其中，为投射辐射力的不准确度，边界辐射热流的不准确度也可表达为和函数。因此，式(9)的边界条件为： 当 (10-a) 当 (10-b) 当 (10-c) 通过求解方程(9)，温度场可获得修正。 1.2 投射辐射力的修正 对于投射辐射力的修正，可以采用扩散综合加速法(Diffusion Synthetic Acceleration, DSA)获得。DSA的核心思想是在每次迭代的过程中，通过低阶的P1扩散方程获得发射项（或散射项）的修正值，从而提升迭代效率，DSA方法的详细介绍可参考文献[1-3]。当和都满足P1扩散近似方程，则由于发射项的不准确导致的投射辐射力修正可由下式(11)求得：