吴子牛编著的《计算流体力学基本原理》中提到以下几种情况会有源项。
1.在非惯性坐标系中，N-S方程右端需添加反映惯性力影响的源项；
2.考虑重力影响时，需在动量方程中添加重力项；
3.考虑两相流时，在N-S方程右端需条件反映各种流体相互干扰的源项（包括质量交换，动量交换，能量交换）；
4.考虑两方程湍流模型时，湍流模型方程右端含有反映湍流生成与消灭的源项；
5.考虑化学反应时，需在组分浓度方程和能量方程添加由化学反应引起的质量变化与能量变化；
6.将高维问题化为低维问题时，也会出现源项。
流体力学中有点源，点汇的概念，我想可能这个与上述提到的应该有些相同吧。
陶文铨老师的书中对源项的描述是这样的：数值传热学书中所指的源项是广义量，它代表那些不能包括到控制方程的非稳态项、对流项与扩散项中所有其它各项之和。书中提到将源项局部线形化，S=Sc+SpTp,这种方法比将源项处理为常数会更合理。
查看了一个用UDF自定义的X方向动量源项
http://caohuali.combust.cn/blog/UploadFiles/2006-3/39619141.gif
经过变形处理后，给出的UDF程序是这样的：
/*******************************************************************/
/* UDF for specifying an x-momentum source term in a
spatially
*/
/* dependent porous
media
*/
/*******************************************************************/
＃i nclude "udf.h"
#define C2 100.0
DEFINE_SOURCE(xmom_source, c, t, dS, eqn)
{
real x[ND_ND];
real con, source;
C_CENTROID(x, c, t);
con = C2*0.5*C_R(c, t)*x[1];
source = -con*fabs(C_U(c, t))*C_U(c,
t);
dS[eqn] = -2.*con*fabs(C_U(c, t));
return source;
}
自己的体会就是：关键在于找到那些不能包括到控制方程的非稳态项、对流项与扩散项中所有项，写出源项的表达式。局部线形化的时候注意Sp要小于等于0.
|
加载中…