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

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

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » fluent中uds的对流项
【调剂】北京石油化工学院2024年16个专业接受调剂
131/2返回列表12下一页
查看: 1676  |  回复: 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美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 科研好难啊，想退学去打工。 +10 byywnyl 2024-04-27 11/550 2024-04-28 11:21 by chemhua
[找工作] 普通院校药学硕士，做合成的，感觉找不到工作 +14 pom戴墨镜 2024-04-24 23/1150 2024-04-28 11:02 by pom戴墨镜
[有机交流] 锂盐的合成 4+3 liurui1203 2024-04-27 3/150 2024-04-28 10:45 by 591950582
[论文投稿] 求论文投稿期刊推荐 15+3 mihudawang 2024-04-26 6/300 2024-04-28 10:16 by 钰波love
[考博] 申博成果界定是根据Jcr分区还是中科院分区 +5 我属驴核动力驴 2024-04-22 6/300 2024-04-28 09:41 by cuan118
[有机交流] 如何分离原料与产品 5+4 qwerasdf587 2024-04-24 16/800 2024-04-28 09:22 by 光脚板bbv
[论文投稿] 关于Journal of Crystal Growth的OA和版面费 50+3 席羊羊 2024-04-26 4/200 2024-04-28 09:06 by 西大刘栋
[硕博家园] 聊天 +13 暮色恋伊人 2024-04-22 14/700 2024-04-27 23:26 by @tanzelin
[硕博家园] 失之东隅 收之桑榆 +5 半世迷离 2024-04-27 10/500 2024-04-27 22:19 by 半世迷离
[硕博家园] 博士白读了 +49 Da_Meng_Zi 2024-04-21 54/2700 2024-04-27 20:12 by 慕容家族
[找工作] 江苏理工学院要慎选 +3 jjchenshui 2024-04-27 3/150 2024-04-27 17:45 by 阿司匹林888
[考博] 真的好想读博！ +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
[论文投稿] 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 我属驴核动力驴
[基金申请] "颜宁:基础研究应顶天立地"能做到基础研究同时顶天立地的才是牛人 +6 zju2000 2024-04-24 6/300 2024-04-27 10:17 by kcmn1000
[考研] 学硕专硕 +5 小蜗牛* 2024-04-26 5/250 2024-04-26 16:43 by 鱼翔浅底1
[考研] 381求调剂 +4 小刺猬987654321 2024-04-25 6/300 2024-04-26 10:57 by czl12138
[考博] 博士招生 +4 zx179 2024-04-24 7/350 2024-04-24 15:01 by H考研成功
[教师之家] 大家访学都是怎么找的啊？ +3 luokereng 2024-04-22 3/150 2024-04-24 11:40 by xuechenli
信息提示
关闭
请填处理意见
关闭