Dissolution of pore-scale soluble substances occurs in applications from environmental hydrogeology to CO2 storage. Development of Darcy-scale models has been widely discussed. This paper proposes an upscaling algorithm to develop large-scale models taking into account Darcy-scale heterogeneities. The theory is based on a comparison between the characteristic length-scale of the Darcy-scale heterogeneities and the length-scale of the dissolution front controlled by a Darcy-scale Damkhler number. The obtained large-scale model is shown to confront favorably to Darcy-scale direct numerical simulations
