[mnolan93]

I am having issues with Population Balance setup. I have three fluids - dissolved paracetemol, IPA, and water. Depending on concentrations in each cell, the nucleation of paracetemol particles will be present. However, how do I get the particles created by the population balance to be pure paracetemol particles? I am using a discrete model with 2 phases (fluids and solid particles). I have written UDFs describing the crystallisation kinetics.

My UDF is also having problems. I am unsure how to fix these. My code is the following:

/************************************************** **********************
UDF for PopBal
************************************************** ***********************/

#include "udf.h"


DEFINE_PB_NUCLEATION_RATE(nuc_rate, cell, thread)
{
real J,C,Cstar,x;
real bo = 6829200; /* constant b0 */
real bi = 149900; /* constant b1 */
real rho = 1260; /* density of paracetemol crystals */
real solute_mass_frac,solvent_mass_frac,antisolvent_mas s_frac;


Thread *tc = THREAD_SUPER_THREAD(thread); /*obtain mixture thread */
Thread **pt = THREAD_SUB_THREADS(tc); /* pointer to sub_threads */
Thread *tp = pt[P_PHASE]; /* primary phase thread */

solute_mass_frac = C_YI(cell,tp,0); /* solute mass fraction */
solvent_mass_frac = C_YI(cell,tp,1); /* solvent mass fraction */
antisolvent_mass_frac = C_YI(cell,tp,2); /* antisolvent mass fraction */


x = antisolvent_mass_frac/(antisolvent_mass_frac+solvent_mass_frac);

Cstar = (0.000001*x*x*x)-(0.0002*x*x)+(0.0097x)+0.1087; /* solubility */

C = solute_mass_frac/solvent_mass_frac; /* concentration */

J = bo*exp(-bi*((log(rho/Cstar))*(log(rho/Cstar)))/(log(C/Cstar)));

return J;

}

DEFINE_PB_GROWTH_RATE(growth_rate, cell, thread)
{
real G,Ko,Ki,C,Cstar,deltaC,x;
real ka = 267; /* constant */
real kb = 89; /* constant */
real kc = -139; /* constant */
real kd = 917500; /* constant */
real ke = 9919300; /* constant */
real solute_mass_frac,solvent_mass_frac,antisolvent_mas s_frac;

Thread *tc = THREAD_SUPER_THREAD(thread); /*obtain mixture thread */
Thread **pt = THREAD_SUB_THREADS(tc); /* pointer to sub_threads */
Thread *tp = pt[P_PHASE]; /* primary phase thread */

solute_mass_frac = C_YI(cell,tp,0); /* solute mass fraction */
solvent_mass_frac = C_YI(cell,tp,1); /* solvent mass fraction */
antisolvent_mass_frac = C_YI(cell,tp,2); /* antisolvent mass fraction */


x = antisolvent_mass_frac/(antisolvent_mass_frac+solvent_mass_frac);

Cstar = (0.000001*x*x*x)-(0.0002*x*x)+(0.0097x)+0.1087; /* solubility */

C = solute_mass_frac/solvent_mass_frac; /* concentration */

Ko = (ka*x*x)+(kb*x)+kc;

Ki = (kd*x)+ke;

deltaC = C-Cstar

if (deltaC <= 0)
{
G = 0;
}
else;
{
G = Ko*pow((deltaC),Ki);
}
return G;
}


H:\\My Documents\\ARP\\UDF\\UDF.c:38: macro `DEFINE_PB_GROWTH_RATE' used with only 3 args
Error: H:\\My Documents\\ARP\\UDF\\UDF.c: line 28: parse error.
Error: H:\\My Documents\\ARP\\UDF\\UDF.c: line 59: parse error.
Error: H:\\My Documents\\ARP\\UDF\\UDF.c: line 69: parse error.
Error: H:\\My Documents\\ARP\\UDF\\UDF.c: line 73: parse error.

Any help is appreciated.

Thanks

[mnolan93]

anyone able to help?

cheers