A bond graph method is described for treating nonlinear diffusion-chemical reaction processes in homogeneous or heterogeneous media. A program has been developed to study steady-state and dynamic response for a system with an arbitrary number of reactants and reactions. An experimental system selected for study is the facilitated diffusion of oxygen through a hemoglobin solution.