The impact of wall shear stress and pressure drop on the stability of the atherosclerotic plaque

Rupture of vulnerable atheromatous plaque in the carotid and coronary arteries often leads to stroke and heart attack respectively. The mechanism of blood flow and plaque rupture in stenotic arteries is still not fully understood. A three dimensional rigid wall model was solved under steady state conditions and unsteady conditions by assuming a time-varying inlet velocity profile to investigate the relative importance of axial forces and pressure drops in arteries with asymmetric stenosis. Flow-structure interactions were investigated for the same geometry and the results were compared with those retrieved with the corresponding 2D cross-section structural models. The Navier-Stokes equations were used as the governing equations for the fluid. The tube wall was assumed hyperelastic, homogeneous, isotropic and incompressible. The analysis showed that the three dimensional behavior of velocity, pressure and wall shear stress is in general very different from that predicted by cross-section models. Pressure drop across the stenosis was found to be much higher than shear stress. Therefore, pressure may be the more important mechanical trigger for plaque rupture other than shear stress, although shear stress is closely related to plaque formation and progression.

Simulation of the interaction between blood flow and atherosclerotic plaque

It has been well accepted that over 50% of cerebral ischemic events are the result of rupture of vulnerable carotid atheroma and subsequent thrombosis. Such strokes are potentially preventable by carotid interventions. Selection of patients for intervention is currently based on the severity of carotid luminal stenosis. It has been, however, widely accepted that luminal stenosis alone may not be an adequate predictor of risk. To evaluate the effects of degree of luminal stenosis and plaque morphology on plaque stability, we used a coupled nonlinear time-dependent model with flow-plaque interaction simulation to perform flow and stress/strain analysis for stenotic artery with a plaque. The Navier-Stokes equations in the Arbitrary Lagrangian-Eulerian (ALE) formulation were used as the governing equations for the fluid. The Ogden strain energy function was used for both the fibrous cap and the lipid pool. The plaque Principal stresses and flow conditions were calculated for every case when varying the fibrous cap thickness from 0.1 to 2mm and the degree of luminal stenosis from 10% to 90%. Severe stenosis led to high flow velocities and high shear stresses, but a low or even negative pressure at the throat of the stenosis. Higher degree of stenosis and thinner fibrous cap led to larger plaque stresses, and a 50% decrease of fibrous cap thickness resulted in a 200% increase of maximum stress. This model suggests that fibrous cap thickness is critically related to plaque vulnerability and that, even within presence of moderate stenosis, may play an important role in the future risk stratification of those patients when identified in vivo using high resolution MR imaging.

Flow of metal into the flash gap in the last stages of a plane-strain closed-die forging operation

In closed-die forging the flash geometry should be such as to ensure that the cavity is completely filled just as the two dies come into contact at the parting plane. If metal is caused to extrude through the flash gap as the dies approach the point of contact — a practice generally resorted to as a means of ensuring complete filling — dies are unnecessarily stressed in a high-stress regime (as the flash is quite thin and possibly cooled by then), which reduces the die life and unnecessarily increases the energy requirement of the operation. It is therefore necessary to carefully determine the dimensions of the flash land and flash thickness — the two parameters, apart from friction at the land, which control the lateral flow. The dimensions should be such that the flow into the longitudinal cavity is controlled throughout the operation, ensuring complete filling just as the dies touch at the parting plane. The design of the flash must be related to the shape and size of the forging cavity as the control of flow has to be exercised throughout the operation: it is possible to do this if the mechanics of how the lateral extrusion into the flash takes place is understood for specific cavity shapes and sizes. The work reported here is part of an ongoing programme investigating flow in closed-die forging. A simple closed shape (no longitudinal flow) which may correspond to the last stages of a real forging operation is analysed using the stress equilibrium approach. Metal from the cavity (flange) flows into the flash by shearing in the cavity in one of the three modes considered here: for a given cavity the mode with the least energy requirement is assumed to be the most realistic. On this basis a map has been developed which, given the depth and width of the cavity as well as the flash thickness, will tell the designer of the most likely mode (of the three modes considered) in which metal in the cavity will shear and then flow into the flash gap. The results of limited set of experiments, reported herein, validate this method of selecting the optimum model of flow into the flash gap.

Effect of a magnetic field on the flow and blood oxygenation in channels of variable cross-section

The effect of a magnetic field on the flow and oxygenation of an incompressible Newtonian conducting fluid in channels with irregular boundaries has been investigated. The geometric parameter δ, which is a ratio of the mean half width of the channel d to the characteristic length λ along the channel over which the significant changes in the flow quantities occur, has been used for perturbing the governing equations. Closed form solutions of the various order equations are presented for the stream function. The equations for oxygen partial pressure remain nonlinear even after perturbation, therefore a numerical solution is presented. The expressions for shear stress at a wall and pressure distributions are derived. Here the separation in the flow occurs at a higher Reynolds number than the corresponding non-magnetic case. It is found that the magnetic field has an effect on local oxygen concentration but has a little effect on the saturation length.

Classical and Cosserat plasticity and viscoplasticity models for slow granular flow

We consider models for the rheology of dense, slowly deforming granular materials based of classical and Cosserat plasticity, and their viscoplastic extensions that account for small but finite particle inertia. We determine the scale for the viscosity by expanding the stress in a dimensionless parameter that is a measure of the particle inertia. We write the constitutive relations for classical and Cosserat plasticity in stress-explicit form. The viscoplastic extensions are made by adding a rate-dependent viscous stress to the plasticity stress. We apply the models to plane Couette flow, and show that the classical plasticity and viscoplasticity models have features that depart from experimental observations; the prediction of the Cosserat viscoplasticity model is qualitatively similar to that of Cosserat plasticity, but the viscosities modulate the thickness of the shear layer.

Laminar flow of a suspension in a curved pipe with varying curvature

The laminar flow of a fairly concentrated suspension (in which the volume fraction Z of the solid particles < 0.4) in a spatially varying periodically curved pipe has been examined numerically. Unlike the case of interacting suspension flows, the particles are found to flow in a well-mixed fashion, altering both the axial and circumferential velocities and consequently the fluid flux in the tube, depending on their diffusivity and inertia. The magnitude of shear stress at the wall is enhanced, suggesting that, if applied to vascular system, the vascular wall could be prone to ulceration during pathological situations like polycythemia. The delay in adaptation of the deviation in Poiseuille flow velocity to the curvature changes is also discussed in detail.

Bed Shear from Velocity Profiles: A New Approach

A new binary law of velocity distribution has been developed to describe the velocity profile for the entire flow region. The law is a combination of logarithmic law, valid in the wall (inner) region, and parabolic law, valid in the core (outer) region of the flow. The validity of the law has been established based on earlier data on flat plates, rough and smooth pipes and experimental data obtained from rigid-walled open channels with plane sand beds. A procedure of estimating bed shear stress from the proposed law of velocity distribution using the measured velocity profile has been evolved. Bed shear estimates made according to this procedure are in agreement with the values obtained from uniform flow analysis in the case of open channel flow over a sediment bed. The proposed method of estimating the bed shear stress from the observed velocity profiles is found to be particularly useful in cases where it is difficult to determine precisely the true bed level, such as in the case of flow over sediment beds.

A simulation study of the mixing, atomic flow and velocity profiles of crystalline materials during sliding

Self-contained Non-Equilibrium Molecular Dynamics (NEMD) simulations using Lennard-Jones potentials were performed to identify the origin and mechanisms of atomic scale interfacial behavior between sliding metals. The mixing sequence and velocity profiles were compared via MD simulations for three cases, viz.: sell-mated, similar and hard-softvcrystal pairs. The results showed shear instability, atomic scale mixing, and generation of eddies at the sliding interface. Vorticity at the interface suggests that atomic flow during sliding is similar to fluid flow under Kelvin-Helmholtz instability and this is supported by velocity profiles from the simulations. The initial step-function velocity profile spreads during sliding. However the velocity profile does not change much at later stages of the simulation and it eventually stops spreading. The steady state friction coefficient during simulation was monitored as a function of sliding velocity. Frictional behavior can be explained on the basis of plastic deformation and adiabatic effects. The mixing layer growth kinetics was also investigated.

Magnetohydrodynamic unsteady flow in a tube of variable cross section in an axial magnetic field

Oscillatory flow in a tube of slowly varying cross section is investigated in the presence of a uniform magnetic field in the axial direction. A perturbation solution including steady streaming is presented. The pressure and shear stress on the wall for various parameters governing the flow are discussed. Physics of Fluids is copyrighted by The American Institute of Physics.

A hydrodynamical study of the flow in renal tubules

The hydrodynamical problem of flow in proximal renal tubule is investigated by considering axisymmetric flow of a viscous, incompressible fluid through a long narrow tube of varying cross-section with reabsorption at the wall. Two cases for reabsorption have been studied (i) when the bulk flow,Q, decays exponentially with the axial distancex, and (ii) whenQ is an arbitrary function ofx such thatQ-Q 0 can be expressed as a Fourier integral (whereQ 0 is the flux atx=0). The analytic expressions for flow variables have been obtained by applying perturbation method in terms of wall parameter ε. The effects of ε on pressure drop across the tube, radial velocity and wall shear have been studied in the case of exponentially decaying bulk flow and it has been found that the results are in agreement with the existing ones for the renal tubules.

Wall Shear Fluctuations in a Turbulent Boundary Layer

IN the last two decades, the instantaneous structure of a turbulent boundary layer has been examined by many in an effort to understand the dynamics of the flow. Distinct and well-defined flow patterns that seem to have great relevance to the turbulence production mechanism have been observed in the wall region.1'2 The flow near the wall is intermittent with periodic eruptions of the fluid, a phenomenon generally termed "bursting process." Earlier investigations in this field were limited to liquid flows at low speeds and the entire flowpattern was observed using flow visualization techniques.Study was later extended to boundary-layer flows in windtunnels at higher speeds and Reynolds numbers using hot-wiresignals for the analysis of the bursting phenomenon.

Pulsatile flow in tubes of varying cross-sections

The pulsatile flow of an incompressible viscous fluid in a cylindrical tube of varying cross section is investigated for small Reynolds numbers. The solutions consist of a stedy and an oscillatory part. The shear stress distribution on the wall is evaluated and discussed in detail for special geometries like tapered tubes, locally constricted tubes and peristaltic tubes. The existence of separation in the flow field is noticed.

Analysis of turbulent confined highly swirling flow

We report here a CFD model of highly swirling flow in a quarl burner using three versions of the k-epsilon model. Results for the recirculating zone, the bounding shear layer and the downstream flow are presented. We discuss, with suitable qualifications, how the model predictions can inform our understanding of this class of flows.

Transition of turbulent pipe flow

The velocity profile in turbulent pipe flow is usually divided into two regions, a wall or inner region and a core or outer region. For the inner region, the viscosity and wall shear stress are the important parameters governing the velocity distribution whereas for the outer region, the wall reduces the velocity below the maximum velocity independent of viscosity. In the present work, a velocity model is proposed for turbulent flow in the wall region of a pipe covering the entire transition from smooth to rough flows. Coupling this model for the wall region with the power law velocity model for the core region, an equation for the friction factor is obtained. The model constants are evaluated by using Nikuradse's experiments in the fully smooth and rough turbulent flows. The model shows good agreement with the friction factor and the velocity profiles obtained by Nikuradse for the transition region of turbulent flow.

Terahertz wave characteristics of a single-walled carbon nanotube containing a fluid flow using the nonlocal Timoshenko beam model

This paper presents the effect of nonlocal scaling parameter on the terahertz wave propagation in fluid filled single walled carbon nanotubes (SWCNTs). The SWCNT is modeled as a Timoshenko beam,including rotary inertia and transverse shear deformation by considering the nonlocal scale effects. A uniform fluid velocity of 1000 m/s is assumed. The analysis shows that, for a fluid filled SWCNT, the wavenumbers of flexural and shear waves will increase and the corresponding wave speeds will decrease as compared to an empty SWCNT. The nonlocal scale parameter introduces certain band gap region in both flexural and shear wave mode where no wave propagation occurs. This is manifested in the wavenumber plots as the region where the wavenumber tends to infinite (or wave speed tends to zero). The frequency at which this phenomenon occurs is called the ``escape frequency''. The effect of fluid density on the terahertz wave propagation in SWCNT is also studied and the analysis shows that as the fluid becomes denser, the wave speeds will decrease. The escape frequency decreases with increase in nonlocal scaling parameter, for both wave modes. We also show that the effect of fluid density and velocity are negligible on the escape frequencies of flexural and shear wave modes. (C) 2010 Elsevier B.V. All rights reserved.