24小时热门版块排行榜    

查看: 809  |  回复: 2

bwj0319

新虫 (小有名气)

[交流] fluen中的UDF程序问题 已有1人参与

模拟往复式压缩机出口管道需要添加一个脉动的压力,脉动形式如图所示,表达式也在下面,周期是0.3s,针对这个问题,由于刚上手fluent不久,对UDF功能也不是很熟,外加C语言也很菜,于是仿照别人的程序编了一段,大家看能实现图示和表达式所需要的脉动压力吗?我的fluent开始算了,看监测点的曲线 好像开始就有脉动了,估计有问题啊 想听听大牛们的意见啊
程序:

/* unsteady.c */
/* UDF for specifying a transient pressure profile boundary condition */
#include "udf.h"
DEFINE_PROFILE(unsteady_pressure,thread,position)
{
face_t f;
begin_f_loop(f,thread)
{
real t=RP_Get_Real("flow-time";

int k;
for(k=0;k<=9;k=k+1)
{
if(t>=0.3*k&&t<0.225+0.3*k)
F_PROFILE(f,thread,position)=310000000;
else
F_PROFILE(f,thread,position)=310000000*(1+0.06*sin(2*(20.93*t-4.71)));
}
}
end_f_loop(f,thread)
}

模拟过程中设置监测点,结果监测点曲线如下附图,监测数据:"Convergence history of Static Pressure on point-1 (in SI units)"
"Flow Time" "Vertex Average Static Pressure"
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991117e+08
0.0060000003 3.0991117e+08
0.0060000003 3.099112e+08
0.0060000003 3.099112e+08
0.0060000003 3.0991123e+08
0.0060000003 3.099113e+08
0.0060000003 3.0991142e+08
0.0060000003 3.0991155e+08
0.0060000003 3.0991168e+08
0.0060000003 3.0991187e+08
0.0070000003 3.0991187e+08
0.0070000003 3.0991216e+08
0.0070000003 3.0991248e+08
0.0070000003 3.0991293e+08
0.0070000003 3.099135e+08
0.0070000003 3.0991421e+08
0.0070000003 3.0991498e+08
0.0070000003 3.0991587e+08
0.0070000003 3.099169e+08
0.0070000003 3.0991805e+08
0.0070000003 3.0991936e+08
0.0070000003 3.099208e+08
0.0070000003 3.0992243e+08
0.0070000003 3.0992416e+08
0.0070000003 3.0992598e+08
0.0070000003 3.0992797e+08
0.0070000003 3.0993011e+08
0.0070000003 3.0993235e+08
0.0070000003 3.0993466e+08
0.0070000003 3.0993709e+08
0.0070000003 3.0993962e+08
0.0080000004 3.0993962e+08
0.0080000004 3.0994237e+08
0.0080000004 3.0994541e+08
0.0080000004 3.099487e+08
0.0080000004 3.0995213e+08
0.0080000004 3.0995571e+08
0.0080000004 3.0995946e+08
0.0080000004 3.0996333e+08
0.0080000004 3.0996723e+08
0.0080000004 3.099712e+08
0.0080000004 3.099752e+08
0.0080000004 3.0997923e+08
0.0080000004 3.0998323e+08
0.0080000004 3.099873e+08
0.0080000004 3.099913e+08
0.0080000004 3.099953e+08
0.0080000004 3.0999933e+08
0.0080000004 3.1000326e+08
0.0080000004 3.1000717e+08
0.0080000004 3.1001104e+08
0.0080000004 3.1001491e+08
0.0090000004 3.1001491e+08
0.0090000004 3.1001882e+08
0.0090000004 3.1002288e+08
0.0090000004 3.1002694e+08
0.0090000004 3.1003114e+08
0.0090000004 3.1003533e+08
0.0090000004 3.1003958e+08
0.0090000004 3.1004384e+08
0.0090000004 3.1004813e+08
0.0090000004 3.1005238e+08
0.0090000004 3.1005664e+08
0.0090000004 3.1006086e+08
0.0090000004 3.1006506e+08
0.0090000004 3.1006918e+08
0.0090000004 3.1007328e+08
0.0090000004 3.1007731e+08
0.0090000004 3.1008128e+08
0.0090000004 3.1008525e+08
0.0090000004 3.1008909e+08
0.0090000004 3.100929e+08
0.0090000004 3.1009664e+08
0.01 3.1009664e+08
0.01 3.1010042e+08
0.01 3.1010413e+08
0.01 3.1010784e+08
0.01 3.1011146e+08
0.01 3.101151e+08
0.01 3.1011869e+08
0.01 3.1012224e+08
0.01 3.1012573e+08
0.01 3.1012918e+08
0.01 3.1013264e+08
0.01 3.1013606e+08
0.01 3.1013946e+08
0.01 3.1014278e+08
0.01 3.1014614e+08
0.01 3.1014944e+08
0.01 3.101527e+08
0.01 3.1015594e+08
0.01 3.1015917e+08
0.01 3.1016234e+08
0.01 3.101655e+08
0.011000001 3.101655e+08
0.011000001 3.1016867e+08
。。。。。。
------------------------------------------------------
结果应该不对,没能实现需要的压力曲线啊 ,有懂UDF程序的大牛们 指导下俺吧 ,愁死了。。。。


压力曲线



压力表达式

监测点压力曲线



[ Last edited by bwj0319 on 2012-7-29 at 21:46 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)

★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
xiegangmai: 金币+1, 谢谢参与! 2012-07-30 21:57:53
for(k=0;k<=9;k=k+1)
{
if(t>=0.3*k&&t<0.225+0.3*k)

k=0,if里面t就是从0开始了
2楼2012-07-29 19:26:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bwj0319

新虫 (小有名气)

引用回帖:
2楼: Originally posted by kylafree at 2012-07-29 19:26:27
for(k=0;k<=9;k=k+1)
{
if(t>=0.3*k&&t<0.225+0.3*k)

k=0,if里面t就是从0开始了

是啊  我就是要从0开始啊,这里是个周期函数,周期为0.3s,每周期里又分为两段了
3楼2012-07-29 22:06:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bwj0319 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 梦想 +5 myrtle 2026-08-26 8/400 2026-09-02 02:42 by Leogzhya
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 为什么国自然不能直接公布 +5 bjdxyxy 2026-08-26 5/250 2026-09-01 17:54 by 小伟大博士
[基金申请] 学科评审组评审是指会评吗? +4 瞬息宇宙 2026-08-31 4/200 2026-09-01 14:58 by jiaoxg
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
信息提示
请填处理意见