From 2a69d99d8320e3073d9b0e5c350af0d34e8017ec Mon Sep 17 00:00:00 2001 From: Ray Speth Date: Sun, 7 Sep 2014 21:17:05 +0000 Subject: [PATCH] [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. --- src/oneD/StFlow.cpp | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) diff --git a/src/oneD/StFlow.cpp b/src/oneD/StFlow.cpp index 56c621ebb..73bfd8f4f 100644 --- a/src/oneD/StFlow.cpp +++ b/src/oneD/StFlow.cpp @@ -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); }