24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2170  |  回复: 5

ltpreason

新虫 (初入文坛)

[求助] udf 根据监测点的压力值设置入口的质量流量

初学udf,请大家帮忙看下哪有问题,或者如果这样思路不对,应该如何做,已被折磨n天,跪求大神指教。
#include "udf.h"
#include "mem.h"
#define m_standard 8000.
#define m_interval 500.
DEFINE_ADJUST(mass_inlet,d)
{
Thread *t1,*t2;
cell_t c;
face_t f;
real xc[ND_ND];
int i;
real p=0.;
thread_loop_c(t1,d)
{
/*坐标为(10,10)的点为监测点,循环找到一些和它接近的点,
把这些点的压力值赋给p,并认为最后一次赋的值即为监测点的值*/
begin_c_loop(c,t1)
{
C_CENTROID(xc,c,t1);
if(sqrt(ND_SUM(pow(xc[0]-10.,2.),
pow(xc[1]-10.,2.)))<0.001)
p=C_P(c,t1);
}
end_c_loop(c,t1)
}
/*根据监测点的压力值,来改变入口的质量流量*/
begin_f_loop(f,t2)
{
if(p<=350000.)
F_PROFILE(f,t2,i) = m_standard+m_interval;
else if(p>350000.&&p<380000.)
F_PROFILE(f,t2,i) = m_standard;
else
F_PROFILE(f,t2,i) = m_standard-m_interval;
}
end_f_loop(f,t2)
}
回复此楼

» 猜你喜欢

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

shuai_1314

捐助贵宾 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
ben_ladeng: 金币+1, 谢谢指点 2013-05-13 19:37:51
你的thread t1和t2没有给出明确的定义,根本不知道是指向哪个区域,看下帮助文档里的data access吧

[ 发自手机版 http://muchong.com/3g ]
率真随性的射手
2楼2013-05-13 07:57:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ltpreason

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by shuai_1314 at 2013-05-13 07:57:05
你的thread t1和t2没有给出明确的定义,根本不知道是指向哪个区域,看下帮助文档里的data access吧

谢谢大神,能再给些具体的修改意见吗?我今天自己看了下那的帮助,还是一头雾水。。
3楼2013-05-13 19:21:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shuai_1314

捐助贵宾 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
xiegangmai: 金币+2, 谢谢参与! 2013-05-15 11:22:00
ltpreason: 金币+5, 有帮助 2013-05-16 20:37:12
引用回帖:
3楼: Originally posted by ltpreason at 2013-05-13 19:21:11
谢谢大神,能再给些具体的修改意见吗?我今天自己看了下那的帮助,还是一头雾水。。...

Thread *thread_name = Lookup_Thread(domain,zone_ID)

用这个命令,zone_ID就是你的各种边界和流动区域对应的那个编号,比如你的fluid编号是2,你的入口边界编号是5

Thread *t1 = Lookup_Thread(d,2);
Thread *t2 = Lookup_Thread(d,5);

应该这样定义,多看看帮助文档搞清楚domain和thread
率真随性的射手
4楼2013-05-14 23:09:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ltpreason

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by shuai_1314 at 2013-05-14 23:09:57
Thread *thread_name = Lookup_Thread(domain,zone_ID)

用这个命令,zone_ID就是你的各种边界和流动区域对应的那个编号,比如你的fluid编号是2,你的入口边界编号是5

Thread *t1 = Lookup_Thread(d,2);
Th ...

非常感谢您!
我加了这两句进去,调用时还是报错。
大神再帮忙看看。。还有点什么问题呢?先谢谢了!
5楼2013-05-16 20:36:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ltpreason

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by shuai_1314 at 2013-05-14 23:09:57
Thread *thread_name = Lookup_Thread(domain,zone_ID)

用这个命令,zone_ID就是你的各种边界和流动区域对应的那个编号,比如你的fluid编号是2,你的入口边界编号是5

Thread *t1 = Lookup_Thread(d,2);
Th ...

非常感谢您!
我加了这两句,可是调用时还是报错。
能再帮忙看看吗?先谢谢您了!
6楼2013-05-16 20:38:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ltpreason 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿安大生物学07初试322、本科二本、调剂求助 +10 李多米lee. 2026-04-12 11/550 2026-04-12 22:58 by yuyin1233
[考研] 电子信息270求调剂 +18 terminal469 2026-04-07 18/900 2026-04-12 16:23 by ajpv风雷
[考研] 一志愿郑州大学 22408 305分求调剂 +5 安小满zzz 2026-04-08 5/250 2026-04-12 00:41 by 蓝云思雨
[考研] 22专硕求调剂 +6 haoyun上岸 2026-04-11 8/400 2026-04-11 23:21 by labixiaoqiao
[考研] 一志愿211,0703化学305分求调剂 +26 严西西戏 2026-04-06 33/1650 2026-04-11 23:01 by 314126402
[考研] 297求调剂 +9 Kwgyz 2026-04-09 9/450 2026-04-11 10:09 by zhq0425
[考研] 22408 366分,本科211,一志愿西工大 +4 Rubt 2026-04-09 4/200 2026-04-10 19:51 by chemisry
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 一志愿沪9,326生物学求相关专业调剂 +4 刘墨墨 2026-04-09 4/200 2026-04-10 12:07 by pengliang8036
[考研] 考研二轮调剂 +8 故人?? 2026-04-09 8/400 2026-04-10 09:44 by 青梅duoduo
[考研] 材料专硕初试分332一志愿西北工业大学, +12 故人?? 2026-04-09 12/600 2026-04-09 18:34 by Ccclqqq
[考研] 化学工程与技术专业一志愿哈工程 291分B区 国家级大创负责人 有一作论文 +13 Emmy~ 2026-04-09 13/650 2026-04-09 14:47 by only周
[考研] 材料307分求大佬组收留 +17 Hll胡 2026-04-07 17/850 2026-04-09 10:53 by liuhuiying09
[考研] 328求调剂 +17 lftmya 2026-04-07 18/900 2026-04-09 08:05 by 5268321
[考研] 机械专硕273请求调剂 +6 庚申壬申 2026-04-07 6/300 2026-04-08 22:41 by bljnqdcc
[考研] 考研求调剂 +4 雯??? 2026-04-08 4/200 2026-04-08 21:44 by 土木硕士招生
[考研] 电子信息346 +4 zuoshaodian 2026-04-08 4/200 2026-04-08 11:54 by zzucheup
[考研] 307求调剂 +14 超级伊昂大王 2026-04-06 14/700 2026-04-08 07:03 by 无际的草原
[考研] 277、学硕,求调剂 数一104, +11 瓶子PZ 2026-04-07 12/600 2026-04-07 23:30 by 一只好果子?
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 09:25 by 小黑不怕难
信息提示
请填处理意见