A pressure drop is specified between the left and right boundaries with no flow at top and bottom. Embedded elliptical regions have permeabilities varying between 0.01 and 50 against background permeability 1. The stream function and velocities are computed based on Darcy's law, using the "Immersed Interface Method" for solving elliptic equations with discontinuous coefficients on uniform Cartesian grids.
Transport of tracer moving in from left boundary (injected for 0 < t < 3.0) is computed using AMRCLAW. Adaptive mesh refinement is partially visible. The regions near the front where no mesh is shown have refinement by an additional factor 2 over the adjacent regions, for an effective 160 x 160 grid.
Results at 4 different times are shown in
