[JeroenH]

I am trying to write a UDF to determine the density as function of composition and temperature using the DEFINE_DPM_PROPERTY udf. Now I am trying to access the species index by name (as I would in a gas mixture) however this does not work for the combusting type DPM particle.

Any advise on how to do this? I have given an example of my current code below:

DEFINE_DPM_PROPERTY(rho_particle, c, t, tp, T)
{
Material *cond_mix = TP_MATERIAL(tp);
int ind_A = mixture_specie_index(cond_mix, "a"); /* specie index A */
int ind_B = mixture_specie_index(cond_mix, "b"); /* specie index B */

int ns;
Material *sp;
real mf[MAX_SPE_EQNS];
mixture_species_loop(cond_mix, sp, ns)
{
mf[ns] = TP_COMPONENT_I(tp, ns); /* component mass fraction */
}

rho_A = ...
rho_B = ...

real rho = 1/(mf[ind_A]/rho_A + mf[ind_B]/rho_B)

return rho;
}


Thank you!