The influence of geometric factors on the wall shear stress distribution in realistic human coronary arteries

Dissertação apresentada na Faculdade de Ciências e Tecnologia da Universidade Nova de Lisboa para a obtenção do grau de Mestre em Engenharia Biomédica. A presente dissertação foi desenvolvida no Erasmus Medical Center em Roterdão, Holanda

Static pressure and wall shear stress distributions in air flow in a seven wire-wrapped rod bundle

An experimental investigation is performed in a turbulent flow in a seven wire-wrapped rod bundle, mounted in an open air facility. Static pressure distributions are measured on central and peripheral rods. By using a Preston tube, the wall shear stress profiles are experimentally obtained along the perimeter of the rods. The geometric parameters of the test section are P/D=1.20 and H/D=15. The measuring section is located at L/D=40 from the air inlet. It is observed that the dimensionless static pressure and wall shear stress profiles are nearly independent of the Reynolds number and strongly dependent of the wire-spacer position, with abrupt variations of the parameters in the neighborhood of the wires.

Does the degree of coarctation of the aorta influence wall shear stress focal heterogeneity?

The development of atherosclerosis in the aorta is associated with low and oscillatory wall shear stress for normal patients. Moreover, localized differences in wall shear stress heterogeneity have been correlated with the presence of complex plaques in the descending aorta. While it is known that coarctation of the aorta can influence indices of wall shear stress, it is unclear how the degree of narrowing influences resulting patterns. We hypothesized that the degree of coarctation would have a strong influence on focal heterogeneity of wall shear stress. To test this hypothesis, we modeled the fluid dynamics in a patient-specific aorta with varied degrees of coarctation. We first validated a massively parallel computational model against experimental results for the patient geometry and then evaluated local shear stress patterns for a range of degrees of coarctation. Wall shear stress patterns at two cross sectional slices prone to develop atherosclerotic plaques were evaluated. Levels at different focal regions were compared to the conventional measure of average circumferential shear stress to enable localized quantification of coarctation-induced shear stress alteration. We find that the coarctation degree causes highly heterogeneous changes in wall shear stress.

Relation between plaque type, plaque thickness, blood shear stress, and plaque stress in coronary arteries assessed by X-ray Angiography and Intravascular Ultrasound

Purpose: Atheromatic plaque progression is affected, among others phenomena, by biomechanical, biochemical, and physiological factors. In this paper, the authors introduce a novel framework able to provide both morphological (vessel radius, plaque thickness, and type) and biomechanical (wall shear stress and Von Mises stress) indices of coronary arteries. Methods: First, the approach reconstructs the three-dimensional morphology of the vessel from intravascular ultrasound(IVUS) and Angiographic sequences, requiring minimal user interaction. Then, a computational pipeline allows to automatically assess fluid-dynamic and mechanical indices. Ten coronary arteries are analyzed illustrating the capabilities of the tool and confirming previous technical and clinical observations. Results: The relations between the arterial indices obtained by IVUS measurement and simulations have been quantitatively analyzed along the whole surface of the artery, extending the analysis of the coronary arteries shown in previous state of the art studies. Additionally, for the first time in the literature, the framework allows the computation of the membrane stresses using a simplified mechanical model of the arterial wall. Conclusions: Circumferentially (within a given frame), statistical analysis shows an inverse relation between the wall shear stress and the plaque thickness. At the global level (comparing a frame within the entire vessel), it is observed that heavy plaque accumulations are in general calcified and are located in the areas of the vessel having high wall shear stress. Finally, in their experiments the inverse proportionality between fluid and structural stresses is observed.

Relation between plaque type, plaque thickness, blood shear stress, and plaque stress in coronary arteries assessed by X-ray Angiography and Intravascular Ultrasound

Purpose: Atheromatic plaque progression is affected, among others phenomena, by biomechanical, biochemical, and physiological factors. In this paper, the authors introduce a novel framework able to provide both morphological (vessel radius, plaque thickness, and type) and biomechanical (wall shear stress and Von Mises stress) indices of coronary arteries. Methods: First, the approach reconstructs the three-dimensional morphology of the vessel from intravascular ultrasound(IVUS) and Angiographic sequences, requiring minimal user interaction. Then, a computational pipeline allows to automatically assess fluid-dynamic and mechanical indices. Ten coronary arteries are analyzed illustrating the capabilities of the tool and confirming previous technical and clinical observations. Results: The relations between the arterial indices obtained by IVUS measurement and simulations have been quantitatively analyzed along the whole surface of the artery, extending the analysis of the coronary arteries shown in previous state of the art studies. Additionally, for the first time in the literature, the framework allows the computation of the membrane stresses using a simplified mechanical model of the arterial wall. Conclusions: Circumferentially (within a given frame), statistical analysis shows an inverse relation between the wall shear stress and the plaque thickness. At the global level (comparing a frame within the entire vessel), it is observed that heavy plaque accumulations are in general calcified and are located in the areas of the vessel having high wall shear stress. Finally, in their experiments the inverse proportionality between fluid and structural stresses is observed.

Impact of local endothelial shear stress on neointima and plaque following stent implantation in patients with ST-elevation myocardial infarction: A subgroup-analysis of the COMFORTABLE AMI-IBIS 4 trial.

BACKGROUND Numerous studies have demonstrated an association between endothelial shear stress (ESS) and neointimal formation after stent implantation. However, the role of ESS on the composition of neointima and underlying plaque remains unclear. METHODS Patients recruited in the Comfortable AMI-IBIS 4 study implanted with bare metal stents (BMS) or biolimus eluting stents (BES) that had biplane coronary angiography at 13month follow-up were included in the analysis. The intravascular ultrasound virtual-histology (IVUS-VH) and the angiographic data were used to reconstruct the luminal surface, and the stent in the stented segments. Blood flow simulation was performed in the stent surface, which was assumed to represent the luminal surface at baseline, to assess the association between ESS and neointima thickness. The predominant ESS was estimated in 3-mm segments and was correlated with the amount of neointima, neointimal tissue composition, and with the changes in the underlying plaque burden and composition. RESULTS Forty three patients (18 implanted with BMS and 25 with BES) were studied. In both stent groups negative correlations were noted between ESS and neointima thickness in BMS (P<0.001) and BES (P=0.002). In BMS there was a negative correlation between predominant ESS and the percentage of the neointimal necrotic core component (P=0.015). In BES group, the limited neointima formation did not allow evaluation of the effect of ESS on its tissue characteristics. ESS did not affect vessel wall remodeling and the plaque burden and composition behind BMS (P>0.10) and BES (P>0.45). CONCLUSIONS ESS determines neointimal formation in both BMS and BES and affects the composition of the neointima in BMS. Conversely, ESS does not impact the plaque behind struts irrespective of stent type throughout 13months of follow-up.

Identification of vascular endothelial genes differentially responsive to fluid mechanical stimuli: cyclooxygenase-2, manganese superoxide dismutase, and endothelial cell nitric oxide synthase are selectively up-regulated by steady laminar shear stress.

Early atherosclerotic lesions develop in a topographical pattern that strongly suggests involvement of hemodynamic forces in their pathogenesis. We hypothesized that certain endothelial genes, which exhibit differential responsiveness to distinct fluid mechanical stimuli, may participate in the atherogenic process by modulating, on a local level within the arterial wall, the effects of systemic risk factors. A differential display strategy using cultured human endothelial cells has identified two genes, manganese superoxide dismutase and cyclooxygenase-2, that exhibit selective and sustained up-regulation by steady laminar shear stress (LSS). Turbulent shear stress, a nonlaminar fluid mechanical stimulus, does not induce these genes. The endothelial form of nitric oxide synthase also demonstrates a similar LSS-selective pattern of induction. Thus, three genes with potential atheroprotective (antioxidant, antithrombotic, and antiadhesive) activities manifest a differential response to distinct fluid mechanical stimuli, providing a possible mechanistic link between endothelial gene expression and early events in atherogenesis. The activities of these and other LSS-responsive genes may have important implications for the pathogenesis and prevention of atherosclerosis.

Critical bed shear stress and threshold of motion of maërl biogenic gravel

Critical bed shear stress for incipient motion has been determined for biogenic free-living coralline algae known as maërl. Maërl from three different sedimentary environments (beach, intertidal, and open marine) in Galway Bay, west of Ireland have been analysed in a rotating annular flume and linear flume. Velocity profile measurements of the benthic boundary layer, using an Acoustic Doppler Velocimeter, have been obtained in four different velocity experiments. The bed shear stress has been determined using three methods: Law of the Wall, Turbulent Kinetic Energy and Reynolds Stress. The critical Shields parameter has been estimated as a non-dimensional mobility number and the results have been compared with the Shields curve for natural sand. Maërl particles fall below this curve because its greater angularity allows grains to be mobilised easier than hydraulically equivalent particles. From previous work, the relationship between grain shape and the settling velocity of maërl suggests that the roughness is greatest for intertidal maërl particles. During critical shear stress determinations, beds of such rough particles exhibited the greatest critical shear stress probably because the particle thalli interlocked and resisted entrainment. The Turbulent Kinetic Energy methodology gives the most consistent results, agreeing with previous comparative studies. Rarely-documented maërl megaripples were observed in the rotating annular flume and are hypothesised to form at velocities ~10 cm s-1 higher than the critical threshold velocity, where tidal currents, oscillatory flow or combined-wave current interaction results in the preferential transport of maërl. A determination of the critical bed shear stress of maërl allows its mobility and rate of erosion and deposition to be evaluated spatially in subsequent applications to biological conservation management.

ACE as a Mechanosensor to Shear Stress Influences the Control of Its Own Regulation via Phosphorylation of Cytoplasmic Ser(1270)

Objectives: We tested whether angiotensin converting enzyme (ACE) and phosphorylation of Ser(1270) are involved in shear-stress (SS)-induced downregulation of the enzyme. Methods and Results: Western blotting analysis showed that SS (18 h, 15 dyn/cm(2)) decreases ACE expression and phosphorylation as well as p-JNK inhibition in human primary endothelial cells (EC). CHO cells expressing wild-type ACE (wt-ACE) also displayed SS-induced decrease in ACE and p-JNK. Moreover, SS decreased ACE promoter activity in wt-ACE, but had no effect in wild type CHO or CHO expressing ACE without either the extra-or the intracellular domains, and decreased less in CHO expressing a mutated ACE at Ser(1270) compared to wt-ACE (13 vs. 40%, respectively). The JNK inhibitor (SP600125, 18 h), in absence of SS, also decreased ACE promoter activity in wt-ACE. Finally, SS-induced inhibition of ACE expression and phosphorylation in EC was counteracted by simultaneous exposure to an ACE inhibitor. Conclusions: ACE displays a key role on its own downregulation in response to SS. This response requires both the extra- and the intracellular domains and ACE Ser(1270), consistent with the idea that the extracellular domain behaves as a mechanosensor while the cytoplasmic domain elicits the downstream intracellular signaling by phosphorylation on Ser(1270).

Shear Stress Induces Nitric Oxide-Mediated Vascular Endothelial Growth Factor Production in Human Adipose Tissue Mesenchymal Stem Cells

It has been demonstrated that human adipose tissue-derived mesenchymal stem cells (hASCs) enhance vascular density in ischemic tissues, suggesting that they can differentiate into vascular cells or release angiogenic factors that may stimulate neoangiogenesis. Moreover, there is evidence that shear stress (SS) may activate proliferation and differentiation of embryonic and endothelial precursor stem cells into endothelial cells (ECs). In this work, we investigated the effect of laminar SS in promoting differentiation of hASCs into ECs. SS (10 dyn/cm(2) up to 96 h), produced by a cone plate system, failed to induce EC markers (CD31, vWF, Flk-1) on hASC assayed by RT-PCR and flow cytometry. In contrast, there was a cumulative production of nitric oxide (determined by Griess Reaction) and vascular endothelial growth factor (VEGF; by ELISA) up to 96 h of SS stimulation ( NO(2)(-) in nmol/10(4) cells: static: 0.20 +/- 0.03; SS: 1.78 +/- 0.38, n = 6; VEGF in pg/10(4) cells: static: 191.31 +/- v35.29; SS: 372.80 +/- 46.74, n = 6, P < 0.05). Interestingly, the VEGF production was abrogated by 5 mM N(G)-L-nitro-arginine methyl ester (L-NAME) treatment (VEGF in pg/10(4) cells: SS: 378.80 +/- 46.74, n = 6; SS + L-NAME: 205.84 +/- 91.66, n = 4, P < 0.05). The results indicate that even though SS failed to induce EC surface markers in hASC under the tested conditions, it stimulated NO-dependent VEGF production.

Bed shear stress measurements in dam break and swash flows

A novel shear plate was used to make direct bed shear stress measurements in laboratory dam break and swash flows on smooth, fixed, impermeable beds. The pressure gradient due to the slope of the fluid free-surface across the plate was measured using pressure transducers. Surface elevation was measured at five locations using acoustic displacement sensors. Flow velocity was measured using an Acoustic-Doppler Velocimeter and calculated using the ANUGA inundation model. The measured bed shear stress at the dam break fluid tip for an initially dry, horizontal bed was close to twice that estimated using steady flow theory. The temporal variation of swash bed shear stress showed a large peak in landward directed stress at the uprush tip, followed by a rapid decay throughout the uprush flow interior. The peak seaward directed stress during the backwash phase was less than half that measured in the uprush. Close to the still water line, in the region of bore collapse and at the time of initial uprush, favourable pressure gradients were measured. In the lower swash region predominately weak adverse pressure gradients were measured.

Effect of shear stress on expression of a recombinant protein by Chinese hamster ovary cells

A flow chamber was used to impart a steady laminar shear stress on a recombinant Chinese hamster ovary (CHO) cell line expressing human growth hormone (hGH). The cells were subjected to shear stress ranging from 0.005 to 0.80 N m(-2). The effect of shear stress on the cell specific glucose uptake, cell specific hGH, and lactate productivity rates were calculated. No morphological changes to the cells were observed over the range of shear stresses examined. When the cells were subjected to 0.10 N m(-2) shear in protein-free media without Pluronic F-68, recombinant protein production ceased with no change in cell morphology, whereas control cultures were expressing hGH at 0.35 mug/10(6) cells/h. Upon addition of the shear protectants, Pluronic F-68 (0.2% [w/v]) or fetal bovine serum (1.0% [v/v] FBS), the productivity of the cells was restored. The effect of increasing shear stress on the cells in protein-free medium containing Pluronic F-68 was also investigated. Cell specific metabolic rates were calculated for cells under shear stress and for no-shear control cultures performed in parallel, with shear stress rates expressed as a percentage of those obtained for control cultures. Upon increasing shear from 0.005 to 0.80 N m(-2), the cell specific hGH productivity decreased from 100% at 0.005 N m(-2) to 49% at 0.80 N m(-2) relative to the no-shear control. A concurrent increase in the glucose uptake rate from 115% at 0.01 N m(-2) to 142% at 0.80 N m(-2), and decreased lactate productivity from 92% to 50%, revealed a change in the yield of products from glucose compared with the static control. It was shown that shear stress, at sublytic levels in medium containing Pluronic F-68, could decrease hGH specific productivity. (C) 2002 Wiley Periodicals, Inc.

Induction of CRP3/MLP expression during vein arterialization is dependent on stretch rather than shear stress

Aims Cysteine- and glycine-rich protein 3/muscle LIM-domain protein (CRP3/MLP) mediates protein-protein interaction with actin filaments in the heart and is involved in muscle differentiation and vascular remodelling. Here, we assessed the induction of CRP3/MLP expression during arterialization in human and rat veins. Methods and results Vascular CRP3/MLP expression was mainly observed in arterial samples from both human and rat. Using quantitative real time RT-PCR, we demonstrated that the CRP3/MLP expression was 10 times higher in smooth muscle cells (SMCs) from human mammary artery (h-MA) vs. saphenous vein (h-SV). In endothelial cells (ECs), CRP3/MLP was scarcely detected in either h-MA or h-SV. Using an ex vivo flow through system that mimics arterial condition, we observed induction of CRP3/MLP expression in arterialized h-SV. Interestingly, the upregulation of CRP3/MLP was primarily dependent on stretch stimulus in SMCs, rather than shear stress in ECs. Finally, using a rat vein in vivo arterialization model, early (1-14 days) CRP3/MLP immunostaining was observed predominantly in the inner layer and later (28-90 days) it appeared more scattered in the vessel layers. Conclusion Here we provide evidence that CRP3/MLP is primarily expressed in arterial SMCs and that stretch is the main stimulus for CRP3/MLP induction in veins exposed to arterial haemodynamic conditions.

Shear stress and sediment transport calculations for swash zone modelling

It is shown that the observed difference in sediment transporting efficiency by the swash uprush, compared with the downrush, could be mainly due to greater bed shear stress for a given velocity in the more abruptly accelerated uprush. The bed shear stress generated by an arbitrary free stream velocity time series is modelled in terms of usual wave boundary layer models plus a phase lead (phi(tau) of the bed shear stress compared with the free stream velocity at the peak frequency. With this approach, the total transport amounts in uprush and downrush can be modelled satisfactorily with the same sediment transport formula, without the need for different uprush and downrush coefficients. While the adaptation of sediment transport formulae from steady flow can thus lead to the right total amounts of sediment moved by this method, the timing of the instantaneous sediment transport rates are probably not accurately modelled due to the highly unsteady nature of the swash and the presence of pre-suspended sediment in the uprush. Nevertheless, the proposed method is a useful intermediate step before we have a complete understanding of sediment transport under very rapid accelerations and of the relative contribution of pre-suspended sediment to the onshore sediment transport in swash zones. (C) 2002 Published by Elsevier Science B.V.

Shear stress and sediment transport calculations for sheet flow under waves

A simple method is provided for calculating transport rates of not too fine (d(50) greater than or equal to 0.20 mm) sand under sheet flow conditions. The method consists of a Meyer-Peter-type transport formula operating on a time-varying Shields parameter, which accounts for both acceleration-asymmetry and boundary layer streaming. While velocity moment formulae, e.g.., = Constant x calibrated against U-tube measurements, fail spectacularly under some real waves (Ribberink, J.S., Dohmen-Janssen, C.M., Hanes, D.M., McLean, S.R., Vincent, C., 2000. Near-bed sand transport mechanisms under waves. Proc. 27th Int. Conf. Coastal Engineering, Sydney, ASCE, New York, pp. 3263-3276, Fig. 12), the new method predicts the real wave observations equally well. The reason that the velocity moment formulae fail under these waves is partly the presence of boundary layer streaming and partly the saw-tooth asymmetry, i.e., the front of the waves being steeper than the back. Waves with saw-tooth asymmetry may generate a net landward sediment transport even if = 0, because of the more abrupt acceleration under the steep front. More abrupt accelerations are associated with thinner boundary layers and greater pressure gradients for a given velocity magnitude. The two real wave effects are incorporated in a model of the form Q(s)(t) = Q(s)[theta(t)] rather than Q(S)(t) = Q(S)[u(infinity)(t)], i.e., by expressing the transport rate in terms of an instantaneous Shields parameter rather than in terms of the free stream velocity, and accounting for both streaming and accelerations in the 0(t) calculations. The instantaneous friction velocities u(*)(t) and subsequently theta(t) are calculated as follows. Firstly, a linear filter incorporating the grain roughness friction factor f(2.5) and a phase angle phi(tau) is applied to u(infinity)(t). This delivers u(*)(t) which is used to calculate an instantaneous grain roughness Shields parameter theta(2.5)(t). Secondly, a constant bed shear stress is added which corresponds to the streaming related bed shear stress -rho ($) over bar((u) over tilde(w) over tilde)(infinity) . The method can be applied to any u(infinity)(t) time series, but further experimental validation is recommended before application to conditions that differ strongly from the ones considered below. The method is not recommended for rippled beds or for sheet flow with typical prototype wave periods and d(50) < 0.20 turn. In such scenarios, time lags related to vertical sediment movement become important, and these are not considered by the present model. (C) 2002 Elsevier Science B.V. All rights reserved.