| 查看: 5420 | 回复: 2 | |||
[交流]
Fluent udf 设置管道入口速度等于出口速度 已有2人参与
|
|
做的一个微管道气-液两相流的模拟，在管道中截取一段流动单元，流动单元的中间有一个气泡。为了维持这个气泡一直在流动单元内，需要控制入口流速，所以打算使用udf使入口速度始终等于出口速度，再加上网格区域的一个平移，应该就能实现这一过程了。写完udf可以完成interpret，但是开始计算时会出现“Error: received a fatal signal (Segmentation fault).”不是很明白为什么 用periodic好像也可以实现这个过程，但只能通过pressure gradient来控制流动状态而不能用mass flow rate。虽然说也能凑合用就是了。。。 文件什么的都传附件了。。跪！求！指！教！ udf如下： #include "udf.h" #define U_MAX 0.5 real store[50]={0.1,0.2,0.3,0.4,0.5}; /* 用数组store保存速度值，考虑过用UDM什么的但是还不清楚应该怎么设置； 随便取了几个初始值，在fluent里面只使用PROFILE而不使用INIT和ADJUST时，可以正常计算，并且能通过fluent里面的plot观察到与初始值对应的速度分布，但是一旦使用INIT就会立刻提示错误； 同时使用PROFILE和ADJUST能够进行一段时间的计算，随后才提示错误，而且此时plot出口速度仍为初始值对应的分布 */ DEFINE_INIT(fully_developed,d) /* 初始化管内流动为充分发展 */ { cell_t c; real x[ND_ND]; real y=0.; Thread *t=Lookup_Thread(d,9); begin_c_loop(c,t) { C_CENTROID(x,c,t); y=x[1]; C_U(c,t)=U_MAX-U_MAX*y*y/(0.000055*0.000055); /* 矩形管道，宽度为0.111mm，因为取对称了所以宽度0.0555mm，速度呈抛物线形分布 */ } end_c_loop(c,t) } DEFINE_ADJUST(outlet_to_inlet,d) /* 将出口thread上的速度赋予数组 */ { int j=0; real u=0.; face_t f; Thread *thread_out=Lookup_Thread(d,12); begin_f_loop(f,thread_out) { u=F_U(f,thread_out); store[j]=u; j++; } end_f_loop(f,thread_out) } DEFINE_PROFILE(udf_inlet,thread,index) /* 将数组中存储的速度值赋予入口 */ { int j=0; real u=0.; real x[ND_ND]; face_t f; begin_f_loop(f,thread) { F_CENTROID(x,f,thread); u=store[j]; F_PROFILE(f,thread,index)=u; j++; } end_f_loop(f,thread); }  unit.jpg |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : unit.zip
2017-05-08 16:37:59, 1.02 M
» 猜你喜欢
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有80人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有3人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
malasong9471
木虫 (正式写手)
- 应助: 58 (初中生)
- 金币: 7278
- 红花: 29
- 帖子: 639
- 在线: 117.1小时
- 虫号: 4100890
- 注册: 2015-09-25
- 性别: GG
- 专业: 工程热力学
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
仅供参考，http://muchong.com/bbs/viewthread.php?tid=10667013&_from=share&u=4100890&h=4285a81d 发自小木虫Android客户端 |
2楼2017-05-17 11:09:18
18093194258
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 50.3
- 帖子: 25
- 在线: 8小时
- 虫号: 4977017
- 注册: 2016-09-03
- 性别: GG
- 专业: 工程热物理相关交叉领域
3楼2019-01-21 20:46:48
10