如图所示,在上面横线上释放固体颗粒,计算颗粒在圆环外表面上的交点和颗粒的流动轨迹。
根据帮助文件里的DPM output输出一直有错误。
是否需要修改里面参数,我是直接用的。
如何求颗粒轨迹与外表面的交点?
#include "udf.h"
DEFINE_DPM_OUTPUT(dpm_output,header,fp,p,thread,plane)
{
char name[100];
if (header)
{
if (NNULLP(thread))
fprintf(fp,"(%s\t)\n",thread->head->dpm_summary.sort_file_name);
else
fprintf(fp,"(%s\t)\n",plane->sort_file_name);
fprintf(fp,"%10s\t%10s\t%10s\t%10s\t%10s\t%10s\t%10s\t%10s\t\n", "dt","diameter","X","Y","Z","U","V","W";
}
else
{
sprintf(name,"%s:%d",p->injection->name,p->part_id);
fprintf(fp, "%10.6g\t%10.6g\t%10.6g\t%10.6g\t%10.6g\t%10.6g\t%10.6g\t%10.6g\t\n",RP_Get_Real("flow-time",P_DIAM(p),P_POS(p)[0],P_POS(p)[1],P_POS(p)[2],P_VEL(p)[0],P_VEL(p)[1],P_VEL(p)[2]);
}
}
|