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+Oc11*n1*n1	             Death rate of class size 1 due to coalesncence  
Bc1	     0	                                                       Birth rate of class size 1 due to coalescence   
Dc2	     Oc12*n2*n1+Oc22*n2*n2+Oc23*n2*n3+Oc22*n2*n2	             Death rate of class size 2 due to coalesncence  
Bc2	     0.5*Oc12*n1*n2+Oc11*n1*n1	                               Birth rate of class size 2 due to coalescence   
Dc3	     Oc13*n1*n3+Oc23*n2*n3+2*Oc33*n3*n3	                                                                      
Bc3	     3/4*Oc13*n3*n1+0.5*Oc12*n2*n1+0.5*Oc23*n3*n2+Oc22*n2*n2	                                                 
Dc4	     0	                                                                                                       
Bc4	     0.25*Oc13*n3*n1+0.5*Oc23*n3*n2+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, m                         
