#include udfh#include mathh#include sg_mphasehDomain D = Get_Domain1;Thread pt;Thread cell_thread;mp_thread_loop_ccell_threadDpt cell_t c; real solid_viscosity; begin_c_loopccell_thread

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

DEFINE_ON_DEMAND(cell_thread_loop) { Domain *D = Get_Domain(1); Thread **pt; Thread *cell_thread;

cell_thread = Lookup_Thread(D, 1); // assuming the cell thread ID is 1

thread_loop_c(pt, cell_thread)
{
    cell_t c;
    
    begin_c_loop(c, cell_thread)
    {
        real solid_viscosity;
        
        real liquid_viscosity = C_MU_L(c, cell_thread);
        real vof_liquid = C_VOF(c, pt[0]);
        real vof_solid = C_VOF(c, pt[1]);
        real shear_rate = C_STRAIN_RATE_MAG(c, cell_thread);
        
        real mix_viscosity = 3.383*pow(shear_rate,-0.3708)+2.64;
        
        solid_viscosity = (mix_viscosity-vof_liquid*liquid_viscosity)/vof_solid;
        
        Message("Solid Viscosity at cell %d: %f\n", c, solid_viscosity);
    }
    end_c_loop(c, cell_thread);
}