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.

How critical is fibrous cap thickness to carotid plaque stability? A flow-plaque interaction model

Background and Purpose Acute cerebral ischemic events are associated with rupture of vulnerable carotid atheroma and subsequent thrombosis. Factors such as luminal stenosis and fibrous cap thickness have been thought to be important risk factors for plaque rupture. We used a flow-structure interaction model to simulate the interaction between blood flow and atheromatous plaque to evaluate the effect of the degree of luminal stenosis and fibrous cap thickness on plaque vulnerability. Methods A coupled nonlinear time-dependent model with a flow-plaque interaction simulation was used to perform flow and stress/strain analysis in a stenotic carotid artery model. The stress distribution within the plaque and the flow conditions within the vessel were calculated for every case when varying the fibrous cap thickness from 0.1 to 2 mm and the degree of luminal stenosis from 10% to 95%. A rupture stress of 300 kPa was chosen to indicate a high risk of plaque rupture. A 1-sample t test was used to compare plaque stresses with the rupture stress. Results High stress concentrations were found in the plaques in arteries with >70% degree of stenosis. Plaque stresses in arteries with 30% to 70% stenosis increased exponentially as fibrous cap thickness decreased. A decrease of fibrous cap thickness from 0.4 to 0.2 mm resulted in an increase of plaque stress from 141 to 409 kPa in a 40% degree stenotic artery. Conclusions There is an increase in plaque stress in arteries with a thin fibrous cap. The presence of a moderate carotid stenosis (30% to 70%) with a thin fibrous cap indicates a high risk for plaque rupture. Patients in the future may be risk stratified by measuring both fibrous cap thickness and luminal stenosis.

Unsteady Free-Convection Flow on a Rotating flat-Plate Under Nonuniform Gravity

The free convection problem with nonuniform gravity finds applications in several fields. For example, centrifugal gravity fieldsarisein many rotating machinery applications. A gravity field is also created artificially in an orbital space station by rotation. The effect of nonuniform gravity due to the rotation of isothermal or nonisothermal plates has been studied by several authors [l-5] using various mathematical techniques.

MHD Flow of a 2nd Grade Fluid in an Orthogonal Rheometer

The magnetohydrodynamics (MHD) flow of a conducting, homogeneous incompressible Rivlin-Ericksen fluid of second grade contained between two infinite, parallel, insulated disks rotating with the same angular velocity about two noncoincident axes, under the application of a uniform transverse magnetic field, is investigated. This model represents the MHD flow of the fluid in the instrument called an orthogonal rheometer, except for the fact that in the rheometer the rotating plates are necessarily finite. An exact solution of the governing equations of motion is presented. The force components in the x and y directions on the disks are calculated. The effects of magnetic field and the viscoelastic parameter on the forces are discussed in detail.

Unsteady incompressible boundary layer flow of a micropolar fluid at a stagnation point

The flow, heat and mass transfer on the unsteady laminar incompressible boundary layer in micropolar fluid at the stagnation point of a 2-dimensional and an axisymmetric body have been studied when the free stream velocity and the wall temperature vary arbitrarily with time. The partial defferential equations governing the flow have been solved numerically using a quasilinear finite-difference scheme. The skin friction, microrotation gradient and heat transfer parameters are found to be strongly dependent on the coupling parameter, mass transfer and time, whereas the effect of the microrotation parameter on the skin friction and heat transfer is rather weak, but microrotation gradient is strongly affected by it. The Prandtl number and the variation of the wall temperature with time affect the heat-transfer very significantly but the skin friction and micrortation gradient are unaffected by them.

Unsteady flow with attenuation in a fluid filled elastic tube with a stenosis

An equation governing the excess pressure has been derived, for an axially tethered and stenosed elastic tube filled with viscous liquid, by introducing the elasticity of the tube through pressure-area relation. This equation is solved numerically for large Womersley parameter and the results are presented for different types of pressure-radius relations and geometries by prescribing an outgoing wave suffering attenuation at some axial point of the tube. For a locally constricted tube it is observed that the pressure oscillates more and generates sound on the down stream side of the constriction.

Compressible boundary-layer flow at a three-dimensional stagnation point with massive blowing

The effect of massive blowing rates on the steady laminar compressible boundary-layer flow with variable gas properties at a 3-dim. stagnation point (which includes both nodal and saddle points of attachment) has been studied. The equations governing the flow have been solved numerically using an implicit finite-difference scheme in combination with the quasilinearization technique for nodal points of attachment but employing a parametric differentiation technique instead of quasilinearization for saddle points of attachment. It is found that the effect of massive blowing rates is to move the viscous layer away from the surface. The effect of the variation of the density- viscosity product across the boundary layer is found to be negligible for massive blowing rates but significant for moderate blowing rates. The velocity profiles in the transverse direction for saddle points of attachment in the presence of massive blowing show both the reverse flow as well as velocity overshoot.

Unsteady Attachment Line Flow On A Flat-Plate With Attached Cylinder

The unsteady laminar incompressible boundary-layer attachment-line flow on a flat plate with attached cylinder with heat and mass transfer has been studied when the free stream velocity, mass transfer and surface wall temperature vary arbitrarily with time. The governing partial differential equations with three independent variables have been solved numerically using an implicit finite-difference scheme. The heat transfer was found to be strongly dependent on the Prandtl number, variation of wall temperature with time and dissipation parameter (for large times). However, the free stream velocity distribution and mass transfer affect both the heat transfer and skin friction.

Flow through a plateau border of cellular foam

Experimentally measured average velocities through plateau borders of stationary cellular foam, when compared with those calculated with the assumption of rigid Plateau Border walls, show that the assumption of rigid walls severely underestimates the velocities. An analysis of the situation wherein plateau border walls have velocities, as decided by the surface viscosity of the system, is presented here. The plateau border is idealized as a pipe of equilateral triangular cross-section with vertices of the triangle having zero velocity. The pertinent form of Navier-Stoke's equations with inhomogeneous boundary conditions and its solution through a procedure of successive approximations is presented in dimensionless form. The solution reduces to the known solution of slow steady flow through a triangular pipe, when surface viscosity is infinite. Results indicate that the assumption of rigid plateau border walls is valid only when value of the inverse of dimensionless surface viscosity is less than 0.044. Beyond that the assumption severely underestimates the flow and the effect of nonrigidity of the wall must be considered.

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.

Microcontinuum approach to the pulsatile flow in tubes with and without longitudinal vibration

A fully developed pulsatile flow in a circular rigid tube is analysed by a microcontinuum approach. Solutions for radial variation of axial velocity and cell rotational velocity across the tube are obtained using the momentum integral method. Simplified forms of the solutions are presented for the relevant physiological data. Marked deviations in the results are observed when compared to a Newtonian fluid model. It is interesting to see that there is sufficient reduction in the mass flow rate, phase lag and friction due to the micropolar character of the fluid.

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.

Optimal hydrothermal load flow: Formulation and a successive approximation solution for fixed head systems

This paper deals with the optimal load flow problem in a fixed-head hydrothermal electric power system. Equality constraints on the volume of water available for active power generation at the hydro plants as well as inequality constraints on the reactive power generation at the voltage controlled buses are imposed. Conditions for optimal load flow are derived and a successive approximation algorithm for solving the optimal generation schedule is developed. Computer implementation of the algorithm is discussed, and the results obtained from the computer solution of test systems are presented.

Time dependent rotational flow of a viscous fluid over an infinite porous disk with a magnetic field

Both the semi-similar and self-similar flows due to a viscous fluid rotating with time dependent angular velocity over a porous disk of large radius at rest with or without a magnetic field are investigated. For the self-similar case the resulting equations for the suction and no mass transfer cases are solved numerically by quasilinearization method whereas for the semi-similar case and injection in the self-similar case an implicit finite difference method with Newton's linearization is employed. For rapid deceleration of fluid and for moderate suction in the case of self-similar flow there exists a layer of fluid, close to the disk surface where the sense of rotation is opposite to that of the fluid rotating far away. The velocity profiles in the absence of magnetic field are found to be oscillatory except for suction. For the accelerating freestream, (semi-similar flow) the effect of time is to reduce the amplitude of the oscillations of the velocity components. On the other hand the effect of time for the oscillating case is just the opposite.