problem in transcient heat flux b.c in UDF
i am simualating solidification and melting problems in the fluent6.0,so for that i have to give heat flux boundary condition varing with time. So i have written the Udf for that but it is diverging after time (15sec) it is diverging in the third condition of the udf.So could u plz help me
#include "udf.h" #define d 0.04 #define C 2.18844 #define D -0.602603 #define A 762.048e04 #define B -0.217412 #define den 8960 #define kc 363 #define cp 423 DEFINE_PROFILE(heat_flux, thread,position) { real alpha, qmax,q; face_t f; alpha=kc/(den*cp); qmax=A*pow(d/alpha,B); begin_f_loop(f,thread) { real t=RP_Get_Real("flow-time"); if(t<=5) F_PROFILE(f,thread,position) = qmax; else if((t>5)&&(t<=15)) F_PROFILE(f,thread,position) = qmax-6.5715e4*(t-5); else F_PROFILE(f,thread,position) =- 2.0152e6*C*pow(t,-D)/(pow(alpha,0.05)); } end_f_loop(f,thread) } |
Re: problem in transcient heat flux b.c in UDF
hello
Couple of mistakes All the number must to have the form 5. or 5.0 not the form 5 in order to assign the decimal number. also you don't need to write this sentences like this else if((t>5)&&(t<=15)) you can write the sentes like this else if(t<=15.0) F_PROFILE(f,thread,position) = qmax-6.5715e4*(t-5.0); Best Regards Alex Munoz |
Re: problem in transcient heat flux b.c in UDF
Hai Alex,
thank's for responding,i did what u said, again the solution is diverging at the third condition i.e [F_PROFILE(f,thread,position) =- 2.0152e6*C*pow(t,-D)/(pow(alpha,0.05));] when the fluent is solving the above equation, and the solution is get diverging. |
Re: problem in transcient heat flux b.c in UDF
Hi
It is correct that you define D=-### and then in your function you writhe pow(t,-D)? I wonder if you are making -*-=+ In addition how do you know that fluent is diverging in the third condition, the code is not able to tell you that type of information. Just it has divergence but how you can ensure that it is due to the UDF?? Regards Alex |
All times are GMT -4. The time now is 05:40. |