| 查看: 3263 | 回复: 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
» 猜你喜欢
假如你的研究生提出不合理要求
已经有10人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 气液两相流问题求助啊
已经有1人回复
- 关于组分传输传质迭代出现致命错误，求指点如何调试udf
已经有2人回复
- 粮仓模拟
已经有1人回复
- 关于C_YI_G（c,t，i）和F_AREA(A,f,t)的使用
已经有1人回复
- 关于udf中多相流的指针设置问题
已经有0人回复
- 气相在液相中的扩散
已经有0人回复
- 关于UDM的使用顺序问题
已经有9人回复
- 关于燃烧反应放出的热量的计算
已经有0人回复
- 关于编写两相流UDF中使用UDM的一些疑问——请大家多交流
已经有2人回复
- fluent UDF的编写 一个发热球体的热量随温度变化，与空气对流散热，请问怎么编写啊
已经有0人回复
mopsite
木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 8682.7
- 散金: 424
- 红花: 2
- 帖子: 2545
- 在线: 615小时
- 虫号: 172296
- 注册: 2006-01-19
- 专业: 传热传质学
6楼2015-09-24 13:44:54
9楼2015-10-29 14:25:20
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
BeObjective
银虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 1475.1
- 红花: 6
- 帖子: 733
- 在线: 82.5小时
- 虫号: 3311402
- 注册: 2014-07-07
- 性别: GG
- 专业: 传热传质学
5楼2015-09-23 10:45:49
mopsite
木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 8682.7
- 散金: 424
- 红花: 2
- 帖子: 2545
- 在线: 615小时
- 虫号: 172296
- 注册: 2006-01-19
- 专业: 传热传质学
7楼2015-10-10 12:59:26
BeObjective
银虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 1475.1
- 红花: 6
- 帖子: 733
- 在线: 82.5小时
- 虫号: 3311402
- 注册: 2014-07-07
- 性别: GG
- 专业: 传热传质学
|
http://muchong.com/bbs/viewthread.php?tid=9465108这是我另外发的一个帖子，这几天也在琢磨这个事，自己想法是问题可能出在两个地方，一个是梯度的处理，另外一个是边界源相位置的处理，我整个流体域分成了两个部分，靠近壁面那层网格被我处理成另外一个流体域，这个zone就成了源相的处理区域，今天在论文看到有人用了其他的办法。谢谢 有兴趣的话，你可以做一做 |
8楼2015-10-11 22:31:06
10楼2015-12-28 17:07:44