壁面热边界条件中辐射边界的问题

wz9533

一般的壁面热边界条件中有辐射传热这一项，如图

                               
登录/注册后可看大图

，
但我没有直接使用辐射传热这一选项，而是按照fluent帮助文件中对于此项的描述

                               
登录/注册后可看大图

编写了如下的udf，并将此udf作为热流密度添加

                               
登录/注册后可看大图

，
但计算发散。

                               
登录/注册后可看大图

程序如下：
# include "udf.h"
DEFINE_PROFILE(Radiation_heatflux,t,i)
{
        face_t f;
        begin_f_loop(f,t)
        {
                F_PROFILE(f,t,i) = 5.67e-8 * (1.e12 - F_T(f,t)*F_T(f,t)*F_T(f,t)*F_T(f,t));
        }
        end_f_loop(f,t)
}

求助各位，是不是我的方法有错误，还是哪里设置错了，我觉得这两种方法应该得到相同的结果啊。

lsyzlq

你少写了ds[eqn]，这个是源项的导数

wz9533

我这个不是源项，是边界的代码

原帖由 lsyzlq 于 2013-9-2 17:03 发表

                               
登录/注册后可看大图

你少写了ds[eqn]，这个是源项的导数