| 查看: 695 | 回复: 2 | |||
[交流]
近壁区域的源项（壁面温度）
|
|
区域如图所示，白色部分是固体，下面是气体。其实想模拟的是固体受热后发生燃烧。思路是取与固体接触的那一层气体单元，判断那一侧气体单元，如果单元内的温度、压力以及与它接触的固体单元内的温度达到一定值后就从这层气体单元内加质量源项和能量源项等。源项的代码（示意代码）如下： 其中的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 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有9人回复
-
时间戳今天，20号变了
已经有5人回复
-
2026国自然函评费到账
已经有9人回复
-
放榜前的不淡定
已经有14人回复
-
建议基金发布提前给出明确的时间点
已经有13人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有20人回复
-
科研孤儿太难了
已经有18人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
-
filecode，4个jtjc了
已经有14人回复
2楼2020-10-27 18:10:07
3楼2020-10-28 23:13:00
40