24小时热门版块排行榜

返回列表

yanjiechao

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3.5
帖子: 17
在线: 19.4小时
虫号: 3836534
注册: 2015-04-27
性别: GG
专业: 工程海洋学

[求助] 用vof 模型，设定第二相（water）的体积分数，需编udf，编译以后出现错误！必重谢

我的研究是多相流，一个水槽设置一个新型齿坝，下边9m是水，上边4.5m是空气，想要用vof 模型捕抓液面，但要设定第二相（water）的体积分数，就需要编写一个udf，编译以后出现错误，希望来这里得到救助！必重谢
#include "udf.h"
#include "sg_mphase.h"
DEFINE_INIT(my_init_function,mixture_domain)
{
int phase_domain_index;
cell_t cell;
Thread *cell_thread;
Domain *subdomain;
real  xc[ND_ND];
   real x;
real y;
real z;
sub_domain_loop(subdomain,mixture_domain,phase_domain_index)
{
if (DOMAIN_ID(subdomain)==2)
   thread_loop_c (cell_thread,subdomain)
{
   begin_c_loop_all (cell,cell_thread)
{
C_CENTROID(xc,cell,cell_thread);
                        x = xc[0];
                        y = xc[1];
                        z = xc[2];
    if ((-37.5.<xc[0]<52.5.)&&(-3.1.<xc[1]<3.1.)&& (0.<xc[2]<9.))
C_VOF(cell,cell_thread)=1.;
else
C_VOF(cell,cell_thread)=0.;
}
end_c_loop_all (cell,cell_thread)
}
}
}

回复此楼