Numerical solution of the advectionereactionediffusion equation at different scales

Solving the transport equation for bimolecular reactive processes in porous media involves several difficulties. The mathematical characteristics of the equation depend on the governing process, for example, when time scales for advection tA, reaction tR and diffusion tD have different orders of magnitude. On the other hand, this equation is based on a continuum model, disregarding inhomogeneities that happen at poral level just where reactions take place. To deal with these problems a different way of modeling the advectionereactionediffusion process is proposed. Based on the Damkohler number an algorithm has been developed to solve the problem for both slow and fast reactions. Spatial segregation of reactants can be incorporated improving the results.