DEFINE_DPM_HEAT_MASS source terms
 

Hi,

I'm trying to write my own mass transfer law for a multicomponent particle injection, but the source terms in the DEFINE_DPM_HEAT_MASS macro really do not do what I want. I've tried to keep it simple by writing a UDF that just transfers mass from the DPM particle to the fluid phase:

#include "udf.h"

/* Species: * Particle: [0]=fdp [1]=plasminogen * Blood: [0]=injected_tpa [1]=fdp_blood [2]=plasma */

/* Source Term Variable Unit

* particle temperature dydt[0] K/s

* particle component mass dydt[1..] kg/s

* gas phase enthalpy dzdt->energy J/s

* gas phase species mass dzdt->species[0..] kg/s */

DEFINE_DPM_HEAT_MASS(simple_source,p,Cp,hgas,hvap, cvap_surf,dydt,dzdt) {

int ns;

int nc = TP_N_COMPONENTS(p); /* number of particle components (2) */

for( ns = 0; ns < nc; ns++ )

{

dydt[1+ns] -= 0.000811578102;

dzdt->species[1+ns] += 0.000811578102;

} }

The evaporating species of fdp is fdb_blood and plasminogen should evaporate (be transferred) to plasma. But for some reason the DPM mass source in the simulation is two magnitudes smaller (2.0344e-06), the particle diameter does not at all getting smaller in the way it should be, and in the contour plot of the species I can't see that any mass is added to my fluid phase. Any ideas what could be wrong??

Thank you very much! Sebastian