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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » FLUENT单机多核并行运算导入udf之后，初始化报错，Received signal SIGSEGV
51/1返回列表
查看: 4967  |  回复: 19
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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;
}
回复此楼

» 猜你喜欢

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

SJTU_SZ

金虫 (小有名气)
引用回帖:
5楼: Originally posted by gafeide at 2015-06-23 19:04:38
但是你之前没有声明这个变量啊
...

比热的函数里面就有这个变量，不用声明的，就跟T和Tref一样。
赞一下 回复此楼
科研！科研！
6楼2015-06-26 08:36:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答

SJTU_SZ

金虫 (小有名气)
真心求高人相助，弄了好几天了，没进展！
赞一下(1人) 回复此楼
科研！科研！
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的回帖
查看全部 20 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿085802 323分求调剂 +12 drizzle_9 2026-04-12 13/650 2026-04-13 01:15 by 幸免 ..
[考研] 322求调剂，08工科 +4 今天是个小号 2026-04-08 4/200 2026-04-13 00:20 by baobaoye
[考研] 22408 352分求调剂 +5 努力的夏末 2026-04-09 5/250 2026-04-12 19:17 by wj165256
[硕博家园] 新一代电子信息294求调剂 不挑学校 +7 Ytyt11 2026-04-09 8/400 2026-04-12 16:57 by ajpv风雷
[考研] 310求调剂 +13 666真好 2026-04-11 14/700 2026-04-12 16:48 by 1005715100
[考研] 085501机械专硕 302分 不挑专业求调剂 +7 汪某. 2026-04-09 7/350 2026-04-11 14:37 by luhong1990
[考研] 085410 273分调剂 +4 X1999 2026-04-09 4/200 2026-04-11 13:05 by pies112
[考研] 还有化工二轮调剂的学校吗 5+14 化工人999 2026-04-09 48/2400 2026-04-11 10:27 by 89436494
[考研] 311求调剂 +13 xyp想读书 2026-04-10 14/700 2026-04-11 09:41 by 猪会飞
[考研] 080500求调剂 +17 黄宇博 2026-04-06 17/850 2026-04-11 08:36 by zhq0425
[考研] 085404 298分求调剂 +10 呼啦呼啦呼呼呼 2026-04-10 11/550 2026-04-10 16:44 by wangy0907
[考研] 273求调剂 +51 麦小叮当 2026-04-06 58/2900 2026-04-10 15:54 by jiajinhpu
[考研] 298求调剂 +13 钉叮咚冬瓜 2026-04-09 13/650 2026-04-10 15:49 by jiajinhpu
[考研] 材料专硕调剂 +16 哈哈哈吼吼吼哈 2026-04-07 17/850 2026-04-09 21:16 by wutongshun
[考研] 085801 总分275 本科新能源 求调剂 +8 bradoner 2026-04-08 9/450 2026-04-09 13:43 by only周
[考研] 327求调剂 +10 Xxjc1107. 2026-04-06 11/550 2026-04-09 01:21 by lature00
[考研] 277求调剂 +4 考研调剂lxh 2026-04-06 6/300 2026-04-08 10:40 by 逆水乘风
[考研] 机械工程264学硕求调剂 +3 qiushangxian 2026-04-06 3/150 2026-04-08 01:53 by Linzejun
[考研] 调剂求助（生物与医药） +6 @6952 2026-04-06 6/300 2026-04-07 23:52 by lys0704
[考研] 11408 325分 +3 jgtxuxgkx 2026-04-07 3/150 2026-04-07 23:10 by lbsjt
信息提示
关闭
请填处理意见
关闭