UDF 出错求助

longzi83

用的non-premixed model模拟雾化燃烧，修改了mixture fraction variance的源项。程序比较简单，不用compile看似没问题，但是迭代时只要有液滴喷入计算区域就会出现 fatal error。用的ANSYS Fluent12.1。请教各位代码是否有什么问题：
#include "udf.h"
#include "dpm.h"
#include "surf.h"  /* for  macros: RP_Cell() & RP_Thread()        */
#include "prop.h"  /* for function: Saturation_Pressure() (of water) */

DEFINE_SOURCE(Var1_source,c,t,dS,eqn)
{
        real source;
          
  if (C_STORAGE_R(c,t,SV_DPMS_MASS) > 0.)
  {
Message("dpm=%f",C_DPMS_PDF_1(c,t));
   }
source = 0.;
dS[eqn] = 0.;
  
  return source;
}

[ 本帖最后由 longzi83 于 2011-2-2 14:05 编辑 ]

invader

udf正确，没检查出问题

longzi83

谢谢楼上。我compile调试一切也看上去都正常，就是一用到source term上去，马上报错。调试了很久，发现如果单独message 代码中的其他各值都正确，唯独C_DPMS_PDF_1(c,t)，即颗粒的蒸发源项一旦message就报错。分析可能错误就是出在这个项上。请问这个参数是否正确？是否有其他办法求得蒸发源项？谢谢

[ 本帖最后由 longzi83 于 2011-1-25 15:44 编辑 ]

invader

C_DPMS_PDF_1(c,t)，
不好意思，我查找了半天都不知道这个宏在帮助文件的什么位置，

意外发现，fluent自带的帮助文档，非常不全面，在网上都有4.6.2，而我的确只有4.4
不知所谓了

longzi83

你好，我在帮助和dpm.h里面也都没有看到这个项，这是以前一个人在老版本中可以调用的东西，不知道为什么现在用不了了。但是如果把它的名称改一下就会报错，貌似定义仍然在，就是用不了。它的定义是SUM(nm)。n 表示液滴个数，m表示质量蒸发率。不知道有没有其他变通的办法定义它。

[ 本帖最后由 longzi83 于 2011-1-26 13:46 编辑 ]

invader

不是版本的事儿，定义肯定还存在，
我看帮助的目的是想找到这个宏的工作机制及其目的，这样的话，帮助修改程序，
不了解程序代码的目的当然看起来就是正确，实际却不知是什么了，呵呵

longzi83

核对了一下，12.0的dpm.h有 64k，比 12.1大多了。各版本sg_mem.h都能看到C_DPMS_PDF_1(c,t)，就是不知道怎么调用的。

invader

呵呵，
我对dpm何pdf都不懂啊，
现在想知道你的程序里，使用了上述函数，那么这个函数在程序里代表的是哪个物理量呢？

invader

DEFINE_DPM_SOURCE
我看帮助文件中，有这个宏，为何你的udf不用这个宏呢，不也是使用了dpm模型的吗？

longzi83

我的意图是要修改flamelet model里的mixture fraction variance源项。其中C_DPMS_PDF_1(c,t)是想调用这个cell中颗粒蒸发的源项，以计算mixture fraction variance源项。尝试了不同的版本，12.0.7可以调用这个参数（可能架构和6.3还差不多，所以能用），但是其他就出错了。好像后面版本PDF的格式不一样了。不知道程序怎么样使用DEFINE_DPM_SOURCE来计算最终的mixture fraction variance源项？谢谢。