[amirhabibdoost]

Dear friends
I face with a problem when I want to reform my UDF. Indeed , I want to add a function ( sin(x)) to my udf , but I don’t know what is necessary corrections.
Please let me any suggestion in order to correct my udf.
Tanks in advance.
Amir

#include "udf.h"

real velocity (real t)
{
real a, b, value;

a = 10.0;
b = 0.34;
value = 0.0;

if ((t >= 0) || (t < 0.2))
{
value = sin (a*t);( here , I want to add sin(x) , but I don’t know what corrections is need.)
}
else
{
value = b;
}

return value;
}

DEFINE_PROFILE(unsteady_velocity, thread, position)
{
face_t f;
real t;

t = CURRENT_TIME;

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

[amirhabibdoost]

Originally Posted by amirhabibdoost

Dear friends
I face with a problem when I want to reform my UDF. Indeed , I want to add a function ( sin(x)) to my udf , but I don’t know what is necessary corrections.
Please let me any suggestion in order to correct my udf.
Tanks in advance.
Amir

#include "udf.h"

real velocity (real t)
{
real a, b, value;

a = 10.0;
b = 0.34;
value = 0.0;

if ((t >= 0) || (t < 0.2))
{
value = sin (a*t);( here , I want to add sin(x) , but I don’t know what corrections is need.)
}
else
{
value = b;
}

return value;
}

DEFINE_PROFILE(unsteady_velocity, thread, position)
{
face_t f;
real t;

t = CURRENT_TIME;

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

urgent help , please

[eng_s_sadeghi]

Hi Dear Mr. Habibdoust
I know you very well
you should use this:

DEFINE_PROFILE(unsteady_velocity, t, i)
{
real pos[ND_ND], x;
face_t f;
begin_f_loop(f, t)
{
F_CENTROID(pos,f,t);
x=pos[0];
F_PROFILE(f, t, i) = sin(x);
}
end_f_loop(f, thread)
}

[amirhabibdoost]

Dear Saeed,
I know you very well too.
Tank for your attention and reply.
you mean that if i use your suggested code, I will have velocity inlet as bellow function:
sin(x).sin(t)
Tanks in advance
Amir