24小时热门版块排行榜    

查看: 1184  |  回复: 2

天堂鱼520

银虫 (小有名气)

[求助] 关于组分传输传质迭代出现致命错误,求指点如何调试udf 已有2人参与

本人在做一个气液两相传质的模拟,编写的udf涉及到物料性质、adjust宏、define mass transfer。
模型类似一个圆柱,但是做的是二维的,也就是直管,液体沿壁面流下,壁面会加热,使液体蒸发,进而两相发生传质
设置是速度入口,自由出口,靠近液体的一侧壁面是加热壁面,另一侧是常温壁面。最初始的状态是直管内充满惰性气体,然后液体开始沿壁面流下,层流
求解的方程有Flow、volume fraction、energy和组分传输
为了使结果收敛,我先只求解了flow和volume fraction,每步都能收敛,2s后流动稳定保存,然后我试了以下尝试:
1)在流动的基础上,直接导入adjust宏,并在phase interaction下导入define mass transfer,开始同时求解Flow、volume fraction、energy和组分传输四个方程,出现致命错误;
2)在流动的基础上,没有导入adjust宏和define mass transfer,不勾选组分传输方程下的两个常用选项,求解Flow、volume fraction和组分传输三个方程,每步都收敛,又是在2s后保存;
3)在2)的基础上,直接导入adjust宏,并在phase interaction下导入define mass transfer,开始同时求解Flow、volume fraction、energy和组分传输四个方程,出现致命错误;
4)在2)的基础上,不导入adjust宏,只导入define mass transfer,开始同时求解Flow、volume fraction、energy和组分传输四个方程,出现致命错误;
5)在2)的基础上,不导入adjust宏,只导入define mass transfer,并将质量源项都改为常数,开始同时求解Flow、volume fraction、energy和组分传输四个方程,出现致命错误;
在代码中,我加入了一些printf("  KL=%10f\n",  KL);这样的输出语言,但是至始至终fluent窗口内都没有出现输出语句

求教各位大神,还有什么方法可以调试?拜托了,各位大神!指点下方向。。。。。
回复此楼
姑娘,没人替你坚强!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

转石千仞

金虫 (正式写手)

击破长空

【答案】应助回帖

妹子,我说过了,调试打印到控制台要用:
Message("the deep value is:%g\n", deep);
deep是参数,如果仅仅测试代码运行,可以用:
Message("the code position A\n";

Message("the code position B\n";
此外,还有:
if()
Error("here is wrong\n";
来使某此参数出问题时跳出计算报错。
调试的路很漫长,你还是慢慢来吧。
------------------------------
http://shop114586048.taobao.com/ ... 1997079397.2.A89lTp
打个广告,欢迎选购。
故善战人之势,如转圆石于千仞之山者,势也
2楼2015-03-20 00:07:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ndzhaolin

金虫 (小有名气)

【答案】应助回帖

应该是mass transfer代码有问题,可以把其改成常数验证
3楼2015-08-16 22:49:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 天堂鱼520 的主题更新
信息提示
请填处理意见