Simulation of haemodynamic flow in head and neck cancer chemotherapy

In silico intravascular optical coherence tomography (IVOCT) for quality assured imaging with reduced intervention

In the clinical application of intravascular optical coherence tomography (IVOCT), it is necessary to flush opaque blood during image acquisition. However, there are no specific standards for how to perform low-dose but effective flushing. In this study, computational fluid dynamics (CFD) and optical models were integrated to numerically simulate the complete process of IVOCT, which includes blood flushing with normal saline followed by image acquisition. Moreover, an intermittent injection scheme was proposed, and its advantages over the conventionally adopted scheme of continuous injection were verified. The results show that intermittent injection can significantly reduce the dosage of normal saline (reduced by 44.4%) with only a slight sacrifice of image quality (reduced by 8.7%, but still acceptable). The developed model and key findings in this work can help surgeons practice optimized IVOCT operations and potentially lead to improved designs of the IVOCT equipment.

Computational fluid dynamics study of intra-arterial chemotherapy for oral cancer

BACKGROUND

Intra-arterial chemotherapy (IAC) for oral cancer can deliver a higher concentration of anticancer agent into a tumor-feeding artery than intravenous systemic chemotherapy. However, distribution of anticancer agent into several branches of the external carotid artery (ECA) in IAC has not demonstrated sufficient treatment efficacy. To improve the effectiveness of IAC, the flow distribution of anticancer agent into the branches of the ECA in several IAC methods was investigated using computational fluid dynamics (CFD).

METHODS

Patient-specific three-dimensional vessel models were created from CT images of 2 patients with tongue cancer. Catheter models were combined with the vessel models. Thirty-two models were generated with varying vertical and horizontal positions of the catheter tip. With the use of a zero-dimensional resistance model of the peripheral vessel network, conventional IAC and superselective IAC were simulated in 30 and 2 models, respectively. The flow distribution of anticancer agent into the branches of the ECA was investigated in 32 models. Additionally, the blood streamline was traced from the inlet of the common carotid artery toward each outlet to examine the flow of anticancer agent in all models, and the wall shear stress of the vessel was calculated for some models.

RESULTS

The CFD simulations could be conducted within a reasonable computational time. In several models, the anticancer agent flowed into the target artery only when the catheter tip was located below the bifurcation of the ECA and each target artery. Furthermore, the anticancer agent tended to flow into the target artery when the catheter tip was shifted toward the target artery. In all ECA branches that had flow of anticancer agent, the blood streamlines to the target arteries contacted the catheter tip. Anticancer agent flowed into only the target artery in patients' models for superselective IAC. However, high wall shear stress was observed at the target artery in one patient's model.

CONCLUSIONS

This CFD study showed that location of the catheter tip was important in controlling the anticancer agent in conventional IAC. The distribution rate of anticancer agent into the tumor-feeding artery tended to increase when the catheter tip was placed below and toward the target artery. Although superselective IAC can reliably supply anticancer agent to the target artery, high wall shear stress at the target artery can occur, depending on vessel geometry of the patient, which may cause serious complications during the treatment.

Investigation of blood flow in the external carotid artery and its branches with a new 0D peripheral model

BACKGROUND

Patient-specific modelling in clinical studies requires a realistic simulation to be performed within a reasonable computational time. The aim of this study was to develop simple but realistic outflow boundary conditions for patient-specific blood flow simulation which can be used to clarify the distribution of the anticancer agent in intra-arterial chemotherapy for oral cancer.

METHODS

In this study, the boundary conditions are expressed as a zero dimension (0D) resistance model of the peripheral vessel network based on the fractal characteristics of branching arteries combined with knowledge of the circulatory system and the energy minimization principle. This resistance model was applied to four patient-specific blood flow simulations at the region where the common carotid artery bifurcates into the internal and external carotid arteries.

RESULTS

Results of these simulations with the proposed boundary conditions were compared with the results of ultrasound measurements for the same patients. The pressure was found to be within the physiological range. The difference in velocity in the superficial temporal artery results in an error of 5.21 ± 0.78 % between the numerical results and the measurement data.

CONCLUSIONS

The proposed outflow boundary conditions, therefore, constitute a simple resistance-based model and can be used for performing accurate simulations with commercial fluid dynamics software.