应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » udf 重力方向变化
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
181/1返回列表
查看: 5468  |  回复: 17
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zcjwyl

金虫 (小有名气)

[交流] udf 重力方向变化

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

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

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

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

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

fluent点滴学习

» 猜你喜欢

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

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

查看全部散金贴
1楼2011-06-13 16:09:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

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 ]
赞一下(8人) 回复此楼
4楼2011-06-13 21:44:52
 已阅   回复此楼   关注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 ]
赞一下(7人) 回复此楼
3楼2011-06-13 20:53:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

heattransfet

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

47425804

木虫 (正式写手)
★
zcjwyl(金币+1):谢谢参与
zcjwyl(金币+2): 你说的很对，谢谢 2011-06-14 10:49:18
重力其实就是个源项，可以用源项UDF实现重力的变化，上面虫友提供的算例非常好，和你的问题实质上是完全一样的。
赞一下(2人) 回复此楼
5楼2011-06-13 23:55:56
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NaOHWTO

金虫 (小有名气)
★ ★
xiegangmai(金币+2): 谢谢应助！ 2011-06-15 19:56:07
zcjwyl(金币+2): 明白了，谢谢。 2011-06-17 10:21:47
重力加速度只是一个方向的线加速度，例如Y向的，所以UDF中用源项加入到Y的动量方程中去，不过，ADJUST宏只是为了调整你给定激励的时间，例如，我的第一个算例，2s后刹车停住了，就没有加速度了，所以这个时候就不要源项了。你不可能去手动改算到2s后停止计算，去掉源项后再算，那样就很麻烦了。
当然，你要一直有加速度，那就只用源项吧。
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

47425804
6楼2011-06-14 20:24:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NaOHWTO

金虫 (小有名气)
引用回帖:
Originally posted by NaOHWTO at 2011-06-13 21:44:52:
这个我自己编了个6自由度的UDF程序，不过这些肯定不会给出来的。
中文文献上的那几篇用FLUENT日志文件+UDF算的有点过于麻烦。总之，
告诉你只要用UDF的源项宏和ADJUST宏就能完成，你想要怎么晃动都行！
附件是 ...

你可以仔细看看动量方程，你在FLUENT面板里面设重力为0，用UDF在源项里面加入 密度*（-9.8m/s), 放到动量方程里面效果是一样的。
赞一下 回复此楼
7楼2011-06-14 20:28:34
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

47425804

木虫 (正式写手)
★
送鲜花一朵
xiegangmai(金币+1): 鼓励讨论交流！ 2011-06-15 19:56:28
引用回帖:
Originally posted by NaOHWTO at 2011-06-14 20:24:21:
重力加速度只是一个方向的线加速度，例如Y向的，所以UDF中用源项加入到Y的动量方程中去，不过，ADJUST宏只是为了调整你给定激励的时间，例如，我的第一个算例，2s后刹车停住了，就没有加速度了，所以这个时候就不 ...

大侠你好！
油箱液面晃动的问题我也做过，不过就是练练手学习学习而已，所以做得很粗糙，远远没有你做出来的这么漂亮！
不知道你是否方便把你做的cas分享下，我不要UDF，只想看看模型的各种设置深入学习下，非常感谢。
如果不方便分享的话，我想请教下，你做这个cas的网格数大约是多少？壁面是怎么设置的？表面张力系数大约是多少？还有我看到你做出来的液面有那种像是波干涉出来的斑纹，很真实，这个是有什么特别的处理方法么？
赞一下(2人) 回复此楼
8楼2011-06-14 20:47:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NaOHWTO

金虫 (小有名气)
把你邮箱给我
回复此楼
9楼2011-06-14 21:55:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lixiaoli_03

新虫 (初入文坛)
★
zcjwyl(金币+1):谢谢参与
引用回帖:
Originally posted by 47425804 at 2011-06-14 20:47:38:
大侠你好！
油箱液面晃动的问题我也做过，不过就是练练手学习学习而已，所以做得很粗糙，远远没有你做出来的这么漂亮！
不知道你是否方便把你做的cas分享下，我不要UDF，只想看看模型的各种设置深入学习下， ...

shixuefeng0808@126.com
赞一下(1人) 回复此楼
11楼2011-07-12 17:26:29
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

模拟船舶晃动

新虫 (初入文坛)
★
zcjwyl(金币+1):谢谢参与
各位前辈好，我也在做晃动模拟，模拟船舶上的油柜晃动。现在只知道个大方向可能要用到6DOF和VOF，恳请指点一二。谢谢
赞一下 回复此楼
12楼2012-02-11 10:45:46
 已阅   回复此楼   关注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的回帖

dubaikobe

新虫 (小有名气)
★
zcjwyl(金币+1): 谢谢参与
不明觉厉，lz能不能留个qq或者群 啥的方便交流啊
赞一下 回复此楼
14楼2013-10-29 10:11:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (初入文坛)
★
zcjwyl(金币+1): 谢谢参与
本帖仅楼主可见
15楼2014-11-20 12:56:01
 已阅   申请仿真EPI   回复此楼   编辑   查看我的主页

bluecore007

新虫 (初入文坛)
★
zcjwyl(金币+1): 谢谢参与
正在做油箱晃动，感谢分享
回复此楼
16楼2015-07-30 10:56:34
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Wwyf

禁虫
★
zcjwyl(金币+1): 谢谢参与
本帖内容被屏蔽
17楼2023-03-12 15:48:46
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Wwyf

禁虫
★
zcjwyl(金币+1): 谢谢参与
本帖内容被屏蔽
18楼2023-03-12 20:25:45
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
huagai10楼
2011-07-01 16:05   回复  
zcjwyl(金币+1):谢谢参与
臭水沟(金币-1): 禁止灌水~~ 2011-07-01 19:59:50
相关版块跳转 我要订阅楼主 zcjwyl 的主题更新
181/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭