[求助]请教高手，udf，谢谢！

wulong2477

在三维的情况下，出口边界条件中，我是想先得到单元的重心坐标，然后再在计算出口的压力，udf这样写的：
     #include "udf.h"
     DEFINE_PROFILE(unsteady_pressure, thread, position)
     {
      real x[ND_ND];
      real y;
      cell_t c;
      begin_c_loop(c, thread)
     {
      C_CENTROID(x,c,thread);
      y = x[1];
      C_P(c, thread, position) =24.01*1000-9.8*1000*y;
      }
      end_c_loop(c, thread)
      }
编译时没有出错，但是在对整个区域初始化时出错了，是不是我的写错了？请各位大虾给与指正，谢谢了！

topshen

第一，出口边界应该是 面 吧，不是格子操作。
   
   
第二，第11行不对：注意C_P(c, thread, position)好像是不可以这样的，这个函数返回值，不能赋值。可以改为F_PROFILE(f，thread,index)=....这样应该ok了。
就是改成
#include "udf.h"
    DEFINE_PROFILE(unsteady_pressure, thread, position)
    {
     real x[ND_ND];
     real y;
     face_t c;
     begin_f_loop(c, thread)
    {
     C_CENTROID(x,c,thread);
     y = x[1];
     F_PROFILE(c, thread, position) =24.01*1000-9.8*1000*y;
     }
     end_f_loop(c, thread)
     }
请试试看！good luck!

topshen

sorry!错了一行！汉！
#include "udf.h"
   DEFINE_PROFILE(unsteady_pressure, thread, position)
   {
    real x[ND_ND];
    real y;
    face_t c;
    begin_f_loop(c, thread)
   {
    F_CENTROID(x,c,thread);
    y = x[1];
    F_PROFILE(c, thread, position) =24.01*1000-9.8*1000*y;
    }
    end_f_loop(c, thread)
    }

wulong2477

谢谢关注，不过还是有点错误，你看一下这是我计算二维时用的出口边界，
     #include "udf.h"
     DEFINE_PROFILE(unsteady_pressure, thread, position)
     {
      real x[ND_ND];
      real y;
      face_t f;
      begin_f_loop(f, thread)
     {
      F_CENTROID(x,f,thread);
      y = x[1];
      F_PROFILE(f, thread, position) =24.01*1000-9.8*1000*y;
      }
      end_f_loop(f, thread)
      }
是不是可以用在三维里？这样对不对啊?

topshen

   
   应该可以用在3维的模型里面，维度ND_ND宏从系统传值过来，二维是2，三维是3。边界的类型都是 面 类型操作，而不管实际是线还是面。
good luck!

topshen

对整个区域（不是边界）初始化时出错的问题：
如果你是要对某个流场区域（局部块或全部）进行初始化操作，第一，简单一点的条件可以在做做网格时候先定义一下，然后在计算的时候用pacth.第二，如果比较复杂，可以用custom fuction来pacth,如果还复杂一点，就编初始化udf吧。

wulong2477

谢谢关注，正在调试中！