高斯热源分布的udf 编程

liyanzz

大家好，我是做焊接熔池温度场模拟的，所用软件为FLUENT。热源选用高斯分布加到壁面上，可我不会对热源进行编程，希望哪位高手能帮帮忙，指点一下，非常感谢！热源公式为
高斯函数分布可表示为 q(r)=（0.65×14×150）÷（π×0.012）×exp(-(x2＋z2)÷0.012)

liyanzz

我仿着帮助文档，编了一个定点热源的例子，我对C语言了解比较少，也不知编的对不对，请各位高手帮我指点一下，万分感谢！
#include "udf.h"
DEFINE_PROFILE(immovable_heat,t,i)
{
real x[ND_ND];
real y,z;
face_t f;
begin_f_loop(f,t)
  {
  F_CENTROID(x,f,t);
  z=x[1];
  y=x[0];
  F_PROFILE(f,t,i)=0.65*16*110/(3.1415926*0.003*0.003)*exp(-(y*y+z*z)/(0.003*0.003);
  }
end_f_loop(f,t)
}