应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » fluent中uds的对流项
【调剂】北京石油化工学院2024年16个专业接受调剂
131/2返回列表12下一页
查看: 1669  |  回复: 12
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者赶紧发sci将赠送您 50 个金币

赶紧发sci

新虫 (小有名气)

[求助] fluent中uds的对流项

小弟最近在做电弧双温模型，两个能量方程都是通过uds写的，程序写好后导入fluent计算几千步就发散了，速度入口附近的速度和压力变得非常大，怀疑是对流项的问题，但是不太会修改，请各位大大帮忙看一下，万分感激！小木虫里不会发图片，发不了公式，这里先附上我的代码：
DEFINE_UDS_FLUX(uds3_flux, f, t, i)
{
        cell_t c0, c1 = -1;
        Thread *t0, *t1 = NULL;
        real NV_VEC(psi_vec);
        real NV_VEC(V);
        real NV_VEC(A), flux = 0.0;
        real ke = 1.38e-23;
        c0 = F_C0(f, t);
        t0 = F_C0_THREAD(f, t);        //内部面，气体流入的为C0，流出为C1。边界面上只有C0没有C1
        F_AREA(A, f, t);
        /* If face lies at domain boundary, use face values; */
        /* If face lies IN the domain, use average of adjacent cells. */
        if (BOUNDARY_FACE_THREAD_P(t)) /*Most face values will be available*/
        {
                real dens;

                /* Depending on its BC, density may not be set on face thread*/
                if (NNULLP(THREAD_STORAGE(t, SV_DENSITY)))
                {
                        dens = edensity(F_UDSI(c0, t0, 3), F_UDSI(c0, t0, 4));  /* Set dens to face value if available */
                        C_UDMI(c0, t0, 16) = edensity(F_UDSI(c0, t0, 3), F_UDSI(c0, t0, 4));
                }
                else
                {
                        dens = edensity(C_UDSI(c0, t0, 3), C_UDSI(c0, t0, 4)); /* else, set dens to cell value */
                        C_UDMI(c0, t0, 17) = edensity(C_UDSI(c0, t0, 3), C_UDSI(c0, t0, 4));
                }
                psi_vec[0] = F_U(f, t)*2.5*dens*ke;
                psi_vec[1] = F_V(f, t)*2.5*dens*ke;
                /*NV_DS(psi_vec, =, F_U(f,t), F_V(f,t), F_W(f,t), *, dens*ke);*/
                flux = NV_DOT(psi_vec, A); /* flux through Face */
               
        }
        else
        {

                c1 = F_C1(f, t);   /* Get cell on other side of face */
                t1 = F_C1_THREAD(f, t);
                psi_vec[0] = C_U(c0, t0)*2.5*ke*edensity(C_UDSI(c0, t0, 3), C_UDSI(c0, t0, 4));
               
                C_UDMI(c0, t0, 18) = C_U(c0, t0)*2.5*ke*edensity(C_UDSI(c0, t0, 3), C_UDSI(c0, t0, 4));
                C_UDMI(c1, t1, 19) = C_V(c1, t1)*2.5*ke*edensity(C_UDSI(c1, t1, 3), C_UDSI(c1, t1, 4));
                psi_vec[1] = C_V(c1, t1)*2.5*ke*edensity(C_UDSI(c1, t1, 3), C_UDSI(c1, t1, 4));
                flux = NV_DOT(psi_vec, A)/2.0; /* Average flux through face */
                C_UDMI(c0, t0, 23) = flux;
               
        }
        /* ANSYS Fluent will multiply the returned value by phi_f (the scalar’s
        value at the face) to get the ‘‘complete’’ advective term. */
        C_UDMI(c0, t0, 24) = flux;
        return flux;
}@wuming524
回复此楼

» 本帖已获得的红花（最新10朵）

whaleof52HZ

» 猜你喜欢
1楼2021-01-24 16:14:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

whaleof52HZ

新虫 (小有名气)
大佬能给个联系方式吗，udf和uds你是咋学的

发自小木虫Android客户端
赞一下(1人) 回复此楼
3楼2021-01-30 23:06:45
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赶紧发sci

新虫 (小有名气)
引用回帖:
3楼: Originally posted by whaleof52HZ at 2021-01-30 23:06:45
大佬能给个联系方式吗，udf和uds你是咋学的

我只是个小菜鸟，，，主要就是看网上流传的udf中文教程，小木虫上很多，然后就论坛上问答中一些对应的程序了解下
赞一下(1人) 回复此楼
4楼2021-02-03 11:41:05
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

whaleof52HZ

新虫 (小有名气)
送红花一朵


发自小木虫Android客户端
2楼2021-01-30 23:05:57
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赶紧发sci

新虫 (小有名气)
有人吗，求助呀
回复此楼
5楼2021-02-03 13:05:57
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赶紧发sci

新虫 (小有名气)
qiuqiuqiu
回复此楼
6楼2021-02-07 16:13:56
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whaleof52HZ

新虫 (小有名气)
引用回帖:
4楼: Originally posted by 赶紧发sci at 2021-02-03 11:41:05
我只是个小菜鸟，，，主要就是看网上流传的udf中文教程，小木虫上很多，然后就论坛上问答中一些对应的程序了解下...

电弧双温模型和超音速等离子弧一样吗

发自小木虫Android客户端
赞一下 回复此楼
7楼2021-02-15 20:11:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赶紧发sci

新虫 (小有名气)
Emmm

发自小木虫IOS客户端
回复此楼
8楼2021-02-21 11:51:20
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赶紧发sci

新虫 (小有名气)
引用回帖:
7楼: Originally posted by whaleof52HZ at 2021-02-15 20:11:36
电弧双温模型和超音速等离子弧一样吗
...

我认为基本方程都是一样的吧，用磁流体模型，就是边界条件之类的区别大一点？

发自小木虫IOS客户端
赞一下 回复此楼
9楼2021-02-21 11:52:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

破下限

新虫 (著名写手)
我用的帮助文档上一模一样的UDF也发散了
赞一下 回复此楼
10楼2021-04-25 22:04:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 赶紧发sci 的主题更新
131/2返回列表12下一页
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 科研好难啊，想退学去打工。 +5 byywnyl 2024-04-27 6/300 2024-04-27 23:13 by byywnyl
[考博] 博士招生-211农业院校 +5 NYC917 2024-04-26 10/500 2024-04-27 22:31 by NYC917
[教师之家] 博士论文被抄袭 +27 和尚敲小木鱼 2024-04-22 44/2200 2024-04-27 21:57 by ay_guobo
[论文投稿] with editor 两个月了，什么原因？ +8 yiersan9 2024-04-24 20/1000 2024-04-27 21:12 by yiersan9
[考博] 真的好想读博！ +16 wangzhe_bs 2024-04-22 23/1150 2024-04-27 17:10 by @tanzelin
[留学DIY] 10043以后咋出去啊 +4 RuiLove 2024-04-22 4/200 2024-04-27 16:03 by 鱼翔浅底1
[论文投稿] 研二光催化6月底四篇二区什么水平 5+5 wjtab 2024-04-22 16/800 2024-04-27 15:47 by cuan118
[论文投稿] Nature一直在编辑手里，考虑好几天了，是悬了吗 +13 彩虹初见 2024-04-24 13/650 2024-04-27 15:24 by cuan118
[考博] 25年博士申请 +7 Changzixuan 2024-04-25 13/650 2024-04-27 14:11 by 我属驴核动力驴
[基金申请] 两类问题算是白选了~ +8 jurkat.1640 2024-04-23 13/650 2024-04-27 12:03 by 淀粉搬运工
[论文投稿] 投稿RSC旗下杂志突然间看不到投稿状态，里面投稿的文章信息不见了，有遇到过的吗？ 50+3 dlying 2024-04-26 3/150 2024-04-27 09:58 by bobvan
[论文投稿] AC审稿 +3 马winner 2024-04-26 6/300 2024-04-27 09:08 by 马winner
[有机交流] 如何分离原料与产品 5+4 qwerasdf587 2024-04-24 13/650 2024-04-26 23:18 by 119966984
[考博] 申博求助 +4 dskabdh 2024-04-24 11/550 2024-04-26 15:54 by dskabdh
[硕博家园] 考研，求职还是考编? +15 xizj 2024-04-21 24/1200 2024-04-26 11:49 by Kan客
[论文投稿] 一直找不到审稿人 +5 lizhengke06 2024-04-21 6/300 2024-04-25 14:01 by chongdong
[基金申请] 国社科项目，你们学校都限额申报吗？ +7 屡战屡败 2024-04-21 10/500 2024-04-25 12:10 by 屡战屡败
[考博] 博士招生 +4 zx179 2024-04-24 7/350 2024-04-24 15:01 by H考研成功
[论文投稿] 期刊推荐 20+4 木颜尘ip 2024-04-22 7/350 2024-04-24 10:06 by bobvan
[考博] 申博成果界定是根据Jcr分区还是中科院分区 +4 我属驴核动力驴 2024-04-22 5/250 2024-04-24 08:47 by 晓目崇
信息提示
关闭
请填处理意见
关闭