u_0 Q_0/((mfmm1.rhoc*(1-phid_0)+mfmm1.rho_s1*phid_0)*pi*r_pipe^2) Inflow velocity
d_ave max(d1,(d1*s1+d2*s2+d3*s3+d4*s4+d5*s5)/(s1+s2+s3+s4+s5)) Average droplet diameter
d_max 1[m]*((sigma/mfmm1.rhoc)*1[s^2/m^3])^0.6/(max(ep,eps)*1[s^3/m^2])^0.4 Maximum stable droplet diameter
R2 mfmm1.rho_s2*spf.gammaT*(s2/d2^3*abs(d2^3-min(d2,d_max)^3)) Break-up rate
R3 mfmm1.rho_s3*spf.gammaT*(s3/d3^3*abs(d3^3-min(d3,d_max)^3)) Break-up rate
R4 mfmm1.rho_s4*spf.gammaT*(s4/d4^3*abs(d4^3-min(d4,d_max)^3)) Break-up rate
R5 mfmm1.rho_s5*spf.gammaT*(s5/d5^3*abs(d5^3-min(d5,d_max)^3)) Break-up rate