[1D] Update thermo properties during Jacobian updates

In combination with the previous commit, this significantly improves the
convergence behavior of the solver when using multicomponent transport
properties. In many cases, the solver is now able to directly solve the
multicomponent problem directly from the mixture-averaged solution without any
timestepping.
This commit is contained in:
Ray Speth 2014-09-07 21:17:05 +00:00
parent 92f8fc4292
commit 2a69d99d83

View file

@ -261,10 +261,9 @@ void StFlow::eval(size_t jg, doublereal* xg,
// update properties
//-----------------------------------------------------
// update thermodynamic and transport properties only if a Jacobian is not
// being evaluated
updateThermo(x, j0, j1);
// update transport properties only if a Jacobian is not being evaluated
if (jg == npos) {
updateThermo(x, j0, j1);
updateTransport(x, j0, j1);
}