Liver cancer arterial perfusion modelling and CFD boundary conditions methodology: a case study of the haemodynamics of a patient-specific hepatic artery in literature-based healthy and tumour-bearing liver scenarios

A Hybrid Particle-Flow CFD Modeling Approach in Truncated Hepatic Arterial Trees for Liver Radioembolization: A Patient-specific Case Study

Hepatocellular carcinoma (HCC) is the most common form of primary liver cancer. At its intermediate, unresectable stage, HCC is typically treated by local injection of embolizing microspheres in the hepatic arteries to selectively damage tumor tissue. Interestingly, computational fluid dynamics (CFD) has been applied increasingly to elucidate the impact of clinically variable parameters, such as injection location, on the downstream particle distribution. This study aims to reduce the computational cost of such CFD approaches by introducing a novel truncation algorithm to simplify hepatic arterial trees, and a hybrid particle-flow modeling approach which only models particles in the first few bifurcations. A patient-specific hepatic arterial geometry was pruned at three different levels, resulting in three trees: Geometry 1 (48 outlets), Geometry 2 (38 outlets), and Geometry 3 (17 outlets). In each geometry, 1 planar injection and 3 catheter injections (each with different tip locations) were performed. For the truncated geometries, it was assumed that, downstream of the truncated outlets, particles distributed themselves proportional to the blood flow. This allowed to compare the particle distribution in all 48 "outlets" for each geometry. For the planar injections, the median difference in outlet-specific particle distribution between Geometry 1 and 3 was 0.21%; while the median difference between outlet-specific flow and particle distribution in Geometry 1 was 0.40%. Comparing catheter injections, the maximum median difference in particle distribution between Geometry 1 and 3 was 0.24%, while the maximum median difference between particle and flow distribution was 0.62%. The results suggest that the hepatic arterial tree might be reliably truncated to estimate the particle distribution in the full-complexity tree. In the resulting hybrid particle-flow model, explicit particle modeling was only deemed necessary in the first few bifurcations of the arterial tree. Interestingly, using flow distribution as a surrogate for particle distribution in the entire tree was considerably less accurate than using the hybrid model, although the difference was much higher for catheter injections than for planar injections. Future work should focus on replicating and experimentally validating these results in more patient-specific geometries.

The Impact of Injection Distance to Bifurcations on Yttrium-90 Distribution in Liver Cancer Radioembolization

PURPOSE

To model the effect of the injection location on the distribution of yttrium-90 (⁹⁰Y) microspheres in the liver during radioembolization using computational simulation and to determine the potential effects of radial movements of the catheter tip.

MATERIALS AND METHODS

Numerical studies were conducted using images from a representative patient with hepatocellular carcinoma. The right hepatic artery (RHA) was segmented from contrast-enhanced cone-beam computed tomography scans. The blood flow was investigated in the trunk of the RHA using numerical simulations for 6 injection position scenarios at 2 sites located at a distance of approximately 5 and 20 mm upstream of the first bifurcation (RHA diameters of approximately 4.6 mm). The ⁹⁰Y delivery to downstream vessels was calculated from the simulated hepatic artery hemodynamics.

RESULTS

Varying the injection location along the RHA and across the vessel cross-section resulted in different simulated microsphere distributions in the downstream vascular bed. When the catheter tip was 5 mm upstream of the bifurcation, ⁹⁰Y distribution in the downstream branches varied by as much as 53% with a 1.5-mm radial movement of the tip. However, the catheter radial movement had a weaker effect on the microsphere distribution when the injection plane was farther from the first bifurcation (20 mm), with a maximum delivery variation of 9% to a downstream branch.

CONCLUSIONS

An injection location far from bifurcations is recommended to minimize the effect of radial movements of the catheter tip on the microsphere distribution.

Larraona G, Sangro B, Bilbao JI, Aramburu J. Computational study of a novel catheter for liver radioembolization

Radioembolization (RE) is a medical treatment for primary and secondary liver cancer that involves the transcatheter intraarterial delivery of micron-sized and radiation-emitting microspheres, with the goal of improving microsphere deposition in the tumoral bed while sparing healthy tissue. An increasing number of in vitro and in silico studies on RE in the literature suggest that the particle injection velocity, spatial location of the catheter tip and catheter type are important parameters in particle distribution. The present in silico study assesses the performance of a novel catheter design that promotes particle dispersion near the injection point, with the goal of generating a particle distribution that mimics the flow split to facilitate tumour targeting. The design is based on two factors: the direction and the velocity at which particles are released from the catheter. A series of simulations was performed with the catheter inserted at an idealised hepatic artery tree with physiologically realistic boundary conditions. Two longitudinal microcatheter positions in the first generation of the tree were studied by analysing the performance of the catheter in terms of the outlet-to-outlet particle distribution and split flow matching. The results show that the catheter with the best performance is one with side holes on the catheter wall and a closed frontal tip. This catheter promotes a flow-split-matching particle distribution, which improves as the injection crossflow increases.

Computational Fluid Dynamics Modeling of Liver Radioembolization: A Review

Yttrium-90 radioembolization (RE) is a widely used transcatheter intraarterial therapy for patients with unresectable liver cancer. In the last decade, computer simulations of hepatic artery hemodynamics during RE have been performed with the aim of better understanding and improving the therapy. In this review, we introduce the concept of computational fluid dynamics (CFD) modeling with a clinical perspective and we review the CFD models used to study RE from the fluid mechanics point of view. Finally, we show what CFD simulations have taught us about the hemodynamics during RE, the current capabilities of CFD simulations of RE, and we suggest some future perspectives.

CFD Simulations of Radioembolization: A Proof-of-Concept Study on the Impact of the Hepatic Artery Tree Truncation

Radioembolization (RE) is a treatment for patients with liver cancer, one of the leading cause of cancer-related deaths worldwide. RE consists of the transcatheter intraarterial infusion of radioactive microspheres, which are injected at the hepatic artery level and are transported in the bloodstream, aiming to target tumors and spare healthy liver parenchyma. In paving the way towards a computer platform that allows for a treatment planning based on computational fluid dynamics (CFD) simulations, the current simulation (model preprocess, model solving, model postprocess) times (of the order of days) make the CFD-based assessment non-viable. One of the approaches to reduce the simulation time includes the reduction in size of the simulated truncated hepatic artery. In this study, we analyze for three patient-specific hepatic arteries the impact of reducing the geometry of the hepatic artery on the simulation time. Results show that geometries can be efficiently shortened without impacting greatly on the microsphere distribution.

A proof-of-concept study of the in-vivo validation of a computational fluid dynamics model of personalized radioembolization

Radioembolization (RE) with yttrium-90 (90Y) microspheres, a transcatheter intraarterial therapy for patients with liver cancer, can be modeled computationally. The purpose of this work was to correlate the results obtained with this methodology using in vivo data, so that this computational tool could be used for the optimization of the RE procedure. The hepatic artery three-dimensional (3D) hemodynamics and microsphere distribution during RE were modeled for six 90Y-loaded microsphere infusions in three patients with hepatocellular carcinoma using a commercially available computational fluid dynamics (CFD) software package. The model was built based on in vivo data acquired during the pretreatment stage. The results of the simulations were compared with the in vivo distribution assessed by 90Y PET/CT. Specifically, the microsphere distribution predicted was compared with the actual 90Y activity per liver segment with a commercially available 3D-voxel dosimetry software (PLANET Dose, DOSIsoft). The average difference between the CFD-based and the PET/CT-based activity distribution was 2.36 percentage points for Patient 1, 3.51 percentage points for Patient 2 and 2.02 percentage points for Patient 3. These results suggest that CFD simulations may help to predict 90Y-microsphere distribution after RE and could be used to optimize the RE procedure on a patient-specific basis.

Multiscale Computational Fluid Dynamics Modeling for Personalized Liver Cancer Radioembolization Dosimetry

Yttrium-90 (90Y) radioembolization is a minimally invasive procedure increasingly used for advanced liver cancer treatment. In this method, radioactive microspheres are injected into the hepatic arterial bloodstream to target, irradiate, and kill cancer cells. Accurate and precise treatment planning can lead to more efficient and safer treatment by delivering a higher radiation dose to the tumor while minimizing the exposure of the surrounding liver parenchyma. Treatment planning primarily relies on the estimated radiation dose delivered to tissue. However, current methods used to estimate the dose are based on simplified assumptions that make the dosimetry results unreliable. In this work, we present a computational model to predict the radiation dose from the 90Y activity in different liver segments to provide a more realistic and personalized dosimetry. Computational fluid dynamics (CFD) simulations were performed in a 3D hepatic arterial tree model segmented from cone-beam CT angiographic data obtained from a patient with hepatocellular carcinoma (HCC). The microsphere trajectories were predicted from the velocity field. 90Y dose distribution was then calculated from the volumetric distribution of the microspheres. Two injection locations were considered for the microsphere administration, a lobar and a selective injection. Results showed that 22% and 82% of the microspheres were delivered to the tumor, after each injection, respectively, and the combination of both injections ultimately delivered 49% of the total administered 90Y microspheres to the tumor. Results also illustrated the nonhomogeneous distribution of microspheres between liver segments, indicating the importance of developing patient-specific dosimetry methods for effective radioembolization treatment.

Transarterial drug delivery for liver cancer: numerical simulations and experimental validation of particle distribution in patient-specific livers

Background: Transarterial therapies are routinely used for the locoregional treatment of unresectable hepatocellular carcinoma (HCC). However, the impact of clinical parameters (i.e. injection location, particle size, particle density etc.) and patient-specific conditions (i.e. hepatic geometry, cancer burden) on the intrahepatic particle distribution (PD) after transarterial injection of embolizing microparticles is still unclear. Computational fluid dynamics (CFD) may help to better understand this impact.Methods: Using CFD, both the blood flow and microparticle mass transport were modeled throughout the 3D-reconstructed arterial vasculature of a patient-specific healthy and cirrhotic liver. An experimental feasibility study was performed to simulate the PD in a 3D-printed phantom of the cirrhotic arterial network.Results: Axial and in-plane injection locations were shown to be effective parameters to steer particles toward tumor tissue in both geometries. Increasing particle size or density made it more difficult for particles to exit the domain. As cancer burden increased, the catheter tip location mattered less. The in vitro study and numerical results confirmed that PD largely mimics flow distribution, but that significant differences are still possible.Conclusions: Our findings highlight that optimal parameter choice can lead to selective targeting of tumor tissue, but that targeting potential highly depends on patient-specific conditions.

Estimation of Yttrium-90 Distribution in Liver Radioembolization using Computational Fluid Dynamics and Deep Neural Networks

Yttrium-90 (⁹⁰Y) radioembolization is a liver cancer therapy based on ⁹⁰Y microspheres injected into the hepatic artery. Current dosimetry methods used to estimate the absorbed dose in order to prescribe the ⁹⁰Y activity to inject are not accurate, which can affect the treatment effectiveness. A new dosimetry based on the hemodynamics simulation of the hepatic arterial tree, CFDose, aimed at overcoming some of the limitations of the current methods. However, due to the expensive computational cost of computational fluid dynamics (CFD) simulations, this method needs to be accelerated before it can be used in real-time during treatment planning. In this paper, we introduce a convolutional neural network model trained with the CFD results of a patient with hepatocellular carcinoma to predict the ⁹⁰Y distribution under different downstream vasculature resistance conditions. The model performance was evaluated using two metrics, the mean squared error and prediction accuracy. The prediction accuracy showed that the average difference between the actual and predicted data was less than 1%. The proposed model could estimate the ⁹⁰Y distribution significantly faster than a CFD simulation.

Roudsari B, Vu CT, Roncali E. Computational Modeling of the Liver Arterial Blood Flow for Microsphere Therapy: Effect of Boundary Conditions

Transarterial embolization is a minimally invasive treatment for advanced liver cancer using microspheres loaded with a chemotherapeutic drug or radioactive yttrium-90 (90Y) that are injected into the hepatic arterial tree through a catheter. For personalized treatment, the microsphere distribution in the liver should be optimized through the injection volume and location. Computational fluid dynamics (CFD) simulations of the blood flow in the hepatic artery can help estimate this distribution if carefully parameterized. An important aspect is the choice of the boundary conditions imposed at the inlet and outlets of the computational domain. In this study, the effect of boundary conditions on the hepatic arterial tree hemodynamics was investigated. The outlet boundary conditions were modeled with three-element Windkessel circuits, representative of the downstream vasculature resistance. Results demonstrated that the downstream vasculature resistance affected the hepatic artery hemodynamics such as the velocity field, the pressure field and the blood flow streamline trajectories. Moreover, the number of microspheres received by the tumor significantly changed (more than 10% of the total injected microspheres) with downstream resistance variations. These findings suggest that patient-specific boundary conditions should be used in order to achieve a more accurate drug distribution estimation with CFD in transarterial embolization treatment planning.

On the importance of spiral-flow inflow boundary conditions when using idealized artery geometries in the analysis of liver radioembolization: A parametric study

In the last decades, the numerical studies on hemodynamics have become a valuable explorative scientific tool. The very first studies were done over idealized geometries, but as numerical methods and the power of computers have become more affordable, the studies tend to be patient specific. We apply the study to the numerical analysis of tumor-targeting during liver radioembolization (RE). RE is a treatment for liver cancer, and is performed by injecting radiolabeled microspheres via a catheter placed in the hepatic artery. The objective of the procedure is to maximize the release of radiolabeled microspheres into the tumor and avoid a healthy tissue damage. Idealized virtual arteries can serve as a generalist approach that permits to separately analyze the effect of a variable in the microsphere distribution with respect to others. However, it is important to use proper physiological boundary conditions (BCs). It is not obvious, the need to account for the effect of tortuosity when using an idealized virtual artery. We study the use of idealized geometry of a hepatic artery as a valid research tool, exploring the importance of using realistic spiral-flow inflow BC. By using a literature-based cancer scenario, we vary two parameters to analyze the microsphere distribution through the outlets of the geometry. The parameters varied are the type of microspheres injected and the microsphere injection velocity. The results with realistic inlet velocity profile showed that the particle distribution in the liver segments is not affected by the analyzed injection velocity values neither by the particle density. NOVELTY STATEMENT: In this article, we assessed the use of idealized geometries as a valid research tool and applied the use of an idealized geometry to the case of an idealized hepatic artery to study the particle-hemodynamics during radioembolization (RE). We studied three different inflow boundary conditions (BCs) to assess the usefulness of the geometry, two types of particle injection velocities and two types of commercially available microspheres for RE treatment. In recent years, the advent in computational resources allowed for more detailed patient-specific geometry generation and discretization and hemodynamics simulations. However, general studies based on idealized geometries can be performed in order to provide medical doctors with some basic and general guidelines when using a given catheter for a given cancer scenario. Moreover, using an idealized geometry can be a reasonable approach which allows us to isolate a given parameter and control other parameters, so that parameters can be independently assessed. Even though an idealized geometry does not match any patient's geometry, the use of an idealized geometry can be valid when drawing general conclusions that may be useful in patient-specific cases. However, we believe that even if an idealized hepatic artery geometry is used for the study, it is necessary to account for the upstream and downstream tortuosity of vessels through the BCs. In this work, we highlighted the need of modeling the tortuosity of upstream and downstream vasculatures through the BCs.

Modeling the hepatic arterial flow in living liver donor after left hepatectomy and postoperative boundary condition exploration

Preoperative and postoperative hepatic perfusion is modeled with one-dimensional (1-D) Navier-Stokes equations. Flow rates obtained from ultrasound (US) data and impedance resulted from structured trees are the inflow and outflow boundary condition (BC), respectively. Structured trees terminate at the size of the arterioles, which can enlarge their size after hepatectomy. In clinical studies, the resistance to pulsatile arterial flow caused by the microvascular bed can be reflected by the resistive index (RI), a frequently used index in assessing arterial resistance. This study uses the RI in a novel manner to conveniently obtain the postoperative outflow impedance from the preoperative impedance. The major emphasis of this study is to devise a model to capture the postoperative hepatic hemodynamics after left hepatectomy. To study this, we build a hepatic network model and analyze its behavior under four different outflow impedance: (a) the same as preoperative impedance; (b) evaluated using the RI and preoperative impedance; (c) computed from structured tree BC with increased radius of terminal vessels; and (d) evaluated using structured tree with both increased radius of root vessel, ie, the outlets of the postoperative hepatic artery, and increased radius of terminal vessels. Our results show that both impedance from (b) and (d) give a physiologically reasonable postoperative hepatic pressure range, while the RI in (b) allows for a fast approximation of postoperative impedance. Since hemodynamics after hepatectomy are not fully understood, the methods used in this study to explore postoperative outflow BC are informative for future models exploring hemodynamic effects of partial hepatectomy.

A methodology for numerically analysing the hepatic artery haemodynamics during B-TACE: a proof of concept

Balloon-occluded transarterial chemoembolisation (B-TACE) is an intraarterial transcatheter treatment for liver cancer. In B-TACE, an artery-occluding microballoon catheter occludes an artery and promotes collateral circulation for drug delivery to tumours. This paper presents a methodology for analysing the haemodynamics during B-TACE, by combining zero-dimensional and three-dimensional modelling tools. As a proof of concept, we apply the methodology to a patient-specific hepatic artery geometry and analyse two catheter locations. Results show that the blood flow redistribution can be predicted in this proof-of-concept study, suggesting that this approach could potentially be used to optimise catheter location.

Numerical zero-dimensional hepatic artery hemodynamics model for balloon-occluded transarterial chemoembolization

Balloon-occluded transarterial chemoembolization (B-TACE) is a valuable treatment option for patients with inoperable malignant tumors in the liver. Balloon-occluded transarterial chemoembolization consists of the transcatheter infusion of an anticancer drug mixture and embolic agents. Contrary to conventional TACE, B-TACE is performed via an artery-occluding microballoon catheter, which makes the blood flow to redistribute due to the intra- and extrahepatic arterial collateral circulation. Several recent studies have stressed the importance of the redistribution of blood flow in enhancing the treatment outcome. In the present study, the geometries of a representative hepatic artery and the communicating arcades (CAs) are modeled. An in silico zero-dimensional hemodynamic model is created by characterizing the geometry and the boundary conditions and then is validated in vitro. The role of CAs is assessed by combining 2 cancer scenarios and 2 catheter locations. The importance of the diameter of the CAs is also studied. Results show that occluding a main artery leads to collateral circulation and CAs start to play a role in blood-flow redistribution. In summary, numerical zero-dimensional simulations permit a fast and reliable approach for exploring the blood-flow redistribution caused by the occlusion of a main artery, and this approach could be used during B-TACE planning.

The role of angled-tip microcatheter and microsphere injection velocity in liver radioembolization: A computational particle-hemodynamics study

Liver radioembolization is a promising treatment option for combating liver tumors. It is performed by placing a microcatheter in the hepatic artery and administering radiation-emitting microspheres through the arterial bloodstream so that they get lodged in the tumoral bed. In avoiding nontarget radiation, the standard practice is to conduct a pretreatment, in which the microcatheter location and injection velocity are decided. However, between pretreatment and actual treatment, some of the parameters that influence the particle distribution in the liver can vary, resulting in radiation-induced complications. The present study aims to analyze the influence of a commercially available microcatheter with an angled tip and particle injection velocity in terms of segment-to-segment particle distribution. Specifically, 4 tip orientations and 2 injection velocities are combined to yield a set of 8 numerical simulations of the particle-hemodynamics in a patient-specific truncated hepatic artery. For each simulation, 4 cardiac pulses are simulated. Particles are injected during the first cycle, and the remaining pulses enable the majority of the injected particles to exit the computational domain. Results indicate that, in terms of injection velocity, particles are more spread out in the cross-sectional lumen areas as the injection velocity increases. The tip's orientation also plays a role because it influences the near-tip hemodynamics, therefore altering the particle travel through the hepatic artery. However, results suggest that particle distribution tries to match the blood flow split, therefore particle injection velocity and microcatheter tip orientation playing a minor role in segment-to-segment particle distribution.

Computational Modeling in Liver Surgery

The need for extended liver resection is increasing due to the growing incidence of liver tumors in aging societies. Individualized surgical planning is the key for identifying the optimal resection strategy and to minimize the risk of postoperative liver failure and tumor recurrence. Current computational tools provide virtual planning of liver resection by taking into account the spatial relationship between the tumor and the hepatic vascular trees, as well as the size of the future liver remnant. However, size and function of the liver are not necessarily equivalent. Hence, determining the future liver volume might misestimate the future liver function, especially in cases of hepatic comorbidities such as hepatic steatosis. A systems medicine approach could be applied, including biological, medical, and surgical aspects, by integrating all available anatomical and functional information of the individual patient. Such an approach holds promise for better prediction of postoperative liver function and hence improved risk assessment. This review provides an overview of mathematical models related to the liver and its function and explores their potential relevance for computational liver surgery. We first summarize key facts of hepatic anatomy, physiology, and pathology relevant for hepatic surgery, followed by a description of the computational tools currently used in liver surgical planning. Then we present selected state-of-the-art computational liver models potentially useful to support liver surgery. Finally, we discuss the main challenges that will need to be addressed when developing advanced computational planning tools in the context of liver surgery.

Prevention and treatment of complications of selective internal radiation therapy: Expert guidance and systematic review

Selective internal radiation therapy (or radioembolization) by intra-arterial injection of radioactive yttrium-90-loaded microspheres is increasingly used for the treatment of patients with liver metastases or primary liver cancer. The high-dose beta-radiation penetrates an average of only 2.5 mm from the source, thus limiting its effects to the site of delivery. However, the off-target diversion of yttrium-90 microspheres to tissues other than the tumor may lead to complications. The most prominent of these complications include radiation gastritis and gastrointestinal ulcers, cholecystitis, radiation pneumonitis, and radioembolization-induced liver disease, which may occur despite careful pretreatment planning. Thus, selective internal radiation therapy demands an expert multidisciplinary team approach in order to provide comprehensive care for patients. This review provides recommendations to multidisciplinary teams on the optimal medical processes in order to ensure the safe delivery of selective internal radiation therapy. Based on the best available published evidence and expert opinion, we recommend the most appropriate strategies for the prevention, early diagnosis, and management of potential radiation injury to the liver and to other organs. (Hepatology 2017;66:969-982).

Computational particle-haemodynamics analysis of liver radioembolization pretreatment as an actual treatment surrogate

Liver radioembolization (RE) is a treatment option for patients with unresectable and chemorefractory primary and metastatic liver tumours. RE consists of intra-arterially administering via catheter radioactive microspheres that locally attack the tumours, sparing healthy tissue. Prior to RE, the standard practice is to conduct a treatment-mimicking pretreatment assessment via the infusion of ^99m Tc-labelled macroaggregated albumin microparticles. The usefulness of this pretreatment has been debated in the literature, and thus, the aim of the present study is to shed light on this issue by numerically simulating the liver RE pretreatment and actual treatment particle-haemodynamics in a patient-specific hepatic artery under two different literature-based cancer scenarios and two different placements of a realistic end-hole microcatheter in the proper hepatic artery. The parameters that are analysed are the following: microagent quantity and size (accounting for RE pretreatment and treatment), catheter-tip position (near the proper hepatic artery bifurcation and away from it), and cancer burden (10% and 30% liver involvement). The conclusion that can be reached from the simulations is that when it comes to mimicking RE in terms of delivering particles to tumour-bearing segments, the catheter-tip position is much more important (because of the importance of local haemodynamic pattern alteration) than the infused microagents (i.e. quantity and size). Cancer burden is another important feature because the increase in blood flow rate to tumour-bearing segments increases the power to drag particles. These numerical simulation-based conclusions are in agreement with clinically observed events reported in the literature. Copyright © 2016 John Wiley & Sons, Ltd.