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

Interpolation #104

Open
ananthanarasimhanj opened this issue May 24, 2023 · 1 comment
Open

Interpolation #104

ananthanarasimhanj opened this issue May 24, 2023 · 1 comment

Comments

@ananthanarasimhanj
Copy link

ananthanarasimhanj commented May 24, 2023

For input of 3 Td, interpolation is used to get rate coefficient from the reported data. Below are the nearby points.

E/N--------- k
2.848036------5.344192e-31
3.274549 ------7.650413e-28

ISSUE:
It is expected that, k should be between e-31 and e-28 for 3 Td. However, when “spline" interpolation was used, the chosen k value for 3 Td was1.18e-24 (spline.i file in in "Interpolation.zip").

ALTERNATE TRIED:
When interpolation was changed to “linear,” the chosen k value was 9.36e-29 (linear.i in "Interpolation.zip").

Interpolation.zip

In another case, with simple Two reaction for N2, for 30 Td, the "spline" picked k as 0, whereas the linear worked fine (N2-Interpolation.zip for reference)

N2-Interpolation.zip

@smpeyres smpeyres added the bug Something isn't working label May 24, 2023
@dcurreli
Copy link
Contributor

dcurreli commented Sep 5, 2023

I suggest to use a linear interpolator in this case. If you really need to spline exponential values (even if I do not think this is the case for this problem), I suggest to add use_log = true in [ChemicalReactions].

@dcurreli dcurreli removed the bug Something isn't working label Sep 5, 2023
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

3 participants