An Integrated Experimental-Computational Study of Vocal Fold Vibration in Type I Thyroplasty.

Subject-specific computational modeling of vocal fold (VF) vibration was integrated with an ex vivo animal experiment of type 1 thyroplasty to study the effect of the implant on the vocal fold vibration. In the experiment, a rabbit larynx was used to simulate type 1 thyroplasty, where one side of the vocal fold was medialized with a trans-muscular suture while the other side was medialized with a silastic implant. Vocal fold vibration was then achieved by flowing air through the larynx and was filmed with a high-speed camera. The three-dimensional computational model was built upon the pre-operative scan of the laryngeal anatomy. This subject-specific model was used to simulate the vocal fold medialization and then the fluid-structure interaction (FSI) of the vocal fold. Model validation was done by comparing the vocal fold displacement with postoperative scan (for medialization), and by comparing the vibratory characteristics with the high-speed images (for vibration). These comparisons showed the computational model successfully captured the effect of the implant and thus has the potential for presurgical planning.

Fluid-Solid Interaction Analysis for Developing In-Situ Strain and Flow Sensors for Prosthetic Valve Monitoring.

Transcatheter aortic valve implantation (TAVI) was initially developed for adult patients, but there is a growing interest to expand this procedure to younger individuals with longer life expectancies. However, the gradual degradation of biological valve leaflets in transcatheter heart valves (THV) presents significant challenges for this extension. This study aimed to establish a multiphysics computational framework to analyze structural and flow measurements of TAVI and evaluate the integration of optical fiber and photoplethysmography (PPG) sensors for monitoring valve function. A two-way fluid-solid interaction (FSI) analysis was performed on an idealized aortic vessel before and after the virtual deployment of the SAPIEN 3 Ultra (S3) THV. Subsequently, an analytical analysis was conducted to estimate the PPG signal using computational flow predictions and to analyze the effect of different pressure gradients and distances between PPG sensors. Circumferential strain estimates from the embedded optical fiber in the FSI model were highest in the sinus of Valsalva; however, the optimal fiber positioning was found to be distal to the sino-tubular junction to minimize bending effects. The findings also demonstrated that positioning PPG sensors both upstream and downstream of the bioprosthesis can be used to effectively assess the pressure gradient across the valve. We concluded that computational modeling allows sensor design to quantify vessel wall strain and pressure gradients across valve leaflets, with the ultimate goal of developing low-cost monitoring systems for detecting valve deterioration.

Research on the Deformation Law of Foundation Excavation and Support Based on Fluid-Solid Coupling Theory.

To investigate the impact of underground water seepage and soil stress fields on the deformation of excavation and support structures, this study initially identified the key influencing factors on excavation deformation. Subsequently, through a finite element simulation analysis using Plaxis, this study explored the effects of critical factors, such as the excavation support form, groundwater lowering depth, permeability coefficient, excavation layer, and sequence on excavation deformation. Furthermore, a comprehensive consideration of various adverse factors was integrated to establish excavation support early warning thresholds, and optimal dewatering strategies. Finally, this study validated the simulation analysis through an on-site in situ testing with wireless sensors in the context of a physical construction site. The research results indicate that the internal support system within the excavation piles exhibited better stability compared to the external anchor support system, resulting in a 34.5% reduction in the overall deformation. Within the internal support system, the factors influencing the excavation deformation were ranked in the following order: water level (35.5%) > permeability coefficient (17.62%) > excavation layer (11.4%). High water levels, high permeability coefficients, and multi-layered soils were identified as the most unfavorable factors for excavation deformation. The maximum deformation under the coupled effect of these factors was established as the excavation support early warning threshold, and the optimal dewatering strategy involved lowering the water level at the excavation to 0.5 m below the excavation face. The on-site in situ monitoring data obtained through wireless sensors exhibited low discrepancies compared to the finite element simulation data, indicating the high precision of the finite element model for considering the fluid-structure interaction.

A mechanically consistent unified formulation for fluid-porous-structure-contact interaction.

Fluid-structure interaction with contact poses profound mathematical and numerical challenges, particularly when considering realistic contact scenarios and the influence of surface roughness. Computationally, contact introduces challenges in altering the fluid domain topology and preserving stress balance. This work introduces a new mathematical framework for a unified continuum description of fluid-porous-structure-contact interaction (FPSCI), leveraging the Navier-Stokes-Brinkman (NSB) equations to incorporate porous effects within the surface asperities in the contact region. Our approach maintains mechanical consistency during contact, circumventing issues associated with contact models and complex interface coupling conditions, allowing for the modeling of tangential creeping flows due to surface roughness. The unified continuum and variational multiscale formulation ensure robustness by enabling stable and unified integration of fluid, porous, and solid sub-problems. Computational efficiency and ease of implementation - key advantages of our approach - are demonstrated by solving two benchmark problems of a falling ball and an idealized heart valve. This research has broad implications for fields reliant on accurate fluid-structure interactions and promising advancements in modeling and numerical simulation techniques.

Mesoscopic ring element growth and deformation induced biofilm streamer evolution in microfluidic channels.

In a fluid environment, biofilms usually form and grow into streamers attached to solid surfaces. Existing research on single streamers studied their formation and failure modes. In the experiment on biofilm growth in a microfluidic channel, we found that rings composed of bacteria and an extracellular matrix are important elements on a mesoscopic scale. In the fluid environment, the failure of these ring elements causes damage to streamers. We simulated the growth and deformation of the ring structure in the micro-channel using multi-agent simulation and fluid-structure coupling of a porous elastic body. Based on this, we simulated the biofilm evolution involving multi-ring deformation, which provides a new length scale to study the biofilm streamer dynamics in fluid environments.

Numerical investigation of quantitative pulmonary pressure ratio in different degrees of stenosis.

BACKGROUND: Pulmonary artery stenosis endangers people's health. Quantitative pulmonary pressure ratio (QPPR) is very important for clinicians to quickly diagnose diseases and develop treatment plans. OBJECTIVE: Our purpose of this paper is to investigate the effects of different degrees (50% and 80%) of pulmonary artery stenosis on QPPR. METHODS: An idealized model is established based on the normal size of human pulmonary artery. The hemodynamic governing equations are solved using fluid-structure interaction. RESULTS: The results show that the QPPR decreases with the increase of stenosis degree, and it is closely related to the pressure drop at both ends of stenosis. Blood flow velocity and wall shear stress are sensitive to the stenosis degree. When the degree of stenosis is 80%, the amplitude of changes of blood flow velocity and wall shear stress at both ends of stenosis is lower. CONCLUSIONS: The results suggest that the degree of pulmonary artery stenosis has a significant impact on QPPR and hemodynamic changes. This study lays a theoretical foundation for further study of QPPR.

Design and Analysis Method of Piezoelectric Liquid Driving Device with Elastic External Displacement.

In piezoelectric drive, resonant drive is an important driving mode in which the external elastic force and electric drive signal are the key factors. In this paper, the effects of the coupling of external elastic force and liquid parameters with the structure on the vibrator resonance frequency and liquid drive are analyzed by numerical simulation. The fluid-structure coupling model for numerical analysis of the elastic force was established, the principle of microdroplet generation and the coupling method of the elastic force were studied, and the changes in the resonant frequency and mode induced by the changes in the liquid parameters in different cavities were analyzed. Through the coupled simulation and calculation of the pressure and deformation of the cavity, the laser vibration measurement test was carried out to test the effect of the vibration mode analysis. The driving model of the fluid jet driven by the elastic force on the piezoelectric drive was further established. The changing shape of the fluid jet under different elastic forces was analyzed, and the influence law of the external elastic force on the change in the droplet separation was determined. It provides reference support for further external microcontrol of droplet motion.

A parametric study regarding structural design of a bioprosthetic aortic valve by 3D fluid-structure interaction simulations.

Since the introduction of transcatheter aortic valve (AV) implantation as a viable option, surgical bioprosthetic AVs have recently started incorporating shorter struts considering future valve-in-valve procedures. However, the effect of leaflet coaptation geometry on the longevity of these valves remains unexplored. To address this gap, we performed a finite element analysis on bioprosthetic AVs with varying strut heights using a two-way fluid-structure interaction method. To establish a baseline, we used a standard height based on a rendered platform image of the CE PERIMOUNT Magna Ease valve from Edward Lifesciences in Irvine, CA. Bovine pericardium properties were assigned to the leaflets, while normal saline properties were used as the recirculating fluid in hemodynamic simulations. The physiological pressure profile of the cardiac cycle was applied between the aorta and left ventricle. We calculated blood flow velocity, effective orifice area (EOA), and mechanical stress on the leaflets. The results reveal that as the strut height increases, the stroke volume increases, leakage volume decreases, and EOA improves. Additionally, the maximum mechanical stress experienced by the leaflet decreases by 62% as the strut height increases to 1.2 times the standard height. This research highlights that a low-strut design in bioprosthetic AVs may negatively affect their durability, which can be useful in design of next-generation bioprosthetic AVs.

Thrombogenic Risk Assessment of Transcatheter Prosthetic Heart Valves Using a Fluid-Structure Interaction Approach.

Background and Objective: Prosthetic heart valve interventions such as TAVR have surged over the past decade, but the associated complication of long-term, life-threatening thrombotic events continues to undermine patient outcomes. Thus, improving thrombogenic risk analysis of TAVR devices is crucial. In vitro studies for thrombogenicity are typically difficult to perform. However, revised ISO testing standards include computational testing for thrombogenic risk assessment of cardiovascular implants. We present a fluid-structure interaction (FSI) approach for assessing thrombogenic risk of prosthetic heart valves. Methods: An FSI framework was implemented via the incompressible computational fluid dynamics multi-physics solver of the Ansys LS-DYNA software. The numerical modeling approach for flow analysis was validated by comparing the derived flow rate of the 29-mm CoreValve device from benchtop testing and orifice areas of commercial TAVR valves in the literature to in silico results. Thrombogenic risk was analyzed by computing stress accumulation (SA) on virtual platelets seeded in the flow fields via Ansys EnSight. The integrated FSI-thrombogenicity methodology was subsequently employed to examine hemodynamics and thrombogenic risk of TAVR devices with two approaches: 1) engineering optimization and 2) clinical assessment. Results: The simulated effective orifice areas of the commercial devices were in the range reported in the literature. The flow rates from the in vitro flow testing matched well with the in silico results. The approach was used to analyze the effect of various TAVR leaflet designs on hemodynamics. Platelets experienced different magnitudes of SA along their trajectories as they flowed past each design. Post-TAVR deployment hemodynamics in patient-specific bicuspid aortic valve anatomies revealed varying degrees of thrombogenic risk for these patients, despite being clinically defined as "mild" paravalvular leak. Conclusions: Our methodology can be used to improve the thromboresistance of prosthetic valves from the initial design stage to the clinic. It allows for unparalleled optimization of devices, uncovering key TAVR leaflet design parameters that can be used to mitigate thrombogenic risk, in addition to patient-specific modeling to evaluate device performance. This work demonstrates the utility of advanced in silico analysis of TAVR devices that can be utilized for thrombogenic risk assessment of other blood recirculating devices.

A hemodynamic model of artery bypass graft considering microcirculation function.

BACKGROUND: The incidence of arterial stenosis is increasing year by year. In order to better diagnose and treat arterial stenosis, numerical simulation technology has become a popular method. OBJECTIVE: A novel model is constructed to investigate the influence of microcirculation on the hemodynamics of artery bypass graft. METHODS: In this paper, a severely narrow artery bypass graft model is considered. The geometric shape includes a narrow artery tube and a bypass graft of the same diameter with a 45° suture angle. The fluid-structure interaction model is considered by finite element numerical calculation, and the flow is simulated with microcirculation as the outlet boundary condition. The changes of blood flow velocity, pressure and wall shear stress are analyzed. RESULTS: The results show that blood almost entirely flows into the graft tube and there is no recirculation area at the anastomosis. CONCLUSION: The artery bypass graft model considering microcirculation function could simulate the physiological characteristics of blood flow more reasonably, and it provide helps for clinicians to diagnose and treat arterial stenosis.

Numerical assessment of using various outlet boundary conditions on the hemodynamics of an idealized left coronary artery model.

Vascular diseases are greatly influenced by the hemodynamic parameters and the accuracy of determining these parameters depends on the use of correct boundary conditions. The present work carries out a two-way fluid-structure interaction (FSI) simulation to investigate the effects of outlet pressure boundary conditions on the hemodynamics through the left coronary artery bifurcation with moderate stenosis (50%) in the left anterior descending (LAD) branch. The Carreau viscosity model is employed to characterise the shear-thinning behaviour of blood. The results of the study reveal that the employment of zero pressure at the outlet boundaries significantly overestimates the values of hemodynamic variables like wall shear stress (WSS), and time-averaged wall shear stress (TAWSS) compared with human healthy and pulsatile pressure outlet conditions. However, the difference between these variables is marginally low for human healthy and pulsatile pressure outlets. The oscillatory shear index (OSI) remains the same across all scenarios, indicating independence from the outlet boundary condition. Furthermore, the magnitude of negative axial velocity and pressure drop across the plaque are found to be higher at the zero pressure outlet boundary condition.

Coronary hemodynamic simulation study.

In this paper, a two-way fluid-structure coupling model is developed to simulate and analyze the hemodynamic process based on dynamic coronary angiography, and examine the influence of different hemodynamic parameters on coronary arteries in typical coronary stenosis lesions. Using the measured FFR pressure data of a patient, the pressure-time function curve is fitted to ensure the accuracy of the boundary conditions. The average error of the simulation pressure results compared to the test data is 6.74%. In addition, the results related to blood flow, pressure contour and wall shear stress contour in a typical cardiac cycle are obtained by simulation analysis. These results are found to be in good agreement with the laws of the real cardiac cycle, which verifies the rationality of the simulation. In conclusion, based on the modeling and hemodynamic simulation analysis process of dynamic coronary angiography, this paper proposes a method to assist the analysis and evaluation of coronary hemodynamic and functional parameters, which has certain practical significance.

Computational fluid-structure interaction analysis of the end-to-side radio-cephalic arteriovenous fistula.

BACKGROUND AND OBJECTIVE: In the current work, we present a descriptive fluid-structure interaction computational study of the end-to-side radio-cephalic arteriovenous fistula. This allows us to account for the different thicknesses and elastic properties of the radial artery and cephalic vein. METHODS: The core of the work consists in simulating different arteriovenous fistula configurations obtained by virtually varying the anastomosis angle, i.e. the angle between the end of the cephalic vein and the side of the radial artery. Since the aim of the work is to understand the blood dynamics in the very first days after the surgical intervention, the radial artery is considered stiffer and thicker than the cephalic vein. RESULTS: Our results demonstrate that both the diameter of the cephalic vein and the anastomosis angle play a crucial role to obtain a blood dynamics without re-circulation regions that could prevent fistula failure. CONCLUSIONS: When an anastomosis angle close to the perpendicular direction with respect to the radial artery is combined with a large diameter of the cephalic vein, the recirculation regions and the low Wall Shear Stress (WSS) zones are reduced. Conversely, from a structural point of view, a low anastomosis angle with a large diameter of the cephalic vein reduces the mechanical stress acting on the vessel walls.

Impact of palatopharyngeal sizes changing on pharyngeal airflow fluctuation and airway vibration in a pediatric airway.

Snoring is common in children and is associated with many adverse consequences. One must study the relationships between pharyngeal morphology and snoring physics to understand snoring progression. Although some model studies have provided fluid-structure interaction dynamic descriptions for the correlation between airway size and snoring physics, the descriptions still need to be further investigated in patient-specific airway models. Fluid-structure interaction studies using patient-specific airway structures complement the above model studies. Based on reported cephalometric measurement methods, this study quantified and preset the size of the palatopharynx airway in a patient-specific airway and investigated how the palatopharynx size affects the pharyngeal airflow fluctuation, soft palate vibration, and glossopharynx vibration with the help of a verified FSI method. The results showed that the stenosis anterior airway of the soft palate increased airway resistance and airway resistance fluctuations, which can lead to increased sleep effort and frequent snoring. Widening of the anterior airway can reduce airflow resistance and avoid obstructing the anterior airway by the soft palate vibration. The pharyngeal airflow resistance, mouth inflow proportion, and soft palate apex displacement have components at the same frequencies in all airway models, and the glossopharynx vibration and instantaneous inflow rate have components at the same frequencies, too. The mechanism of this same frequency fluctuation phenomenon can be explained by the fluid-structure interaction dynamics of an ideal coupled model consisting of a flexible plate model and a collapsible tube model. The results of this study demonstrate the potential of FSI in studying snoring physics and clarify to some degree the mechanism of airway morphology affecting airway vibration physics.

Investigation of the ultrasound-induced collapse of air bubbles near soft materials.

A numerical investigation into the ultrasound-induced collapse of air bubbles near soft materials, utilizing a novel multi-material diffuse interface method (DIM) model with block-structured adaptive mesh refinement is presented. The present work expands from a previous five-equation DIM by incorporating Eulerian hyperelasticity. The model is applicable to any arbitrary number of interacting fluid and solid material. A single conservation law for the elastic stretch tensor enables tracking the deformations for all the solid materials. A series of benchmark cases are conducted, and the solution is found to be in excellent agreement against theoretical data. Subsequently, the ultrasound-induced bubble-tissue flow interactions are examined. The bubble radius was found to play a crucial role in dictating the stresses experienced by the tissue, underscoring its significance in medical applications. The results reveal that soft tissues primarily experience tensile forces during these interactions, suggesting potential tensile-driven injuries that may occur in relevant treatments. Moreover, regions of maximal tensile forces align with tissue elongation areas. It is documented that while early bubble dynamics remain relatively unaffected by changes in shear modulus, at later stages of the penetration processes and the deformation shapes, exhibit notable variations. Lastly, it is demonstrated that decreasing standoff distances enhances the interaction between bubbles and tissue, thereby increasing the stress levels in the tissue, although the behavior of the bubble dynamics remains largely unchanged.

Analysis of Intracranial Aneurysm Haemodynamics Altered by Wall Movement.

Computational fluid dynamics is intensively used to deepen our understanding of aneurysm growth and rupture in an attempt to support physicians during therapy planning. Numerous studies assumed fully rigid vessel walls in their simulations, whose sole haemodynamics may fail to provide a satisfactory criterion for rupture risk assessment. Moreover, direct in vivo observations of intracranial aneurysm pulsation were recently reported, encouraging the development of fluid-structure interaction for their modelling and for new assessments. In this work, we describe a new fluid-structure interaction functional setting for the careful evaluation of different aneurysm shapes. The configurations consist of three real aneurysm domes positioned on a toroidal channel. All geometric features, employed meshes, flow quantities, comparisons with the rigid wall model and corresponding plots are provided for the sake of reproducibility. The results emphasise the alteration of flow patterns and haemodynamic descriptors when wall deformations were taken into account compared with a standard rigid wall approach, thereby underlining the impact of fluid-structure interaction modelling.

Patient-specific analysis of bicuspid aortic valve hemodynamics using a fully coupled fluid-structure interaction model.

Bicuspid aortic valve (BAV), the most common congenital heart disease, is prone to develop significant valvular dysfunction and aortic wall abnormalities such as ascending aortic aneurysm. Growing evidence has suggested that abnormal BAV hemodynamics could contribute to disease progression. In order to investigate BAV hemodynamics, we performed 3D patient-specific fluid-structure interaction (FSI) simulations with fully coupled blood flow dynamics and valve motion throughout the cardiac cycle. Results showed that the hemodynamics during systole can be characterized by a systolic jet and two counter-rotating recirculation vortices. At peak systole, the jet was usually eccentric, with asymmetric recirculation vortices and helical flow motion in the ascending aorta. The flow structure at peak systole was quantified using the vorticity, flow rate reversal ratio and local normalized helicity (LNH) at four locations from the aortic root to the ascending aorta. The systolic jet was evaluated with the peak velocity, normalized flow displacement, and jet angle. It was found that peak velocity and normalized flow displacement (rather than jet angle) gave a strong correlation with the vorticity and LNH in the ascending aorta, which suggests that these two metrics could be used for clinical noninvasive evaluation of abnormal blood flow patterns in BAV patients.

Study on the optimal elastic modulus of flexible blades for right heart assist device supporting patients with single-ventricle physiologies.

Background: Patients with single-ventricle physiologies continue to experience insufficient circulatory power after undergoing palliative surgeries. This paper proposed a right heart assist device equipped with flexible blades to provide circulatory assistance for these patients. The optimal elastic modulus of the flexible blades was investigated through numerical simulation. Methods: A one-way fluid-structure interaction (FSI) simulation was employed to study the deformation of flexible blades during rotation and its impact on device performance. The process began with a computational fluid dynamics (CFD) simulation to calculate the blood pressure rise and the pressure on the blades' surface. Subsequently, these pressure data were exported for finite element analysis (FEA) to compute the deformation of the blades. The fluid domain was then recreated based on the deformed blades' shape. Iterative CFD and FEA simulations were performed until both the blood pressure rise and the blades' shape stabilized. The blood pressure rise, hemolysis risk, and thrombosis risk corresponding to blades with different elastic moduli were exhaustively evaluated to determine the optimal elastic modulus. Results: Except for the case at 8,000 rpm with a blade elastic modulus of 40 MPa, the pressure rise associated with flexible blades within the studied range (rotational speeds of 4,000 rpm and 8,000 rpm, elastic modulus between 10 MPa and 200 MPa) was lower than that of rigid blades. It was observed that the pressure rise corresponding to flexible blades increased as the elastic modulus increased. Additionally, no significant difference was found in the hemolysis risk and thrombus risk between flexible blades of various elastic moduli and rigid blades. Conclusion: Except for one specific case, deformation of the flexible blades within the studied range led to a decrease in the impeller's functionality. Notably, rotational speed had a more significant impact on hemolysis risk and thrombus risk compared to blade deformation. After a comprehensive analysis of blade compressibility, blood pressure rise, hemolysis risk, and thrombus risk, the optimal elastic modulus for the flexible blades was determined to be between 40 MPa and 50 MPa.

Research of biomimetic corrugation on the blade flutter suppression in large-scale wind turbine systems.

Aiming at the blade flutter of large horizontal-axis wind turbines, a method by utilizing biomimetic corrugation to suppress blade flutter is first proposed. By extracting the dragonfly wing corrugation, the biomimetic corrugation airfoil is constructed, finding that mapping corrugation to the airfoil pressure side has better aerodynamic performance. The influence of corrugation type, amplitudeλ, and intensity on airfoil flutter is analyzed using orthogonal experiment, which determines that theλhas the greatest influence on airfoil flutter. Based on the fluctuation range of the moment coefficient ΔCm, the optimal airfoil flutter suppression effect is obtained when the type is III,λ= 0.6, and intensity is denser (n= 13). The effective corrugation layout area in the chord direction is determined to be the leading edge, and the ΔCmof corrugation airfoil is reduced by 7.405%, compared to the original airfoil. The application of this corrugation to NREL 15 MW wind turbine 3D blades is studied, and the influence of corrugation layout length in the blade span direction on the suppressive effect is analyzed by fluid-structure interaction. It is found that when the layout length is 0.85 R, the safety marginSfreaches a maximum value of 0.3431 Hz, which is increased 2.940%. The results show that the biomimetic corrugated structure proposed in this paper can not only improve the aerodynamic performance by changing the local flow field on the surface of the blade, but also increase the structural stiffness of the blade itself, and achieve the effect of flutter suppression.

In silico analysis of embolism in cerebral arteries using fluid-structure interaction method.

Ischemic stroke, particularly embolic stroke, stands as a significant global contributor to mortality and long-term disabilities. This paper presents a comprehensive simulation of emboli motion through the middle cerebral artery (MCA), a prevalent site for embolic stroke. Our patient-specific computational model integrates major branches of the middle cerebral artery reconstructed from magnetic resonance angiography images, pulsatile flow dynamics, and emboli of varying geometries, sizes, and material properties. The fluid-structure interactions method is employed to simulate deformable emboli motion through the middle cerebral artery, allowing observation of hemodynamic changes in artery branches upon embolus entry. We investigated the impact of embolus presence on shear stress magnitude on artery walls, analyzed the effects of embolus material properties and geometries on embolus trajectory and motion dynamics within the middle cerebral artery. Additionally, we evaluated the non-Newtonian behavior of blood, comparing it with Newtonian blood behavior. Our findings highlight that embolus geometry significantly influences trajectory, motion patterns, and hemodynamics within middle cerebral artery branches. Emboli with visco-hyperelastic material properties experienced higher stresses upon collision with artery walls compared to those with hyperelastic properties. Furthermore, considering blood as a non-Newtonian fluid had notable effects on emboli stresses and trajectories within the artery, particularly during collisions. Notably, the maximum von Mises stress experienced in our study was 21.83 kPa, suggesting a very low probability of emboli breaking during movement, impact, and after coming to a stop. However, in certain situations, the magnitude of shear stress on them exceeded 1 kPa, increasing the likelihood of cracking and disintegration. These results serve as an initial step in anticipating critical clinical conditions arising from arterial embolism in the middle cerebral artery. They provide insights into the biomechanical parameters influencing embolism, contributing to improved clinical decision-making for stroke management.