| 查看: 3431 | 回复: 12 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
BeObjective银虫 (正式写手)
|
[求助]
关于UDF源项的编写遇到了一下问题，有大神做过类似的东西的望指点 已有1人参与
|
||
|
首先介绍一下我的问题吧，做混合气体（空气和水蒸气）的冷凝，网格如下图所示： 二维网格，上部是进口，下面是出口，左边是恒温壁面，右边是绝热壁面，蒸汽会在左边的壁面上进行冷凝，使用的是species模型，冷凝的过程是通过UDF添加质量、能量源项的方式来实现的。公式如下图： 这几个是原论文中的质量源项公式： 我自己的想法是这样的：质量源项出现的区域肯定是在靠近壁面的地方，而在这里的模型中忽略了液膜的厚度，那么靠近壁面的那一层网格就是源项需要处理的区域，于是通过FLUENT的Adapt>mark，标记那一层网格，然后separate zone将整个流体域分成了两个zone，这样udf源项只在靠近壁面的那一层网格中去做。为了简化问题，我先只做了质量源项的部分，能量源项暂时不考虑，写完程序如下： #include "udf.h" DEFINE_SOURCE(wall_msource,c,t,dS,eqn) { real source,temp,pres,dens; real m_f,diff; int i=0; temp=C_T(c,t); pres=C_P(c,t); dens=C_R(c,t); m_f=C_YI(c,t,i); C_UDMI(c,t,0)=C_YI_G(c,t,i)[0]; diff=(2.775e-6+4.479e-8*temp+1.656e-10*temp*temp)*pres/100; source=dens*diff*C_UDMI(c,t,0)/(1-m_f); return source; } 关于这个程序，有一点说明，论文中公式后面几个是计算质量分数的，但是，udf能够存储水蒸汽质量分数，为什么论文中还多此一举去求质量分数？这是我的一点疑问，于是我这里就直接获得求解器的质量分数，并且获得了质量分数的梯度，这样问题就来了，在这个区域只有一层网格，获得质量分数没问题，但是这个质量分数梯度的获得有没有问题，按理说质量分数的梯度至少是两个网格单元做的差分，一层网格获得梯度貌似存在问题。这个程序写完我导入FLUENT进行计算，就出现问题了，在残差收敛后，我监控了进出口的质量流量，发现这个源项产生的作用几乎为0，进出口完全持平，于是我为了排除是不是源项太小产生的误差直接把质量源项部分改写成-0.4，得到非常显著的进出口流量差异。以下两个图就是UDF源项和constant-0.4源项所产生的水蒸气质量分数分布云图： 图一是UDF的云图，而且放大了很多： 图二是常数源项的云图： 两者对比非常的明显。图一貌似没有任何变化，图二的分布很明显。这两个计算的CASE都已经上传上去了。 最近也是做这个东东做的有点麻木了，其中也做了很多交叉对比，没有实质性的进展，看看有没有有经验的大神一眼看穿问题出在哪里？只要确实能点到问题所在，不会吝啬金币，谢了！  2.jpg  3.jpg  4.jpg  5.jpg  无标题.jpg  1.jpg  2.jpg  3.jpg  4.jpg  5.jpg |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : fluent111.cas
- 附件 2 : fluent112.cas
2015-09-21 22:26:13, 694.51 K
2015-09-21 22:26:14, 694.61 K
» 猜你喜欢
又一批高校组建人工智能学院 师资行吗 不是骗人吗
已经有7人回复
-
国自然面上和省基金B类撒花
已经有18人回复
-
有没有学校收留
已经有3人回复
-
312求调剂
已经有3人回复
-
华师大读博
已经有5人回复
-
急需审稿人！！！
已经有3人回复
-
申博/考博
已经有8人回复
-
295分求调剂
已经有6人回复
-
085600材料与化工调剂
已经有6人回复
-
有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业
已经有7人回复
|
本帖内容被屏蔽 |
11楼2016-09-05 13:08:35
BeObjective
银虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 1475.1
- 红花: 6
- 帖子: 733
- 在线: 82.5小时
- 虫号: 3311402
- 注册: 2014-07-07
- 性别: GG
- 专业: 传热传质学
2楼2015-09-22 12:04:47
BeObjective
银虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 1475.1
- 红花: 6
- 帖子: 733
- 在线: 82.5小时
- 虫号: 3311402
- 注册: 2014-07-07
- 性别: GG
- 专业: 传热传质学
3楼2015-09-22 12:07:32
4楼2015-09-22 15:23:31
5