24小时热门版块排行榜    

查看: 1148  |  回复: 1
【悬赏金币】回答本帖问题,作者tsv249052将赠送您 500 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

tsv249052

新虫 (著名写手)

[求助] 求助,UDS源项,一计算就发散是什么问题已有1人参与

求助,我自己编的udf程序,使用了define adjust和uds源项,在fluent中编译build和加载load都没有问题,但是计算一开始就显示uds-0发散了,不知道是什么原因,部分代码如下:
#include "udf.h"

#define PI 3.14159 /* 定义π*/
#define K 1.38e-23 /*定义玻尔兹曼常数 */
#define M 2.99e-26 /*定义水分子质量*/
#define R 461.5 /*定义水蒸气气体常数*/
#define f_v 0.4 /*相对湿度*/

real S1, G, J, pv1, ps1, T, P, sur_tension, density_c, density_v;/*不饱和度,吉布斯焓,成核率,水蒸气分压力,水蒸气饱和压力,温度,压力*/

/*计算成核率,液滴生长速度,临界半径*/
DEFINE_ADJUST(chenghelv, d)
{
        Thread *t;
        cell_t c;

        thread_loop_c (t,d)
        {
        begin_c_loop_all(c,t)
        {
        T = C_T(c,t);
        P = C_P(c,t);
        ps1 = -2.88934 + 0.37966 * T - 0.00925 * pow(T, 2) + 1.54e-4 * pow(T, 3);/*饱和压力拟合*/
        density_v = -0.00134 + 0.00136 * T - 3.27e-5 * pow(T, 2) + 7.76e-7 * pow(T, 3);/*蒸气密度拟合*/
        density_c = 1000.07904 + 0.01231 * T - 0.00586 * pow(T,2) + 1.6e-5 * pow(T,3);/*水密度拟合*/
        sur_tension = 0.07566 - 1.41e-4 * T - 2.68e-7 * pow(T,2);/*表面张力拟合*/
        S1 = P / ps1;
        pv1 = P * f_v;
        G = 16 / 3 * PI * pow((M / (density_c * log(S1) * K * T)),2) * pow(sur_tension, 3); /*计算吉布斯自由焓*/
        C_UDMI(c,t,0) = sqrt(2 / PI * sur_tension) * pow(M ,-3/2) * pow(density_v, 2) / density_c * exp( -G / K * T);/*成核率*/
        end_c_loop_all(c,t)
        }
                }
}

/*第一个输送方程Q0源项*/
DEFINE_SOURCE(q0_source,c,t,dS,eqn)
{
        return C_UDMI(c,t,0);
}

/*第二个输送方程Q1源项*/
DEFINE_SOURCE(q1_source,c,t,dS,eqn)
{
        return  C_UDMI(c,t,0) * pow(C_UDMI(c,t,2),2) + C_R(c,t) * C_UDSI(c,t,0) * C_UDMI(c,t,1);
}

/*第三个输送方程Q2源项*/
DEFINE_SOURCE(q2_source,c,t,dS,eqn)
{
        return  C_UDMI(c,t,0) * pow(C_UDMI(c,t,2),2) + 2 * C_R(c,t) * C_UDSI(c,t,1) * C_UDMI(c,t,1);
}

/*第四个输送方程g的源项*/
DEFINE_SOURCE(g_source,c,t,dS,eqn)
{
        return 4 * PI * density_c * (1/3 * C_UDMI(c,t,0) * pow(C_UDMI(c,t,2),3) + C_R(c,t) * C_UDSI(c,t,2) * C_UDMI(c,t,1));
}

/*能量方程源项*/
DEFINE_SOURCE(energy_source,c,t,dS,eqn)
{
        real dt, temp, L;
        dt = RP_Get_Real("physical-time-step";
        temp = C_T(c,t);
        L = 2502 - 2.42169 * temp;
        return L * C_UDSI(c,t,3) / dt;
}
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Diligent雨

新虫 (初入文坛)

【答案】应助回帖

把udf简化尝试,运用排除法试试;可能是fluent内部设置问题,查看与uds0所有的设置;减小松弛因子。希望有帮助
焊接fluent熔池
2楼2020-03-23 22:01:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tsv249052 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 数理的人才答辩通知发了吗 +8 yzy3327 2024-05-30 8/400 2024-06-03 18:19 by xuxuxuxujie
[论文投稿] 希望嫡长子顺利出场! +6 C刊霸王 2024-05-31 7/350 2024-06-03 09:51 by 畅21
[论文投稿] 《控制理论与应用》期刊收版面费吗? +5 ygj2015 2024-05-29 5/250 2024-06-03 09:19 by lmpmsm
[论文投稿] 选期刊 5+3 jfdhj 2024-05-29 5/250 2024-06-03 07:47 by 黑大环境队长
[论文投稿] 求助大神,Fe和Al离子对MOF都有淬灭,当两种离子共存时,怎么区分两种离子? 10+3 maoxiao 2024-06-02 3/150 2024-06-03 01:43 by nono2009
[论文投稿] 编辑对文章提的意见 40+3 qa8386136 2024-06-01 3/150 2024-06-02 13:28 by bobvan
[论文投稿] 没收到邮件 10+4 荣小撇 2024-05-31 9/450 2024-06-02 13:25 by bobvan
[考博] 导师不让硕转博,让我去国外读博,能理解吗? +12 萧山幽谷 2024-05-29 20/1000 2024-06-02 12:01 by yuan0806
[考博] 申请2024或2025年博士研究生 +5 嘟噜嘟1 2024-05-29 12/600 2024-06-01 22:36 by 嘟噜嘟1
[教师之家] 在大地上我们只过一生---看完我的阿勒泰上头了好几天,完结那天晚上几乎失眠 +11 瞬息宇宙 2024-05-27 13/650 2024-06-01 22:15 by otani
[考博] 24or25材料专业申博 +4 农夫三拳有点痛 2024-05-30 11/550 2024-06-01 14:45 by Napoleonsky
[论文投稿] 求Sci期刊推荐 10+4 甄小鱼 2024-05-30 7/350 2024-06-01 10:41 by bobvan
[高分子] MMA预聚体光固化发雾问题求助 +3 惠亚金总 2024-05-29 10/500 2024-05-31 14:59 by 惠亚金总
[硕博家园] 求问endnote上Wiley引用格式在哪下载 +3 Elon. 2024-05-30 3/150 2024-05-31 10:13 by crazy peng
[材料综合] 真空封石英管 北京 +4 dessha 2024-05-29 5/250 2024-05-30 16:40 by mpdfwxgui
[基金申请] 信息学部函评结束了吗? +6 ducan21 2024-05-28 7/350 2024-05-29 12:10 by WORLD0256
[论文投稿] 真急着毕业,CPB主编终审17天了,邮件催稿了两次,就是一点动静没有 5+3 kkkk夏 2024-05-28 6/300 2024-05-29 11:18 by hitsdu
[论文投稿] 高手朋友推荐比较容易投稿和录用的SCI期刊,不投稿收费SC,对分区没有要求 5+3 xintangren 2024-05-28 4/200 2024-05-29 10:46 by xintangren
[论文投稿] 核心初审被拒,理由是“选题的意义不明确,文章写得不像是科技论文”,怎么改 5+3 工藤雷花樱 2024-05-27 8/400 2024-05-29 10:09 by topedit
[基金申请] 如果您是国自然评审专家 +3 丁香园账户 2024-05-28 3/150 2024-05-29 06:44 by gaohui8888
信息提示
请填处理意见