24小时热门版块排行榜    

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

SJTU_SZ

金虫 (小有名气)

[求助] FLUENT单机多核并行运算导入udf之后,初始化报错,Received signal SIGSEGV 已有1人参与

用FLUENT的圆管内超临界水流动传热,采用SST模型,用udf在material物性设置里面定义超临界水的物性随温度的变化,然后初始化的时候报错,错误信息如下:
==============================================================================

Node 4: Process 19756: Received signal SIGSEGV.

==============================================================================
MPI Application rank 0 exited before MPI_Finalize() with status 2
The fl process could not be started.
下面是我的udf代码,udf直接在记事本里写的,然后把扩展名改成.c   ,请问是udf的问题,还是我别的什么地方弄错了?PS,我串行运算的时候也报错Error:  received a fatal signal (Segmentation fault).
Error Object: #f,求大神指教!

#include "udf.h"
DEFINE_PROPERTY(viscosity_polynomial,cell,thread)
{
real viscosity;
real temperature;
temperature=C_T(cell,thread);
if(temperature>=540&&temperature<670)
  viscosity=-0.02974+(2.06393E-4)*temperature-(5.32888E-7)*pow(temperature,2)+(6.09625E-10)*pow(temperature,3)-(2.61186E-3)*pow(temperature,4);
else if(temperature>=670&&temperature<720)
  viscosity=1.71089-0.00967*temperature+(2.05003E-5)*pow(temperature,2)-(1.93154E-8)*pow(temperature,3)+(6.82497E-12)*pow(temperature,4);
else
  viscosity=4.60518E-4-(1.89262E-6)*temperature+(3.04391E-9)*pow(temperature,2)-(2.12799E-12)*pow(temperature,3)+(5.55983E-16)*pow(temperature,4);
return viscosity;
}
DEFINE_PROPERTY(density_polynomial,cell,thread)
{
real density;
real temperature;
temperature=C_T(cell,thread);
if(temperature>=540&&temperature<650)
  density=-102064.64697+719.52081*temperature-1.88404*pow(temperature,2)+0.00219*pow(temperature,3)-(9.59734E-7)*pow(temperature,4);
else if(temperature>=650&&temperature<680)
  density=1.05818E8-637783.40085*temperature-1441.25876*pow(temperature,2)-1.44725*pow(temperature,3)+(5.44861E-4)*pow(temperature,4);
else if(temperature>=680&&temperature<720)
  density=1.45053E7-81799.75493*temperature+173.01793*pow(temperature,2)-0.16267*pow(temperature,3)+(5.7363E-5)*pow(temperature,4);
else
  density=9785.72376-39.65886*temperature+0.06117*pow(temperature,2)-(4.2103E-5)*pow(temperature,3)+(1.08805E-8)*pow(temperature,4);
return density;
}
DEFINE_PROPERTY(thermal_conductivity_polynomial,cell,thread)
{
real tc;
real temperature;
temperature=C_T(cell,thread);
if(temperature>=540&&temperature<650)
  tc=-46.31273+0.31446*temperature-(7.86857E-4)*pow(temperature,2)+(8.75221E-7)*pow(temperature,3)-(3.6729E-10)*pow(temperature,4);
else if(temperature>=650&&temperature<680)
  tc=-10013.68647+61.91813*temperature-0.14356*pow(temperature,2)+(1.47931E-4)*pow(temperature,3)-(5.7166E-8)*pow(temperature,4);
else if(temperature>=680&&temperature<750)
  tc=4105.15866-22.5438*temperature+0.04644*pow(temperature,2)-(4.25217E-5)*pow(temperature,3)+(1.46033E-8)*pow(temperature,4);
else
  tc=7.23134-0.02932*temperature+(4.49844E-5)*pow(temperature,2)-(3.0539E-8)*pow(temperature,3)+(7.76698E-12)*pow(temperature,4);
return tc;
}
DEFINE_SPECIFIC_HEAT(specific_heat_polynomial,T,Tref,h,yi)
{
real cp;
if(T>=540&&T<650)
  cp=1.02221E7-70713.70545*T+183.47921*pow(T,2)-0.21157*pow(T,3)+(9.15037E-5)*pow(T,4);
else if(T>=650&&T<675)
  cp=-2.81076E10+(1.69883E8)*T-384994.92162*pow(T,2)+387.72315*pow(T,3)-0.14641*pow(T,4);
else if(T>=675&&T<710)
  cp=-5.27703E9+3.05849E7*T-66448.74541*pow(T,2)+64.13917*pow(T,3)-0.02321*pow(T,4);
else if(T>=710&&T<900)
  cp=6.4899E6-31048.84474*T+55.78053*pow(T,2)-0.04456*pow(T,3)+(1.33526E-5)*pow(T,4);
else
  cp=117954.93316-415.67508*T+0.56772*pow(T,2)-(3.47613E-4)*pow(T,3)+(8.04047E-8)*pow(T,4);
*h=cp*(T-Tref);
return cp;
}
回复此楼
科研!科研!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿_凡

木虫 (小有名气)

引用回帖:
15楼: Originally posted by 叶子闹 at 2018-04-04 15:13:41
请问你的问题解决了吗?我最近在编译导入UDF初始化之后也遇到了同样的问题...

请问您的问题解决了吗?我最近也碰到同样的问题
19楼2020-03-04 20:48:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答

SJTU_SZ

金虫 (小有名气)

真心求高人相助,弄了好几天了,没进展!
科研!科研!
2楼2015-06-23 08:43:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gafeide

铜虫 (初入文坛)

倒数第三行,*h是什么意思啊

[ 发自小木虫客户端 ]
3楼2015-06-23 13:00:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

SJTU_SZ

金虫 (小有名气)

引用回帖:
3楼: Originally posted by gafeide at 2015-06-23 13:00:38
倒数第三行,*h是什么意思啊

显焓的值,我看fluent帮助里面这么写的,要求定义这么一个形式。
科研!科研!
4楼2015-06-23 14:20:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-25 01:47 by BZKMTicpDhFj
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +7 scalable 2026-08-24 7/350 2026-08-24 23:01 by anata1209
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见