| 查看: 1840 | 回复: 4 | |||
[求助]
udf程序改为并行
|
|
程序的目的是根据监测点的压力调节入口的质量流量，维持监测点的压力近似恒定，单核计算时程序可以用，一旦并行就出错，自己试着按帮助里的方法改一直没通（比如加#if !HOST #endif之类），恳请各位指点，先谢过了！把自己的金币全部悬赏！ DEFINE_ADJUST(get_p_monitor,d) { Thread *t; cell_t c; real xc[ND_ND];/*网格质心坐标*/ thread_loop_c(t,d) { begin_c_loop(c,t) { C_CENTROID(xc,c,t); if (sqrt(ND_SUM(pow(xc[0] + 2.,2.), pow(xc[1] -1.3，2.)， pow(xc[2]-0.,2.)))<0.01 pressure=C_P(c,t);/*点(-2,1.3,0)为监测点，当网格质心与监测点的距离小于0.01时，认为该网格的压力就是监测点的压力*/ ｝ end_c_loop(c,t) } printf("\n pressure=%g\n",pressure); } DEFINE_PROFILE(mass_inlet_as_p,t,i) { face_t f; real mass=0.5; if(pressure<=480000.) mass=mass+0.01; else if(pressure>480000.&&.pressure<520000.)/*根据监测点的压力调节入口质量流量*/ mass=mass; else mass=mass-0.01; begin_f_loop(f,t) { F_PROFILE(f,t,i)=mass; } end_f_loop(f,t) printf("\n pressure=%g\n",pressure); } 如果一并行，每计算一步会出现四个压力值（我用四核并行），其中一个为真实值，另三个是0，然后到下面用压力判断质量流量时用的一直是那个0的压力，所以就流量一直升高。 |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 动网格的UDF问题
已经有21人回复
- udf在非host节点下并行计算message
已经有11人回复
- FLUENT中编写UDF时指针变量是否要初始化
已经有4人回复
- 单机多核并行计算下UDF的问题
已经有9人回复
- 一个UDF并行报错
已经有17人回复
- fluent UDF代码调用预定义程序的问题（如fortran程序）
已经有15人回复
- 进行单机多核udf的并行计算，需要怎么设置？
已经有9人回复
- 求助udf的输入输出
已经有16人回复
- UDF 改为并行计算
已经有9人回复
- 并行计算没有输出？怎么回事
已经有16人回复
- UDF 热源随温度变化的处理
已经有4人回复
- 单机多核一个cpu的串行并行求助
已经有26人回复
- 激光焊接udf，加载温度场不对
已经有7人回复
- FLUENT读取并行case自动删除文件
已经有3人回复
- 关于MPI并行编程的基本问题
已经有8人回复
- 【求助】单机多核fluent并行udf
已经有13人回复
- 【求助】单机多核UDF怎么并行计算
已经有15人回复
2楼2013-11-06 18:00:08
3楼2013-11-07 20:23:45
4楼2013-11-14 19:19:26
【答案】应助回帖
★ ★ ★ ★ ★
1592203609: 金币+5, 谢谢回帖 2013-11-15 21:29:15
1592203609: 金币+5, 谢谢回帖 2013-11-15 21:29:15
|
DEFINE_ADJUST(get_p_monitor,d) { Thread *t; cell_t c; real xc[ND_ND]; thread_loop_c(t,d) { begin_c_loop(c,t) { C_CENTROID(xc,c,t); if (sqrt(ND_SUM(pow(xc[0] + 2.,2.), pow(xc[1] -1.3，2.)， pow(xc[2]-0.,2.)))<0.01 pressure=C_P(c,t); /*你这里前面网格的压力会不断的被后面网格所覆盖，Fluent会循环所有网格，所以最后输出的pressure变量只是最后一个网格的压力。可把C_P(c,t)放在一个UDM里面（自己查查，不难）*/ ｝ end_c_loop(c,t) } printf("\n pressure=%g\n",pressure); /*由于并行时流场会被分为几份，输出时有的区域的压力就会输出为零。pressure初始化为零？省略了？*/ } DEFINE_PROFILE(mass_inlet_as_p,t,i) { face_t f; real mass=0.5; if(pressure<=480000.) mass=mass+0.01; else if(pressure>480000.&&.pressure<520000.)/*根据监测点的压力调节入口质量流量*/ mass=mass; else mass=mass-0.01; begin_f_loop(f,t) { F_PROFILE(f,t,i)=mass; } end_f_loop(f,t) printf("\n pressure=%g\n",pressure); /*这种写法输出不会崩溃？这里不是每个网格都输出？我记得DEFINE_PROFILE会循环每个网格，也就是每个网格应该都会执行printf啊？*/ } |
5楼2013-11-15 17:55:42
25