Although a full understanding of the hepatic circulation is one of the keys to successfully perform liver surgery and to elucidate liver pathology, relatively little is known about the functional organization of the liver vasculature. Therefore, we materialized and visualized the human hepatic vasculature at different scales, and performed a morphological analysis by combining vascular corrosion casting with novel micro‐computer tomography (CT) and image analysis techniques. A human liver vascular corrosion cast was obtained by simultaneous resin injection in the hepatic artery (HA) and portal vein (PV). A high resolution (110 μm) micro‐CT scan of the total cast allowed gathering detailed macrovascular data. Subsequently, a mesocirculation sample (starting at generation 5; 88 × 68 × 80 mm³) and a microcirculation sample (terminal vessels including sinusoids; 2.0 × 1.5 × 1.7 mm³) were dissected and imaged at a 71‐μm and 2.6‐μm resolution, respectively. Segmentations and 3D reconstructions allowed quantifying the macro‐ and mesoscale branching topology, and geometrical features of HA, PV and hepatic venous trees up to 13 generations (radii ranging from 13.2 mm to 80 μm; lengths from 74.4 mm to 0.74 mm), as well as microvascular characteristics (mean sinusoidal radius of 6.63 μm). Combining corrosion casting and micro‐CT imaging allows quantifying the branching topology and geometrical features of hepatic trees using a multiscale approach from the macro‐ down to the microcirculation. This may lead to novel insights into liver circulation, such as internal blood flow distributions and anatomical consequences of pathologies (e.g. cirrhosis).

The perfusion of the liver microcirculation is often analyzed in terms of idealized functional units (hexagonal liver lobules) based on a porous medium approach. More elaborate research is essential to assess the validity of this approach and to provide a more adequate and quantitative characterization of the liver microcirculation. To this end, we modeled the perfusion of the liver microcirculation using an image-based three-dimensional (3D) reconstruction of human liver sinusoids and computational fluid dynamics techniques. After vascular corrosion casting, a microvascular sample (±0.134 mm3) representing three liver lobules, was dissected from a human liver vascular replica and scanned using a high resolution (2.6 μm) micro-CT scanner. Following image processing, a cube (0.15 × 0.15 × 0.15 mm3) representing a sample of intertwined and interconnected sinusoids, was isolated from the 3D reconstructed dataset to define the fluid domain. Three models were studied to simulate flow along three orthogonal directions (i.e., parallel to the central vein and in the radial and circumferential directions of the lobule). Inflow and outflow guidances were added to facilitate solution convergence, and good quality volume meshes were obtained using approximately 9 × 106 tetrahedral cells. Subsequently, three computational fluid dynamics models were generated and solved assuming Newtonian liquid properties (viscosity 3.5 mPa s). Post-processing allowed to visualize and quantify the microvascular flow characteristics, to calculate the permeability tensor and corresponding principal permeability axes, as well as the 3D porosity. The computational fluid dynamics simulations provided data on pressure differences, preferential flow pathways and wall shear stresses. Notably, the pressure difference resulting from the flow simulation parallel to the central vein (0–100 Pa) was clearly smaller than the difference from the radial (0–170 Pa) and circumferential (0–180 Pa) flow directions. This resulted in a higher permeability along the central vein direction (kd,33 = 3.64 × 10−14 m2) in comparison with the radial (kd,11 = 1.56 × 10−14 m2) and circumferential (kd,22 = 1.75 × 10−14 m2) permeabilities which were approximately equal. The mean 3D porosity was 14.3. Our data indicate that the human hepatic microcirculation is characterized by a higher permeability along the central vein direction, and an about two times lower permeability along the radial and circumferential directions of a lobule. Since the permeability coefficients depend on the flow direction, (porous medium) liver microcirculation models should take into account sinusoidal anisotropy.

Hypothermic machine perfusion (HMP) is experiencing a revival in organ preservation due to the limitations of static cold storage and the need for better preservation of expanded criteria donor organs. For livers, perfusion protocols are still poorly defined, and damage of sinusoidal endothelial cells and heterogeneous perfusion are concerns. In this study, an electrical model of the human liver blood circulation is developed to enlighten internal pressure and flow distributions during HMP. Detailed vascular data on two human livers, obtained by combining vascular corrosion casting, micro-CT-imaging and image processing, were used to set up the electrical model. Anatomical data could be measured up to 5-6 vessel generations in each tree and showed exponential trend lines, used to predict data for higher generations. Simulated flow and pressure were in accordance with literature data. The model was able to simulate effects of pressure-driven HMP on liver hemodynamics and reproduced observations such as flow competition between the hepatic artery and portal vein. Our simulations further indicate that, from a pure biomechanical (shear stress) standpoint, HMP with low pressures should not result in organ damage, and that fluid viscosity has no effect on the shear stress experienced by the liver microcirculation in pressure-driven HMP.

The intraperitoneal (IP) administration of chemotherapy is an alternative treatment for peritoneal carcinomatosis, allowing for higher intratumor concentrations of the cytotoxic agent compared to intravenous administration. Nevertheless, drug penetration depths are still limited to a few millimeters. It is thus necessary to better understand the limiting factors behind this poor penetration in order to improve IP chemotherapy delivery. By developing a three-dimensional computational fluid dynamics (CFD) model for drug penetration in a tumor nodule, we investigated the impact of a number of key parameters on the drug transport and penetration depth during IP chemotherapy. Overall, smaller tumors showed better penetration than larger ones, which could be attributed to the lower IFP in smaller tumors. Furthermore, the model demonstrated large improvements in penetration depth by subjecting the tumor nodules to vascular normalization therapy, and illustrated the importance of the drug that is used for therapy. Explicitly modeling the necrotic core had a limited effect on the simulated penetration. Similarly, the penetration depth remained virtually constant when the Darcy permeability of the tissue changed. Our findings illustrate that the developed parametrical CFD model is a powerful tool providing more insight in the drug transport and penetration during IP chemotherapy.

The hepatic blood circulation is complex, particularly at the microcirculatory level. Previously, 2D liver lobule models using porous media and a 3D model using real sinusoidal geometries have been developed. We extended these models to investigate the role of vascular septa (VS) and anisotropic permeability. The lobule was modelled as a hexagonal prism (with or without VS) and the tissue was treated as a porous medium (isotropic or anisotropic permeability). Models were solved using computational fluid dynamics. VS inclusion resulted in more spatially homogeneous perfusion. Anisotropic permeability resulted in a larger axial velocity component than isotropic permeability. A parameter study revealed that results are most sensitive to the lobule size and radial pressure drop. Our model provides insight into hepatic microhaemodynamics, and suggests that inclusion of VS in the model leads to perfusion patterns that are likely to reflect physiological reality. The model has potential for applications to unphysiological and pathological conditions.