The focus of this paper is the numerical solution of space-fractional heat conduction equations. Typical numerical treatments for fractional diffusion equations (FDEs) are constructed to function on structured grids and cannot be readily applied to problems in arbitrarily shaped domains. Presented is an unstructured control volume finite element method for the solution of FDEs in non-Cartesian domains. This method approximates the fractional derivative flux at given integration points with a weighted sum of gradients along orthogonal spines. Accuracy and utility of the approach are demonstrated by solution of one- and two-sided FDEs in rectangular and cylindrical domains.