| 查看: 1832 | 回复: 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; } |
回复此楼» 猜你喜欢
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有198人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有11人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有12人回复
-
桂林理工大学物理学专业招收调剂
已经有13人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
美国顶级物理期刊《应用物理快报》（APL）的编辑欺骗和歧视作者及AIP的官僚主义傲慢
已经有6人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
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
