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;