24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 5471  |  回复: 17
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zcjwyl

金虫 (小有名气)


[交流] udf 重力方向变化

Hi, 各位虫友。在玩一个典型问题for fun:2D 矩形容器内部一半water一半air,容器绕中心轴线左右晃动,角速度和摆幅已知道,VOF计算自由液面。

现在用动网格可以实现容器的俯仰,不过计算消耗大了些。好点的办法是假设容器是不动的,而让重力的方向做相应变化,就不涉及到动网格了。

现在问题是Fluent里面似乎没有什么地方可以编译重力,operation condition panel里面只能定义一个定常的重力。看到有文章立报道过fluent里变重力的,不过作者没有公布详细的做法。

哪位大侠有做个类似问题的,请指教。
回复此楼

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

fluent点滴学习

» 猜你喜欢

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

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

查看全部散金贴

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

丹丹丫头

新虫 (初入文坛)



zcjwyl(金币+1): 谢谢参与
引用回帖:
1120833楼: Originally posted by NaOHWTO at 2011-06-14 20:28:34:
你可以仔细看看动量方程,你在FLUENT面板里面设重力为0,用UDF在源项里面加入 密度*(-9.8m/s), 放到动量方程里面效果是一样的。

请教一个问题,如果添加一个倾斜圆管的重力源项,且密度是按boussinesq函数变化的,UDF该如何添加呢,我自己写了一个UDF,但是计算结果感觉不太对,能帮忙看一下吗,三角函数这样写对吗
#include "udf.h"
#define g 9.8
#define a 30
#define b 1000
#define d 293.15
DEFINE_SOURCE(ymom_source, c, t, dS, eqn)
{
real m,con, source;
m=1/d;
con=m*b*g*cos(a);
source= con *(C_T(c, t)-d);
dS[eqn] = con;
return source;
}
13楼2012-04-30 14:25:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

heattransfet

新虫 (职业作家)


★ ★
zcjwyl(金币+1):谢谢参与
臭水沟(金币+1): 谢谢应助~~ 2011-06-13 23:11:48
用UDS编写新的方程
2楼2011-06-13 18:18:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NaOHWTO

金虫 (小有名气)


★ ★ ★
zcjwyl(金币+1):谢谢参与
臭水沟(金币+2): 谢谢应助~~ 2011-06-13 23:12:09
要用UDF哦,
附件算例,是汽车油箱直线60m/s跑, 突然刹车,2s后停下来的油箱晃动。

[ Last edited by NaOHWTO on 2011-6-13 at 21:52 ]
3楼2011-06-13 20:53:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NaOHWTO

金虫 (小有名气)


★ ★ ★
臭水沟(金币+3): 给力!谢谢应助,希望常来交流! 2011-06-13 23:13:21
zcjwyl(金币+3): 谢谢你的算例,不过直接用Adjust宏可以修改重力吗,Fluent代码里重力的定义是怎么写的?M_gravity。我现在的做法只是直接修改动量源项,不过还是麻烦点,最好直接编辑重力。 2011-06-14 10:47:03
这个我自己编了个6自由度的UDF程序,不过这些肯定不会给出来的。
中文文献上的那几篇用FLUENT日志文件+UDF算的有点过于麻烦。总之,
告诉你只要用UDF的源项宏和ADJUST宏就能完成,你想要怎么晃动都行!
附件是一个6自由度的算例[ Last edited by NaOHWTO on 2011-6-13 at 21:50 ]
4楼2011-06-13 21:44:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见