[Nicolastheterminator]

Hi everyone,

I have a problem with my porous flow UDFs. My geometry is a simple 2D channel with a porous region blocking flow. I specify a pressure inlet and outlet.

When solving my problem with the porous flow solver in Fluent, I get a nice reasonable flow. But when I do it with a UDF defining porous flow with a source UDF, I get the same result except some very ugly redlines going into and out of the porous region that are entirely unexpected (see image). My UDFs should be fine as they are simple Darcian flow so I wonder what is wrong?

I suspect it is the mesh, but I would really appreciate advice on what to do.

Thanks,
Nick

Here is my Porous flow UDF

#include "udf.h"
#define urf 0.1 // under-relaxation factor
real s1 = 0. , s2 = 0.;

DEFINE_SOURCE(Darcian_x,c,t,dS,eqn)
{
//k is the permeability, defined in m^2

real source, k, mu;

mu = C_MU_L(c,t); k = 1.4E-11;

source = - C_U(c,t)*mu/k;

dS[eqn] = - mu/k;

s1 = s1*(1-urf) + source*urf;

return s1;

}

DEFINE_SOURCE(Darcian_y,c,t,dS,eqn)
{
//k is the permeability, defined in m^2

real source, k, mu;
mu = C_MU_L(c,t); k = 1E-20;

source = - C_V(c,t)*mu/k;

dS[eqn] = - mu/k;

s2 = s2*(1-urf) + source*urf;

return s2;

}