research/cantera
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Fixed a problem with the VCS equilibrium solver at constant HP or UP

Browse source 
A sign error when limiting the temperature step caused the solver to fail if
the initial temperature was higher than sqrt(Tlow*Thigh).
This commit is contained in:
Ray Speth 2012-10-12 20:34:59 +00:00
parent e9026fb683
commit a6dc994439
1 changed files with 1 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
src/equil/vcs_MultiPhaseEquil.cpp
View file

@ -294,7 +294,7 @@ int vcs_MultiPhaseEquil::equilibrate_HP(doublereal Htarget,
Tnew = sqrt(Tlow*Thigh);
dT = Tnew - Tnow;
if (dT < -200.) {
dT = 200;
dT = -200;
}
if (dT > 200.) {
dT = 200.;