| 查看: 4974 | 回复: 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; } |
回复此楼» 猜你喜欢
085600材料与化工349分求调剂
已经有16人回复
-
085801电气专硕272求调剂
已经有14人回复
-
271求调剂
已经有27人回复
-
调剂求收留
已经有33人回复
-
300分求调剂 （085501机械专硕，本科扬大）
已经有10人回复
-
求调剂
已经有15人回复
-
材料085601调剂
已经有34人回复
-
机械还有还有名额吗？太难了
已经有8人回复
-
288求调剂，一志愿华南理工大学071005
已经有24人回复
-
085600材料与化工329分求调剂
已经有23人回复
15楼2018-04-04 15:13:41
SJTU_SZ
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 579.5
- 帖子: 54
- 在线: 32.7小时
- 虫号: 2098217
- 注册: 2012-10-31
- 性别: GG
- 专业: 多相流热物理学
2楼2015-06-23 08:43:25
3楼2015-06-23 13:00:38
SJTU_SZ
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 579.5
- 帖子: 54
- 在线: 32.7小时
- 虫号: 2098217
- 注册: 2012-10-31
- 性别: GG
- 专业: 多相流热物理学
4楼2015-06-23 14:20:42
5