| 查看: 3707 | 回复: 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) } |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
国自然面上复盘~欢迎讨论
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
基金不中，共勉
已经有5人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有3人回复
-
我不理解！
已经有23人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有10人回复
-
系统查不到
已经有10人回复
-
怎么查啊
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
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
ykxyer6楼2017-07-04 11:20:59
ykxyer
铁杆木虫 (文坛精英)
- 应助: 2 (幼儿园)
- 金币: 9436.9
- 散金: 8350
- 红花: 71
- 沙发: 188
- 帖子: 17480
- 在线: 672.6小时
- 虫号: 1673326
- 注册: 2012-03-07
- 性别: GG
- 专业: 金属功能材料
7楼2017-07-04 13:32:43
ykxyer
铁杆木虫 (文坛精英)
- 应助: 2 (幼儿园)
- 金币: 9436.9
- 散金: 8350
- 红花: 71
- 沙发: 188
- 帖子: 17480
- 在线: 672.6小时
- 虫号: 1673326
- 注册: 2012-03-07
- 性别: GG
- 专业: 金属功能材料
8楼2017-07-04 13:34:16
ykxyer
铁杆木虫 (文坛精英)
- 应助: 2 (幼儿园)
- 金币: 9436.9
- 散金: 8350
- 红花: 71
- 沙发: 188
- 帖子: 17480
- 在线: 672.6小时
- 虫号: 1673326
- 注册: 2012-03-07
- 性别: GG
- 专业: 金属功能材料
9楼2017-07-04 16:54:57
10楼2019-04-17 13:07:34
投票:
