:( 紧急求助：如何通过能量源相给耦合壁面添加热量（附程序）

hf77105

任务说明：在燃气与壁面的耦合传热计算过程中，另外给壁面加入一定的热量,要求沿着壁面分段加入不同的热量。
已知：沿壁面 当x(X轴坐标)<=-35.58    加入的热流密度为q=Qmax=21.6e6 W/M*M;
               x=0         q=0.5*Qmax;
               x=59.6      q=0.1*Qmax;
               x=145.5     q=0.02*Qmax;
               x=159.45    壁面沿X轴终点坐标
               其它各处的热流密度插值求出
本人想通过对燃气的能量源项添加内热源来实现这一目的,以下是我写的UDF程序，但使用后壁面的温度和热流密度基本没有变化，这是为何，
请高手指点！
                 
&#35;include "udf.h"
DEFINE_SOURCE(cell_heat_source, cell, thread, dS, eqn)
{
  real xc[ND_ND];
  real source;
  C_CENTROID(xc, cell, thread);

      if(xc[0]<=-35.58)
      {source = 21.6e6;
        dS[eqn] = 0.;}
      else if(xc[0]<=0.)
{source = 21.6e6-0.5*21.6e6*(xc[0]+35.58)/(0.+35.58);
         dS[eqn]=0.;}
      else if(xc[0]<=59.6)
{source = 0.5*21.6e6 - 0.4*21.6e6*(xc[0]-0.)/(59.6-0.);
        dS[eqn] = 0;}
      else if(xc[0]<=145.5)
{source = 0.1*21.6e6 - 0.08*21.6e6*(xc[0]-59.6)/(145.5-59.6);
        dS[eqn] = 0;}
      else
        source = 0.02*21.6e6 - 0.02*21.6e6*(xc[0]-145.5)/(159.45-145.5);
        dS[eqn]=0;
   
   return source;
}

hf77105

高手帮忙顶一下好吗？

tmxk100

还不深明白你说的意思，可以在qq上进行讨论，qq：76396250

waitingforu

我怎么觉得是第二类边界条件阿？？define_profile

wearessweare

就不必用源项，边界条件的宏DEFINE PROFILE就可以修改热流密度，你写的能量源项相当于常数不利于迭代收敛，传热的过程也是从燃料到壁面，传热过程也反了

hf77105

如果用 define profile 修改后的热流密度,是通过壁面边界条件加进去的吗?