Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Formulas used to calculate kf #28

Open
wittregr opened this issue May 12, 2020 · 3 comments
Open

Formulas used to calculate kf #28

wittregr opened this issue May 12, 2020 · 3 comments
Assignees
Labels
question Further information is requested

Comments

@wittregr
Copy link
Contributor

What are the formulas used to calculate kf? Please add this information to the documentation for kf.out and include the units kf will have for each formula.

@MaxRCohen
Copy link

Hi @mbkumar, looks like we still need this information. Could you please add it where the documentation says "1. Adsorption reaction:" and "2. Surface reaction:"?

@skasiraj skasiraj added the question Further information is requested label Oct 18, 2021
@skasiraj skasiraj self-assigned this Oct 19, 2021
@skasiraj
Copy link
Contributor

The formulas used for gas and surface reactions are implemented by Cantera using the Arrhenius and Modified-Arrhenius forms. If this information is consistent with your experience we can update the documentation and close this issue.

  • For surface reaction without sticking coefficients it uses this expression. Here A, b, and E_a are the modified Arrhenius parameters. If activation energy is defined as coverage dependent then a_k, mk, Ek, are the coverage dependencies from species k. For pMuTT generated input I am not sure if we ever use this form. If we do not provide values for a_k. m_k and E_k, then they will be set to a default of 0 and the expression reduced to the typical modified Arrhenius form.
    image
  • For sticking type reactions this is the expression for the forward rate constant. Where gamma_tot is the total molar site density, m is the sum of all the surface reactant stoichiometric coefficients, and W is the molecular weight of the gas phase species
    image
  • The The unit system is (kmol, m, s). The actual units depend on the reaction order and the dimensionality (surface or bulk).

@MaxRCohen
Copy link

To confirm, is the surface reaction term A equal to (k_BT/h)(1/gamma_tot^(m-1))? k_B is the Boltzman constant, h is Planck's constant, and the other terms are the same as in the adsorption section. For the adsorption section, is the lower case gamma equal to the sticking coefficient? If so, then these formulas are consistent with what I know for Chemkin and should be included in the documentation. However, @wittregr should confirm the coverage dependent formulation for the surface reactions since I less familiar with it.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
question Further information is requested
Projects
None yet
Development

No branches or pull requests

3 participants