In the present study, a dynamic mathematical model for the growth of granulocyte progenitor cells in the hematopoietic process is developed based on the principles of diffusion and chemical reaction. This model simulates granulocyte progenitor cell growth and oxygen consumption in a three-dimensional (3-D) perfusion bio-reactor. Material balances on cells are coupled to the nutrient balances in 3-D matrices to determine the effects of transport limitations on cell growth. The method of volume averaging is used to formulate the material balances for the cells and the nutrients in the porous matrix containing the cells. All model parameters are obtained from the literature. The maximum cell volume fraction reached when oxygen is depleted in the cell layer at 15 days and is nearly 0.63, corresponding to a cell density of 2.25 x 10(8) cells/mL. The substrate inhibition kinetics for cell growth lead to complex effects with respect to the roles of oxygen concentration and supply by convection and diffusion on cell growth. Variation in the height of the liquid layer above the cell matrix where nutrient supply is introduced affected the relative and absolute amounts of oxygen supply by hydrodynamic flow and by diffusion across a gas permeable FEP membrane. Mass transfer restrictions of the FEP membrane are considerable, and the supply of oxygen by convection is essential to achieve higher levels of cell growth. A maximum growth rate occurs at a specific flow rate. For flow rates higher than this optimal, the high oxygen concentration led to growth inhibition and for lower flow rates growth limitations occur due to insufficient oxygen supply. Because of the non-linear effects of the autocatalytic substrate inhibition growth kinetics coupled to the convective transport, the rate of growth at this optimal flow rate is higher than that in a corresponding well-mixed reactor where oxygen concentration is set at the maximum indicated by the inhibitory kinetics.(C) 2005 Wiley Periodicals, Inc.