CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > ANSYS > FLUENT

Non-eq. porous model

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   August 10, 2020, 13:41
Default Non-eq. porous model
New Member
Franz Roman
Join Date: Jun 2009
Posts: 28
Rep Power: 15
franzdrs is on a distinguished road
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
real rho_air, mu, k_air, cp_air, va, Re, Pr, Nu, ht;
cell_t c;

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;
franzdrs is offline   Reply With Quote


non-equilibrium thermal, porous medium, udf

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Non-equilibrium thermal model for porous medium mmunige FLUENT 1 July 14, 2016 07:03
Use of non-thermal equilibrium in porous media model highhopes FLUENT 0 November 9, 2014 21:40
local thermal non equilibrium model in porous media adi.ptb Fluent UDF and Scheme Programming 0 November 1, 2014 16:00
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32
porous filtration membrane area in model catherina CFX 3 October 27, 2011 07:05

All times are GMT -4. The time now is 02:14.