Non-Newtonian Rheology in Blood Circulation

Blood is a complex suspension that demonstrates several non-Newtonian rheological characteristics such as deformation-rate dependency, viscoelasticity and yield stress. In this paper we outline some issues related to the non-Newtonian effects in blood circulation system and present modeling approaches based mostly on the past work in this field.Comment: 26 pages, 5 figures, 2 table

A mathematical model of Doxorubicin treatment efficacy on non-Hodgkin’s lymphoma: Investigation of current protocol through theoretical modelling results

Doxorubicin treatment outcomes for non-Hodgkin’s lymphomas (NHL) are mathematically modelled and computationally analyzed. The NHL model includes a tumor structure incorporating mature and immature vessels, vascular structural adaptation and NHL cell-cycle kinetics in addition to Doxorubicin pharmacokinetics (PK) and pharmacodynamics (PD). Simulations provide qualitative estimations of the effect of Doxorubicin on high-grade (HG), intermediate-grade (IG) and low-grade (LG) NHL. Simulation results imply that if the interval between successive drug applications is prolonged beyond a certain point, treatment will be inefficient due to effects caused by heterogeneous blood flow in the system

Computational Simulations for Aortic Coarctation: Representative Results From a Sampling of Patients

Treatments for coarctation of the aorta (CoA) can alleviate blood pressure (BP) gradients(D), but long-term morbidity still exists that can be explained by altered indices of hemodynamics and biomechanics. We introduce a technique to increase our understanding of these indices for CoA under resting and nonresting conditions, quantify their contribution to morbidity, and evaluate treatment options. Patient-specific computational fluid dynamics (CFD) models were created from imaging and BP data for one normal and four CoA patients (moderate native CoA: D12 mmHg, severe native CoA: D25 mmHg and postoperative end-to-end and end-to-side patients: D0 mmHg). Simulations incorporated vessel deformation, downstream vascular resistance and compliance. Indices including cyclic strain, time-averaged wall shear stress (TAWSS), and oscillatory shear index (OSI) were quantified. Simulations replicated resting BP and blood flow data. BP during simulated exercise for the normal patient matched reported values. Greatest exercise-induced increases in systolic BP and mean and peak DBP occurred for the moderate native CoA patient (SBP: 115 to 154 mmHg; mean and peak DBP: 31 and 73 mmHg). Cyclic strain was elevated proximal to the coarctation for native CoA patients, but reduced throughout the aorta after treatment. A greater percentage of vessels was exposed to subnormal TAWSS or elevated OSI for CoA patients. Local patterns of these indices reported to correlate with atherosclerosis in normal patients were accentuated by CoA. These results apply CFD to a range of CoA patients for the first time and provide the foundation for future progress in this area

Disambiguating the role of blood flow and global signal with partial information decomposition

Global signal (GS) is an ubiquitous construct in resting state functional magnetic resonance imaging (rs-fMRI), associated to nuisance, but containing by definition most of the neuronal signal. Global signal regression (GSR) effectively removes the impact of physiological noise and other artifacts, but at the same time it alters correlational patterns in unpredicted ways. Performing GSR taking into account the underlying physiology (mainly the blood arrival time) has been proven to be beneficial. From these observations we aimed to: 1) characterize the effect of GSR on network-level functional connectivity in a large dataset; 2) assess the complementary role of global signal and vessels; and 3) use the framework of partial information decomposition to further look into the joint dynamics of the global signal and vessels, and their respective influence on the dynamics of cortical areas. We observe that GSR affects intrinsic connectivity networks in the connectome in a non-uniform way. Furthermore, by estimating the predictive information of blood flow and the global signal using partial information decomposition, we observe that both signals are present in different amounts across intrinsic connectivity networks. Simulations showed that differences in blood arrival time can largely explain this phenomenon, while using hemodynamic and calcium mouse recordings we were able to confirm the presence of vascular effects, as calcium recordings lack hemodynamic information. With these results we confirm network-specific effects of GSR and the importance of taking blood flow into account for improving de-noising methods. Additionally, and beyond the mere issue of data denoising, we quantify the diverse and complementary effect of global and vessel BOLD signals on the dynamics of cortical areas

Outflow boundary conditions for 3D simulations of non-periodic blood flow and pressure fields in deformable arteries

The simulation of blood flow and pressure in arteries requires outflow boundary conditions that incorporate models of downstream domains. We previously described a coupled multidomain method to couple analytical models of the downstream domains with 3D numerical models of the upstream vasculature. This prior work either included pure resistance boundary conditions or impedance boundary conditions based on assumed periodicity of the solution. However, flow and pressure in arteries are not necessarily periodic in time due to heart rate variability, respiration, complex transitional flow or acute physiological changes. We present herein an approach for prescribing lumped parameter outflow boundary conditions that accommodate transient phenomena. We have applied this method to compute haemodynamic quantities in different physiologically relevant cardiovascular models, including patient-specific examples, to study non-periodic flow phenomena often observed in normal subjects and in patients with acquired or congenital cardiovascular disease. The relevance of using boundary conditions that accommodate transient phenomena compared with boundary conditions that assume periodicity of the solution is discussed

Structural Adaptation and Heterogeneity of Normal and Tumor Microvascular Networks

Relative to normal tissues, tumor microcirculation exhibits high structural and functional heterogeneity leading to hypoxic regions and impairing treatment efficacy. Here, computational simulations of blood vessel structural adaptation are used to explore the hypothesis that abnormal adaptive responses to local hemodynamic and metabolic stimuli contribute to aberrant morphological and hemodynamic characteristics of tumor microcirculation. Topology, vascular diameter, length, and red blood cell velocity of normal mesenteric and tumor vascular networks were recorded by intravital microscopy. Computational models were used to estimate hemodynamics and oxygen distribution and to simulate vascular diameter adaptation in response to hemodynamic, metabolic and conducted stimuli. The assumed sensitivity to hemodynamic and conducted signals, the vascular growth tendency, and the random variability of vascular responses were altered to simulate ‘normal’ and ‘tumor’ adaptation modes. The heterogeneous properties of vascular networks were characterized by diameter mismatch at vascular branch points (d3var) and deficit of oxygen delivery relative to demand (O2def). In the tumor, d3var and O2def were higher (0.404 and 0.182) than in normal networks (0.278 and 0.099). Simulated remodeling of the tumor network with ‘normal’ parameters gave low values (0.288 and 0.099). Conversely, normal networks attained tumor-like characteristics (0.41 and 0.179) upon adaptation with ‘tumor’ parameters, including low conducted sensitivity, increased growth tendency, and elevated random biological variability. It is concluded that the deviant properties of tumor microcirculation may result largely from defective structural adaptation, including strongly reduced responses to conducted stimuli