Arora, GouravGouravAroraHussain, SaddamSaddamHussainWalker, GavinGavinWalkerCattani, CarloCarloCattaniHeinrich, StefanStefanHeinrichSingh, MehakpreetMehakpreetSingh2026-02-272026-02-272026-01-29Chemical Engineering Science 325: 123400 (2026)https://hdl.handle.net/11420/61790Particulate processes such as bubble breakage, granule preparation, and ore extraction are often modeled using population balance equations. While most studies rely on linear fragmentation models, these fail to capture realistic dynamics where breakage occurs due to particle collisions, leading to nonlinear collisional breakage equations. The mathematical complexity of nonlinear kernels and integrals has limited progress in this direction, with very few analytical solutions available. In this work, we develop an accurate approach based on a modified variational iteration method for solving multidimensional collisional breakage equations. The method allows flexible operator selection and efficient determination of the Lagrange multiplier. A rigorous convergence and error analysis is carried out using the fixed-point theorem. The approach is further extended to two- and three-dimensional models, including collisional fragmentation cases. Numerical experiments confirm the accuracy and efficiency of the method compared with existing approaches.en0009-25092026Elsevier34A3435Q7045K0545L10Convergence analysisFinite volume schemeLagrange multiplierMathematics Subject Classification (AMC) MSC 34A12Modified variational iteration methodNonlinear equationParticle dynamicsTechnology::600: TechnologyTechnology::660: Chemistry; Chemical EngineeringJournal Article10.1016/j.ces.2026.123400Journal Article