Pb in drinking water is one of the main sources of Pb in human bodies. Besides excess lead incidents caused by lead service lines, Pb contamination in non-lead pipe systems (e.g. copper and PVC) is also on the rise. Brass fixtures and lead-solder connections are the primary sources of lead in non-lead systems, which cause intermittent peaks of Pb at the consumer tap. The concentrations of lead in tap water depend on pipeline geometry, sampling methods, and the characteristic of the pipe flow. Using a 3-dimensional computational fluid dynamics model, this study analyzes the Pb concentration variations at taps of copper water supply systems. The turbulent diffusion and shear flow dispersion are fully simulated in the model to provide the most accurate prediction. Water parcels containing lead (clouds) are formed adjacent to lead sources during stagnation and are then dispersed with the flow of water through the pipe when the tap is opened. The geometry of the pipeline has a significant impact on the monitored Pb levels. The complex flow condition in elbow areas leads to a more intense mixing of lead ions. Therefore, the Pb levels at consumer taps in complex plumbing systems in high-rise buildings are higher than in straight pipelines. When the sampling flow rate is large, the peak height of Pb is slightly higher due to higher turbulent intensity. Pb concentrations in sequential samples are predicted with sample volumes of 50, 250, 500, and 1000 mL. Lead levels may be diluted significantly when samples are taken in large volumes (e.g. 500 mL or 1 L), resulting in an underestimation of Pb levels at taps. A sample volume of 250–500 mL is recommended in sequential sampling programs on non-lead service lines in order to identify Pb sources.