A simple numerical model to simulate athermal avalanches is presented. The model is inspired by the "porous collapse" process where the compression of porous materials generates collapse cascades, leading to power law distributed avalanches. The energy (E), amplitude (A_{max}), and size (S) exponents are derived by computer simulation in two approximations. Time-dependent "jerk" spectra are calculated in a single avalanche model where each avalanche is simulated separately from other avalanches. The average avalanche profile is parabolic, the scaling between energy and amplitude follows E∼A_{max}^{2}, and the energy exponent is ε = 1.33. Adding a general noise term in a continuous event model generates infinite avalanche sequences which allow the evaluation of waiting time distributions and pattern formation. We find the validity of the Omori law and the same exponents as in the single avalanche model. We then add spatial correlations by stipulating the ratio G/N between growth processes G (linked to a previous event location) and nucleation processes N (with new, randomly chosen nucleation sites). We found, in good approximation, a power law correlation between the energy exponent ε and the Hausdorff dimension H_{D} of the resulting collapse pattern H_{D}-1∼ɛ^{-3}. The evolving patterns depend strongly on G/N with the distribution of collapse sites equally power law distributed. Its exponent ɛ_{topo} would be linked to the dynamical exponent ε if each collapse carried an energy equivalent to the size of the collapse. A complex correlation between ɛ,ɛ_{topo}, and H_{D} emerges, depending strongly on the relative occupancy of the collapse sites in the simulation box.