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	     2*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	     Oc12*n2*n1+2*Oc22*n2*n2+Oc23*n2*n3	       Death rate of class size 2 due to coalesncence  
Bc2	     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*Oc12*n2*n1+5/4*Oc13*n1*n3+Oc22*n2*n2	                                                 
Dc4	     0	                                                                                        
Bc4	     3/4*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 for all discrete phase, m  
