Research from this group has recently led to the development of a genipin-crosslinked gelatin gel that may be used as a controlled release matrix for bioactive compounds. This study presents a model that simulates the release of entrapped caffeine from the hydrogel and the ingress of water into these gels. Fick's second law Abbreviations: A, Dimensionless front-fixing radius; B, Dimensionless front-fixing height; C 0 , Caffeine concentration at initial conditions in the cylindrical matrix; C 1,eq , Equilibrium water concentration in the fully swollen cylindrical matrix; C 1 , Water concentration in the cylindrical matrix; C 2 , Caffeine concentration in the cylindrical matrix;C 1 , Dimensionless water concentration in the cylindrical matrix; C 2 , Dimensionless caffeine concentration in the cylindrical matrix; D 1 , Water diffusion coefficient in the cylindrical matrix; D 2 , Caffeine diffusion coefficient in the cylindrical matrix; D 1,eq , Diffusion coefficient of caffeine in the fully swollen cylindrical matrix; D 2,eq , Diffusion coefficient of water in the fully swollen cylindrical matrix; r , Radial coordinate in the cylindrical matrix; r 0 , Initial radius of the cylindrical matrix;r , Dimensionless radius of the cylindrical matrix; r N , New radius of the cylindrical matrix; r t , (New radius of the cylindrical matrix)/(initial radius of the cylindrical matrix); t, Time; z, Axial coordinate in the cylindrical matrix; z 0 , Initial half-height of the cylindrical matrix;z, Dimensionless half-height of the cylindrical matrix; z N , New half height of cylindrical matrix; z t , (New half-height of the cylindrical matrix)/(initial half-height of the cylindrical matrix); V p , Dry matrix volume; β 1 , Dimensionless constant for the diffusion of water; β 2 , Dimensionless constant for the diffusion of caffeine; ρ 1 , Density of water; ρ 2 , Density of caffeine; τ , Dimensionless time.of diffusion is used to describe the water penetration and bioactive release in the system. To validate the model, caffeine release experiments from the gel were carried out. The predicted bioactive release profiles are in very good agreement with experimental data at different gel compositions. The model may also be used for a wide range of bioactive molecules and hydrogels with different cylindrical dimensions.