4. METHOD OF SOLUTION
The numerical technique used is based upon dividing the particle size domain into m sections and imposing the condition of mass conservation for each chemical component for the processes considered. Aerosol mass concentrations are grouped into sections (i.e., size classes) for which an average composition is determined. For m sections, a set of 2m(m+2) sectional coefficients must be calculated before integrating in time. These coefficients are determined from the basic coagulation, condensation, and deposition coefficients. Since the sectional coefficients depend on the physical properties of the containment chamber (e.g., temperature, pressure, chamber volume, and deposition surface area), they will generally need to be recalculated for a particular application. However, for a given containment chamber, the sectional coefficients will probably vary only with temperature and pressure. Consequently, the code has been developed so that sectional coefficients are stored at a user-specified upper and lower bound for both termperature and pressure, and linear interpolation is used to determine the appropriate sectional coefficients for a given temperature and pressure. A Runge-Kutta-Fehlberg method is used to integrate in time.
