| 查看: 4534 | 回复: 4 | |||
[求助]
质量源项、动量源项
|
|
source = -mass/C_VOLUME(cell_gas,thread_gas)/CURRENT_TIMESTEP; /* 传质源项，单位 kg/m3-s */ mass具体是如何定义的？相减符号如何？ 还有，谁有动量源项，我已经在UDF里面写好质量源项，怎样加动量源项？在UDF里面 |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教设置源项的一个疑问
已经有6人回复
- FLUENT源项UDF非稳态迭代的一个疑问请教
已经有8人回复
- 时间函数的能量方程源项怎么表达
已经有3人回复
- UDF源项的一个疑问
已经有14人回复
- 求助气固两相流的边界条件问题
已经有8人回复
- 【讨论】Fluent边界条件及虚拟质量力
已经有6人回复
- 【求助】求助或讨论一下关于能量传递源项中各个焓值的含义
已经有3人回复
- 中科院高能物理所，东莞散列中子源项目的offer，这样的国家实验室的到底怎么样？
已经有9人回复
- 【求助】相变模拟
已经有6人回复
- 【求助】请问Fluent中如何简单地添加流动区域各坐标下的动量源项？
已经有16人回复
- 【讨论】请教一个气-液-固三相流模拟可能性的问题
已经有6人回复
2楼2012-06-21 22:35:29
3楼2012-06-27 09:50:57
|
DEFINE_SOURCE(gas_source, cell_gas, thread_gas, dS, eqn) /* gas相o2组分的传质源项 */ { real source; Domain *domain_liquid = Get_Domain(3); /* liquid相domain指针 */ Thread *thread_mix = THREAD_SUPER_THREAD(thread_gas); /* mixture的fluid-air区域指针 */ Thread *thread_liquid = Lookup_Thread(domain_liquid, 2); /* liquid相的fluid-water区域指针 */ cell_t cell_liquid; real m0_gas, m0_liquid, m1_gas, m1_liquid; if(C_UDMI(cell_gas, thread_mix, 0 ) == -1) /* 非相界面cell */ { source = 0.0; /* 传质源项为0 */ } else /* 相界面cell */ { 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)*kc*C_R(cell_liquid, thread_liquid)*C_VOLUME(cell_liquid, thread_liquid)/(kc*C_R(cell_liquid, thread_liquid)*C_VOLUME(cell_liquid, thread_liquid)+C_R(cell_gas, thread_gas)*C_VOLUME(cell_gas, thread_gas));/* 分配后液相中o2组分质量，单位 kg */ m1_gas = m1_liquid*C_R(cell_gas, thread_gas)*C_VOLUME(cell_gas, thread_gas)/(C_R(cell_liquid, thread_liquid)*C_VOLUME(cell_liquid, thread_liquid)*kc); /* 分配后气相中o2组分质量，单位 kg */ source = (m1_gas-m0_gas)/C_VOLUME(cell_gas, thread_gas)/(CURRENT_TIMESTEP+0.023); /* 传质源项，单位 kg/m3-s */ } dS[eqn] = 0.0; return source; } DEFINE_SOURCE(liquid_source, cell_liquid, thread_liquid, dS, eqn) /* liquid相o2组分的传质源项 */ { real source; Domain *domain_gas = Get_Domain(2); /* gas相domain指针 */ Thread *thread_mix = THREAD_SUPER_THREAD(thread_liquid); /* mixture的fluid-water区域指针 */ Thread *thread_gas = Lookup_Thread(domain_gas, 3); /* gas相的fluid-air区域指针 */ cell_t cell_gas; real m0_gas, m0_liquid, m1_gas, m1_liquid; if(C_UDMI(cell_liquid, thread_mix, 0 ) == -1) /* 非相界面cell */ { source = 0.0; /* 传质源项为0 */ } else /* 相界面cell */ { cell_gas = C_UDMI(cell_liquid, 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)*kc*C_R(cell_liquid, thread_liquid)*C_VOLUME(cell_liquid, thread_liquid)/(kc*C_R(cell_liquid, thread_liquid)*C_VOLUME(cell_liquid, thread_liquid)+C_R(cell_gas, thread_gas)*C_VOLUME(cell_gas, thread_gas));/* 分配后液相中o2组分质量，单位 kg */ m1_gas = m1_liquid*C_R(cell_gas, thread_gas)*C_VOLUME(cell_gas, thread_gas)/(C_R(cell_liquid, thread_liquid)*C_VOLUME(cell_liquid, thread_liquid)*kc); /* 分配后气相中o2组分质量，单位 kg */ source = (m1_liquid-m0_liquid)/C_VOLUME(cell_liquid, thread_liquid)/(CURRENT_TIMESTEP+0.023); /* 传质源项，单位 kg/m3-s */ } dS[eqn] = 0.0; return source; } DEFINE_SOURCE(gas_source1, cell_gas, thread_gas, dS, eqn) /* gas相o2组分的传质源项 */ { real source; Domain *domain_liquid = Get_Domain(3); /* liquid相domain指针 */ Thread *thread_mix = THREAD_SUPER_THREAD(thread_gas); /* mixture的fluid-air区域指针 */ Thread *thread_liquid = Lookup_Thread(domain_liquid, 2); /* liquid相的fluid-water区域指针 */ cell_t cell_liquid; real m0_gas, m0_liquid, m1_gas, m1_liquid; if(C_UDMI(cell_gas, thread_mix, 1 ) == -1) /* 非相界面cell */ { source = 0.0; /* 传质源项为0 */ } else /* 相界面cell */ { 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)*kc*C_R(cell_liquid, thread_liquid)*C_VOLUME(cell_liquid, thread_liquid)/(kc*C_R(cell_liquid, thread_liquid)*C_VOLUME(cell_liquid, thread_liquid)+C_R(cell_gas, thread_gas)*C_VOLUME(cell_gas, thread_gas));/* 分配后液相中o2组分质量，单位 kg */ m1_gas = m1_liquid*C_R(cell_gas, thread_gas)*C_VOLUME(cell_gas, thread_gas)/(C_R(cell_liquid, thread_liquid)*C_VOLUME(cell_liquid, thread_liquid)*kc); /* 分配后气相中o2组分质量，单位 kg */ source = (m1_gas-m0_gas)/C_VOLUME(cell_gas, thread_gas)/CURRENT_TIMESTEP; /* 传质源项，单位 kg/m3-s */ } dS[eqn] = 0.0; return source; } 这是我写的源项 -----但是我觉得我的周期性浓度变化没有写进来----求助啊 |
4楼2012-07-08 22:13:12
★
xiegangmai: 金币+1, 鼓励讨论交流！ 2012-07-08 22:35:08
xiegangmai: 金币+1, 鼓励讨论交流！ 2012-07-08 22:35:08
|
DEFINE_SOURCE(gas_source, cell_gas, thread_gas, dS, eqn) /* gas相o2组分的传质源项 */ { real source; Domain *domain_liquid = Get_Domain(3); /* liquid相domain指针 */ Thread *thread_mix = THREAD_SUPER_THREAD(thread_gas); /* mixture的fluid-air区域指针 */ Thread *thread_liquid = Lookup_Thread(domain_liquid, 2); /* liquid相的fluid-water区域指针 */ cell_t cell_liquid; real m0_gas, m0_liquid, m1_gas, m1_liquid; if(C_UDMI(cell_gas, thread_mix, 0 ) == -1) /* 非相界面cell */ { source = 0.0; /* 传质源项为0 */ } else /* 相界面cell */ { 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)*kc*C_R(cell_liquid, thread_liquid)*C_VOLUME(cell_liquid, thread_liquid)/(kc*C_R(cell_liquid, thread_liquid)*C_VOLUME(cell_liquid, thread_liquid)+C_R(cell_gas, thread_gas)*C_VOLUME(cell_gas, thread_gas));/* 分配后液相中o2组分质量，单位 kg */ m1_gas = m1_liquid*C_R(cell_gas, thread_gas)*C_VOLUME(cell_gas, thread_gas)/(C_R(cell_liquid, thread_liquid)*C_VOLUME(cell_liquid, thread_liquid)*kc); /* 分配后气相中o2组分质量，单位 kg */ source = (m1_gas-m0_gas)/C_VOLUME(cell_gas, thread_gas)/(CURRENT_TIMESTEP+0.023); /* 传质源项，单位 kg/m3-s */ } dS[eqn] = 0.0; return source; } 我这个因为是周期型传质---两边都传---然后右边的传质实际上是下一个长度周期的呀---所以我不知道怎样弄啊---我的udf 关于质量分数哪儿写的不对啊 |
5楼2012-07-08 22:28:44
25