24小时热门版块排行榜    

查看: 3884  |  回复: 14

malu430

金虫 (小有名气)

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

在用vof模型计算气液两相逆流
想在气液界面加上气液剪切力动量源项
现在可以找到气液界面
并且也可以找到离界面3--4个网格大小的网格的坐标
但是怎么通过这个坐标反过来找到网格,并且获取网格中的一些信息呢,比如速度之类
感谢大侠帮忙
回复此楼
已阅   回复此楼   关注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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +6 scalable 2026-08-24 6/300 2026-08-24 20:28 by 苦寒来香
[基金申请] 我面上完蛋了 +9 且听虎啸 2026-08-20 10/500 2026-08-24 20:08 by 苦寒来香
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 没有任何消息-是不是就凉了 +6 图啦图啦 2026-08-24 6/300 2026-08-24 17:44 by 小小红帆
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 明天应该可查了!? +4 chengyan1220 2026-08-23 4/200 2026-08-24 08:47 by gloomy6159
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见