A model is formulated to describe dissolution of naphthalene from an insoluble nonaqueous phase liquid (NAPL) and its subsequent biodegradation in the aqueous phase in completely mixed batch reactors. The physicochemical processes of equilibrium partitioning and mass transfer of naphthalene between the NAPL and aqueous phases were incorporated into the model. Biodegradation kinetics were described by Monod＇s microbial growth kinetic model, modified to account for the inhibitory effects of 1,2-naphthoquinone formed during naphthalene degradation under certain conditions. System parameters and biokinetic coefficients pertinent to the NAPL-water systems were determined either by direct measurement or from nonlinear regression of the naphthalene mineralization profiles obtained from batch reactor tests with two-component NAPLs comprised of naphthalene and heptamethylnonane. The NAPLs contained substantial mass of naphthalene, and naphthalene biodegradation kinetics were evaluated over the time required for near complete depletion of naphthalene from the NAPL. Model predictions of naphthalene mineralization time profiles compared favorably to the general trends observed in the data obtained from laboratory experiments with the two-component NAPL, as well as with two coal tars obtained from the subsurface at contaminated sites and composed of many different PAHs (polycyclic aromatic hydrocarbon compounds). The effects of varying the NAPL mass and the naphthalene mole fractions in the NAPL are discussed. It was observed that the time to achieve a given percent removal of naphthalene does not change significantly with the initial mass of naphthalene in a fixed volume of the NAPL. Significant changes in the mineralization profiles are observed when the volume (and mass) of NAPL in the system is changed.