What is the udf for permeability or viscous resistance in Fluent 6.1?

#include "udf.h" #include "sg.h"

#define R1soil 8e13 /* Resistances */ #define R2soil 2e13

/* Define Resistences */

DEFINE_PROFILE(one_res, thread, np) { cell_t c; real x[ND_ND]; real y;

begin_c_loop (c,thread)

{ C_CENTROID(x,c,thread);

y =x[1];

if ((y<h1 && y>=h2) || (y<h3 && y>=h4) || (y<h5 && y>=h6 ) )

F_PROFILE(c,thread,np) = R1soil; else if ((y<h2 && y>=h3) || (y<h4 && y>=h5) || (y<h6))

F_PROFILE(c,thread,np) = R1graval;

} end_c_loop (c,thread) }

What is the h's (h2,h3,h4,h5,h6) stand for?

Thank you

THis is UDF for a Packed column.. h1..h2.are column heights..i.e. distance along height...the udf specifies resistances for two layers of packings..

i had used it successfully

hope this helps.. by -Umesh

is there some syntax error in the last 2-3 lines of the code given by Umesh?

hi the code is correct... but some symbols were not displyed... just refer to this one

--The Code-----------

#include "udf.h"

#include "sg.h"

#define R1soil 8e13 /* Resistances */

#define R2soil 2e13

#define h1 1.8

#define h2 0.3

/* Define Resistences */

DEFINE_PROFILE(one_res, thread, np) { cell_t c;

real x[ND_ND];

real y;

begin_c_loop (c,thread)

{ C_CENTROID(x,c,thread);

y =x[1];

if (y less than h1 and and y greater than= h2)

F_PROFILE(c,thread,np) = R1soil;

else if (y less than h2)

F_PROFILE(c,thread,np) = R1graval;

} end_c_loop (c,thread) }

[reynolds039]

hello
what is the viscous resistance in multiphase flow throw a porous medium for each phase
(oil-water) , give the udf for each of them please