Much of the algorithm for relaxing the element potentials
to equilibrium employs the Brinkley algorithm (const T, P)
(p. 129 Smith and Missen). This algorithm does better at
relaxing the system towards steady state from a bad initial
guess.
The old algorithm now does the last little bit and it does
the calculation for fixed pairs other than const T, P, using
the Brinkley algorithm's initial guess.
Modified the damping algorithm in the old method. It now doesn't
create false convergence failures when it is close to the
end of the calculation.
The new solver has been tested against a matrix of initial conditions
for the gri30 mechanism and a 200 species jp8 mechanism.
Yet to do:
This routine should work but has not been checked out against:
phases with electric charge
nonideal phases.
This routine should probably not yet work with
molality based phases.
More work may have to be done with the Brinkley algorithm for
problems which suffer from stoichiometric constraints. A full
Gauss-Jordon elimination step should probably be added to
flag cases where there are effective stoichiometric constraints.
solver.
This test does a matrix of T, pres, and mixture fraction for the
gri30 mechanism. It does a total of 1500 calculation starting from
a bad initial guess.
it sums to one. Previously, it was set uniformly to zero. The mass
fraction of species 0 is set to 1.
This was noticed by a user, and may have resulted in seg faults before
the ThermoPhase object was properly initialized.
These are needed in a few places so make doesn't do any extra work
when it doesn't need to.
Changed a few of the "install multiple files at once" cases into
installations where one file is installed at a time. Some install
commands, including the one that comes with cantera can only
install one file at a time.
