| 查看: 2034 | 回复: 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; } |
回复此楼» 猜你喜欢
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有149人回复
-
课题组招1-2名硕士
已经有2人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有11人回复
3楼2019-04-22 12:15:38
2楼2019-04-11 19:38:26
4楼2019-04-24 23:02:21
