| 查看: 1870 | 回复: 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人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有218人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
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