Direct infusions into brain parenchyma of biological therapeutics for serious brain diseases have been, and are being, considered. However, individual brains, as well as distinct cytoarchitectural regions within brains, vary in their response to fluid flow and pressure. Further, the tissue responds dynamically to these stimuli, requiring a nonlinear treatment of equations that would describe fluid flow and drug transport in brain. We here report in detail on an individual–specific model and a comparison of its prediction with simulations for living porcine brains. Two critical features we introduced into our model — absent from previous ones, but requirements for any useful simulation — are the infusion-induced interstitial expansion and the backflow. These are significant determinants of the flow. Another feature of our treatment is the use of cross–property relations to obtain individual–specific parameters that are coefficients in the equations. The quantitative results are at least encouraging, showing a high fraction of overlap between the computed and measured volumes of distribution of a tracer molecule, and are potentially clinically useful. Several improvements are called for; principally a treatment of the interstitial expansion more fundamentally based on poroelasticity, and a better delineation of the diffusion tensor of a particle confined to the interstitial spaces.