In vivo mimicking model for solid tumor towards hydromechanics of tissue deformation and creation of necrosis

The present work addresses transvascular and interstitial fluid transport inside a solid tumor surrounded by normal tissue (close to an in vivo mimicking setup). In general, biological tissues behave like a soft porous material and show mechanical behavior towards the fluid motion through the interstitial space. In general, forces like viscous drag that are associated with the fluid flow may compress the tissue material. On the macroscopic level, we try to model the motion of fluids and macromolecules through the interstitial space of solid tumor and the normal tissue layer. The transvascular fluid transport is assumed to be governed by modified Starling’s law. The poroelastohydrodynamics (interstitial hydrodynamics and the deformation of tissue material) inside the tumor and normal tissue regions is modeled using linearized biphasic mixture theory. Correspondingly, the velocity distribution of fluid is coupled to the displacement field of the solid phase (mainly cellular phase and extracellular matrix) in both the normal and tumor tissue regions. The corresponding velocity field is used within the transport reaction equation for fluids and macromolecules through interstitial space to get the overall solute (e.g., nutrients, drug, and other macromolecules) distribution. This study justifies that the presence of the normal tissue layer plays a significant role in delaying/assisting necrosis inside the tumor tissue. It is observed that the exchange process of fluids and macromolecules across the interface of the tumor and normal tissue affects the effectiveness factor corresponding to the tumor tissue.     

Evaluation of a Voxelized Model Based on DCE-MRI for Tracer Transport in Tumor

Recent advances in the treatment of cancer involving therapeutic agents have shown promising results. However, treatment efficacy can be limited due to inadequate and uneven uptake in solid tumors, thereby making the prediction of drug transport important for developing effective therapeutic strategies. In this study, a patient-specific computational porous media model (voxelized model) was developed for predicting the interstitial flow field and distribution of a systemically delivered magnetic resonance (MR) visible tracer in a tumor. The benefits of a voxel approach include less labor and less computational time (approximately an order of magnitude reduction compared to the traditional computational fluid dynamics (CFD) approach developed earlier by our group). The model results were compared with that obtained from a previous approach based on unstructured meshes along with MR-measured tracer concentration data within tumors, using statistical analysis and qualitative representations. The statistical analysis indicated the similarity between the structured and unstructured models' results with a low root mean square error (RMS) and a high correlation coefficient. The voxelized model captured features of the flow field and tracer distribution such as high interstitial fluid pressure inside the tumor and the heterogeneous distribution of the tracer. Predictions of tracer distribution by the voxelized approach also resulted in low RMS error when compared with MR-measured data over a 1?h time course. The similarity in the voxelized model results with experiment and the nonvoxelized model predictions were maintained across three different tumors. Overall, the voxelized model serves as a reliable and swift alternative to approaches using unstructured meshes in predicting extracellular transport within tumors.     

Numerical modeling of fluid flow in solid tumors

A mathematical model of interstitial fluid flow is developed, based on the application of the governing equations for fluid flow, i.e., the conservation laws for mass and momentum, to physiological systems containing solid tumors. The discretized form of the governing equations, with appropriate boundary conditions, is developed for a predefined tumor geometry. The interstitial fluid pressure and velocity are calculated using a numerical method, element based finite volume. Simulations of interstitial fluid transport in a homogeneous solid tumor demonstrate that, in a uniformly perfused tumor, i.e., one with no necrotic region, because of the interstitial pressure distribution, the distribution of drug particles is non-uniform. Pressure distribution for different values of necrotic radii is examined and two new parameters, the critical tumor radius and critical necrotic radius, are defined. Simulation results show that: 1) tumor radii have a critical size. Below this size, the maximum interstitial fluid pressure is less than what is generally considered to be effective pressure (a parameter determined by vascular pressure, plasma osmotic pressure, and interstitial osmotic pressure). Above this size, the maximum interstitial fluid pressure is equal to effective pressure. As a consequence, drugs transport to the center of smaller tumors is much easier than transport to the center of a tumor whose radius is greater than the critical tumor radius; 2) there is a critical necrotic radius, below which the interstitial fluid pressure at the tumor center is at its maximum value. If the tumor radius is greater than the critical tumor radius, this maximum pressure is equal to effective pressure. Above this critical necrotic radius, the interstitial fluid pressure at the tumor center is below effective pressure. In specific ranges of these critical sizes, drug amount and therefore therapeutic effects are higher because the opposing force, interstitial fluid pressure, is low in these ranges.     

Computational Modeling of 3D Tumor Growth and Angiogenesis for Chemotherapy Evaluation

Solid tumors develop abnormally at spatial and temporal scales, giving rise to biophysical barriers that impact anti-tumor chemotherapy. This may increase the expenditure and time for conventional drug pharmacokinetic and pharmacodynamic studies. In order to facilitate drug discovery, we propose a mathematical model that couples three-dimensional tumor growth and angiogenesis to simulate tumor progression for chemotherapy evaluation. This application-oriented model incorporates complex dynamical processes including cell- and vascular-mediated interstitial pressure, mass transport, angiogenesis, cell proliferation, and vessel maturation to model tumor progression through multiple stages including tumor initiation, avascular growth, and transition from avascular to vascular growth. Compared to pure mechanistic models, the proposed empirical methods are not only easy to conduct but can provide realistic predictions and calculations. A series of computational simulations were conducted to demonstrate the advantages of the proposed comprehensive model. The computational simulation results suggest that solid tumor geometry is related to the interstitial pressure, such that tumors with high interstitial pressure are more likely to develop dendritic structures than those with low interstitial pressure.     

Drug delivery to brain tumors

We develop a macroscopic model for delivering drug to brain tumors. The model accounts for bulk convective and diffusive transport across the blood-brain barrier and through the interstitial space. Through mathematical analysis and simulations, we assess the effects of changing parameters (within physiological bounds) on drug delivery. We find that there is an optimal treatment for convective drug delivery to the center of the tumor. We interpret this phenomenon in terms of traffic flow. The implications of our analyses on existing chemotherapeutic protocols are discussed.     

A Mathematical Model of the Enhanced Permeability and Retention Effect for Liposome Transport in Solid Tumors

The discovery of the enhanced permeability and retention (EPR) effect has resulted in the development of nanomedicines, including liposome-based formulations of drugs, as cancer therapies. The use of liposomes has resulted in substantial increases in accumulation of drugs in solid tumors; yet, significant improvements in therapeutic efficacy have yet to be achieved. Imaging of the tumor accumulation of liposomes has revealed that this poor or variable performance is in part due to heterogeneous inter-subject and intra-tumoral liposome accumulation, which occurs as a result of an abnormal transport microenvironment. A mathematical model that relates liposome accumulation to the underlying transport properties in solid tumors could provide insight into inter and intra-tumoral variations in the EPR effect. In this paper, we present a theoretical framework to describe liposome transport in solid tumors. The mathematical model is based on biophysical transport equations that describe pressure driven fluid flow across blood vessels and through the tumor interstitium. The model was validated by direct comparison with computed tomography measurements of tumor accumulation of liposomes in three preclinical tumor models. The mathematical model was fit to liposome accumulation curves producing predictions of transport parameters that reflect the tumor microenvironment. Notably, all fits had a high coefficient of determination and predictions of interstitial fluid pressure agreed with previously published independent measurements made in the same tumor type. Furthermore, it was demonstrated that the model attributed inter-subject heterogeneity in liposome accumulation to variations in peak interstitial fluid pressure. These findings highlight the relationship between transvascular and interstitial flow dynamics and variations in the EPR effect. In conclusion, we have presented a theoretical framework that predicts inter-subject and intra-tumoral variations in the EPR effect based on fundamental properties of the tumor microenvironment and forms the basis for transport modeling of liposome drug delivery.     

Drug delivery patterns for different stenting techniques in coronary bifurcations: a comparative computational study

The treatment of coronary bifurcation lesions represents a challenge for the interventional cardiologists due to the lower rate of procedural success and the higher risk of restenosis. The advent of drug-eluting stents (DES) has dramatically reduced restenosis and consequently the request for re-intervention. The aim of the present work is to provide further insight about the effectiveness of DES by means of a computational study that combines virtual stent implantation, fluid dynamics and drug release for different stenting protocols currently used in the treatment of a coronary artery bifurcation. An explicit dynamic finite element model is developed in order to obtain realistic configurations of the implanted devices used to perform fluid dynamics analysis by means of a previously developed finite element method coupling the blood flow and the intramural plasma filtration in rigid arteries. To efficiently model the drug release, a multiscale strategy is adopted, ranging from lumped parameter model accounting for drug release to fully 3-D models for drug transport to the artery. Differences in drug delivery to the artery are evaluated with respect to local drug dosage. This model allowed to compare alternative stenting configurations (namely the Provisional Side Branch, the Culotte and the Inverted Culotte techniques), thus suggesting guidelines in the treatment of coronary bifurcation lesions and addressing clinical issues such as the effectiveness of drug delivery to lesions in the side branch, as well as the influence of incomplete strut apposition and overlapping stents.     

Analysis of resonance chemotherapy in leukemia treatment via multi-staged population balance models

An age-structured population balance model that explicitly models cell cycle phases is developed to investigate the effects of cell cycle specific (CCS) drugs. In particular, the benefits of timing CCS drug treatments in resonance chemotherapy are predicted and measured directly in vitro before evaluating likely in vivo scenarios. The phase transition rates are measured in vitro for the HL60 leukemia cell line and are used to predict the transient phase dynamics after exposure to the S phase specific drug, camptothecin. The phase oscillations predicted by the model are observed experimentally and the timing of a second camptothecin pulse is shown to significantly alter the overall treatment effectiveness. To explore the feasibility of designing resonance chemotherapeutic treatments to preferentially eliminate one cell type over another, Jurkat and HL60 leukemia cells are exposed to the same dual-pulse camptothecin treatment regimen. With the model framework validated for simplified cases, the model is used to extrapolate the effectiveness of resonance chemotherapy considering in vivo effects such as quiescence, drug metabolism, drug properties, and transport considerations that were not included in the in vitro experiments. While resonance chemotherapy is intuitive and looks promising in vitro, when in vivo considerations are included in the model, the phenomenon is dampened and the window of applicability becomes narrower.     

Effect of tumor shape,size, and tissue transport properties on drug delivery to solid tumors

Background

The computational methods provide condition for investigation related to the process of drug delivery, such as convection and diffusion of drug in extracellular matrices, drug extravasation from microvessels or to lymphatic vessels. The information of this process clarifies the mechanisms of drug delivery from the injection site to absorption by a solid tumor. In this study, an advanced numerical method is used to solve fluid flow and solute transport equations simultaneously to investigate the effect of tumor shape and size on drug delivery to solid tumor.

Methods

The advanced mathematical model used in our previous work is further developed by adding solute transport equation to the governing equations. After applying appropriate boundary and initial conditions on tumor and surrounding tissue geometry, the element-based finite volume method is used for solving governing equations of drug delivery in solid tumor. Also, the effects of size and shape of tumor and some of tissue transport parameters such as effective pressure and hydraulic conductivity on interstitial fluid flow and drug delivery are investigated.

Results

Sensitivity analysis shows that drug delivery in prolate shape is significantly better than other tumor shapes. Considering size effect, increasing tumor size decreases drug concentration in interstitial fluid. This study shows that dependency of drug concentration in interstitial fluid to osmotic and intravascular pressure is negligible.

Conclusions

This study shows that among diffusion and convection mechanisms of drug transport, diffusion is dominant in most different tumor shapes and sizes. In tumors in which the convection has considerable effect, the drug concentration is larger than that of other tumors at the same time post injection.

     

Aerosol deposition characteristics in distal acinar airways under cyclic breathing conditions

Although the major mechanisms of aerosol deposition in the lung are known, detailed quantitative data in anatomically realistic models are still lacking, especially in the acinar airways. In this study, an algorithm was developed to build multigenerational three-dimensional models of alveolated airways with arbitrary bifurcation angles and spherical alveolar shape. Using computational fluid dynamics, the deposition of 1- and 3-μm aerosol particles was predicted in models of human alveolar sac and terminal acinar bifurcation under rhythmic wall motion for two breathing conditions (functional residual capacity = 3 liter, tidal volume = 0.5 and 0.9 liter, breathing period = 4 s). Particles entering the model during one inspiration period were tracked for multiple breathing cycles until all particles deposited or escaped from the model. Flow recirculation inside alveoli occurred only during transition between inspiration and expiration and accounted for no more than 1% of the whole cycle. Weak flow irreversibility and convective transport were observed in both models. The average deposition efficiency was similar for both breathing conditions and for both models. Under normal gravity, total deposition was ~33 and 75%, of which ~67 and 96% occurred during the first cycle, for 1- and 3-μm particles, respectively. Under zero gravity, total deposition was ~2-5% for both particle sizes. These results support previous findings that gravitational sedimentation is the dominant deposition mechanism for micrometer-sized aerosols in acinar airways. The results also showed that moving walls and multiple breathing cycles are needed for accurate estimation of aerosol deposition in acinar airways.     

Numerical Modeling of Interstitial Fluid Flow Coupled with Blood Flow through a Remodeled Solid Tumor Microvascular Network

Modeling of interstitial fluid flow involves processes such as fluid diffusion, convective transport in extracellular matrix, and extravasation from blood vessels. To date, majority of microvascular flow modeling has been done at different levels and scales mostly on simple tumor shapes with their capillaries. However, with our proposed numerical model, more complex and realistic tumor shapes and capillary networks can be studied. Both blood flow through a capillary network, which is induced by a solid tumor, and fluid flow in tumor’s surrounding tissue are formulated. First, governing equations of angiogenesis are implemented to specify the different domains for the network and interstitium. Then, governing equations for flow modeling are introduced for different domains. The conservation laws for mass and momentum (including continuity equation, Darcy’s law for tissue, and simplified Navier–Stokes equation for blood flow through capillaries) are used for simulating interstitial and intravascular flows and Starling’s law is used for closing this system of equations and coupling the intravascular and extravascular flows. This is the first study of flow modeling in solid tumors to naturalistically couple intravascular and extravascular flow through a network. This network is generated by sprouting angiogenesis and consisting of one parent vessel connected to the network while taking into account the non-continuous behavior of blood, adaptability of capillary diameter to hemodynamics and metabolic stimuli, non-Newtonian blood flow, and phase separation of blood flow in capillary bifurcation. The incorporation of the outlined components beyond the previous models provides a more realistic prediction of interstitial fluid flow pattern in solid tumors and surrounding tissues. Results predict higher interstitial pressure, almost two times, for realistic model compared to the simplified model.     

Evolution of osmotic pressure in solid tumors

The mechanical microenvironment of solid tumors includes both fluid and solid stresses. These stresses play a crucial role in cancer progression and treatment and have been analyzed rigorously both mathematically and experimentally. The magnitude and spatial distribution of osmotic pressures in tumors, however, cannot be measured experimentally and to our knowledge there is no mathematical model to calculate osmotic pressures in the tumor interstitial space. In this study, we developed a triphasic biomechanical model of tumor growth taking into account not only the solid and fluid phase of a tumor, but also the transport of cations and anions, as well as the fixed charges at the surface of the glycosaminoglycan chains. Our model predicts that the osmotic pressure is negligible compared to the interstitial fluid pressure for values of glycosaminoglycans (GAGs) taken from the literature for sarcomas, melanomas and adenocarcinomas. Furthermore, our results suggest that an increase in the hydraulic conductivity of the tumor, increases considerably the intratumoral concentration of free ions and thus, the osmotic pressure but it does not reach the levels of the interstitial fluid pressure.     

Multidrug resistance in lactic acid bacteria: molecular mechanisms and clinical relevance

The active extrusion of cytotoxic compounds from the cell by multidrug transporters is one of the major causes of failure of chemotherapeutic treatment of tumor cells and of infections by pathogenic microorganisms. The secondary multidrug transporter LmrP and the ATP-binding cassette (ABC) type multidrug transporter LmrA in Lactococcus lactis are representatives of the two major classes of multidrug transporters found in pro- and eukaryotic organisms. Therefore, knowledge of the molecular properties of LmrP and LmrA will have a wide significance for multidrug transporters in all living cells, and may enable the development of specific inhibitors and of new drugs which circumvent the action of multidrug transporters. Interestingly, LmrP and LmrA are transport proteins with very different protein structures, which use different mechanisms of energy coupling to transport drugs out of the cell. Surprisingly, both proteins have overlapping specificities for drugs, are inhibited by t he same set of modulators, and transport drugs via a similar transport mechanism. The structure-function relationships that dictate drug recognition and transport by LmrP and LmrA will represent an intriguing new area of research.     

Understanding the role of the tumour vasculature in the transport of drugs to solid cancer tumours

OBJECTIVES: The vasculature of tumours imposes certain barriers that transport of anti-cancer drugs must overcome. Here follows an account of development of a general computational model that describes the mechanisms of drug transport to a solid tumour, with an emphasis on modelling the vasculature using solute transport concepts. MATERIALS AND METHODS: Investigation into the biological parameters that enhance/prevent anticancer drug transport to the tumour provides a means to evaluate the effects of these parameters on the treatment process. Sensitivity analysis of these provides useful insights concerning anticancer drug transport mechanisms from the vasculature to the solid tumour for a non-specified drug and non-specified solid tumour by revealing the conditions that promote or prevent effective drug transport. The effect of the vasculature on transport efficiency is studied using a parametric analysis of some of the transport and biological parameters. Understanding the various transport mechanisms provides a basis to evaluate the effectiveness of the drug treatment a priori. RESULTS: It was found that increases in the capillary hydrostatic pressure, diffusive permeability coefficient and hydraulic conductivity all result in a decrease in tumour size. Similarly, decreases in the interstitium hydrostatic pressure and filtration constant result in a decrease in tumour size. Dependence of the change in the tumour size to changes in these parameters is non-linear. These results demonstrate the potential of the integrated computational model of the tumour and its vasculature to estimate efficacy of a particular treatment process. Regardless of the dependency of the outcome on the assumed model parameters and the assumed kinetics, mathematical models of this type can provide more explanation on the issues related to the transport barriers, the efficacy of the treatment, and the development of effective anticancer drugs. A case study also is presented to demonstrate the model's flexibility to accommodate a two-cell-glioma population.     

Influence of interstitial fluid dynamics on growth and therapy of angiogenic tumor. Analysis by mathematical model

We have developed a spatially distributed mathematical model of angiogenic tumor growth in tissue with account of interstitial fluid dynamics and bevacizumab monotherapy. In this model the process of neovascularization is initiated by tumor cells in a state of metabolic stress, vascular endothelial growth factor (VEGF) being its main mediator. The model takes into consideration the convection flows arising in dense tissue due to active proliferation and migration of tumor cells as well as interstitial fluid inflow from blood vascular system, its outflow through lymphatic system and redistribution in the area of tumor growth. The work considers the diffusive approximation of interstitial fluid dynamics in tumor and normal tissue. Numerical study of the model showed that in absence of therapy a peritumoral edema is formed due to the increase of interstitial fluid inflow from angiogenic capillaries. In the case of rapid interstitial fluid outflow through lymphatic system and its fast transport from necrotic zone to normal tissue the regimes of full growth stop are observed in case of low-invasive tumor. Under bevacizumab monotherapy the peritumoral edema vanishes and low-invasive tumor may not only decelerate its growth, but also start shrinking for a large range of parameters.     

Interstitial Fluid Flow and Drug Delivery in Vascularized Tumors: A Computational Model

Interstitial fluid is a solution that bathes and surrounds the human cells and provides them with nutrients and a way of waste removal. It is generally believed that elevated tumor interstitial fluid pressure (IFP) is partly responsible for the poor penetration and distribution of therapeutic agents in solid tumors, but the complex interplay of extravasation, permeabilities, vascular heterogeneities and diffusive and convective drug transport remains poorly understood. Here we consider–with the help of a theoretical model–the tumor IFP, interstitial fluid flow (IFF) and its impact upon drug delivery within tumor depending on biophysical determinants such as vessel network morphology, permeabilities and diffusive vs. convective transport. We developed a vascular tumor growth model, including vessel co-option, regression, and angiogenesis, that we extend here by the interstitium (represented by a porous medium obeying Darcy''s law) and sources (vessels) and sinks (lymphatics) for IFF. With it we compute the spatial variation of the IFP and IFF and determine its correlation with the vascular network morphology and physiological parameters like vessel wall permeability, tissue conductivity, distribution of lymphatics etc. We find that an increased vascular wall conductivity together with a reduction of lymph function leads to increased tumor IFP, but also that the latter does not necessarily imply a decreased extravasation rate: Generally the IF flow rate is positively correlated with the various conductivities in the system. The IFF field is then used to determine the drug distribution after an injection via a convection diffusion reaction equation for intra- and extracellular concentrations with parameters guided by experimental data for the drug Doxorubicin. We observe that the interplay of convective and diffusive drug transport can lead to quite unexpected effects in the presence of a heterogeneous, compartmentalized vasculature. Finally we discuss various strategies to increase drug exposure time of tumor cells.     

Systematic interaction of plasma albumin with the efficacy of chemotherapeutic drugs

Albumin, as the most abundant plasma protein, plays an integral role in the transport of a variety of exogenous and endogenous ligands in the bloodstream and extravascular spaces. For exogenous drugs, especially chemotherapeutic drugs, binding to and being delivered by albumin can significantly affect their efficacy. Meanwhile, albumin can also bind to many endogenous ligands, such as fatty acids, with important physiological significance that can affect tumor proliferation and metabolism. In this review, we summarize how albumin with unique properties affects chemotherapeutic drugs efficacy from the aspects of drug outcome in blood, toxicity, tumor accumulation and direct or indirect interactions with fatty acids, plus application of albumin-based carriers for anti-tumor drug delivery.     

Influence of interstitial bone microcracks on strain-induced fluid flow

It is well known that microcracks act as a stimulus for bone remodelling, initiating resorption by osteoclasts and new bone formation by osteoblasts. Moreover, microcracks are likely to alter the fluid flow and convective transport through the bone tissue. This paper proposes a quantitative evaluation of the strain-induced interstitial fluid velocities developing in osteons in presence of a microcrack in the interstitial bone tissue. Based on Biot theory in the low-frequency range, a poroelastic model is carried out to study the hydro-mechanical behaviour of cracked osteonal tissue. The finite element results show that the presence of a microcrack in the interstitial osteonal tissue may drastically reduce the fluid velocity inside the neighbouring osteons. This fluid inactive zone inside osteons can cover up to 10% of their surface. Consequently, the fluid environment of bone mechano-sensitive cells is locally modified.     

A new method for isolation of interstitial fluid from human solid tumors applied to proteomic analysis of ovarian carcinoma tissue

Major efforts have been invested in the identification of cancer biomarkers in plasma, but the extraordinary dynamic range in protein composition, and the dilution of disease specific proteins make discovery in plasma challenging. Focus is shifting towards using proximal fluids for biomarker discovery, but methods to verify the isolated sample's origin are missing. We therefore aimed to develop a technique to search for potential candidate proteins in the proximal proteome, i.e. in the tumor interstitial fluid, since the biomarkers are likely to be excreted or derive from the tumor microenvironment. Since tumor interstitial fluid is not readily accessible, we applied a centrifugation method developed in experimental animals and asked whether interstitial fluid from human tissue could be isolated, using ovarian carcinoma as a model. Exposure of extirpated tissue to 106 g enabled tumor fluid isolation. The fluid was verified as interstitial by an isolated fluid:plasma ratio not significantly different from 1.0 for both creatinine and Na(+), two substances predominantly present in interstitial fluid. The isolated fluid had a colloid osmotic pressure 79% of that in plasma, suggesting that there was some sieving of proteins at the capillary wall. Using a proteomic approach we detected 769 proteins in the isolated interstitial fluid, sixfold higher than in patient plasma. We conclude that the isolated fluid represents undiluted interstitial fluid and thus a subproteome with high concentration of locally secreted proteins that may be detected in plasma for diagnostic, therapeutic and prognostic monitoring by targeted methods.     

Numerical modelling of mass transport in an arterial wall with anisotropic transport properties

Coronary artery disease results in blockages or narrowing of the artery lumen. Drug eluting stents (DES) were developed to replace bare metal stents in an effort to combat re-blocking of the diseased artery following treatment. The numerical models developed within this study focus on representing the changing trends of drug delivery from an idealised DES in an arterial wall with an anisotropic ultra-structure. To reduce the computational burden of solving coupled physics problems, a model reduction strategy was adopted. Particular focus has been placed upon adequately modelling the influence of strut compression as there is a paucity of numerical studies that account for changes in transport properties in compressed regions of the arterial wall due to stent deployment. This study developed an idealised numerical modelling framework to account for the changes in the directionally dependent porosity and tortuosities of the arterial wall as a result of radial strut compression. The results show that depending on the degree of strut compression, trends in therapeutic drug delivery within the arterial wall can be either increased or decreased. The study highlights the importance of incorporating compression into numerical models to better represent transport within the arterial wall and suggests an appropriate numerical modelling framework that could be utilised in more realistic patient specific arterial geometries.