Chemical reaction didn't make sense

MFiX
#1

Hi all,

I tried to run a reaction of sorbent being heated and releases CO2 and freeing surface sites. Namely PR → 2 ML + CO2.

In my usr_rates.f file, the reaction is linear proportional to the PR concentration. However, even when the PR concentration is very low, there are still a lot of CO2 being output. Meanwhile, the CO2 measured at the outlet was almost constant, but the PR concentration did not decrease at all. Did I misunderstand how the reaction rate work?

I attached my code in the thread.
Stack_Desorp (5).mfx (15.7 KB)
usr_rates (6).f (2.8 KB)

Thanks,
Jack

#2

I think you need to use the bulk density ROP_s instead of the density RO_s when you compute c_PR :

c_PR  = (ROP_s(IJK,1)*X_s(IJK,1,PR)/MW_s(1,PR))
2 Likes
#3

Thanks Jeff! It’s working now.