24小时热门版块排行榜    

查看: 3877  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 能否退出参与的面上项目解除限项 +15 koalala 2026-08-24 17/850 2026-08-24 14:16 by hhs666
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 18/900 2026-08-24 12:03 by iaeyuan
[基金申请] 范进中举一文的中心思想 +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
[基金申请] 什么时候开奖? +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
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +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
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见