| 查看: 3260 | 回复: 9 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
ykxyer铁杆木虫 (文坛精英)
|
[求助]
fluent UDF如何实现瞬态，壁面heat flux为 上一时刻温度与此时温度差的函数的？
|
|
|
壁面热流边界满足公式heat flux=C_all*Density_all*t_all*(T1-T0))/delta_t；其中 C_all为总的热容； Density_all为总的密度； t_all为总的厚度； T1为前一时刻的壁面温度； T0为此时的壁面温度； delta_t为两温度时间的差值； 壁面边界处，材料会随着温度的升高而吸收一部分热量，而这部分传热以heat flux形式表示，fluent中heat flux的单位为W/m2， 因此就简化了公式Q=C*m*(T1-T2)，m=density*Volume近似等于：密度*厚度*面积，这部分热量不是随着时间变化的，而是随着温度变化的， fluent中的单位为W/m2，在计算热量 时会乘以时间，因此在UDF中除了时间。 /////////////////////////////////////////////////////////////////////////////// 下面是我自己写得UDF,但是初始化的时候有如下提示的错误： Error: FLUENT received fatal signal (ACCESS_VIOLATION) 1. Note exact events leading to error. 2. Save case/data under new name. 3. Exit program and restart to continue. 4. Report error to your distributor. Error Object: #f /////////////////////////////////////////////////////////////////////////////// #include "udf.h" /*5-the wall heat flux*/ /*-----------------------------------------------------------------------*/ DEFINE_PROFILE(wall_heat_flux, thread, position) { /*--------------壁面材料储热部分--------------------*/ real Density_Al=2700; /*铝的密度*/ real Density_CRFP=1570; /*CRFP的密度*/ real C_Al=902; /*铝的热容*/ real C_CRFP=840; /*CRFP的热容*/ real t_Al=0.0048; /*铝的厚度*/ real t_CRFP=0.0025; /*CRFP的厚度*/ real w1,w2; /*w1为铝的质量分数，w2为CRFP的质量分数*/ real Density_all,C_all; /*壁面材料总的密度，热容*/ real heat_flux; /*壁面材料储热部分近似热流量*/ real current_time step; /*时间步长*/ real t; cell_t c; Thread *t; face_t f; begin_f_loop(f,thread) { t_wall=F_T(f,thread); heat_flux_1=(1/(1/h_hydrogen+thickness_all/k_wall+1/h_air))*(t_air-t_wall); /*对流交换单位面积氢气热流密度*/ w1=Density_Al*t_Al/(Density_Al*t_Al+Density_CRFP*t_CRFP); w2=Density_CRFP*t_CRFP/(Density_Al*t_Al+Density_CRFP*t_CRFP); Density_all=1/(w1/Density_Al+w2/Density_CRFP); C_all=C_Al*w1+C_CRFP*w2; current_timestep = RP_Get_Real("physical-time-step"  ;C_UDMI(c,t,0) = C_T_M1(c,t); /*获得前一时刻温度值*/ C_UDMI(c,t,1)= C_T_M2(c,t);/*获得前前时刻温度值*/ heat_flux=C_all*Density_all*(t_Al+t_CRFP)*(C_UDMI(c,t,0) - C_UDMI(c,t,1))/current_timestep; F_PROFILE(f,thread,position)=heat_flux; } end_f_loop(f,thread) } |
回复此楼» 猜你喜欢
博士读完未来一定会好吗
已经有24人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
ykxyer
铁杆木虫 (文坛精英)
- 应助: 2 (幼儿园)
- 金币: 9436.9
- 散金: 8350
- 红花: 71
- 沙发: 188
- 帖子: 17480
- 在线: 672.6小时
- 虫号: 1673326
- 注册: 2012-03-07
- 性别: GG
- 专业: 金属功能材料
9楼2017-07-04 16:54:57
2楼2017-07-04 10:06:22
ykxyer
铁杆木虫 (文坛精英)
- 应助: 2 (幼儿园)
- 金币: 9436.9
- 散金: 8350
- 红花: 71
- 沙发: 188
- 帖子: 17480
- 在线: 672.6小时
- 虫号: 1673326
- 注册: 2012-03-07
- 性别: GG
- 专业: 金属功能材料
3楼2017-07-04 10:16:11
4楼2017-07-04 10:37:44