Solving a steady flow in a subdomain of a mesh using FiPy

Question

I\'m fairly new to FiPy and I\'m currently facing an issue of which I am sure can be solved easily: I want to solve a 3D steady flow of the form:

    eq = ( Diffu