Computational exploration of injection strategies for improving medicine distribution in the liver

The liver, a vital metabolic organ, is always susceptible to various diseases that ultimately lead to fibrosis, cirrhosis, acute liver failure, chronic liver failure, and even cancer. Optimal and specific medicine delivery in various diseases, hepatectomy, shunt placement, and other surgical interventions to reduce liver damage, transplantation, optimal preservation, and revival of the donated organ all rely on a complete understanding of perfusion and mass transfer in the liver. This study aims to simulate the computational fluid dynamics of perfusion and the temporal-spatial distribution of a medicine in a healthy liver to evaluate the hemodynamic characteristics of flow and medicine transport with the purpose of more effective liver treatment.

Patient-specific geometries of parenchyma and hepatic artery, portal vein, and hepatic vein vessels of a healthy liver were segmented and reconstructed from the abdominal computed tomography scan images. Mesh was generated for the comprehensive combined model using unstructured tetrahedral elements. Transient pressure values were applied as boundary conditions at the portal vein and hepatic artery inlets, and pressure outlet boundary condition was assumed at the hepatic vein outlet. Medicine injection was done through the portal vein. The liver parenchyma was assumed to be a porous medium. Finally, computational fluid dynamics (CFD) simulation was performed to investigate blood perfusion, medicine distribution, and saturation time.

The velocity parameter values calculated for the hepatic artery, portal vein, and hepatic vein vessels were consistent with the physiological ranges. The mass flow rate was higher in the portal vein than in the hepatic artery, which was consistent with high perfusion through the portal vein. The portal pressure gradient was 8.53 mmHg. From a pharmacokinetic viewpoint, medicine distribution in porous tissue was a heterogeneous process. Medicine distribution was higher at end-diastolic pressure than at peak-systolic pressure which showed the influence of hepatic artery flow. The tissue was saturated faster at first 40 s and with decreasing porosity, saturation time decreased, and distribution improved.

The right lobe included a higher number of vascular terminals due to its larger volume, and the flow rate was higher in this lobe compared to the left lobe. This showed the significant effect of the right lobe on the overall function of the body. Recirculation flow zones along hepatic artery and portal vein branches emphasized the sensitivity of downstream vessels. Rotational flow and potential vortex formation at the hepatic vein outlet may indicate a risk of plaque and clot formation in this region. The heterogeneous distribution of medicine indicated the importance of injection time in treating liver diseases. The percentage of tissue porosity affected the saturation time, so adjusting the medicine dose and injection time could be challenging in treatments.