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

Derivatives of Net Production Rates w.r.t. Temperature #129

Open
koendev opened this issue Jun 3, 2020 · 2 comments
Open

Derivatives of Net Production Rates w.r.t. Temperature #129

koendev opened this issue Jun 3, 2020 · 2 comments

Comments

@koendev
Copy link

koendev commented Jun 3, 2020

Hi,

The Kinetics header file provides the Jacobian of the species net production rates w.r.t. the species densities, which is very useful. Does something similar exist for the derivatives w.r.t. the temperature? If not, are there any plans to implement this in the future?

Thank you for your help!
Koen Devesse

@grgbellasvki
Copy link
Collaborator

Dear Koen,

The derivatives of the Arrhenius kf with respect to temperature is already available in the code for the one-temperature model. You can find it in RateLaws.h:77 and the function is called derivative. If this is what you need, let us know if you need help on how to access this function.

For multi-temperature model it is not available, but it may be available in the future.

@koendev
Copy link
Author

koendev commented Jul 30, 2020

Dear,

Thank you for your reply! I'm sorry it took me so long to respond.

Having the derivatives for the kf would indeed be enough for me to do what I need! If you could help me on how to use that function, that would be great, as I don't understand what the different inputs are. I guess the kf themselves can be found through the functions forward- and backwardRateCoefficients, with multidimensional double arrays as inputs, correct?

Finally, is there also a way to easily access all the stochiometric coefficients for all the species and all the reactions? I would also need those, but I don't see an easy function that outputs those.

Thank you very much!
Koen Devesse

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

No branches or pull requests

2 participants