  Many processes can correspond to reactive impregnation in porous solids.
These processes are usually numerically computed by classical methods like
finite element method, finite volume method, etc. The disadvantage of these
methods remains in the computational time. The convergence and accuracy require
a small step-time and a small mesh size, which is expensive in computational
time and can induce a spurious oscillation. In order to avoid this problem, we
propose a Self-organized Gradient Percolation algorithm. This method permits to
reduce the CPU time drastically.
