24小时热门版块排行榜    

查看: 5096  |  回复: 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本身就包含了气液界面啊。。。
回复此楼

» 收录本帖的淘帖专辑推荐

fluent气液两相udf

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liangjj626

新虫 (初入文坛)

你好,传质的udf编译好了加载到哪个边界上的?气泡patch上不知道这个udf加载到哪儿
5楼2017-03-17 01:14:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 1036985870 的主题更新
信息提示
请填处理意见