This Gibbs energy loss is linear in the molarflow with the linear factor kC:
du := n_flow/kC;
function linearPotentialDiff
  extends partialPotentialDiff;
  input Real kC(unit = "mol2/(J.s)") = 1 "Linear factor" annotation(
    Dialog(enable = true));
end linearPotentialDiff;