24小时热门版块排行榜

>论坛更新日志 (4131)
>博后之家 (204)
>导师招生 (201)
>论文投稿 (187)
>文献求助 (165)
>虫友互识 (145)
>有机交流 (80)
>硕博家园 (75)
>考博 (57)
>招聘信息布告栏 (43)
>学术会议 (27)
>休闲灌水 (27)
>找工作 (16)
>考研 (14)
>基金申请 (13)
>教师之家 (11)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xuxia8495

银虫 (小有名气)

应助: 16 (小学生)
金币: 403.6
帖子: 226
在线: 84.9小时
虫号: 1607257

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

使用DEFINE_INIT宏初始化两相体积，出现了问题：

如图所示：
我的模型是一个长方体，沿z轴方向将模型分成两部分，z向坐标范围为（-0.1,0.3）
0 初始化后，检查相分布，看第二相体积分数与定义的一致（图2），而查看第一相体积分布，发现整个区域全为第一相（图1）

UDF如下：
#include "udf.h"

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

thread_loop_c(t,d)
{
pri_th = THREAD_SUB_THREAD(t,0);
sec_th = THREAD_SUB_THREAD(t,1);
begin_c_loop_all (c,t)
{
C_CENTROID(xc,c,t);
if(xc[2]>=0&&xc[2]<=0.3)
{
C_VOF(c,sec_th)=1;
C_T(c,t)=1820;}
else
{
C_VOF(c,sec_th)=0;
C_T(c,t)=300;
}

}
end_c_loop_all(c,t)
}
}

请高手帮忙看看是哪出现问题了，谢谢！

图片1.png

图片2.png

回复此楼

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2012-11-06 17:32:40

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

bamboo_621

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 7.8
帖子: 1
在线: 75.1小时
虫号: 2138886

★ ★
xuxia8495(金币+1): 谢谢参与
xiegangmai: 金币+1, 谢谢参与！ 2013-07-10 23:03:12

引用回帖:

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

patch的话可以写个journal文件，每次计算之前导入journal文件就好了吧

赞一下

回复此楼

10楼2013-07-10 15:31:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 14 个回答

fangdian

木虫 (正式写手)

应助: 201 (大学生)
金币: 2836.4
帖子: 818
在线: 206.9小时
虫号: 1926821

★ ★
xuxia8495(金币+1): 谢谢参与
xiegangmai: 金币+1, 谢谢参与 2012-11-12 19:36:57

实在不行的话，对基本相也初始化吧。

赞一下

回复此楼

2楼2012-11-06 20:06:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

转石千仞

金虫 (正式写手)

应助: 24 (小学生)
金币: 941
帖子: 673
在线: 161.6小时
虫号: 1163587

★ ★ ★
xuxia8495(金币+1): 谢谢参与
xiegangmai: 金币+2, 谢谢参与 2012-11-07 20:19:01

这是我写的，你看看罢。

CODE:

#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)

        }

}