[iiigoriii]

Dear OpenFOAM users,

I made my own reaction rate function (in according to hints from this forum) and everything works perfectly. I have total match between manual calculation, DWSIM and OpenFOAM. If somebody will be interested in this in future, little hint -> check units. It is biggest mind****.

I created reaction rate as polynomial, dependent on temperature. I would like to add also mass fraction dependence. My reaction rate is created from LangmuirHinshelwoodReactionRate. This reaction rate can use variables: pressure (P), temperature (T) and molar concentration (c).

For example in StandardChemistryModel.C i see for molar mass, but I don't know how to define specieThermo_.

scalar W = specieThermo_[i].W()

Is there a simple (or relative) way to add mass fraction "Y" to LangmuirHinshelwoodReactionRateI.H ? (or actual "rho"? from "c" it is simple to calculate "Y"...)

thanks, BR, ii