[交流] 使用DEFINE_INIT宏初始化两相体积

#include "udf.h"   /*初始化体积分数与温度*/
DEFINE_INIT(my_init_function,domain)
{
Domain *d=Get_Domain(1);
Domain *sd=Get_Domain(2);
Domain *th=Get_Domain(3);
Thread *t;
cell_t c;
real xc[ND_ND];
thread_loop_c(t,d)
        {
                begin_c_loop(c,t)
                {
                        C_CENTROID(xc,c,t);
                        if(xc[1]<0.5)
                        C_T(c,t)=400.;
                        else
                        C_T(c,t)=300.;
                }
                end_c_loop(c,thread)
        }
thread_loop_c(t,sd)
        {
                begin_c_loop(c,t)
                {
                        C_CENTROID(xc,c,t);
                        if(xc[1]<1.5)
                        C_VOF(c,t)=0.;
                        else
                        C_VOF(c,t)=1.;
                }
                end_c_loop(c,thread)
        }
}

3楼: Originally posted by 转石千仞 at 2012-11-07 08:39:45
这是我写的，你看看罢。
#include "udf.h"   /*初始化体积分数与温度*/
DEFINE_INIT(my_init_function,domain)
{
Domain *d=Get_Domain(1);
Domain *sd=Get_Domain(2);
Domain *th=Get_Domain(3);
...

8楼: Originally posted by 18662515085 at 2012-11-09 07:34:51
建议你在划分网格的时候定义好区域
然后再FLUENT里面初始化以后，直接PATCH
没必要使用UDF定义
QQ：1442083080

9楼: Originally posted by xuxia8495 at 2012-11-09 08:38:12
谢谢啊，感觉patch没有用UDF方便。。。...

10楼: Originally posted by bamboo_621 at 2013-07-10 15:31:41
patch的话可以写个journal文件，每次计算之前导入journal文件就好了吧...

7楼: Originally posted by xuxia8495 at 2012-11-08 09:47:20
这样修改一下就可以了

#include "udf.h"

DEFINE_INIT(my_init,d)
{
cell_t c;
Thread *t;
real xc;
Thread *pri_th,*sec_th;

thread_loop_c(t,d)
{
pri_th = THREAD_SUB_THREAD(t,0);
s ...