[lehoanganh07]

Hi, I posted one thread about this problem. Now I solved my problem and could compile to Fluent. This UDF file used to midify the turbulent viscosity for k epsilon model by Reboud formular. when I ran solution, althout, not diver but the residua become unstable which is reach to 10^3 for continuity residual. So I think have some thing wrong in my code, but I cant find out.
someone please help me! Thank you very much!


#include "udf.h"
#include "math.h"

#define pc 1944610000.0
#define t0 3836.94
#define ck 472.26
#define R 8.31

DEFINE_TURBULENT_VISCOSITY(reboud_mut,cell,thread)
{
real mut;
real rhol;
real rhov;
real cmu = 0.09;
real damp;
real voidf;
real k =C_K(cell,thread);
real d =C_D(cell,thread);

int liquid=0;
int vapor =1;

Thread *phase1;
Thread *phase2;

phase1=THREAD_SUB_THREAD(thread,liquid);
phase2=THREAD_SUB_THREAD(thread,vapor );

rhol=C_R(cell,phase1);
rhov=C_R(cell,phase2);
voidf=C_VOF(cell,phase2);

damp=rhov+(rhol-rhov)*pow((1.0-voidf),10);
mut=damp*cmu*k*k/d;

return (mut);
}