The prediction of soot emissions is a requirement for the development of gas turbine combustors. Apart from global quantities related to soot mass, future regulations also call for the control of particle number. Therefore, theoretical models for soot from combustion devices must include various nucleation, growth, and oxidation mechanisms and aerosol physics in order to predict the soot particle number distribution. This paper introduces an approach based on Incompletely Stirred Reactor Network (ISRN) modeling that simplifies calculations and facilitates parametric analyses with very complex soot models. The method is based on the Conditional Moment Closure (CMC) combustion model and Incompletely Stirred Reactor (ISR) theory, which is here extended to a reactor network formulation. An ISR is a volume that is inhomogeneous in terms of mixture fraction but is characterized by homogeneous conditional averages, such as temperature conditioned on the mixture fraction having a particular value. A network of ISRs can then be deployed to separately capture soot production and oxidation regions exhibiting different degrees of micro-mixing rates and residence times, as typically observed in practical combustors. The ISRN approach is demonstrated in this paper for a model aero-engine combustor for which detailed CFD and experimental data exist, and it is found that reasonable accuracy is produced at a significantly reduced computational cost.