| 查看: 6335 | 回复: 9 | |||||
[交流]
UDF中如何编写体积分数的梯度 已有6人参与
|
| 在气液两相流UDF中，需要编写气相和液相体积分数的梯度，应该用什么宏呢？ |
回复此楼» 收录本帖的淘帖专辑推荐
 传热传质 | 博士科研新路历程 |
» 猜你喜欢
听说今天filecode变了
已经有7人回复
-
基金中了
已经有3人回复
-
filecode
已经有10人回复
-
filecode
已经有12人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有5人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有7人回复
-
影响面上的因素
已经有11人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 编UDF将一个边界上的温度作为另一个边界条件的温度值
已经有13人回复
- 关于 Get_Domain()的问题
已经有5人回复
- DEFINE_ADJUST中的
已经有18人回复
- FLUENT中如何设置进料的浓度？
已经有7人回复
- Fluent欧拉模型中相的体积分数的设定
已经有16人回复
- udf获得多相流中某一项中某一组分的含量
已经有14人回复
- 如何在udf里面 判断界面的位置啊（vof模型）
已经有9人回复
- fluent 6.3 13.0版本中多相流的比较
已经有20人回复
- vof模型中获取第二项的体积分数
已经有6人回复
- vof模型算气液两相逆流 使用udf加动量源项 碰到问题 请高手解答
已经有14人回复
- 【求助】FLUENT UDF 气液两相流 汽蚀 汽化潜热(问题详见内容)
已经有21人回复
- VOF体积分数的梯度为什么为0呢？？
已经有10人回复
- 请教指定浓度的法向梯度为0的UDF如何写
已经有5人回复
- UDF中浓度梯度怎么写？
已经有3人回复
- udf中怎样获取流体流速？
已经有3人回复
- 介个UDF有啥问题
已经有10人回复
- 【求助】用vof模型计算出的界面，有木有udf语句可以找得到
已经有6人回复
- 【求助】udf读取混合物密度失败、及平均流速读取？
已经有8人回复
- 【求助】请教Fluent多相流流模型中的一个问题——相之间互相贯穿
已经有12人回复
chenjian5076
木虫 (小有名气)
老和山少侠
- 应助: 20 (小学生)
- 金币: 1635.5
- 红花: 10
- 帖子: 218
- 在线: 248.7小时
- 虫号: 1907370
- 注册: 2012-07-23
- 性别: GG
- 专业: 多相流热物理学
2楼2013-12-13 17:03:15
|
我编了这样一段程序。解释的时候出现：line 15: invalid type for binary expression: float * array 2 of float.的错误。line 15 以下是程序： #include "udf.h" #include "sg_mphase.h" #define T_SAT 319 #define LAT_HT 2.392e6 DEFINE_SOURCE(vap_src, cell, pri_th, dS, eqn) { Thread * mix_th, *sec_th; real m_dot_v,q; mix_th = THREAD_SUPER_THREAD(pri_th); sec_th = THREAD_SUB_THREAD(mix_th, 1); if(C_T(cell, mix_th)<=T_SAT) { q=-2*C_K_L(cell,pri_th)*C_T_G(cell,pri_th)*C_VOF(cell,pri_th); /*此行提示错误*/ m_dot_v=q/LAT_HT; dS[eqn] = -2*C_K_L(cell,pri_th)*C_T_G(cell,pri_th)/LAT_HT; } else { q=C_VOF_G(cell,sec_th); m_dot_v=q/LAT_HT; dS[eqn] = 0; } return m_dot_v; } DEFINE_SOURCE(liq_src, cell, sec_th, dS, eqn) { Thread *mix_th, *pri_th; real m_dot_l, q; mix_th = THREAD_SUPER_THREAD(sec_th); pri_th = THREAD_SUB_THREAD(mix_th,0); if(C_T(cell, mix_th)<=T_SAT) { q=2*C_K_L(cell,sec_th)*C_T_G(cell,sec_th)*C_VOF_G(cell,sec_th); m_dot_l=q/LAT_HT;; dS[eqn] =2*C_K_L(cell,pri_th)*C_T_G(cell,pri_th)/LAT_HT; } else { q=-2*C_K_L(cell,pri_th)*C_T_G(cell,pri_th)*C_VOF_G(cell,pri_th); m_dot_l=q/LAT_HT; dS[eqn]=0.; } return m_dot_l; } DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn) { Thread *pri_th, *sec_th; real m_dot, q; pri_th = THREAD_SUB_THREAD(mix_th, 0); sec_th = THREAD_SUB_THREAD(mix_th, 1); if(C_T(cell, mix_th)<=T_SAT) { q=2*C_K_L(cell,pri_th)*C_T_G(cell,pri_th)*C_VOF_G(cell,pri_th);; dS[eqn] = 0.; } else { q=-2*C_K_L(cell,sec_th)*C_T_G(cell,sec_th)*C_VOF_G(cell,sec_th); dS[eqn]=0.; } return q; } 我的udf编的不咋样。你帮我看看。我定义的是气液相变的源项。多谢啦 |
3楼2013-12-13 18:47:56
chenjian5076
木虫 (小有名气)
老和山少侠
- 应助: 20 (小学生)
- 金币: 1635.5
- 红花: 10
- 帖子: 218
- 在线: 248.7小时
- 虫号: 1907370
- 注册: 2012-07-23
- 性别: GG
- 专业: 多相流热物理学
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
这个函数在UDF帮助文件中叫做单元格函数，意思就是通过单元格的数值的差求得梯度，但是开始算的时候其实第一步是没有梯度的，所以这个地方会报错，具体的改进方法你可以看以下链接：http://www.cfd-online.com/Forums ... eulerian-model.html |
4楼2013-12-16 21:20:05
5楼2014-06-16 10:27:06
6楼2014-12-05 09:19:51
7楼2015-06-24 15:08:05
8楼2016-03-07 19:02:20
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
本帖内容被屏蔽 |
9楼2017-08-15 16:31:17
|
本帖内容被屏蔽 |
10楼2017-08-15 16:31:52
50