| 查看: 1266 | 回复: 6 | ||
[求助]
udf致命错误 已有1人参与
|
|
我仿真的模型是表层土壤浅埋换热管，研究非饱和土壤中热湿耦合 迁移，现阶段只考虑液态水的迁移。模型有混合mixture模型，porous模型，利用udf编写方程，程序如下： #include"udf.h" #include<math.h> DEFINE_UDS_UNSTEADY(myunsteady,c,t,i,apu,su) { real physical_dt,vol,phi_old; physical_dt=RP_Get_Real("physical_time_step"  ;vol=C_VOLUME(c,t); *apu=-vol/physical_dt; phi_old=C_STORAGE_R(c,t,SV_UDSI_M1(i)); *su=vol*phi_old/physical_dt; } DEFINE_SOURCE(source_moving,cell,sec_th,dS,eqn) { real A[ND_ND]; real source; real temp0,temp1,temp2; real vol; real vof; real a0[ND_ND],a1[ND_ND],a2[ND_ND]; real diff; real xv[ND_ND],yv[ND_ND],zv[ND_ND]; Domain*d; Thread*mix_th; face_t f; d=Get_Domain(2); mix_th=THREAD_SUPER_THREAD(sec_th); //指针指向混合项的第二项即液态水 temp0=C_T_G(cell,mix_th)[0]; temp1=C_T_G(cell,mix_th)[1]; temp2=C_T_G(cell,mix_th)[2]; vol=C_VOLUME(cell,mix_th); //控制体体积； vof=C_UDSI(cell,mix_th,0); //控制体 内含湿量； diff=-2.3495*pow(10,-5)*pow(vof,6)+2.0370*pow(10,-5)*pow(vof,5)-6.6946*pow(10,-6)*pow(vof,4)+1.0438*pow(10,-6)*pow(vof,3)-7.8699*pow(10,-8)*pow(vof,2)+2.7350*pow(10,-9)*vof+5.4039*pow(10,-12); //等温扩散系数 xv[0]=1;xv[1]=0;xv[2]=0; yv[0]=0;yv[1]=1;yv[2]=0; zv[0]=0;zv[1]=0;zv[2]=1; begin_f_loop(f,mix_th) //历遍混合相的 所有面； { F_AREA(A,f,mix_th); a0[ND_ND]+=temp0*NV_DOT(A,xv); //x方向的温度梯度与控制体x方向 投影面积的乘积； a1[ND_ND]+=temp1*NV_DOT(A,yv); a2[ND_ND]+=temp2*NV_DOT(A,zv); } end_f_loop(f,mix_th) source=diff*(a0[0]+a1[1]+a2[2])/vol; dS[eqn]=0; return source; } DEFINE_DIFFUSIVITY(MyDiff,c,t,i) { real vof; real con_zl; real con_nd; real MyDiff; vof=C_UDSI(c,t,0); con_zl=0.075636667-0.000114214*C_T(c,t)-0.000000279*C_T(c,t)*C_T(c,t); con_nd=0.00041212*exp(-C_T(c,t)/12.37666)+0.00119916*exp(-C_T(c,t)/41.70239)+0.00018058; MyDiff=9.15003*pow(vof,3.613655)*con_zl*0.0014953452/(con_nd*0.073240827); return MyDiff; } 在进行计算的时候显示致命错误，求大神指点  @wuming524 |
回复此楼» 猜你喜欢
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有14人回复
-
网上报道青年教师午睡中猝死、熬夜猝死的越来越多，主要哪些原因引起的？
已经有4人回复
-
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有10人回复
-
版面费该交吗
已经有12人回复
-
面上可以超过30页吧？
已经有5人回复
-
“人文社科而论，许多学术研究还没有达到民国时期的水平”
已经有5人回复
-
什么是人一生最重要的？
已经有4人回复
2楼2016-04-19 17:27:41
3楼2016-04-20 10:30:39
4楼2016-04-25 22:40:03
5楼2017-03-12 10:02:17
6楼2019-01-10 13:54:53
7楼2019-01-10 14:07:12