| 查看: 1753 | 回复: 10 | ||
xxkobett新虫 (初入文坛)
|
[求助]
气液相变模型，可是每次初始化都提示Segmentation Fault，新人求解答
|
|
#include "udf.h" #include "sg_mphase.h" #include "mem.h" #define T_SAT 373.15 #define LAT_HT 2.257e6 DEFINE_SOURCE(vap_src, cell, pri_th, dS, eqn) { Thread *mix_th, *sec_th; real m_dot_v; mix_th = THREAD_SUPER_THREAD(pri_th); sec_th = THREAD_SUB_THREAD(mix_th, 1); if(C_T(cell, mix_th)<=T_SAT) { m_dot_v = -2*C_K_L(cell,sec_th)*NV_DOT(C_VOF_G(cell,sec_th),C_T_G(cell,mix_th))/LAT_HT; dS[eqn] = 0; } else { m_dot_v = 2*C_K_L(cell,pri_th)*NV_DOT(C_VOF_G(cell,pri_th),C_T_G(cell,mix_th))/LAT_HT; dS[eqn] = 0; } return m_dot_v; } DEFINE_SOURCE(lip_src, cell, sec_th, dS, eqn) { Thread * mix_th, *pri_th; real m_dot_l; mix_th = THREAD_SUPER_THREAD(sec_th); pri_th = THREAD_SUB_THREAD(mix_th, 0); if(C_T(cell, mix_th)<=T_SAT) { m_dot_l = 2*C_K_L(cell,sec_th)*NV_DOT(C_VOF_G(cell,sec_th),C_T_G(cell,mix_th))/LAT_HT; dS[eqn] = 0; } else { m_dot_l = -2*C_K_L(cell,pri_th)*NV_DOT(C_VOF_G(cell,pri_th),C_T_G(cell,mix_th))/LAT_HT; dS[eqn] = 0; } return m_dot_l; } DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn) { Thread *pri_th, *sec_th; real m_dot; pri_th = THREAD_SUB_THREAD(mix_th, 0); sec_th = THREAD_SUB_THREAD(mix_th, 1); if(C_T(cell, mix_th)<=T_SAT){ m_dot = 2*C_K_L(cell,sec_th)*NV_DOT(C_VOF_G(cell,sec_th),C_T_G(cell,mix_th))/LAT_HT; dS[eqn] = 0; } else { m_dot = -2*C_K_L(cell,pri_th)*NV_DOT(C_VOF_G(cell,pri_th),C_T_G(cell,mix_th))/LAT_HT; dS[eqn] = 0; } return LAT_HT*m_dot; }@wuming524 |
回复此楼» 猜你喜欢
青基如何破局
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
基金不中，共勉
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有11人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有7人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有6人回复
-
怎么看青基中了没有啊
已经有5人回复
-
梦想
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【趣闻】杨振宁撰文：《李政道传》有多处失实(1)
已经有41人回复
xxkobett
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 14.8
- 散金: 10
- 帖子: 13
- 在线: 3.2小时
- 虫号: 5129008
- 注册: 2016-10-17
- 专业: 工程热物理相关交叉领域
2楼2016-12-06 11:35:23
xxkobett
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 14.8
- 散金: 10
- 帖子: 13
- 在线: 3.2小时
- 虫号: 5129008
- 注册: 2016-10-17
- 专业: 工程热物理相关交叉领域
3楼2016-12-06 12:38:01
xxkobett
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 14.8
- 散金: 10
- 帖子: 13
- 在线: 3.2小时
- 虫号: 5129008
- 注册: 2016-10-17
- 专业: 工程热物理相关交叉领域
4楼2016-12-07 08:39:30
5楼2016-12-07 10:14:29
xxkobett
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 14.8
- 散金: 10
- 帖子: 13
- 在线: 3.2小时
- 虫号: 5129008
- 注册: 2016-10-17
- 专业: 工程热物理相关交叉领域
6楼2016-12-09 16:46:32
7楼2016-12-09 17:51:02
|
老版本软件里，并没有c_vof_g这个宏，需编程自己求解。其他几个_g 老版本里同样是求解器默认不保存这些数据的，以节省内存，不过可以在solver/set/expert 下用命令设置保存这几个梯度宏 发自小木虫Android客户端 |
8楼2016-12-09 22:13:19
wendong1993
新虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 59.5
- 散金: 19
- 红花: 1
- 帖子: 49
- 在线: 10.5小时
- 虫号: 1725971
- 注册: 2012-03-29
- 性别: GG
- 专业: 工程热物理相关交叉领域
9楼2016-12-12 10:55:00
10楼2017-08-07 00:24:00
10