Choice of boundary condition for lattice-Boltzmann simulation of moderate-Reynolds-number flow in complex domains

Uncertainty quantification of the lattice Boltzmann method focussing on studies of human-scale vascular blood flow

Uncertainty quantification is becoming a key tool to ensure that numerical models can be sufficiently trusted to be used in domains such as medical device design. Demonstration of how input parameters impact the quantities of interest generated by any numerical model is essential to understanding the limits of its reliability. With the lattice Boltzmann method now a widely used approach for computational fluid dynamics, building greater understanding of its numerical uncertainty characteristics will support its further use in science and industry. In this study we apply an in-depth uncertainty quantification study of the lattice Boltzmann method in a canonical bifurcating geometry that is representative of the vascular junctions present in arterial and venous domains. These campaigns examine how quantities of interest-pressure and velocity along the central axes of the bifurcation-are influenced by the algorithmic parameters of the lattice Boltzmann method and the parameters controlling the values imposed at inlet velocity and outlet pressure boundary conditions. We also conduct a similar campaign on a set of personalised vessels to further illustrate the application of these techniques. Our work provides insights into how input parameters and boundary conditions impact the velocity and pressure distributions calculated in a simulation and can guide the choices of such values when applied to vascular studies of patient specific geometries. We observe that, from an algorithmic perspective, the number of time steps and the size of the grid spacing are the most influential parameters. When considering the influence of boundary conditions, we note that the magnitude of the inlet velocity and the mean pressure applied within sinusoidal pressure outlets have the greatest impact on output quantities of interest. We also observe that, when comparing the magnitude of variation imposed in the input parameters with that observed in the output quantities, this variability is particularly magnified when the input velocity is altered. This study also demonstrates how open-source toolkits for validation, verification and uncertainty quantification can be applied to numerical models deployed on high-performance computers without the need for modifying the simulation code itself. Such an ability is key to the more widespread adoption of the analysis of uncertainty in numerical models by significantly reducing the complexity of their execution and analysis.

Uncertainty quantification of the impact of peripheral arterial disease on abdominal aortic aneurysms in blood flow simulations

Peripheral arterial disease (PAD) and abdominal aortic aneurysms (AAAs) often coexist and pose significant risks of mortality, yet their mutual interactions remain largely unexplored. Here, we introduce a fluid mechanics model designed to simulate the haemodynamic impact of PAD on AAA-associated risk factors. Our focus lies on quantifying the uncertainty inherent in controlling the flow rates within PAD-affected vessels and predicting AAA risk factors derived from wall shear stress. We perform a sensitivity analysis on nine critical model parameters through simulations of three-dimensional blood flow within a comprehensive arterial geometry. Our results show effective control of the flow rates using two-element Windkessel models, although specific outlets need attention. Quantities of interest like endothelial cell activation potential (ECAP) and relative residence time are instructive for identifying high-risk regions, with ECAP showing greater reliability and adaptability. Our analysis reveals that the uncertainty in the quantities of interest is 187% of that of the input parameters. Notably, parameters governing the amplitude and frequency of the inlet velocity exert the strongest influence on the risk factors' variability and warrant precise determination. This study forms the foundation for patient-specific simulations involving PAD and AAAs which should ultimately improve patient outcomes and reduce associated mortality rates.

High resolution simulation of basilar artery infarct and flow within the circle of Willis

On a global scale, cerebro- and cardiovascular diseases have long been one of the leading causes of death and disability and their prevalence appears to be increasing in recent times. Understanding potential biomarkers and risk factors will help to identify individuals potentially at risk of suffering an ischemic stroke. However, the widely variable construction of the cerebral vasculature makes it difficult to provide a specific assessment without the knowledge of a patient's physiology. In this paper we use the 3D blood flow simulator HemeLB to study flow within three common structural variations of the circle of Willis during and in the moments after a blockage of the basilar artery. This tool, based on the lattice Boltzmann method, allows the 3D flow entering the basilar artery to be finely controlled to replicate the cessation of blood feeding this particular vessel-we demonstrate this with several examples including a sudden halt to flow and a gradual loss of flow over three heartbeat cycles. In this work we start with an individualised 3D representation of a full circle of Willis and then construct two further domains by removing the left or right posterior communicating arteries from this geometry. Our results indicate how, and how quickly, the circle of Willis is able to redistribute flow following such a stroke. Due to the choice of infarct, the greatest reduction in flow was observed in the posterior cerebral arteries where flow was reduced by up to 70% in some cases. The high resolution domains used in this study permit the velocity magnitude and wall shear stress to be analysed at key points during and following the stroke. The model we present here indicates how personalised vessels are required to provide the best insight into stroke risk for a given individual.

Parametric analysis of an efficient boundary condition to control outlet flow rates in large arterial networks

Substantial effort is being invested in the creation of a virtual human-a model which will improve our understanding of human physiology and diseases and assist clinicians in the design of personalised medical treatments. A central challenge of achieving blood flow simulations at full-human scale is the development of an efficient and accurate approach to imposing boundary conditions on many outlets. A previous study proposed an efficient method for implementing the two-element Windkessel model to control the flow rate ratios at outlets. Here we clarify the general role of the resistance and capacitance in this approach and conduct a parametric sweep to examine how to choose their values for complex geometries. We show that the error of the flow rate ratios decreases exponentially as the resistance increases. The errors fall below 4% in a simple five-outlets model and 7% in a human artery model comprising ten outlets. Moreover, the flow rate ratios converge faster and suffer from weaker fluctuations as the capacitance decreases. Our findings also establish constraints on the parameters controlling the numerical stability of the simulations. The findings from this work are directly applicable to larger and more complex vascular domains encountered at full-human scale.

An efficient, localised approach for the simulation of elastic blood vessels using the lattice Boltzmann method

Many numerical studies of blood flow impose a rigid wall assumption due to the simplicity of its implementation compared to a full coupling with a solid mechanics model. In this paper, we present a localised method for incorporating the effects of elastic walls into blood flow simulations using the lattice Boltzmann method implemented by the open-source code HemeLB. We demonstrate that our approach is able to more accurately capture the flow behaviour expected in elastic walled vessels than ones with rigid walls. Furthermore, we show that this can be achieved with no loss of computational performance and remains strongly scalable on high performance computers. We finally illustrate that our approach captures the same trends in wall shear stress distribution as those observed in studies using a rigorous coupling between fluid dynamics and solid mechanics models to solve flow in personalised vascular geometries. These results demonstrate that our model can be used to efficiently and effectively represent flows in elastic blood vessels.

High fidelity blood flow in a patient-specific arteriovenous fistula

An arteriovenous fistula, created by artificially connecting segments of a patient’s vasculature, is the preferred way to gain access to the bloodstream for kidney dialysis. The increasing power and availability of supercomputing infrastructure means that it is becoming more realistic to use simulations to help identify the best type and location of a fistula for a specific patient. We describe a 3D fistula model that uses the lattice Boltzmann method to simultaneously resolve blood flow in patient-specific arteries and veins. The simulations conducted here, comprising vasculatures of the whole forearm, demonstrate qualified validation against clinical data. Ongoing research to further encompass complex biophysics on realistic time scales will permit the use of human-scale physiological models for basic and clinical medicine.

Enhanced single-node lattice Boltzmann boundary condition for fluid flows

We propose a procedure to implement Dirichlet velocity boundary conditions for complex shapes that use data from a single node only, in the context of the lattice Boltzmann method. Two ideas are at the base of this approach. The first is to generalize the geometrical description of boundary conditions combining bounce-back rule with interpolations. The second is to enhance them by limiting the interpolation extension to the proximity of the boundary. Despite its local nature, the resulting method exhibits second-order convergence for the velocity field and shows similar or better accuracy than the well-established Bouzidi's scheme for curved walls [M. Bouzidi, M. Firdaouss, and P. Lallemand, Phys. Fluids 13, 3452 (2001)]PHFLE61070-663110.1063/1.1399290. Among the infinite number of possibilities, we identify several meaningful variants of the method, discerned by their approximation of the second-order nonequilibrium terms and their interpolation coefficients. For each one, we provide two parametrized versions that produce viscosity independent accuracy at steady state. The method proves to be suitable to simulate moving rigid objects or surfaces moving following either the rigid body dynamics or a prescribed kinematic. Also, it applies uniformly and without modifications in the whole domain for any shape, including corners, narrow gaps, or any other singular geometry.

Towards blood flow in the virtual human: efficient self-coupling of HemeLB

Many scientific and medical researchers are working towards the creation of a virtual human—a personalized digital copy of an individual—that will assist in a patient’s diagnosis, treatment and recovery. The complex nature of living systems means that the development of this remains a major challenge. We describe progress in enabling the HemeLB lattice Boltzmann code to simulate 3D macroscopic blood flow on a full human scale. Significant developments in memory management and load balancing allow near linear scaling performance of the code on hundreds of thousands of computer cores. Integral to the construction of a virtual human, we also outline the implementation of a self-coupling strategy for HemeLB. This allows simultaneous simulation of arterial and venous vascular trees based on human-specific geometries.

Consistent vortex initialization for the athermal lattice Boltzmann method

A barotropic counterpart of the well-known convected vortex test case is rigorously derived from the Euler equations along with an athermal equation of state. Starting from a given velocity distribution corresponding to an intended flow recirculation, the athermal counterpart of the Euler equations are solved to obtain a consistent density field. The present initialization is assessed on a standard lattice Boltzmann solver based on the D2Q9 lattice. Compared to the usual isentropic initialization, a much lower spurious relaxation toward the targeted solution is observed, which is due to the spatial resolution rather than approximated macroscopic quantities. The amplitude of the spurious waves can be further reduced by including an off-equilibrium part in the initial distribution functions.

A voxel image-based pulmonary airflow simulation method with an automatic detection algorithm for airway outlets

Investigations of pulmonary airflows in respiratory systems are important for the diagnostics and treatment of pulmonary diseases. For accurate prediction of the flow field in an airway, a numerical simulation must be conducted using the true geometry from computed tomography (CT) data. Flow simulation is still a difficult task because of the mesh generation process and preprocessing setup procedures. In this study, we developed a voxel image-based simulation method using an automatic detection algorithm for airway outlets to simplify the simulation process and improve its applicability in the medical field. Our approach is based on the lattice Boltzmann method with a topology analysis algorithm, which can preserve all raw information from the original CT images and give an accurate flow field inside the airways. Our method can reproduce the essential flow features inside airways, is highly efficient, and decreases the overall processing time. Thus, it has a great potential for future real-time airflow analyses to provide airflow information to medical experts. HIGHLIGHTS: This paper proposed a voxel image-based simulation method with a novel automatic outlet-selecting algorithm to calculate the velocity and pressure of physiological flows in multi-generation-branched airways. Our approach simplifies the simulation process by automatically applying the boundary conditions to large numbers of outlets and minimizes the time-consuming mesh generation process. Our proposed method has considerable potential for real-time simulations improving the applicability to patient-specific medical diagnostics and treatments.

Suitability of lattice Boltzmann inlet and outlet boundary conditions for simulating flow in image-derived vasculature

The lattice Boltzmann method (LBM) is a popular alternative to solving the Navier-Stokes equations for modeling blood flow. When simulating flow using the LBM, several choices for inlet and outlet boundary conditions exist. While boundary conditions in the LBM have been evaluated in idealized geometries, there have been no extensive comparisons in image-derived vasculature, where the geometries are highly complex. In this study, the Zou-He (ZH) and finite difference (FD) boundary conditions were evaluated in image-derived vascular geometries by comparing their stability, accuracy, and run times. The boundary conditions were compared in four arteries: a coarctation of the aorta, dissected aorta, femoral artery, and left coronary artery. The FD boundary condition was more stable than ZH in all four geometries. In general, simulations using the ZH and FD method showed similar convergence rates within each geometry. However, the ZH method proved to be slightly more accurate compared with experimental flow using three-dimensional printed vasculature. The total run times necessary for simulations using the ZH boundary condition were significantly higher as the ZH method required a larger relaxation time, grid resolution, and number of time steps for a simulation representing the same physiological time. Finally, a new inlet velocity profile algorithm is presented for complex inlet geometries. Overall, results indicated that the FD method should generally be used for large-scale blood flow simulations in image-derived vasculature geometries. This study can serve as a guide to researchers interested in using the LBM to simulate blood flow.

Choice of no-slip curved boundary condition for lattice Boltzmann simulations of high-Reynolds-number flows

Various curved no-slip boundary conditions available in literature improve the accuracy of lattice Boltzmann simulations compared to the traditional staircase approximation of curved geometries. Usually, the required unknown distribution functions emerging from the solid nodes are computed based on the known distribution functions using interpolation or extrapolation schemes. On using such curved boundary schemes, there will be mass loss or gain at each time step during the simulations, especially apparent at high Reynolds numbers, which is called mass leakage. Such an issue becomes severe in periodic flows, where the mass leakage accumulation would affect the computed flow fields over time. In this paper, we examine mass leakage of the most well-known curved boundary treatments for high-Reynolds-number flows. Apart from the existing schemes, we also test different forced mass conservation schemes and a constant density scheme. The capability of each scheme is investigated and, finally, recommendations for choosing a proper boundary condition scheme are given for stable and accurate simulations.

Validation of Patient-Specific Cerebral Blood Flow Simulation Using Transcranial Doppler Measurements

We present a validation study comparing results from a patient-specific lattice-Boltzmann simulation to transcranial Doppler (TCD) velocity measurements in four different planes of the middle cerebral artery (MCA). As part of the study, we compared simulations using a Newtonian and a Carreau-Yasuda rheology model. We also investigated the viability of using downscaled velocities to reduce the required resolution. Simulations with unscaled velocities predict the maximum flow velocity with an error of less than 9%, independent of the rheology model chosen. The accuracy of the simulation predictions worsens considerably when simulations are run at reduced velocity, as is for example the case when inflow velocities from healthy individuals are used on a vascular model of a stroke patient. Our results demonstrate the importance of using directly measured and patient-specific inflow velocities when simulating blood flow in MCAs. We conclude that localized TCD measurements together with predictive simulations can be used to obtain flow estimates with high fidelity over a larger region, and reduce the need for more invasive flow measurement procedures.

Clinical validation and assessment of aortic hemodynamics using computational fluid dynamics simulations from computed tomography angiography

Background

Hemodynamic information including peak systolic pressure (PSP) and peak systolic velocity (PSV) carry an important role in evaluation and diagnosis of congenital heart disease (CHD). Since MDCTA cannot evaluate hemodynamic information directly, the aim of this study is to provide a noninvasive method based on a computational fluid dynamics (CFD) model, derived from multi-detector computed tomography angiography (MDCTA) raw data, to analyze the aortic hemodynamics in infants with CHD, and validate these results against echocardiography and cardiac catheter measurements.

Methods

This study included 25 patients (17 males, and 8 females; a median age of 2 years, range: 4 months–4 years) with CHD. All patients underwent both transthoracic echocardiography (TTE) and MDCTA within 2 weeks prior to cardiac catheterization. CFD models were created from MDCTA raw data. Boundary conditions were confirmed by lumped parameter model and transthoracic echocardiography (TTE). Peak systolic velocity derived from CFD models (PSV_CFD) was compared to TTE measurements (PSV_TTE), while the peak systolic pressure derived from CFD (PSP_CFD) was compared to catheterization (PSP_CC). Regions with low and high peak systolic wall shear stress (PSWSS) were also evaluated.

Results

PSV_CFD and PSP_CFD showed good agreements between PSV_TTE (r = 0.968, p < 0.001; mean bias = − 7.68 cm/s) and PSP_CC (r = 0.918, p < 0.001; mean bias = 1.405 mmHg). Regions with low and high PSWSS) can also be visualized. Skewing of velocity or helical blood flow was also observed at aortic arch in patients.

Conclusions

Our result demonstrated that CFD scheme based on MDCTA raw data is an accurate and convenient method in obtaining the velocity and pressure from aorta and displaying the distribution of PSWSS and flow pattern of aorta. The preliminary results from our study demonstrate the capability in combining clinical imaging data and novel CFD tools in infants with CHD and provide a noninvasive approach for diagnose of CHD such as coarctation of aorta in future.

Modeling Patient-Specific Magnetic Drug Targeting Within the Intracranial Vasculature

Drug targeting promises to substantially enhance future therapies, for example through the focussing of chemotherapeutic drugs at the site of a tumor, thus reducing the exposure of healthy tissue to unwanted damage. Promising work on the steering of medication in the human body employs magnetic fields acting on nanoparticles made of paramagnetic materials. We develop a computational tool to aid in the optimization of the physical parameters of these particles and the magnetic configuration, estimating the fraction of particles reaching a given target site in a large patient-specific vascular system for different physiological states (heart rate, cardiac output, etc.). We demonstrate the excellent computational performance of our model by its application to the simulation of paramagnetic-nanoparticle-laden flows in a circle of Willis geometry obtained from an MRI scan. The results suggest a strong dependence of the particle density at the target site on the strength of the magnetic forcing and the velocity of the background fluid flow.

Transitional hemodynamics in intracranial aneurysms - Comparative velocity investigations with high resolution lattice Boltzmann simulations, normal resolution ANSYS simulations, and MR imaging

PURPOSE

Blood flow in intracranial aneurysms has, until recently, been considered to be disturbed but still laminar. Recent high resolution computational studies have demonstrated, in some situations, however, that the flow may exhibit high frequency fluctuations that resemble weakly turbulent or transitional flow. Due to numerous assumptions required for simplification in computational fluid dynamics (CFD) studies, the occurrence of these events, in vivo, remains unsettled. The detection of these fluctuations in aneurysmal blood flow, i.e., hemodynamics by CFD, poses additional challenges as such phenomena cannot be captured in clinical data acquisition with magnetic resonance (MR) due to inadequate temporal and spatial resolutions. The authors' purpose was to address this issue by comparing results from highly resolved simulations, conventional resolution laminar simulations, and MR measurements, identify the differences, and identify their causes.

METHODS

Two aneurysms in the basilar artery, one with disturbed yet laminar flow and the other with transitional flow, were chosen. One set of highly resolved direct numerical simulations using the lattice Boltzmann method (LBM) and another with adequate resolutions under laminar flow assumption were conducted using a commercially available ANSYS Fluent solver. The velocity fields obtained from simulation results were qualitatively and statistically compared against each other and with MR acquisition.

RESULTS

Results from LBM, ANSYS Fluent, and MR agree well qualitatively and quantitatively for one of the aneurysms with laminar flow in which fluctuations were <80 Hz. The comparisons for the second aneurysm with high fluctuations of > ∼ 600 Hz showed vivid differences between LBM, ANSYS Fluent, and magnetic resonance imaging. After ensemble averaging and down-sampling to coarser space and time scales, these differences became minimal.

CONCLUSIONS

A combination of MR derived data and CFD can be helpful in estimating the hemodynamic environment of intracranial aneurysms. Adequately resolved CFD would suffice gross assessment of hemodynamics, potentially in a clinical setting, and highly resolved CFD could be helpful in a detailed and retrospective understanding of the physiological mechanisms.

Prediction of the moments in advection-diffusion lattice Boltzmann method. II. Attenuation of the boundary layers via double-Λ bounce-back flux scheme

Impact of the unphysical tangential advective-diffusion constraint of the bounce-back (BB) reflection on the impermeable solid surface is examined for the first four moments of concentration. Despite the number of recent improvements for the Neumann condition in the lattice Boltzmann method-advection-diffusion equation, the BB rule remains the only known local mass-conserving no-flux condition suitable for staircase porous geometry. We examine the closure relation of the BB rule in straight channel and cylindrical capillary analytically, and show that it excites the Knudsen-type boundary layers in the nonequilibrium solution for full-weight equilibrium stencil. Although the d2Q5 and d3Q7 coordinate schemes are sufficient for the modeling of isotropic diffusion, the full-weight stencils are appealing for their advanced stability, isotropy, anisotropy and anti-numerical-diffusion ability. The boundary layers are not covered by the Chapman-Enskog expansion around the expected equilibrium, but they accommodate the Chapman-Enskog expansion in the bulk with the closure relation of the bounce-back rule. We show that the induced boundary layers introduce first-order errors in two primary transport properties, namely, mean velocity (first moment) and molecular diffusion coefficient (second moment). As a side effect, the Taylor-dispersion coefficient (second moment), skewness (third moment), and kurtosis (fourth moment) deviate from their physical values and predictions of the fourth-order Chapman-Enskog analysis, even though the kurtosis error in pure diffusion does not depend on grid resolution. In two- and three-dimensional grid-aligned channels and open-tubular conduits, the errors of velocity and diffusion are proportional to the diagonal weight values of the corresponding equilibrium terms. The d2Q5 and d3Q7 schemes do not suffer from this deficiency in grid-aligned geometries but they cannot avoid it if the boundaries are not parallel to the coordinate lines. In order to vanish or attenuate the disparity of the modeled transport coefficients with the equilibrium weights without any modification of the BB rule, we propose to use the two-relaxation-times collision operator with free-tunable product of two eigenfunctions Λ. Two different values Λ_{v} and Λ_{b} are assigned for bulk and boundary nodes, respectively. The rationale behind this is that Λ_{v} is adjustable for stability, accuracy, or other purposes, while the corresponding Λ_{b}(Λ_{v}) controls the primary accommodation effects. Two distinguished but similar functional relations Λ_{b}(Λ_{v}) are constructed analytically: they preserve advection velocity in parabolic profile, exactly in the two-dimensional channel and very accurately in a three-dimensional cylindrical capillary. For any velocity-weight stencil, the (local) double-Λ BB scheme produces quasi-identical solutions with the (nonlocal) specular-forward reflection for first four moments in a channel. In a capillary, this strategy allows for the accurate modeling of the Taylor-dispersion and non-Gaussian effects. As illustrative example, it is shown that in the flow around a circular obstacle, the double-Λ scheme may also vanish the dependency of mean velocity on the velocity weight; the required value for Λ_{b}(Λ_{v}) can be identified in a few bisection iterations in given geometry. A positive solution for Λ_{b}(Λ_{v}) may not exist in pure diffusion, but a sufficiently small value of Λ_{b} significantly reduces the disparity in diffusion coefficient with the mass weight in ducts and in the presence of rectangular obstacles. Although Λ_{b} also controls the effective position of straight or curved boundaries, the double-Λ scheme deals with the lower-order effects. Its idea and construction may help understanding and amelioration of the anomalous, zero- and first-order behavior of the macroscopic solution in the presence of the bulk and boundary or interface discontinuities, commonly found in multiphase flow and heterogeneous transport.

Computational fluid dynamics assisted characterization of parafoveal hemodynamics in normal and diabetic eyes using adaptive optics scanning laser ophthalmoscopy

Diabetic retinopathy (DR) is the leading cause of visual loss in working-age adults worldwide. Previous studies have found hemodynamic changes in the diabetic eyes, which precede clinically evident pathological alterations of the retinal microvasculature. There is a pressing need for new methods to allow greater understanding of these early hemodynamic changes that occur in DR. In this study, we propose a noninvasive method for the assessment of hemodynamics around the fovea (a region of the eye of paramount importance for vision). The proposed methodology combines adaptive optics scanning laser ophthalmoscopy and computational fluid dynamics modeling. We compare results obtained with this technique with in vivo measurements of blood flow based on blood cell aggregation tracking. Our results suggest that parafoveal hemodynamics, such as capillary velocity, wall shear stress, and capillary perfusion pressure can be noninvasively and reliably characterized with this method in both healthy and diabetic retinopathy patients.

MRI-based computational hemodynamics in patients with aortic coarctation using the lattice Boltzmann methods: Clinical validation study

Purpose

To introduce a scheme based on a recent technique in computational hemodynamics, known as the lattice Boltzmann methods (LBM), to noninvasively measure pressure gradients in patients with a coarctation of the aorta (CoA). To provide evidence on the accuracy of the proposed scheme, the computed pressure drop values are compared against those obtained using the reference standard method of catheterization.

Materials and Methods

Pre‐ and posttreatment LBM‐based pressure gradients for 12 patients with CoA were simulated for the time point of peak systole using the open source library OpenLB. Four‐dimensional (4D) flow‐sensitive phase‐contrast MRI at 1.5 Tesla was used to acquire flow and to setup the simulation. The vascular geometry was reconstructed using 3D whole‐heart MRI. Patients underwent pre‐ and postinterventional pressure catheterization as a reference standard.

Results

There is a significant linear correlation between the pretreatment catheter pressure drops and those computed based on the LBM simulation, r=.85, P<.001. The bias was ‐0.58 ± 4.1 mmHg and was not significant ( P=0.64) with a 95% confidence interval (CI) of ‐3.22 to 2.06. For the posttreatment results, the bias was larger and at ‐2.54 ± 3.53 mmHg with a 95% CI of ‐0.17 to ‐4.91 mmHg.

Conclusion

The results indicate a reasonable agreement between the simulation results and the catheter measurements. LBM‐based computational hemodynamics can be considered as an alternative to more traditional computational fluid dynamics schemes for noninvasive pressure calculations and can assist in diagnosis and therapy planning.

Level of Evidence: 3

J. Magn. Reson. Imaging 2017;45:139–146.

Computer simulations reveal complex distribution of haemodynamic forces in a mouse retina model of angiogenesis

There is currently limited understanding of the role played by haemodynamic forces on the processes governing vascular development. One of many obstacles to be overcome is being able to measure those forces, at the required resolution level, on vessels only a few micrometres thick. In this paper, we present an in silico method for the computation of the haemodynamic forces experienced by murine retinal vasculature (a widely used vascular development animal model) beyond what is measurable experimentally. Our results show that it is possible to reconstruct high-resolution three-dimensional geometrical models directly from samples of retinal vasculature and that the lattice-Boltzmann algorithm can be used to obtain accurate estimates of the haemodynamics in these domains. We generate flow models from samples obtained at postnatal days (P) 5 and 6. Our simulations show important differences between the flow patterns recovered in both cases, including observations of regression occurring in areas where wall shear stress (WSS) gradients exist. We propose two possible mechanisms to account for the observed increase in velocity and WSS between P5 and P6: (i) the measured reduction in typical vessel diameter between both time points and (ii) the reduction in network density triggered by the pruning process. The methodology developed herein is applicable to other biomedical domains where microvasculature can be imaged but experimental flow measurements are unavailable or difficult to obtain.

Flexible composition and execution of high performance, high fidelity multiscale biomedical simulations

Multiscale simulations are essential in the biomedical domain to accurately model human physiology. We present a modular approach for designing, constructing and executing multiscale simulations on a wide range of resources, from laptops to petascale supercomputers, including combinations of these. Our work features two multiscale applications, in-stent restenosis and cerebrovascular bloodflow, which combine multiple existing single-scale applications to create a multiscale simulation. These applications can be efficiently coupled, deployed and executed on computers up to the largest (peta) scale, incurring a coupling overhead of 1-10% of the total execution time.

Impact of blood rheology on wall shear stress in a model of the middle cerebral artery

Perturbations to the homeostatic distribution of mechanical forces exerted by blood on the endothelial layer have been correlated with vascular pathologies, including intracranial aneurysms and atherosclerosis. Recent computational work suggests that, in order to correctly characterize such forces, the shear-thinning properties of blood must be taken into account. To the best of our knowledge, these findings have never been compared against experimentally observed pathological thresholds. In this work, we apply the three-band diagram (TBD) analysis due to Gizzi et al. (Gizzi et al. 2011 Three-band decomposition analysis of wall shear stress in pulsatile flows. Phys. Rev. E 83, 031902. (doi:10.1103/PhysRevE.83.031902)) to assess the impact of the choice of blood rheology model on a computational model of the right middle cerebral artery. Our results show that, in the model under study, the differences between the wall shear stress predicted by a Newtonian model and the well-known Carreau-Yasuda generalized Newtonian model are only significant if the vascular pathology under study is associated with a pathological threshold in the range 0.94-1.56 Pa, where the results of the TBD analysis of the rheology models considered differs. Otherwise, we observe no significant differences.