  The multiple-network poroelasticity (MPET) equations describe deformation and
pressures in an elastic medium permeated by interacting fluid networks. In this
paper, we (i) place these equations in the theoretical context of coupled
elliptic-parabolic problems, (ii) use this context to derive residual-based a
posteriori error estimates and indicators for fully discrete MPET solutions and
(iii) evaluate the performance of these error estimators in adaptive algorithms
for a set of test cases: ranging from synthetic scenarios to physiologically
realistic simulations of brain mechanics.
