A ‘Variable Blend and Rates’ model is developed for SAG mills which extends upon the JKSimMet® Variable Rates Model. Similarly, a ‘Variable Blend Residence Constrained Perfect Mixing’ Ball Mill model is also introduced, also extending upon the existing JKSimMet® model. These models are used in conjunction with JKSimMet® while employing external technical programs. The technical model is then used to inform an economic model described elsewhere in Hanhiniemi and Powell (2022), which calculate operating and unit costs based on operating conditions, design, and ore characterisation. The SAG and ball mill models include constraints beyond those in the JKSimMet® models. The SAG model constraints include Size Specific Energy (SSE) verification, charge weight limit, and assessment of the impact of media to ore composition in the charge using the SAG Charge SSE Model developed in this paper. The ball mill model includes residence time, volumetric capacity, and SSE verification constraints. These models are applied to Cortez Gold Mine in the processing of two very different ore types and their blends; ‘Horse Canyon’ which is softer in coarser sizes and competent is fine sizes, and ‘Wenban’ which is competent at coarser sizes but softer in fine sizes than Horse Canyon. Horse Canyon also has a finer run-of-mine size distribution then Wenban. The application of the model at this site included undertaking five surveys at varying blends. The behaviours of the blends and individual components is then studied in the SABC circuit. A non_x0002_linearity is observed with throughputs increasing with blending. The newly developed modelling approach captures these responses. Controlled blending based on ore competence-by-size and size distribution was implemented by the site, which then realised significant improvements in throughput without any capital cost (Ntiamoah et al, 2020). The modelling developed in this paper is then compared to existing practice, and the new method was found to provide improved fidelity and utility. A three-dimensional surface is then developed using the model describing the SAG mill throughput capacity at varying conditions (speed-filling) versus varying blend. The change in throughput as a function of only the changing charge composition is presented. The result is nonlinear, and explanations are provided to describe this behaviour.