A new method for the numerical simulation of the radiation chemistry of aqueous polymer solutions is introduced. The method makes use of a deterministic approach combining the conventional homogeneous radiation chemistry of water with the chemistry of polymer radicals and other macromolecular species. The method is applied on single-pulse irradiations of aqueous polymer solutions. The speciation of macromolecular species accounts for the variations in the number of alkyl radicals per chain, molecular weight, and number of internal loops (as a consequence of an intramolecular radical-radical combination). In the simulations, the initial polymer molecular weight, polymer concentration, and dose per pulse (function of pulse length and dose rate during the pulse) were systematically varied. In total, 54 different conditions were simulated. The results are well in line with the available experimental data for similar systems. At a low polymer concentration and a high dose per pulse, the kinetics of radical decay is quite complex for the competition between intra- and intermolecular radical-radical reactions, whereas at a low dose per pulse the kinetics is purely second-order. The simulations demonstrate the limitations of the polymer in scavenging all the radicals generated by water radiolysis when irradiated at a low polymer concentration and a high dose per pulse. They also show that the radical decay of lower-molecular-weight chains is faster and to a larger extent dominated by intermolecular radical-radical reactions, thus explaining the mechanism behind the experimentally observed narrowing of molecular weight distributions.