Dear altruists,
I am trying to simulate a transient fluent model with a wall heat source. When I had a udf for heat source, the solution becomes divergent, but when I put constant value, the simulation works. Is there any recommended solution?
#include "udf.h"
DEFINE_PROFILE(Frictional_heat_flux,thread, index)
{
real x[ND_ND];
real mue=0.3;
real omegaV;
real temp;
real sticking=0.65;
real theta;
real R;
real pressure=40e06;
real yieldstress;
real cf=0.95;
real U=0.004656;
real omega=32.5;
face_t f;
begin_f_loop(f, thread)
{
F_CENTROID(x,f, thread);
temp=F_T(f, thread);
if (temp