Numerical Simulation of Wall Shear Stress and Particle-Based Hemodynamic Parameters in Pre-Cuffed and Streamlined End-to-Side Anastomoses

Computer Modeling of Carotid Endarterectomy With the Different Shape Patches and Prediction of the Atherosclerotic Plaque Formation Zones

The article describes a method for constructing geometric models of the carotid bifurcation and computer simulation of endarterectomy surgery with the patches of various configurations. The purpose of this work is to identify the areas of the greatest risk of restenosis in the constructed models and to conduct a comparative analysis of risk factors when using the patches of different widths and shapes. The method is demonstrated on a reconstructed model of a healthy vessel. Its building is based on a preoperative computed tomography study of a particular patient's affected vessel. The flow in the vessel is simulated by computational fluid dynamics using data from the patient's ultrasound Doppler velocimetry. Risk factors are assessed through the hemodynamic indices on the vessel wall associated with Wall Shear Stress. The distribution of risk zones in the healthy vessel, presumably leading to its observed lesion (plaque), is analyzed. Comparative evaluation of 10 various patches implantation results is carried out and the optimal variant is determined. The proposed method can be used to predict the hemodynamic results of surgery using patches of various sizes and shapes.

Patient-Specific Hemodynamics of New Coronary Artery Bypass Configurations

PURPOSE

This study aims to quantify the patient-specific hemodynamics of complex conduit routing configurations of coronary artery bypass grafting (CABG) operation which are specifically suitable for off-pump surgeries. Coronary perfusion efficacy and local hemodynamics of multiple left internal mammary artery (LIMA) with sequential and end-to-side anastomosis are investigated. Using a full anatomical model comprised of aortic arch and coronary artery branches the optimum perfusion configuration in multi-vessel coronary artery stenosis is desired.

METHODOLOGY

Two clinically relevant CABG configurations are created using a virtual surgical planning tool where for each configuration set, the stenosis level, anastomosis distance and angle were varied. A non-Newtonian computational fluid dynamics solver in OpenFOAM incorporated with resistance boundary conditions representing the coronary perfusion physiology was developed. The numerical accuracy is verified and results agreed well with a validated commercial cardiovascular flow solver and experiments. For segmental performance analysis, new coronary perfusion indices to quantify deviation from the healthy scenario were introduced.

RESULTS

The first simulation configuration set;-a CABG targeting two stenos sites on the left anterior descending artery (LAD), the LIMA graft was capable of 31 mL/min blood supply for all the parametric cases and uphold the healthy LAD perfusion in agreement with the clinical experience. In the second end-to-side anastomosed graft configuration set;-the radial artery graft anastomosed to LIMA, a maximum of 64 mL/min flow rate in LIMA was observed. However, except LAD, the obtuse marginal (OM) and second marginal artery (m2) suffered poor perfusion. In the first set, average wall shear stress (WSS) were in the range of 4 to 35 dyns/cm² for in LAD. Nevertheless, for second configuration sets the WSS values were higher as the LIMA could not supply enough blood to OM and m2.

CONCLUSION

The virtual surgical configurations have the potential to improve the quality of operation by providing quantitative surgical insight. The degree of stenosis is a critical factor in terms of coronary perfusion and WSS. The sequential anastomosis can be done safely if the anastomosis angle is less than 90 degrees regardless of degree of stenosis. The smaller proposed perfusion index value, O(0.04 - 0) × 10², enable us to quantify the post-op hemodynamic performance by comparing with the ideal healthy physiological flow.

NUMERICAL INVESTIGATION OF THE EFFECTS OF BED SHAPE ON THE END-TO-SIDE CABG HEMODYNAMICS

Disrupted flow initiates and aggravates intimal thickening in the end-to-side (ETS) coronary artery bypass grafting (CABG), which may lead to failure. To enhance the post-intervention hemodynamics, the geometry is either optimized or totally reconfigured. Majority of configurations proposed by researchers have not suited CABG surgery, for they entailed rigorous manipulation on conventional grafts in situ, which was neither swift nor straightforward. The aim of the present study is, thus, to introduce a slight, yet effective, modification to a conventional ETS CABG configuration, and numerically investigate its effects on updated hemodynamic and structural environment, anticipating the longevity of proposed configuration and CABG success. This fairly simple modification may easily be made positioning a pre-designed anastomotic device between the bed of host artery in the conventional ETS CABG and its surrounding tissues. Conducting comprehensive numerical simulations, performance of the proposed configuration was assessed using idealized and patient-specific geometries of the conventional ETS CABG. Blood flow was simulated in a conventional and an updated CABG configuration considering 2-way fluid–structure interaction. Results revealed that, although the proposed configuration may induce higher structural stresses in vessels walls, it may improve important hemodynamic metrics such as wall shear stress gradient, oscillatory shear index, and relative residence time on host artery bed reducing disruption of flow. This study may also set the stage for design engineers and regulatory officials to evolve ETS CABG toward more hemodynamics-friendly approaches. Further in vitro, preclinical, and clinical experiments are, yet, entailed to accomplish ideal designs of procedural guidelines/grafts.

Analysis of Computational Fluid Dynamics and Particle Image Velocimetry Models of Distal-End Side-to-Side and End-to-Side Anastomoses for Coronary Artery Bypass Grafting in a Pulsatile Flow

BACKGROUND

Intimal hyperplasia (IH) is a major cause of graft failure. Hemodynamic factors such as stagnation and disturbed blood flow are involved in IH formation. The aim of this study is to perform a comparative analysis of distal-end side-to-side (deSTS) and end-to-side (ETS) anastomoses using computational fluid dynamics (CFD) after validating the results via particle image velocimetry (PIV).

METHODS AND RESULTS

We investigated the characteristics of our target flow fields using CFD under steady and pulsatile flows. CFD via PIV under steady flow in a 10-times-actual-size model was validated. The CFD analysis revealed a recirculation zone in the heel region in the deSTS and ETS anastomoses and at the distal end of the graft, and just distal to the toe of the host artery in the deSTS anastomoses. The recirculation zone sizes changed with the phase shift. We found regions of low wall shear stress and high oscillating shear index in the same areas. The PIV and CFD results were similar.

CONCLUSIONS

It was demonstrated that the hemodynamic characteristics of CFD and PIV is the difference between the deSTS and ETS anastomoses; that is, the deSTS flow peripheral to the distal end of the graft, at the distal end and just distal to the toe of the host artery is involved in the IH formation.

Associations between Carotid Artery Plaque Score, Carotid Hemodynamics and Coronary Heart Disease

Background: The carotid artery plaque score (PS) is an independent predictor of Coronary Heart Disease (CHD). This study aims to evaluate the combination of PS and carotid hemodynamics to predict CHD. Methods: A total of 476 patients who underwent carotid ultrasonography and coronary angiography were divided into two groups depending on the presence of CHD. PS, carotid intima-media thickness, and carotid blood flow were measured. Receiver operating characteristic curve analysis was performed to establish the best prediction model for CHD presence. Results: Age, sex, carotid intima-media thickness of internal carotid artery and carotid bifurcation, PS, peak systolic velocity (PSA) of right internal carotid artery (RICA), and most resistance index data were significantly related with the presence of CHD. The area under the curve for a collective model, which included factors of the PS, carotid hemodynamics and age, was significantly higher than the other model. Age, PS, and PSA of RICA were significant contributors for predicting CHD presence. Conclusions: The model of PS and PSA of RICA has greater predictive value for CHD than PS alone. Adding age to PS and PSA of RICA further improves predictive value over PS alone.

Effect of swirling blood flow on vortex formation at post-stenosis

Various clinical observations reported that swirling blood flow is a normal physiological flow pattern in various vasculatures. The swirling flow has beneficial effects on blood circulation through the blood vessels. It enhances oxygen transfer and reduces low-density lipoprotein concentration in the blood vessel by enhancing cross-plane mixing of the blood. However, the fluid-dynamic roles of the swirling flow are not yet fully understood. In this study, inhibition of material deposition at the post-stenosis region by the swirling flow was observed. To reveal the underlying fluid-dynamic characteristics, pathline flow visualization and time-resolved particle image velocimetry measurements were conducted. Results showed that the swirling inlet flow increased the development of vortices at near wall region of the post-stenosis, which can suppress further development of stenosis by enhancing transport and mixing of the blood flow. The fluid-dynamic characteristics obtained in this study would be useful for improving hemodynamic characteristics of vascular grafts and stents in which the stenosis frequently occurred. Moreover, the time-resolved particle image velocimetry measurement technique and vortex identification method employed in this study would be useful for investigating the fluid-dynamic effects of the swirling flow on various vascular environments.

Lagrangian postprocessing of computational hemodynamics

Recent advances in imaging, modeling, and computing have rapidly expanded our capabilities to model hemodynamics in the large vessels (heart, arteries, and veins). This data encodes a wealth of information that is often under-utilized. Modeling (and measuring) blood flow in the large vessels typically amounts to solving for the time-varying velocity field in a region of interest. Flow in the heart and larger arteries is often complex, and velocity field data provides a starting point for investigating the hemodynamics. This data can be used to perform Lagrangian particle tracking, and other Lagrangian-based postprocessing. As described herein, Lagrangian methods are necessary to understand inherently transient hemodynamic conditions from the fluid mechanics perspective, and to properly understand the biomechanical factors that lead to acute and gradual changes of vascular function and health. The goal of the present paper is to review Lagrangian methods that have been used in post-processing velocity data of cardiovascular flows.

Coronary artery bypass grafting hemodynamics and anastomosis design: a biomedical engineering review

In this paper, coronary arterial bypass grafting hemodynamics and anastomosis designs are reviewed. The paper specifically addresses the biomechanical factors for enhancement of the patency of coronary artery bypass grafts (CABGs). Stenosis of distal anastomosis, caused by thrombosis and intimal hyperplasia (IH), is the major cause of failure of CABGs. Strong correlations have been established between the hemodynamics and vessel wall biomechanical factors and the initiation and development of IH and thrombus formation. Accordingly, several investigations have been conducted and numerous anastomotic geometries and devices have been designed to better regulate the blood flow fields and distribution of hemodynamic parameters and biomechanical factors at the distal anastomosis, in order to enhance the patency of CABGs. Enhancement of longevity and patency rate of CABGs can eliminate the need for re-operation and can significantly lower morbidity, and thereby reduces medical costs for patients suffering from coronary stenosis. This invited review focuses on various endeavors made thus far to design a patency-enhancing optimized anastomotic configuration for the distal junction of CABGs.

Balloon-assisted maturation (BAM) of the arteriovenous fistula: the good, the bad, and the ugly

Arteriovenous fistula (AVF) maturation failure is currently a huge clinical problem. One approach to enhance the AVF maturation rate is an aggressive sequence of balloon angioplasty procedures, often known as balloon-assisted maturation. The goal of the current paper is to explore the pros and cons of this procedure and to try and better identify its impact on AVF maturation.

A novel method for the generation of multi-block computational structured grids from medical imaging of arterial bifurcations

In this study a description of a new approach, for the generation of multi-block structured computational grids on patient-specific bifurcation geometries is presented. The structured grid generation technique is applied to data obtained by medical imaging examination, resulting in a surface conforming, high quality, multi-block structured grid of the branching geometry. As a case study application a patient specific abdominal aorta bifurcation is selected. For the evaluation of the grid produced by the novel method, a grid convergence study and a comparison between the grid produced by the method and unstructured grids produced by commercial meshing software are carried out.

Full-hexahedral structured meshing for image-based computational vascular modeling

Image-based computational modeling offers a virtual access to spatially and temporally high resolution flow and structural mechanical data in vivo. Due to inter-subject morphological variability, mesh generation represents a critical step in modeling the patient-specific geometry and is usually performed using unstructured tetrahedral meshing algorithms. Although hexahedral structured meshes are known to provide higher accuracy and reduce the computational costs both for Finite Element Analysis and Computational Fluid Dynamics, their application in computational cardiovascular studies is challenging due to the complex 3D and branching topology of vascular territories. In this study, we propose a robust procedure for structured mesh generation, tailoring the mesh structure to the subject-specific vessel topology. The proposed methodology is based on centerline-based synthetic descriptors (i.e. centerlines, radii and centerlines' normals) which are used to solve the meshing problem following a bottom-up approach. First, topologically equivalent block-structures are placed inside and outside the lumen domain. Then, a projection operation is performed, returning a parametric volume mesh which fits the original triangulated model with sub-micrometric accuracy. Additionally, a three-layered arterial wall (resembling the intima, media and adventitia) is artificially generated, with the possibility of setting variable thickness (e.g. proximal-to-distal tapering) and material anisotropy (e.g. position-dependent collagen-fibers' orientation). This new meshing procedure, implemented using open-source software packages only, is demonstrated on two challenging human cases, being an aortic arch and an abdominal aortic aneurysm. High-quality meshes are generated in both cases, according to shape-quality metrics. By increasing the computation accuracy, the developed meshing tool has the potential to further add "confidence" to the use of computational methods in vascular applications.

Real-time measurement of solute transport within the lacunar-canalicular system of mechanically loaded bone: direct evidence for load-induced fluid flow

Since proposed by Piekarski and Munro in 1977, load-induced fluid flow through the bone lacunar-canalicular system (LCS) has been accepted as critical for bone metabolism, mechanotransduction, and adaptation. However, direct unequivocal observation and quantification of load-induced fluid and solute convection through the LCS have been lacking due to technical difficulties. Using a novel experimental approach based on fluorescence recovery after photobleaching (FRAP) and synchronized mechanical loading and imaging, we successfully quantified the diffusive and convective transport of a small fluorescent tracer (sodium fluorescein, 376 Da) in the bone LCS of adult male C57BL/6J mice. We demonstrated that cyclic end-compression of the mouse tibia with a moderate loading magnitude (-3 N peak load or 400 µε surface strain at 0.5 Hz) and a 4-second rest/imaging window inserted between adjacent load cycles significantly enhanced (+31%) the transport of sodium fluorescein through the LCS compared with diffusion alone. Using an anatomically based three-compartment transport model, the peak canalicular fluid velocity in the loaded bone was predicted (60 µm/s), and the resulting peak shear stress at the osteocyte process membrane was estimated (∼5 Pa). This study convincingly demonstrated the presence of load-induced convection in mechanically loaded bone. The combined experimental and mathematical approach presented herein represents an important advance in quantifying the microfluidic environment experienced by osteocytes in situ and provides a foundation for further studying the mechanisms by which mechanical stimulation modulates osteocytic cellular responses, which will inform basic bone biology, clinical understanding of osteoporosis and bone loss, and the rational engineering of their treatments.

An extended convection diffusion model for red blood cell-enhanced transport of thrombocytes and leukocytes

Transport phenomena of platelets and white blood cells (WBCs) are fundamental to the processes of vascular disease and thrombosis. Unfortunately, the dilute volume occupied by these cells is not amenable to fluid-continuum modeling, and yet the cell count is large enough that modeling each individual cell is impractical for most applications. The most feasible option is to treat them as dilute species governed by convection and diffusion; however, this is further complicated by the role of the red blood cell (RBC) phase on the transport of these cells. We therefore propose an extended convection-diffusion (ECD) model based on the diffusive balance of a fictitious field potential, Psi, that accounts for the gradients of both the dilute phase and the local hematocrit. The ECD model was applied to the flow of blood in a tube and between parallel plates in which a profile for the RBC concentration field was imposed and the resulting platelet concentration field predicted. Compared to prevailing enhanced-diffusion models that dispersed the platelet concentration field, the ECD model was able to simulate a near-wall platelet excess, as observed experimentally. The extension of the ECD model depends only on the ability to prescribe the hematocrit distribution, and therefore may be applied to a wide variety of geometries to investigate platelet-mediated vascular disease and device-related thrombosis.

Mechanical analysis of end-to-end silk-sutured anastomosis for robot-assisted surgery

BACKGROUND

Robot-assisted anastomosis holds great promise for the future. To secure surgery quality, some key process factors, such as the force arrangement of sutures, should be provided because of the lack of haptic feedback in robotics systems

METHODS

A model of anastomosis is presented to establish the mechanical relationship between vessel and sutures. Stress distribution of the vessel loaded by the suture was then achieved through finite-element simulations, based on the material property test results. Further, experiments were performed to validate the reliability of the FEM simulation of the anastomosis process.

RESULTS

To avoid blood osmosis, the allowable lower limit of the suture tension was 0.05 N. To keep the tissue free from injury, the allowable upper limit of tension on the suture was 0.4 N.

CONCLUSIONS

The study provided meaningful results for directing the robot-assisted anastomosis procedure and design of the surgical tools.

Flow analyses in the lower airways: patient-specific model and boundary conditions

Computational fluid dynamics (CFD) is increasingly applied in the respiratory domain. The ability to simulate the flow through a bifurcating tubular system has increased the insight into the internal flow dynamics and the particular characteristics of respiratory flows such as secondary motions and inertial effects. The next step in the evolution is to apply the technique to patient-specific cases, in order to provide more information about pathological airways. This study presents a patient-specific approach where both the geometry and the boundary conditions (BC) are based on individual imaging methods using computed tomography (CT). The internal flow distribution of a 73-year-old female suffering from chronic obstructive pulmonary disease (COPD) is assessed. The validation is performed through the comparison of lung ventilation with gamma scintigraphy. The results show that in order to obtain agreement within the accuracy limits of the gamma scintigraphy scan, both the patient-specific geometry and the BC (driving pressure) play a crucial role. A minimal invasive test (CT scan) supplied enough information to perform an accurate CFD analysis. In the end it was possible to capture the pathological features of the respiratory system using the imaging and computational fluid dynamics techniques. This brings the introduction of this new technique in the clinical practice one step closer.