找回密码
 注册
查看: 1811|回复: 6

求各位大神指点迷津,小弟感激涕零啊

[复制链接]
发表于 2013-3-31 10:40:27 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?注册

x
各位大神,小弟要编一个关于进口mass flow 的UDF,其中的质量流量用进口平均密度乘以随时间变化的速度表示,一下是小弟编的udf代码,求大神们指点,万分感激啊
#include"udf.h"
real NV_VEC(A);
real avg_dest;
real sum_D_A=0.0;
real sum_A=0.0;
real T;

real unsteady_velocity;
Thread*thread_out;
face_t f;
Domain*domain;
DEFINE_ADJUST(adjust,d)
{
domain=Get_Domain(1);
thread_out=Lookup_Thread(domain,4);
begin_f_loop(f,thread_out)
   real T=CURRENT_TIME;
   {
     F_AREA(A,f,thread_out);
     sum_A+=NV_MAG(A);
     sum_D_A+=NV_MAG(A)*F_P(f,thread_out)/F_T(f,thread_out)*(1/287.1);
    }
end_f_loop(f,thread_out)
   unsteady_velocity=157.1*sin(3140.*T);
   avg_dest=sum_D_A/sum_A;
}
DEFINE_PROFILE(Mass_flux,t,i)
{
  begin_f_loop(f,t)
   {
     F_PROFILE(f,t,i)=avg_dest*unsteady_velocity;
   }
  end_f_loop(f,t)
}
 楼主| 发表于 2013-4-1 15:20:22 | 显示全部楼层
有人吗
 楼主| 发表于 2013-4-1 15:21:06 | 显示全部楼层
又能解答问题的大神吗
发表于 2013-4-2 09:42:08 | 显示全部楼层

回复 3# 李智慧19900624 的帖子

我是新手,请问楼主该程序中thread_out=Lookup_Thread(domain,4);中ID4是怎么得到的?我想要将指针指向具体一个区域中的某一个面上,ID值如何知道呢?
 楼主| 发表于 2013-4-2 11:19:17 | 显示全部楼层
这个好像是固定值,具体你把函数挂到哪里就是哪个区域的指针吧
发表于 2013-4-2 14:00:11 | 显示全部楼层

回复 5# 李智慧19900624 的帖子

谢谢楼主
 楼主| 发表于 2013-4-2 22:27:55 | 显示全部楼层
有能帮帮俺的大神出现吗
您需要登录后才可以回帖 登录 | 注册

本版积分规则

快速回复 返回顶部 返回列表