| 查看: 1657 | 回复: 3 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 Akun 的 10 个金币 | |||
[交流]
有关动网格面被分割，该怎么实现并行下使用
|
|||
|
最近在做颗粒沉积方面的模拟，用到了define_dpm_erosion这个宏，模拟的情况是圆管壁面沉积。 1.模拟过程包含两部分：a.erosion宏的沉积判定；b.管壁面根据沉积高度进行define_grid_motion的动网格执行； 2.因为流域颗粒比较多，所以计划采用多核进行计算，缩短计算时间； 3.用八核计算时，发现管壁面被分成了4各部分，因为grid_motion中要用到begin_f_loop，是否需要进行并行化语句； 4.目前serial模式没有问题，并行会因为动网格更新时出现负体积而报错，我猜是因为并行将壁面分割的原因（问题可能就出现在分割面重合的那几个节点上） 以下是我grid_motion中的部分语言：中间我还用到了数组计算每个节点的位移高度 DEFINE_GRID_MOTION(dongwangge,domain,dt,time,dtime) { Thread *tf=DT_THREAD(dt); face_t f; int i,j; int m,n,l; real cos_theta,sin_theta; real sum=0; int step; int cha; FILE *fp; real vf[1000][2]; real vn[1000][2]; real angle[1000]; real node[1000][2]; real A[ND_ND]; real A1[ND_ND]; Node *v; Node *vs; Thread *tf0; cell_t c0; step=N_TIME; cha=step-(step/200)*200; if(cha==0) { .......(利用数组计算了每个节点应该位移的高度) SET_DEFORMING_THREAD_FLAG(THREAD_T0(tf));// 开始动网格 i=0; n=0; begin_f_loop(f,tf) { f_node_loop(f,tf,n) { v=F_NODE(f,tf,n); vs=F_NODE_SHADOW(f,tf,n); if(NODE_POS_NEED_UPDATE(v)) { NODE_POS_UPDATED(v); cos_theta=vn[0]/NV_MAG(vn); sin_theta=vn[1]/NV_MAG(vn); if(angle<0.965926) { NODE_X(v)=(node[i-1][0]+node[0]+node[i+1][0])/3; NODE_Y(v)=(node[i-1][1]+node[1]+node[i+1][1])/3;//把节点位置移动至三角形重心位置 NODE_X(vs)=(node[i-1][0]+node[0]+node[i+1][0])/3; NODE_Y(vs)=(node[i-1][1]+node[1]+node[i+1][1])/3; } NODE_X(v)=NODE_X(v)-HEIGHT_NODE*cos_theta; NODE_Y(v)=NODE_Y(v)-HEIGHT_NODE*sin_theta; //网格位置更新 NODE_X(vs)=NODE_X(vs)-HEIGHT_NODE*cos_theta; NODE_Y(vs)=NODE_Y(vs)-HEIGHT_NODE*sin_theta; i++; } } } end_f_loop(f,tf) begin_f_loop(f,tf) { c0 = F_C0(f, tf); tf0 = F_C0_THREAD(f,tf); for (i=1; i<9; i++) { F_UDMI(f,tf,i)=0.; C_UDMI(c0,tf0,i)=0.; } } end_f_loop(f, tf) } else return; }  八核分布图.png  壁面四核分布图.png |
回复此楼» 猜你喜欢
求Sci期刊推荐
已经有6人回复
-
在大地上我们只过一生---看完我的阿勒泰上头了好几天，完结那天晚上几乎失眠
已经有11人回复
-
导师不让硕转博，让我去国外读博，能理解吗？
已经有18人回复
-
研0二导师分到新来的博士后靠谱吗
已经有4人回复
-
化学口B0109(高分子合成)，拿青年基金一般需要怎样的文章水平？
已经有18人回复
-
申请2024或2025年博士研究生
已经有9人回复
-
纠结选哪一个期刊，电化学领域
已经有10人回复
-
MMA预聚体光固化发雾问题求助
已经有10人回复
-
苯酚和炔溴反应成醚
已经有12人回复
-
双非本科毕业论文，气人
已经有4人回复
2楼2018-10-10 12:08:32
3楼2022-05-02 21:08:39
4楼2022-05-02 21:19:36
5