24小时热门版块排行榜

返回列表

nba24

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 413.7
红花: 4
帖子: 130
在线: 59.5小时
虫号: 1570800
注册: 2012-01-09
专业: 化学反应工程

[求助] UDF ---传质

高手好，
      液液两相泰勒流传质，大家看看我的各个源项设置对不对啊！
      不要打酱油啊
1、我只设置了X方向动量源项，没有设置Y方向动量源项
2、质量源项，我既设置了mass  ，又设置了C4H8O2（乙酸乙酯）的质量源项。
不知道我这样设置对不对啊，各位评论一下
动量源项UDF
cell_liquid = C_UDMI(cell_gas, thread_mix, 0); /* 获取相界面另一侧的cell编号 */
/* 下面就可通过thread_mix, cell_gas, thread_gas, cell_liquid, thread_liquid五个变量获取需要的数据计算动量源项 */
m0_gas = C_YI(cell_gas, thread_gas, 0)*C_R(cell_gas, thread_gas)*C_VOLUME(cell_gas, thread_gas); /* 气相中o2组分质量，单位 kg */
m0_liquid = C_YI(cell_liquid, thread_liquid, 0)*C_R(cell_liquid, thread_liquid)*C_VOLUME(cell_liquid, thread_liquid); /* 液相中o2组分质量，单位 kg */
m1_liquid = (m0_gas+m0_liquid)/(1.0+1/kc); /* 分配后液相中o2组分质量，单位 kg */
m1_gas = m1_liquid/kc; /* 分配后气相中o2组分质量，单位 kg */
source = (m1_gas-m0_gas)*C_U(cell_gas, thread_mix)/C_VOLUME(cell_gas, thread_gas)/CURRENT_TIMESTEP; /* 动量源项，单位 N/m3 */

质量源项UDF：
cell_liquid = C_UDMI(cell_gas, thread_mix, 1); /* 获取相界面另一侧的cell编号 */
/* 下面就可通过thread_mix, cell_gas, thread_gas, cell_liquid, thread_liquid五个变量获取需要的数据计算传质源项 */
m0_gas = C_YI(cell_gas, thread_gas, 0)*C_R(cell_gas, thread_gas)*C_VOLUME(cell_gas, thread_gas); /* 气相中o2组分质量，单位 kg */
m0_liquid = C_YI(cell_liquid, thread_liquid, 0)*C_R(cell_liquid, thread_liquid)*C_VOLUME(cell_liquid, thread_liquid); /* 液相中o2组分质量，单位 kg */
m1_liquid = (m0_gas+m0_liquid)/(1.0+1/kc); /* 分配后液相中o2组分质量，单位 kg */
m1_gas = m1_liquid/kc; /* 分配后气相中o2组分质量，单位 kg */
source = (m1_gas-m0_gas)/C_VOLUME(cell_gas, thread_gas)/CURRENT_TIMESTEP; /* 传质源项，单位 kg/m3-s */

大家对我做的点评一下啊，急用

左油相，右水相；初始酯质量分数，油相0.1，水相0

动量源项设置

质量源项设置

回复此楼