Udf for heat source
A porous medium is included in a domain. I want to define a heat source to the porous medium which will be a function of the temperature difference between the porous temperature and the domain temperature. Thanks for any help, Thomas

Re: Udf for heat source
It's possible. The DEFINE_SOURCE macro is available for your problem.
Keil 
Re: Udf for heat source
If you're using the standard porous media model in Fluent 5 or 6 you won't be able to model interphase heat transfer or have a model that assumes nonequilibrium between the two phases. Fluent assumes, a single energy equation for solid and fluid as described in the manual.
To model a porous media without thermal equilibirum you need to use udfs/udss... I have presented some code for this in 2000 and 2001. Its also in the udf repository....search the archives for the posts/code. Greg 
Re: Udf for heat source
Interesting,
where to find this udf repository? thanks in advance, Laika, still orbiting 
Re: Udf for heat source
Repository at
ftp://ftp.cfdonline.com/fluent and my code at ftp://ftp.cfdonline.com/fluent/udf_scalar_dos.c Greg 
Dear Greg Perkin,
This is to let you know that I am working with heat transfer in porous media and have found your replies on this site very helpful. I wanted to download "udf_scalar_dos.c" but unable to download. I am using my user name and password of this site to download but it showing unable to login. Please let me know how to download it, or if possible,please email the file to me. Thanking you Sajid 
(define_source) heat generation changs with location
hello everyone!
I write a udf for heat generation with different location(xyz);),but only a part of fluid has heat generation:confused:. project:a solid (0.08X0.08X0.1 to xyz) with 2 cylinder(r=0.01),cylinder is fluid with heat generation.:) the code:#include "udf.h" DEFINE_SOURCE(heat_flux, cell, t, dS, eqn) { real a, b, c; real x[2] = { 0.02, 0.02 }, y[2] = { 0,0 }; //x、y source location real p[8] = { 0.025, 500000000, 0.025, 500000000, 0.075, 500000000, 0.075, 500000000 };//0.025、0.075 is Z location，source value follows. real pos[ND_ND]; real source=0; int i, j = 0; C_CENTROID(pos, cell, t); a = pos[0]; //abc acquire xyz b = pos[1]; c = pos[2]; if ((c >= p[0]  0.025) && (c <= p[0] + 0.025)) //judge Z location { i=0; while(i<2){ if ((a <= x[i] + 0.01 && a >= x[i]  0.01) && (b <= y[i] + 0.01 && b >= y[i]  0.01)) //judge xy，if meeting the range,give source right value {source = p[0 + i+1]; dS[eqn] = 0; i=2;} i=i+1; } } else if ((c >= p[4]  0.025) && (c <= p[4] + 0.025)) { i=0; while(i<2){ if ((a <= x[i] + 0.01 && a >= x[i]  0.01) && (b <= y[i] + 0.01 && b >= y[i]  0.01)) {source = p[4+2*i+1]; dS[eqn] = 0; i=2;} i=i+1; } } else { source=0; dS[eqn]=0; } return source; } the result shows the right cylinder has heat generation:) ,but the other one has zero.:( thanks for your kindness.contact me:956400837@qq.com 
(define_source) heat generation changs with location
hello everyone!
I write a udf for heat generation with different location(xyz);),but only a part of fluid has heat generation:confused:. project:a solid (0.08X0.08X0.1 to xyz) with 2 cylinder(r=0.01),cylinder is fluid with heat generation.:) the code:#include "udf.h" DEFINE_SOURCE(heat_flux, cell, t, dS, eqn) { real a, b, c; real x[2] = { 0.02, 0.02 }, y[2] = { 0,0 }; //x、y source location real p[8] = { 0.025, 500000000, 0.025, 500000000, 0.075, 500000000, 0.075, 500000000 };//0.025、0.075 is Z location，source value follows. real pos[ND_ND]; real source=0; int i, j = 0; C_CENTROID(pos, cell, t); a = pos[0]; //abc acquire xyz b = pos[1]; c = pos[2]; if ((c >= p[0]  0.025) && (c <= p[0] + 0.025)) //judge Z location { i=0; while(i<2){ if ((a <= x[i] + 0.01 && a >= x[i]  0.01) && (b <= y[i] + 0.01 && b >= y[i]  0.01)) //judge xy，if meeting the range,give source right value {source = p[0 + i+1]; dS[eqn] = 0; i=2;} i=i+1; } } else if ((c >= p[4]  0.025) && (c <= p[4] + 0.025)) { i=0; while(i<2){ if ((a <= x[i] + 0.01 && a >= x[i]  0.01) && (b <= y[i] + 0.01 && b >= y[i]  0.01)) {source = p[4+2*i+1]; dS[eqn] = 0; i=2;} i=i+1; } } else { source=0; dS[eqn]=0; } return source; } the result shows the right cylinder has heat generation:) ,but the other one has zero.:( thanks for your kindness. 
All times are GMT 4. The time now is 20:41. 