VOF-UDF-UDS-Floating point error 问题

guohf

我用VOF模式模拟液体在高电压下喷射到空气中的流动，想得到泰勒锥，x轴为水平轴， y轴为垂直向下。用水做第二相，空气为第一相，入口速度为抛物形（其值很小）， 其UDF如下：
#include "udf.h"
DEFINE_PROFILE(velocity_inlet,thread,position)
{
real x[ND_ND];
real y;
face_t f;
begin_f_loop(f,thread)
    {
      F_CENTROID(x,f,thread);
      y=x[0];
      F_PROFILE(f,thread,position)=0.00001*(3.0-1.5*y*y/8.0);
    }
  end_f_loop(f,thread)
}

电压是14kv,电场满足 div(grad(V))=0, 其中V为电压。

模式选用定常、层流。条件全部设好后，进行迭代时，不能运行，出错提示：
Floating point error: invalid number

附件是我的cas文件，请高手帮我看一下，出错在什么地方，我刚接触这方面的问题，是不是对问题的理解有误。

先谢谢各位。

guohf

我用定值速度入口也试过，还是一样的问题。

我的工程背景就是高压下的溶液射流，主要想看一下在高电场作用下射流规迹或气液界面的变化。

昨天没把cas文件传上去，今天补上了，请各位帮忙。做了很长时间了点头绪都没有，很急。

谢谢

注：因cas文件太大，压缩了，上传时将文件扩展名改为了.cas，请下载后将文件扩展名改回 ".rar"

cone.cas (42.37 KB, 下载次数: 2)

2010-5-20 11:37 上传

点击文件名下载附件