This study is a collaborative effort among lactation specialists and fluid dynamic engineers. The paper presents clinical results for suckling pressure pattern in lactating human breast as well as a 3D computational fluid dynamics (CFD) modeling of milk flow using these clinical inputs. The investigation starts with a careful, statistically representative measurement of suckling vacuum pressure, milk flow rate, and milk intake in a group of infants. The results from clinical data show that suckling action does not occur with constant suckling rate but changes in a rhythmic manner for infants. These pressure profiles are then used as the boundary condition for the CFD study using commercial ansys fluent software. For the geometric model of the ductal system of the human breast, this work takes advantage of a recent advance in the development of a validated phantom that has been produced as a ground truth for the imaging applications for the breast. The geometric model is introduced into CFD simulations with the aforementioned boundary conditions. The results for milk intake from the CFD simulation and clinical data were compared and cross validated. Also, the variation of milk intake versus suckling pressure are presented and analyzed. Both the clinical and CFD simulation show that the maximum milk flow rate is not related to the largest vacuum pressure or longest feeding duration indicating other factors influence the milk intake by infants.