-
-
Notifications
You must be signed in to change notification settings - Fork 346
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
Fix some 1D solver errors #1103
Conversation
Codecov Report
@@ Coverage Diff @@
## main #1103 +/- ##
=======================================
Coverage 73.49% 73.49%
=======================================
Files 364 364
Lines 47809 47808 -1
=======================================
Hits 35138 35138
+ Misses 12671 12670 -1
Continue to review full report at Codecov.
|
Fixed an error in the error handler for singular Jacobians where the wrong domain was being identified as containing the problematic row, that ultimately resulted in an integer divide by zero error instead of a CanteraError identifying which domain, grid point, and solution component the row corresponded to. This fixes a regression introduced in e8194af.
15e368a
to
7267fed
Compare
A reasonable initial guess for lambda can be obtained by noting the proportionality between lambda and rho*V**2 in the radial momentum equation. This fixes some cases where the correct solution involves large values of lambda and the initial guess of lambda = 0 was leading to singular Jacobians. See Cantera#1067
7267fed
to
f51f58a
Compare
This is analogous to the issue with multicomponent transport / Soret diffusion that was reported in Cantera#615 and fixed in Cantera#631.
8bb2141
to
f4dff38
Compare
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks good to me - thanks for tracking down the issue!
Changes proposed in this pull request
CanteraError
identifying which domain, grid point, and solution component the row corresponded to. This fixes a regression introduced in e8194af.FreeFlame
If applicable, fill in the issue number this pull request is fixing
Resolves #1067
If applicable, provide an example illustrating new features this pull request is introducing
Checklist
scons build
&scons test
) and unit tests address code coverage