应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » udf编写多孔介质气固换热系数遇到问题，求助~~~
101/1返回列表
查看: 3116  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

高歌旅客Les

银虫 (小有名气)

[求助] udf编写多孔介质气固换热系数遇到问题，求助~~~ 已有4人参与

在做多孔介质局部非热平衡模型，设置系数h时用udf编程，udf可以编译，但是初始化的时候总是报错，求助大神，到底怎么回事呢？？？？
#include "udf.h"
#include "metric.h"
#include "sg.h"
#include "mem.h"
#include "math.h"

#define Dp 0.00006
#define zeta 0.6
DEFINE_PROFILE(Velocity_inlet,thread,position)
{
        real u;
        real t;
        face_t f;
        begin_f_loop(f,thread)
        {
                t=RP_Get_Real("flow-time";
                F_PROFILE(f,thread,position)=5*sin(188.5*t);
        }
        end_f_loop(f,thread)
}
DEFINE_PROFILE(Pressure_outlet,thread,position)
{
        real p;
        real t;
        face_t f;
        begin_f_loop(f,thread)
        {
                t=RP_Get_Real("flow-time";
                F_PROFILE(f,thread,position)=4500000+1170000*sin(188.5*t);
        }
        end_f_loop(f,thread)
}
/*3.
This module calculates the Heat transfer coefficient h(loc) for both the solid
and the fluid phase based on the particle Reynolds Number Rep.

*/

double heat_t_coeff(cell_t c,Thread *t)
{
        //FILE *t;
        //*HTC = (kf*(2+1.1Pr^0.333*Re^0.6))/dp;*/
        double Pr,Re,HTC,U,dv,Rep,Area,Ux;
        //Interstial area between the solid and the Fluid
        Area=4.0*(1-zeta)/Dp;
        //void diameter
        dv=(4.0*zeta)/Area;
        //Average velocity
        U=sqrt(pow(C_U(c,t),2.)+pow(C_V(c,t),2.));
        // modified Reynolds Number based on the void diameter
        Re=(U*dv*C_R(c,t))/C_MU_EFF(c,t);
        //Particle Reynolds Number based on Diamter of the particle(Dp).
//        Rep = (U*Dp*C_R(c,t))/C_MU_EFF(c,t);
        Rep=(U*Dp*C_R(c,t))/C_MU_L(c,t);

        //Fluid Prandtl Number
//        Pr = (C_MU_EFF(c,t)*C_CP(c,t))/C_K_L(c,t);
//        Pr = (C_MU_L(c,t)*C_CP(c,t))/therm_conduct(C_T(c,t)) ;
        Pr=(C_MU_L(c,t)*C_CP(c,t))/C_K_L(c,t) ;


//        C_UDMI(c,t,1)=therm_conduct(C_T(c,t));  
       
        /*Use this part of the code for E2 model of HTC*/
//    printf("hello";
//        HTC = (C_K_L(c,t)*(2.0 +(1.1*pow(Pr,0.333)*pow(Rep,0.6))))/Dp;
//        HTC = (therm_conduct(C_T(c,t)) *(2.0 +(1.1*pow(Pr,0.333)*pow(Rep,0.6))))/Dp;
//        HTC = (C_K_L(c,t)*(2.0 +(1.1*pow(Pr,0.333)*pow(Rep,0.6))))/Dp;
        HTC=(C_K_L(c,t)*(2.0 +(1.1*pow(Pr,0.333)*pow(Rep,0.6))))/Dp;
        /*Use this part of code for E2 model of HTC

if(Rep < 100)
                HTC=0.004*(C_K_EFF(c,t)/dv)*pow(Pr,0.3333)*pow(Re,1.35);
        else
                HTC =0.0156*(C_K_EFF(c,t)/dv)*pow(Pr,0.3333)*pow(Re,1.04);
        */
        return HTC;
}

初始化的时候报错如下：
Error:  received a fatal signal (Segmentation fault).
Error Object: #f
chip-exec: heat_t_coeff: wrong return type: void udf function expected
Error:  received a fatal signal (Segmentation fault).

求助大神！！！！！！@wuming524
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2016-09-06 21:31:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

【答案】应助回帖

感谢参与，应助指数 +1
http://www.cfd-online.com/Forums ... entation-fault.html
The segmentation fault means that you try to access a variable, which is not saved.
赞一下(1人) 回复此楼
2楼2016-09-06 22:13:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mycc

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
高歌旅客Les: 金币+5, ★★★很有帮助 2016-09-07 10:31:06
你这个宏都没有用到嘛，double heat_t_coeff(cell_t c,Thread *t)
初始化错误的话，可以尝试一个个加载宏，看哪个宏加载之后出现错误，那么问题就出现在那个宏之中
赞一下 回复此楼
项目汇总：http://www.cfluid.com/forum.php?mod=viewthread&amp;tid=114340&amp;extra=
3楼2016-09-07 08:09:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dai26520

专家顾问 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
heat_t_coeff(cell_t c,Thread *t)我觉得这个地方直接用一个DEFINE-PROFILE宏就行，直接定义气固对流换热系数，然后在GUI里面设置对流换热系数时加载DEFINE-PROFILE宏就行了。另外你这个可能还需要一个DEFINE-ADJUST宏。
赞一下(2人) 回复此楼
个性只是一个代号，做好自己就行
4楼2016-09-07 08:35:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

高歌旅客Les

银虫 (小有名气)
引用回帖:
3楼: Originally posted by mycc at 2016-09-07 08:09:55
你这个宏都没有用到嘛，double heat_t_coeff(cell_t c,Thread *t)
初始化错误的话，可以尝试一个个加载宏，看哪个宏加载之后出现错误，那么问题就出现在那个宏之中

应该是最后那个对流换热系数的问题，因为之前的速度和压力我都加载过，没什么问题~~~
然后我是新手，想请教一下，这个宏我都没用到是啥意思呢？？？
赞一下 回复此楼
5楼2016-09-07 10:18:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mycc

专家顾问 (正式写手)
引用回帖:
5楼: Originally posted by 高歌旅客Les at 2016-09-07 10:18:42
应该是最后那个对流换热系数的问题，因为之前的速度和压力我都加载过，没什么问题~~~
然后我是新手，想请教一下，这个宏我都没用到是啥意思呢？？？...

UDF 编译好之后，需要加载到相应的位置。 比如DEFINE-PROFILE宏加载到边界条件上
你这个heat_t_coeff是一个自定义的宏，你无法加载到fluent中，需要通过fluent已经存在define宏来调用。
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

高歌旅客Les
项目汇总：http://www.cfluid.com/forum.php?mod=viewthread&amp;tid=114340&amp;extra=
6楼2016-09-07 11:54:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

高歌旅客Les

银虫 (小有名气)
送红花一朵
引用回帖:
6楼: Originally posted by mycc at 2016-09-07 11:54:00
UDF 编译好之后，需要加载到相应的位置。 比如DEFINE-PROFILE宏加载到边界条件上
你这个heat_t_coeff是一个自定义的宏，你无法加载到fluent中，需要通过fluent已经存在define宏来调用。...

感谢您的回答，我懂了！
所以我要用heat_t_coeff的话，要先用fluent自带的define来调用对吧，那么我选用哪个呢？DEFINE_property？
麻烦您了，感谢！
赞一下 回复此楼
7楼2016-09-07 22:01:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dai26520

专家顾问 (著名写手)

【答案】应助回帖

引用回帖:
7楼: Originally posted by 高歌旅客Les at 2016-09-07 22:01:56
感谢您的回答，我懂了！
所以我要用heat_t_coeff的话，要先用fluent自带的define来调用对吧，那么我选用哪个呢？DEFINE_property？
麻烦您了，感谢！...

DEFINE-PROFILE宏
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

高歌旅客Les
个性只是一个代号，做好自己就行
8楼2016-09-08 08:38:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

高歌旅客Les

银虫 (小有名气)
送红花一朵
引用回帖:
8楼: Originally posted by dai26520 at 2016-09-08 08:38:10
DEFINE-PROFILE宏...

感谢您的帮助~~~
我有个问题，profile不是定义边界条件的宏么？为什么可以用profile呢？？？
麻烦啦~
赞一下 回复此楼
9楼2016-09-08 09:25:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奔向大海

新虫 (初入文坛)

【答案】应助回帖

您好！我想请问您一下，为什么设置4个HTC呢？而且其中三个表达式完全一样？您能解释一下吗？灰常谢谢！！！
赞一下 回复此楼
10楼2018-06-27 13:53:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 高歌旅客Les 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 13:54 by GkCVyODlx4mR
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 jCd0dEvKHShX 2026-08-29 6/300 2026-08-30 13:54 by GkCVyODlx4mR
[基金申请] 有没有仍没收到信息的 +6 德尚中行 2026-08-27 7/350 2026-08-30 12:53 by Bxndjjd
[基金申请] 国社科又开始会评了，不知道这次命运如何 +3 雨打竹帘 2026-08-30 4/200 2026-08-30 12:35 by 白大多亮
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:39 by ZPa0EcMwuECS
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 06:57 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 06:22 by ZPa0EcMwuECS
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 06:00 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 13/650 2026-08-30 05:38 by ZPa0EcMwuECS
[公派出国] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:45 by ZPa0EcMwuECS
[基金申请] 为什么到现在没收到通知？ +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
关闭
请填处理意见
关闭