| 查看: 3444 | 回复: 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) } |
回复此楼» 猜你喜欢
288资源与环境专硕求调剂，不限专业，有学上就行
已经有6人回复
-
生物学学硕，一志愿湖南大学，初试成绩338
已经有8人回复
-
266分，求材料相关专业调剂
已经有6人回复
-
304材料求调剂
已经有5人回复
-
359求调剂
已经有5人回复
-
085600 286分 材料求调剂
已经有12人回复
-
各位老师好，我的一志愿为北京科技大学085601材料专硕
已经有10人回复
-
一志愿：西北大学，英一数一408-284分求调剂
已经有5人回复
-
291求调剂
已经有10人回复
-
一志愿郑州大学，080500学硕，总分317分求调剂
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
ykxyer
铁杆木虫 (文坛精英)
- 应助: 2 (幼儿园)
- 金币: 9436.9
- 散金: 8350
- 红花: 71
- 沙发: 188
- 帖子: 17480
- 在线: 672.6小时
- 虫号: 1673326
- 注册: 2012-03-07
- 性别: GG
- 专业: 金属功能材料
3楼2017-07-04 10:16:11
2楼2017-07-04 10:06:22
4楼2017-07-04 10:37:44
ykxyer
铁杆木虫 (文坛精英)
- 应助: 2 (幼儿园)
- 金币: 9436.9
- 散金: 8350
- 红花: 71
- 沙发: 188
- 帖子: 17480
- 在线: 672.6小时
- 虫号: 1673326
- 注册: 2012-03-07
- 性别: GG
- 专业: 金属功能材料
|
我现在把程序改成下面的了，但是还是有报错：Error: Update_Time_Level: invalid data，百度了一下，也没找到类似错误的是啥原因？ #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) { 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"; heat_flux=C_all*Density_all*(t_Al+t_CRFP)*(C_T_M1(c,t) - C_T_M2(c,t))/current_timestep; F_PROFILE(f,thread,position)=heat_flux; } end_f_loop(f,thread) } |
5楼2017-07-04 10:42:39
