压力边界条件是该边界出口速度的函数，请问该函数哪里错了？

易阳

大家好，我用FLUENT模拟三维低速压气机，出口压力边界条件是该边界出口速度的函数，编好UDF如下，先静态计算1000步，然后加载UDF，算了500步正常，然后100步，进口速度没变，出口速度，压力，流量直线上升，超过计算极限，请大家帮忙看看，是UDF的问题吗？
#include "udf.h"

DEFINE_PROFILE(unsteady_pressure, thread, position)
{
  face_t f;
  real k;
  


  begin_f_loop(f, thread)
  {
    k = F_U(f,thread);
    F_PROFILE(f, thread, position) = 101325.0 + k*k;
  }
  end_f_loop(f, thread)

}

另，我开始想直接用UDF算，但初始化时用进口速度，提示出现错误，我才不得不先计算静态，问一下，这个错误是什么原因？

易阳

刚看了一下论坛帖子，，发现了一个问题，，就是我是按照FLUENT帮助文档的例子写的，，那个例子是根据位置来给定压力，，所以用的是，，position，而我这里应该是速度，，，知道哪里错了，，但不知道怎么改啊，，加油，，

易阳

很显然，直接把，位置POSITION改成，速度VELOCITY是不行的，，请帮个忙啊，，

易阳

怎么没个人回复呢？
看了一天，，发现直接用UDF无法初始化的问题，，应该是没有使用能量方程的缘故，，，使用了能量方程之后可以计算，，，至于结果，，好象没有明显的问题，，
而那个先用稳态值做初始化，然后 再用UDF算到了，500步就直线上升，，有资料介绍，，好象是因为进出口区域太短了的缘故，，有非物理反射，，，奇怪，，我已经延长了双倍弦长啊，，过会试验一下延长5倍，，，，，