Delivery of drug using nanocarriers tethered with vasculature-targeting epitopes aims to maximize the therapeutic efficacy of the drug while minimizing the drug side effects. Circadian rhythm which is governed by the central nervous system has implications for targeted drug delivery due to sleep-wake cycle changes in blood flow dynamics. This paper presents an advanced fluid dynamics modeling method that is based on viscous incompressible shear-rate fluid (blood) coupled with an advection-diffusion equation to simulate the formation of drug concentration gradients in the blood stream and buildup of concentration at the targeted site. The method is equipped with an experimentally calibrated nanoparticle-endothelial cell adhesion model that employs Robin boundary conditions to describe nanoparticle retention based on probability of adhesion, a friction model accounting for surface roughness of endothelial cell layer, and a dispersion model based on Taylor-Aris expression for effective diffusion in the boundary layer. The computational model is first experimentally validated and then tested on engineered bifurcating arterial systems where impedance boundary conditions are applied at the outflow to account for the downstream resistance at each outlet.It is then applied to a virtual geometric model of anin vivoarterial tree developed through MRI-based image processing techniques. These simulations highlight the potential of the computational model for drug transport, adhesion, and retention at multiple sites in virtualin vivomodels. The model provides a virtual platform for exploring circadian rhythm modulated blood flow for targeted drug delivery while minimizing thein vivoexperimentation.