[huela]

Hello,
I have a problem with my udf. I want if the Pressure falls under 0 that the value of Pressure would be set 0.

#include <udf.h>
DEFINE_ON_DEMAND(hallo)
{
Domain *d;

real patm = 101325.0;
real pmin = 0.0;
real pressure = 0;
Thread *t;
cell_t c;
d = Get_Domain(1);

thread_loop_c(t,d)
{
begin_c_loop(c,t)
{

pressure = C_P(c,t);

if (pressure < 0.)
pressure = 0.0;
}
end_c_loop(c,t)

begin_c_loop(c,t)
{
pressure = C_P(c,t);
C_UDMI(c,t,0) = pressure;
}
end_c_loop(c,t)
}

}


Thank you for help!!

[coglione]

DEFINE_PROFILE() is the way to go to fix a flow variable to a certain value for any cell in your grid.
Have a look at the UDF-Userguid section 2.3.15, there is also an example that might be helpful.

BTW, what is the reason for your intention to prevent negative pressure values? If it seems to you that this is unphysical, keep in mind that in Fluent pressures are always defined as gauge pressure to an operating pressure defined by you and static pressure may drop significantly for high speed flows and may become negative which does not necessarily mean that your solution is wrong.

cheers

[chaozhong.qin]

Hi,

I found a small mistake in your UDF, as follows:

begin_c_loop(c,t)
{

pressure = C_P(c,t);

if (pressure < 0.)
pressure = 0.0;
(note that: C_P(c,t) = pressure;
return the value to C_P(c,t)
)

}
end_c_loop(c,t)

[huela]

Thank you for answer!

But what I have to change in my UDF?

Thank you!

[sozen]

Originally Posted by huela

Thank you for answer!

But what I have to change in my UDF?

Thank you!

Could you solve the problem?