| 查看: 1805 | 回复: 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; } |
回复此楼» 猜你喜欢
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有80人回复
-
物理学I论文润色/翻译怎么收费?
已经有89人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
2楼2019-04-11 19:38:26
3楼2019-04-22 12:15:38
4楼2019-04-24 23:02:21
刘广义
铁杆木虫 (著名写手)
- 应助: 60 (初中生)
- 金币: 4971.7
- 散金: 25
- 红花: 18
- 帖子: 1972
- 在线: 572小时
- 虫号: 2889824
- 注册: 2013-12-23
- 性别: GG
- 专业: 燃烧学
5楼2021-10-03 10:20:47
