[rakadit]

Dear all!

Please help me to write UDF for wall shear stress equation: tau=constant*viscosity*velocitygradient.

I have tried with define profile UDF using cell loop but it is resulting in acess violation error?

thanks

[Amir]

Hi,
when you want to use gradients, you have to prevent the solver from freeing up memory by issuing the text command solve/set/expert and then answering yes to the question Keep temporary solver memory from being freed?
if it doesn't work, post your UDF and I'll correct it for you.

[rakadit]

Thanks for your suggestion! The UDF which I had written for calculating wall shear is like thislease correct the UDF and send the corrected one to me.

#include "udf.h"
DEFINE_PROFILE(stressx,t,i)
{
cell_t c;
real dvdy_r;
real alpha=0.7;
real dvdy=C_DVDY(c,t);
real mu_a=C_MU_L(c,t);
dvdy_r=alpha*dvdy;
begin_c_loop(c,t)
{
F_PROFILE(c,t,i)=mu_a*dvdy_r ;
}
end_c_loop(c,t)
}

[Amir]

It seems that you need to focus on UDF concepts more carefully.

#include "udf.h"
DEFINE_PROFILE(stressx,t,i)
{
face_t f;
cell_t c;
Thread *tc;
real dvdy_r;
real alpha=0.7;
real dvdy;
real mu_a;
begin_f_loop(f,t)
{
    c=F_C0(f,t);
    tc=THREAD_T0(t);
    mu_a=C_MU_L(c,tc);
    dvdy=C_DVDY(c,tc);
    dvdy_r=alpha*dvdy;
    F_PROFILE(f,t,i)=mu_a*dvdy_r;
}
end_f_loop(f,t)
}

Regards,

Amir

[rakadit]

Dear Mr. Amir,

Thanks for correcting UDF.Yes,you are right I am new to UDF writing .I will soon give u feedback on corrected UDF.

with regards,


Rakesh

[rakadit]

Dear Mr. Amir,

Sorry for late response.I had gone to some remote place.I am working with the UDF.it is working fine.Results are coming but correctness of results is to be confirmed.
Thanks again!


Rakesh