Three-dimensional hemodynamics analysis of the circle of Willis in the patient-specific nonintegral arterial structures

Multi-Dimensional Modeling of Cerebral Hemodynamics: A Systematic Review

The Circle of Willis (CoW) describes the arterial system in the human brain enabling the neurovascular blood supply. Neurovascular diseases like intracranial aneurysms (IAs) can occur within the CoW and carry the risk of rupture, which can lead to subarachnoid hemorrhage. The assessment of hemodynamic information in these pathologies is crucial for their understanding regarding detection, diagnosis and treatment. Multi-dimensional in silico approaches exist to evaluate these hemodynamics based on patient-specific input data. The approaches comprise low-scale (zero-dimensional, one-dimensional) and high-scale (three-dimensional) models as well as multi-scale coupled models. The input data can be derived from medical imaging, numerical models, literature-based assumptions or from measurements within healthy subjects. Thus, the most realistic description of neurovascular hemodynamics is still controversial. Within this systematic review, first, the models of the three scales (0D, 1D, 3D) and second, the multi-scale models, which are coupled versions of the three scales, were discussed. Current best practices in describing neurovascular hemodynamics most realistically and their clinical applicablility were elucidated. The performance of 3D simulation entails high computational expenses, which could be reduced by analyzing solely the region of interest in detail. Medical imaging to establish patient-specific boundary conditions is usually rare, and thus, lower dimensional models provide a realistic mimicking of the surrounding hemodynamics. Multi-scale coupling, however, is computationally expensive as well, especially when taking all dimensions into account. In conclusion, the 0D-1D-3D multi-scale approach provides the most realistic outcome; nevertheless, it is least applicable. A 1D-3D multi-scale model can be considered regarding a beneficial trade-off between realistic results and applicable performance.

Effects of abnormal vertebral arteries and the circle of Willis on vertebrobasilar dolichoectasia: A multi-scale simulation study

BACKGROUND

Vertebrobasilar dolichoectasia is a rare cerebrovascular disease characterized by obvious extension, dilation and tortuosity of vertebrobasilar artery, and its pathophysiological mechanism is not clear. This study focused on local hemodynamic changes in basilar arteries with typical vertebrobasilar dolichoectasia, together with unbalanced vertebral arteries and abnormal structures of the circle of Willis, through multi-scale modeling.

METHODS

Three-dimensional models of 3 types of vertebrobasilar arteries were constructed from magnetic resonance images. The first type has no vertebrobasilar dolichoectasia, the second type has vertebrobasilar dolichoectasia and balanced vertebral arteries, and the third type has vertebrobasilar dolichoectasia and unbalanced vertebral arteries. A lumped parameter model of the circle of Willis was established and coupled to these three-dimensional models.

FINDINGS

The results showed that unbalanced bilateral vertebral arteries, especially single vertebral artery deletion mutation, might associate with higher wall shear stress on anterior wall of basilar artery in patients with vertebrobasilar dolichoectasia. And unbalanced bilateral vertebral arteries would increase the blood pressure in basilar artery. Meanwhile, missing communicating arteries in the circle of Willis, especially bilateral posterior communicating arteries absences, would significantly increase blood pressure in basilar artery. The unilateral absence of posterior communicating arteries would increase differences in blood flow between the left and right posterior cerebral arteries.

INTERPRETATION

This study provided a multi-scale modeling method and some preliminary results for helping understand the role of hemodynamics in occurrence and development of vertebrobasilar dolichoectasia.

Determining Clinically-Viable Biomarkers for Ischaemic Stroke Through a Mechanistic and Machine Learning Approach

Assessment of distal cerebral perfusion after ischaemic stroke is currently only possible through expensive and time-consuming imaging procedures which require the injection of a contrast medium. Alternative approaches that could indicate earlier the impact of blood flow occlusion on distal cerebral perfusion are currently lacking. The aim of this study was to identify novel biomarkers suitable for clinical implementation using less invasive diagnostic techniques such as Transcranial Doppler (TCD). We used 1D modelling to simulate pre- and post-stroke velocity and flow wave propagation in a typical arterial network, and Sobol’s sensitivity analysis, supported by the use of Gaussian process emulators, to identify biomarkers linked to cerebral perfusion. We showed that values of pulsatility index of the right anterior cerebral artery > 1.6 are associated with poor perfusion and may require immediate intervention. Three additional biomarkers with similar behaviour, all related to pulsatility indices, were identified. These results suggest that flow pulsatility measured at specific locations could be used to effectively estimate distal cerebral perfusion rates, and ultimately improve clinical diagnosis and management of ischaemic stroke.

Middle cerebral artery pressure laterality in patients with symptomatic ICA stenosis

An internal carotid artery (ICA) stenosis can potentially decrease the perfusion pressure to the brain. In this study, computational fluid dynamics (CFD) was used to study if there was a hemispheric pressure laterality between the contra- and ipsilateral middle cerebral artery (MCA) in patients with a symptomatic ICA stenosis. We further investigated if this MCA pressure laterality (ΔPMCA) was related to the hemispheric flow laterality (ΔQ) in the anterior circulation, i.e., ICA, proximal MCA and the proximal anterior cerebral artery (ACA). Twenty-eight patients (73±6 years, range 59-80 years, 21 men) with symptomatic ICA stenosis were included. Flow rates were measured using 4D flow MRI data (PC-VIPR) and vessel geometries were obtained from computed tomography angiography. The ΔPMCA was calculated from CFD, where patient-specific flow rates were applied at all input- and output boundaries. The ΔPMCA between the contra- and ipsilateral side was 6.4±8.3 mmHg (p<0.001) (median 3.9 mmHg, range -1.3 to 31.9 mmHg). There was a linear correlation between the ΔPMCA and ΔQICA (r = 0.85, p<0.001) and ΔQACA (r = 0.71, p<0.001), respectively. The correlation to ΔQMCA was weaker (r = 0.47, p = 0.011). In conclusion, the MCA pressure laterality obtained with CFD, is a promising physiological biomarker that can grade the hemodynamic disturbance in patients with a symptomatic ICA stenosis.

Validation and Diagnostic Performance of a CFD-Based Non-invasive Method for the Diagnosis of Aortic Coarctation

Purpose: The clinical diagnosis of aorta coarctation (CoA) constitutes a challenge, which is usually tackled by applying the peak systolic pressure gradient (PSPG) method. Recent advances in computational fluid dynamics (CFD) have suggested that multi-detector computed tomography angiography (MDCTA)-based CFD can serve as a non-invasive PSPG measurement. The aim of this study was to validate a new CFD method that does not require any medical examination data other than MDCTA images for the diagnosis of CoA.

Materials and methods: Our study included 65 pediatric patients (38 with CoA, and 27 without CoA). All patients underwent cardiac catheterization to confirm if they were suffering from CoA or any other congenital heart disease (CHD). A series of boundary conditions were specified and the simulated results were combined to obtain a stenosis pressure-flow curve. Subsequently, we built a prediction model and evaluated its predictive performance by considering the AUC of the ROC by 5-fold cross-validation.

Results: The proposed MDCTA-based CFD method exhibited a good predictive performance in both the training and test sets (average AUC: 0.948 vs. 0.958; average accuracies: 0.881 vs. 0.877). It also had a higher predictive accuracy compared with the non-invasive criteria presented in the European Society of Cardiology (ESC) guidelines (average accuracies: 0.877 vs. 0.539).

Conclusion: The new non-invasive CFD-based method presented in this work is a promising approach for the accurate diagnosis of CoA, and will likely benefit clinical decision-making.

Longitudinal computational fluid dynamics study of stenosis and aneurysmal degeneration of an aortorenal bypass

Saphenous vein graft (SVG) bypass placement is regarded as the optimal option for renal artery stenosis, which usually causes secondary hypertension and poor renal perfusion. Using computational fluid dynamics, this study aimed to investigate the underlying hemodynamic mechanism of the vein aneurysm and stenosis after aortorenal bypass surgery. Three-dimensional models were reconstructed based on computed tomographic angiography images of a 20-year-old female patient who suffered from uncontrollable hypertension using the image processing package Mimics (Materialise). The morphology and hemodynamic parameters in the healthy state, at initial presentation and at post-operative 9-month and 2-year follow-ups after surgery were analysed. The hemodynamic parameters became normal in the left and right renal arteries after bypass surgery. However, flow separation and stagnation occurred at the post-operative 9-month aorta-vein anastomosis, which caused asymmetrical flow and extremely high wall shear stress (WSS) and WSS gradients at the outflow vein tract, where the stenosis occurred 2 years later. In addition, the graft bending produced an asymmetrical flow pattern downstream. This research revealed that the abnormal hemodynamics, including flow separation and extremely high WSS values and gradients, caused by the retrograde flow of aortorenal bypass may be responsible for the SVG degeneration. In addition, flow asymmetry due to vessel bending is a potential risk factor for SVG aneurysm dilation.

PIV investigation of the flow fields in subject-specific vertebro-basilar (VA-BA) junction

Background

As the only arterial structure of which two main arteries merged into one, the vertebro-basilar (VA-BA) system is one of the favorite sites of cerebral atherosclerotic plaques. The aim of this study was to investigate the detailed hemodynamics characteristics in the VA-BA system.

Methods

A scale-up subject-specific flow phantom of VA-BA system was fabricated based on the computed tomography angiography (CTA) scanning images of a healthy adult. Flow fields in eight axial planes and six radial planes were measured and analyzed by using particle image velocimetry (PIV) under steady flow conditions of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Re}=300$$\end{document}Re=300, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Re}=500$$\end{document}Re=500. A water–glycerin mixture was used as the working fluid.

Results

The flow in the current model exhibited highly three-dimensional characteristics. The confluence of VAs flow formed bimodal velocity distribution near the confluence apex. Due to the asymmetrical structural configuration, the bimodal velocity profile skewed towards left, and sharper peaks were observed under higher Reynolds condition. Secondary flow characterized by two vortices formed in the radial planes where 10 mm downstream the confluence apex and persists along the BA under both Reynolds numbers. The strength of secondary flow under \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Re}=500$$\end{document}Re=500 is around 8% higher than that under \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Re}=300$$\end{document}Re=300, and decayed nonlinearly along the flow direction. In addition, a low momentum recirculation region induced by boundary layer separation was observed near the confluence apex. The wall shear stress (WSS) in the recirculation area was found to be lower than 0.4 Pa. This region coincides well with the preferential site of vascular lesions in the VA-BA system.

Conclusions

This preliminary study verified that the subject-specific in-vitro experiment is capable of reflecting the detailed flow features in the VA-BA system. The findings from this study may help to expand the understanding of the hemodynamics in the VA-BA system, and further clarifying the mechanism that underlying the localization of vascular lesions.

Intelligent Diabetes Detection System based on Tongue Datasets

Background:

Scanning Electron Microscope (SEM) Camera Imaging shows and helps analyze hidden organs in the human body. SEM image analysis provides in-depth and critical details of organ abnormalities. Similarly, the human tongue finds use in the detection of organ dysfunction with tongue reflexology.

Objective:

To detect diabetes at an early stage using a non-invasive method of diabetes detection through tongue images and to utilize the reasonable cost of modality (SEM camera) for capturing the tongue images instead of the existing and expensive imaging modalities like X-ray, Computed Tomography, Magnetic Resonance Imaging, Positron Emission Tomography, Single-Photon Emission Computed Tomography etc.

Methods:

The tongue image is captured via SEM camera, it is preprocessed to remove noise and resize the tongue such that it is suitable for segmentation. Greedy Snake Algorithm (GSA) is used to segment the tongue image. The texture features of the tongue are analyzed and finally it is classified as diabetic or normal.

Results:

Failure of organs stomach, intestine, liver and pancreas results in change of the color of the tongue, coating thickness and cracks on the tongue. Changes in pancreas proactive behavior also reflect on tongue coating. The tongue coating texture varies from white or vanilla to yellow also the tongue coating thickness also increases.

Conclusion:

In this paper, the author proposes to diagnose Diabetes Type2 (DT2) at an early stage from tongue digital image. The tongue image is acquired and processed with Greedy Snake Algorithm (GSA) to extract edge and texture features.

Imaging modalities to diagnose carotid artery stenosis: progress and prospect

In the past few decades, imaging has been developed to a high level of sophistication. Improvements from one-dimension (1D) to 2D images, and from 2D images to 3D models, have revolutionized the field of imaging. This not only helps in diagnosing various critical and fatal diseases in the early stages but also contributes to making informed clinical decisions on the follow-up treatment profile. Carotid artery stenosis (CAS) may potentially cause debilitating stroke, and its accurate early detection is therefore important. In this paper, the technical development of various CAS diagnosis imaging modalities and its impact on the clinical efficacy is thoroughly reviewed. These imaging modalities include duplex ultrasound (DUS), computed tomography angiography (CTA) and magnetic resonance angiography (MRA). For each of the imaging modalities considered, imaging methodology (principle), critical imaging parameters, and the extent of imaging the vulnerable plaque are discussed. DUS is usually the initial recommended CAS diagnostic examination. However, for the therapeutic intervention, either MRA or CTA is recommended for confirmation, and for added information on intracranial cerebral circulation and aortic arch condition for procedural planning. Over the past few decades, the focus of CAS diagnosis has also shifted from pure stenosis quantification to plaque characterization. This has led to further advancement in the existing imaging tools and development of other potential imaging tools like Optical coherence tomography (OCT), photoacoustic tomography (PAT), and infrared (IR) thermography.

A novel vision-based system for quantitative analysis of abdominal aortic aneurysm deformation

BACKGROUND

In clinical diagnostics, combination of different imaging techniques is applied to assess spatial configuration of the abdominal aortic aneurysm (AAA) and deformation of its wall. As deformation of aneurysm wall is crucial parameter in assessing wall rupture, we aimed to develop and validate a Non-Invasive Vision-Based System (NIVBS) for the analysis of 3D elastic artificial abdominal aortic models. 3D-printed elastic AAA models from four patients were applied for the reconstruction of real hemodynamic. During experiments, the inlet boundary conditions included the injection volume and frequency of pulsation averaged from electrocardiography traces. NIVBS system was equipped with nine cameras placed at a constant distance to record wall movement from 360^o angle and a dedicated set of artificial lights providing coherent illumination. Additionally, self-prepared algorithms for image acquisition, processing, segmentation, and contour detection were used to analyze wall deformation. Finally, the shape deformation factor was applied to evaluate aorta's deformation. Experimental results were confronted with medical data from AngioCT and 2D speckle-tracking echocardiography (2DSTE).

RESULTS

Image square analyses indicated that the optimal distance between the camera's lens and the investigated object was in the range of 0.30-0.35 m. There was approximately 1.44% difference observed in aneurysm diameters between NIVBS (86.57 ± 5.86 mm) and AngioCT (87.82 ± 6.04 mm) (p = 0.7764). The accuracy of developed algorithm for the reconstruction of the AAA deformation was equal to 98.56%. Bland-Altman analysis showed that the difference between clinical data (2DSTE) and predicted wall deformation (NIVBS) for all patients was 0.00 mm (confidence interval equal to 0.12 mm) for aneurysm size, 0.01 mm (confidence interval equal to 0.13 mm) and 0.00 mm (confidence interval equal to 0.09 mm) for the anterior and posterior side, as well as 0.01 mm (confidence interval equal to 0.18 mm) and 0.01 mm (confidence interval equal to 0.11 mm) for the left and right side. The optimal range of camera's lens did not affect acquired values.

CONCLUSIONS

The NIVBS with proposed algorithm that reconstructs the pressure from surrounding organs is appropriate to analyze the AAAs in water environment. Moreover, NIVBS allowed detailed quantitative analysis of aneurysm sac wall deformation.

COMPUTATIONAL INVESTIGATION OF THROMBIN CONCENTRATION IN CEREBRAL ANEURYSMS TREATED WITH FLOW-DIVERTING STENTS

Flow-diverting stent is an ongoing embolization device to treat cerebral aneurysms, and it diverts the flow direction to reduce the flow velocity inside the aneurysmal sacs and promote the thrombus formation. However, its effect for aneurysm embolization is controversial. A hemodynamic-biomedical coupling model was constructed to describe the generation and transport of thrombin in arteries, and the model was applied to investigate the variation of thrombin concentration, which plays a key role in thrombus formation, in two patient-specific cerebral aneurysm models when they are treated with Pipeline flow diverting stents. It is observed from computational fluid dynamics simulations that thrombin concentration in the aneurysmal sac without collateral artery increases significantly after Pipeline implantation, however, it has hardly any variation in the aneurysmal sac without collateral artery or in the giant aneurysmal sac after Pipeline implantation. Therefore, we believe that single Pipeline is very effective to embolize a small aneurysm without collateral artery, but cannot embolize a giant aneurysm or a small aneurysm with a collateral artery on its sac effectively.

Double stent insertion for combined malignant airway and superior vena cava obstruction

The present study sought to assess the feasibility and effectiveness of double stent insertion as a means of managing combined malignant airway and superior vena cava (SVC) obstruction (CMASO).From July 2010 to January 2018, twelve consecutive patients with CMASO were treated by double stent insertion (airway and SVC stents) in our centers. We assessed data pertaining to both technical and clinical success rates, as well as to long-term patient outcomes.The use of double stents (12 airway stents and 19 SVC stents) was technically successful in all study subjects, with a 0 to 92 days period between the 2 stent insertions (mean 27.8 days). Patients did not show evidence of any procedure-related complications. Mean patient Hugh-Jones grades improved from 4.4 ± 0.5 before inserting the airway stent down to 1.2 ± 0.4 following this insertion (P < .001). Mean SVC pressure was reduced from 17.5 ± 2.8 mm Hg before the stent insertion down to 6.7 ± 1.4 mm Hg following this insertion (P < .001). Fifty-six days after insertion, a single patient experienced re-obstruction of their SVC stent. All patients died within the follow-up period, with a median survival time of 113 days for these 12 patients.This double stent insertion protocol is both effective and safe as a means of offering palliative care to those with CMASO.

Hydraulic conductivity and low-density lipoprotein transport of the venous graft wall in an arterial bypass

Background

Blood flow condition may have influence upon the hydraulic conductivity of venous graft (L_p,vein) in an arterial bypass, then affecting the accumulation of low-density lipoproteins (LDLs) within the graft wall. To probe this possibility, we first measured in vitro the filtration rates of swine lateral saphenous vein segments under different flow rates, and the correlation of L_p,vein with wall shear stress (WSS) was then obtained.

Results

The experimental results showed that when WSS was very low, L_p,vein would increase drastically with WSS from 1.16 ± 0.15 × 10⁻¹¹ m/s Pa at 0 dyn/cm² to 2.17 ± 0.20 × 10⁻¹¹ m/s Pa at 0.7 dyn/cm², then became constant of approximately 2.33 × 10⁻¹¹ m/s Pa as the WSS increased further. Based on the experimental results, we assumed three different cases of L_p,vein and numerically simulated the LDLs transport in an arterial bypass model with venous graft. Case A: L_p,vein = 2.33 × 10⁻¹¹ m/s Pa; Case B: L_p,vein = 1.16 × 10⁻¹¹ m/s Pa (static condition with WSS of 0); Case C: L_p,vein was shear dependent. The simulation showed that the deposition/accumulation of LDLs within the venous graft wall in Case A was greatly enhanced when compared with that in Case B. However, the LDL accumulation in the graft wall was similar for Case A and Case C.

Conclusions

Our study, therefore, indicates that when the venous graft was implanted as a bypass graft, the L_p,vein might remain nearly constant along its whole length except for very few areas where the value of WSS was extremely low (less than 0.7 dyn/cm²) and the effects of L_p,vein modulated by blood flow on LDL transport may be neglected.

Investigating heartbeat-related in-plane motion and stress levels induced at the aortic root

Background

The axial motion of aortic root (AR) due to ventricular traction was previously suggested to contribute to ascending aorta (AA) dissection by increasing its longitudinal stress, but AR in-plane motion effects on stresses have never been studied. The objective is to investigate the contribution of AR in-plane motion to AA stress levels.

Methods

The AR in-plane motion was assessed on magnetic resonance imagining data from 25 healthy volunteers as the movement of the AA section centroid. The measured movement was prescribed to the proximal AA end of an aortic finite element model to investigate its influences on aortic stresses. The finite element model was developed from a patient-specific geometry using LS-DYNA solver and validated against the aortic distensibility. Fluid–structure interaction (FSI) approach was also used to simulate blood hydrodynamic effects on aortic dilation and stresses.

Results

The AR in-plane motion was 5.5 ± 1.7 mm with the components of 3.1 ± 1.5 mm along the direction of proximal descending aorta (PDA) to AA centroid and 3.0 ± 1.3 mm perpendicularly under the PDA reference system. The AR axial motion elevated the longitudinal stress of proximal AA by 40% while the corresponding increase due to in-plane motion was always below 5%. The stresses at proximal AA resulted approximately 7% less in FSI simulation with blood flow.

Conclusions

The AR in-plane motion was comparable with the magnitude of axial motion. Neither axial nor in-plane motion could directly lead to AA dissection. It is necessary to consider the heterogeneous pressures related to blood hydrodynamics when studying aortic wall stress levels.

Electronic supplementary material

The online version of this article (10.1186/s12938-019-0632-7) contains supplementary material, which is available to authorized users.

Effects on the pulmonary hemodynamics and gas exchange with a speed modulated right ventricular assist rotary blood pump: a numerical study

Rotary blood pumps (RBPs) are the newest generation of ventricular assist devices. Although their continuous flow characteristics have been accepted widely, more and more research has focused on the pulsatile modulation of RBPs in an attempt to provide better perfusion. In this study, we investigated the effects of an axial RBP serving as the right ventricular assist device on pulmonary hemodynamics and gas exchange using a numerical method with a complete cardiovascular model along with airway mechanics and a gas exchange model. The RBP runs in both constant speed and synchronized pulsatile modes using speed modulation. Hemodynamics and airway O₂ and CO₂ partial pressures were obtained under normal physiological conditions, and right ventricle failure conditions with or without RBP. Our results showed that the pulsatile mode of the RBP could support right ventricular assist to restore most hemodynamics. Using speed modulation, both pulmonary arterial pressure and flow pulsatility were increased, while there was only very little effect on alveolar O₂ and CO₂ partial pressures. This study could provide basic insight into the influence of pulmonary hemodynamics and gas exchange with speed modulated right ventricular assist RBPs, which is concerned when designing their pulsatile control methods.

Using the Pulse Contour Method to Measure the Changes in Stroke Volume during a Passive Leg Raising Test

The pulse contour method is often used with the Windkessel model to measure stroke volume. We used a digital pressure and flow sensors to detect the parameters of the Windkessel model from the pulse waveform. The objective of this study was to assess the stability and accuracy of this method by making use of the passive leg raising test. We studied 24 healthy subjects (40 ± 9.3 years), and used the Medis^® CS 1000, an impedance cardiography, as the comparing reference. The pulse contour method measured the waveform of the brachial artery by using a cuff. The compliance and resistance of the peripheral artery was detected from the cuff characteristics and the blood pressure waveform. Then, according to the method proposed by Romano et al., the stroke volume could be measured. This method was implemented in our designed blood pressure monitor. A passive leg raising test, which could immediately change the preloading of the heart, was done to certify the performance of our method. The pulse contour method and impedance cardiography simultaneously measured the stroke volume. The measurement of the changes in stroke volume using the pulse contour method had a very high correlation with the Medis^® CS 1000 measurement, the correlation coefficient of the changed ratio and changed differences in stroke volume were r² = 0.712 and r² = 0.709, respectively. It was shown that the stroke volume measured by using the pulse contour method was not accurate enough. But, the changes in the stroke volume could be accurately measured with this pulse contour method. Changes in stroke volume are often used to understand the conditions of cardiac preloading in the clinical field. Moreover, the operation of the pulse contour method is easier than using impedance cardiography and echocardiography. Thus, this method is suitable to use in different healthcare fields.

Blood homocysteine levels could predict major adverse cardiac events in patients with acute coronary syndrome: A STROBE-compliant observational study

The Global Registry of Acute Coronary Events (GRACE) risk score independently predicts major adverse cardiac events (MACEs) in patients with acute coronary syndrome (ACS). This study aims to evaluate whether the level of plasma homocysteine in addition to the GRACE score enhances the predictive value for MACEs in patients with acute coronary syndrome.A total of 361 patients with ACS evaluated at our hospital were included in the study and tested for blood homocysteine levels. We recorded 40 (11.1%) instances of MACE during a median follow-up of 43.3 months (quartile 40.6-44.4 months), including 29 cases (8.0%) of all-cause death and 11 cases (3.1%) of nonfatal myocardial infarction.The GRACE score was significantly associated with homocysteine levels, and multivariate Cox regression analysis showed that both the GRACE risk score and homocysteine content were independent predictors of MACEs (HR 2.63; 95% confidence interval (CI) 1.54 to 4.49; P < .001 and 2.27; 1.06 to 4.86; P = .035, respectively). Moreover, meta-analysis showed that as the homocysteine level increased, the incidence of MACEs also increased (log-rank 8.41; P = .015). GRACE scores adjusted by homocysteine level increased the area under the curve (AUC) from 0.78 to 0.83 (P = 0.006).Blood homocysteine levels are significantly associated with the GRACE risk score, and using both parameters can further improve risk stratification in patients with acute coronary syndrome.

Quantifying Hemodynamic Changes in Moyamoya Disease Based on Two-Dimensional Cine Phase-Contrast Magnetic Resonance Imaging and Computational Fluid Dynamics

OBJECTIVE

To investigate hemodynamic changes in moyamoya disease (MMD) via two-dimensional cine phase-contrast magnetic resonance imaging and computational fluid dynamics.

METHODS

In 18 patients with MMD and 10 healthy control subjects, phase-contrast magnetic resonance imaging was performed to quantify flow rate of main supplying arteries, including internal carotid arteries (ICAs) and vertebral arteries. Mean flow rate in these vessels was adopted as the patient-specific boundary condition for computational fluid dynamics simulation of the circle of Willis in MMD and control groups. Pressure drop in both ICAs and their difference, wall shear stress and secondary flow in the carotid siphon of ICAs, and flow rate and size of posterior communicating arteries (PComAs) were compared between MMD and control groups. Four patients with MMD underwent follow-up scans for longitudinal comparison.

RESULTS

Phase-contrast magnetic resonance imaging data revealed significantly different flow rate in the left ICA and right vertebral arteries between MMD and control groups. Computational fluid dynamics simulation demonstrated similar wall shear stress and similar secondary flow of both ICAs but significantly higher pressure drop in left ICA, higher pressure drop difference between left ICA and right ICA, and higher flow rate in PComAs in patients with MMD compared with control subjects. Significantly increased size of left PComA in patients with MMD was also found. Follow-up results confirmed that the combination of pressure drop difference, flow rate, and size of PComAs can potentially assist long-term prognosis after surgery.

CONCLUSIONS

Pressure drop difference, flow rate, and size of PComAs can be used to evaluate impairments in cerebrovascular reserve and indicate long-term prognosis in MMD.

Pulse wave response characteristics for thickness and hardness of the cover layer in pulse sensors to measure radial artery pulse

Background

Piezo-resistive pressure sensors are widely used for measuring pulse waves of the radial artery. Pulse sensors are generally fabricated with a cover layer because pressure sensors without a cover layer are fragile when they come into direct contact with the skin near the radial artery. However, no study has evaluated the dynamic pulse wave response of pulse sensors depending on the thickness and hardness of the cover layer. This study analyzed the dynamic pulse wave response according to the thickness and hardness of the cover layer and suggests an appropriate thickness and hardness for the design of pulse sensors with semiconductor device-based pressure sensors.

Methods

Pulse sensors with 6 different cover layers with various thicknesses (0.8 mm, 1 mm, 2 mm) and hardnesses (Shore type A; 30, 43, 49, 71) were fabricated. Experiments for evaluating the dynamic pulse responses of the fabricated sensors were performed using a pulse simulator to transmit the same pulse wave to each of the sensors. To evaluate the dynamic responses of the fabricated pulse sensors, experiments with the pulse sensors were conducted using a simulator that artificially generated a constant pulse wave. The pulse wave simulator consisted of a motorized cam device that generated the artificial radial pulse waveform by adjusting the stroke of the cylindrical air pump and an air tube that conveyed the pulse to the artificial wrist.

Results

The amplitude of the measured pulse pressure decreased with increasing thickness and hardness of the cover layer. Normalized waveform analysis showed that the thickness rather than the hardness of the cover layer contributed more to waveform distortion. Analysis of the channel distribution of the pulse sensor with respect to the applied constant dynamic pressure showed that the material of the cover layer had a large effect.

Conclusions

In this study, in-line array pulse sensors with various cover layers were fabricated, the dynamic pulse wave responses according to the thickness and the hardness of the cover layer were analyzed, and an appropriate thickness and hardness for the cover layer were suggested. The dynamic pulse wave responses of pulse sensors revealed in this study will contribute to the fabrication of improved pulse sensors and pulse wave analyses.

The feasibility of 1-stop examination of coronary CT angiography and abdominal enhanced CT

This study aims to evaluate the feasibility of performing coronary computed tomography angiography (CCTA) and abdominal enhanced computed tomography (CT) with 1-time injection of the agent.CCTA images (right coronary artery, left anterior descending coronary artery, and left circumflex coronary artery) were collected from 20 patients who completed a 1-stop combined examination of CCTA and abdominal enhanced CT (group A), 20 patients who only underwent abdominal enhanced CT (group B1), and 20 patients who only underwent CCTA (group B2). These images were interpreted using the 5-point Likert scale system by 2 experienced radiologists, and abdominal images were observed for breathing artifact. CT value, signal-to-noise ratio (SNR), and CTDI were recorded and compare among the 3 groups.The difference in image quality of the coronary and total volume of the contrast agent between group A and group B1 was not statistical significant (P > .05). The CT value and SNR in group B1 (CCTA) (CT: 394.65 ± 59.23, SNR: 17.38 ± 4.13) increased, compare with Group A (CT: 360.35 ± 34.16, SNR: 13.76 ± 1.84, P = .03, .01), while CTDI was undifferentiated between group A (17.14 ± 6.20) and group B1 (18.38 ± 9.79) (P = .64). The difference in CT value and SNR at the arterial phase and CT value at the venous phase between group A (abdomen) and group B2 were statistically significant, the CTDI in group A (9.09 ± 1.05) increased, compared with group B2 (8.23 ± 1.33) (P = .03), and SNR at the venous phase in group B2 (12.50 ± 2.43) increased, compared with group A (10.89 ± 2.03) (P = .03).Revolution CT can capture full images and very rapidly switch to the scan mode, enabling a 1-stop axial CCTA and enhanced helical abdominal scan. The 1-stop combined scan resulted in a satisfactory image quality, which reduced the contrast agent dose and simplified the workflow.The 1-stop combined scan allows for the high success rate of the examination, reduces the number of examinations, and decreases the dose and risk of injection of the contrast agent. This would be helpful for patients to obtain diagnostic images in time.

Comparisons of the accuracy of radiation diagnostic modalities in brain tumor: A nonrandomized, nonexperimental, cross-sectional trial

Tumor morphology improved sensitivity, accuracy, and specificity of the diagnosis, but all diagnostic techniques have attenuation correction issues.To compare computed tomographic (CT), positron emission tomographic (PET), and magnetic resonance imaging (MRI) characteristics of patients with brain tumor in a Chinese setting.A nonrandomized, nonexperimental, cross-sectional trial.Jining No. 1 People's Hospital, China.In total, 127 patients who had clinically confirmed a brain tumor were included in the cross-sectional study. Patients were subjected to brain CT, MRI, and PET. The tumors resected after brain surgery were subjected to morphological diagnosis. Statistical analysis of data of surgically removed tumor and that of different methods of diagnosis was performed using Wilcoxon test following Tukey-Kramer test. Spearmen correlation was performed between diagnostic modalities and in vivo morphology. Results were considered significant at 99% of confidence level.The data of diameter and volume of tumor derived from CT (Spearman r = 0.9845 and 0.9706), and MRI (Spearman r = 0.955 and 0.2378) were failed to correlate with that of that of the surgically removed tumor. However, prediction of diameter and volume of the tumor by PET (Spearman r = 0.9922 and 0.9921) were correlated with that of the surgically removed tumor. CT and MRI were failed to quantified pituitary adenomas.The study was recommended PET for assessment of brain tumor.

L1 incidence reflects pelvic incidence and lumbar lordosis mismatch in sagittal balance evaluation

Retrospective study.To investigate the radiologic and geometrical association between L1 incidence (L1I) with pelvic incidence/lumbar lordosis (PI/LL) mismatch and T1 incidence (T1I) with PI/LL/thoracic kyphosis (TK) mismatch.The relationship between PI and LL is not clear, and it might be because of the absence of a direct radiologic parameter to represent PI/LL mismatch. To the best of our knowledge, this is the first report on a direct radiologic parameter for representing PI/LL mismatch.This study is a retrospective review of 146 patients who underwent anteroposterior and lateral standing radiographs of the whole spine. L1I was defined as the angle between the line perpendicular to the L1 upper endplate and the line connecting the midpoint of the sacral endplate to the center of both femoral heads. T1I was defined as the angle between the line perpendicular to the T1 upper endplate and the line connecting the midpoint of the sacral endplate to the center of both femoral heads. Both were validated using the Pearson correlation coefficient and linear regression analysis.Radiologically measured L1I and T1I were coterminous with calculated measurements of ΔPI/LL and ΔPI/LL/TK in terms of means and standard deviations, respectively. Excellent correlations were found between L1I and ΔPI/LL, and T1I and ΔPI/LL/TK (R = 0.997, P < .01; R = 0.981, P < .01, respectively). In linear regression analysis, the slope and intercept of L1I were 0.991 and -0.041, with a predictability of 99.4% (R = 0.994), and those of T1I were 0.990 and -0.026, with a predictability of 99.0% (R = 0.990), respectively.L1I and T1I were strongly correlated with PI/LL mismatch and PI/LL/TK mismatch, respectively. L1I and T1I are direct parameters that represent PI/LL mismatch and PI/LL/TK mismatch. They would be useful in analyzing sagittal balance.

LEVEL OF EVIDENCE

Level 3.

The effect of inlet and outlet boundary conditions in image-based CFD modeling of aortic flow

Background

Computational modeling of cardiovascular flow is a growing and useful field, but such simulations usually require the researcher to guess the flow’s inlet and outlet conditions since they are difficult and expensive to measure. It is critical to determine the amount of uncertainty introduced by these assumptions in order to evaluate the degree to which cardiovascular flow simulations are accurate. Our work begins to address this question by examining the sensitivity of flow to several different assumed velocity inlet and outlet conditions in a patient-specific aorta model.

Methods

We examined the differences between plug flow, parabolic flow, linear shear flows, skewed cubic flow profiles, and Womersley flow at the inlet. Only the shape of the inlet velocity profile was varied—all other parameters were identical among these simulations. Secondary flow in the form of a counter-rotating pair of vortices was also added to parabolic axial flow to study its effect on the solution. In addition, we examined the differences between two-element Windkessel, three element Windkessel and the outflow boundary conditions. In these simulations, only the outlet boundary condition was varied.

Results

The results show axial and in-plane velocities are considerably different close to the inlet for the cases with different inlet velocity profile shapes. However, the solutions are qualitatively similar beyond 1.75D, where D is the inlet diameter. This trend is also observed in other quantities such as pressure and wall shear stress. Normalized root-mean-square deviation, a measure of axial velocity magnitude differences between the different cases, generally decreases along the streamwise coordinate. The linear shear inlet velocity boundary condition and plug velocity boundary condition solution exhibit the highest time-averaged wall shear stress, approximately \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$8\%$$\end{document}8% higher than the parabolic inlet velocity boundary condition. Upstream of 1D from the inlet, adding secondary flow has a significant impact on temporal wall shear stress distributions. This is especially observable during diastole, when integrated wall shear stress magnitude varies about \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$26\%$$\end{document}26% between simulations with and without secondary flow. The results from the outlet boundary condition study show the Windkessel models differ from the outflow boundary condition by as much as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$18\%$$\end{document}18% in terms of time-averaged wall shear stress. Furthermore, normalized root-mean-square deviation of axial velocity magnitude, a measure of deviation between Windkessel and the outflow boundary condition, increases along the streamwise coordinate indicating larger variations near outlets.

Conclusion

It was found that the selection of inlet velocity conditions significantly affects only the flow region close to the inlet of the aorta. Beyond two diameters distal to the inlet, differences in flow solution are small. Although additional studies must be performed to verify this result, the data suggest that it is important to use patient-specific inlet conditions primarily if the researcher is concerned with the details of the flow very close to the inlet. Similarly, the selection of outlet conditions significantly affects the flow in the vicinity of the outlets. Upstream of five diameters proximal to the outlet, deviations between the outlet boundary conditions examined are insignificant. Although the inlet and outlet conditions only affect the flow significantly in their respective neighborhoods, our study indicates that outlet conditions influence a larger percentage of the solution domain.

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.

Hemodynamics of cerebral bridging veins connecting the superior sagittal sinus based on numerical simulation

Background

The physiological and hemodynamic features of bridging veins involve wall shear stress (WSS) of the cerebral venous system. Based on the data of cadavers and computational fluid dynamics software pack, the hemodynamic physical models of bridging veins (BVs) connecting superior sagittal sinus (SSS) were established.

Results

A total of 137 BVs formed two clusters along the SSS: anterior group and posterior group. The diameters of the BVs in posterior group were larger than of the anterior group, and the entry angle was smaller. When the diameter of a BV was greater than 1.2 mm, the WSS decreased in the downstream wall of SSS with entry angle less than 105°, and the WSS also decreased in the upstream wall of BVs with entry angle less than 65°. The minimum WSS in BVs was only 63% of that in SSS. Compared with the BVs in anterior group, the minimum WSS in the posterior group was smaller, and the distance from location of the minimum WSS to the dural entrance was longer.

Conclusion

The cerebral venous thrombosis occurs more easily when the diameter of a BV is greater than 1.2 mm and the entry angle is less than 65°. The embolus maybe form earlier in the upstream wall of BVs in the posterior part of SSS.

Electronic supplementary material

The online version of this article (10.1186/s12938-018-0466-8) contains supplementary material, which is available to authorized users.

Carotid DSA based CFD simulation in assessing the patient with asymptomatic carotid stenosis: a preliminary study

Background

Cerebrovascular events are frequently associated with hemodynamic disturbance caused by internal carotid artery (ICA) stenosis. It is challenging to determine the ischemia-related carotid stenosis during the intervention only using digital subtracted angiography (DSA). Inspired by the performance of well-established FFRct technique in hemodynamic assessment of significant coronary stenosis, we introduced a pressure-based carotid arterial functional assessment (CAFA) index generated from computational fluid dynamic (CFD) simulation in DSA data, and investigated its feasibility in the assessment of hemodynamic disturbance preliminarily using pressure-wired measurement and arterial spin labeling (ASL) MRI as references.

Methods

The cerebral multi-delay multi-parametric ASL-MRI and carotid DSA including trans-stenotic pressure-wired measurement were implemented on a 65-year-old man with asymptomatic unilateral (left) ICA stenosis. A CFD simulation using simplified boundary condition was performed in DSA data to calculate the CAFA index. The cerebral blood flow (CBF) and arterial transit time (ATT) of ICA territories were acquired.

Results

CFD simulation showed good correlation (r = 0.839, P = 0.001) with slight systematic overestimation (mean difference − 0.007, standard deviation 0.017) compared with pressure-wired measurement. No significant difference was observed between them (P = 0.09). Though the narrowing degree of in the involved ICA was about 70%, the simulated and measured CAFA (0.942/0.937) revealed a functionally nonsignificant stenosis which was also verified by a compensatory final CBF (fronto-temporal/fronto-parietal region: 51.58/45.62 ml/100 g/min) and slightly prolonged ATT (1.23/1.4 s) in the involved territories, together with a normal left–right percentage difference (2.1–8.85%).

Conclusions

The DSA based CFD simulation showed good consistence with invasive approach and could be used as a cost-saving and efficient way to study the relationship between hemodynamic disorder caused by ICA stenosis and subsequent perfusion variations in brain. Further research should focus on the role of noninvasive pressure-based CAFA in screening asymptomatic ischemia-causing carotid stenosis.

A patient-specific lumped-parameter model of coronary circulation

A new lumped-parameter model for coronary hemodynamics is developed. This model is developed for the whole coronary network based on CT scans of a patient-specific geometry including the right coronary tree, which is absent in many previous mathematical models. The model adopts the structured tree model boundary conditions similar to the work of Olufsen et al., thus avoiding the necessity of invasive perfusion measurements. In addition, we also incorporated the effects of the head loss at the two inlets of the large coronary arteries for the first time. The head loss could explain the phenomenon of a sudden increase of the resistance at the inlet of coronary vessel. The estimated blood pressure and flow rate results from the model agree well with the clinical measurements. The computed impedances also match the experimental perfusion measurement. The effects of coronary arterial stenosis are considered and the fractional flow reserve and relative flow in the coronary vessels for a stenotic vessel computed in this model show good agreement with published experimental data. It is believed that the approach could be readily translated to clinical practice to facilitate real time clinical diagnosis.

Patch-based local learning method for cerebral blood flow quantification with arterial spin-labeling MRI

Arterial spin-labeling (ASL) perfusion MRI is a non-invasive method for quantifying cerebral blood flow (CBF). Standard ASL CBF calibration mainly relies on pair-wise subtraction of the spin-labeled images and controls images at each voxel separately, ignoring the abundant spatial correlations in ASL data. To address this issue, we previously proposed a multivariate support vector machine (SVM) learning-based algorithm for ASL CBF quantification (SVMASLQ). But the original SVMASLQ was designed to do CBF quantification for all image voxels simultaneously, which is not ideal for considering local signal and noise variations. To fix this problem, we here in this paper extended SVMASLQ into a patch-wise method by using a patch-wise classification kernel. At each voxel, an image patch centered at that voxel was extracted from both the control images and labeled images, which was then input into SVMASLQ to find the corresponding patch of the surrogate perfusion map using a non-linear SVM classifier. Those patches were eventually combined into the final perfusion map. Method evaluations were performed using ASL data from 30 young healthy subjects. The results showed that the patch-wise SVMASLQ increased perfusion map SNR by 6.6% compared to the non-patch-wise SVMASLQ.

Large-scale subject-specific cerebral arterial tree modeling using automated parametric mesh generation for blood flow simulation

In this paper, we present a novel technique for automatic parametric mesh generation of subject-specific cerebral arterial trees. This technique generates high-quality and anatomically accurate computational meshes for fast blood flow simulations extending the scope of 3D vascular modeling to a large portion of cerebral arterial trees. For this purpose, a parametric meshing procedure was developed to automatically decompose the vascular skeleton, extract geometric features and generate hexahedral meshes using a body-fitted coordinate system that optimally follows the vascular network topology. To validate the anatomical accuracy of the reconstructed vasculature, we performed statistical analysis to quantify the alignment between parametric meshes and raw vascular images using receiver operating characteristic curve. Geometric accuracy evaluation showed an agreement with area under the curves value of 0.87 between the constructed mesh and raw MRA data sets. Parametric meshing yielded on-average, 36.6% and 21.7% orthogonal and equiangular skew quality improvement over the unstructured tetrahedral meshes. The parametric meshing and processing pipeline constitutes an automated technique to reconstruct and simulate blood flow throughout a large portion of the cerebral arterial tree down to the level of pial vessels. This study is the first step towards fast large-scale subject-specific hemodynamic analysis for clinical applications.

Minimizing the blood velocity differences between phase-contrast magnetic resonance imaging and computational fluid dynamics simulation in cerebral arteries and aneurysms

The integration of phase-contrast magnetic resonance images (PC-MRI) and computational fluid dynamics (CFD) is a way to obtain detailed information of patient-specific hemodynamics. This study proposes a novel strategy for imposing a pressure condition on the outlet boundary (called the outlet pressure) in CFD to minimize velocity differences between the PC-MRI measurement and the CFD simulation, and to investigate the effects of outlet pressure on the numerical solution. The investigation involved ten patient-specific aneurysms reconstructed from a digital subtraction angiography image, specifically on aneurysms located at the bifurcation region. To evaluate the effects of imposing the outlet pressure, three different approaches were used, namely: a pressure-fixed (P-fixed) approach; a flow rate control (Q-control) approach; and a velocity-field-optimized (V-optimized) approach. Numerical investigations show that the highest reduction in velocity difference always occurs in the V-optimized approach, where the mean of velocity difference (normalized by inlet velocity) is 19.3%. Additionally, the highest velocity differences appear near to the wall and vessel bifurcation for 60% of the patients, resulting in differences in wall shear stress. These findings provide a new methodology for PC-MRI integrated CFD simulation and are useful for understanding the evaluation of velocity difference between the PC-MRI and CFD.