|
楼主 |
发表于 2009-9-18 11:12:21
|
显示全部楼层
我自己刚编了一段,恳请大家指正!
#include <udf.h>
DEFINE_PROFILE(air_velocity,thread,position)
{
face_t f;
real P_0=0;
real P_C;
real deltp=P_C-P_0;
real r;
real zeta;
begin_f_loop(f, thread)
{
P_C=F_P(f,t);
if (p_C<p_0)
vel=sqrt(2*deltp/(r*zeta));
else
vel=0;
return vel;
}
end_f_loop(f, thread)
} |
|