请教各位大神，有什么命令能把DPM里面的incomplete值输出

fty0083

我现在做dpm，发现有很多颗粒的FATE是incomplete，看过了帮助，我理解这些颗粒基本上应该是悬浮了。但是考虑到数值收敛性，需要检测一下这个值。我看过了dpm.h的内容，基本上incomplete、trap、abort这些都是一个叫Particle_FATE的值。所以我想可能会有什么命令能返回这个东西，我能把它记录到文件里面去。

下面内容引自文件dpm.h：

typedef enum
{
  FATE_ABORTED=0,
  FATE_INCOMPLETE,
  FATE_EVAPORATED,
  FATE_ESCAPED,    /* also includes trapped */
  FATE_MAX
} Particle_Fate;

另外，这里还有一个结构体：

struct global_tracked_particle_vars_struct
{
**********************
float intFact[MAX_CELL_NODES];

  int n_tracked;
  int n_escaped;
  int n_aborted;
  int n_trapped;
  int n_evaporated;
  int n_incomplete;
  int n_coalesced;
  int n_splashed;
  int n_shed;
}
我估计，fluent自己生成的报告可能就是从这些数里面输出的，但是，我想把这些数记录下来。

[ 本帖最后由 fty0083 于 2011-9-14 09:12 编辑 ]

fty0083

呵呵，今天自己研究了一下，问题已经解决了。下面说一下解决办法和心得。
这个数可以用下面的命令返回
p->gvtp.n_incomplete
简单解释一下，p是Tracked_Particle类型，所以这个命令必须在带有*tp的宏里面用。另外，这个语句是从指针p的结构体里面指向了所引用的global_tracked_particle_vars_struct类，原结构体定义了这个类被称为gvtp，这个类是全局追中参数类，里面有n_incomplete，所以用结构体表示法p->gvtp.n_incomplete指向这个值。值得提的是，fluent的DPM提供P_VEL(p)表示当前颗粒的速度，也可以用p->state.V，这个值的引用可以和p->gvtp.n_incomplete一起理解。可见，fluent能返回的信息有很多，实在找不到，就去找头文件吧，好用！

yiteng

原帖由 fty0083 于 2011-9-14 18:57 发表

                               
登录/注册后可看大图

呵呵，今天自己研究了一下，问题已经解决了。下面说一下解决办法和心得。
这个数可以用下面的命令返回
p->gvtp.n_incomplete
简单解释一下，p是Tracked_Particle类型，所以这个命令必须在带有*tp的宏里面用。另外 ...

楼主很强。

小博

不知道楼主现在还有没有用这个论坛，我现在在做气液固多相流下管道的颗粒沉积模拟研究，有些问题想问问楼主，能不能QQ交流下呢？