sasa 发表于 2008-9-24 09:54

求助：动量源项dS[eqn]怎么写？body force

当source=(C_U(cell, thread)*C_V_G(cell, thread)+C_V(cell, thread)*C_V_G(cell, thread))
dS怎么对C_V(cell, thread)求导？
是否可以只给source 不给dS
body force 相可以作为动量源项添加吗？是否有专门的body force命令？

请不吝赐教

cora 发表于 2008-9-24 14:28

一个例子，自己照着改吧


/*******************************************************************/
/* 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;
real con, source;

C_CENTROID(x, c, t);
con = C2*0.5*C_R(c, t)*x;

source = -con*fabs(C_U(c, t))*C_U(c, t);
dS = -2.*con*fabs(C_U(c, t));

return source;
}

sasa 发表于 2008-9-24 15:24

谢谢回复

这个例子的source可以直接对x方向的速度求导得到dS【eqn】
(C_U(cell, thread)*C_U_G(cell, thread)+C_V(cell, thread)*C_U_G(cell, thread))项中含有速度对x的求导和速度对y的求导，要得到dS【eqn】如何对这个求导进行x方向的速度的求导呢？

chbo76 发表于 2008-9-25 10:58

可以只给source 不给dS

sasa 发表于 2008-9-25 11:41

对于连续性方程的源项，也可以只给source，不给dS吗？不影响计算结果？

chbo76 发表于 2008-9-27 10:53

我没试过。
你可以试试，应该就是收敛慢点吧。

cora 发表于 2008-9-28 07:25

原帖由 sasa 于 2008-9-25 11:41 发表 http://www.chinavib.com/forum/images/common/back.gif
对于连续性方程的源项，也可以只给source，不给dS吗？不影响计算结果？

只影响收敛性，不影响结果

sasa 发表于 2008-9-28 09:14

谢谢谢谢～

zhuimeng163 发表于 2008-11-29 08:59

也影响结果吧好像

chengzinj 发表于 2011-3-23 15:34

现在遇到同样的问题，也是动量源项不能对速度求导。。。请问有人做过吗？
页: [1]
查看完整版本: 求助：动量源项dS[eqn]怎么写？body force