This thesis is concerned with the simulation of particle physics processes involving strong interactions in modern event generators. New algorithms to reinstate colour in colour- ordered amplitudes through colour dressing are presented and their analytical and numerical properties are discussed in detail. The colour-dressed Berends-Giele recursive relations are extended to the full Standard Model Lagrangian and implemented into the numerical program Comix for large multiplicity matrix element computation. New algorithms for phase space integration are proposed, whereof one is capable to effectively couple colour and momentum sampling. Comparisons to other high-multiplicity generators are shown. QCD parton evolution and the CKKW algorithm to...