| 查看: 922 | 回复: 1 | |||
[求助]
如何用UDF定义流体通过一个cell所用时间已有1人参与
|
|
double gasDissolvedRat(cell_t cell, Thread *cell_t) { double criNum, factor_index, J_ng; double vol_cell, mass_bubble, mass_max; double mass_flux, time_cell, pressure; double sol_CH4, mol_con; int n; face_t f; Thread *liq,*tf; criNum = gasMolNum_criticalcluster(cell, cell_t); factor_index = solubility_in_water_CH4(cell, cell_t); J_ng = factor_index*exp(-0.5*pow(criNum,2.0/3.0)); if (J_ng < 1.0) { J_ng = 0.0; } mass_bubble = mol_mass_CH4*criNum*J_ng/Avo_constant/1.0e12; vol_cell = C_VOLUME(cell, cell_t); mass_flux = 0; c_face_loop(cell,cell_t,n) /* loops over all faces of a cell */ { f = C_FACE(cell,cell_t,n); tf=C_FACE_THREAD(cell,cell_t,n); mass_flux = mass_flux + F_FLUX(f,tf); } liq =THREAD_SUB_THREAD(cell_t, 0); pressure = (C_P(cell, cell_t) + pressure_atm)/1.0e5; time_cell = vol_cell*C_R(cell, liq)/mass_flux ; sol_CH4 = -3.918E-06*pressure*pressure + 1.419E-03*pressure - 2.160E-04; mol_con = C_YI(cell,liq,0)*1000/mol_mass_CH4; if (sol_CH4<mol_con) { mass_max = (mol_con - sol_CH4)*mol_mass_CH4*C_R(cell,liq)/time_cell/1.0e3; }else { mass_max = 0.0; } if (mass_max>mass_bubble) { return mass_bubble; }else { return mass_max; } } 上述是关于含有CH4气体的水溶液在压力降低的情况下，释放CH4的函数。程序一旦运行，马上就会floating point exception。 1) time_cell 是流体流过一个网格所用的时间，想知道关于time_cell 的描述是否正确; 2) 大家是如何定义流体流过一个网格所用的时间; 3) 如何获取一个内部面上的速度向量； 4）程序还有什么其它错误。 谢谢谢。。。。  |
回复此楼» 收录本帖的淘帖专辑推荐
 udf相关 |
» 猜你喜欢
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有3人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有12人回复
-
博士读完未来一定会好吗
已经有27人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
怎么用udf读取在fluent中设置的surface-point？
已经有7人回复- dpm中用udf修改drag force，怎么获得颗粒所在cell流体的信息？
已经有3人回复
- fluent多孔介质孔隙率是位置的函数问题求助
已经有9人回复
- 流体物性随温度变化如何在fluent中设置
已经有11人回复
- 解释型UDF如何调用温度梯度
已经有9人回复
- UDF中的入门级问题
已经有13人回复
- fluent UDF 修改物性
已经有13人回复
- UDF中如何编写体积分数的梯度
已经有9人回复
- 用udf实现coupled边界固体区域中温度最值提取
已经有8人回复
- udf中如何编译壁面温度为X坐标的函数！
已经有8人回复
- fluent中流体出口在液面以下，怎么设置
已经有5人回复
- 利用UDF提取与壁面边界相邻的一层cell的方法
已经有6人回复
- fluent udf请教
已经有16人回复
- fluent 流体入口的热源加载问题
已经有5人回复
- UDF的编译通不过
已经有8人回复
- udf中在单元格的所有面上的循环
已经有4人回复
- 新手求救啊！Fluent非稳态计算 怎么样才能在每一个时间步上调用一次UDF啊？
已经有15人回复
- 在ANSYS14中的FLuent中，写关于流体比热物性的UDF
已经有14人回复
- 【求助】请问Fluent中如何简单地添加流动区域各坐标下的动量源项？
已经有16人回复
2楼2015-05-20 00:45:33