(急)請教各位高手關於使用udf加入熱源的問題(附程序)

blackcup

[这个贴子最后由blackcup在 2006/01/23 06:09pm 第 2 次编辑]

我的問題是要在一個區域內加入一隨距離 y 遞減之熱源
source term為C1*sqrt(C2*decay)  其中decay=exp(-EPSON*S*y)
程序如下:
#include "udf.h"
#define DHP 1.66e-5
#define KP 0.6
#define M 3.0
#define FI 0.1
#define EPSON 200.0
#define S 0.15
#define I0 20000.0
#define KT 0.5
DEFINE_SOURCE(laser_source, cell, thread, dS, eqn)
{
  real xc[ND_ND], source;
  real C1,C2;
  real decay,y;
  begin_c_loop(cell, thread) {
  C_CENTROID(xc, cell, thread);
  C1=DHP*KP*M;
  C2=FI*EPSON*S*I0/KT;
  y = xc[1];
  decay=exp(-EPSON*S*y);
  source=C1*sqrt(C2*decay);
  dS[eqn]=0.0;
  }end_c_loop(cell,thread);
  return source;
  }

interprete沒有問題,但帶入後跑出來的結果似乎不正確,
請各位高手可以幫我看看這樣寫是否符合我所要定義的source term

feich

去掉loop循环试试看

popolar

begin_c_loop(cell, thread)
end_c_loop(cell,thread);
应该没有必要使用这个循环吧