eta_phi d(mfmm1.mumixture,phid) Derivative of viscosity with respect to concentration
J1 -rho_s*(phid*D1*(phid*gamma_dotx+phidx*gamma_dot)+phid*phid*gamma_dot*D2*eta_phi*phidx/mfmm1.mumixture) Particle flux, x direction
f_hinder eta_f*(1-phi0)/mfmm1.mumixture Hindering function
J2 -rho_s*(phid*D1*(phid*gamma_doty+phidy*gamma_dot)+phid*phid*gamma_dot*D2*eta_phi*phidy/mfmm1.mumixture)+rho_s*f_hinder*phid*u_st Particle flux, y direction
cd rho_s*phid/mfmm1.rho Mass fraction of dispersed phase