Computational fluid dynamics of the right ventricular outflow tract and of the pulmonary artery: a bench model of flow dynamics.

OBJECTIVES: The reconstruction of the right ventricular outflow tract (RVOT) with valved conduits remains a challenge. The reoperation rate at 5 years can be as high as 25% and depends on age, type of conduit, conduit diameter and principal heart malformation. The aim of this study is to provide a bench model with computer fluid dynamics to analyse the haemodynamics of the RVOT, pulmonary artery, its bifurcation, and left and right pulmonary arteries that in the future may serve as a tool for analysis and prediction of outcome following RVOT reconstruction. METHODS: Pressure, flow and diameter at the RVOT, pulmonary artery, bifurcation of the pulmonary artery, and left and right pulmonary arteries were measured in five normal pigs with a mean weight of 24.6 ± 0.89 kg. Data obtained were used for a 3D computer fluid-dynamics simulation of flow conditions, focusing on the pressure, flow and shear stress profile of the pulmonary trunk to the level of the left and right pulmonary arteries. RESULTS: Three inlet steady flow profiles were obtained at 0.2, 0.29 and 0.36 m/s that correspond to the flow rates of 1.5, 2.0 and 2.5 l/min flow at the RVOT. The flow velocity profile was constant at the RVOT down to the bifurcation and decreased at the left and right pulmonary arteries. In all three inlet velocity profiles, low sheer stress and low-velocity areas were detected along the left wall of the pulmonary artery, at the pulmonary artery bifurcation and at the ostia of both pulmonary arteries. CONCLUSIONS: This computed fluid real-time model provides us with a realistic picture of fluid dynamics in the pulmonary tract area. Deep shear stress areas correspond to a turbulent flow profile that is a predictive factor for the development of vessel wall arteriosclerosis. We believe that this bench model may be a useful tool for further evaluation of RVOT pathology following surgical reconstructions.

Low-pressure environment and remodelling of the forearm vein in Brescia-Cimino haemodialysis access.

BACKGROUND: The aim of the study was to determine which, and to what extent, haemodynamic parameters contribute to the remodelling of the venous limb of the Brescia-Cimino haemodialysis access. METHODS: The dimensions of the radial artery and the venous limb of the haemodialysis access were measured by an echo-tracking technique. In six ESRD patients undergoing primary arteriovenous fistula (AVF) formation, vessel diameter, wall thickness, blood pressure and blood flow were measured after the operation, and at 1 and 3 months follow-up. The contralateral forearm vessels in their native position served as baseline values for comparison. RESULTS: The diameter of the proximal antecubital vein progressively increased over the study period without reaching significant differences (4430, 5041 and 6620 microm at weeks 1, 4 and 12 respectively), whereas the intima-media thickness remained unchanged. The venous dilatation was associated with a reduction of the mean shear stress that culminated after the operation and progressively returned to normal venous values at 3 months (24.5 vs 10.4 dyne/cm(2), P<0.043). Thus the venous limb of the AVF undergoes eccentric hypertrophy as demonstrated by the increase in wall cross-sectional area (4.42 vs 6.32 mm(2) at week 1 vs week 12, P<0.028). At the time of the operation, the blood pressure in the AVF was 151+/-14/92.4+/-11 mmHg vs 49+/-19/24.5+/-6 mmHg (means+/-SEM) for the radial artery and the venous limb of the vascular access, respectively. One year after the operation the blood pressure in the venous limb had not changed: 42+/-14/25.3+/-7 mmHg (means+/-SEM). Under these conditions, the systolo-diastolic diameter changes observed in the radial artery and the antecubital vein were within a similar range at all time points: 56+/-17 vs 90+/-26 microm (means+/-SEM) at week 12. CONCLUSIONS: The increased circumferential stress resulting from the flow-mediated dilatation rather than the elevation of blood pressure appears to represent the main contributing factor to the eccentric hypertrophy of the venous limb of Brescia-Cimino haemodialysis access.

Role of hemodynamic forces in the ex vivo arterialization of human saphenous veins.

BACKGROUND: Human saphenous vein grafts are one of the salvage bypass conduits when endovascular procedures are not feasible or fail. Understanding the remodeling process that venous grafts undergo during exposure to arterial conditions is crucial to improve their patency, which is often compromised by intimal hyperplasia. The precise role of hemodynamic forces such as shear stress and arterial pressure in this remodeling is not fully characterized. The aim of this study was to determine the involvement of arterial shear stress and pressure on vein wall remodeling and to unravel the underlying molecular mechanisms. METHODS: An ex vivo vein support system was modified for chronic (up to 1 week), pulsatile perfusion of human saphenous veins under controlled conditions that permitted the separate control of arterial shear stress and different arterial pressure (7 mm Hg or 70 mm Hg). RESULTS: Veins perfused for 7 days under high pressure (70 mm Hg) underwent significant development of a neointima compared with veins exposed to low pressure (7 mm Hg). These structural changes were associated with altered expression of several molecular markers. Exposure to an arterial shear stress under low pressure increased the expression of matrix metalloproteinase (MMP)-2 and MMP-9 and tissue inhibitor of metalloproteinase (TIMP)-1 at the transcript, protein, and activity levels. This increase was enhanced by high pressure, which also increased TIMP-2 protein expression despite decreased levels of the cognate transcript. In contrast, the expression of plasminogen activator inhibitor-1 increased with shear stress but was not modified by pressure. Levels of the venous marker Eph-B4 were decreased under arterial shear stress, and levels of the arterial marker Ephrin-B2 were downregulated under high-pressure conditions. CONCLUSIONS: This model is a valuable tool to identify the role of hemodynamic forces and to decipher the molecular mechanisms leading to failure of human saphenous vein grafts. Under ex vivo conditions, arterial perfusion is sufficient to activate the remodeling of human veins, a change that is associated with the loss of specific vein markers. Elevation of pressure generates intimal hyperplasia, even though veins do not acquire arterial markers. CLINICAL RELEVANCE: The pathological remodeling of the venous wall, which leads to stenosis and ultimately graft failure, is the main limiting factor of human saphenous vein graft bypass. This remodeling is due to the hemodynamic adaptation of the vein to the arterial environment and cannot be prevented by conventional therapy. To develop a more targeted therapy, a better understanding of the molecular mechanisms involved in intimal hyperplasia is essential, which requires the development of ex vivo models of chronic perfusion of human veins.

Procedure for human saphenous veins ex vivo perfusion and external reinforcement.

The mainstay of contemporary therapies for extensive occlusive arterial disease is venous bypass graft. However, its durability is threatened by intimal hyperplasia (IH) that eventually leads to vessel occlusion and graft failure. Mechanical forces, particularly low shear stress and high wall tension, are thought to initiate and to sustain these cellular and molecular changes, but their exact contribution remains to be unraveled. To selectively evaluate the role of pressure and shear stress on the biology of IH, an ex vivo perfusion system (EVPS) was created to perfuse segments of human saphenous veins under arterial regimen (high shear stress and high pressure). Further technical innovations allowed the simultaneous perfusion of two segments from the same vein, one reinforced with an external mesh. Veins were harvested using a no-touch technique and immediately transferred to the laboratory for assembly in the EVPS. One segment of the freshly isolated vein was not perfused (control, day 0). The two others segments were perfused for up to 7 days, one being completely sheltered with a 4 mm (diameter) external mesh. The pressure, flow velocity, and pulse rate were continuously monitored and adjusted to mimic the hemodynamic conditions prevailing in the femoral artery. Upon completion of the perfusion, veins were dismounted and used for histological and molecular analysis. Under ex vivo conditions, high pressure perfusion (arterial, mean = 100 mm Hg) is sufficient to generate IH and remodeling of human veins. These alterations are reduced in the presence of an external polyester mesh.

Ex vivo Pulsatile Perfusion of Human Saphenous Veins Induces Intimal Hyperplasia and Increased Levels of the Plasminogen Activator Inhibitor 1.

Vessel wall trauma induces vascular remodeling processes including the development of intimal hyperplasia (IH). To assess the development of IH in human veins, we have used an ex vivo vein support system (EVVSS) allowing the perfusion of freshly isolated segments of saphenous veins in the presence of a pulsatile flow which reproduced arterial conditions regarding shear stress, flow rate and pressure during a period of 7 and 14 days. Compared to the corresponding freshly harvested human veins, histomorphometric analysis showed a significant increase in the intimal thickness which was already maximal after 7 days of perfusion. Expression of the endothelial marker CD31 demonstrated the presence of endothelium up to 14 days of perfusion. In our EVVSS model, the activity as well as the mRNA and protein expression levels of plasminogen activator inhibitor 1, the inhibitor of urokinase-type plasminogen activator (uPA) and tissue-type plasminogen activator (tPA), were increased after 7 days of perfusion, whereas the expression levels of tPA and uPA were not altered. No major change was observed between 7 and 14 days of perfusion. These data show that our newly developed EVVSS is a valuable setting to study ex vivo remodeling of human veins submitted to a pulsatile flow.

The use of external mesh reinforcement to reduce intimal hyperplasia and preserve the structure of human saphenous veins.

The saphenous vein is the conduit of choice in bypass graft procedures. Haemodynamic factors play a major role in the development of intimal hyperplasia (IH), and subsequent bypass failure. To evaluate the potential protective effect of external reinforcement on such a failure, we developed an ex vivo model for the perfusion of segments of human saphenous veins under arterial shear stress. In veins submitted to pulsatile high pressure (mean pressure at 100 mmHg) for 3 or 7 days, the use of an external macroporous polyester mesh 1) prevented the dilatation of the vessel, 2) decreased the development of IH, 3) reduced the apoptosis of smooth muscle cells, and the subsequent fibrosis of the media layer, 4) prevented the remodelling of extracellular matrix through the up-regulation of matrix metalloproteinases (MMP-2, MMP-9) and plasminogen activator type I. The data show that, in an experimental ex vivo setting, an external scaffold decreases IH and maintains the integrity of veins exposed to arterial pressure, via increase in shear stress and decrease wall tension, that likely contribute to trigger selective molecular and cellular changes.

Fluid dynamics simulation of right ventricular outflow tract oversizing.

OBJECTIVES: Repair of the right ventricular outflow tract (RVOT) in paediatric cardiac surgery remains challenging due to the high reoperation rate. Intimal hyperplasia and consequent arteriosclerosis is one of the most important limitation factors for graft durability. Since local shear stress and pressure are predictive elements for intimal hyperplasia and wall degeneration, we sought to determine in an oversized 12-mm RVOT model, with computed fluid dynamics simulation, the local haemodynamical factors that may explain intimal hyperplasia. This was done with the aim of identifying the optimal degree of oversizing for a 12-mm native RVOT. METHODS: Twenty domestic pigs, with a weight of 24.6 ± 0.89 kg and a native RVOT diameter of 12 ± 1.7 mm, had valve conduits of 12, 16, 18 and 20 mm implanted. Pressure and flow were measured at 75, 100 and 125% of normal flow at RVOT at the pulmonary artery, pulmonary artery bifurcation and at the left and right pulmonary arteries. Three-dimensional computed fluid dynamics (CFD) simulation in all four geometries in all flow modalities was performed. Local shear stress and pressure conditions were investigated. RESULTS: Corresponding to 75, 100 and 125% of steady-state flow, three inlet velocity profiles were obtained, 0.2, 0.29 and 0.36 m/s, respectively. At inflow velocity profiles, low shear stress areas, ranged from 0 to 2 Pa, combined with high-pressure areas ranging from 11.5 to 12.1 mmHg that were found at distal anastomosis, at bifurcation and at the ostia of the left and right pulmonary arteries in all geometries. CONCLUSIONS: In all three oversized geometries, the local reparation of shear stress and pressure in the 16-mm model showed a similar local profile as in the native 12 mm RVOT. According to these findings, we suggest oversizing the natural 12-mm RVOT by not more than 4 mm. The elements responsible for wall degeneration and intimal hyperplasia remain very similar to the conditions present in native RVOT.

Transient shear banding in time-dependent fluids

We study the dynamics of shear-band formation and evolution using a simple rheological model. The description couples the local structure and viscosity to the applied shear stress. We consider in detail the Couette geometry, where the model is solved iteratively with the Navier-Stokes equation to obtain the time evolution of the local velocity and viscosity fields. It is found that the underlying reason for dynamic effects is the nonhomogeneous shear distribution, which is amplified due to a positive feedback between the flow field and the viscosity response of the shear thinning fluid. This offers a simple explanation for the recent observations of transient shear banding in time-dependent fluids. Extensions to more complicated rheological systems are considered.

Understanding the mechanisms of lung mechanical stress

Physical forces affect both the function and phenotype of cells in the lung. Bronchial, alveolar, and other parenchymal cells, as well as fibroblasts and macrophages, are normally subjected to a variety of passive and active mechanical forces associated with lung inflation and vascular perfusion as a result of the dynamic nature of lung function. These forces include changes in stress (force per unit area) or strain (any forced change in length in relation to the initial length) and shear stress (the stress component parallel to a given surface). The responses of cells to mechanical forces are the result of the cell's ability to sense and transduce these stimuli into intracellular signaling pathways able to communicate the information to its interior. This review will focus on the modulation of intracellular pathways by lung mechanical forces and the intercellular signaling. A better understanding of the mechanisms by which lung cells transduce physical forces into biochemical and biological signals is of key importance for identifying targets for the treatment and prevention of physical force-related disorders.

Cyclic extrusion of a lava dome based on a stick-slip mechanism

Lava dome eruptions are sometimes characterised by large periodic fluctuations in extrusion rate over periods of hours that may be accompanied by Vulcanian explosions and pyroclastic flows. We consider a simple system of nonlinear equations describing a 1D flow of lava extrusion through a deep elastic dyke feeding a shallower cylindrical conduit in order to simulate this short-period cyclicity. Stick-slip conditions depending on a critical shear stress are assumed at the wall boundary of the cylindrical conduit. By analogy with the behaviour of industrial polymers in a plastic extruder, the elastic dyke acts like a barrel and the shallower cylindrical portion of the conduit as a die for the flow of magma acting as a polymer. When we applied the model to the Soufrière Hills Volcano, Montserrat, for which the key parameters have been evaluated from previous studies, cyclic extrusions with periods from 3 to 30 h were readily simulated, matching observations. The model also reproduces the reduced period of cycles observed when a major unloading event occurs due to lava dome collapse.

Polymer brushes under shear: molecular dynamics simulations compared to experiments

Surfaces coated with polymer brushes in a good solvent are known to exhibit excellent tribological properties. We have performed coarse-grained equilibrium and nonequilibrium molecular dynamics (MD) simulations to investigate dextran polymer brushes in an aqueous environment in molecular detail. In a first step, we determined simulation parameters and units by matching experimental results for a single dextran chain. Analyzing this model when applied to a multichain system, density profiles of end-tethered polymer brushes obtained from equilibrium MD simulations compare very well with expectations based on self-consistent field theory. Simulation results were further validated against and correlated with available experimental results. The simulated compression curves (normal force as a function of surface separation) compare successfully with results obtained with a surface forces apparatus. Shear stress (friction) obtained via nonequilibrium MD is contrasted with nanoscale friction studies employing colloidal-probe lateral force microscopy. We find good agreement in the hydrodynamic regime and explain the observed leveling-off of the friction forces in the boundary regime by means of an effective polymer–wall attraction.

Simulating startup shear of entangled polymer melts

Start-up shear rheology is a standard experiment used for characterizing polymer flow, and to test various models of polymer dynamics. A rich phenomenology is developed for behavior of entangled monodisperse linear polymers in such tests, documenting shear stress overshoots as a function of shear rates and molecular weights. A tube theory does a reasonable qualitative job at describing these phenomena, although it involves several drastic approximations and the agreement can be fortuitous. Recently, Lu and coworkers published several papers [e.g. Lu {\it et al.} {\it ACS Macro Lett}. 2014, 3, 569-573] reporting results from molecular dynamics simulations of linear entangled polymers, which contradict both theory and experiment. Based on these observations, they made very serious conclusions about the tube theory, which seem to be premature. In this letter, we repeat simulations of Lu {\it et al.} and systematically show that neither their simulation results, nor their comparison with theory are confirmed.

Effect of fiber orientation on the shear behavior of glass fiber/epoxy composites

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Static analysis of functionally graded cylindrical and conical shells or panels using the generalized unconstrained third order theory coupled with the stress recovery

A 2D Unconstrained Third Order Shear Deformation Theory (UTSDT) is presented for the evaluation of tangential and normal stresses in moderately thick functionally graded conical and cylindrical shells subjected to mechanical loadings. Several types of graded materials are investigated. The functionally graded material consists of ceramic and metallic constituents. A four parameter power law function is used. The UTSDT allows the presence of a finite transverse shear stress at the top and bottom surfaces of the graded shell. In addition, the initial curvature effect included in the formulation leads to the generalization of the present theory (GUTSDT). The Generalized Differential Quadrature (GDQ) method is used to discretize the derivatives in the governing equations, the external boundary conditions and the compatibility conditions. Transverse and normal stresses are also calculated by integrating the three dimensional equations of equilibrium in the thickness direction. In this way, the six components of the stress tensor at a point of the conical or cylindrical shell or panel can be given. The initial curvature effect and the role of the power law functions are shown for a wide range of functionally conical and cylindrical shells under various loading and boundary conditions. Finally, numerical examples of the available literature are worked out.

Analysis of coronary bifurcations by intravascular ultrasound and virtual histology

Background Regional differences in shear stress have been identified as reason for early plaque formation in vessel bifurcations. We aimed to investigate regional plaque morphology and composition using intravascular ultrasound (IVUS) and virtual histology (IVUS–VH) in coronary artery bifurcations. Methods We performed IVUS and IVUS–VH studies at coronary bifurcations to analyze segmental plaque burden and composition of different segments in relation to their orientation to the bifurcation. Results A total of 236 patients with a mean age of 59 ± 11 years (69% male) were analyzed. Plaque burden was higher at the contralateral vessel wall facing the bifurcation compared to the ipsilateral vessel wall and this difference was true for proximal and distal segments (proximal: 37 ± 12% and 45 ± 15% for segments at the ipsilateral and contralateral vessel wall, respectively, p < 0.001; distal: 37 ± 10% and 47 ± 15% for segments at the ipsilateral and contralateral vessel wall, respectively, p < 0.001). In addition, these segments exhibited a higher proportion of dense calcium and a lower proportion of fibrous tissue and fibro fatty tissue. Conclusions Segments on the contralateral wall of the bifurcation which have previously been identified as regions with low shear stress not only exhibited a higher plaque burden, but also a higher degree of calcification.