Diffuse interface approach to oxygen transport and metabolism under blood flow dynamics in microcirculations

The relationship between the spatiotemporal distribution of oxygen transport and blood flow dynamics, accounting for the motion and deformation of individual red blood cells (RBCs), is of fundamental importance for understanding microcirculation systems. Three-dimensional (3D) modeling is indispensable for addressing complex oxygen transport and cellular behaviors in capillary networks; however, the computational approach is formidable for enforcing interface (or jump) conditions on largely moving and deforming interfaces. In this paper, we propose a diffuse interface approach for the oxygen transport using a mixture formulation. We formulate oxygen transport using an advection-diffusion-reaction equation and rewrite all governing equations in mixture forms using phase indicator functions, where all the interface conditions are included in the governing equations. This innovation avoids the complexity associated with discontinuities for largely moving interfaces in highly dense RBC conditions. We model cellular flow as a fluid-membrane interaction problem using the immersed boundary method (IBM). The method allows the seamless calculation of coupling problems for cellular flows and oxygen transports in the cytoplasm (internal fluid) of the RBC, plasma (external fluid), and tissue regions using a fixed Cartesian coordinate mesh. The proposed method accurately captures the analytical solution for spherically symmetric diffusion, and successfully demonstrates oxygen transport in both straight capillaries and their networks. These rigorous analyses suggest that RBCs can autonomously regulate the oxygen supply to tissues in response to the local tissue oxygenation level, resulting in the establishment of homogeneous tissue oxygenation.