The current study was aimed at developing a package of “model + algorithm” for the design of resin in pulp (RIP) and carbon in pulp (CIP) processes of gold, uranium, and base metals. For this purpose, we first formulated a double‐resistance model for irreversible adsorption (accompanied by chemical reaction) in CSTRs, and modified the McKay et al. semi‐analytical model for reversible uptake in a similar system. We then devised two algorithms for the design and optimization of reversible and irreversible RIP and CIP cascades. The developed algorithms were applied on the extraction of copper, uranium, and gold. The packages are able to specify the optimum number of stages, reactor volume (V), resin flow rate (m.), and resin hold up (φs). This shows an evident advantage over the McCabe‐Thiele method, whose only result is the number of stages. The new methods have only two adjustable parameters of diffusivity (Deff) and liquid film coefficient (k), whose determination can be carried out through simple batch experiments. As an additional work, this study also presents a method for mathematical troubleshooting for the Mc Kay et al. model.