24小时热门版块排行榜    

查看: 3829  |  回复: 14

malu430

金虫 (小有名气)

[求助] vof模型算气液两相逆流 使用udf加动量源项 碰到问题 请高手解答

在用vof模型计算气液两相逆流
想在气液界面加上气液剪切力动量源项
现在可以找到气液界面
并且也可以找到离界面3--4个网格大小的网格的坐标
但是怎么通过这个坐标反过来找到网格,并且获取网格中的一些信息呢,比如速度之类
感谢大侠帮忙
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

malu430

金虫 (小有名气)

xiegangmai: 回帖置顶 2012-10-26 19:59:51
引用回帖:
10楼: Originally posted by luowenxu at 2012-09-01 14:05:50
你的界面是取体积分数0.5得到的吗?
你可以使用export输出界面物理量

嗯 是取0.5得到的
这个我试过了 哈哈 的确可以输出
11楼2012-09-01 16:57:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

fangdian

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
malu430: 金币+2 2012-09-01 12:02:45
能不能通过循环去找:如果网格点的中心点坐标和你找到的那个网格坐标基本相近,就采用那个网格。比方说类似代码:
thread_loop_c(t,d)
{
begin_c_loop(c,t)
{C_CENTROID(x,c,t);
if(fabs(x[0]-你找到的x坐标)<=1e-6&&fabs(x[0]-你找到的y坐标)<=1e-6&&fabs(x[0]-你找到的z坐标)<=1e-6)
{加入输出单元值的代码.......
    }
}
end_c_loop(c,t)
}
6楼2012-09-01 08:55:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangdian

木虫 (正式写手)

上面代码中if写错,应该为:if(fabs(x[0]-你找到的x坐标)<=1e-6&&fabs(x[1]-你找到的y坐标)<=1e-6&&fabs(x[2]-你找到的z坐标)<=1e-6)
7楼2012-09-01 08:56:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

gyctju

金虫 (著名写手)

很难,几乎是mission impossible
2楼2012-08-30 23:53:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ya634984094

铜虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
malu430: 金币+2 2012-08-31 09:54:37
是不是要加表面张力啊??你用的神马模型?再看怎样加力啊~
3楼2012-08-31 00:19:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

malu430

金虫 (小有名气)

引用回帖:
2楼: Originally posted by gyctju at 2012-08-30 23:53:43
很难,几乎是mission impossible

啊 这样啊 真的吗?
4楼2012-08-31 09:54:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

malu430

金虫 (小有名气)

引用回帖:
3楼: Originally posted by ya634984094 at 2012-08-31 00:19:09
是不是要加表面张力啊??你用的神马模型?再看怎样加力啊~

不是表面张力 是要加气液剪切力
我用的vof模型啊
看了很多文献 都是写个udf导进去就可以实现
谢谢
5楼2012-08-31 09:56:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

malu430

金虫 (小有名气)

引用回帖:
6楼: Originally posted by fangdian at 2012-09-01 08:55:41
能不能通过循环去找:如果网格点的中心点坐标和你找到的那个网格坐标基本相近,就采用那个网格。比方说类似代码:
thread_loop_c(t,d)
{
begin_c_loop(c,t)
{C_CENTROID(x,c,t);
if(fabs(x-你找到的x坐标)< ...

恩。这个办法不错
会不会计算量有点大。。。
8楼2012-09-01 12:03:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

malu430

金虫 (小有名气)

引用回帖:
7楼: Originally posted by fangdian at 2012-09-01 08:56:49
上面代码中if写错,应该为:if(fabs(x-你找到的x坐标)<=1e-6&&fabs(x-你找到的y坐标)<=1e-6&&fabs(x-你找到的z坐标)<=1e-6)

谢谢
9楼2012-09-01 12:03:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luowenxu

木虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
malu430: 金币+2 2012-09-01 16:56:47
你的界面是取体积分数0.5得到的吗?
你可以使用export输出界面物理量

[ 发自手机版 http://muchong.com/3g ]
若想人前显贵,必先人后受罪
10楼2012-09-01 14:05:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 malu430 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 13/650 2026-08-04 07:30 by k7OM8YghWbkC
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +9 锦衣卫寒战 2026-07-28 13/650 2026-08-04 07:21 by Ermito
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jRl3mE6ddZGq 2026-08-03 8/400 2026-08-04 06:15 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 06:07 by k7OM8YghWbkC
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jRl3mE6ddZGq 2026-08-03 4/200 2026-08-04 04:54 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:46 by k7OM8YghWbkC
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 9/450 2026-08-04 03:56 by kf2781974
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见