24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3356  |  回复: 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 的主题更新
信息提示
请填处理意见