udf of the "particle-body-force" for help

zhaoh

Dear all:

My udf of "particle-body-force" is as follows: #include "udf.h" #include "dpm.h" #include "surf.h" #define C1 2.2 #define C2 1.5 DEFINE_DPM_BODY_FORCE(p_b_force,p,i) { real bforce; if(i==1) bforce=0; else if(i==0) bforce=C1*P_POS(p)[0]+C2*P_POS(p)[1]; return (bforce/P_MASS(p)); } When i compiled it and iterate it, there is a mistake: divergence detected - temporarily reducing Courant number to 0.5 and trying again... divergence detected - temporarily reducing Courant number to 0.05 and trying again... divergence detected - temporarily reducing Courant number to 0.005 and trying again... divergence detected - temporarily reducing Courant number to 0.0005 and trying again... divergence detected - temporarily reducing Courant number to 5e-05 and trying again... Error: Internal error at line 872 in file 'rp_mstage.c'. divergence detected in AMG solver Error Object: ()

and there are often a mistake as follows: Can't figure which cell particle is going to -- ABORT Can't figure which cell particle is going to -- ABORT Can't figure which cell particle is going to -- ABORT Can't figure which cell particle is going to -- ABORT Can't figure which cell particle is going to -- ABORT Can't figure which cell particle is going to -- ABORT Can't figure which cell particle is going to -- ABORT Can't figure which cell particle is going to -- ABORT Can't figure which cell particle is going to -- ABORT\

I don`t know how to solve these questions. Who can help me? Thanks alot! Expecting for your reply!