<p>Subcutaneous administration has become a common approach for drug delivery of biotherapeutics, such as monoclonal antibodies, which is achieved mainly by absorption through the lymphatic system. This dissertation focuses on the computational modeling of the fluid flow and solute transport in the skin tissue and the quantitative investigation of lymphatic uptake. First, the various mechanisms governing drug transport and lymphatic uptake of biotherapeutics through subcutaneous injection are investigated quantitatively through high-fidelity numerical simulations, including lymphatic drainage, blood perfusion, binding, and metabolism. The tissue is modeled as a homogeneous porous medium using both a single-layered domain and a multi-layered domain, which includes the epidermis, dermis, hypodermis (subcutaneous tissue), and muscle layers. A systematic parameter study is conducted to understand the roles of different properties of the tissue in terms of permeability, porosity, and vascular permeability. The role of binding and metabolism on drug absorption is studied by varying the binding parameters for different macromolecules after coupling the transport equation with a pharmacokinetic equation. The interstitial pressure plays an essential role in regulating the absorption of unbound drug proteins during the injection, while the binding and metabolism of drug molecules reduce the total free drugs. </p> <p>  </p> <p>The lymphatic vessel network is essential to achieve the functions of the lymphatic system. Thus, the drug transport and lymphatic uptake through a three-dimensional hybrid discrete-continuum vessel network in the skin tissue are investigated through high-fidelity numerical simulations. The explicit heterogeneous vessel network is embedded into the continuum model to investigate the role of explicit heterogeneous vessel network in drug transport and absorption. The solute transport across the vessel wall is investigated under various transport conditions. The diffusion of the drug solutes through the explicit vessel wall affects the drug absorption after the injection, while the convection under large interstitial pressure dominates the drug absorption during the injection. The effect of diffusion cannot be captured by the previously developed continuum model. Furthermore, the effects of injection volume and depth on the lymphatic uptake are investigated in a multi-layered domain. The injection volume significantly affects lymphatic uptake through the heterogeneous vessel network, while the injection depth has little influence. At last, the binding and metabolism of drug molecules are studied to bridge the simulation to the experimentally measured drug clearance. </p> <p><br></p> <p>Convective transport of drug solutes in biological tissues is regulated by the interstitial fluid pressure, which plays a crucial role in drug absorption into the lymphatic system through the subcutaneous (SC) injection.  An approximate continuum poroelasticity model is developed to simulate the pressure evolution in the soft porous tissue during an SC injection. This poroelastic model mimics the deformation of the tissue by introducing the time variation of the interstitial fluid pressure. The advantage of this method lies in its computational time efficiency and simplicity, and it can accurately model the relaxation of pressure. The interstitial fluid pressure obtained using the proposed model is validated against both the analytical and the numerical solution of the poroelastic tissue model. The decreasing elasticity elongates the relaxation time of pressure, and the sensitivity of pressure relaxation to elasticity decreases with the hydraulic permeability, while the increasing porosity and permeability due to deformation alleviate the high pressure. </p> <p><br></p> <p>At last, an improved Kedem-Katchalsky model is developed to study solute transport across the lymphatic vessel network, including convection and diffusion in the multi-layered poroelastic tissue with a hybrid discrete-continuum vessel network embedded inside. The effect of different drug solutes with different Stokes radii and different structures of the lymphatic vessel network, such as fractal trees and Voronoi structure, on the lymphatic uptake is investigated. The drug solute with a small size has a larger partition coefficient and diffusivity across the openings of the lymphatic vessel wall, which favors drug absorption. The Voronoi structure is found to be more efficient in lymphatic uptake. </p> <p><br></p> <p>The systematic and quantitative investigation of subcutaneous absorption based on high-fidelity numerical simulations can provide guidance on the optimization of drug delivery systems and is valuable for the translation of bioavailability from the pre-clinical species to humans. We provide a novel approach to studying the diffusion and convection of drug molecules into the lymphatic system by developing the hybrid discrete-continuum vessel network. The study of the solute transport across the discrete lymphatic vessel walls further improves our understanding of lymphatic uptake. The novel and time-efficient computational model for solute transport across the lymphatic vasculature connects the microscopic properties of the lymphatic vessel membrane to macroscopic drug absorption. The comprehensive hybrid vessel network model developed here can be further used to improve our understanding of the diseases caused by the disturbed lymphatic system, such as lymphedema, and provide insights into the treatment of diseases caused by the malfunction of lymphatics.</p>