[anubhavd]

Dear All

I am writing a udf in which I need to call the volume fraction of primary and secondary phase to calculate some quantity. My case has two phases and I am using Eulerian multi phase. Please help in calling the Volume fractions I am attaching the body of the written UDF. The part which seems confusing to me is marked as red in the attached udf. Is the attached code is correct??

#include "udf.h"
#include "sg_pb.h"
#include "sg_mphase.h"

DEFINE_PB_NUCLEATION_RA_rate_J, cell, thread)
{


double J, rho, A, mol_wt, P, K, C, VF_w, VF_b, T;


Thread *tc = THREAD_SUPER_THREAD(thread);
Thread **pt = THREAD_SUB_THREADS(tc);
Thread *tp = pt[P_PHASE];


T = 300;
mol_wt = 0.018; /* ......unit 'Kg/mol'*/
rho = 1000; /* ......unit 'Kg/m^3)*/
A = 6.023*pow(10,23); /* ......dimensionless*/
K = 1.3807*pow(10,-23); /* ......J K-1*/

VF_w = C_VOF(cell,pt[0]);

VF_b = C_VOF(cell,pt[1]);

C = (rho * VF_w* A/mol_wt);
P = C*K*T;


J = 150000;

C_UDMI(cell, thread, 0) = VF_w;
C_UDMI(cell, thread, 1) = VF_b;
C_UDMI(cell, thread, 2) = C;
C_UDMI(cell, thread, 3) = P;


return J;

}