This commit cleans up the incompressible solver and adds a DarcyFlux solver in model module. Some documentation for both classes has been added. The convection code is only linear at the moment.