VOF-UDF-UDS-Floating point error 问题

guohf00001

我用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文件，请高手帮我看一下，出错在什么地方，我刚接触这方面的问题，是不是对问题的理解有误。

先谢谢各位。

jialixujia

我也老中碰到这种问题。。。。哎。。情况不一样。。。

znon

没细看，
其中，对于边界，Axial应该是X-轴方向。
Apply an axis boundary on the centerline (x-axis).
For geometry that is axisymmetric or axisymmetric swirl , the centerline (x-axis) boundary type should be set to axis. See Section  7.3.17.