add preliminary Darcy solver functions