Dissolution of lead, and its subsequent transport from source to the point of consumption in a household drinking water system is modeled. Radial diffusion and convective mass transfer are used to determine the increase in lead concentration in water under stagnant and flow conditions, respectively. Subsequent transport to the point of consumption is modeled using axial convection. A computer code has been developed to predict the instantaneous and mean lead concentrations at the faucet for a given water consumption pattern, and locations of multiple lead sources in the pipeline. The model has been used to analyze the effect of protocol design on measured variability in lead levels.