A two-dimensional numerical analysis of unsteady flow in the carotid artery bifurcation. A comparison with three-dimensional in-vitro measurements and the influence of minor stenoses

In the present study a two-dimensional finite element model for incompressible Newtonian flow is applicated to the modelling of carotid artery flow. In earlier studies, the numerical model was validated experimentally for several flow configurations. In general the pulsatile flow is characterized by reversed flow regions at the non-divider side walls of both the internal and external carotid arteries. The unsteadiness of the flow is associated with rather complex spatial and temporal velocity distributions and leads to temporal variations of the location and length of the reversed flow regions. As a consequence, pronounced spatial and temporal variations in the wall shear stresses are found. At the non-divider side walls, wall shear stresses are relatively low and exhibits an oscillatory behaviour in space and time. At the divider side walls, wall shear stresses are relatively high and approximately follow the flow rate distribution in time. The aim of this study is not only to present two-dimensional calculations but also to compare the calculated two-dimensional velocity profiles with those from three-dimensional experiments. It is observed that in the common carotid artery and in the proximal parts of the internal and external carotid arteries, the two-dimensional numerical model provides valuable information with respect to the three-dimensional configuration. In the more distal parts of especially the internal carotid artery, deviations are found between the two-dimensional numerical and three-dimensional experimental model. These deviations can mainly be attributed to the neglect of the secondary velocity distribution in the two-dimensional model. In the two-dimensional numerical model the influence of a minor stenosis in the internal carotid artery is hardly distinguishable from a minor geometrical variation without stenosis. Full three-dimensional analyses of the influence of minor stenoses are needed to prove numerically whether in-vivo measurements of the axial velocity distribution are useful in the detection of minor stenoses.     

Flow pattern analysis in a highly stenotic patient-specific carotid bifurcation model using a turbulence model

The aim of this study is to investigate the blood flow pattern in carotid bifurcation with a high degree of luminal stenosis, combining in vivo magnetic resonance imaging (MRI) and computational fluid dynamics (CFD). A newly developed two-equation transitional model was employed to evaluate wall shear stress (WSS) distribution and pressure drop across the stenosis, which are closely related to plaque vulnerability. A patient with an 80% left carotid stenosis was imaged using high resolution MRI, from which a patient-specific geometry was reconstructed and flow boundary conditions were acquired for CFD simulation. A transitional model was implemented to investigate the flow velocity and WSS distribution in the patient-specific model. The peak time-averaged WSS value of approximately 73 Pa was predicted by the transitional flow model, and the regions of high WSS occurred at the throat of the stenosis. High oscillatory shear index values up to 0.50 were present in a helical flow pattern from the outer wall of the internal carotid artery immediately after the throat. This study shows the potential suitability of a transitional turbulent flow model in capturing the flow phenomena in severely stenosed carotid arteries using patient-specific MRI data and provides the basis for further investigation of the links between haemodynamic variables and plaque vulnerability. It may be useful in the future for risk assessment of patients with carotid disease.     

A model for shear stress sensing and transmission in vascular endothelial cells

Arterial endothelial cell (EC) responsiveness to flow is essential for normal vascular function and plays a role in the development of atherosclerosis. EC flow responses may involve sensing of the mechanical stimulus at the cell surface with subsequent transmission via cytoskeleton to intracellular transduction sites. We had previously modeled flow-induced deformation of EC-surface flow sensors represented as viscoelastic materials with standard linear solid behavior (Kelvin bodies). In the present article, we extend the analysis to arbitrary networks of viscoelastic structures connected in series and/or parallel. Application of the model to a system of two Kelvin bodies in parallel reveals that flow induces an instantaneous deformation followed by creeping to the asymptotic response. The force divides equally between the two bodies when they have identical viscoelastic properties. When one body is stiffer than the other, a larger fraction of the applied force is directed to the stiffer body. We have also probed the impact of steady and oscillatory flow on simple sensor-cytoskeleton-nucleus networks. The results demonstrated that, consistent with the experimentally observed temporal chronology of EC flow responses, the flow sensor attains its peak deformation faster than intracellular structures and the nucleus deforms more rapidly than cytoskeletal elements. The results have also revealed that a 1-Hz oscillatory flow induces significantly smaller deformations than steady flow. These results may provide insight into the mechanisms behind the experimental observations that a number of EC responses induced by steady flow are not induced by oscillatory flow.     

DSA-Based Quantitative Assessment of Cerebral Hypoperfusion in Patients with Asymmetric Carotid Stenosis

Digital subtraction angiography (DSA) is often used to evaluate the morphological and pathological changes of cerebral arteries in clinical practice. This study aims to explore the possibility of assessing cerebral hypoperfusion with DSA in patients with carotid stenosis. Thirty patients with a mild to severe stenosis on one side, and a mild stenosis on the other side of the carotid artery were recruited. Frontal, parietal, temporal and occipital lobes were chosen as regions of interest for measuring the quantitative perfusion parameters from their time-density curves (TDCs) of DSA images. The perfusion parameters were compared between the two hemispheres by using paired t-test. In addition, the bilateral asymmetry of these parameters was calculated and its correlation with the bilateral asymmetry in stenosis was analyzed. The parameters included mean transit time (MTT), time of contrast uptake (TU), time taken to the half peak value (1/2TMAX), area under the curve (AUC) were significantly prolonged at the severe stenosis side than those at the mild stenosis side in frontal lobe (P=0.013; P=0.041; P=0.009; P=0.027) and parietal lobe (P=0.008; P=0.041; P=0.002; P=0.012). The asymmetric ratios of MTT and AUC showed statistically significant correlations with stenosis asymmetry in all four lobes. MTT, TU, 1/2TMAX and AUC could reflect the bilateral asymmetry of the cerebral perfusion. These DSA parameters, therefore, may be used for the evaluation of cerebral hypoperfusion caused by carotid stenosis.     

Turbulence detection in a stenosed artery bifurcation by numerical simulation of pulsatile blood flow using the low-Reynolds number turbulence model

The pulsatile blood flow in a partially blocked artery is significantly altered as the flow regime changes through the cardiac cycle. This paper reports on the application of a low-Reynolds turbulence model for computation of physiological pulsatile flow in a healthy and stenosed carotid artery bifurcation. The human carotid artery was chosen since it has received much attention because atherosclerotic lesions are frequently observed. The Wilcox low-Re k-omega turbulence model was used for the simulation since it has proven to be more accurate in describing transition from laminar to turbulent flow. Using the FIDAP finite element code a validation showed very good agreement between experimental and numerical results for a steady laminar to turbulent flow transition as reported in a previous publication by the same authors. Since no experimental or numerical results were available in the literature for a pulsatile and turbulent flow regime, a comparison between laminar and low-Re turbulent calculations was made to further validate the turbulence model. The results of this study showed a very good agreement for velocity profiles and wall shear stress values for this imposed pulsatile laminar flow regime. To explore further the medical aspect, the calculations showed that even in a healthy or non-stenosed artery, small instabilities could be found at least for a portion of the pulse cycle and in different sections. The 40% and 55% diameter reduction stenoses did not significantly change the turbulence characteristics. Further results showed that the presence of 75% stenoses changed the flow properties from laminar to turbulent flow for a good portion of the cardiac pulse. A full 3D simulation with this low-Re-turbulence model, coupled with Doppler ultrasound, can play a significant role in assessing the degree of stenosis for cardiac patients with mild conditions.     

Weightlifting performance is related to kinematic and kinetic patterns of the hip and knee joints

The purpose of this study was to investigate the correlations between biomechanical outcome measures and weightlifting performance. Joint kinematics and kinetics of the hip, knee, and ankle were calculated while 10 subjects performed a clean at 85% of 1 repetition maximum (1RM). Kinematic and kinetic time-series patterns were extracted with principal components analysis. Discrete scores for each time-series pattern were calculated and used to determine how each pattern was related to body mass-normalized 1RM. Two hip kinematic and 2 knee kinetic patterns were significantly correlated with relative 1RM. The kinematic patterns captured hip and trunk motions during the first pull and hip joint motion during the movement transition between the first and second pulls. The first kinetic pattern captured a peak in the knee extension moment during the second pull. The second kinetic pattern captured a spatiotemporal shift in the timing and amplitude of the peak knee extension moment. The kinematic results suggest that greater lift mass was associated with steady trunk position during the first pull and less hip extension motion during the second-knee bend transition. Further, the kinetic results suggest that greater lift mass was associated with a smaller knee extensor moments during the first pull, but greater knee extension moments during the second pull, and an earlier temporal transition between knee flexion-extension moments at the beginning of the second pull. Collectively, these results highlight the importance of controlled trunk and hip motions during the first pull and rapid employment of the knee extensor muscles during the second pull in relation to weightlifting performance.     

Symptomatic carotid ischaemic events: safest and most cost effective way of selecting patients for angiography,before carotid endarterectomy.

OBJECTIVE: To determine the safest, least costly, and most effective way to select patients with symptomatic carotid ischaemic events for carotid angiography before carotid endarterectomy. DESIGN: Prospective cohort study. SETTING: University departments of clinical neurosciences and clinical neurology. PATIENTS: 485 Patients with carotid territory transient ischaemic attacks of the brain (n = 224) or eye (n = 162) or retinal infarction (n = 99) were referred to a single neurologist between 1976 and 1986. INTERVENTIONS: Clinical examination by auscultation over the precordium, supraclavicular fossae, and neck vessels (all patients). Cerebral angiography of patients suitable for carotid endarterectomy. MAIN OUTCOME MEASURES: Financial cost and number of disabling strokes after angiography. RESULTS: 296 Patients were investigated by cerebral angiography. Ischaemic symptoms had occurred in the distribution of 298 internal carotid arteries (symptomatic) that were imaged, two patients having bilateral symptoms. The presence or absence of a carotid bruit and the maximum percentage diameter stenosis of the origin of the symptomatic internal carotid artery were correlated. The prevalence of mild disease (diameter stenosis greater than or equal to 25%) of the symptomatic internal carotid artery was 57%, and if an ipsilateral carotid bruit was heard the probability of mild stenosis rose to 92%. The prevalence of moderate disease of the symptomatic internal carotid artery (stenosis greater than or equal to 50%) was 39%, and if a bruit was heard the probability doubled to 78%. The prevalence of severe internal carotid disease (stenosis greater than or equal to 75%) was 22%, and if a bruit was heard the probability was more than double, at 49%. The direct cost to both the NHS and the private health sector of investigating patients with symptomatic carotid ischaemia was estimated for several strategies of carotid artery imaging and expressed in terms of financial cost and number of strokes after angiography incurred in detecting all patients with diameter stenosis of the symptomatic internal carotid artery of greater than or equal to 25%, 50%, or 75%. To detect diameter stenosis of the internal carotid artery of greater than or equal to 25% it is most cost effective to proceed directly to cerebral angiography in patients with a carotid bruit over the symptomatic carotid bifurcation and to screen patients without a carotid bruit by duplex carotid ultrasonography; patients in whom duplex ultrasonography discloses stenosis of greater than or equal to 25% are then referred for cerebral angiography. To detect only more severe internal carotid disease (stenosis of greater than or equal to 50%) the same policy applies, unless the local duplex ultrasonographic service is particularly efficient and reliable, when it is probably most cost effective and safer to screen all patients by this method irrespective of the findings on cervical auscultation. To detect stenosis of 75% or greater it is most cost effective to screen all patients with duplex ultrasonography, whether a carotid bruit is present or not, because this approach reduces the number of angiograms required, is the least expensive, and results in the least number of strokes after angiography. CONCLUSIONS: Patients selection for cerebral angiography before carotid endarterectomy needs to be appropriate and cost effective. Sound clinical evaluation and duplex carotid ultrasound are required. The findings of this study should not be applied to other medical centres without first considering possible differences in the prevalence of carotid artery disease, the efficiency and reliability of duplex ultrasonography, the local complication rates of cerebral angiography, and the local costs of the imaging procedures.     

Dependence of central airway resistance on frequency and tidal volume: a model study

The resistance of a hollow cast of human central airways was measured during true sinusoidal airflow oscillations over a wide range of frequencies (0.5-40 Hz) and for various flow amplitudes up to 8 l/s. Pressure and flow were measured in the trachea with high-performance transducers, digitized and averaged over 100 cycles. Data were studied at two points in the flow cycle: at peak inspiratory and expiratory flows and in the two neighborhoods around zero flow where airway resistance (Rv approximately equal to o) was taken as the average slope of the pressure-flow (P-V) curve in each zone. When data obtained near peak flow were plotted in terms of dimensionless pressure drop vs. peak Reynolds number (Rem) and compared with steady-state data, we found no difference up to 2 Hz as previously reported (Isabey and Chang, J. Appl. Physiol. 51: 1338-1348, 1981), a slight decay in pressure drop between 4 and 8 Hz, a frequency-dependent increase in peak flow resistance at high frequencies (10-40 Hz) governed by the Strouhal number alpha 2/Rem beyond alpha 2/Rem = 0.5. On the other hand RV approximately equal to o was found to increase relative to steady state as local acceleration increases, e.g., as peak flow increases at a fixed frequency; this differs from the classical linear theory of oscillatory flow in a long straight tube. To explain these results, we had to use, as in our previous study, an alternative expression for the Strouhal number, i.e., epsilon = L X A X (dV/dt)/V2 (where L and A are the length and cross-sectional area of the trachea and V is a constant flow range over which resistance around flow reversal was computed), which accurately reflects the ratio of local acceleration [d(V/A)/dt)] to convective acceleration [(V/A)2/L] in developing branching flow. Finally, to delineate the regions of dominance of each of the dimensionless parameters, we compiled frequency-tidal volume diagrams for peak flows as well as for reversal. Epsilon, which is negligible near peak flows, appeared to govern the oscillatory P-V relationship near flow reversal in a transitional region of the diagram located between regions of steadiness, or moderate unsteadiness, and a region of dominant unsteadiness governed by alpha.     

Numerical investigation of physiologically realistic pulsatile flow through arterial stenosis

Numerical simulations of pulsatile blood flow in straight tube stenosis models were performed to investigate the poststenotic flow phenomena. In this study, three axisymmetrical and three asymmetrical stenosis models with area reduction of 25%, 50% and 75% were constructed. A measured human common carotid artery blood flow waveform was used as the upstream flow condition which has a mean Reynold's number of 300. All calculations were performed with high spatial and temporal resolutions. Flow features such as velocity profiles, flow separation zone (FSZ), and wall shear stress (WSS) distributions in the poststenotic region for all models are presented. The results have demonstrated that the formation and development of FSZs in the poststenotic region are very complex, especially in the flow deceleration phase. In axisymmetric stenoses the poststenotic flow is more sensitive to changes in the degree of stenosis than in asymmetric models. For severe stenoses, the stenosis influence length is shorter in asymmetrical models than in axisymmetrical cases. WSS oscillations (between positive and negative values) have been observed at various downstream locations in some models. The amplitude of the oscillation depends strongly on the axial location and the degree of stenosis.     

Flow field and oscillatory shear stress in a tuning-fork-shaped model of the average human carotid bifurcation

The oscillatory shear index (OSI) was developed based on the hypothesis that intimal hyperplasia was correlated with oscillatory shear stresses. However, the validity of the OSI was in question since the correlation between intimal thickness and the OSI at the side walls of the sinus in the Y-shaped model of the average human carotid bifurcation (Y-AHCB) was weak. The objectives of this paper are to examine whether the reason for the weak correlation lies in the deviation in geometry of Y-AHCB from real human carotid bifurcation, and whether this correlation is clearly improved in the tuning-fork-shaped model of the average human carotid bifurcation (TF-AHCB). The geometry of the TF-AHCB model was based on observation and statistical analysis of specimens from 74 cadavers. The flow fields in both models were studied and compared by using flow visualization methods under steady flow conditions and by using laser Doppler anemometer (LDA) under pulsatile flow conditions. The TF-shaped geometry leads to a more complex flow field than the Y-shaped geometry. This added complexity includes strengthened helical movements in the sinus, new flow separation zone, and directional changes in the secondary flow patterns. The results show that the OSI-values at the side walls of the sinus in the TF-shaped model were more than two times as large as those in the Y-shaped model. This study confirmed the stronger correlation between the OSI and intimal thickness in the tuning-fork geometry of human carotid bifurcation, and the TF-AHCB model is a significant improvement over the traditional Y-shaped model.     

Turbulence modeling in three-dimensional stenosed arterial bifurcations

Under normal healthy conditions, blood flow in the carotid artery bifurcation is laminar. However, in the presence of a stenosis, the flow can become turbulent at the higher Reynolds numbers during systole. There is growing consensus that the transitional k-omega model is the best suited Reynolds averaged turbulence model for such flows. Further confirmation of this opinion is presented here by a comparison with the RNG k-epsilon model for the flow through a straight, nonbifurcating tube. Unlike similar validation studies elsewhere, no assumptions are made about the inlet profile since the full length of the experimental tube is simulated. Additionally, variations in the inflow turbulence quantities are shown to have no noticeable affect on downstream turbulence intensity, turbulent viscosity, or velocity in the k-epsilon model, whereas the velocity profiles in the transitional k-omega model show some differences due to large variations in the downstream turbulence quantities. Following this validation study, the transitional k-omega model is applied in a three-dimensional parametrically defined computer model of the carotid artery bifurcation in which the sinus bulb is manipulated to produce mild, moderate, and severe stenosis. The parametric geometry definition facilitates a powerful means for investigating the effect of local shape variation while keeping the global shape fixed. While turbulence levels are generally low in all cases considered, the mild stenosis model produces higher levels of turbulent viscosity and this is linked to relatively high values of turbulent kinetic energy and low values of the specific dissipation rate. The severe stenosis model displays stronger recirculation in the flow field with higher values of vorticity, helicity, and negative wall shear stress. The mild and moderate stenosis configurations produce similar lower levels of vorticity and helicity.     

Numerical analysis of flow through a severely stenotic carotid artery bifurcation

The results of computational simulations may supplement MR and other in vivo diagnostic techniques to provide an accurate picture of the hemodynamics in particular vessels, which may help demonstrate the risks of embolism or plaque rupture posed by particular plaque deposits. In this study, a model based on an endarterectomy specimen of the plaque in a carotid bifurcation was examined. The flow conditions include steady flow at Reynolds numbers of 300, 600, and 900 as well as unsteady pulsatile flow. Both dynamic pressure and wall shear stress are very high, with shear values up to 70 N/m2, proximal to the stenosis throat in the internal carotid artery, and both vary significantly through the flow cycle. The wall shear stress gradient is also strong along the throat. Vortex shedding is observed downstream of the most severe occlusion. Two turbulence models, the Chien and Goldberg varieties of k-epsilon, are tested and evaluated for their relevance in this geometry. The Chien model better captures phenomena such as vortex shedding. The flow distal to stenosis is likely transitional, so a model that captures both laminar and turbulent behavior is needed.     

In vivo demonstration of flow recirculation and turbulence downstream of graded stenoses in canine arteries

The velocity field around arterial stenoses was investigated using a pulsed doppler velocimeter in vivo. Asymmetric zones of recirculation were identified by systolic flow reversal in the post-stenotic field in carotid and iliac arteries of anesthetised dogs. There was a close correlation between shear intensity and turbulence as estimated by the velocity difference between the jet and the recirculation zone and by maximum spectral width respectively. Under the conditions of these experiments, stenosis grade (% diameter reduction) dominated hemodynamic variables such as Reynolds number, oscillation and pulsatility in determining the intensity of turbulence. The method used does not appear to have sufficient resolution to distinguish between turbulence and discrete oscillating velocities (vorticity) nor to allow determination of wall shear stress though the pattern of change of the latter is similar to that found downstream of axisymmetric stenosis in models using steady flow.     

A model for shear stress-induced deformation of a flow sensor on the surface of vascular endothelial cells

Fluid mechanical shear stress elicits humoral, metabolic, and structural responses in vascular endothelial cells (ECs); however, the mechanisms involved in shear stress sensing and transduction remain incompletely understood. Beyond being responsive to shear stress, ECs distinguish among and respond differently to different types of shear stress. Recent observations suggest that endothelial shear stress sensing may occur through direct interaction of the flow with cell-surface structures that act as primary flow sensors. This paper presents a mathematical model for the shear stress-induced deformation of a flow sensor on the EC surface. The sensor is modeled as a cytoskeleton-coupled viscoelastic structure exhibiting standard linear solid behavior. Since ECs respond differently to different types of flow, the deformation and resulting velocity of the sensor in response to steady, non-reversing pulsatile, and oscillatory flow have been studied. Furthermore, the sensitivity of the results to changes in various model parameters including the magnitude of applied shear stress, the constants that characterize the viscoelastic behavior, and the pulsatile flow frequency (f) has been investigated. The results have demonstrated that in response to a suddenly applied shear stress, the sensor exhibits a level of instantaneous deformation followed by gradual creeping to the long-term response. The peak deformation increases linearly with the magnitude of the applied shear stress and decreases for viscoelastic constants that correspond to stiffer sensors. While the sensor deformation depends on f for low f values, the deformation becomes f -independent above a critical threshold frequency. Finally, the peak sensor deformation is considerably larger for steady and non-reversing pulsatile flow than for oscillatory flow. If the extent of sensor deformation correlates with the intensity of flow-mediated endothelial signaling, then our results suggest possible mechanisms by which ECs distinguish among steady, non-reversing pulsatile, and oscillatory shear stress.     

A theory of blood flow in skeletal muscle

A theoretical analysis of blood flow in the microcirculation of skeletal muscle is provided. The flow in the microvessels of this organ is quasi steady and has a very low Reynolds number. The blood is non-Newtonian and the blood vessels are distensible with viscoelastic properties. A formulation of the problem is provided using a viscoelastic model for the vessel wall which was recently derived from measurements in the rat spinotrapezius muscle (Skalak and Schmid-Sch?nbein, 1986b). Closed form solutions are derived for several physiologically important cases, such as perfusion at steady state, transient and oscillatory flows. The results show that resting skeletal muscle has, over a wide range of perfusion pressures an almost linear pressure-flow curve. At low flow it exhibits nonlinearities. Vessel distensibility and the non-Newtonian properties of blood both have a strong influence on the shape of the pressure-flow curve. During oscillatory flow the muscle exhibits hysteresis. The theoretical results are in qualitative agreement with experimental observations.     

Flow in a catheterized curved artery with stenosis

The fluid mechanics of blood flow in a catheterized curved artery with stenosis is studied through a mathematical analysis. Blood is modelled as an incompressible Newtonian fluid and the flow is assumed to be steady and laminar. An approximate analytic solution to the problem is obtained through a double series perturbation analysis for the case of small curvature and mild stenosis. The effect of catheterization on various physiologically important flow characteristics (i.e. the pressure drop, impedance and the wall shear stress) is studied for different values of the catheter size and Reynolds number of the flow. It is found that all these flow characteristics vary markedly across a stenotic lesion. Also, increase in the catheter size leads to a considerable increase in their magnitudes. These results are used to obtain the estimates of increased pressure drop across an arterial stenosis when a catheter is inserted into it. Our calculations, based on the geometry and flow conditions existing in coronary arteries, suggest that, in the presence of curvature and stenosis, and depending on the value of k (ratio of catheter size to vessel size) ranging from 0.1 to 0.4, the pressure drop increases by a factor ranging from 1.60 to 5.16. But, in the absence of curvature and stenosis, with the same range of catheter size, this increased factor is about 1.74-4.89. These estimates for the increased pressure drop can be used to correct the error involved in the measured pressure gradients using catheters. The combined effects of stenosis and curvature on flow characteristics are also studied in detail. It is found that the effect of stenosis is more dominant than that of the curvature. Due to the combined effect of stenosis, curvature and catheterization, the secondary streamlines are modified in a cross-sectional plane. The insertion of a catheter into the artery leads to the formation of increased number of secondary vortices.     

Velocity profiles and streamlines of a revolution post-stenotic flow.

Velocity fields have been measured in models simulating arterial stenoses for continuous and revolution flows. A pulsed Doppler velocimeter allows for velocity readings in the entire tube and in the wall area. Streamlines are determined by numerical solving of the system of equations defining the current function. Velocity profiles and streamlines are presented and discussed either for steady or for unsteady flows, with different Reynolds numbers and variable degrees of stenosis. There is, in the wall area, a recirculating zone made of a well-defined rouleau. Its length varies increasingly according to the increasing severity of the stenosis. The stability of axial flow depends on the input profile, the degree of stenosis and the Reynolds number. Plotting streamlines allows to describe accurately the flow; its quantitative aspect offers advantages with respect to conventional visualization mode.     

Flow patterns and wall shear stress distribution in human internal carotid arteries: the geometric effect on the risk for stenoses

It has been widely observed that atherosclerotic stenosis occurs at sites with complex hemodynamics, such as arteries with high curvature or bifurcations. These regions usually have very low or highly oscillatory wall shear stress (WSS). In the present study, 3D sinusoidally pulsatile blood flow through the models of internal carotid artery (ICA) with different geometries was investigated with computational simulation. Three preferred sites of stenoses were found along the carotid siphon with low and highly oscillatory WSS. The risk for stenoses at these sites was scaled with the values of time-averaged WSS and oscillating shear index (OSI). The local risk for stenoses at every preferred site of stenoses was found different between 3 types of ICA, indicating that the geometry of the blood vessel plays significant roles in the atherogenesis. Specifically, the large curvature and planarity of the vessel were found to increase the risk for stenoses, because they tend to lower WSS and elevate OSI. Therefore, the geometric study makes it possible to estimate the stenosis location in the ICA siphon as long as the shape of ICA was measured.     

颈动脉支架置入术后患者认知功能的变化及CT灌注成像技术的疗效评估价值研究

目的:探讨颈动脉狭窄患者在行颈动脉支架置入术后认知功能的变化情况,并分析CT灌注成像(CTP)对手术疗效的评估价值。方法:选取2015年10月到2018年3月在济宁医学院附属日照市人民医院接受治疗的颈动脉狭窄患者80例,其中有症状性颈动脉狭窄患者49例作为有症状组,无症状性颈动脉狭窄患者31例作为无症状组。所有患者均接受颈动脉支架置入手术及CTP检查,采用蒙特利尔认知评估量表(MoCA)、简易智能精神状况量表(MMSE)、搭火柴测验(Stick Test)综合评价颈动脉狭窄患者的术前、术后1周、术后3个月、术后6个月认知功能的变化情况,比较有症状组和无症状组患者的CTP相对灌注参数。结果:术后1周,颈动脉狭窄患者的MoCA总分、视空间/执行能力、注意力、延迟回忆以及MMSE总分、Stick Test总分较术前有所降低(P0.05);术后3个月、术后6个月,颈动脉狭窄患者的MoCA总分、视空间/执行能力、注意力、延迟回忆以及MMSE总分、Stick Test总分较术前有所升高(P0.05)。术前,有症状组的相对血流达峰时间、相对平均通过时间长于无症状组,相对脑血流量低于无症状组(P0.05);术后1周,无症状组的相对血流达峰时间较术前有所缩短,且短于有症状组(P0.05);术后1周,有症状组的相对血流达峰时间、相对平均通过时间较术前有所缩短,相对脑血流量较术前有所升高(P0.05);两组术前、术后1周相对脑血容量比较均无统计学差异(P0.05)。结论:颈动脉支架置入术后患者会出现暂时的、可逆的认知功能恶化,但最终认知功能会得到明显的改善。CTP可发现异常的脑灌注情况,同时能够较好地评价颈动脉支架置入术治疗颈动脉狭窄患者的疗效。     

Secondary velocity fields in the conducting airways of the human lung

An understanding of flow and dispersion in the human respiratory airways is necessary to assess the toxicological impact of inhaled particulate matter as well as to optimize the design of inhalable pharmaceutical aerosols and their delivery systems. Secondary flows affect dispersion in the lung by mixing solute in the lumen cross section. The goal of this study is to measure and interpret these secondary velocity fields using in vitro lung models. Particle image velocimetry experiments were conducted in a three-generational, anatomically accurate model of the conducting region of the lung. Inspiration and expiration flows were examined under steady and oscillatory flow conditions. Results illustrate secondary flow fields as a function of flow direction, Reynolds number, axial location with respect to the bifurcation junction, generation, branch, phase in the oscillatory cycle, and Womersley number. Critical Dean number for the formation of secondary vortices in the airways, as well as the strength and development length of secondary flow, is characterized. The normalized secondary velocity magnitude was similar on inspiration and expiration and did not vary appreciably with generation or branch. Oscillatory flow fields were not significantly different from corresponding steady flow fields up to a Womersley number of 1 and no instabilities related to shear were detected on flow reversal. These observations were qualitatively interpreted with respect to the simple streaming, augmented dispersion, and steady streaming convective dispersion mechanisms.