| 查看: 5657 | 回复: 10 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
1036985870金虫 (小有名气)
|
[求助]
fluent模拟传质过程 已有2人参与
|
|
|
各位做过传质模拟的坛友们大家好，我模拟的是氮气和二氧化碳在水中上升过程，气泡是我用patch实现的。以下是我编写的二氧化碳溶解过程的udf: #include "udf.h" #define H 0.39 /*定义二氧化碳溶解度系数*/ #define module_weight 44. /*定义二氧化碳的摩尔质量*/ #define density_CO2 1.7878 /*定义二氧化碳的密度*/ DEFINE_MASS_TRANSFER(liq_gas_transfer, cell, thread, from_index,from_species_index, to_index, to_species_index) { real m_lg,pressure_gas; Thread *gas, *liq; gas = THREAD_SUB_THREAD(thread, from_index); /*获取气相指针*/ liq = THREAD_SUB_THREAD(thread, to_index); /*获取液相指针*/ pressure_gas =C_R(cell,gas)/density_O2*C_YI(cell,gas,0)*C_P(cell,gas); /*获取气相中二氧化碳的压力*/ m_lg = module_weight*(H*pressure_gas-C_YI(cell,liq,0)*C_R(cell,liq)/module_weight)*log(H*pressure_gas/(H*pressure_gas-C_YI(cell,liq,0) *C_R(cell,liq)/module_weight))/(CURRENT_TIMESTEP+0.000000000001); /*获取氧气传质速率，单位为kg/m3s*/ return m_lg; } 迭代开始出现Error: Floating point error: invalid number Error Object: () 请问这种问题该怎么解决，我怀疑是存在0项，但是不知道问题出现在哪里？还请各位大神指点指点！另外，我在论坛看到有关传质过程的模拟，追踪气液界面的C_VOF(c,t)本来的范围就是0到1之间，有人说寻找两相间的界面的话，直接采用C_VOF(c,t)=0.5来判断，但是感觉是不对的。还有就是是不是mass transfer本身就包含了气液界面啊。。。 |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有18人回复
-
基金系统什么内容也没有
已经有9人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有11人回复
-
怎么查啊
已经有6人回复
-
国自然面上复盘~欢迎讨论
已经有12人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有6人回复
-
面上合作单位盖章
已经有5人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
申请删除本帖
已经有7人回复
-
基金不中，共勉
已经有10人回复
ljqmmljq
铁虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1613.3
- 帖子: 317
- 在线: 12.4小时
- 虫号: 13164584
- 注册: 2018-12-18
- 专业: 工程热物理与能源利用
11楼2019-01-17 10:55:05
1036985870
金虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 330.6
- 红花: 3
- 帖子: 113
- 在线: 76.9小时
- 虫号: 2355149
- 注册: 2013-03-17
- 性别: GG
- 专业: 化学反应工程
2楼2016-01-11 10:21:15
songsny
铜虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 112
- 帖子: 39
- 在线: 28.2小时
- 虫号: 2381059
- 注册: 2013-03-27
- 性别: MM
- 专业: 工程热物理相关交叉领域
3楼2016-01-28 11:22:31
4楼2016-04-05 20:27:08
10