Question

fluent中udf编程提示出错 Error: received a fatal signal (Segmentation fault). Error Object: #f

需要给图中空心桨盘施加动量源，以测得该桨盘转动对计算域流场的影响，程序如下#include "udf.h" /*Define the constant variables*/#define rho 1.225 #define pi 3.14159 #define omega 1164 /*For certain propeller with a 10000 RPM*/ #define theta0 11 #define dr 0.07 /*This should based on the mesh we generated*/ #define dophy 9.0 #define r0 0.225#define r1 0.925#define c 0.1 DEFINE_SOURCE(sourcex,cell,thread,dS,eqn){ real x[ND_ND], x1,y1,z1,r, v1,v2,v3,vx,vy,vz,v; real phi,alpha,ophy,cy,cx,l,d,Fn,Ft,Fx,Fy,Fz; real sx,sy,sz,source; begin_c_loop(cell,thread){ C_CENTROID(x,cell,thread); x1 = x[0]; y1 = x[1];if(sqrt(x1*x1+y1*y1)>0&&sqrt(x1*x1+y1*y1)<r0)r=0;else if(sqrt(x1*x1+y1*y1)==r0&&sqrt(x1*x1+y1*y1)<r1&&sqrt(x1*x1+y1*y1)==r1)r=sqrt(x1*x1+y1*y1);else if(sqrt(x1*x1+y1*y1)>r1)r=0; if (x1>0.&&y1>0.)ophy=asin(fabs(y1/r)); else if (x1<0.&&y1>0.)ophy=asin(fabs(x1/r))+pi/2.; else if (x1<0.&&y1<0.)ophy=asin(fabs(y1/r))+pi; else if (x1>0.&&y1<0.)ophy=asin(fabs(x1/r))+3.*pi/2.; else if (x1==0.&&y1<0.)ophy=3./2.*pi; else if (x1==0.&&y1>0.)ophy=0.5*pi; else if (x1>0.&&y1==0.)ophy=0.; else if (x1<0.&&y1==0.)ophy=pi; v1=C_U(cell,thread); v2=C_V(cell,thread); v3=C_W(cell,thread); vx=v1+omega*y1;vy=v2-omega*x1;vz=v3;v=sqrt(vx*vx+vy*vy+vz*vz);phi=atan(-vz/(vx*vx+vy*vy));alpha=theta0-phi; cy=0.11*alpha*180.0/pi; cx=0.001*alpha*180.0/pi; l=0.5*rho*v*v*cy*c*dr; d=0.5*rho*v*v*cx*c*dr; Fn=l*cos(phi)-d*sin(phi); Ft=l*sin(phi)-d*cos(phi);Fx=Ft*sin(ophy);Fy=-Ft*cos(ophy);Fz=Fn; sx=-dophy*Fx/pi;sy=-dophy*Fy/pi;sz=-dophy*Fz/pi;source=sx; C_UDMI(cell,thread,0)=source;} end_c_loop(cell,thread) dS[eqn]=0.0; return source; } 请问这里是有什么问题吗？编译udf时没有报错，一初始化运行就报错

Answer 1

1.先初始化，迭代一步后，再导入udf--问题依旧
2.先保存退出软件后，再重新打开--问题依旧
3.将4个udf分别拆开后，分别导入--没有出现问题
4.将其中的两个udf组合起来，导入--问题依旧（据此推断是udf的问题），下面是我的udf；将DEFINE_SOURCE(liq_src…）和DEFINE_SOURCE(vap_src…）这两个udf组合在一起，出现了开始的同样问题。请大神帮忙看下，问题出在哪？滴水之恩，涌泉相报，谢谢大神们！！！
附件截图是我引用udf的三个位置，附件文件是我的udf。

#include "udf.h" 
#include "sg_mphase.h"
#define T_SAT 278 
#define LAT_HT 2260000
DEFINE_SPECIFIC_HEAT(my_user_cp, T, Tref, h, yi) 
{ 
real cp=1200.176+2164.7398*pow(0.35,1.16317)+0.0229*pow(T,1.306682)+0.3285*pow(0.35,1.16317)*pow(T,1.306682); 
*h = cp*(T-Tref); 
return cp; 
}


DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn)//与下面组合后出现问题
{
Thread *mix_th,*sec_th;
real m_dot_l;
mix_th=THREAD_SUPER_THREAD(pri_th);
sec_th=THREAD_SUB_THREAD(mix_th,1);
if(C_T(cell,mix_th)>=T_SAT)
{
m_dot_l=-0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;
dS[eqn]=-0.1*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;
}
else
{
m_dot_l=0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
dS[eqn]=0.0;
}
return m_dot_l;
}

DEFINE_SOURCE(vap_src,cell,sec_th,dS,eqn)//与上面组合后出现问题 
{
Thread *mix_th,*pri_th;
real m_dot_v;
mix_th=THREAD_SUPER_THREAD(sec_th);
pri_th=THREAD_SUB_THREAD(mix_th,0);
if(C_T(cell, mix_th)>=T_SAT)
{
m_dot_v=0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;
dS[eqn]=0.0;
}
else
{
m_dot_v=-0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
dS[eqn]=-0.1*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;
}
return m_dot_v;
}

DEFINE_SOURCE(enrg_src,cell,mix_th,dS,eqn)
{
Thread *pri_th,*sec_th;
real m_dot;
pri_th=THREAD_SUB_THREAD(mix_th,0);
sec_th=THREAD_SUB_THREAD(mix_th,1);
if(C_T(cell, mix_th)>=T_SAT)
{
m_dot=-0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;
dS[eqn]=-0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)/T_SAT;
} 
else
{
m_dot=0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
dS[eqn]=-0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)/T_SAT;
}
return LAT_HT*m_dot;  
}


QQ1.png



QQ2.png

Answer 2

需要把fluent/define/memory开了，你这个没开