I want to write a UDF for heat transfer coefficient which changes as a function of the surface temperature as follows:

htc=24.1e-4*Epsilon*pow(T, 2.0).

So, anyone tell me how to write the UDF with which DEFINE macros?

Thanks in advance!

#include "udf.h"

DEFINE_PROFILE(pressure_profile,t,i) {

face_t f; double htc; real Epsilon=10.;

begin_f_loop(f,t)

{

htc=24.1e-4*Epsilon*pow(F_T(f,t), 2.0);

F_PROFILE(f,t,i) = htc;

} end_f_loop(f,t) }

Thank you for your reply!