Residence times and basins of attraction for a realistic right internal carotid artery with two aneurysms

We are presenting computational fluid dynamics simulation results for the flow in an anatomically accurate right internal carotid artery, exhibiting two saccular aneurysms close to each other. Our study focuses on the investigation of passage times for blood cells through the two-aneurysm malformation. We construct residence time maps that exhibit strong non-uniformity, linked to the entry of fluid in only the first, only the second, or in both aneurysms. An entrance index is computed, showing qualitatively the regions at an arterial section upstream of the aneurysms, where cells following one of these scenarios emanate. The significance of the residence time profiles and entry scenarios obtained is discussed with respect to thrombosis and pharmacokinetics. Preliminary evidence that the inflow-outflow patterns of the two aneurysms may be leading to particularly complex flow and to chaotic mixing is discussed.

Three-dimensional coupled fluid-structure simulation of pericardial bioprosthetic aortic valve function

A computational, three-dimensional coupled fluid-structure dynamics model was developed for a generic pericardial aortic valve in a rigid aortic root graft with physiologic sinuses. Valve geometry was based on that of the natural valve. Blood flow was modeled as pulsatile, laminar, Newtonian, incompressible flow. The structural model accounted for material and geometric nonlinearities and also simulated leaflet coaptation. A body fitted grid was used to subdivide the flow domain into computational finite volume cells. Shell finite elements were used to discretize the leaflet volume. A finite volume computational fluid dynamics code and finite element structure dynamics code were used to solve the flow and structure equations, respectively. The fluid flow and structural equations were coupled using an implicit "influence coefficient" technique. Physiologic ventricular and aortic pressure waveforms were prescribed as the flow boundary conditions. The aortic flow field, valve structural configuration, and leaflet stresses were computed at 2 msec intervals. Model predictions on aortic flow and transient variation in valve orifice area were in close agreement with corresponding experimental in vitro data. These findings suggest that the computer model has potential for being a powerful design tool for bioprosthetic aortic valves.

Doppler velocimetry in cavitating media

A numerical model has been developed to simulate propagation of ultrasonic beams in inhomogeneous moving media. The model is based on the ray theory of propagating waves, valid in the limit of high frequencies. The resulting equations depend only on local values of the velocity field and the speed of sound. In its implementation, the model assumes that the interactions of sound with the surrounding flow field are decoupled. This allows for applying the model in a post processing mode to flows computed by other means. The model was used to investigate beam behavior in unsteady cavitating flows. The study was motivated by reports of cavitation occurring in mitral bi-leaflet mechanical heart valves. The flow field and cavitation physics were simulated using a general purpose computer code, CFD-ACE. The ultrasonic beam model was then used to calculate the beam path, orientation, and frequency changes in the transient cavitating region. Results show that the presence of cavitation can fundamentally alter the beam propagation characteristics. Simple models that assume rectilinear propagation cannot, by definition, handle such flows. Cavitation incurs very large variations in the local sound speed, which in turn can induce very large distortions in the beam. This fact has strong ramifications regarding the accuracy of ultrasonic velocimetry when simple models are used to interpret Doppler data gathered under such flow conditions.

Application of computational fluid dynamics techniques to blood pumps

Present-day computational fluid dynamics (CFD) techniques can be used to analyze the behavior of fluid flow in a variety of pumps. CFD can be a powerful tool during the design stage for rapid virtual prototyping of different designs, analyzing performance parameters, and making design improvements. Computational flow solutions provide information such as the location and size of stagnation zones and the local shear rate. These parameters can be correlated to the extent of hemolysis and thrombus formation and are critical to the success of a blood pump. CFD-ACE, an advanced commercial CFD code developed by CFD Research Corporation, has been applied to fluid flows in rotary machines, such as axial flow pumps and inducers. Preprocessing and postprocessing tools for efficient grid generation and advanced graphical flow visualization are integrated seamlessly with CFD-ACE. The code has structured multiblock grid capability, non-Newtonian fluid treatment, a variety of turbulence models, and an Eulerian-Langrangian particle tracking model. CFD-ACE has been used successfully to study the flow characteristics in an axial flow blood pump. An unstructured flow solver that greatly automates the process of grid generation and speeds up the flow simulation is under development.

Numerical study of squeeze-flow in tilting disc mechanical heart valves

BACKGROUND AND AIM OF THE STUDY

The squeeze-flow that develops during valve closure is believed to cause cavitation in mitral mechanical heart valves (MHVs).

METHODS

Squeeze-flow was studied in tilting disc MHVs using two different numerical approaches. In the first decoupled analysis, experimental measurements of valve closing velocities were input into a computer model which simulated the resulting flow field. The second coupled approach involved simulation of the occluder motion and housing deformation in response to the surrounding flow.

RESULTS

Both models predicted the likelihood of cavitation during the squeeze-flow phase of valve closure. They also indicated that valve mounting compliance influences the squeeze-flow field. The coupled analysis also showed that squeeze-flow is influenced by fluid viscosity and the geometry of the contact region.

CONCLUSIONS

These parameters could therefore influence the inception of cavitation in tilting disc mitral MHVs.

An experimental-computational analysis of MHV cavitation: effects of leaflet squeezing and rebound

A combined experimental-computational study was performed to investigate the flow mechanics which could cause cavitation during the squeezing and rebounding phases of valve closure in the 29 mm mitral bileaflet Edwards-Duromedics (ED) mechanical heart valve (MHV). Leaflet closing motion was measured in vitro, and input into a computational fluid mechanics software package, CFD-ACE, to compute flow velocities and pressures in the small gap space between the occluder tip and valve housing. The possibility of cavitation inception was predicted when fluid pressures dropped below the saturated vapor pressure for blood plasma. The computational analysis indicated that cavitation is more likely to be induced during valve rebound rather than the squeezing phase of valve closure in the 29 mm ED-MHV. Also, there is a higher probability of cavitation at lower values of the gap width at the point of impact between the leaflet tip and housing. These predictions of cavitation inception are not likely to be significantly influenced by the water-hammer pressure gradient that develops during valve closure.

Dynamics of oxygen unloading from sickle erythrocytes

The objective of this work is to theoretically model oxygen unloading in sickle red cells. This has been done by combining into a single model diffusive transport mechanisms, which have been well-studied for normal red cells, and the hemoglobin polymerization process, which has been previously been studied for deoxyhemoglobin-S solutions and sickle cells in near-equilibrium situations. The resulting model equations allow us to study the important processes of oxygen delivery and polymerization simultaneously. The equations have been solved numerically by a finite-difference technique. The oxygen unloading curve for sickle erythrocytes is biphasic in nature. The rate of unloading depends in a complicated way on (a) the kinetics of hemoglobin S polymerization, (b) the kinetics of hemoglobin deoxygenation, and (c) the diffusive transport of both free oxygen and oxy-hemoglobin. These processes interact. For example, the hemoglobin S polymer interferes with the transport of both free oxygen and unpolymerized oxy-hemoglobin, and this is accounted for in the model by diffusivities which depend on the polymer and solution hemoglobin concentration. Other parameters which influence the interaction of these processes are the concentration of 2,3-diphosphoglycerate and total hemoglobin concentration. By comparing our model predictions for oxygen unloading with simpler predictions based on equilibrium oxygen affinities, we conclude that the relative rate of oxygen unloading of cells with different physical properties cannot be correctly predicted from the equilibrium affinities. To describe the unloading process, a kinetic calculation of the sort we give here is required.