  In this paper, we propose a new finite element solution approach to the
multi-compartmental Darcy equations describing flow and interactions in a
porous medium with multiple fluid compartments. We introduce a new numerical
formulation and a block-diagonal preconditioner. The robustness with respect to
variations in material parameters is demonstrated by theoretical considerations
and numerical examples.
