| 查看: 1772 | 回复: 4 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
添加连续性方程质量源项UDF编译无错误，进行计算时出现错误 已有3人参与
|
|||
|
本人目前所做内容为研究地埋管周围土壤热湿耦合传递特性，现添加连续性方程源项UDF编译无错误，进行计算时出现错误： Error: received a fatal signal (Segmentation fault). Error: received a fatal signal (Segmentation fault). Error Object: #f 哪位大神能帮忙看看是什么地方得错误，万分感谢！！！ 连续性方程源项UDF如下： #include"udf.h" DEFINE_SOURCE(source_mass,cell,t,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]; real x; Domain *domain; Thread *c_thread; Thread *subthread; cell_t c; face_t f; Domain *d=Get_Domain(2); //第一相对应的子域 domain_id 是 2 Thread *mixture_thread = THREAD_SUPER_THREAD(subthread); temp0=C_T_G(cell,c_thread)[0]; temp1=C_T_G(cell,c_thread)[1]; temp2=C_T_G(cell,c_thread)[2]; vol=C_VOLUME(cell,c_thread); //控制体体积； vof=C_UDSI(cell,c_thread,0); //控制体内含湿量； x=vof*12.53; diff=(7*(-2.41-0.002*pow(vof,-1.75))/0.0718)*(-1.425*pow(10,-4)-4.76*pow(10,-7)*C_T(c,t))*9.48*pow(10,x)/pow(10,10); //等温扩散系数 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,c_thread) //历遍混合相的所有面 { F_AREA(A,f,c_thread); //F_AREA可以被用于返回一个实数的面积向量 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,c_thread) source=diff*(a0[0]+a1[1]+a2[2])/vol; dS[eqn]=0; return source; } |
回复此楼» 猜你喜欢
重庆交通大学光子学微结构与器件课题组招收2026年硕士研究生信息
已经有2人回复
-
一志愿郑大材料学硕298分，求调剂
已经有6人回复
-
物理学I论文润色/翻译怎么收费?
已经有200人回复
-
寻合作：应力腐蚀多尺度模拟
已经有5人回复
-
考研交流
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
帮我的英语口语老师找学生
已经有1人回复
3楼2019-04-22 12:15:38
2楼2019-04-11 19:38:26
4楼2019-04-24 23:02:21
