[eddiegolpesar]

Hi
I prepared the following UDF for Darcy term in DC casting, for velocity in x direction. It was interpreted successfully. I want use some other UDFs like this for velocity in y direction, k and epsilon. But I see only the last interpreted UDF and I can't see all of them in the graphical panel in source terms part. Please help me....

/*EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE EEEEEEEEEEEEEEEEE

UDF that adds Darcy source term for U in the X direction

DDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD DDDDDDDDDDDDDDDDD*/

#include "udf.h"

#define FST 0.3
#define TL 903
#define TS 743
#define TF 933
#define kp 0.16
#define shi 1.0e-15
#define visc 0.001
#define K0 2.0e-11

real TEMP, DEPVAR, SDEPVAR, PERM, FL, FS;
/*real DEPVAR;
real SDEPVAR;
real PERM;
real FL;
real FS;*/

DEFINE_SOURCE(darcy_ux, c, t, dS, eqn)
{
real source;

TEMP= C_T(c,t);
DEPVAR= C_U(c,t);
SDEPVAR= 0. ;

/***** cacculation of FL & FS *****/
if (TEMP >= TL)
{
FL=1. ;
}
if ((TEMP > TS) && (TEMP < TL))
{
FL= 1.-1./(1.-kp)*(TL-TEMP)/(TF-TEMP);
}
if (TEMP <= TS)
{
FL= 0. ;
}

/***** cacculation Permeability & Darcy term *****/
if (FS==0.)
{
source = dS[eqn] = 0. ;
}
if ((FS>0.) && (FS<=FST))
{
source = dS[eqn] = 0. ;
}
if ((FS>=FST) && (FS<1.))
{
PERM=K0*pow((1.-FS),3.)/pow(FS,2.);
source = -visc/(PERM+shi)*DEPVAR + visc/(PERM+shi)*SDEPVAR ;
dS[eqn]= visc/(PERM+shi) ;
}
return source;
}