A random packing algorithm is implemented to construct numerically models o
f heterogeneous propellants. These models consist of random distributions o
f spheres in a periodic cube. The spheres can hal-e various sizes, randomly
assigned if desired. Packing fractions are calculated for a bimodal model
(spheres of two different sizes) and compare vi ell with old experimental d
ata for steel shot. Slices through the cube define surfaces that represent
propellant surfaces, and the stoichiometry of these surfaces is a stochasti
c variable. Several realizations of the propellant cube are constructed and
variations in surface stoichiometry are examined using histograms and Four
ier series. A sample dame calculation is presented and solved using data de
fined by a single surface realization.