请教：在DPM模型中采用UDF加入电磁力后运行出错

yf1978

我的计算网格模型是一个二维轴对称模型。加入UDF后计算几步后总是出现如下错误：
Error:
FLUENT received fatal signal (ACCESS_VIOLATION)
1. Note exact events leading to error.
2. Save case/data under new name.
3. Exit program and restart to continue.
4. Report error to your distributor.
Error Object: ()
UDF程序如下：
#include "udf.h"
#define Q 1.0         /* particle electric charge      */
#define BZ 1.0
DEFINE_DPM_BODY_FORCE(particle_body_force, p, i)
{
        Thread *t;
        cell_t c;
        Domain *d;
        real bforce;
        real xc[ND_ND];
        real x, y;
/* loop over all cell threads in the domain  */
  thread_loop_c (t,d)
  {
/* loop over all cells  */
   begin_c_loop_all (c,t)
    {
          C_CENTROID(xc,c,t);
          x=xc[0];
  y=xc[1];
        if((y<=0.03) & (x<=0.35) & (x>=0.15))
        {
           if(i==0) bforce=Q*BZ*P_VEL(p)[1];      
           else if(i==1) bforce=-Q*BZ*P_VEL(p)[0];
         
        }
        else
           bforce=0.0;
        /* an acceleration should be returned */
        return (bforce/P_MASS(p));
    }
   end_c_loop_all (c,t)
  }
请各位多多指教，不胜感谢！

a2883524

不知道是不是你拷贝错误，就你上面程序，后面少一个大括号。第二个就是一般return语句在最后的大括号前面，不是在循环里面，我想放在里面可能会导致循环出错。

yf1978

谢谢斑竹。
少一个括号是拷贝的问题。
我把程序改了一下，还是不行，收敛10步左右后出现错误提示。请再帮我看看原因是什么？错误信息还是
Error:
FLUENT received fatal signal (ACCESS_VIOLATION)
1. Note exact events leading to error.
2. Save case/data under new name.
3. Exit program and restart to continue.
4. Report error to your distributor.
Error Object: ()
程序如下：
DEFINE_DPM_BODY_FORCE(particle_body_force, p, i)
{
        Thread *t;
        cell_t c;
        Domain *d;
        real bforce;
        real xc[ND_ND];
        real x, y;
          C_CENTROID(xc,c,t);
          x=xc[0];
  y=xc[1];
        if((x<=0.35) & (x>=0.15))
        {
           if(i==0) bforce=Q*BZ*P_VEL(p)[1];      
           else if(i==1) bforce=-Q*BZ*P_VEL(p)[0];
         
        }
        else
           bforce=0.0;
        /* an acceleration should be returned */
        return (bforce/P_MASS(p));
}

yf1978

另外，我想请教一下如何能得到离散相颗粒的坐标值？我使用了提取离散相颗粒质量和速度值的宏，如上。
我查了一下帮助文档，里面有很多宏可以得到离散相的各种参数，唯独没有如何提取离散相颗粒的坐标的宏。

hq111

我也是出现这样的问题呢，不知怎么解决的？呵呵，我的qq号是715176266，你若是解决的请指点！谢谢

octen-ol

我也想知道

LHZ

离散相颗粒坐标，P_POS(p)