volatile organic compounds (vocs) are one of the major concerns for indoor air quality. a new method for treating vocs in the air is through an advance oxidation process (aop) known as photocatalytic oxidation (pco). this method uses ultraviolet (uv) radiant power to activate the surface of a catalyst (e.g. ti02). from the past experimental work that has been done, it can be concluded that pco reactions follow the langmuir-hinshelwood (l-h) kinetic model. however, the complexity of the l-h kinetic model is difficult to simulate using existing computation fluid dynamic (cfd) software. in this thesis, a new method for modeling the l-h surface reaction kinetics is proposed. the focus of this work is on the development of a novel approach to model the complex surface reaction rate expressions in order to define pco reaction rates on the photocatalyst surfaces. a new approach is developed to adapt the overall experimental reaction rates, which are in terms of the total system volume. this adaptation will help in deriving the actual rate of reaction happening on the catalyst surface in terms of catalyst surface area. two cases were studied in order to demonstrate how this new approach can be used to accurately model the complex reaction kinetics of pco systems. in each case, an integrated cfd model was developed to accurately predict the rate ofvoc decomposition based on the work conducted by shiraishi et al., (2005b) and brosillon et al., (2008). in the first case, the experimental kinetic model for formaldehyde decomposition was adapted in order to describe a surface reaction based on formaldehyde concentration on the catalyst surface using three different approaches. it was determined that a two part polynomial rate expression was the most accurate one, as it was able to account for the higher initial rate of reaction. however, the exponential model did give reasonable results as well. in the second case, the reaction rate model was able to predict the rate of decomposition for butyric acid in the air for a variety of initial concentrations and uv irradiance levels at the catalyst surface. the developed cfd model results also discredited the assumptions made in a number of published papers that uv irradiance levels are uniform across a catalyst surface. finally, a simple case study was developed in order to demonstrate how the novel approach to reaction modeling could be used to predict pco system performance treating air in a close system. in this case, the pco system was capable of treating air contaminated with butyric acid, as well as quickly reducing the concentration below the odor threshold.