|
发表于 2013-5-21 11:53:11
|
显示全部楼层
回复 1# ltpreason 的帖子
#include "udf.h"
#define m_standard 8000.
#define m_interval 500.
real pressure;
DEFINE_ADJUST(mass_inlet_Get_P,d)
{
Thread *t;
cell_t c;
real xc[ND_ND];
thread_loop_c(t,d)
{
begin_c_loop(c,t)
{
C_CENTROID(xc,c,t);
if(sqrt(ND_SUM(pow(xc[0]-10.,2.), pow(xc[1]-10.,2.), pow(xc[2]-10.,2.)))<0.001)
pressure=C_P(c,t);
}
end_c_loop(c,t)
}
}
DEFINE_PROFILE(mass_inlet_as_P,t,i)
{
face_t f;
begin_f_loop(f,t)
{
/*根据监测点的压力值,来改变入口的质量流量*/
if(pressure<=350000.)
F_PROFILE(f,t,i) = m_standard+m_interval;
else if(pressure>350000.&&pressure<380000.)
F_PROFILE(f,t,i) = m_standard;
else
F_PROFILE(f,t,i) = m_standard-m_interval;
}
end_f_loop(f,t)
} |
|