版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

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

回复此楼

» 收录本帖的淘帖专辑推荐

fluent点滴学习

fluent软件学习

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

udf中的面宏如何使用，例如F_UDMI(f,t,i),和C_UDMI(c,t,i)的使用有什么区别？已经有0人回复
FLUENT UDF求助已经有2人回复
UDF宏DEFINE_DPM_INJECTION_INIT不能正常喷淋已经有2人回复
关于Lookup_Thread（domain，id）宏求教各位前辈已经有2人回复
关于C_YI_G（c,t，i）和F_AREA(A,f,t)的使用已经有1人回复
fluent udf帮助文件里面可压缩液体的计算公式源自哪里？已经有6人回复
DPM使用UDF时提示 P_POS: undeclared variable 已经有3人回复
汽蚀（空化）UDF编写中遇到的问题已经有13人回复
[已完结]UDF计算species 模型出错已经有3人回复
壁面热边界条件中辐射边界的问题已经有11人回复
《跟我一起学C++》目录及视频下载地址已经有116人回复
如何用Compute_Force_And_Moment函数求流体对运动翼型的俯仰力矩已经有12人回复
[已完结]使用fprintf显示不了需要的变量值已经有17人回复
ANYSY UDF使用問題已经有1人回复
导入表面张力系数的udf文件后报错已经有3人回复
AMBER使用antechamber 准备配体小分子出问题，求指导方法，谢谢已经有9人回复
udf编程通过，但是在FLUENT中运行出错，麻烦udf高手帮忙看看已经有4人回复
[Free]C语言深度解剖－解开程序员面试笔试的秘密.陈正冲.2008 已经有36人回复

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

查看全部散金贴

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

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

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)

        }

}

赞一下

回复此楼

» 本帖已获得的红花（最新10朵）

xuxia8495

3楼2012-11-07 08:39:45

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

转石千仞

金虫 (正式写手)

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

★ ★
xiegangmai: 金币+2, 辛苦了！ 2012-11-07 20:19:10

你的问题好像是：在混合域中循环多相线的问题，这个要用到特殊的宏，好像是：
mp_thread_loop_c
具体的你要自己去看了。
你可以画网格时定义区域，再patch值，也非常方便，后处理也好做些。

赞一下

回复此楼

4楼2012-11-07 08:44:43

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

xuxia8495

银虫 (小有名气)

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

送鲜花一朵

引用回帖:

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);
...

非常谢谢！
我运行了一下这种写法，发现问题并没有解决：检查相分布，查看基本相体积分数，与定义的一致，而查看第二相体积分布，发现整个区域全为第二相-------看来需要分别定义基本相和第二项

赞一下

回复此楼

5楼2012-11-08 09:37:24

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

xuxia8495

银虫 (小有名气)

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

#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,pri_th)=0;
C_T(c,t)=1820;}
else
{
C_VOF(c,pri_th)=1;
C_T(c,t)=300;
}

}
end_c_loop_all(c,t)
}
}

赞一下

回复此楼

6楼2012-11-08 09:46:38

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

xuxia8495

银虫 (小有名气)

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

★ ★
xiegangmai: 金币+2, 鼓励讨论交流 2012-11-08 19:39:00

这样修改一下就可以了

#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,pri_th)=0;
C_T(c,t)=1820;}
else
{
C_VOF(c,pri_th)=1;
C_T(c,t)=300;
}

}
end_c_loop_all(c,t)
}
}

赞一下

回复此楼

7楼2012-11-08 09:47:20

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

18662515085

银虫 (正式写手)

应助: 142 (高中生)
金币: 419.8
帖子: 487
在线: 220.4小时
虫号: 1955839

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

建议你在划分网格的时候定义好区域
然后再FLUENT里面初始化以后，直接PATCH
没必要使用UDF定义
QQ：1442083080

赞一下

回复此楼

8楼2012-11-09 07:34:51

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

xuxia8495

银虫 (小有名气)

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

引用回帖:

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

谢谢啊，感觉patch没有用UDF方便。。。

赞一下

回复此楼

9楼2012-11-09 08:38:12

已阅回复此楼关注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的回帖

meiguimimi

银虫 (小有名气)

应助: 1 (幼儿园)
金币: 338.8
帖子: 53
在线: 51.3小时
虫号: 2807908

★
xuxia8495(金币+1): 谢谢参与

引用回帖:

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

请问，journal文件需要怎么编写呢？可否提供一个参考方向？

赞一下

回复此楼

11楼2015-01-16 16:54:41

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

死天

专家顾问 (著名写手)

专家经验: +60
应助: 399 (硕士)
金币: 3648.9
帖子: 1336
在线: 443.8小时
虫号: 2643057

★
xuxia8495(金币+1): 谢谢参与

建议用patch ，不要用init，因为用UDF是和初始化一起工作，而Patch相当于算一步。只初始化的话，体积分数没有算，即各个相的体积分数之和为1，这步没有算。所以建议用patch

赞一下

回复此楼

12楼2015-01-16 17:09:19

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

DGE旅行者

新虫 (小有名气)

应助: 3 (幼儿园)
金币: 370.8
帖子: 150
在线: 47.3小时
虫号: 3604073

★
xuxia8495(金币+1): 谢谢参与

引用回帖:

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 ...

这个还是不行吧

回复此楼

13楼2015-01-21 10:32:19

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

KEENWANG

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 14.8
帖子: 3
在线: 5.8小时
虫号: 8327991

★
xuxia8495(金币+1): 谢谢参与

楼主解决了吗？可以分享一下结果吗？

赞一下

回复此楼

14楼2018-12-11 17:04:56

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

相关版块跳转我要订阅楼主 xuxia8495 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 一志愿985，本科211，0817化学工程与技术319求调剂 +5	Liwangman 2026-03-15	5/250	2026-03-16 17:10 by 我的船我的海
[考研] 0703一志愿211 285分求调剂 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 308求调剂 +3	是Lupa啊 2026-03-16	3/150	2026-03-16 10:07 by 求调剂zz
[教师之家] 焦虑 +7	水冰月月野兔 2026-03-13	9/450	2026-03-16 10:00 by Quakerbird
[考研] 290求调剂 +5	孔志浩 2026-03-12	10/500	2026-03-16 09:01 by 余晖&
[考研] 326求调剂 +4	上岸的小葡 2026-03-15	5/250	2026-03-16 08:39 by Linda Hu
[考研] 机械专硕调剂 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 266求调剂 +4	学员97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 307求调剂 +7	超级伊昂大王 2026-03-10	7/350	2026-03-14 00:49 by JourneyLucky
[考研] 材料工程,326分，求调剂 +6	KRSLSR 2026-03-10	6/300	2026-03-13 23:47 by JourneyLucky
[考研] 材料371求调剂 +9	鳄鱼? 2026-03-11	11/550	2026-03-13 22:53 by JourneyLucky
[考研] 求调剂（材料与化工327） +4	爱吃香菜啦 2026-03-11	4/200	2026-03-13 22:11 by JourneyLucky
[考研] 293求调剂 +3	世界首富 2026-03-11	3/150	2026-03-13 16:27 by JourneyLucky
[考研] 【0856】化学工程（085602）313 分，本科学科评估A类院校化学工程与工艺，诚求调剂 +7	小刘快快上岸 2026-03-11	7/350	2026-03-13 16:06 by ruiyingmiao
[考研] 307求调剂 +5	超级伊昂大王 2026-03-12	5/250	2026-03-13 15:56 by 棒棒球手
[考研] 308求调剂 +3	是Lupa啊 2026-03-12	3/150	2026-03-13 14:30 by 求调剂zz
[考研] 0817化学工程与技术考研312分调剂 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu
[考博] 读博申请 +5	感dd 2026-03-10	7/350	2026-03-11 17:02 by QGZDSYS
[考研] 化工0817调剂 +8	灿若星晨 2026-03-10	8/400	2026-03-10 22:44 by 星空星月
[考研] 085602化工求调剂 +7	董boxing 2026-03-10	7/350	2026-03-10 17:07 by BruceLiu320

24小时热门版块排行榜

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

» 收录本帖的淘帖专辑推荐

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到: 查看全部散金贴

» 本帖已获得的红花（最新10朵）

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

查看全部散金贴