Oc11	    2.62e-7 [m^3/s]	                               Coalescence Rate                                
Oc12	    3.48e-7[m^3/s]	                                Coalescence Rate between class 1 and 2          
Oc13	    4.76e-7[m^3/s]	                                                                                
Oc14	    6.7e-7[m^3/s]	                                                                                 
Oc22	    4.49e-7[m^3/s]	                                                                                
Oc23	    5.96e-7[m^3/s]	                                                                                
Oc24	    8.16e-7[m^3/s]	                                                                                
Oc33	    7.69e-7[m^3/s]	                                                                                
Oc34	    1.02e-6[m^3/s]	                                                                                
Oc44	    1.32e-6[m^3/s]	                                                                                
Dc1	     Oc11*n1*n1+Oc12*n1*n2+Oc13*n1*n3	              Death rate of class size 1 due to coalesncence  
Bc1	     0	                                             Birth rate of class size 1 due to coalescence   
Dc2	     0.5*Oc12*n2*n1+Oc22*n2*n2+Oc23*n2*n3	          Death rate of class size 2 due to coalesncence  
Bc2	     0.5*Oc11*n1*n1	                                Birth rate of class size 2 due to coalescence   
Dc3	     0.25*Oc13*n1*n3+0.5*Oc23*n2*n3+Oc33*n3*n3	                                                     
Bc3	     0.5*Oc12*n2*n1+0.5*Oc22*n2*n2	                                                                 
Dc4	     0	                                                                                             
Bc4	     0.25*Oc13*n3*n1+0.5*Oc23*n3*n2+0.5*Oc33*n3*n3	                                                 
R1	      Bc1-Dc1	                                                                                       
R2	      Bc2-Dc2	                                                                                       
R3	      Bc3-Dc3	                                                                                       
R4	      Bc4-Dc4	                                                                                       
v_total	 n1*v1+n2*v2+n3*v3+n4*v4	                       Total volume of discrete phase                  
n_total	 n1+n2+n3+n4	                                   total number of particle density                
m_total	 v_total*rhod	                                  total mass of discrete phase                    
a_total	 n1*a1+n2*a2+n3*a3+n4*a4	                       total area of discrete phase                    
d32	     6*v_total/a_total/n_total	                     mean Sauter diameter for all discrete phase, m  
