大家帮忙给看看这个入口边界udf

悠悠过客

系统入口有个止回阀，阀的开启由阀两侧的压力差控制，阀门只能向系统内开启。阀门上游的压力定位大气压，系统内的压力是随时间变化的。
我自己写了一段程序来定义入口的速度。大家帮忙给看有什么毛病。注：压力为表压值
#include "udf.h"
face_t f;
real r;
real zeta;
DEFINE_PROFILE(velocity,thread,position)
{
        real P_0=0; /*大气压*/
        real P_c; /*系统内的压力*/
        real delt_P=P_0-P_c;
        begin_f_loop(f,thread)
        {
                P_c=F_P(f,thread); /*对进气口截面压力取值*/
                if (P_0>P_c)
                        F_PROFILE(f,thread,position)=sqrt(2*delt_P/(r*zeta));/*给定一个压差与速度的关系式*/
                else
                        F_PROFILE(f,thread,position)=0;
        }
        end_f_loop(f,thread)       
}

xueyeguchou

把下面这三个放在DEFINE_PROFILE(velocity,thread,position)下面试试，你这几个参数在速度profile里面貌似没有啊
face_t f;
real r;
real zeta;