[SANDIPSINH JODDHA]

#include "udf.h"

#include "mem.h"



#define Q 145000 //J/mol

#define A 8.86e+6

#define n 3.55

#define R 8.314 //J/molK

#define alpha 0.045

DEFINE_PROPERTY(cell_viscosity, c, t)

{

real temp = C_T(c, t);

real mu_lam;

real Z;

real Flow_stress;

real x;

real strain;

if (strain < 0)

{

strain = 10;

}

else


{
strain = C_STRAIN_RATE_MAG(c, t);

}


Z = (strain)*exp(Q / (R*temp));

x = pow((Z / A), (1. / n));

Flow_stress = (1. / alpha)*(log(x + sqrt(pow(x, 2) + 1.)));

mu_lam = Flow_stress / (3.*strain);

if (mu_lam > 0)

return mu_lam;

}

we are trying to define the viscosity in term of strain rate and temperature using UDF code but, we are not getting value as expected.

[AlexanderZ]

you should specify what exactly is different from your expectations
compile the code and you will get following output:

warning C4715: 'cell_viscosity' : not all control paths return a value
warning C4700: uninitialized local variable 'strain' used

fix that

take into account that DEFINE_PROPERTY macro is a loop over all cells in domain, but your variable strain contains only 1 value for the last cell, it overwrites each other cell. You may use used defined memory to store strain on each cell

Also I would initialize strain for the 0-time moment

[SANDIPSINH JODDHA]

Can you please show us how to do that because I am new to this and when I try this code every cell get's same viscosity value.