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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » udf 根据监测点的压力值设置入口的质量流量
61/1返回列表
查看: 2061  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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)
}
回复此楼

» 猜你喜欢
1楼 2013-05-12 20:40:20
已阅   回复此楼   关注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 ]
赞一下(1人) 回复此楼
率真随性的射手
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 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭