|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Franz Roman
Join Date: Jun 2009
Posts: 28
Rep Power: 17 ![]() |
Hello,
I have a problem with the non-eq. porous medium. To calculate the heat transfer coefficient between the solid material in the porous zone and the air I want to use a user defined profile. I want to access the air properties like density, specific heat, thermal conductivity and viscosity, and also the velocity magnitude in each element of the porous zone, from within the UDF. I included some UDMs to store the values of these variables and check that they are correct, and I saw that the density and the viscosity of the fluid agree with the set up in the Materials Panel, but for specific heat and thermal conductivity the UDMs have values of almost zero (i.e. -1.42e37), although in the Materials Panel they have correct constant values. Also, the velocity is displayed as uniform at 1.0 m/s in a plane through the porous media, although the contours of velocity magnitude show a variable velocity, which had been calculated previously. The UDF is below, it seems ok to me. Or is accessing some variables not allowed from withing a DEFINE_PROFILE? Another issue is that, when Fluent created the solid porous zone, the two open faces of the porous zone (though which air could go in and out) were converted to walls, and it is not possible to convert them to interior. I would appreciate some assistance. #include "udf.h" /* Heat and moisture transfer in tomato cooling */ /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ #define dp 0.07 #define porosity 0.44 DEFINE_PROFILE(heat_transf_coeff,t,i) { real rho_air, mu, k_air, cp_air, va, Re, Pr, Nu, ht; cell_t c; begin_c_loop(c,t) { rho_air = C_R(c,t); mu = C_MU_L(c,t); k_air = C_K_L(c,t); cp_air = C_CP(c,t); va = pow(pow(C_U(c,t),2) + pow(C_V(c,t),2) + pow(C_W(c,t),2),(1/2)); /* velocity magnitude at cell */ Re = (va/porosity)*dp/(mu/rho_air); Pr = mu*cp_air/k_air; Nu = 2.0 + 0.552*pow(Re, 0.53)*pow(Pr, 0.33); ht = Nu*k_air/dp; C_UDMI(c,t,0) = rho_air; C_UDMI(c,t,1) = mu; C_UDMI(c,t,2) = k_air; C_UDMI(c,t,3) = cp_air; C_UDMI(c,t,4) = va; F_PROFILE(c,t,i) = ht; } end_c_loop(c,t) } |
|
![]() |
![]() |
![]() |
Tags |
non-equilibrium thermal, porous medium, udf |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Non-equilibrium thermal model for porous medium | mmunige | FLUENT | 1 | July 14, 2016 06:03 |
Use of non-thermal equilibrium in porous media model | highhopes | FLUENT | 0 | November 9, 2014 20:40 |
local thermal non equilibrium model in porous media | adi.ptb | Fluent UDF and Scheme Programming | 0 | November 1, 2014 15:00 |
Water subcooled boiling | Attesz | CFX | 7 | January 5, 2013 03:32 |
porous filtration membrane area in model | catherina | CFX | 3 | October 27, 2011 06:05 |