| 查看: 2106 | 回复: 9 | |||
[交流]
关于输出捕获的颗粒物的坐标问题
|
|||
|
最近在做颗粒物的捕获。即当颗粒物碰到壁面后判定颗粒物速度，然后进行捕获或者反射，并将捕获的颗粒物坐标输出到文件，UDF如下 #include "udf.h" #include "dpm.h" FILE *fp; DEFINE_DPM_BC(JIEZHA,p,t,f,f_normal,dim) { fp=fopen("d:\\out.txt","a+"  ;if(NV_MAG(P_VEL(p))<8) { fprintf(fp," %g %g %g\n",P_POS(p)[0],P_POS(p)[1],P_POS(p)[2]); return PATH_ABORT; fclose(fp); } else return PATH_ACTIVE; } 然而在进行非稳态计算后，计算总是会卡在Advancing DPM injections ....，out.txt文件会被创建，但是里面并没有数据，请问是什么原因造成的，我的UDF有什么问题吗 |
回复此楼» 猜你喜欢
遇见不省心的家人很难过
已经有16人回复
-
退学或坚持读
已经有25人回复
-
博士延得我，科研能力直往上蹿
已经有4人回复
-
免疫学博士有名额，速联系
已经有14人回复
-
面上基金申报没有其他的参与者成吗
已经有4人回复
-
多组分精馏求助
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 三、关于“天然动力能量场核”构造推想
已经有0人回复
» 抢金币啦！回帖就可以得到:
-
南方医科大学中药学院 申请考核博士一名 （天然药化方向，天然产物分离经验优先）
+1/272
-
中国石油大学（华东）吴传德教授团队（国家杰青）2026硕、博招生
+2/260
-
香港岭南大学郭瑛课题组博士招生 (2026年秋季入学）
+1/185
-
黄汉民团队联合淮北师范大学招聘师资博士后（年薪30-40万）
+1/84
-
深圳大学柔性电子材料方向“申请-考核制”博士生招生
+2/66
-
娃娃们今儿考试喽。。。。
+1/61
-
香港科技大学计算物理及流体力学课题组招收全奖博士后及博士生（2026年9月入学）
+1/41
-
香港科技大学计算物理及流体力学课题组招收全奖博士后及博士生（2026年9月入学）
+1/38
-
香港科技大学计算物理及流体力学课题组招收全奖博士后及博士生（2026年9月入学）
+1/36
-
教育部重点实验室和清华大学某国家重点实验室，联合培养硕生、博生，并长期招博士后
+1/30
-
山东科技大学招聘化学化工博士博士后
+1/28
-
德国Karlsruhe Institute of Technology招收电化学储能及联合培养CSC博士
+1/15
-
中科院深圳先进院-免疫治疗方向-招收1名博士生(26年9月入学)
+1/13
-
考博求助
+1/13
-
南方科技大学仉晶晶课题组诚聘有机方向博士后、科研助理
+1/8
-
上海理工大学“新能源材料”专业-赵斌教授招收申请考核制博士生【能源催化方向】
+1/2
-
苏州大学医学院纳米生物医学方向招收申请考核制博士生1名
+1/1
-
澳科大诚招2026年秋季全奖博士研究生（药剂学/生物材料方向）
+1/1
-
211大学-课题组招收学术博士生
+1/1
-
接理论计算，主要包含第一性原理、分子动力学、机器学习、有限元模拟等，欢迎学术交流
+1/1
2楼2017-04-13 11:32:25
4楼2017-04-14 13:33:22
5楼2017-04-27 14:16:36
★
苦逼锅炉工(金币+1): 谢谢参与
苦逼锅炉工(金币+1): 谢谢参与
|
要我说吧，这个问题应该这么看: 1. 能生成out.txt文件，那么这个UDF能够正常调用并至少能够运行至 fp=xxxxx 这一行； 2.能运行fp=xxxxx这一行照理来说没理由不能继续运行if这一行； 3.然而运行了if这一行之后没有记录结果，说明没有任何一次if条件判断是True ； 4.于是乎得出结论，你的判断条件有问题，没有任何一句 NV_MAG(P_VEL(p))<8； 5.怎么验证？把NV_MAG(P_VEL(p))<8中的这个8设定成88888888888888 6.如果这样修改后out.txt中依旧没有数据，那么，我上面说的都是废话(′;ω;｀) 7.再次ps：fclos(fp）这一句永远不会被执行；只要运行到Return就表示本次对该程序的调用结束，并返回Return中的值，Return之后的语句不会被执行； 8.在另一个帖子里回复你关于输出好多点的问题。 9.对8的补充：首先还是确认一下你的case是只在某一个时间点射入制定数量的颗粒，有可能你的case一直在入射颗粒所以才会一直trap。 |
7楼2017-04-27 16:41:04
8楼2017-06-08 09:45:26
9楼2017-06-08 16:12:55
10楼2022-05-23 10:52:00
简单回复
苦逼锅炉工3楼
2017-04-13 15:42
回复
2017-04-27 14:52
回复
苦逼锅炉工(金币+1): 谢谢参与
顶 发自小木虫IOS客户端