Equilibrating a mixture that contains an interface phase fails. I don't think there's any fundamental reason why such phases shouldn't be able to be included in a multiphase equilibrium calculation.
File ~/src/cantera/build/python/cantera/mixture.pyx:330, in cantera.mixture.Mixture.equilibrate()
328 successively more verbose information.
329 """
--> 330 self.mix.equilibrate(stringify(XY.upper()), stringify(solver), rtol,
331 max_steps, max_iter, estimate_equil, log_level)
CanteraError:
*******************************************************************************
CanteraError thrown by Phase::assignDensity:
density must be positive. density = nan
*******************************************************************************
Problem description
Equilibrating a mixture that contains an interface phase fails. I don't think there's any fundamental reason why such phases shouldn't be able to be included in a multiphase equilibrium calculation.
Steps to reproduce
Behavior
System information
main
branch at 4565a55e4Additional context
Originally reported on the Users' Group.