关于udf小程序
#include"udf.h"
#include"mem.h"
#define p_pretotal 350000
#define pcr 185000
#define R 287.06
#define k 1.4
#define T 300
DEFINE_ADJUST(on_demamd_avg,d)
{
real p_sum = 0;
real p;
real v = 0;
real r = 0,tem;
real V;
real xc[ND_ND];
Thread *t;
face_t f;
Domain *d;
d = Get_Domain(1); /*获取液态区域的网格 */
t=Lookup_Thread(d,12);/* 找到出口截面的线指针(ID在fluent中的phase找的) */
/* 单元循环 */
begin_f_loop(f,t)
{
F_CENTROID(xc,f,t); /* 将网格坐标存入数组xc中 */
p_sum = F_P(f, t);
if(p_sum<=p_pretotal)
{
if(p_sum>=pcr)
{
F_U(f,t)=sqrt((2*k*R*T/(k-1))*pow((1-(p_sum/p_pretotal)),(k-1)/k));
}
else
F_U(f,t)=sqrt(k*R*T*pow((p_sum/p_pretotal),(k-1)/k));
}
else
F_U(f,t)=0;
}
end_f_loop(f,t);
}
通过获取压力 比较压力得到相应的速度公式 求高手指点