Abstract. This study established a multiscale computational model for intraperitoneal (IP) chemotherapy, to depict the time-dependent and spatial-dependent drug concentrations in peritoneal tumors as functions of drug properties (size, binding, diffusivity, permeability), transport mechanisms (diffusion, convection), spatial-dependent tumor heterogeneities (vessel density, cell density, pressure gradient), and physiological properties (peritoneal pressure, peritoneal fluid volume). Equations linked drug transport and clearance on three scales (tumor, IP cavity, whole organism). Paclitaxel was the test compound. The required model parameters (tumor diffusivity, tumor hydraulic conductivity, vessel permeability and surface area, microvascular hydrostatic pressure, drug association with cells) were obtained from literature reports, calculation, and/or experimental measurements. Drug concentration-time profiles in peritoneal fluid and plasma were the boundary conditions for tumor domain and blood vessels, respectively. The finite element method was used to numerically solve the nonlinear partial differential equations for fluid and solute transport. The resulting multiscale model accounted for intratumoral spatial heterogeneity, depicted diffusive and convective drug transport in tumor interstitium and across blood vessels, and provided drug flux and concentration as a function of time and spatial position in the tumor. Comparison of model-predicted tumor spatiokinetics with experimental results (autoradiographic data of 3H-paclitaxel in IP ovarian tumors in mice, 6 h posttreatment) showed good agreement (1% deviation for area under curve and 23% deviations for individual data points, which were several-fold lower compared to the experimental intertumor variations). The computational multiscale model provides a tool to quantify the effects of drug-, tumor-, and host-dependent variables on the concentrations and residence time of IP therapeutics in tumors.