High Shear Stress Regions Research Articles

Modeling Fibrin Accumulation on Flow-Diverting Devices for Intracranial Aneurysms.

The mechanisms leading to aneurysm occlusion after treatment with flow-diverting devices are not fully understood. Flow modification induces thrombus formation within the aneurysm cavity, but fibrin can simultaneously accumulate and cover the device scaffold, leading to further flow modification. However, the interplay and relative importance of these processes are not clearly understood. A computational model of fibrin accumulation and flow modification after flow diversion treatment of cerebral aneurysms has been developed under the guidance of invitro experiments and observations. The model is based on the loose coupling of flow and transport-reaction equations that are solved separately by independent codes. Interaction or reactive terms account for thrombin production from prothrombin stimulated by thrombogenic metallic wires and inhibition by antithrombin as well as fibrin production from fibrinogen stimulated by thrombin and flow shear stress, and fibrin adhesion to device wires and already attached fibrin. The computational model was demonstrated and tested on idealized vessel and aneurysm geometries. The model was able to reproduce the salient features of fibrin accumulation after the deployment of flow-diverting devices in idealized invitro models of cerebral aneurysms. Namely, fibrin production in regions of high shear stress, initial accumulation at the inflow zone, and progressive occlusion of the device and corresponding flow attenuation. The computational model linking flow dynamics to fibrin production, transport, and adhesion can be used to investigate and better understand the effects that lead to fibrin accumulation and the resulting aneurysm inflow reduction and intra-aneurysmal flow modulation.

Alterations in nasal airflow and air conditioning after endoscopic nasopharyngectomy for recurrent nasopharyngeal carcinoma: a pilot computational fluid dynamics study

Endoscopic nasopharyngectomy represents a significant intervention for recurrent nasopharyngeal carcinoma (NPC). Various surgical techniques, including transnasal and transoral approaches, are employed. However, the impact of these procedures on nasal airflow dynamics is not well understood. This computational fluid dynamics (CFD) study aimed to investigate alterations in nasal airflow and air conditioning following endoscopic nasopharyngectomy. A 55-year-old male patient with recurrent NPC was selected, whose CT data were utilized for image reconstruction. A preoperative model and two postoperative models, including the transnasal and transoral approach models, were established. The airflow patterns and various CFD parameters were analyzed. In the postoperative models, the high-speed airflow went along the soft palate and into the nasopharyngeal outlet, and there was the low-speed turbulence in the expanded nasopharyngeal cavity. Compared to the preoperative model, the postoperative models exhibited reductions in surface-to-volume ratio, nasal resistance, airflow velocity and proportion of high wall shear stress regions in nasopharynx. The changing trends of nasopharyngeal air temperature and humidity in the preoperative and transoral models were consistent. The heating and humidification efficiency decreased in the transnasal model compared to the transoral model. The endoscopic nasopharyngectomy for recurrent NPC affects the nasal airflow and warming and humidification function. The transoral approach has less influence on aerodynamics of the upper airway compared to the transnasal approach. From a CFD perspective, the endoscopic nasopharyngectomy does not increase the risk of postoperative complications, including the empty nose syndrome and the carotid blowout syndrome.

Computational model-based hemodynamic comparisons of traditional and modified idealized models of autologous radiocephalic fistula.

Autologous arteriovenous fistula (AVF) is a commonly used vascular access (VA) for hemodialysis, and hemodynamic changes are one of the main factors for its failure. To explore the effect of geometry on the hemodynamics in the AVF, a modified model is built with a gradual and smooth turn at the anastomosis and is compared with the traditional model, which has an abrupt sharp turn at the anastomisis. Transient computational fluid dynamics (CFD) simulations were performed for the comparison and analysis of the hemodynamic fields of the two models at different stages of the pulse cycle. The results showed that the low shear stress region and high oscillatory shear stress region in the modified AVF model coincided with regions of intimal hyperplasia that have been identified by previous studies. A comparison with the blood flow velocities measured in vivo was performed, and the error between the simulation results and the medical data was reduced by 22% in the modified model, which verifies the rationality and utility of the modified model.

Ultrasonic irrigation flows in root canals: effects of ultrasound power and file insertion depth

Ultrasonic irrigation during root canal treatment can enhance biofilm disruption. The challenge is to improve the fluid flow so that the irrigant reaches areas inaccessible to hand instrumentation. The aim of this study is to experimentally investigate how the flow field and hydrodynamic forces induced by ultrasonic irrigation are influenced by the ultrasound power and file insertion depth. A root canal phantom was 3D printed and used as a mold for the fabrication of a PDMS channel. An ultrasonic instrument with a #15K-file provided the irrigation. The flow field was studied by means of Particle Image Velocimetry (PIV). The time averaged velocity and shear stress distributions were found to vary significantly with ultrasound power. Their maximum values increase sharply for low powers and up to a critical power level. At and above this setting, the flow pattern changes, from the high velocity and shear stress region confined in the vicinity of the tip, to one covering the whole root canal domain. Exceeding this threshold also induces a moderate increase in the maximum velocities and shear stresses. The insertion depth was found to have a smaller effect on the measured velocity and shear stresses. Due to the oscillating nature of the flow, instantaneous maximum velocities and shear stresses can reach much higher values than the mean, especially for high powers. Ultrasonic irrigation will benefit from using a higher power setting as this does produce greater shear stresses near the walls of the root canal leading to the potential for increased biofilm removal.

Development and validation of a mathematical model for evaluating shear-induced damage of von Willebrand factor

PurposeTo develop a mathematical model for predicting shear-induced von Willebrand factor (vWF) function modification which can be used to guide ventricular assist devices (VADs) design, and evaluate the damage of high molecular weight multimers (HMWM)-vWF in VAD patients for reducing clinical complications. MethodsMathematical models were constructed based on three morphological variations (globular vWF, unfolded vWF and degraded vWF) of vWF under shear stress conditions, in which parameters were obtained from previous studies or fitted by experimental data. Different clinical support modes (pediatric vs. adult mode), different VAD operating states (pulsation vs. constant mode) and different clinical VADs (HeartMate II, HeartWare and CentriMag) were utilized to analyze shear-induced damage of HMWM-vWF based on our vWF model. The accuracy and feasibility of the models were evaluated using various experimental and clinical cases, and the biomechanical mechanisms of HMWM-vWF degradation induced by VADs were further explained. ResultsThe mathematical model developed in this study predicted VAD-induced HMWM-vWF degradation with high accuracy (correlation with experimental data r2 > 0.99). The numerical results showed that VAD in the pediatric mode resulted in more HMWM-vWF degradation per unit time and per unit flow rate than in the adult mode. However, the total degradation of HMWM-vWF is less in the pediatric mode than in the adult mode because the pediatric mode has fewer times of blood circulation than the adult mode in the same amount of time. The ratio of HMWM-vWF degradation was lower in the pulsation mode than in the constant mode. This is due to the increased flushing of VADs in the pulsation mode, which avoids prolonged stagnation of blood in high shear regions. This study also found that the design feature, rotor size and volume of the VADs, and the superimposed regions of high shear stress and long residence time inside VADs affect the degradation of HMWM-vWF. The axial flow VADs (HeartMate II) showed higher degradation of HMWM-vWF compared to centrifugal VADs (HeartWare and CentriMag). Compared to fully magnetically suspended VADs (CentriMag), hydrodynamic suspended VADs (HeartWare) produced extremely high degradation of HWMW-vWF in its narrow hydrodynamic clearance. Finally, the study used a mathematical model of HMWM-vWF degradation to interpret the clinical statistics from a biomechanical perspective and found that minimizing the rotating speed of VADs within reasonable limits helps to reduce HWMW-vWF degradation. All predicted conclusions are supported by the experimental and clinical data. ConclusionThis study provides a validated mathematical model to assess the shear-induced degradation of HMWM-vWF, which can help to evaluate the damage of HMWM-vWF in patients implanted with VADs for reducing clinical complications, and to guide the optimization of VADs for improving hemocompatibility.

Hostile Hemodynamics in Distal Stent Graft-Induced New Entry Prior to Aortic Rupture: A Comparison of Transient versus Steady-State CFD Simulations.

Computational fluid dynamics (CFD) simulations model blood flow in aortic pathologies. The aim of our study was to understand the local hemodynamic environment at the site of rupture in distal stent graft-induced new entry (dSINE) after frozen elephant trunk with a clinically time efficient steady-flow simulation versus transient simulations. Steady-state simulations were performed for dSINE, prior and after its development and prior to aortic rupture. To account for potential turbulences due geometric changes at the dSINE location, Reynolds-averaged Navier-Stokes equations with the realizable k-ε model for turbulences were applied. Transient simulations were performed for comparison. Hemodynamic parameters were assessed at various locations of the aorta. Post-dSINE, jet-like flow due to luminal narrowing was observed which increased prior to rupture and resulted in focal neighbored regions of high and low wall shear stress (WSS). Prior to rupture, aortic diameter at the rupture site increased lowering WSS at the entire aortic circumference. Concurrently, WSS and turbulence increased locally above the entry tear at the inner aortic curvature. Turbulent kinetic energy and WSS elevation in the downstream aorta demonstrated enhanced stress on the native aorta. Results of steady-state simulations were in good qualitative agreement with transient simulations. Steady-flow CFD simulations feasible at clinical time scales prior to aortic rupture reveal a hostile hemodynamic environment at the dSINE rupture site in agreement with lengthy transient simulations. Consequently, our developed approach may be of value in treatment planning where a fast assessment of the local hemodynamic environment is essential.

Active flow control of the airflow of a ship at yaw

This paper implements the steady Coanda effect active flow control (AFC) on the Chalmers ship model (CSM) to study its influence on the ship’s side force and airwake under the yaw effect. The study is conducted numerically using Large Eddy Simulation (LES) with Wall-Adapting Local-Eddy Viscosity (WALE) model. Numerical methods are validated by the experimental data acquired from the baseline CSM under 10∘ port-side wind. The model with AFC is created by modifying the square-shaped hanger base to the Coanda surface and added with injection slots along the base’s roof edge and two side edges. The results show that the base-shape modification significantly alters the vortex structure on deck from z-direction vortex (ZV) to streamwise vortex (SV), and the steady Coanda effect with a momentum coefficient (Cμ) of 0.02 further enhances the SV with the removal of port-side vortex (PV). The side force and yaw moment are reduced by 5.27% and 7.97%, respectively in the AFC case due to the reduction of port-side (windward) ship-surface pressure. Furthermore, the current AFC can suppress the low-speed region and alleviate the velocity gradient in the lateral direction, which mitigates the regions of high TKE (turbulent kinetic energy) and high shear stress along the port-side deck.

Cleaning of simple cohesive soil layers in a radial flow cell

A radial flow cell (RFC, inlet tube radius 0.95 mm and gaps of 1–4.2 mm) was used to investigate the removal of two thin (80–230 µm thick) model soil layers from glass and 304 stainless steel substrates by the flow of water at 20 °C. Under the flow conditions employed (Reynolds numbers 200–1400), inertial effects give rise to recirculation zones and regions of high shear stress on the bottom, soiled plate. The soils were dried layers of (i) instant coffee and (ii) a domestic abrasive cleaning suspension comprising fine particulates in a soluble matrix. Cleaning data exhibited a constant local erosion rate which varied strongly with radial position. For both soils, cleaning involved the growth of a circular cleaned region and redeposition of particulate matter in a ring at locations close to the foot of the secondary recirculation zone predicted by 2D axisymmetric CFD simulations. Removal beyond this location was observed with the coffee layers, indicating that cleaning for this soil was controlled primarily by simple diffusion mechanisms. The effect of channel aspect ratio and flow rate on the location of recirculation zones and shear stress distributions was investigated. The local cleaning rate in these steady flows was not linked simply to local wall shear stress.

Fluorescent Nanodiamonds for Detecting Free-Radical Generation in Real Time during Shear Stress in Human Umbilical Vein Endothelial Cells.

Free-radical generation is suspected to play a key role in cardiovascular diseases. Another crucial factor is shear stress. Human umbilical vein endothelial cells (HUVECS), which form the lining of blood vessels, require a physiological shear stress to activate many vasoactive factors. These are needed for maintaining vascular cell functions such as nonthrombogenicity, regulation of blood flow, and vascular tone. Additionally, blood clots form at regions of high shear stress within a blood vessel. Here, we use a new method called diamond magnetometry which allows us to measure the dynamics of free-radical generation in real time under shear stress. This quantum sensing technique allows free-radical detection with nanoscale resolution at the single-cell level. We investigate radical formation in HUVECs in a microfluidic environment under different flow conditions typically found in veins and arteries. Here, we looked into free-radical formation before, during, and after flow. We found that the free-radical production varied depending on the flow conditions. To confirm the magnetometry results and to differentiate between radicals, we performed conventional fluorescent reactive oxygen species (ROS) assays specific for superoxide, nitric oxide, and overall ROS.

Interplay between particulate fouling and its flow disturbance: Numerical and experimental studies

In this study, we examine the complex interplay between particulate fouling and flow disturbance via particle simulation and microfluidic experiments. From the experiment, it is found that a particle is transported and deposited following the streamline, which is affected by the existing particle deposit. Inspired by this microfluidic observation, we develop a novel coarse-grained particle simulation, which continuously updates the flow field by considering the particle deposits, to investigate particulate fouling more intensively. Four physical factors are proposed to describe quantitatively the particulate fouling: effective collision, hydrodynamic boundary layer, inaccessible zone, and high shear stress region. In the initial stage of fouling, a number of particles are deposited within the growing hydrodynamic boundary layer, where effective collisions are readily achieved. In the latter stage of fouling, the inaccessible zone and high shear stress region are expanded, which causes the particle deposition rate to decrease. This result provides a comprehensive understanding of the particulate fouling through the combination of the four physical factors.

Refining a numerical model for device-induced thrombosis and investigating the effects of non-Newtonian blood models

Thrombosis is one of the main causes of failure in device implantation. Computational thrombosis simulation is a convenient approach to evaluate the risk of thrombosis for a device. However, thrombosis is a complicated process involving multiple species and reactions. Application of a macroscopic, single-scale computational model for device-induced thrombosis is a cost-effective approach. The current study has refined an existing thrombosis model, which simulated thrombosis by tracing four species in blood: non-activated platelets, activated platelets, surface adherent platelets, and ADP. Platelets are activated mechanically by shear stress, and chemically by ADP. Platelet adhesion occurs on surfaces with low wall shear stress with platelet aggregation inhibited in regions of high shear stress. The study improves the existing thrombosis model by: 1) Modifying the chemical platelet activation function so that ADP activates platelets; 2) Modifying the function describing thrombus deposition and growth to distinguish between thrombus deposition on wall surfaces and thrombus growth on existing thrombus surfaces; 3) Modifying the thrombus breakdown function to allow for thrombus breakdown by shear stress; 4) Modeling blood flow as non-Newtonian. The results show that the inclusion of ADP and the use of a non-Newtonian model improve agreement with experimental data.

Mechanical Forces Impacting Cleavage of Von Willebrand Factor in Laminar and Turbulent Blood Flow

Von Willebrand factor (VWF) is a large multimeric hemostatic protein. VWF is critical in arresting platelets in regions of high shear stress found in blood circulation. Excessive cleavage of VWF that leads to reduced VWF multimer size in plasma can cause acquired von Willebrand syndrome, which is a bleeding disorder found in some heart valve diseases and in patients receiving mechanical circulatory support. It has been proposed that hemodynamics (blood flow) found in these environments ultimately leads to VWF cleavage. In the context of experiments reported in the literature, scission theory, developed for polymers, is applied here to provide insight into flow that can produce strong extensional forces on VWF that leads to domain unfolding and exposure of a cryptic site for cleavage through a metalloproteinase. Based on theoretical tensile forces, laminar flow only enables VWF cleavage when shear rate is large enough (>2800 s−1) or when VWF is exposed to constant shear stress for nonphysiological exposure times (>20 min). Predicted forces increase in turbulence, increasing the chance for VWF cleavage. These findings can be used when designing blood-contacting medical devices by providing hemodynamic limits to these devices that can otherwise lead to acquired von Willebrand syndrome.

Development of a test method to investigate mode II fracture and dissection of arteries

The current study presents the development and implementation of a bespoke experimental technique to generate and characterise mode II crack initiation and propagation in arterial tissue. The current study begins with a demonstration that lap-shear testing of arterial tissue results in mixed mode fracture, rather than mode II. We perform a detailed computational design of a bespoke experimental method (which we refer to as a shear fracture ring test (SFRT)) to robustly and repeatably generate mode II crack initiation and propagation in arteries. This method is based on generating a localised region of high shear adjacent to a cylindrical loading bar. Placement of a radial notch in this region of high shear stress is predicted to result in a kinking of the crack during a mode II initiation and propagation of the crack over a long distance in the circumferential (c)-direction along the circumferential-axial (c-a) plane. Fabrication and experimental implementation of the SFRT on excised ovine aorta specimens confirms that the bespoke test method results in pure mode II initiation and propagation. We demonstrate that the mode II fracture strength along the c-a plane is eight times higher than the corresponding mode I strength determined from a standard peel test. We also calibrate the mode II fracture energy based on our measurement of crack propagation rates. The mechanisms of fracture uncovered in the current study, along with our quantification of mode II fracture properties have significant implications for current understanding of the biomechanical conditions underlying aortic dissection.

Near- and far-field structure of shallow mixing layers between parallel streams

Abstract

Macrophage accumulation within coronary arterial wall in diabetic patients with acute coronary syndrome: a study with in-vivo intravascular imaging modalities

Background and aimsMacrophage accumulation in arteriosclerotic plaque of coronary arteries is involved in plaque destabilization. Atherosclerosis has been known to be progressive in patients with type 2 diabetes mellitus (DM). This study compared the features of 3-dimensional (3D) spatial distribution of macrophage accumulation within coronary artery wall between acute coronary syndrome (ACS) patients with DM (n = 20) and those without (non-DM, n = 20) by using intravascular ultrasound (IVUS) and optical coherence tomography (OCT).MethodsThe OCT-derived macrophage accumulation was measured within the proximal left anterior-descending artery. This measurement was performed for the whole vessel segment of interest, higher shear stress region (flow divider side) and lower shear stress region (the opposite side).ResultsNormalized macrophage accumulation per unit length of the whole segment of interest was significantly larger in ACS patients with DM than without. In non-DM patients, macrophage density per IVUS-derived plaque volume was significantly higher in high shear stress region compared to low shear stress region, however, there was no significant difference between the two regions in DM patients. The macrophage density in the low shear stress region was significantly higher in the DM group than in the non-DM group. A multivariate analysis showed that the presence of DM was a major determinant for macrophage distribution.ConclusionsMacrophage accumulation was more abundant and homogeneous within coronary arterial wall in DM patients with ACS compared to non-DM patients, suggesting that plaque destabilization may occur more widely throughout coronary wall in DM patients.

Influence of long-term thermal ageing (＞10,000h) on low stress creep deformation of P92 heat-resistant steels

ABSTRACTMicrostructure evolution of P92 heat-resistant steel during long-term service has high influence on the creep property. In this paper, P92 heat-resistant steel samples were aged at 898 K for 10,000 h, 15,000 h and 25,000 h. The effect of thermal ageing on microstructure was analysed through optical microscopy and scanning electron microscopy. The creep behaviour at specific shear stresses (14.7 MPa, 20 MPa, 25.6 MPa, 41 MPa) was investigated through helicoid-spring creep testing. The relationship between microstructure evolution and creep properties was discussed. The results demonstrated that the basic morphology of P92 heat-resistant steel samples following long-term thermal ageing was martensite strip structure, while the precipitated phases were MX, Laves and M23C6. The coarsening of M23C6 carbide was slower than Laves phase during thermal ageing of 15,000 ~ 25,000 h. Precipitation strengthening for the grain interior was corresponding to the Orowan mechanism. In the regions of high shear stresses (14.7 MPa, 20 MPa), the Laves phases played an important role in dislocation motion inhibition at grain boundaries, while in the regions of low shear stresses (25.6 MPa, 41 MPa), the M23C6 phases played a significant role in dislocation motion inhibition at grain boundaries. Moreover, the main precipitate which affected the creep of P92 heat-resistant steels under the lower shear stresses of 14.7 MPa and 20 MPa was Laves phase and the higher shear stresses of 25.6 MPa and 41 MPa was mainly M23C6 phase.

Hemolysis associated with Impella heart pump positioning: In vitro hemolysis testing and computational fluid dynamics modeling.

Short-term mechanical circulatory support devices provide temporary hemodynamic support in heart failure and are increasingly used to enable recovery or as a bridge to decision. Blood damage with mechanical circulatory support devices is influenced by many factors, including the magnitude and duration of shear stress and obstruction to blood flow. This study aimed to evaluate the effects of the Impella CP® heart pump positioning on hemolysis using in vitro hemolysis testing and computational fluid dynamics modeling. The in vitro hemolysis testing was conducted per the recommended Food and Drug Administration and American Society for Testing and Materials guidelines. The bench hemolysis testing and computational fluid dynamics simulation analysis were performed for both normal operating (outlet unobstructed) and outlet-obstructed condition of Impella CP (mimicking outlet on the aortic valve due to improper positioning). The modified index of hemolysis was 2.78 ± 0.69 at normal operating conditions compared to 18.7 ± 7.8 when the Impella CP outlet was obstructed (p = 0.002). Computational fluid dynamics modeling showed about three times increase in exposure time to regions of high shear stress when the Impella CP outlet was obstructed compared to unobstructed condition, thus supporting the experimental observations. Based on these results, it is recommended to ensure proper placement of Impella CP via regular monitoring using echocardiographic guidance or other methods to minimize the risk of hemolysis associated with an obstructed outflow.

Computational fluid dynamic analysis of physical forces playing a role in brain organoid cultures in two different multiplex platforms

BackgroundOrganoid cultivation in suspension culture requires agitation at low shear stress to allow for nutrient diffusion, which preserves tissue structure. Multiplex systems for organoid cultivation have been proposed, but whether they meet similar shear stress parameters as the regularly used spinner flask and its correlation with the successful generation of brain organoids has not been determined.ResultsHere we used computational fluid dynamics (CFD) to simulate two multiplex culture conditions: steering plates on an orbital shaker and the use of a previously described bioreactor. The bioreactor had low speed and high shear stress regions that may affect cell aggregate growth, depending on volume, whereas the computed variables of the steering plates were closer to those of the spinning flask.ConclusionOur protocol improves the initial steps of the standard brain organoid formation, and the produced organoids displayed regionalized brain structures, including retinal pigmented cells. Overall, we conclude that suspension culture on orbital steering plates is a cost-effective practical alternative to previously described platforms for the cultivation of brain organoids for research and multiplex testing.

Computational Fluid Dynamics Characterization of Two Patient-Specific Systemic-to-Pulmonary Shunts before and after Operation.

Studying the haemodynamics of the central shunt (CS) and modified Blalock–Taussig shunt (MBTS) benefits the improvement of postoperative recovery for patients with an aorta-pulmonary shunt. Shunt configurations, including CS and MBTS, are virtually reconstructed for infants A and B based on preoperative CT data, and three-dimensional models of A, 11 months after CS, and B, 8 months after MBTS, are reconstructed based on postoperative CT data. A series of parameters including energy loss, wall shear stress, and shunt ratio are computed from simulation to analyse the haemodynamics of CS and MBTS. Our results showed that the shunt ratio of the CS is approximately 30% higher than the MBTS and velocity distribution in the left pulmonary artery (LPA) and right pulmonary artery (RPA) was closer to a natural development in the CS than the MBTS. However, energy loss of the MBTS is lower, and the MBTS can provide more symmetric pulmonary artery (PA) flow than the CS. With the growth of infants A and B, the shunt ratio of infants was decreased, but maximum wall shear stress and the distribution region of high wall shear stress (WSS) were increased, which raises the probability of thrombosis. For infant A, the preoperative abnormal PA structure directly resulted in asymmetric growth of PA after operation, and the LPA/RPA ratio decreased from 0.49 to 0.25. Insufficient reserved length of the MBTS led to traction phenomena with the growth of infant B; on the one hand, it increased the eddy current, and on the other hand, it increased the flow resistance of anastomosis, promoting asymmetric PA flow.

Strength of Strained Two‐Phase Mixtures: Application to Rapid Creep and Stress Amplification in Subduction Zone Mélange

AbstractAseismic creep may occur by distributed deformation in mélange shear zones comprising weak matrix and stronger clast materials. Slow slip events or steady tectonic displacement can be distributed over <100‐m thick shear zones if weak matrix controls bulk shear zone deformation. We use 2‐D numerical models to quantify the rheology of moderately strained (shear strain <1.75) mélange for various volumetric proportions of competent clasts. Mélange deformation with <50% clasts is matrix dominated and can accommodate steady creep. At higher clast proportions mélange viscosity increases more than tenfold after small strains, because strong clasts interact and form force chains. Clast shear stress is amplified above the imposed shear stress, by a factor of <14 where force chains develop. Slow slip events may occur due to a temporary absence of force chains, while localized regions of high shear stress generate coincident fracturing and potentially tremor events.