[burak_ekmekci]

Hi guys ;

- How can ı create this udf (below)? and is it possible to create a udf for discharge, not velocity ?

between time =0 second, t=360 discharge=14
between t=360, t= 690 d=0.14444(x-360)+14
between t=690, t=1020 d=61.67-0.14444(x-690)
between t=1020, t=1320 d=14


thanks ...

[burak_ekmekci]

is there any comment or suggestion ?

I have created this codes that is below. it is working but there is a problem with increasing velocity. when time gets 360 s, velocity does not increase ? what might be problem ?

please help.

thanks.

#include "udf.h"



real Dt1=360;
real Dt2=330;
real Dt3=330;
real Dt4=360;


real Tevapfz1=0.182;
real Tevapfz2=0.182;
real Tevapfz3=0.4;
real Tevapfz4=0.182;
real Tevapfz5=0.182;


real A=360;
real B=0.182;

real c_time=0.0;
real aralik=0.0;


{
c_time=CURRENT_TIME;

if(c_time<Dt1)
{
A=0;
B=Tevapfz1;
aralik=1.0;
}

if((c_time<(Dt1+Dt2))&&(c_time>=Dt1))
{
A=(Tevapfz2-Tevapfz3)/(-Dt2);
B=(Tevapfz3*Dt1-Tevapfz2*(Dt1+Dt2))/(-Dt2);
aralik=2.0;
}

if((c_time<(Dt1+Dt2+Dt3))&&(c_time>=(Dt1+Dt2)))
{
A=(Tevapfz3-Tevapfz4)/(-Dt3);
B=(Tevapfz4*(Dt1+Dt2)-Tevapfz3*(Dt1+Dt2+Dt3))/(-Dt3);
aralik=3.0;
}

if((c_time<(Dt1+Dt2+Dt3+Dt4))&&(c_time>(Dt1+Dt2+Dt 3)))
{
A=0;
B=(Tevapfz5*(Dt1+Dt2+Dt3)-Tevapfz4*(Dt1+Dt2+Dt3+Dt4))/(-Dt4);
aralik=4.0;
}



printf("VELOCITY=%f\n\n", (A*CURRENT_TIME+B));
printf("ARALIK=%f\n\n", aralik);

}


DEFINE_PROFILE(evapfz_temp, thread, index)
{
real x[ND_ND];
real y;
face_t f;

begin_f_loop(f, thread)
{
F_CENTROID(x,f,thread);
F_PROFILE(f, thread, index) = A*c_time+B;
}
end_f_loop(f, thread)

}