Wall Shear Stress Distribution Research Articles

Numerical simulation of shelter effect assessment for single-row windbreaks on the periphery of oasis farmland

Dust storms, resulting from aeolian erosion, pose significant environmental hazards, while farmland is one of the main sources of dust storm release. An effective strategy to mitigate surface wind speed and curb dust emissions involves the establishment of windbreaks on the periphery of oasis farmland. This study conducts a series of numerical simulations using computational fluid dynamics (CFD) to explore the airflow fields around windbreaks with diverse systematic structural parameters, encompassing porosity, planting spacing, and fence effects. The main findings are as follows: (1) When the vegetation porosity is consistent (e.g., porosity α = 0.7, 0.8, and 0.9), exact geometry results can effectively reflect the distribution of wall shear stress, while the porous medium model overlooks these details. (2) The “Venturi effect” contributes to the acceleration of surface erosion and improper planting spacing results in an elevation of near-surface velocity. Planting spacing of 0.5 m demonstrates superior wind speed reduction performance, mitigating aeolian erosion and accumulation. (3) When the fence is positioned at l = 5h (h represents the height of the windbreaks), the flow field around the windbreaks is minimally influenced. The optimal placement distance for fences should be close to the windbreaks, featuring minimal porosity (l = 0 h, α = 0.1), extending the shelter distance from 3 to 4 h to 5–6 h. The research results can provide a theoretical basis for optimizing the plant configuration of biological desertification control and soil erosion control measures.

Experimental study on hemodynamics of an end-to-side anastomosis

A three-dimensional and three-component velocity measurement on the flow field in a 45° end-to-side anastomosis model is conducted to investigate the hemodynamics, which is an important factor to the intimal hyperplasia formation and graft failure after surgery. Thanks to the advanced volumetric measurement technology of tomographic particle image velocimetry, the recirculation zone, low-speed region, and the spiral flow structures can be visualized. As a result, the flow field of three cases with the local maximum velocity of 0.15, 0.8, and 1.4 m/s are visible and the inlet velocity profile tends to be skewed as the flow rate increases. The mean vorticity contours indicate that the positive vortex center rotates 6.47°, 50.23°, and 90.4° and the negative vortex center rotates 20.44°, 15.73°, and 68.47°, respectively, in three cases. The instantaneous vortex structures identified by the λci criterion demonstrate two large-scale vortex structures in the distal section. The two vortices have the tendency to intertwine while one of them decays earlier. The wall shear stress (WSS) distributions on the entire model with the local maximum of 0.8, 5.8, and 13.8 Pa in three cases have been quantitatively achieved. The abnormal WSS and WSS gradient can help localize risk areas and understand the intimal hyperplasia formation. A detailed illustration of hemodynamics inside the 45° end-to-side anastomosis model has been provided, which demonstrates more comprehensive large-scale flow structures and abnormal WSS regions. Combined with the information of flow structures and WSS distribution, the understanding of the hemodynamics in the anastomosis can be strengthened.

Enhancing Wind Turbine Blade Preventive Maintenance Procedure through Computational Fluid Dynamics-Based Prediction of Wall Shear Stress

Wind turbine blades are essential parts of wind energy systems and are frequently exposed to harsh environmental elements, such as strong winds, turbulence, and corrosive atmospheric elements. Over time, these circumstances may result in serious harm to blades, such as delamination and erosion, which may negatively affect the wind turbine’s functionality and durability. Accurate prediction of various types of damage is crucial to improve the toughness and lifespan of wind turbine blades and to maximize the overall effectiveness of wind energy systems. This article presents a novel computational fluid dynamics (CFDs)-based method for analyzing the distribution of wall shear stress on turbine blades, aimed at publicizing the yearly maintenance procedure. The investigation results from the CFDs, when compared with the current situation in a wind turbine farm in Thailand, confirmed that our wall shear stress modeling accurately predicted wind turbine damage. A maximum wall shear stress level higher than 5.00 Pa in the case of PA 90°, incoming air velocity 10.00 m/s, and 15 rpm was the main contribution to presenting the erosion and delamination from current drone inspection in wind turbine farms. In conclusion, these findings demonstrated the potential of using CFDs to predict wind turbine blade delamination and erosion, thereby significantly contributing to the development of specific and accurate yearly preventive maintenance. The proposed CFDs-based approach should serve as a sustainability tool for local human development, benefiting wind turbine engineers and operating technicians by providing them with a deeper understanding of the local flow conditions and wall shear stress distribution along wind turbine blades. This enables them to make informed decisions regarding blade design and maintenance.

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.

Investigation of the Influence of Flow Pattern on the Internal Corrosion in Horizontal–Vertical Downward Pipe With Oil–Water Flow

Abstract In the petroleum and petrochemical industries, oil–water flow is widespread inside the pipes. The existence of water often results in internal corrosion in the horizontal–vertical downward pipe when water contacts the pipe wall. Surface wetting behavior and wall shear stress (WSS) are two important factors affecting corrosion procedure, which are governed by the flow patterns. With the propose to mitigate corrosion, focus shall be concentrated on the impact of flow pattern toward corrosion. In this work, the flow regime with oil–water flow in the horizontal–vertical pipe is investigated by computational fluid dynamics simulations. The cases with different mixture velocity (0.1 m/s–2.2 m/s) and different water cut (3%–40%) are investigated. The key discovery in this paper is that five types of flow patterns can be identified based on the multiphase flow in the horizontal–vertical pipe, which is rarely reported in recent work. According to the results of the surface wetting status and wall shear stress distribution, the severe corrosion area is predicted and classified into five types. The inside wall of elbows and the outside wall of vertical pipes are the area's most susceptible to corrosion, and the results are well in line with the on-site data.

Investigation of Rupture Risk of Thoracic Aortic Aneurysms via Fluid–Structure Interaction and Artificial Intelligence Method

AbstractPatient-specific studies on vascular flows have significantly increased for hemodynamics due to the need for different observation techniques in clinical practice. In this study, we investigate aortic aneurysms in terms of deformation, stress, and rupture risk. The effect of Ascending Aortic Diameter (AAD) was investigated in different aortic arches (19.81 mm, 42.94 mm, and 48.01 mm) via Computational Fluid Dynamics (CFD), Two-way coupling Fluid–Structure Interactions (FSI) and deep learning. The non-newtonian Carreau viscosity model was utilized with patient-specific velocity waveform. Deformations, Wall Shear Stresses (WSSs), von Mises stress, and rupture risk were presented by safety factors. Results show that the WSS distribution is distinctly higher in rigid cases than the elastic cases. Although WSS values rise with the increase in AAD, aneurysm regions indicate low WSS values in both rigid and elastic artery solutions. For the given AADs, the deformations are 2.75 mm, 6. 82 mm, and 8.48 mm and Equivalent von Mises stresses are 0.16 MPa, 0.46 MPa, and 0.53 MPa. When the rupture risk was evaluated for the arteries, the results showed that the aneurysm with AAD of 48.01 mm poses a risk up to three times more than AAD of 19.81 mm. In addition, an Artificial neural network (ANN) method was developed to predict the rupture risk with a 98.6% accurate prediction by numerical data. As a result, FSI could indicate more accurately the level of rupture risk than the rigid artery assumptions to guide the clinical assessments and deep learning methods could decrease the computational costs according to CFD and FSI.

Numerical modeling of the two-dimensional free and impinging jets with solid and erosion bed using FLOW-3D

ABSTRACT This study simulates the flow properties of a plane-layered jet in free and impinging jets on a horizontal surface using Flow-3D and validates them with previous research and experimental data. The flow characteristics, such as velocity and shear stress distribution, are compared with theoretical relationships. The simulations vary the outlet velocities for the free jet and the surface types and slip conditions for the impinging jet. The results show that the free jet flow follows the self-similar theory, except for a faster decay of numerical velocities at the end of the jet’s path. The results also reveal a virtual origin (λ) before the nozzle outlet, ranging from 7.6 to 8.4 mm. The impinging jet results match the experimental wall shear stress distribution diagram, with the peak at y/H = 0.14. The transverse velocity profiles and the wall shear stress distribution of the impinging jet are influenced by the slip condition and the surface type, but the difference between smooth and granular surfaces is less than 3%. This study offers valuable insights into the behaviour of free and impinging jets on a horizontal surface and shows the effectiveness of the Flow-3D numerical model.

The interplay between bacterial biofilms, encrustation, and wall shear stress in ureteral stents: a review across scales

Ureteral stents are hollow tubes that are inserted into the ureter to maintain the flow of urine from the kidney to the bladder. However, the use of these indwelling stents is associated with potential complications. Biofilm, an organized consortium of bacterial species embedded within a self-producing extracellular matrix, can attach to the outer and inner surfaces of ureteral stents. Furthermore, encrustation - defined as the buildup of mineral deposits on the stent surface - can occur independently or in parallel with biofilm formation. Both phenomena can cause stent obstruction, which can lead to obstructive pyelonephritis and make stent removal difficult. Understanding the influence of flow on the development of biofilm and encrustation and the impact of small mechanical environmental changes (e.g., wall shear stress distribution) is key to improve the long-term performance of stents. Identifying the optimal stent properties to prevent early bacterial attachment and/or crystal deposition and their growth, would represent a breakthrough in reducing biofilm-/encrustation-associated complications. This review identifies the most prevalent bacterial strains and crystal types associated with ureteral stents, and the process of their association with the stent surface, which often depends on patient comorbidities, stent material, and indwelling time. Furthermore, we focus on the often-overlooked role of fluid dynamics on biofilm and encrustation development in ureteral stents, across a range of physical scales (i.e., from micro- to macro-scale) with the aim of providing a knowledge base to inform the development of safer and more effective ureteral stents.

Influence of Wavy Arteries and Veins on Hemodynamic Characteristics: A Numerical Study

The present work is focused on the study of hemodynamic characteristics for tortuous arteries/veins. Tortuosity in arteries/veins is defined by introducing waviness in the wall of the tube. Analysis is further extended for bifurcated veins with and without wavy walls. Waviness is defined by two geometric parameters; pitch and depth of the wave. Four different combinations of pitch and depth are studied and compared with a plain straight wall. The present study is carried out numerically by using a computational fluid dynamics tool. Hemodynamics for a steady flow of blood is investigated through pressure, velocity, and wall shear stress distribution. Waviness in the wall of arteries/veins creates a recirculation zone at the crest and trough of the wall. Occurrence of the recirculation zone leads to reduction in velocity which in turn reduces wall shear stress. Variation in the magnitude of the velocity and corresponding wall shear stress at the crest and trough of the wavy wall depends on the pitch and depth of the artery/veins (tube).

Effects of hematocrit levels on flow structures and stress levels in the healthy and diseased carotid arteries

Flow structures and wall shear stress (WSS) distribution in human arteries depend strongly on blood rheology. Blood viscosity undergoes significant alterations due to varying hematocrit (Hct) levels. Hematocrit levels can change from the normal level (45%) due to numerous disorders such as leukemia, lymphoma, and congenital heart disease. However, the effects of Hct levels on the flow features and WSS levels in large stenotic arteries are not well reported in literature. The present study computationally assesses the variation in flow features in healthy and stenosed carotid arteries due to changes in blood Hct levels. We have reported variations in axial and secondary flow characteristics, vortex structures, and WSS-based parameters for different Hct levels. Specifically, we consider high (65%) and low (25%) Hct levels in addition to a normal Hct level. Our results reveal that both high and low Hct levels have adverse effects on the hemodynamic features inside a healthy carotid artery model. The low Hct level displays a much more widespread low time-averaged WSS (TAWSS) area, while the high Hct shows an increased oscillatory shear index (OSI) zone at the carotid sinus wall. In stenotic arteries, low Hct enhances the spread of the low TAWSS zone in the post-stenotic carotid sinus, whereas high Hct increases the risk of plaque rupture. The high Hct level also increases the spread of high OSI zone at the inner wall of internal carotid artery for the moderately stenosed artery.

Computational fluid dynamics study on three‐dimensional polymeric scaffolds to predict wall shear stress using machine learning models for bone tissue engineering applications

AbstractGeometrical patterns and dimensions of the polymeric scaffold play a major role in controlling the degradation and mechanical stimuli for osteogenic differentiation. Wall shear stress (WSS) analysis of scaffold provides a better understanding of the body fluid flow dynamics. A computational fluid dynamics (CFD) study was carried out to understand velocity profile and WSS distribution when the strands are arranged in rectangular and triangular pitch for the different strand diameters and spacing. The number of scaffold surfaces with less than 30 mPa and maximum and average WSS was estimated to check the suitability of the scaffold for loading stem cells. This situation is favorable to induce osteogenic activity and cell viability. Higher spacing/pitch between the strands increases the chances of scaffold surface having WSS less than 30 mPa. When the spacing and diameter are smaller, there is no significant variation in WSS and pressure drop between rectangular and triangular pitch arrangement is observed. Machine learning (ML) models were developed to predict WSS distribution and to reduce the computational cost involved in solving the Navier–Stokes equation. XG Boost and support vector machine (SVM) models outperform the other models in predicting the WSS with high R2 and five‐fold cross‐validation accuracy and are helpful in predicting the optimal design parameters of a three‐dimensional scaffold.

Positive Pressure and Negative Pressure Irrigation Dynamics with Different Needle Designs Using Computational Fluid Dynamics

This study aimed to investigate the irrigation dynamics of the positive pressure side-vented (SV) needle, EndoVac (micropores) needle and modified apical negative pressure (mANP) open-ended needle using computational fluid dynamics (CFD). A simulation of a prepared root canal (conical frustum) of 15 mm length with an apical diameter of 0.40 mm following Protaper F4 apical preparation was created using three-dimensional (3D) CAD software. The 3D simulated needle of SV 30G needle, EndoVac with micropores needle and mANP, 30G flat open-ended needle were also created. The irrigation dynamics were evaluated through transient CFD simulations. In addition, the irrigation dynamics of mANP at 0.2 mm, 0.5 mm, and 1.0 mm short from the working length were also assessed. The EndoVac and mANP showed negative apical static pressure and streamline patterns able to reach the apical region, thus indicating negligible extrusion. Meanwhile, SV showed positive apical static pressure and almost nonexistent streamlines beyond the needle tip. The SV showed the highest wall shear stress (WSS) magnitude of 1030Pa whereas Endovac (161 Pa) and mANP1 (258 Pa). However, SV revealed lower average WSS (10 Pa) compared to mANP1 (13 Pa) and mANP2 (11 Pa). This is due to SV developed a localised maximum WSS opposite the open vent area only therefore, uneven distribution of WSS. The EndoVac system developed a localised maximum WSS in the pair of micropores furthest away from the apical. CFD analysis of the EndoVac, mANP and SV showed different technique approach, needle design and needle depths insertion affect the irrigation dynamics pattern and magnitude.

Numerical simulation of pulsatile blood flow through eccentric double stenosed carotid artery

Purpose This study aims to perform three-dimensional numerical computations for blood flow through a double stenosed carotid artery. Pulsatile flow with Womersley number (Wo) of 4.65 and Reynolds number (Re) of 425, based on the diameter of normal artery and average velocity of inlet pulse, was considered. Design/methodology/approach Finite volume method based ANSYS Fluent 20.1 was used for solving the governing equations of three-dimensional, laminar, incompressible and non-Newtonian blood flow. A high-quality grid with sufficient refinement was generated using ICEM CFD 20.1. The time-averaged flow field was captured to investigate the effect of severity and eccentricity on the lumen flow characteristics. Findings The results show that an increase in interspacing between blockages brings shear layer instability within the region between two blockages. The velocity profile and wall shear stress distribution are found to be majorly influenced by eccentricity. On the other hand, their peak magnitude is found to be primarily influenced by severity. Results have also demonstrated that the presence of eccentricity in stenosis would assist in flow development. Originality/value Variation in severity and interspacing was considered with a provision of eccentricity equal to 10% of diameter. Eccentricity refers to the offset between the centreline of stenosis and the centreline of normal artery. For the two blockages, severity values of 40% and 60% based on diameter reduction were permuted, giving rise to four combinations. For each combination, three values of interspacing in the multiples of normal artery diameter (D), viz. 4D, 6D and 8D were considered.

Optical measurement of the shear stress and velocity distribution in an idealized deglutition process

The rheological properties of food have an impact on physical interaction with various haptic sensors in the mouth during the chewing and swallowing process. Therefore, the investigations on food rheology are dedicated to improving the texture or mouthfeel of a product. It is advantageous to perform such rheological measurements in geometries and under conditions that reflect the flow conditions during the food distortion in the mouth cavity. Such investigations have been limited primarily to viscous or lumpy foodstuffs. Our study is dedicated to the compression of liquid foams in an idealized deglutition process. To reflect the flow conditions in such a process, we developed an experimental setup consisting of a quasi-two-dimensional replica of the palate and a movable tongue. During the idealized deglutition process, the optical measurements of bubble displacement and deformation were performed. The deglutition process was modelled by varying the compression speed of the tongue towards the static palate. The variation in mean bubble size and liquid fraction of foam for two tongue geometries with different roughness are available. It is observed that the velocity distribution of bubble motion during the deglutition corresponds to a rigid body flow, notwithstanding the surface roughness. The resulting wall-shear stress, however, was higher for the rough tongue plate for all foam bubble sizes and liquid fraction. The resulting velocity and wall-shear stress distributions can provide information about the haptic perception in the mouth during the swallowing process.

Impact of Cross-Limb Stent-Graft Configuration on Hemodynamics in Abdominal Aortic Aneurysm Interventional Therapy.

The cross-limb (CL) technique is a commonly used endovascular treatment for addressing unfavorable anatomical features in patients with abdominal aortic aneurysm (AAA). The configuration of CL stent-graft plays a critical role in determining the postoperative hemodynamic properties and physiological behaviors, which ultimately impact the efficacy and safety of endovascular AAA treatment. This study aims to investigate the relationship between hemodynamics and CL stent-graft configuration from a hemodynamic perspective. Five distinct geometric models of cross-limb (CL) stent-graft configurations were constructed by optimizing the real clinical computed tomography angiography (CTA) data. These models varied in main body lengths and cross angles and were used to perform numerical simulations to analyze various hemodynamic parameters. Flow pattern, distribution of wall shear stress (WSS)-related parameters, localized normalized helicity (LNH), pressure drop, and the displacement force of all models were examined in this paper. In patient-specific cases, helical flow and WSS increase with the main body. However, it also generated secondary flow in localized areas, leading to increased oscillation in the WSS direction. Notably, increasing the stent graft's main body length or decreasing the cross angle reduced the displacement force exerted on the stent-graft. Reducing the cross angle did not significantly alter the hemodynamic characteristics. In the clinical practice of CL deployment, it is crucial to carefully consider the stent-graft configuration and the patient specific to achieve optimal postoperative outcomes. This study provides valuable insights for guiding stent selection and treatment planning in patients with abdominal aortic aneurysm undergoing CL techniques, from a hemodynamic perspective.

Inter-corelab agreement of angiography-derived wall shear stress distribution

Abstract Aims Wall shear stress (WSS) estimated in models reconstructed from intravascular imaging and 3D-quantitative coronary angiography (QCA) data provides important prognostic information and enables the identification of high-risk lesions. However, these analyses are time-consuming and require expertise limiting WSS applications and adoption in clinical practice. Recently, a novel software has been developed for real-time computation of WSS and multidirectional WSS distribution. The aim of this study is to examine its inter-corelab reproducibility. Methods and results Sixty lesions (20 in coronary bifurcations) with a borderline negative fractional flow reserve were selected and processed using the CAAS Workstation WSS prototype (version 6.0 PIE MEDICAL) to estimate WSS and multi-directional WSS. Analysis was performed by two Core laboratories; the estimations of the two corelabs for the WSS in 3 mm segments across each reconstructed vessel were extracted and compared. To examine the concordance between inter-observer analysis of clinically relevant values, the scatterplot was divided into 4 quadrants with AS 61.3% and maximum WSS 8.24 Pa based on the previous publications and Cohen’s Kappa value was presented using these cut-off values. In total 700 segments (256 located in bifurcated vessels) were included in the present analysis. Average analysis time on the software was 219±45s in straight model and 291±61s in bifurcation model. A high intraclass correlation was noted for all the WSS metrics for the estimations of the two core labs irrespective of the presence (range: 0.90 – 0.92) or absence (range: 0.89 - 0.90) of a coronary bifurcation, while the ICC was good-moderate for the multidirectional WSS (range: 0.72-0.86). In addition, the bias between the WSS estimations was low and the limits of agreement were narrow in Bland-Altman analyses (Figure 1). When using the cut-off value of maximum WSS and % AS, a substantial agreement was noted between the two estimations of the two core labs; κ=0.77 for identifying maximum WSS >8.24Pa and κ=0.71 for identifying AS >61.3% (Figure 2). Conclusion The CAAS workstation enables reproducible 3D-QCA reconstruction and computation of WSS metrics. Further research is needed to explore its efficacy in predicting cardiovascular events.Figure 1Figure 2

Conditional analysis on extreme wall shear stress and heat flux events in compressible turbulent boundary layers

This study presents a comprehensive analysis on the extreme positive and negative events of wall shear stress and heat flux fluctuations in compressible turbulent boundary layers (TBLs) solved by direct numerical simulations. To examine the compressibility effects, we focus on the extreme events in two representative cases, i.e. a supersonic TBL of Mach number $M=2$ and a hypersonic TBL of $M=8$ , by scrutinizing the coherent structures and their correlated dynamics based on conditional analysis. As characterized by the spatial distribution of wall shear stress and heat flux, the extreme events are indicated to be closely related to the structural organization of wall streaks, in addition to the occurrence of the alternating positive and negative structures (APNSs) in the hypersonic TBL. These two types of coherent structures are strikingly different, namely the nature of wall streaks and APNSs are shown to be related to the solenoidal and dilatational fluid motions, respectively. Quantitative analysis using a volumetric conditional average is performed to identify and extract the coherent structures that directly account for the extreme events. It is found that in the supersonic TBL, the essential ingredients of the conditional field are hairpin-like vortices, whose combinations can induce wall streaks, whereas in the hypersonic TBL, the essential ingredients become hairpin-like vortices as well as near-wall APNSs. To quantify the momentum and energy transport mechanisms underlying the extreme events, we proposed a novel decomposition method for extreme skin friction and heat flux, based on the integral identities of conditionally averaged governing equations. Taking advantage of this decomposition method, the dominant transport mechanisms of the hairpin-like vortices and APNSs are revealed. Specifically, the momentum and energy transports undertaken by the hairpin-like vortices are attributed to multiple comparable mechanisms, whereas those by the APNSs are convection dominated. In that, the dominant transport mechanisms in extreme events between the supersonic and hypersonic TBLs are indicated to be totally different.

Implementing new computational methods for the study of JCT and SC inner wall basement membrane biomechanics and hydrodynamics

PurposeThe conventional aqueous outflow pathway, which includes the trabecular meshwork (TM), juxtacanalicular tissue (JCT), and the inner wall endothelium of Schlemm's canal (SC), regulates intraocular pressure (IOP) by controlling the aqueous humor outflow resistance. Despite its importance, our understanding of the biomechanics and hydrodynamics within this region remains limited. Fluid-structure interaction (FSI) offers a way to estimate the biomechanical properties of the JCT and SC under various loading and boundary conditions, providing valuable insights that are beyond the reach of current imaging techniques. MethodsIn this study, a normal human eye was fixed at a pressure of 7 mm Hg, and two radial wedges of the TM tissues, which included the SC inner wall basement membrane and JCT, were dissected, processed, and imaged using 3D serial block-face scanning electron microscopy (SBF-SEM). Four different sets of images were used to create 3D finite element (FE) models of the JCT and inner wall endothelial cells of SC with their basement membrane. The outer JCT portion was carefully removed as the outflow resistance is not in that region, leaving only the SCE inner wall and a few µm of the tissue, which does contain the resistance. An inverse iterative FE algorithm was then utilized to calculate the unloaded geometry of the JCT/SC complex at an aqueous humor pressure of 0 mm Hg. Then in the model, the intertrabecular spaces, pores, and giant vacuole contents were replaced by aqueous humor, and FSI was employed to pressurize the JCT/SC complex from 0 to 15 mm Hg. ResultsIn the JCT/SC complex, the shear stress of the aqueous humor is not evenly distributed. Areas proximal to the inner wall of SC experience larger stresses, reaching up to 10 Pa, while those closer to the JCT undergo lower stresses, approximately 4 Pa. Within this complex, giant vacuoles with or without I-pore behave differently. Those without I-pores experience a more significant strain, around 14%, compared to those with I-pores, where the strain is roughly 9%. ConclusionsThe distribution of aqueous humor wall shear stress is not uniform within the JCT/SC complex, which may contribute to our understanding of the underlying selective mechanisms in the pathway.

Numerical study of irreversibility analysis on pulsatile flow of blood through a ω-shaped stenotic artery

The mathematical modeling of blood flow in a stenotic blood vessel is very important for understanding the effects of various geometric and rheological parameters on the normal flow of blood. The present study is concerned with the heat and mass transfer in the blood flow of a non-Newtonian fluid through a ω-shaped stenotic arterial segment. The non-Newtonian behavior of blood is described mathematically by employing the constitutive equation of the Herschel-Bulkley fluid. The wall of the artery is considered to be rigid having ω-shaped constriction in its lumen. The mild stenosis assumption is used to simplify the fully coupled momentum, energy, and concentration equations as well as the constitutive equation of the Herschel-Bulkley fluid. Numerical solutions for fluid velocity, temperature, mass concentration, and entropy generation are evaluated by using the explicit finite difference scheme. The effects of various emerging parameters on the velocity, wall shear stress, flow rate, temperature, concentration, and entropy distributions have been analyzed in detail. A comparison between different fluid models has also been presented. It is found that the temperature in the Herschel-Bulkley fluid is marginally lower than that of the Newtonian, Power law, and Bingham fluids whereas an opposite trend is noticed in the case of mass concentration.

Wall shear stress and pressure patterns in aortic stenosis patients with and without aortic dilation captured by high-performance image-based computational fluid dynamics.

Spatial patterns of elevated wall shear stress and pressure due to blood flow past aortic stenosis (AS) are studied using GPU-accelerated patient-specific computational fluid dynamics. Three cases of moderate to severe AS, one with a dilated ascending aorta and two within the normal range (root diameter less than 4cm) are simulated for physiological waveforms obtained from echocardiography. The computational framework is built based on sharp-interface Immersed Boundary Method, where aortic geometries segmented from CT angiograms are integrated into a high-order incompressible Navier-Stokes solver. The key question addressed here is, given the presence of turbulence due to AS which increases wall shear stress (WSS) levels, why some AS patients undergo much less aortic dilation. Recent case studies of AS have linked the existence of an elevated WSS hotspot (due to impingement of AS on the aortic wall) to the dilation process. Herein we further investigate the WSS distribution for cases with and without dilation to understand the possible hemodynamics which may impact the dilation process. We show that the spatial distribution of elevated WSS is significantly more focused for the case with dilation than those without dilation. We further show that this focal area accommodates a persistent pocket of high pressure, which may have contributed to the dilation process through an increased wall-normal forcing. The cases without dilation, on the contrary, showed a rather oscillatory pressure behaviour, with no persistent pressure "buildup" effect. We further argue that a more proximal branching of the aortic arch could explain the lack of a focal area of elevated WSS and pressure, because it interferes with the impingement process due to fluid suction effects. These phenomena are further illustrated using an idealized aortic geometry. We finally show that a restored inflow eliminates the focal area of elevated WSS and pressure zone from the ascending aorta.