24小时热门版块排行榜    

查看: 187  |  回复: 0

tsduyikang

新虫 (小有名气)

[求助] 对一个边界面上所有节点上的压力值进行操作,若出现负压,对负压力值置0。

不确定是对流体域的cell压力进行操作还是对边界面上的压力进行操作?若程序有错误,请指正,谢谢

1、对流体域cell压力进行操作程序
#include "udf.h"
DEFINE_ADJUST(adjust_pressure_cell,domain)
{
#if RP_NODE
Thread *c_thread;
cell_t c;
domain=Get_Domain(1);
c_thread=Lookup_Thread(domain,3);
begin_c_loop(c,t)
{
        if(C_P(c,t)<0)
        {
                C_P(c,t)=0;
        }
        else
        {
C_P(c,t)=0;
        }

}
end_c_loop
#endif
}
2、对边界面压力进行操作程序
DEFINE_ADJUST(adjust_pressure_face,domain)
{
#if RP_NODE
Thread *c_thread;
cell_t c;
domain=Get_Domain(1);
c_thread=Lookup_Thread(domain,4);
begin_c_loop(c,t)
{
        if(F_P(c,t)<0)
        {
                F_P(c,t)=0;
        }
        else
        {
F_P(c,t)=0;
        }

}
end_c_loop
#endif
}
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tsduyikang 的主题更新
信息提示
请填处理意见