24小时热门版块排行榜    

查看: 973  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

王聪920704

新虫 (初入文坛)

[求助] UDF求助 已有2人参与

本人是UDF菜鸟,想请各位大神帮忙指导。我想提取近壁面第一层网格处的速度,之前从木虫里也找到一个类似的程序,现贴出如下:
#include "udf.h"
DEFINE_ADJUST(adjust,d)
{
face_t f;
cell_t c0;
real u_cell,n=0;
real u_face;
int ID=3;
Thread *thread_wall,*tc0;
thread_wall=Lookup_Thread(d,ID);
tc0=THREAD_T0(thread_wall);
begin_f_loop(f,thread_wall)
{
c0=F_C0(f,thread_wall);
C_UDMI(c0,tc0,0)=C_U(c0,tc0)/F_T(f,thread_wall);
}
end_f_loop(f,thread_wall)
}
我用了一个简单的圆柱试算了一下,发现速度仍为零,期待各位大神的指导,谢谢。
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gyctju

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
王聪920704: 金币+5, ★★★★★最佳答案 2016-11-27 20:36:21
目测是需要把thread_wall换成tc0

或者把“tc0= .。。。” 放到循环里面去

另外, 是此 “王聪”是我认识的那个“王聪”吗?
4楼2016-11-25 20:41:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

mycc

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
王聪920704: 金币+5, ★★★★★最佳答案 2016-11-27 20:36:30
应该没有问题,你算了多少步?另外/F_T(f,thread_wall);是为什么?
项目汇总:http://www.cfluid.com/forum.php?mod=viewthread&tid=114340&extra=
2楼2016-11-25 08:39:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

王聪920704

新虫 (初入文坛)

没说清楚,不好意思,这里面要用到近壁区第一层网格处的速度除以温度,所以有分母这一项。我只算了200步左右,因为很好收敛所以就没继续往下算,我在想是不是近壁区网格需要加密?
3楼2016-11-25 08:57:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

王聪920704

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by gyctju at 2016-11-25 20:41:37
目测是需要把thread_wall换成tc0

或者把“tc0= .。。。” 放到循环里面去

另外, 是此 “王聪”是我认识的那个“王聪”吗?

您好,刚尝试了一下,速度还是零,不知道问题出在哪里,您那边的话能帮我测试下这个程序吗?
5楼2016-11-27 21:33:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见