| 查看: 1727 | 回复: 12 | ||
| 【悬赏金币】回答本帖问题，作者赶紧发sci将赠送您 50 个金币 | ||
[求助]
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» 猜你喜欢
数理学部函评几号结束？
已经有4人回复
-
有博导缺学生吗，我可以当牛做马
已经有53人回复
-
浙江的高校现在门槛都这么高吗
已经有6人回复
-
民办高校入职后稳定吗？
已经有6人回复
-
南京林业大学老师首聘期考核不过，服药自杀
已经有33人回复
-
投稿成功后又想撤回
已经有3人回复
-
博士去高校就是为了有寒暑假吗？
已经有14人回复
-
求期刊推荐
已经有4人回复
-
一个审稿人已返回，也没邀请审稿人，也没人审稿中，编辑要闹啥？
已经有8人回复
-
基金委也是打工的
已经有6人回复
whaleof52HZ
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 703.1
- 散金: 150
- 沙发: 2
- 帖子: 270
- 在线: 10.9小时
- 虫号: 10258417
- 注册: 2018-09-12
- 性别: GG
- 专业: 可再生与替代能源利用中的
3楼2021-01-30 23:06:45
4楼2021-02-03 11:41:05
whaleof52HZ
新虫 (小有名气)
2楼2021-01-30 23:05:57
5楼2021-02-03 13:05:57
6楼2021-02-07 16:13:56
whaleof52HZ
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 703.1
- 散金: 150
- 沙发: 2
- 帖子: 270
- 在线: 10.9小时
- 虫号: 10258417
- 注册: 2018-09-12
- 性别: GG
- 专业: 可再生与替代能源利用中的
7楼2021-02-15 20:11:36
8楼2021-02-21 11:51:20
9楼2021-02-21 11:52:40
破下限
新虫 (著名写手)
- 应助: 3 (幼儿园)
- 金币: 6422.2
- 散金: 86
- 红花: 6
- 帖子: 1152
- 在线: 367.3小时
- 虫号: 3639947
- 注册: 2015-01-11
- 性别: GG
- 专业: 可再生与替代能源利用中的
10楼2021-04-25 22:04:49
4