Modeling of the Ion Exchange in Pulp Suspensions by Gibbs Energy Minimization