[sdsuronnie]

#include "udf.h"

DEFINE_PROFILE(velocity_magnitude, t, i)
{
real velocity;
real the_current_time;
face_t f;

the_current_time = CURRENT_TIME;



if ((the_current_time>=0.0) && (the_current_time<0.25))
{
F_PROFILE(f, thread, position) =-39.04*t*t+9.76*t;
}
if ((the_current_time>=0.25) && (the_current_time<0.3))
{
F_PROFILE(f, thread, position) =89.6*t*t-49.28*t+6.72;
}
if ((the_current_time>=0.3) && (the_current_time<=1))
{
F_PROFILE(f, thread, position) =0;
}

begin_f_loop(f,t)
{
F_PROFILE(f,t,i) = velocity;
}
end_f_loop(f, t)
}

Tha above is the code , however when i try to interpret it, it has err
hope experts here can help me


best regard
Ronnie
sdsuronnie@gmail.com

thanks

[sdsuronnie]

it is an inlet velocity boundary condition

[mprinkey]

DEFINE_PROFILE is called for an entire face thread...the "t" in the DEFINE_PROFILE macro at the top of the UDF. You *must* loop over all of the faces on the thread and set F_PROFILE for each face. That is happening in what appears to be vestigial code at the bottom there, but not in the conditional test.

Inside of your if statements checking the_current_time, you need to have those loops. And your are using the thread pointer t in your piecewise definition of the velocity, which makes no possible sense. I assume that should be the_current_time.

...
if ((the_current_time>=0.0) && (the_current_time<0.25))
{
 begin_f_loop(f,t)
 {
  F_PROFILE(f, thread, position) =-39.04*the_current_time*the_current_time+9.76*the_current_time;
 }
end_f_loop(f, t)
}
...

[sdsuronnie]

Thank you for your comments, really help me a lot

best regard

Ronnie