+ remove initial equation (initial purpose was to help
                initialization by provided inverse function of homotopy
        + enable user to set T/deltaP ref used to compute damper Kv
        + fix homotopy simplified expr (to compute mass flow based on
            Kv)