24小时热门版块排行榜

>论坛更新日志 (3395)
>虫友互识 (1225)
>休闲灌水 (247)
>导师招生 (187)
>文献求助 (115)
>论文投稿 (81)
>考博 (46)
>博后之家 (44)
>硕博家园 (38)
>基金申请 (35)
>考研 (34)
>招聘信息布告栏 (28)
>找工作 (28)
>教师之家 (24)
>论文道贺祈福 (18)
>公派出国 (18)

返回列表

【悬赏金币】回答本帖问题，作者Kellyshuo将赠送您 35 个金币

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

Kellyshuo

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 532
帖子: 6
在线: 17.3小时
虫号: 10175522
注册: 2018-09-06
专业: 工程热物理相关交叉领域

[求助] fluent UDF源项编译通过后，初始化闪退，求大佬指点一二已有1人参与

源项文件编译通过后，fluent初始化闪退代码见下：
#include "udf.h"
#include "math.h"
#define pre_exp 1.12e-05
#define active 29.16e+03
#define r 8.314
#define l 30e-06
#define mh2 2
define_source(h2_membrane_mass,c,t,ds,eqn)
{
real a[nd_nd];
real source;
real pr,pm,t;
real vol0,area0;
real vol1,area1;
face_t f;
real id = 3.0 ;
domain *d = get_domain(1);
thread *thread_mem_f = lookup_thread(d,id);

cell_t c0 = f_c0(f,thread_mem_f);
cell_t c1 = f_c1(f,thread_mem_f);
thread *tc0 = thread_t0(thread_mem_f);
thread *tc1 = thread_t1(thread_mem_f);
vol0 = c_volume(c0,tc0);
f_area(a,f,tc0);
area0 = nv_mag(a);
t = c_t(c0,tc0);
pr = c_p(c0,tc0) + 101325;
pm = c_p(c1,tc1) + 101325;
source = - area0*mh2*pre_exp*exp(-active/t)/(vol0*l)*(sqrt(pr)-sqrt(pm));
return source;

}

回复此楼