Allow BodyForce to use variables as arguments
This MR adds the ability to use variables in BodyForce
function definitions, so that certain simulations are now enabled with the use of IMEX time stepping -- in particular, buoyancy terms in IncNavierStokesSolver
for convective flows and simple reaction-diffusion equations in ADRSolver
(through the use of a new equation system that this MR introduces).