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

sasa

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

请不吝赐教

cora

一个例子，自己照着改吧

1. 
   
2. /*******************************************************************/
   
3. /* UDF for specifying an x-momentum source term in a spatially     */
   
4. /* dependent porous media                                          */
   
5. /*******************************************************************/
   
6. 
   
7. ＃i nclude "udf.h"
   
8. 
   
9. #define C2 100.0
   
10. 
    
11. DEFINE_SOURCE(xmom_source, c, t, dS, eqn)
    
12. {
    
13.   real x[ND_ND];
    
14.   real con, source;
    
15. 
    
16.   C_CENTROID(x, c, t);
    
17.   con = C2*0.5*C_R(c, t)*x[1];
    
18. 
    
19.   source = -con*fabs(C_U(c, t))*C_U(c, t);
    
20.   dS[eqn] = -2.*con*fabs(C_U(c, t));
    
21. 
    
22.   return source;
    
23. }

复制代码

sasa

谢谢回复

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

chbo76

可以只给source 不给dS[eqn]

sasa

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

chbo76

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

cora

原帖由 sasa 于 2008-9-25 11:41 发表

                               
登录/注册后可看大图

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

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

sasa

谢谢谢谢～

zhuimeng163

也影响结果吧好像

chengzinj

现在遇到同样的问题，也是动量源项不能对速度求导。。。请问有人做过吗？