应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 关于UDF源项的编写遇到了一下问题，有大神做过类似的东西的望指点
131/2返回列表12下一页
查看: 3202  |  回复: 12
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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都已经上传上去了。
最近也是做这个东东做的有点麻木了，其中也做了很多交叉对比，没有实质性的进展，看看有没有有经验的大神一眼看穿问题出在哪里？只要确实能点到问题所在，不会吝啬金币，谢了！关于UDF源项的编写遇到了一下问题，有大神做过类似的东西的望指点
2.jpg


关于UDF源项的编写遇到了一下问题，有大神做过类似的东西的望指点-1
3.jpg


关于UDF源项的编写遇到了一下问题，有大神做过类似的东西的望指点-2
4.jpg


关于UDF源项的编写遇到了一下问题，有大神做过类似的东西的望指点-3
5.jpg


关于UDF源项的编写遇到了一下问题，有大神做过类似的东西的望指点-4
无标题.jpg


关于UDF源项的编写遇到了一下问题，有大神做过类似的东西的望指点-5
1.jpg


关于UDF源项的编写遇到了一下问题，有大神做过类似的东西的望指点-6
2.jpg


关于UDF源项的编写遇到了一下问题，有大神做过类似的东西的望指点-7
3.jpg


关于UDF源项的编写遇到了一下问题，有大神做过类似的东西的望指点-8
4.jpg


关于UDF源项的编写遇到了一下问题，有大神做过类似的东西的望指点-9
5.jpg
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : fluent111.cas
    • 2015-09-21 22:26:13, 694.51 K
  • 附件 2 : fluent112.cas
    • 2015-09-21 22:26:14, 694.61 K

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
liveinday-tightcompartments.
1楼2015-09-21 22:38:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

mopsite

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
BeObjective: 金币+10, ★★★很有帮助, 不好意思，回复慢了，确实是饱和湿度的求解公式，谢谢！UDF你用的多吗？最近老出现segementation fault，而且有时候计算完以后吧，保存case，重新导入没法进行后处理，unavailable。 2015-09-30 07:43:03
那个计算的 是饱和湿度吧，不是直接提取的。。。
赞一下(2人) 回复此楼
6楼2015-09-24 13:44:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

唐小虎

新虫 (初入文坛)
我在做多孔介质，其中多孔介质的界面要模拟蒸发，用到UDF，也卡主好久了。。。希望以后多多交流，你问题解决了也可以交流下经验
赞一下(1人) 回复此楼
9楼2015-10-29 14:25:20
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

BeObjective

银虫 (正式写手)
没人看嘛 ，自己顶一下了
回复此楼
liveinday-tightcompartments.
2楼2015-09-22 12:04:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

BeObjective

银虫 (正式写手)
会的不要藏着嘛

发自小木虫Android客户端
回复此楼
liveinday-tightcompartments.
3楼2015-09-22 12:07:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

830706ljx

新虫 (初入文坛)
不会也帮顶下贴吧，最近也在做udf，现在连编译也编译不成

发自小木虫Android客户端
赞一下 回复此楼
4楼2015-09-22 15:23:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

BeObjective

银虫 (正式写手)
卡住了，真蛋疼，一个人做东西有时候真的很麻烦，复杂一点就容易出错，谢谢顶一顶

发自小木虫Android客户端
赞一下 回复此楼
liveinday-tightcompartments.
5楼2015-09-23 10:45:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mopsite

木虫 (著名写手)
将你的case和udf发来看看，是什么情况导致segementation fault的
赞一下 回复此楼
7楼2015-10-10 12:59:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

BeObjective

银虫 (正式写手)
引用回帖:
7楼: Originally posted by mopsite at 2015-10-10 12:59:26
将你的case和udf发来看看，是什么情况导致segementation fault的

http://muchong.com/bbs/viewthread.php?tid=9465108这是我另外发的一个帖子，这几天也在琢磨这个事，自己想法是问题可能出在两个地方，一个是梯度的处理，另外一个是边界源相位置的处理，我整个流体域分成了两个部分，靠近壁面那层网格被我处理成另外一个流体域，这个zone就成了源相的处理区域，今天在论文看到有人用了其他的办法。谢谢 有兴趣的话，你可以做一做
赞一下 回复此楼
liveinday-tightcompartments.
8楼2015-10-11 22:31:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

下雨天呼呼

新虫 (初入文坛)
内容已删除
回复此楼
10楼2015-12-28 17:07:44
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 BeObjective 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭