A method of forming a wormhole in a porous medium, comprising: running a stimulation simulator comprising: gridding a treatment region of the porous medium into a plurality of cells comprising a first portion designated as matrix cells and a second portion designated as wormhole cells; populating the simulator with static properties of the porous medium and reaction kinetic properties for reaction of the porous medium with a reactant in a treatment fluid; modeling the matrix cells wherein a medium of the matrix cells comprises matrix material behaving as a single permeability, single porosity system; modeling the wormhole cells in a wormhole initiation stage wherein a medium of the respective wormhole initiation stage cells has a solid saturation above a respective critical solid saturation and is comprised of the matrix material behaving as a single permeability, single porosity system; modeling at least a portion of the wormhole cells in a wormhole growth stage wherein the respective wormhole cells have a solid saturation equal to or less than the respective critical sold saturation, and wherein the wormhole growth stage cells comprise two different interconnected media comprised respectively of the matrix material and a wormhole material having a fluid mobility as a function of solid saturation; and obtaining optimized treatment fluid injection parameters; and injecting the treatment fluid into the treatment region of the porous medium according to the optimized treatment fluid injection parameters to form the wormhole. running a stimulation simulator comprising: gridding a treatment region of the porous medium into a plurality of cells comprising a first portion designated as matrix cells and a second portion designated as wormhole cells; populating the simulator with static properties of the porous medium and reaction kinetic properties for reaction of the porous medium with a reactant in a treatment fluid; modeling the matrix cells wherein a medium of the matrix cells comprises matrix material behaving as a single permeability, single porosity system; modeling the wormhole cells in a wormhole initiation stage wherein a medium of the respective wormhole initiation stage cells has a solid saturation above a respective critical solid saturation and is comprised of the matrix material behaving as a single permeability, single porosity system; modeling at least a portion of the wormhole cells in a wormhole growth stage wherein the respective wormhole cells have a solid saturation equal to or less than the respective critical sold saturation, and wherein the wormhole growth stage cells comprise two different interconnected media comprised respectively of the matrix material and a wormhole material having a fluid mobility as a function of solid saturation; and obtaining optimized treatment fluid injection parameters; and gridding a treatment region of the porous medium into a plurality of cells comprising a first portion designated as matrix cells and a second portion designated as wormhole cells; populating the simulator with static properties of the porous medium and reaction kinetic properties for reaction of the porous medium with a reactant in a treatment fluid; modeling the matrix cells wherein a medium of the matrix cells comprises matrix material behaving as a single permeability, single porosity system; modeling the wormhole cells in a wormhole initiation stage wherein a medium of the respective wormhole initiation stage cells has a solid saturation above a respective critical solid saturation and is comprised of the matrix material behaving as a single permeability, single porosity system; modeling at least a portion of the wormhole cells in a wormhole growth stage wherein the respective wormhole cells have a solid saturation equal to or less than the respective critical sold saturation, and wherein the wormhole growth stage cells comprise two different interconnected media comprised respectively of the matrix material and a wormhole material having a fluid mobility as a function of solid saturation; and obtaining optimized treatment fluid injection parameters; and injecting the treatment fluid into the treatment region of the porous medium according to the optimized treatment fluid injection parameters to form the wormhole.