24小时热门版块排行榜    

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

白白的大萝卜

新虫 (小有名气)

[求助] UFD调试问题!! 已有1人参与

在Fluent 中加入源项的UDF  一初始化就开始报错  Error: received a fatal signal(Segmentation fault).,UDF是照着别人的文献写得,没有书写错误,现在不知道问题出在哪里,本人对UDF和C也是刚入门的新手,在这里 希望好心人能帮忙调试一下,有偿也可以的!!!救救孩子吧  万分感谢!
#include "udf.h"
/*Define which user-defined scalars to use.*/
enum
{
   p,ads,d_ads,q_ads
};
DEFINE_ADJUST(p_absolute,domain)
{
        Thread*t;
        cell_t c;
        real p_operating;
        real p_static;
        /*solve the absolute pressure.*/
        thread_loop_c(t,domain)
        {
        if(NULL!=THREAD_STORAGE(t,SV_UDS_I(p)))
          {
                begin_c_loop(c,t)
                {
                p_static=C_P(c,t);
                p_operating=RP_Get_Real("operating-pressure";
                C_UDSI(c,t,p)=p_static+p_operating;
                }
            end_c_loop(c,t)
      }
        }
}
DEFINE_ADJUST(adsorption,domain)
{
        Thread*t;
        cell_t c;
        real n_o;
        real E_t;
        real x_p;
        real A;
        real x_b;
        real n;
        /*solve the adsorption equation.*/
        thread_loop_c(t,domain)
        {
                if(NULL!=THREAD_STORAGE(t,SV_UDS_I(ads))&&NULL!=T_STORAGE_R_NV(t,SV_UDSI_M1(ads)))
         {
          begin_c_loop(c,t)
            {
        real tem=C_T(c,t);
        n_o=71.6;
        E_t=3080+18.9*tem;
        x_p=1.47E9/C_UDSI(c,t,p);
        A=8.31429*tem*log(x_p);
        x_b=-A*A/(E_t*E_t);
        C_UDSI(c,t,ads)=n_o*exp(x_b);
        }
        end_c_loop(c,t)
     }
    }
}

DEFINE_ADJUST(q_adsorption,domain)
{
        Thread*t;
        cell_t c;
        /*solve the absolute pressure.*/
        thread_loop_c(t,domain)
        {
                if(NULL!=THREAD_STORAGE(t,SV_UDS_I(q_ads))&&NULL!=T_STORAGE_R_NV(t,SV_UDSI_M1(q_ads)))
                {
                        begin_c_loop(c,t)
                        {
                real physical_dt;
                physical_dt=RP_Get_Real("physical-time-step";
                C_UDSI(c,t,q_ads)=C_UDSI_M1(c,t,q_ads)+0.001*physical_dt*(C_UDSI_M1(c,t,ads)-C_UDSI_M1(c,t,q_ads));
                        }
                        end_c_loop(c,t)
                }
        }
}
DEFINE_ADJUST(d_adsorption,domain)
{
        Thread*t;
        cell_t c;
        /*solve the energy source.*/
        thread_loop_c(t,domain)
        {
           if(NULL!=THREAD_STORAGE(t,SV_UDS_I(d_ads))&&NULL!=T_STORAGE_R_NV(t,SV_UDSI_M1(d_ads)))
           {
                   begin_c_loop(c,t)
            {
           C_UDSI(c,t,d_ads)=-0.5422633*0.001*(C_UDSI(c,t,ads)-C_UDSI(c,t,q_ads));
        }   
        end_c_loop(c,t)
       }
    }
}
DEFINE_SOURCE(mass_source,c,t,dS,eqn)
{
        real m_source;
        m_source=C_UDSI_M1(c,t,d_ads);
        dS[eqn]=0.0;
        return m_source;
}
DEFINE_SOURCE(enery_source,c,t,dS,eqn)
{
        real e_source;
        e_source=-C_UDSI_M1(c,t,d_ads)*4*1.58E6;
        dS[eqn]=0.0;
        return e_source;
}
回复此楼

» 猜你喜欢

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

saw8124

新虫 (初入文坛)

【答案】应助回帖

大佬,请问您吸附调试出来了吗?我做的时候遇到了一些小问题,想请教
4楼2021-01-23 10:02:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

SuperUDF

银虫 (小有名气)

UDS没有看到对流项,扩散项,源项也不是UDS的,里面问题很大着呢。调试可以用VC++ UDF Studio软件(自己百度),一步一步跟踪变量的值。找人代做估计不靠谱,UDS没有几千搞不定。你有钱除外。

» 本帖已获得的红花(最新10朵)

2楼2019-12-02 13:54:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

白白的大萝卜

新虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by SuperUDF at 2019-12-02 13:54:29
UDS没有看到对流项,扩散项,源项也不是UDS的,里面问题很大着呢。调试可以用VC++ UDF Studio软件(自己百度),一步一步跟踪变量的值。找人代做估计不靠谱,UDS没有几千搞不定。你有钱除外。

好的,谢谢大佬!我去试一试吧
3楼2019-12-02 15:10:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 复试调剂 +6 Copy267 2026-03-10 6/300 2026-03-13 09:00 by Linda Hu
[考研] 0703化学求调剂 +3 绿豆芹菜汤 2026-03-12 3/150 2026-03-13 00:01 by 万物之源~氢
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-12 21:53 by 飞仙戟羽
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
[考研] 285求调剂 +3 柴郡猫_ 2026-03-12 3/150 2026-03-12 15:39 by 布丁aaa
[考研] 085701环境工程专硕求调剂 +8 幼稚园大班打手 2026-03-08 8/400 2026-03-12 15:12 by 每天只摆一小会
[考研] 271求调剂 +7 生如夏花… 2026-03-11 7/350 2026-03-12 09:32 by allen-yin
[考研] 材料工程330分求调剂,一志愿985 +5 小材化本科 2026-03-07 5/250 2026-03-11 16:51 by 535743368
[考研] 083000环境科学与工程调剂 +8 mingmingry 2026-03-09 9/450 2026-03-11 10:23 by 沙漠之狐994
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
[考研] 337一志愿华南理工材料求调剂 +5 mysdl 2026-03-07 5/250 2026-03-10 14:54 by ruiyingmiao
[考博] 找博导 +4 小呆呆熊 2026-03-07 4/200 2026-03-10 11:08 by 天道酬勤2026686
[基金申请] 面上项目还需要AI说明吗? +3 liyundong 2026-03-08 3/150 2026-03-09 22:30 by kingkocxr
[考研] 310 070300化学求调剂 +4 扑风铃的猫 2026-03-08 5/250 2026-03-09 19:46 by 30660438
[考研] 求0856调剂 +5 小力气珂珂 2026-03-08 5/250 2026-03-09 06:42 by houyaoxu
[考研] 求调剂,一志愿华中科大0702,数一英一,293 +4 小罗露一二 2026-03-07 4/200 2026-03-08 16:36 by 星空星月
[考研] 081700学硕一志愿北京化工大学数二英一过六级有竞赛求调剂 +5 galaxary 2026-03-07 7/350 2026-03-08 04:37 by wang_dand
[考研] 材料化工求调剂 +7 博斯特525 2026-03-06 8/400 2026-03-08 03:51 by lfhuang
[考研] 化学290求调剂 +3 两颗 西柚 2026-03-07 4/200 2026-03-07 20:55 by macy2011
[考研] 求调剂 +4 呼呼?~+123456 2026-03-06 4/200 2026-03-06 23:11 by L135790
信息提示
请填处理意见