[perun333]

Hello guys,

I have problem when Im starting my computing with udf (fatal error). But compiling without problems. I have 2D axisymmetric model, where I want to keeping time and space temperature gradients only in some part of the 2D model. Here is my code:

#include "udf.h"

DEFINE_EXECUTE_AT_END(gradienty)
{
Domain *domain;
cell_t c;
Thread *ct;
face_t f;

FILE *fw1;
FILE *fw2;
int ID = 1;

real gradient_prostorovy;
real gradient_casovy;


gradient_prostorovy = 0.;
gradient_casovy = 0.;
domain = Get_Domain(1);


ct = Lookup_Thread(domain,ID);

begin_f_loop(f,ct){

{
gradient_prostorovy = MAX(gradient_prostorovy,C_UDMI(c,ct,0) = sqrt(C_T_G(c,ct)[1]*C_T_G(c,ct)[1]));


gradient_casovy = MAX(gradient_casovy,C_UDMI(c,ct,1) = 3600 * (C_T(c,ct) - C_T_M1(c,ct)) / CURRENT_TIMESTEP);

}
end_f_loop(f,ct)
}

fw1 = fopen("gradient_prostorovy.txt","w");
fprintf(fw1,"%lf ", CURRENT_TIME);
fprintf(fw1,"%lf\n", gradient_prostorovy);
fclose(fw1);

fw2 = fopen("gradient_casovy.txt","w");
fprintf(fw2,"%lf ", CURRENT_TIME);
fprintf(fw2,"%lf\n", gradient_casovy);
fclose(fw2);

}


PS: For next step I would like to make udf for changing heat flux on the wall. If will be gradient_casovy (time_gradient) more then 19 I would like to use heat flux - 10 else +10.

So I would like to do it by this way

DEFINE_PROFILE(muj_profil,t,i)
{
face_t f;

begin_f_loop(f,t)
{

if( gradient_casovy > 19. ) {
F_PROFILE(f,t,i)=.....-10.;
}

else {
F_PROFILE(f,t,i)=....+10.;
}

}
end_f_loop(f,t)
}

I would be happy for help