24小时热门版块排行榜    

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

菲之小部落

木虫 (小有名气)

[求助] fluent中嵌入UDF后,初始化后,出现分段性错误是怎么回事? 已有1人参与

我的UDF如下:
#include "udf.h"

#define o2water_coeff 2.08e-9;
#define o2air_coeff 0.178e-4;
#define h_coeff 0.04375;
#define air_density 1.3;
#define water_density 998;
#define water_viscosity 0.001;
#define air_viscosity 0.000017;
#define surface_tension 0.0719404;

DEFINE_PROPERTY(cell2_diameter,c,t)
{
float d;
float mu_dia;
float eps=C_D(c,t);
d=C_VOF(c,t);
mu_dia=5.12*pow((0.0719404/(1.3*998)),0.6)*pow((0.000017/0.001),0.25)*pow(d,0.65)*pow(eps,-0.4);
return mu_dia;
}
DEFINE_SOURCE(ca_source,c,t,dS,eqn)
{
float con;
float d;
float a;
float b;
float y;
float eps=C_D(c,t);
d=C_VOF(c,t);
a=1.3*2*sqrt(2.08e-9/3.14);
b=pow(eps*998/0.001,0.25);
y=6*d/(5.12*pow((0.0719404/(1.3*998)),0.6)*pow((0.000017/0.001),0.25)*pow(d,0.65)*pow(eps,-0.4));
con=a*b*y;
dS[eqn]=con*h_coeff;
return con;
}
DEFINE_SOURCE(cw_source,c,t,dS,eqn)
{
float con1;
float a;
float b;
float y;
float d;
float eps=C_D(c,t);
d=C_VOF(c,t);
a=998*2*sqrt(2.08e-9/3.14);
b=pow(C_D(c,t)*998/0.001,0.25);
y=6*d/(5.12*pow((0.0719404/(1.3*998)),0.6)*pow((0.000017/0.001),0.25)*pow(d,0.65)*pow(eps,-0.4));
con1=a*b*y;
dS[eqn]=-con1;
return con1;
}
DEFINE_DIFFUSIVITY(ca_diffusivity,c,t,i)
{
return (0.178e-4*1.3+C_MU_T(c,t)/0.7);
}
DEFINE_DIFFUSIVITY(cw_diffusivity,c,t,i)
{
return(2.08e-9*998+C_MU_T(c,t)/0.7);
}
设为VOF,三相,空气,水,氧气,其中水为主相,空气和氧气均为第二相,设置了氧气在空气中的传质系数,氧气在水中的传质系数,就是UDF中设置的,在边界条件中,气体进口处,设置空气的体积系数为1,然后以气体进口为初始化条件进行初始化,结果出现如下error:received a fatal signal(segmentation fault)
求助于各位大神,哪里出了错??
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

死天

专家顾问 (著名写手)

【答案】应助回帖

氧气本来就是和空气混溶啊
9楼2014-09-07 13:57:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

ziyun899

禁虫 (初入文坛)

本帖内容被屏蔽

2楼2014-07-14 21:19:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲之小部落

木虫 (小有名气)

引用回帖:
2楼: Originally posted by ziyun899 at 2014-07-14 21:19:44
楼主,我也遇到了这个问题,请问你的问题解决了么

没有,一直找不到解决方法,你呢?
tomorrowisanotherday
3楼2014-07-15 11:01:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ziyun899

禁虫 (初入文坛)

本帖内容被屏蔽

4楼2014-07-15 15:59:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 03:55 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 03:40 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 11/550 2026-08-23 02:55 by OEbVnUOu01ol
[教师之家] 跳槽后在研项目怎么办? +3 简单化xn 2026-08-22 7/350 2026-08-22 19:46 by 简单化xn
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见