We extend estimation of range flow to handle brightness changes in image data caused by inhomogeneous illumination. Standard range flow computes 3D velocity fields using both range and intensity image sequences. Towards this end, range flow estimation combines a depth change model with a brightness constancy model. However, local brightness is generally not preserved when object surfaces rotate relative to the camera or the light sources, or when surfaces move in inhomogeneous illumination. We describe and investigate different approaches to handle such brightness changes. A straightforward approach is to prefilter the intensity data such that brightness changes are suppressed, for instance by a highpass or a homomorphic filter. Such prefiltering may, though, reduce the signal to noise ratio. An alternative novel approach is to replace the brightness constancy model by (1) a gradient constancy model, or (2) by a combination of gradient and brightness constancy constraints used earlier successfully for optical flow, or (3) by a physics-based brightness change model. In performance tests, the standard version and the novel versions of range flow estimation are investigated using prefiltered or non-prefiltered synthetic data with available ground truth. Furthermore, the influences of additive Gaussian noise and simulated shot noise are investigated. We finally compare all range flow estimators on real data.