Computational Fluid Dynamics Analysis and Empirical Evaluation of Carboxymethylcellulose/Alginate 3D Bioprinting Inks for Screw-Based Microextrusion

Three-dimensional microextrusion bioprinting technology uses pneumatics, pistons, or screws to transfer and extrude bioinks containing biomaterials and cells to print biological tissues and organs. Computational fluid dynamics (CFD) analysis can simulate the flow characteristics of bioinks in a control volume, and the effect on cell viability can be predicted by calculating the physical quantities. In this study, we developed an analysis system to predict the effect of a screw-based dispenser system (SDS) on cell viability in bioinks through rheological and CFD analyses. Furthermore, carboxymethylcellulose/alginate-based bioinks were used for the empirical evaluation of high-viscous bioinks. The viscosity of bioinks was determined by rheological measurement, and the viscosity coefficient for the CFD analysis was derived from a correlation equation by non-linear regression analysis. The mass flow rate derived from the analysis was successfully validated by comparison with that from the empirical evaluation. Finally, the cell viability was confirmed after bioprinting with bioinks containing C2C12 cells, suggesting that the developed SDS may be suitable for application in the field of bioengineering. Consequently, the developed bioink analysis system is applicable to a wide range of systems and materials, contributing to time and cost savings in the bioengineering industry.

Evaluation of the shrinkage process of a neck remnant after stent-coil treatment of a cerebral aneurysm using silent magnetic resonance angiography and computational fluid dynamics analysis: illustrative case

BACKGROUND

Silent magnetic resonance angiography (MRA) mitigates metal artifacts, facilitating clear visualization of neck remnants after stent and coil embolization of cerebral aneurysms. This study aims to scrutinize hemodynamics at the neck remnant by employing silent MRA and computational fluid dynamics.

OBSERVATIONS

The authors longitudinally tracked images of a partially thrombosed anterior communicating artery aneurysm's neck remnant, which had been treated with stent-assisted coil embolization, using silent MRA over a decade. Computational fluid dynamics delineated the neck remnant's reduction process, evaluating hemodynamic parameters such as flow rate, wall shear stress magnitude and vector, and streamlines. The neck remnant exhibited diminishing surface area, volume, neck size, dome depth, and aspect ratio. Its reduction correlated with a decline in the flow rate ratio of the remnant dome to the inflow parent artery. Analysis delineated regions within the contracting neck remnant characterized by consistently low average wall shear stress magnitude and variation, accompanied by notable variations in wall shear stress vector directionality.

LESSONS

Evaluation of neck remnants after stent-coil embolization is possible through silent MRA and computational fluid dynamics. Predicting the neck remnant reduction may be achievable through hemodynamic parameter analysis.

Effect of blood viscosity on the hemodynamics of arteriovenous fistulae based on numerical investigation

Arteriovenous fistula (AVF) is the most commonly used vascular access for hemodialysis in patients with end-stage renal disease. Vascular diseases such as atherosclerosis and thrombosis, triggered by altered hemodynamic conditions, are the main causes of access failure. Changes in blood viscosity accelerate access dysfunction by affecting local velocities and wall shear stress (WSS) distribution in the circulation. Numerical simulation was employed to analyze and compare the hemodynamic behavior of AVF under different blood viscosities (0.001-0.012 Pa∙s). An idealized three-dimensional model with end-to-side anastomosis was established. Transient simulations were conducted using pulsatile inlet velocity and outflow as boundary conditions. The simulation results reveal the blood flow state of AVF under different viscosity physiological conditions and derive the rule of change. When blood viscosity increases, the local velocity in the disturbed region slows down and the stagnation time becomes longer, resulting in increased deposition of substances. As blood viscosity increases, the level of shear stress on the entire wall of the fistula increases accordingly. WSS values at high viscosities above 0.007 Pa∙s showed significantly larger low-shear regions near the anastomosis and increased chances of inducing atheromatous plaques. This research has revealed the correlation between blood dynamic viscosity and the hemodynamic behavior of AVF. Elevated whole blood viscosity increases the incidence of access obstruction and vascular disease leading to fistula failure. The study provides a basis for optimizing the distribution of hemodynamic parameters in the fistula for hemodialysis patients.

CFD analysis of the hyper-viscous effects on blood flow across abdominal aortic aneurysm in COVID patients: multiphysics approach

Recent research has shown that individuals suffering from COVID-19 are accommodating an elevated level of blood viscosity due to the morphological changes in blood cells. As viscosity is a major flow parameter influencing the flow across a stenosis or an aneurysm, the examination of the significance of hyperviscosity in COVID patients is imperative in arterial pathologies. In this research, we have considered a patient-specific case in which the aneurysm is located along the abdominal aortal walls. Recent research on the side effects of COVID-19 voiced out the various effects on the circulatory system of humans. Also, as abdominal aneurysms exist very often among individuals, causing the death of 150-200 million every year, the hyper-viscous effects of blood on the flow across the diseased aorta are explored by considering the elevated viscosity levels. In vitro explorations contribute considerably to the clinical methods and treatments to be regarded. The objective of the present inquest is to research the flow field in aneurysmatic-COVID-affected patients considering the elastic nature of vessel walls, using the arbitrary Lagrangian-Eulerian approach. The study supports the various clinical findings that voiced the detrimental effects associated with blood hyperviscosity. The simulation results obtained, by solving the fluid mechanics' equations coupled with the solid mechanics' equations, employing a FEM solver suggest that the elevated stress imparted by the hyper-viscous flows on the walls of the aneurysmal aorta can trigger the fastening of the aneurysmal sac enlargement or rupture.

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.

A real-time patient-specific treatment strategy for enhanced external counterpulsation

Diastolic/systolic blood pressure ratio (D/S) ≥ 1.2 is the gold standard of enhanced external counterpulsation (EECP) treatment, but it does not show a clear clinical correspondence with the configuration of the EECP mode. As such, a single target results in different treatment effects in different individuals. The local haemodynamic effect (wall shear stress, WSS) of EECP on vascular endothelial cells is conducive to promote the growth of collateral circulation vessels and restore the blood supply distal to the stenosis lesion. Considering the haemodynamic effects of WSS on human arteries, this study developed a real-time patient-specific treatment strategy of EECP for patients with cardio-cerebrovascular diseases. Based on patient-specific haemodynamic data from 113 individuals, an optimization algorithm was developed to achieve the individualization of a 0D lumped-parameter model of the human circulatory system, thereby simulating the patient-specific global haemodynamic effects. 0D/3D coupled cardio-cerebrovascular models of two subjects were established to simulate the local WSS. We then established statistical models to evaluate clinically unmeasurable WSS based on measurable global haemodynamic indicators. With the aim of attaining appropriate area- and time-averaged WSS (ATAWSS, 4-7 Pa), as evaluated by global haemodynamic indicators, a closed-loop feedback tuning method was developed to provide patient-specific EECP treatment strategies. Results showed that for clinical data collected from 113 individuals, the individualized 0D model can accurately simulate patient-specific global haemodynamic effects (average error <5%). Based on two subjects, the statistical models can be used to evaluate local ATAWSS (error <6%) for coronary arteries and for cerebral arteries. An EECP mode planned by the patient-specific treatment strategy can promote an appropriate ATAWSS within a 16 s calculation time. The real-time patient-specific treatment strategy of EECP is expected to improve the long-term outcome for each patient and have potential clinical significance.

Design of experiment approach to identify the dominant geometrical feature of left coronary artery influencing atherosclerosis

Background and Objective. Coronary artery geometry heavily influences local hemodynamics, potentially leading to atherosclerosis. Consequently, the unique geometrical configuration of an individual by birth can be associated with future risk of atherosclerosis. Although current researches focus on exploring the relationship between local hemodynamics and coronary artery geometry, this study aims to identify the order of influence of the geometrical features through systematic experiments, which can reveal the dominant geometrical feature for future risk assessment.Methods. According to Taguchi's method of design of experiment (DoE), the left main stem (LMS) length (lLMS), curvature (kLMS), diameter (dLMS) and the bifurcation angle between left anterior descending (LAD) and left circumflex (LCx) artery (αLAD-LCx) of two reconstructed patient-specific left coronary arteries (LCA) were varied in three levels to create L9 orthogonal array. Computational fluid dynamic (CFD) simulations with physiological boundary conditions were performed on the resulting eighteen LCA models. Average helicity intensity (h2) and relative atheroprone area (RAA) of near-wall hemodynamic descriptors were analyzed.Results. The proximal LAD (LADproximal) was identified to be the most atheroprone region of the left coronary artery due to higherh2,large RAA of time averaged wall shear stress (TAWSS < 0.4 Pa), oscillatory shear index (OSI ∼ 0.5) and relative residence time (RRT > 4.17 Pa-1). In both patient-specific cases, based onh2and TAWSS,dlmsis the dominant geometric parameter while based on OSI and RRT,αLAD-LCxis the dominant one influencing hemodynamic condition in proximal LAD (p< 0.05). Based on RRT, the rank of the geometrical factors is:αLAD-LCx>dLMS>lLMS>kLMS, indicating thatαLAD-LCxis the most dominant geometrical factor affecting hemodynamics at proximal LAD which may influence atherosclerosis.Conclusion. The proposed identification of the rank of geometrical features of LCA and the dominant feature may assist clinicians in predicting the possibility of atherosclerosis, of an individual, long before it will occur. This study can further be translated to be used to rank the influence of several arterial geometrical features at different arterial locations to explore detailed relationships between the arterial geometrical features and local hemodynamics.

A vertebrobasilar junction aneurysm successfully treated with a combination of surgical clipping and flow diverter placement based on the results of computational fluid dynamics analysis: illustrative case

BACKGROUND

The treatment of vertebrobasilar junction (VBJ) aneurysms is challenging. Although flow diverters (FDs) are a possible treatment option, geometrical conditions hinder intervention. VBJ aneurysms possess dual inflow vessels from the bilateral vertebral arteries (VAs), one of which is ideally occluded prior to FD treatment. However, it remains unclear which VA should be occluded.

OBSERVATIONS

A 75-year-old male with a growing VBJ complex aneurysm exhibiting invagination toward the brainstem and causing perifocal edema required intervention. Preoperative computational fluid dynamics (CFD) analysis demonstrated that left VA occlusion would result in more stagnant flow and less impingement of flow than right VA occlusion. According to the simulated strategy, surgical clipping of the left VA just proximal to the aneurysm was performed, followed by FD placement from the basilar artery trunk to the right VA. The patient demonstrated tolerance of the VA occlusion, and follow-up computed tomography angiography at 18 months after FD treatment confirmed the disappearance of the aneurysm.

LESSONS

Preoperative flow dynamics simulations using CFD analysis can reveal an optimal treatment strategy involving a hybrid surgery that combines FD placement and direct surgical occlusion for a VBJ complex aneurysm.

Multiyear Interval Changes in Aortic Wall Shear Stress in Patients with Bicuspid Aortic Valve Assessed by 4D Flow MRI

BACKGROUND

In patients with bicuspid aortic valve (BAV), 4D flow MRI can quantify regions exposed to abnormal aortic hemodynamics, including high wall shear stress (WSS), a known stimulus for arterial wall dysfunction. However, the long-term multiscan reproducibility of 4D flow MRI-derived hemodynamic parameters is unknown.

PURPOSE

To investigate the long-term stability of 4D flow MRI-derived peak velocity, WSS, and WSS-derived heatmaps in patients with BAV undergoing multiyear surveillance imaging.

STUDY TYPE

Retrospective.

POPULATION

20 BAV patients (mean age 48.4 ± 13.9 years; 14 males) with five 4D flow MRI scans, with intervals of at least 6 months between scans, and 125 controls (mean age: 50.7 ± 15.8 years; 67 males).

FIELD STRENGTH/SEQUENCE

1.5 and 3.0T, prospectively ECG and respiratory navigator-gated aortic 4D flow MRI.

ASSESSMENT

Automated AI-based 4D flow analysis pipelines were used for data preprocessing, aorta 3D segmentation, and quantification of ascending aorta (AAo) peak velocity, peak systolic WSS, and heatmap-derived relative area of elevated WSS compared to WSS ranges in age and sex-matched normative control populations. Growth rate was derived from the maximum AAo diameters measured on the first and fifth MRI scans.

STATISTICAL TESTS

One-way repeated measures analysis of variance. P < 0.05 indicated significance.

RESULTS

One hundred 4D flow MRI exams (five per patient) were analyzed. The mean total follow-up duration was 5.5 ± 1.1 years, and the average growth rate was 0.3 ± 0.2 mm/year. Peak velocity, peak systolic WSS, and relative area of elevated WSS did not change significantly over the follow-up period (P = 0.64, P = 0.69, and P = 0.35, respectively). The patterns and areas of elevated WSS demonstrated good reproducibility on semiquantitative assessment.

CONCLUSION

4D flow MRI-derived peak velocity, WSS, and WSS-derived heatmaps showed good multiyear and multiscan stability in BAV patients with low aortic growth rates. These findings underscore the reliability of these metrics in monitoring BAV patients for potential risk of dilation.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY: Stage 1.

Stress Load and Ascending Aortic Aneurysms: An Observational, Longitudinal, Single-Center Study Using Computational Fluid Dynamics

Ascending aortic aneurysm (AAoA) is a silent disease with high mortality; however, the factors associated with a worse prognosis are not completely understood. The objective of this observational, longitudinal, single-center study was to identify the hemodynamic patterns and their influence on AAoA growth using computational fluid dynamics (CFD), focusing on the effects of geometrical variations on aortic hemodynamics. Personalized anatomic models were obtained from angiotomography scans of 30 patients in two different years (with intervals of one to three years between them), of which 16 (53%) showed aneurysm growth (defined as an increase in the ascending aorta volume by 5% or more). Numerically determined velocity and pressure fields were compared with the outcome of aneurysm growth. Through a statistical analysis, hemodynamic characteristics were found to be associated with aneurysm growth: average and maximum high pressure (superior to 100 Pa); average and maximum high wall shear stress (superior to 7 Pa) combined with high pressure (>100 Pa); and stress load over time (maximum pressure multiplied by the time interval between the exams). This study provides insights into a worse prognosis of this serious disease and may collaborate for the expansion of knowledge about mechanobiology in the progression of AAoA.

Oscillatory shear stress is elevated in patients with bicuspid aortic valve and aortic regurgitation: a 4D flow cardiovascular magnetic resonance cross-sectional study

AIMS

Patients with bicuspid aortic valve (BAV) and aortic regurgitation have higher rate of aortic complications compared with patients with BAV and stenosis, as well as BAV without valvular disease. Aortic regurgitation alters blood haemodynamics not only in systole but also during diastole. We therefore sought to investigate wall shear stress (WSS) during the whole cardiac cycle in BAV with aortic regurgitation.

METHODS AND RESULTS

Fifty-seven subjects that underwent 4D flow cardiovascular magnetic resonance imaging were included: 13 patients with BAVs without valve disease, 14 BAVs with aortic regurgitation, 15 BAVs with aortic stenosis, and 22 normal controls with tricuspid aortic valve. Peak and time averaged WSS in systole and diastole and the oscillatory shear index (OSI) in the ascending aorta were computed. Student's t-tests were used to compare values between the four groups where the data were normally distributed, and the non-parametric Wilcoxon rank sum tests were used otherwise. BAVs with regurgitation had similar peak and time averaged WSS compared with the patients with BAV without valve disease and with stenosis, and no regions of elevated WSS were found. BAV with aortic regurgitation had twice as high OSI as the other groups (P ≤ 0.001), and mainly in the outer mid-to-distal ascending aorta.

CONCLUSION

OSI uniquely characterizes altered WSS patterns in BAVs with aortic regurgitation, and thus could be a haemodynamic marker specific for this specific group that is at higher risk of aortic complications. Future longitudinal studies are needed to verify this hypothesis.

Analysis and numerical investigation of bile flow dynamics within the strictured biliary duct

The mechanics of bile flow in the biliary system plays an important role in studying bile stasis and gallstone formation. Bile duct stricture is an abnormal phenomenon that refers to the bile duct getting smaller or narrower. The main objective of this study is to study the influence of stricture on bile flow dynamics using numerical methods. We employed a numerical Computational Fluid Dynamics model of the bile flow within a strictured hepatic duct. We studied and compared the influence of stricture severity, stricture length, eccentricity, and bile flow property on the bile flow dynamics. The bile flow velocity, pressure distribution, pressure drop, and wall shear stress are provided in detail. The stricture alters the normal bile flow pattern and increases flow resistance. At the location upstream and downstream of the stricture, bile flow slows down. In the area of the stricture throat, bile flow is accelerated, and recirculation forms behind the stricture. The maximum pressure drop of the biliary system increases with the stricture length. The eccentricity makes the flow deflect away from the duct's centerline. The behavior of the deflected flow is significantly altered downstream of the stricture. Such bile flow behavior as deceleration and recirculation may lead to cholestasis. Stricture alters bile flow in the biliary tract, causing changes in biliary hydrodynamic indexes, which could potentially serve as an omen for gallstone formation and other related diseases. The consideration of the bile duct stricture could lead to better patient stratification.

Evaluation of air impingement for dry-cleaning nonfat dry milk residues on a stainless-steel surface

The use of air jet impingement to remove residues from surfaces in food manufacturing operations offers an alternative to the use of water and liquid cleaning agents. During this investigation, air impingement was used to remove nonfat dry milk (NFDM) residues from a stainless-steel surface. The influence of the water activity (aw ) of the residue, the time after the residue reached an equilibrium water activity, and the thickness of residue at the time of removal from the surface have been investigated. All three factors had a significant effect on the time for removal. An increase in the water activity, the time at equilibrium, the sample thickness, or a combination of all three resulted in an increase in the time required to remove the deposits. Visible changes in the structure of deposits were observed as NFDM samples equilibrated to water activities above 0.43. NFDM residues with water activities less than 0.33 were removed within 1 s of using air impingement regardless of wall shear stress. When the water activities were greater than 0.50, the thickness of deposit was greater than 1 mm, and the time after reaching an equilibrium water activity was over 7 days, more than 5 min of air impingement with wall shear stress over 9.48 Pa was required to remove the residue. The results from these experiments indicated that air impingement has the potential to provide effective cleaning in manufacturing facilities for low-moisture foods. PRACTICAL APPLICATION: The introduction of water in low-moisture food environments is often undesirable due to the possibility of pathogenic microorganism growth. The normal cleaning operations in the food industry use water as a cleaning agent. This study evaluates the application of air impingement technology as a dry-cleaning method.

Impact of hematocrit on pulsatile blood flow in stenosed arteries: a computational study in healthy, diabetic, and anemic models

In this study, researchers aim to enhance the realism of circulatory system simulations, focusing on factors affecting flow variations, particularly in stenotic arteries of individuals with altered hematocrit levels. Through extensive data collection and varied conditions, the goal is to attain more precise and valid results. The study conducts approximate simulations to comprehensively describe the dynamic motion of pulsatile flow. Different values of inlet velocity (UDF) are introduced, considering potential arterial distortion or occlusion due to plaque deposition, along with variations in hematocrit (Hct) levels commonly observed in patients. Three distinct types of pulsatile blood flow, corresponding to diabetes (Hct 65%), healthy (Hct 45%), and anemia (Hct 25%), are studied and compared. The research illuminates that stenosis in arteries with varying hematocrit levels significantly impacts hydrodynamic features, potentially predisposing individuals to cardiovascular diseases. Through meticulous analysis, several conclusions about hemodynamic characteristics are drawn. It is observed that both velocity and wall shear stress exhibit variation along the affected artery, influenced by stenosis and changes in hematocrit levels. Notably, the highest influence on velocity and wall shear stress is observed with Hct 65%, compared to Hct 45% and Hct 25% at the moment of stenosis. These findings hold substantial practical implications for the field of cardiovascular health, providing valuable insights into blood flow behavior in stenotic arteries with diverse hematocrit levels. Ultimately, this research contributes to more effective clinical interventions.

A computational model-based study on the feasibility of predicting post-splenectomy thrombosis using hemodynamic metrics

For portal hypertensive patients with splenomegaly and hypersplenism, splenectomy is an effective surgery to relieve the complications. However, patients who have undergone splenectomy often suffer from portal venous system thrombosis, a sequela that requires prophylaxis and timely treatment to avoid deterioration and death. The aim of this study is to investigate the feasibility of predicting post-splenectomy thrombosis using hemodynamic metrics based on computational models. First, 15 portal hypertensive patients who had undergone splenectomy were enrolled, and their preoperative clinical data and postoperative follow-up results were collected. Next, computational models of the portal venous system were constructed based on the preoperative computed tomography angiography images and ultrasound-measured flow velocities. On this basis, splenectomy was mimicked and the postoperative area of low wall shear stress (ALWSS) was simulated for each patient-specific model. Finally, model-simulated ALWSS was statistically compared with the patient follow-up results to investigate the feasibility of predicting post-splenectomy thrombosis using hemodynamic metrics. Results showed that ALWSS could predict the occurrence of post-splenectomy thrombosis with the area under the receiver operating characteristic curve (AUC) equal to 0.75. Moreover, statistical analysis implied that the diameter of the splenic vein is positively correlated with ALWSS (r = 0.883, p < 0.0001), and the anatomical structures of the portal venous system also influence the ALWSS. These findings demonstrated that the computational model-based hemodynamic metric ALWSS, which is associated with the anatomorphological features of the portal venous system, is capable of predicting the occurrence of post-splenectomy thrombosis, promoting better prophylaxis and postoperative management for portal hypertensive patients receiving splenectomy.

Hemodynamic Parameters in the Parent Arteries of Unruptured Intracranial Aneurysms Depend on Aneurysm Size and Are Different Compared to Contralateral Arteries: A 7 Tesla 4D Flow MRI Study

BACKGROUND

Different Circle of Willis (CoW) variants have variable prevalences of aneurysm development, but the hemodynamic variation along the CoW and its relation to presence and size of unruptured intracranial aneurysms (UIAs) are not well known.

PURPOSE

Gain insight into hemodynamic imaging markers of the CoW for UIA development by comparing these outcomes to the corresponding contralateral artery without an UIA using 4D flow magnetic resonance imaging (MRI).

STUDY TYPE

Retrospective, cross-sectional study.

SUBJECTS

Thirty-eight patients with an UIA, whereby 27 were women and a mean age of 62 years old.

FIELD STRENGTH/SEQUENCE

Four-dimensional phase-contrast (PC) MRI with a 3D time-resolved velocity encoded gradient echo sequence at 7 T.

ASSESSMENT

Hemodynamic parameters (blood flow, velocity pulsatility index [vPI], mean velocity, distensibility, and wall shear stress [peak systolic (WSSMAX ), and time-averaged (WSSMEAN )]) in the parent artery of the UIA were compared to the corresponding contralateral artery without an UIA and were related to UIA size.

STATISTICAL TESTS

Paired t-tests and Pearson Correlation tests. The threshold for statistical significance was P < 0.05 (two-tailed).

RESULTS

Blood flow, mean velocity, WSSMAX , and WSSMEAN were significantly higher, while vPI was lower, in the parent artery relative to contralateral artery. The WSSMAX of the parent artery significantly increased linearly while the WSSMEAN decreased linearly with increasing UIA size.

CONCLUSIONS

Hemodynamic parameters and WSS differ between parent vessels of UIAs and corresponding contralateral vessels. WSS correlates with UIA size, supporting a potential hemodynamic role in aneurysm pathology.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 2.

Treatment of aneurysmal artery with PED: A case report

RATIONALE

Pipeline embolization devices are a new treatment for intracranial aneurysms, especially irregular, giant aneurysms.

PATIENT CONCERNS

A 48-year-old female patient presented with a headache in the frontal part for 3 weeks.

DIAGNOSES

Cerebral computed tomographic angiography and magnetic resonance angiography showed a saccular aneurysm in the right internal carotid artery that wrapped the parent artery. Digital subtraction angiography provided images with 1 large saccular internal carotid aneurysm.

INTERVENTIONS

The patient was treated by flow diverter stent and coil embolization and dual antiplatelet therapy with aspirin and ticagrelor in the ICU and was discharged after 10 days without complications.

OUTCOMES

One year after interventional therapy, repeated digital subtraction angiography showed no recurrence of aneurysm and embolization well, and the patient reported improvement in neurological symptoms.

LESSONS

Aneurysmal parent artery is a rare phenomenon. The combination of flow diverter stents and coil embolization to treat cases with large saccular aneurysms has important clinical significance and this may provide a reference for clinical treatment of aneurysmal parent artery.

Four-Dimensional Flow MRI-Derived Hemodynamics in Abdominal Aortic Aneurysms: Reproducibility and Associations With Diameter, Intraluminal Thrombus Volume, and Vorticity

BACKGROUND

Maximum diameter measurements are used to assess the rupture risk of abdominal aortic aneurysms (AAAs); however, these are not precise enough to predict all ruptures. Four-dimensional (4D) flow MRI-derived parameters provide additional information by visualizing hemodynamics in AAAs but merit further investigation before they are clinically applicable.

PURPOSE

To assess the reproducibility of 4D flow MRI-derived hemodynamics, to investigate possible correlations with lumen and maximum diameter, and to explore potential relationships with vorticity and aneurysm growth.

STUDY TYPE

Prospective single-arm study.

POPULATION

A total of 22 (71.5 ± 6.1 years, 20 male) asymptomatic AAA patients with a maximum diameter of at least 30 mm.

FIELD STRENGTH/SEQUENCE

A 3.0 T/Free-breathing 4D flow MRI phase-contrast acquisition with retrospective ECG-gating.

ASSESSMENT

Patients underwent two consecutive 4D flow MRI scans 1-week apart. Aortic volumes were segmented from time-averaged phase contrast magnetic resonance angiographies. Reproducibility was assessed by voxelwise analysis after registration. Mean flow velocity, mean wall shear stress (WSS), mean lumen diameter, and qualitative vorticity scores were assessed. In addition, Dixon MRI and retrospective surveillance data were used to study maximum diameter (including thrombus), intraluminal thrombus volume (ILT), and growth rate.

STATISTICAL TESTS

For reproducibility assessment, Bland-Altman analyses, Pearson correlation, Spearman's correlation, and orthogonal regression were conducted. Potential correlations between hemodynamics and vorticity scores were assessed using linear regression. P < 0.05 was considered statistically significant.

RESULTS

Test-retest median Pearson correlation coefficients for flow velocity and WSS were 0.85 (IQR = 0.08) m/sec and 0.82 (IQR = 0.10) Pa, respectively. Mean WSS significantly correlated with mean flow velocity (R = 0.75) and inversely correlated with mean lumen diameter (R = -0.73). No significant associations were found between 4D flow MRI-derived hemodynamic parameters and maximum diameter (flow velocity: P = 0.98, WSS: P = 0.22).

DATA CONCLUSION

A 4D flow MRI is robust for assessing the hemodynamics within AAAs. No correlations were found between hemodynamic parameters and maximum diameter, ILT volume and growth rate.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 2.

Numerical investigation of unruptured middle cerebral artery bifurcation aneurysms: influence of aspect ratio

An aneurysm is a disease condition, which is due to the pathological weakening of an arterial wall. These aneurysms are often found in various branch points and bifurcations of an artery in the cerebral circulation. Most aneurysms come to medical attention, either due to brain hemorrhages caused by rupture or found unruptured. To consider surgically invasive treatment modalities, clinicians need scientific methods such as, hemodynamic analysis to assess rupture risk. The arterial wall loses its structural integrity when wall shear stress (WSS) and other hemodynamic parameters exceed a certain threshold. In the present study, numerical simulations are carried out for unruptured middle cerebral artery (MCA) aneurysms. Three distinct representative sizes are chosen from a larger patient pool of 26 MCA aneurysms. Logically, these aneurysms represent three growth stages of any patient with similar anatomical structure. Simulations are performed to compare the three growth phases (with different aspect ratios) of an aneurysm and correlate their hemodynamic parameters. Simulations with patient specific boundary conditions reveal that, aneurysms with a higher aspect ratio (AR) correspond to an attendant decrease in both time-averaged wall shear stress (TAWSS) and spatial wall shear stress gradients (WSSG). Smaller MCAs were observed to have higher positive wall shear stress divergence (WSSD), exemplifying the tensile nature of arterial wall stretching. Present study identifies positive wall shear stress divergence (PWSSD) to be a potential biomarker for evaluating the growth of an aneurysm.

First experiences of ultrasound vector flow imaging at the femoropopliteal artery in peripheral arterial disease

Background: The femoropopliteal artery (FPA) plays a central role in diagnosing and treating peripheral arterial disease (PAD). FPA lesions are the most frequent cause of intermittent claudication, and no other artery of the lower extremities is recanalised more frequently. Generally, ultrasound is the primary imaging tool in PAD, particularly FPA. With the development of high-frame-rate ultrasound technology in addition to traditional ultrasound modes, vector flow imaging (VFI) has provided deeper haemodynamic insights when used in the carotid artery. Here, we report the use of VFI at the FPA level in routine PAD examinations. Patients and methods: In this single-centre prospective study, we evaluated consecutive patients with PAD using B-mode imaging, colour Doppler, pulsed wave Doppler (PW) and vector flow. Hemodynamic parameters at predefined locations at the carotid artery and FPA were compared. Results: Qualitatively adequate VFI at all sites was possible in 76% of the patients with PAD. With decreasing volume flow from the common carotid artery to the internal carotid artery and from the common femoral artery via the superficial femoral artery to the popliteal artery, the correlation between VFI- and PW-derived-volume flow was high at every site. Based on different techniques, the VFI-derived values were significantly lower than the PW-derived values. The mean wall shear stress was significantly lower at all femoropopliteal sites than at the carotid sites, whereas the oscillatory shear index at the femoral site was higher than that at the carotid sites rather than at the popliteal location. Conclusions: Our findings suggest that vector flow data acquisition in the FPA is feasible in most patients with PAD. Therefore, with knowledge of the method and its limitations, VFI provides haemodynamic information beyond traditional ultrasound techniques and is a promising new tool for flow analysis in PAD.

Orthopedic Scaffolds: Evaluation of Structural Strength and Permeability of Fluid Flow via an Open Cell Neovius Structure for Bone Tissue Engineering

The ability of bone to regenerate itself through mechanobiological responses is its dynamic property. Mechanical cues from a neighboring environment produce the structural strain to promote blood flow and bone marrow mobility that in turn aids the bone regeneration process. Occurrences of these phenomena are crucial for the success of metallic scaffolds implanted in the host bone tissue. Thus, permeability and fluid flow-induced wall shear stress (WSS) are two parameters that directly influence cell bioactivities inside a scaffold and are crucial for effective bone tissue regeneration. Given that the scaffolds shall be implanted in the body, permeability assessment was carried out using non-Newtonian fluid. In this work, the triply periodic minimal surface scaffolds with Neovius architectures were fabricated by using selective laser melting technology. The estimation of fluid flow was carried out using computational fluid dynamics (CFD) analysis with a non-Newtonian blood fluid model. Further, the structural strength of various open cell Neovius lattices was evaluated using a static compression test, and in vitro cell culture using Alamar blue assay was evaluated. Results revealed that the values of intrinsic blood flow permeability of the three-dimensional (3D)-printed open cell porous scaffold with Neovius architecture were of the same order of magnitude as those of human bone, ranging from 0.0025 × 10-9 to 0.0152 × 10-9 m2. The structural elastic modulus and compressive strength of NOCL40, NOCL50, and NOCL60 lattices range from 3.27 to 3.71 GPa and 194 to 205 MPa, respectively. All of the values are comparable to the human bone, thus making these lattices a suitable alternative for orthopedic applications.

Influence of Aortic Valve Stenosis and Wall Shear Stress on Platelets Function

Aortic valve stenosis (AS) is a common heart valve disease in the elderly population, and its pathogenesis remains an interesting area of research. The degeneration of the aortic valve leaflets gradually progresses to valve sclerosis. The advanced phase is marked by the presence of extracellular fibrosis and calcification. Turbulent, accelerated blood flow generated by the stenotic valve causes excessive damage to the aortic wall. Elevated shear stress due to AS leads to the degradation of high-molecular weight multimers of von Willebrand factor, which may involve bleeding in the mucosal tissues. Conversely, elevated shear stress has been associated with the release of thrombin and the activation of platelets, even in individuals with acquired von Willebrand syndrome. Moreover, turbulent blood flow in the aorta may activate the endothelium and promote platelet adhesion and activation on the aortic valve surface. Platelets release a wide range of mediators, including lysophosphatidic acid, which have pro-osteogenic effects in AS. All of these interactions result in blood coagulation, fibrinolysis, and the hemostatic process. This review summarizes the current knowledge on high shear stress-induced hemostatic disorders, the influence of AS on platelets and antiplatelet therapy.

Cerebral artery signal intensity gradient from Time-of-Flight Magnetic Resonance Angiography and clinical outcome in lenticulostriate infarction: a retrospective cohort study

Purpose

Lenticulostriate infarction requires further research of arterial hemodynamic factors, as the disease is diagnosed in the absence of major arterial stenosis or cardioembolism.

Methods

In this multicenter retrospective cohort study, we included patients who were hospitalized for lenticulostriate infarction from January 2015 to March 2021 at three stroke centers in South Korea. We obtained hemodynamic information on cerebral arteries using signal intensity gradient (SIG), an in-vivo approximated wall shear stress (WSS) derived from Time-of-Flight Magnetic Resonance Angiography (TOF-MRA). A favorable outcome was defined as a modified Rankin Scale of 0 to 2 at hospital discharge.

Results

A total of 294 patients were included, of whom 146 (49.7%) had an unfavorable outcome. The unfavorable outcome group showed significantly lower SIG in both middle cerebral arteries (MCAs) than the favorable group (5.2 ± 1.2 SI/mm vs. 5.9 ± 1.2, p < 0.001), and similar findings were observed in other cerebral arteries. The SIGs in both MCAs were independently associated with favorable outcome, with an odds ratio of 1.42 (95% confidence interval, 1.11-1.80; p = 0.005) for the right MCA and 1.49 (95% CI, 1.15-1.93; p = 0.003) for the left MCA, after adjusting for potential confounders. Similar findings were observed in other cerebral artery SIGs.

Conclusion

Cerebral artery SIG from TOF-MRA was significantly associated with short-term functional outcomes in patients with lenticulostriate infarction. Further studies are needed to investigate the temporal relationships of SIG in patients with cerebral infarction.

In-silico enhanced animal study of pulmonary artery pressure sensors: assessing hemodynamics using computational fluid dynamics

To assess whether in-silico models can be used to predict the risk of thrombus formation in pulmonary artery pressure sensors (PAPS), a chronic animal study using pigs was conducted. Computed tomography (CT) data was acquired before and immediately after implantation, as well as one and three months after the implantation. Devices were implanted into 10 pigs, each one in the left and right pulmonary artery (PA), to reduce the required number of animal experiments. The implantation procedure aimed at facilitating optimal and non-optimal positioning of the devices to increase chances of thrombus formation. Eight devices were positioned non-optimally. Three devices were positioned in the main PA instead of the left and right PA. Pre-interventional PA geometries were reconstructed from the respective CT images, and the devices were virtually implanted at the exact sites and orientations indicated by the follow-up CT after one month. Transient intra-arterial hemodynamics were calculated using computational fluid dynamics. Volume flow rates were modelled specifically matching the animals body weights. Wall shear stresses (WSS) and oscillatory shear indices (OSI) before and after device implantation were compared. Simulations revealed no relevant changes in any investigated hemodynamic parameters due to device implantation. Even in cases, where devices were implanted in a non-optimal manner, no marked differences in hemodynamic parameters compared to devices implanted in an optimal position were found. Before implantation time and surface-averaged WSS was 2.35 ± 0.47  Pa, whereas OSI was 0.08 ± 0.17 , respectively. Areas affected by low WSS magnitudes were 2.5 ± 2.7  cm2 , whereas the areas affected by high OSI were 18.1 ± 6.3  cm2 . After device implantation, WSS and OSI were 2.45 ± 0.49  Pa and 0.08 ± 0.16 , respectively. Surface areas affected by low WSS and high OSI were 2.9 ± 2.7  cm2 , and 18.4 ± 6.1  cm2 , respectively. This in-silico study indicates that no clinically relevant differences in intra-arterial hemodynamics are occurring after device implantation, even at non-optimal positioning of the sensor. Simultaneously, no embolic events were observed, suggesting that the risk for thrombus formation after device implantation is low and independent of the sensor position.

The Stress Phase Angle Measurement Using Spectral Domain Optical Coherence Tomography

The stress phase angle (SPA), defined as the temporal phase angle between circumferential stress (CS) in the arterial wall and wall shear stress (WSS), is utilized to investigate the interactions between CS and WSS. SPA serves as an important parameter for the early diagnosis of cardiovascular disease. In this study, we proposed a novel method for measuring SPA using spectral domain optical coherence tomography (SD-OCT). The multi-M-mode scan strategy is adopted for interference spectrum acquisition. The phases of CS and WSS are extracted from the corresponding structural and flow velocity images of SD-OCT. The method is validated by measuring SPA in the outflow tract (OFT) of chick embryonic hearts and the common carotid artery of mice. To the best of our knowledge, this is the first time that OCT has been used for SPA measurement.

The site-specific distribution of atheromatous plaques in the coronary arteries

The etiology of atherosclerosis is still unknown, but there are several hypotheses trying to explain this complex disease. Most consider atherosclerosis as a cholesterol storage disease. However, hypercholesterolemia is not a cause but a risk factor. Besides, like other well-known systemic risk factors, it does not explain the uneven distribution of atheromatous plaques in the vasculature. Atherosclerotic lesions develop mainly at vulnerable "risk points" of the arterial wall such as curvatures and near side branches, and predominantly in the left anterior descending (LAD), while the left circumflex (LCx) artery is relatively spared. Furthermore, atheromatous plaques are present mainly in the proximal segments in the LAD and LCx, in contrast to the right coronary artery (RCA), where plaques are more evenly distributed. The hemodynamic theory explains to some extent the distribution of atherosclerotic lesions and considers atherosclerosis as a reactive biological response of endothelial cells to wall shear stress. In this review, we discuss the interplay of concentration of low-density lipoproteins at the luminal surface and local hemodynamic forces (disturbed flow) that reduce wall shear stress in the process of plaque formation. Moreover, we present the distribution of atheromatous plaques in the coronary arteries in autopsy studies and imaging methods such as cardiac computed tomography angiography and invasive coronary angiography.

The influencing mechanism of iliac vein stent implantation for hemodynamics at the bifurcation

In the intervention with stent implantation for iliac vein compression syndrome (IVCS), it remains unclear about the influencing mechanism of the structure and implantation position of the stent for the hemodynamics of the affected site. In this paper, an iliac vein model was established. Besides, the computational fluid dynamics (CFD) was utilized to analyze the time-averaged wall shear stress (TAWSS) and oscillatory shear index (OSI) in a sine period after stent implantation based on the three different implantation positions of two iliac vein stents (the left branch outlet, contralateral disturbed flow and main iliac vein). The influence of the structure and implantation position of the stent on blood flow was revealed. These findings were verified by the particle image velocimetry (PIV) experiment. The results indicated that the maximum blood flow velocity of the iliac vein decreased after the stent implantation. Among the three positions, the influence of stent implantation on the iliac vein blood flow was the least when the stent implantation was performed at the left branch outlet; the influence of stent implantation on the iliac vein blood flow was the greatest when the stent implantation was performed at the contralateral disturbed flow. Moreover, there was little influence of Venastent implantation on the blood flow. These results could provide a scientific foundation for implantation treatment and stent design related to IVCS.

Computational fluid dynamics as supporting technology for coronary artery disease diagnosis and treatment: an international survey

Background

Computational fluid dynamics (CFD) is emerging as an effective technology able to improve procedural outcomes and enhance clinical decision-making in patients with coronary artery disease (CAD). The present study aims to assess the state of knowledge, use and clinical acceptability of CFD in the diagnosis and treatment of CAD.

Methods

We realized a 20-questions international, anonymous, cross-sectional survey to cardiologists to test their knowledge and confidence on CFD as a technology applied to patients suffering from CAD. Responses were recorded between May 18, 2022, and June 12, 2022.

Results

A total of 466 interventional cardiologists (mean age 48.4 ± 8.3 years, males 362), from 42 different countries completed the survey, for a response rate of 45.9%. Of these, 66.6% declared to be familiar with the term CFD, especially for optimization of existing interventional techniques (16.1%) and assessment of hemodynamic quantities related with CAD (13.7%). About 30% of respondents correctly answered to the questions exploring their knowledge on the pathophysiological role of some CFD-derived quantities such as wall shear stress and helical flow in coronary arteries. Among respondents, 85.9% would consider patient-specific CFD-based analysis in daily interventional practice while 94.2% declared to be interested in receiving a brief foundation course on the basic CFD principles. Finally, 87.7% of respondents declared to be interested in a cath-lab software able to conduct affordable CFD-based analyses at the point-of-care.

Conclusions

Interventional cardiologists reported to be profoundly interested in adopting CFD simulations as a technology supporting decision making in the treatment of CAD in daily practice.

Effects of size and shape of the side holes of a double J stent on the ureter fluid flow after stenosis

The effect of side holes morphology changes in double J stent (DJS) on encrustation was analyzed using computational fluid dynamics (CFD). We analyzed DJS side holes with inner diameter of 1 mm and outer diameters of 1 (type A), 1.2 (type B) and 1.4 (type C) mm, respectively. Concentric stenosis with three intraureteral degree (0%, 12%, and 88%) was analyzed. The flow rate, shear stress and wall shear stress (WSS) distribution were investigated. Urine flow through SH1 before the ureteropelvic junction (UPJ) differed based on the ureteral stenosis degree. The sum of flow rates through the SHs increased with diameter. In the stented ureter with 12% stenosis, the flow rate through SH1 approximately doubled than that without ureteral stenosis, and the flow rate through SH1 was maximal for the type 'C' stent in both 12% and 88% ureteral stenosis. The mean shear stress in the SHs increased with the degree of stenosis. The WSS around the SHs was higher for type 'C' than types A and B. From the flow rates and shear stresses in and around the SHs, the larger SH diameter of the DJS from the UPJ to mid-ureter is expected to induce encrustation reduction, especially in patients with urinary lithiasis.

Angiogenic Microvascular Wall Shear Stress Patterns Revealed Through Three-dimensional Red Blood Cell Resolved Modeling

The wall shear stress (WSS) exerted by blood flowing through microvascular capillaries is an established driver of new blood vessel growth, or angiogenesis. Such adaptations are central to many physiological processes in both health and disease, yet three-dimensional (3D) WSS characteristics in real angiogenic microvascular networks are largely unknown. This marks a major knowledge gap because angiogenesis, naturally, is a 3D process. To advance current understanding, we model 3D red blood cells (RBCs) flowing through rat angiogenic microvascular networks using state-of-the-art simulation. The high-resolution fluid dynamics reveal 3D WSS patterns occurring at sub-endothelial cell (EC) scales that derive from distinct angiogenic morphologies, including microvascular loops and vessel tortuosity. We identify the existence of WSS hot and cold spots caused by angiogenic surface shapes and RBCs, and notably enhancement of low WSS regions by RBCs. Spatiotemporal characteristics further reveal how fluctuations follow timescales of RBC "footprints." Altogether, this work provides a new conceptual framework for understanding how shear stress might regulate EC dynamics in vivo.

Biomechanical determinants of endothelial permeability assessed in standard and modified hollow-fibre bioreactors

Effects of mechanical stress on the permeability of vascular endothelium are important to normal physiology and in the development of atherosclerosis. Here we elucidate novel effects using commercially available and modified hollow-fibre bioreactors, in which endothelial cells form confluent monolayers lining plastic capillaries with porous walls, contained in a cartridge. The capillaries were perfused with a near-aortic waveform, and permeability was assessed by the movement of rhodamine-labelled albumin from the intracapillary to the extracapillary space. Permeability was increased by acute application of shear stress and decreased by chronic shear stress compared with a static control: this has previously been shown only for multidirectional flows. Increasing viscosity reduced permeability under both acute and chronic shear; since shear rate remained unchanged, these effects resulted from altered shear stress. Reducing pulsatility increased permeability, contrary to the widely held assumption that flow which is highly oscillatory causes endothelial dysfunction. Chronic convection across the monolayer increased effective permeability more than could be explained by the addition of advective transport, contrary to results from previous acute experiments. The off-the-shelf and modified bioreactors provide an excellent tool for investigating the biomechanics of endothelial permeability and have revealed novel effects of flow duration, viscosity, pulsatility and transmural flow.

Effect of Particle Migration on the Stress Field in Microfluidic Flows of Blood Analog Fluids at High Reynolds Numbers

In the present paper, we investigate how the reductions in shear stresses and pressure losses in microfluidic gaps are directly linked to the local characteristics of cell-free layers (CFLs) at channel Reynolds numbers relevant to ventricular assist device (VAD) applications. For this, detailed studies of local particle distributions of a particulate blood analog fluid are combined with wall shear stress and pressure loss measurements in two complementary set-ups with identical flow geometry, bulk Reynolds numbers and particle Reynolds numbers. For all investigated particle volume fractions of up to 5%, reductions in the stress and pressure loss were measured in comparison to a flow of an equivalent homogeneous fluid (without particles). We could explain this due to the formation of a CFL ranging from 10 to 20 μm. Variations in the channel Reynolds number between Re = 50 and 150 did not lead to measurable changes in CFL heights or stress reductions for all investigated particle volume fractions. These measurements were used to describe the complete chain of how CFL formation leads to a stress reduction, which reduces the apparent viscosity of the suspension and results in the Fåhræus-Lindqvist effect. This chain of causes was investigated for the first time for flows with high Reynolds numbers (Re∼100), representing a flow regime which can be found in the narrow gaps of a VAD.

Numerical investigation of the effect of carotid bifurcation stenosis degree on pulsatility characteristics

Arterial bifurcations are regions that are susceptible to hemodynamic effects and thrombus formation. In the current study, the hemodynamic effects of a simplified 3D model of an arterial bifurcation were simulated using the commercial computational fluid dynamics software FLUENT. The non-Newtonian properties of blood were modeled using the Carreau model, and the pulsation dynamics and heat transfer characteristics of blood at different degrees of stenosis in the arterial bifurcation were analyzed. The results indicate that arterial stenosis caused by a thrombus when the pulsation velocity reaches its peak has an essential impact on blood transport. The stenosis of the bifurcation increases the peak pulsatile flow pressure drop, and each 0.5 mm stenosis of the arterial bifurcation increases the mean wall shear stress of the bifurcated segment by approximately 0.25 Pa. From the heat transfer perspective, arterial stenosis has little effect on the heat transfer coefficient. The heat transfer coefficient measured inside the bifurcation is much larger than that measured outside the bifurcation. The stenosis of the arterial bifurcation causes an increase in the mean velocity of the arterial cross-section, and the volume-averaged absolute vorticity is introduced to quantify the secondary flow effect during the pulsation cycle, where the arterial stenosis causes an increase in the mean absolute vorticity at pulsation velocity and accelerates the decay of the vorticity at uniform velocity. In this paper, the hemodynamics of carotid bifurcation pulsation is analyzed in conjunction with flow field properties to reveal the flow field dynamics factors and heat transfer characteristics of local stenosis of the carotid bifurcation and to conduct an exploratory study for the diagnosis and treatment of carotid bifurcation thrombosis.

Transport of nitrite from large arteries modulates regional blood flow during stress and exercise

Background

Acute cardiovascular stress increases systemic wall shear stress (WSS)-a frictional force exerted by the flow of blood on vessel walls-which raises plasma nitrite concentration due to enhanced endothelial nitric oxide synthase (eNOS) activity. Upstream eNOS inhibition modulates distal perfusion, and autonomic stress increases both the consumption and vasodilatory effects of endogenous nitrite. Plasma nitrite maintains vascular homeostasis during exercise and disruption of nitrite bioavailability can lead to intermittent claudication.

Hypothesis

During acute cardiovascular stress or strenuous exercise, we hypothesize enhanced production of nitric oxide (NO) by vascular endothelial cells raises nitrite concentrations in near-wall layers of flowing blood, resulting in cumulative NO concentrations in downstream arterioles sufficient for vasodilation.

Confirmation and implications

Utilizing a multiscale model of nitrite transport in bifurcating arteries, we tested the hypothesis for femoral artery flow under resting and exercised states of cardiovascular stress. Results indicate intravascular transport of nitrite from upstream endothelium could result in vasodilator-active levels of nitrite in downstream resistance vessels. The hypothesis could be confirmed utilizing artery-on-a-chip technology to measure NO production rates directly and help validate numerical model predictions. Further characterization of this mechanism may improve our understanding of symptomatic peripheral artery occlusive disease and exercise physiology.

Fluid-structure interaction study for biomechanics and risk factors in Stanford type A aortic dissection

Aortic dissection is a life-threatening condition with a rising prevalence in the elderly population, possibly as a consequence of the increasing population life expectancy. Untreated aortic dissection can lead to myocardial infarction, aortic branch malperfusion or occlusion, rupture, aneurysm formation and death. This study aims to assess the potential of a biomechanical model in predicting the risks of a non-dilated thoracic aorta with Stanford type A dissection. To achieve this, a fully coupled fluid-structure interaction model was developed under realistic blood flow conditions. This model of the aorta was developed by considering three-dimensional artery geometry, multiple artery layers, hyperelastic artery wall, in vivo-based physiological time-varying blood velocity profiles, and non-Newtonian blood behaviours. The results demonstrate that in a thoracic aorta with Stanford type A dissection, the wall shear stress (WSS) is significantly low in the ascending aorta and false lumen, leading to potential aortic dilation and thrombus formation. The results also reveal that the WSS is highly related to blood flow patterns. The aortic arch region near the brachiocephalic and left common carotid artery is prone to rupture, showing a good agreement with the clinical reports. The results have been translated into their potential clinical relevance by revealing the role of the stress state, WSS and flow characteristics as the main parameters affecting lesion progression, including rupture and aneurysm. The developed model can be tailored for patient-specific studies and utilised as a predictive tool to estimate aneurysm growth and initiation of wall rupture inside the human thoracic aorta.

Development and Evaluation of a Multifrequency Ultrafast Doppler Spectral Analysis (MFUDSA) Algorithm for Wall Shear Stress Measurement: A Simulation and In Vitro Study

Cardiovascular pathology is the leading cause of death and disability in the Western world, and current diagnostic testing usually evaluates the anatomy of the vessel to determine if the vessel contains blockages and plaques. However, there is a growing school of thought that other measures, such as wall shear stress, provide more useful information for earlier diagnosis and prediction of atherosclerotic related disease compared to pulsed-wave Doppler ultrasound, magnetic resonance angiography, or computed tomography angiography. A novel algorithm for quantifying wall shear stress (WSS) in atherosclerotic plaque using diagnostic ultrasound imaging, called Multifrequency ultrafast Doppler spectral analysis (MFUDSA), is presented. The development of this algorithm is presented, in addition to its optimisation using simulation studies and in-vitro experiments with flow phantoms approximating the early stages of cardiovascular disease. The presented algorithm is compared with commonly used WSS assessment methods, such as standard PW Doppler, Ultrafast Doppler, and Parabolic Doppler, as well as plane-wave Doppler. Compared to an equivalent processing architecture with one-dimensional Fourier analysis, the MFUDSA algorithm provided an increase in signal-to-noise ratio (SNR) by a factor of 4-8 and an increase in velocity resolution by a factor of 1.10-1.35. The results indicated that MFUDSA outperformed the others, with significant differences detected between the typical WSS values of moderate disease progression (p = 0.003) and severe disease progression (p = 0.001). The algorithm demonstrated an improved performance for the assessment of WSS and has potential to provide an earlier diagnosis of cardiovascular disease than current techniques allow.

Trigeminal neuralgia caused by a persistent primitive trigeminal artery: preoperative three-dimensional multifusion imaging and computational fluid dynamics analysis. Illustrative case

BACKGROUND

Trigeminal neuralgia (TN) is caused by trigeminal nerve compression by colliding vessels. Preoperative three-dimensional (3D) multifusion images are useful for surgical simulations. Moreover, computational fluid dynamics (CFD) analysis of colliding vessels may be useful for hemodynamic evaluation at the site of neurovascular contact (NVC).

OBSERVATIONS

A 71-year-old woman had TN due to compression of the trigeminal nerve by the superior cerebellar artery (SCA) fused with the persistent primitive trigeminal artery (PTA). Preoperative 3D multifusion simulation images of silent magnetic resonance (MR) angiography and MR cisternography depicted the NVC, including the trigeminal nerve, SCA, and PTA. CFD analysis revealed the hemodynamic condition of the NVC, including the SCA and PTA. The wall shear stress magnitude (WSSm) at the NVC showed a local elevation due to flow confluence from the SCA and PTA. High WSSm was observed in the NVC.

LESSONS

Preoperative simulation images of MR angiography and MR cisternography may depict the NVC. CFD analysis can provide the hemodynamic condition at the NVC.

Modification of the swirling well cell culture model to alter shear stress metrics

Effects of hemodynamic shear stress on endothelial cells have been extensively investigated using the "swirling well" method, in which cells are cultured in dishes or multiwell plates placed on an orbital shaker. A wave rotates around the well, producing complex patterns of shear. The method allows chronic exposure to flow with high throughput at low cost but has two disadvantages: a number of shear stress characteristics change in a broadly similar way from the center to the edge of the well, and cells at one location in the well may release mediators into the medium that affect the behavior of cells at other locations, exposed to different shears. These properties make it challenging to correlate cell properties with shear. The present study investigated simple alterations to ameliorate these issues. Flows were obtained by numerical simulation. Increasing the volume of fluid in the well-altered dimensional but not dimensionless shear metrics. Adding a central cylinder to the base of the well-forced fluid to flow in a square toroidal channel and reduced multidirectionality. Conversely, suspending a cylinder above the base of the well made the flow highly multidirectional. Increasing viscosity in the latter model increased the magnitude of dimensional but not dimensionless metrics. Finally, tilting the well changed the patterns of different wall shear stress metrics in different ways. Collectively, these methods allow similar flows over most of the cells cultured and/or allow the separation of different shear metrics. A combination of the methods overcomes the limitations of the baseline model.

Computational Fluid Dynamics (CFD) Model for Analysing the Role of Shear Stress in Angiogenesis in Rheumatoid Arthritis

Rheumatoid arthritis (RA) is an autoimmune disease characterised by an attack on healthy cells in the joints. Blood flow and wall shear stress are crucial in angiogenesis, contributing to RA's pathogenesis. Vascular endothelial growth factor (VEGF) regulates angiogenesis, and shear stress is a surrogate for VEGF in this study. Our objective was to determine how shear stress correlates with the location of new blood vessels and RA progression. To this end, two models were developed using computational fluid dynamics (CFD). The first model added new blood vessels based on shear stress thresholds, while the second model examined the entire blood vessel network. All the geometries were based on a micrograph of RA blood vessels. New blood vessel branches formed in low shear regions (0.840-1.260 Pa). This wall-shear-stress overlap region at the junctions was evident in all the models. The results were verified quantitatively and qualitatively. Our findings point to a relationship between the development of new blood vessels in RA, the magnitude of wall shear stress and the expression of VEGF.

Feasibility study of combining wall shear stress and elastography to assess the vascular status of carotid artery

INTRODUCTION

At present, early detection of the potential risk of atherosclerosis and prevention is of great significance to reduce the occurrence of stroke.

AIM

This study aims to explore the value of combining the wall shear stress measured by ultrasound vector flow imaging technique and sound touch elastography of common carotid artery in normal adults using the Mindray Resona 7 ultrasound system.

METHODS

Forty volunteers (mean age 39.5 y, 23 females, 17 males) were divided into four groups according to their age. All volunteers underwent ultrasound carotid artery examination, and the values of wall shear stress and elasticity on the posterior wall of the common carotid artery were measured using advanced imaging functions, vector flow imaging technique, and sound touch elastography.

RESULTS

Different cut-off values of wall shear stress were used to investigate the significance between two groups with corresponding sound touch elastography values. It can be seen that the statistical difference could be found when the mean wall shear stress was larger than 1.5 Pa approximately (statistical significance was defined when P < 0.05), and the sound touch elastography value was positively correlated with the wall shear stress value.

CONCLUSION

This study reveals that the combination of wall shear stress and sound touch elastography is an effective and feasible method for assessing carotid artery health. When the mean wall shear stress value is over 1.5 Pa, the corresponding sound touch elastography value increases significantly. The risk of atherosclerosis increases with the stiffness of blood vessel walls.

Association of local solid mechanical, hemodynamic and morphological characteristics with ruptured intracranial aneurysm

The rupture of intracranial aneurysms (IAs) is a complicated phenomenon of which the mechanism is not fully understood. The purpose of this study is to associate local solid mechanical, hemodynamic, and morphological characteristics with rupture regions through statistical means, in an attempt to identify the parameters that are indicative of rupture propensity for IAs. Twenty patient-specific ruptured IA models were reconstructed from digital subtraction angiography (DSA), and applied in the analysis of wall tension, wall shear stress (WSS) and curvature. The precise rupture locations were marked out through intraoperative videos. Pearson correlation analysis was employed to investigate the correlations of these three parameters with patient characteristics and global geometric features. Univariate and multivariate logistic regression analysis were further performed on wall tension, WSS and curvature with regards to rupture and nonrupture regions. Receiver operating characteristic (ROC) analysis defining area under the curve (AUC) was performed on these three parameters. The univariate model of wall tension (AUC, 0.9750), WSS (AUC, 0.9300), curvature (0.8150) and their combined multivariate model (AUC, 0.9875) all present high AUC values. The wall tension, WSS and curvature are acceptable parameters relating to rupture regions. The rupture odd is more sensitive to the wall tension and WSS than curvature. Each logistic model is capable in discriminating ruptures from nonrupture regions, while the multivariate model is the most efficient.

Effect of temperature, wall shear stress, and NaOH concentration on cleaning effectiveness

The objective of this research was to evaluate the effect of various operational parameters and their interaction on cleaning rate of NaOH. The parameters include wall shear stress, temperatures, and NaOH concentrations of cleaning fluid. The higher cleaning effectiveness for proteinaceous deposits was achieved at higher wall shear stress. The wall shear stress of 2.42 Pa removed over 90% of foulant after 10 min of operation, which was significantly higher than the removal of lower wall shear stress (0.84 and 0.39 Pa). Similarly, the cleaning rate increased with increase in temperature and concentration of cleaning solution. The use of cleaning solution (0.05% NaOH) at 65°C provided significantly higher cleaning rate than 25°C. A cleaning solution concentration of 0.5% NaOH provided significantly higher removal of foulant than 0% or 0.05% concentrations. However, analysis on the interaction between temperature and the wall shear stress suggested that the temperature above 45°C and wall shear stress above 0.65 Pa did not provide significant improvement in cleaning efficacy. When a cleaning solution temperature was maintained at 45°C, higher wall shear stress provided more rapid removal of the foulant when the cleaning agent concentration was 0.05% or 0.5%. In water rinse conducted without chemical agents, no much improvements in foulant removal were observed with increase in wall shear stress. The change of activation energy (E_a ) indicates that water rinse was sensitive to temperature change at higher wall shear stress. However, cleaning with cleaning agents was less sensitive to temperature variation compared to the water rinse. PRACTICAL APPLICATION: Considering the intense use of water and chemical compounds during cleaning operation, finding operating condition to ensure acceptable cleanness with minimized input is desired. Temperature, cleaning agent concentration, and wall shear stress are the major components to influence the cleaning efficacy. The result of the study demonstrates the potential for optimization of in-place cleaning by appropriate adjustments of the CIP operating parameters.

Magnetic resonance imaging-based hemodynamic wall shear stress alters aortic wall tissue biomechanics in bicuspid aortic valve patients

OBJECTIVE

In this study we aimed to conclusively determine whether altered aortic biomechanics are associated with wall shear stress (WSS) independent of region of tissue collection. Elevated WSS in the ascending aorta of patients with bicuspid aortic valve has been shown to contribute to local maladaptive aortic remodeling and might alter biomechanics.

METHODS

Preoperative 4-dimensional flow magnetic resonance imaging was performed on 22 patients who underwent prophylactic aortic root and/or ascending aorta replacement. Localized elevated WSS was identified in patients using age-matched healthy atlases (n = 60 controls). Tissue samples (n = 78) were collected and categorized according to WSS (elevated vs normal) and region. Samples were subjected to planar biaxial testing. To fully quantify the nonlinear biomechanical response, the tangential modulus (local stiffness) at a low-stretch (LTM) and high-stretch (HTM) linear region and the onset (TZo) and end stress of the nonlinear transition zone were measured. A linear mixed effect models was implemented to determine statistical relationships.

RESULTS

A higher LTM in the circumferential and axial direction was associated with elevated WSS (P = .007 and P = .018 respectively) independent of collection region. Circumferential TZo and HTM were higher with elevated WSS (P = .024 and P = .003); whereas the collection region was associated with variations in axial TZo (P = .013), circumferential HTM (P = .015), and axial HTM (P = .001).

CONCLUSIONS

This study shows strong evidence that biomechanical changes in the aorta are strongly associated with hemodynamics, and not region of tissue collection for bicuspid valve aortopathy patients. Elevated WSS is associated with tissue behavior at low stretch ranges (ie, LTM and TZo).

4D-flow MRI derived wall shear stress for the risk stratification of bicuspid aortic valve aortopathy: A systematic review

Purpose

Current intervention guidelines for bicuspid aortic valve (BAV) associated ascending aorta (AAo) dilatation are suboptimal predictors of clinical outcome. There is growing interest in identifying better biomarkers such as wall shear stress (WSS) to help risk stratify BAV aortopathy. The aim of the systematic review is to synthesize existing evidence of the relationship between WSS and aortopathy in the BAV population.

Methods

A comprehensive literature search of available major databases was performed in May 2022 to include studies that used four-dimensional flow cardiac magnetic resonance (4D-flow) MRI to quantify WSS in the AAo in adult BAV populations. Summary results and statistical analysis were provided for key numerical results. A narrative summary was provided to assess similarities between studies.

Results

A total of 26 studies that satisfied selection criteria and quality assessment were included in the review. The presence of BAV resulted in significantly elevated WSS magnitude and circumferential WSS, but not axial WSS. The presence of aortic stenosis had additional impact on WSS and flow alterations. BAV phenotypes were associated with different WSS distributions and flow profiles. Altered protein expression in the AAo wall associated with WSS supported the contribution of altered hemodynamics to aortopathy in addition to genetic factors.

Conclusion

WSS has the potential to be a valid biomarker for BAV aortopathy. Future work would benefit from larger study cohorts with longitudinal evaluations to further characterize WSS association with aortopathy, mortality, and morbidities.

Systematic review registration

https://www.crd.york.ac.uk/prospero/display_record.php?ID=CRD42022337077, identifier CRD42022337077.

Adaptive Wall Shear Stress Imaging in Phantoms, Simulations and In Vivo

WSS measurement is challenging since it requires sensitive flow measurements at a distance close to the wall. The aim of this study is to develop an ultrasound imaging technique which combines vector flow imaging with an unsupervised data clustering approach that automatically detects the region close to the wall with optimally linear flow profile, to provide direct and robust WSS estimation. The proposed technique was evaluated in phantoms, mimicking normal and atherosclerotic vessels, and spatially registered Fluid Structure Interaction (FSI) simulations. A relative error of 6.7% and 19.8% was obtained for peak systolic (WSSPS) and end diastolic (WSSED) WSS in the straight phantom, while in the stenotic phantom, a good similarity was found between measured and simulated WSS distribution, with a correlation coefficient, R, of 0.89 and 0.85 for WSSPS and WSSED, respectively. Moreover, the feasibility of the technique to detect pre-clinical atherosclerosis was tested in an atherosclerotic swine model. Six swines were fed atherogenic diet, while their left carotid artery was ligated in order to disturb flow patterns. Ligated arterial segments that were exposed to low WSSPS and WSS characterized by high frequency oscillations at baseline, developed either moderately or highly stenotic plaques (p < 0.05). Finally, feasibility of the technique was demonstrated in normal and atherosclerotic human subjects. Atherosclerotic carotid arteries with low stenosis had lower WSSPS as compared to control subjects (p < 0.01), while in one subject with high stenosis, elevated WSS was found on an arterial segment, which coincided with plaque rupture site, as determined through histological examination.

Correlation between aortic valve protein levels and vector flow mapping of wall shear stress and oscillatory shear index in patients supported with continuous-flow left ventricular assist devices

BACKGROUND

Continuous-flow left ventricular assist devices commonly lead to aortic regurgitation, which results in decreased pump efficiency and worsening heart failure. We hypothesized that non-physiological wall shear stress and oscillatory shear index alter the abundance of structural proteins in aortic valves of left ventricular assist device (LVAD) patients.

METHODS

Doppler images of aortic valves of patients undergoing heart transplants were obtained. Eight patients had been supported with LVADs, whereas 10 were not. Aortic valve tissue was collected and protein levels were analyzed using mass spectrometry. Echocardiographic images were analyzed and wall shear stress and oscillatory shear index were calculated. The relationship between normalized levels of individual proteins and in vivo echocardiographic measurements was evaluated.

RESULTS

Of the 57 proteins of interest, there was a strong negative correlation between levels of 15 proteins and the wall shear stress (R < -0.500, p ≤ 0.05), and a moderate negative correlation between 16 proteins and wall shear stress (R -0.500 to -0.300, p ≤ 0.05). Gene ontology analysis demonstrated clusters of proteins involved in cellular structure. Proteins negatively correlated with WSS included those with cytoskeletal, actin/myosin, cell-cell junction and extracellular functions. C: In aortic valve tissue, 31 proteins were identified involved in cellular structure and extracellular junctions with a negative correlation between their levels and wall shear stress. These findings suggest an association between the forces acting on the aortic valve (AV) and leaflet protein abundance, and may form a mechanical basis for the increased risk of aortic leaflet degeneration in LVAD patients.

Changes in Internal Thoracic Artery Blood Flow According to the Degree of Stenosis of the Anterior Descending Branch of the Left Coronary Artery

PURPOSE

Computational fluid dynamics has enabled the evaluation of coronary flow reserve. The purpose of this study was to clarify the hemodynamic variation and reserve potential of the left internal thoracic artery (LITA).

METHODS

Four patients were selected on the basis of various native coronary stenosis patterns and graft design. The wall shear stress and oscillatory shear index were measured, and one patient was selected. Next, we created three hypothetical lesions with 75%, 90%, and 99% stenosis in front of the graft anastomosis, and compared the changes in LITA blood flow and coronary flow distribution.

RESULTS

In the 75% to 90% stenosis model, blood flow was significantly higher in the native coronary flow proximal to the coronary artery bypass anastomosis regardless of time phase. In the 99% stenosis model, blood flow from the LITA was significantly dominant compared to native coronary flow at the proximal site of anastomosis. The range of LITA flow variability was the largest at 99% stenosis, with a difference of 70 ml/min.

CONCLUSION

The 99% stenosis model showed the highest LITA flow. The range of LITA flow variability is large, suggesting that it may vary according to the rate of native coronary stenosis.

Computational Fluid Dynamics Analysis of Blood Flow Changes during the Growth of Saccular Abdominal Aortic Aneurysm

Computational fluid dynamics analysis of the growth process of saccular abdominal aortic aneurysm was performed. A 3D model of aortic aneurysm was created based on CT images. Properties in terms of wall shear stress, mean flow velocity, mean pressure, energy loss, and pressure loss coefficient were calculated using thermal fluid analysis software "ANSYS CFX." As the aneurysm expanded, the mean flow velocity decreased and the wall shear stress, mean pressure, energy loss, and pressure loss coefficient increased. Wall shear stress increased when the aneurysm was small, suggesting that is related to the development and growth of the aneurysm. (This is secondary publication from J Jpn Coll Angiol 2021; 61: 3-10.).

Patient-specific fluid simulation of transcatheter mitral valve replacement in mitral annulus calcification

Introduction

Transcatheter mitral valve replacement is a promising alternative to open-heart surgery in elderly patients. Patients with severe mitral annulus calcification (MAC) are a particularly high-risk population, where postprocedural complications can have catastrophic effects. Amongst these, obstruction of the left ventricular outflow tract can lead to ventricular hypertrophic remodeling and subsequent heart failure, while subclinical valve thrombosis can result in early bioprosthetic valve failure.

Methods

To elucidate the mechanisms of left ventricular outflow tract obstruction and valve thrombosis following valve-in-MAC procedures, we used image processing and Computational Fluid Dynamics (CFD) software to generate patient- and device-specific models based on preprocedural CT data. Personalized computer simulations were performed to predict the left ventricular haemodynamics after implantation in three patients with severe MAC.

Results

The simulations have successfully captured the increased pressure gradient in the left ventricular outflow tract as a result of the partial obstruction due to the implanted valve. Regions of wall shear stress above the threshold value for platelet activation were also observed on the bioprosthetic frame as a result of the reduced outflow tract area, which led to increases in flow resistance and blood residence time inside the ventricle. Consistent with these findings, areas of slow recirculating flow and blood stasis formed near the valve frame, creating potential pro-thrombotic conditions.

Discussion

This study provides insight into the relationship between size and shape of the outflow tract post-implantation, pressure gradients and pro-thrombotic flow metrics such as wall shear stress and blood residence time. Results show the potential of CFD modeling to bring key functional metrics into preprocedural assessment for a comprehensive evaluation of post-procedural risks beyond anatomical factors. Following further validation and extension to the atrial chamber, this approach can provide an in-depth analysis of the likelihood of valvular thrombosis.

Modified no-touch technique for radio-cephalic arteriovenous fistula increases primary patency and decreases juxta-anastomotic stenosis

OBJECTIVE

Low primary patency rate is a major problem of radio-cephalic arteriovenous fistula (RC-AVF) creation. Radial artery deviation and reimplantation (RADAR) is associated with low juxta-anastomotic stenosis rate. However, inflow artery stenosis is prominent with RADAR. To further reduce injury to veins and arteries during operation, a modified no-touch technique (MNTT) was used to create RC-AVF.

METHODS

We retrospectively reviewed our prospectively maintained database of patients with end-stage renal disease (ESRD)s undergoing RC-AVF creation for hemodialysis using either the MNTT between January 2021 and January 2022 (MNTT group) or conventional surgical procedure ( end-to-side vein-to-artery anastomosis) between October 2016 and October 2017 (Control group). Patients who chose to undergo RC-AVF surgery underwent standardized preoperative mapping and postoperative fistula evaluations using duplex ultrasound. Additionally, 4D flow MRI data were used to visualize and quantify the hemodynamics of one RC-AVF by MNTT. Outcomes included primary patency, juxta-anastomotic stenosis, and maturation rates.

RESULTS

Forty patients underwent RC-AVFs by MNTT, compared to 60 patients in the control group. The MNTT group had a higher primary unassisted patency rate than the control group (p = 0.038). Juxta-anastomotic stenosis (all on the cephalic vein) occurred in 4 (10%) patients who underwent MNTT. RC-AVF maturation rates after 3 months were not different between both groups (maturation rate: 90% and 81.7% in the MNTT and control groups, respectively, p = 0.253). COX regression showed that both conventional AVF surgery (p = 0.031) and smaller cephalic vein diameter (p = 0.034) were associated with higher odds of RC-AVF failure. The AVF flow within the proximal vein remained helical during cardiac cycle. The distribution of wall shear stress (WSS) and oscillatory shear index (OSI) differed from that of conventional surgical AVF.

CONCLUSION

RC-AVF by MNTT increases primary patency rate and decreases juxta-anastomotic stenosis rate. The improvement in hemodynamics may be one of the important reasons for the better patency rate of in the RC-AVF by MNTT group.