24小时热门版块排行榜

返回列表

summer同学

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 118.8
散金: 30
帖子: 49
在线: 23.8小时
虫号: 4026814
注册: 2015-08-17
性别: MM
专业: 多相流热物理学

[求助] udf求助已有1人参与

#include "udf.h"
#define MW_W 0.018
#define pi 3.14159
#define a_c 0.04
#define R  8.314
#define N 185.371
#define L 22.6
DEFINE_PROFILE(hflux,thread,index)
{
      face_t f;
real sumh=0.0;
      real A[ND_ND];
      Thread *threadh;
begin_f_loop(f,thread)
{
      real Ti;
   Ti=F_T(f,thread);
   F_AREA(A,f,thread);
   sumh+=(2*a_c/(2-a_c))*sqrt(MW_W/(2*pi*R))*NV_MAG(A)*(611.2*(pow(2.718,(17.62*(Ti-273.15)/

(Ti-30.03))))/sqrt(Ti)-N)*L;
}
end_f_loop(f,thread)
      F_PROFILE(f,thread,index)=sumh;
}
定义了一个heat flux的UDF，与界面温度有关，编译没错，但是怀疑没起到实际作用，求各位大神指点指点。

回复此楼

» 猜你喜欢

拟解决的关键科学问题还要不要写已经有8人回复
最失望的一年已经有12人回复
存款400万可以在学校里躺平吗已经有29人回复
求推荐英文EI期刊已经有5人回复
请教限项目规定已经有4人回复
国自然申请面上模板最新2026版出了吗？已经有20人回复
26申博已经有3人回复
基金委咋了？2026年的指南还没有出来？已经有10人回复
基金申报已经有6人回复
疑惑？已经有5人回复

1楼 2017-09-02 11:21:42

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

tumu_tu

金虫 (小有名气)

应助: 5 (幼儿园)
金币: 2650
散金: 37
红花: 2
帖子: 167
在线: 60.6小时
虫号: 3769411
注册: 2015-03-27
性别: GG
专业: 建筑环境与结构工程

F_PROFILE(f,thread,index)=sumh;没有在面上循环，这样只在一个网格面上定义了值，其他没定义。

赞一下

回复此楼

2楼2017-09-02 22:05:35

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

summer同学

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 118.8
散金: 30
帖子: 49
在线: 23.8小时
虫号: 4026814
注册: 2015-08-17
性别: MM
专业: 多相流热物理学

引用回帖:

2楼: Originally posted by tumu_tu at 2017-09-02 22:05:35
F_PROFILE(f,thread,index)=sumh;没有在面上循环，这样只在一个网格面上定义了值，其他没定义。

感谢你的回复，F_PROFILE(f,thread,index)=sumh目的在于返回循环里面求和的最终值sumh，你的意思是把这一语句放入循环里面吗？

赞一下

回复此楼

3楼2017-09-03 08:34:13

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

tumu_tu

金虫 (小有名气)

应助: 5 (幼儿园)
金币: 2650
散金: 37
红花: 2
帖子: 167
在线: 60.6小时
虫号: 3769411
注册: 2015-03-27
性别: GG
专业: 建筑环境与结构工程

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
summer同学: 金币+20, ★★★很有帮助 2017-09-04 08:36:29

引用回帖:

3楼: Originally posted by summer同学 at 2017-09-03 08:34:13
感谢你的回复，F_PROFILE(f,thread,index)=sumh目的在于返回循环里面求和的最终值sumh，你的意思是把这一语句放入循环里面吗？...

再写一个begin...end_f_loop(f,thread)，把F_PROFILE放进去。

CODE:

#include "udf.h"

#define MW_W 0.018

#define pi 3.14159

#define a_c 0.04

#define R 8.314

#define N 185.371

#define L 22.6

DEFINE_PROFILE(hflux, thread, index)

{

        face_t f;

        real sumh = 0.0;

        real A[ND_ND];

        Thread *threadh;

        begin_f_loop(f, thread)

        {

                real Ti;

                Ti = F_T(f, thread);

                F_AREA(A, f, thread);

                sumh += (2 * a_c / (2 - a_c))*sqrt(MW_W / (2 * pi*R))*NV_MAG(A)*(611.2*(pow(2.718, (17.62*(Ti - 273.15) / (Ti - 30.03)))) / sqrt(Ti) - N)*L;

        }

        end_f_loop(f, thread)

        begin_f_loop(f, thread)

        {

                F_PROFILE(f, thread, index) = sumh;

        }

        end_f_loop(f, thread)

}