对流换热边界UDF问题

wenchang

我的问题是一个壁面是对流换热边界条件
其中对流换热系数用下式计算 h=17.6*(0.024-0.86*y)^0.32；
Y是沿壁面的高度。程序如下，但是编译不通，哪位高手指点一下，感激不尽。
#include "udf.h"
&#35;include <math.h>
DEFINE_PROFILE(Coefficient,t,i)
{
  real x[ND_ND];             /* this will hold the position vector */
  real y;
  face_t f;
  
  
  begin_f_loop(f,t)
    {
      F_CENTROID(x,f,t);
      y=x[1];
      F_PROFILE(f,t,i)=17.6*pow((0.024-0.86*y),0.32);
    }
  end_f_loop(f,t)
}

kevinxiao

错误提示是什么？

wenchang

iter     time/iter
!    1!    1
Error: Floating point error: invalid number
Error Object: ()
Error: FLUENT received a fatal signal (SEGMENTATION VIOLATION).
Error Object: ()

kevinxiao

编译不能通过还是不能计算？code本身好像没什么问题。
检查一下你的unit，确定0.024-0.86*y>0