Diffusion in pure gels and gels with immobilized cells was analyzed. A model of diffusion assuming a homogeneous cell distribution in gel was improved by introducing a tortuosity value. By theoretical analysis and numerical modeling it was shown that the tortuosity of a gel with immobilized cells is the product of two factors: (1) tortuosity generated by the cells, T c , and (2) tortuosity of the gel matrix, T g , both variables being a function of cell volume fraction, φ c . Total tortuosity is thus T Σ ) T c T g . On the basis of this approach, it was possible to analyze diffusivity data for gels with immobilized cells. It was shown that, in these systems, the diffusivity η ) D e /D 0 is a complex function of (1) diffusivity in the gel, η g , and (2) diffusivity in immobilized cells, η c . The developed model allowed for the description of the dependence of D e /D 0 on φ c . Comparison with numerous published experimental data showed a good fit. Observed deviations might be explained by nonhomogeneous cell distributions inside the gel matrix.