24小时热门版块排行榜

返回列表

木风10086

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 117
红花: 5
帖子: 53
在线: 40.6小时
虫号: 1266670
注册: 2011-04-15
性别: GG
专业: 流体力学

[交流] 近壁区域的源项（壁面温度）

区域如图所示，白色部分是固体，下面是气体。其实想模拟的是固体受热后发生燃烧。思路是取与固体接触的那一层气体单元，判断那一侧气体单元，如果单元内的温度、压力以及与它接触的固体单元内的温度达到一定值后就从这层气体单元内加质量源项和能量源项等。源项的代码（示意代码）如下：
其中的C_T(cell,thread和C_P(cell,thread)都是判断气体层单元的温度和压力，而F_T(f,t)本来是想着无法获得旁边固体单元的温度，所以用这层边界面的温度代替，但是计算后发现这个F_T(f,t)变量根本没起到作用。虽然源项能计算，但是明明表面温度还很低源项就加进去了。请问应该如果处理？或者能不能直接得到与气体相邻的固体单元内的温度？
DEFINE_SOURCE(c_m, cell, thread, dS, eqn)
{
real source;
if (C_T(cell,thread) >= 320.&&C_P(cell,thread) >= 120000.&&F_T(f,t)>=1000.)
{
   source =1
   dS[eqn]=0;
   }
else
{
   source =0;
   dS[eqn] = 0;
   }
return source;
}

DEFINE_SOURCE(c_e, cell, thread, dS, eqn)
{
real source;
if (C_T(cell,thread) >= 320.&&C_P(cell,thread) >= 120000.&&F_T(f,t)>=1000.)
{
   source =1;
   dS[eqn] = 0;
   }
else
{
   source = 0;
   dS[eqn] = 0;
   }
return source;
}

123321.png

回复此楼