| 查看: 1725 | 回复: 5 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
模拟同轴套管式换热器，编UDF将一个出口边界上的温度作为另一个进口边界条件的温度值 已有2人参与
|
||||
|
大家好，我现在正在利用Fluent模拟地下同轴套管式换热器的各参数对其性能的影响，类似于模拟地源热泵的U型管换热器，但我的钻孔深度在2000米左右，远超地源热泵技术的钻孔深度。 我的模型如下图所示，我取模型的1/4进行三维研究，环腔和内管流体为水，两者之间是内管璧，环腔外为钻孔区域，钻孔外为岩石区域，共五个体。现在模拟同轴套管式换热器的换热量，流体水采用外进内出，且是稳态模拟。 我现在遇到了一些问题，我把套管式换热器的环腔进口设置为速度进口1（inlet1）、环腔出口设置为压力出口1（outlet1）、内管下面设置为速度进口2（inlet2）、内管上面设置为压力出口2（outlet2），我所关注及要求的就是内管出口即压力出口2（outlet2）的平均温度。 问题是环腔出口即压力出口1（outlet1）的平均温度，必须与内管下面即速度进口2（inlet2）的平均温度保持一致，即必须把环腔出口即压力出口1（outlet1）的边界条件设置为内管下面即速度进口2（inlet2）的边界条件，也就是两者温度保持一致，请问怎么用UDF编写来实现这个过程？最好有详细编程，我是Fluent新手，再此谢谢大家了！  图片1.png |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有41人回复
-
filecode
已经有9人回复
-
请教兼职经验
已经有3人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有5人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有13人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有11人回复
-
读博的好处
已经有4人回复
-
filecode
已经有17人回复
-
FileCode能看出啥？
已经有32人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有15人回复
dai26520
专家顾问 (著名写手)
-
专家经验: +219 - 仿真EPI: 1
- 应助: 224 (大学生)
- 金币: 4097
- 散金: 4143
- 红花: 107
- 帖子: 1850
- 在线: 529.6小时
- 虫号: 1619980
- 注册: 2012-02-16
- 专业: 高分子合成化学
- 管辖: 仿真模拟
3楼2017-02-04 18:17:23
2楼2017-01-30 22:11:32
xingfuww
专家顾问 (正式写手)
-
专家经验: +78 - 仿真EPI: 2
- 应助: 46 (小学生)
- 金币: 6637.4
- 散金: 52
- 红花: 29
- 帖子: 542
- 在线: 284小时
- 虫号: 1776059
- 注册: 2012-04-24
- 性别: GG
- 专业: 化学工程及工业化学
- 管辖: 仿真模拟
4楼2017-02-05 10:48:39
xingfuww
专家顾问 (正式写手)
-
专家经验: +78 - 仿真EPI: 2
- 应助: 46 (小学生)
- 金币: 6637.4
- 散金: 52
- 红花: 29
- 帖子: 542
- 在线: 284小时
- 虫号: 1776059
- 注册: 2012-04-24
- 性别: GG
- 专业: 化学工程及工业化学
- 管辖: 仿真模拟
【答案】应助回帖
|
我这个是以前帮别人写的一个 可以同步的速度入口 你可以参考一下。 一侧的出口 变为另外一侧入口。 # include "udf.h" #include <stdio.h> #include <stdlib.h> face_t f; Domain *domain; Thread *t_out; DEFINE_EXECUTE_AT_END(flux_end) { real outlet_flux; real flow_time = RP_Get_Real("flow-time"  ;FILE *f_frc; real flow=0.0; domain = Get_Domain(1); t_out=Lookup_Thread(domain,2); begin_f_loop(f,t_out) { flow+=F_FLUX(f,t_out); } end_f_loop(f,t_out) outlet_flux=flow; /*printf("%f,%f\n",flow_time, outlet_flux); */ f_frc=fopen("face_interior1_flux.txt","w" ;fprintf(f_frc,"%f %f\n",flow_time, outlet_flux); fclose(f_frc); } DEFINE_PROFILE(symmetry_outlet_flux,t,i) { real flow_time,outlet_flux; FILE *f_frc; f_frc=fopen("face_interior1_flux.txt","r" ;fscanf(f_frc,"%lf %lf",&flow_time,&outlet_flux);/*flow_time也读入了但是这里没用到*/ fclose(f_frc); begin_f_loop(f,t) { F_PROFILE(f,t,i) = -outlet_flux; } end_f_loop(f,t) } DEFINE_PROFILE(aa,t,i) { face_t f; real time; int a; time=CURRENT_TIME; begin_f_loop(f,t) { a=(int)(time/10); if(time>=(10*a+0)&&time<=(10*a+5)) F_PROFILE(f,t,i) =-(time-10*a)*(time-10*a) ;/*周期10秒*/ else F_PROFILE(f,t,i) =0 ; } end_f_loop(f,t) } |
5楼2017-02-05 11:02:40
