Multiscale modeling of lymphatic drainage from tissues using homogenization theory

Lymphatic capillary drainage of interstitial fluid under both steady-state and inflammatory conditions is important for tissue fluid balance, cancer metastasis, and immunity. Lymphatic drainage function is critically coupled to the fluid mechanical properties of the interstitium, yet this coupling is poorly understood. Here we sought to effectively model the lymphatic-interstitial fluid coupling and ask why the lymphatic capillary network often appears with roughly a hexagonal architecture. We use homogenization method, which allows tissue-scale lymph flow to be integrated with the microstructural details of the lymphatic capillaries, thus gaining insight into the functionality of lymphatic anatomy. We first describe flow in lymphatic capillaries using the Navier–Stokes equations and flow through the interstitium using Darcy's law. We then use multiscale homogenization to derive macroscale equations describing lymphatic drainage, with the mouse tail skin as a basis. We find that the limiting resistance for fluid drainage is that from the interstitium into the capillaries rather than within the capillaries. We also find that between hexagonal, square, and parallel tube configurations of lymphatic capillary networks, the hexagonal structure is the most efficient architecture for coupled interstitial and capillary fluid transport; that is, it clears the most interstitial fluid for a given network density and baseline interstitial fluid pressure. Thus, using homogenization theory, one can assess how vessel microstructure influences the macroscale fluid drainage by the lymphatics and demonstrate why the hexagonal network of dermal lymphatic capillaries is optimal for interstitial tissue fluid clearance.