[求助] 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) }

回复此楼

» 猜你喜欢

• 085602化工求调剂（331分） 已经有8人回复
• 307求调剂 已经有18人回复
• 085701环境工程求调剂 已经有12人回复
• 288资源与环境专硕求调剂，不限专业，有学上就行 已经有11人回复
• 085600材料与化工调剂 已经有9人回复
• 一志愿华东师范大学有机化学专业，初试351分，复试被刷求调剂! 已经有8人回复
• 266分，求材料相关专业调剂 已经有7人回复
• 材料求调剂 一志愿哈工大总分298分，前三科223分 已经有7人回复
• 303求调剂 已经有7人回复
• 083000学硕274求调剂 已经有9人回复

» 本主题相关商家推荐: (我也要在这里推广)