volScalarField& e = thermo.he();
const volScalarField& psi = thermo.psi();