如何使用UDF添加颗粒间相互作用力（dpm模型）

lgsh

请教：如何使用udf添加颗粒间相互作用力的项呢？有一个关于体积力的宏（DEFINE_DPM_BODY_FORCE）,但是如何考虑被跟踪颗粒周围的多个颗粒间的相互作用呢，有些纠结~请问该从哪方面入手呢？

lgsh

可有人做过附加体积力算例的？给点意见呗~~~~~~~
#include "udf.h"
#include "surf.h"
#define A  4.e-19               
DEFINE_DPM_BODY_FORCE(Van_der_waals_force,p,i)
{
       
        cell_t c = RP_CELL(&(p -> cCell));
        Thread *t = RP_THREAD(&(p -> cCell));
        Particle *pi;  

        float bforce=0, bforce_x=0, bforce_y=0, bforce_z=0, coo=0, s=0;
        int   j=0, num=0;
        float pos[4]={0};

        begin_particle_cell_loop(pi,c,t)
        {                       
                num++;  
        }       
        end_particle_cell_loop(pi,c,t)

        if(num>0)
        {
                begin_particle_cell_loop(pi,c,t)
                {
                        for(j=0; j<3; j++)
                                pos[j] = P_POS(pi)[j] - P_POS(p)[j];
                        coo = sqrt(pos[0]*pos[0] + pos[1]*pos[1] + pos[2]*pos[2]);
                        s = coo-0.5*(P_DIAM(p)+P_DIAM(pi));

               
                        if(s <= 2.e-6)
                        {
                                bforce = A*P_DIAM(p)/(24*s*s);
                                bforce_x += bforce*pos[0]/coo;
                                bforce_y += bforce*pos[1]/coo;
                                bforce_z += bforce*pos[2]/coo;
                        }
                }
                end_particle_cell_loop(pi,c,t)
        }

    if(i==0) bforce = bforce_x;      
    else if(i==1) bforce = bforce_y;
    else
        bforce = bforce_z;
        /* an acceleration should be returned */
    return (bforce/P_MASS(p));
}
这是我写的求单个网格内颗粒间相互作用力的UDF 不知道有没有问题 请高人看一下 帮帮小弟吧

Conany

这个帖子咋没人顶嘞？

83136715

顶一下，我也想知道

他山有爱

我也用到           begin_particle_cell_loop(pi,c,t)
但是很容易出错。你解决了吗

newcfd123

他山有爱 发表于 2013-5-20 14:49
我也用到           begin_particle_cell_loop(pi,c,t)
但是很容易出错。你解决了吗

你好，请问解决了吗