VOF模型，密度UDF设置问题

ayeelan

请教下面的程序哪里有问题？模型是VOF模型，液态相水和气体air，想设置新的液态的密度，关于压力、温度、Xw的，现在interpreted之后，在materials里面选不中这个density01，先谢谢指教。

#include "udf.h"

DEFINE_PROPERTY(density01,cell,thread)
{
        real T=C_T(cell,thread);
        real P=C_P(cell,thread);
        real cw=C_YI(cell,thread,0);

        real R=8.314472;
        real Mh2o=18.;
        real Rw=1.;
        real Mba=70.1;
        real Rbaliq=741.8;
        real Rpu=1100.;

        real xw0=0.02;
        real cw0=0.02;
        real Xw=(cw0-cw)/cw0;
        real mba0=13.;
        real mpu0=260.;
        real Lliq0=mba0/mpu0;
        real Lgas=0.005;
       
        real density01;
        density01=(1+xw0+Lliq0)/((Xw*R*T*xw0/(P*Mh2o))+((1-Xw)*xw0/Rw)+
                (Lgas*R*T/(P*Mba))+((Lliq0-Lgas)/(Rbaliq))+(1/Rpu));
       
        return density01;
}

zzt_1234

建议你先看看一些例子，你这个程序肯定不行。你都没有F_PROFILE命令。

兰式折射

可以编译后试试

兰式折射

interpreted之后也可以选啊

zhangtao2030

浪子轻狂

你这个udf书写错误

ayeelan

real T=C_T(cell,thread);
        real P=C_P(cell,thread);
        real cw=C_YI(cell,thread,0);
这三行去掉就可以用，随便加哪一行都不能用，请教从网格中获取信息该怎么弄啊？

znon

不知道这个能不能管用，你可以试一下。

把 你的常量提前定义和赋值，然后给温度、压力和质量分数赋值。

ayeelan

谢谢你。

我把VOF模型去掉，就是算一个简单的流域，C_T就可以用了，但是VOF就不行，这个有什么细节设置么？还是程序里面要分清取的是哪个区域的温度？

invader

楼主，你骗人，下面的跟帖的人也骗人，
这个udf明明编译后就能选中，也能用，

200902303

你好，能请教你个问题嘛？

ayeelan

元老兄弟，，，我当时不懂啊，这个经过改动一下就可以用了，初学者总是碰到问题搞不定，有些还是元老们看来很简单的~~多多关照

ayeelan

另外，还有两个问题想请教一下
请教：VOF模型中想实现这样一个现象：两相中只有一相可以通过出口边界，而另一相不能流过去，这样的边界条件该如何实现呢？

还有就是怎么在UDF中写一个温度对时间的求导呢？       dT/dt     还有 dP/dt
谢谢斑竹~

mayo_0

还有就是怎么在UDF中写一个温度对时间的求导呢？       dT/dt     还有 dP/dt
谢谢斑竹~

这个FLUENT6.2 UDF的帮助文件里就有一个示例文件，还就是温度的。

[ 本帖最后由 mayo_0 于 2011-3-8 06:04 编辑 ]

ayeelan

谢谢
在帮助里还是没找到啊，请问下是这样写么？运行的时候提示错误，一步都不能计算。

real physical_dt=RP_Get_Real("physical-time-step");
real dTdt=(C_T(cell,thread)-C_T_M1(cell,thread))/physical_dt;


报的错：
Error:
FLUENT received fatal signal (ACCESS_VIOLATION)
1. Note exact events leading to error.
2. Save case/data under new name.
3. Exit program and restart to continue.
4. Report error to your distributor.
Error Object: ()