Shmuel Einav

A laser-Doppler velocimetry study of ensemble-averaged characteristics of the turbulent near wake of a square cylinder

Ensemble-averaged statistics at constant phase of the turbulent near-wake flow (Reynolds number ≈ 21400 around a square cylinder have been obtained from two-component laser-Doppler measurements. Phase was defined with reference to a signal taken from a pressure sensor located at the midpoint of a cylinder sidewall. The distinction is drawn between the near wake where the shed vortices are ‘mature’ and distinct and a base region where the vortices grow to maturity and are then shed. Differences in length and velocity scales and vortex celerities between the flow around a square cylinder and the more frequently studied flow around a circular cylinder are discussed. Scaling arguments based on the circulation discharged into the near wake are proposed to explain the differences. The relationship between flow topology and turbulence is also considered with vorticity saddles and streamline saddles being distinguished. While general agreement with previous studies of flow around a circular cylinder is found with regard to essential flow features in the near wake, some previously overlooked details are highlighted, e.g. the possibility of high Reynolds shear stresses in regions of peak vorticity, or asymmetries near the streamline saddle. The base region is examined in more detail than in previous studies, and vorticity saddles, zero-vorticity points, and streamline saddles are observed to differ in importance at different stages of the shedding process.

A hypothesis for vulnerable plaque rupture due to stress-induced debonding around cellular microcalcifications in thin fibrous caps

In this article, we advance a hypothesis for the rupture of thin fibrous cap atheroma, namely that minute (10-μm-diameter) cellular-level microcalcifications in the cap, which heretofore have gone undetected because they lie below the visibility of current in vivo imaging techniques, cause local stress concentrations that lead to interfacial debonding. New theoretical solutions are presented for the local stress concentration around these minute spherical inclusions that predict a nearly 2-fold increase in interfacial stress that is relatively insensitive to the location of the hypothesized microinclusions in the cap. To experimentally confirm the existence of the hypothesized cellular-level microcalcifications, we examined autopsy specimens of coronary atheromatous lesions using in vitro imaging techniques whose resolution far exceeds conventional magnetic resonance imaging, intravascular ultrasound, and optical coherence tomography approaches. These high-resolution imaging modalities, which include confocal microscopy with calcium-specific staining and micro-computed tomography imaging, provide images of cellular-level calcifications within the cap proper. As anticipated, the minute inclusions in the cap are very rare compared with the numerous calcified macrophages observed in the necrotic core. Our mathematical model predicts that inclusions located in an area of high circumferential stress (>300 kPa) in the cap can intensify this stress to nearly 600 kPa when the cap thickness is <65 μm. The most likely candidates for the inclusions are either calcified macrophages or smooth muscle cells that have undergone apoptosis.

Influence of microcalcifications on vulnerable plaque mechanics using FSI modeling

Sudden heart attacks remain one of the primary causes of premature death in the developed world. Asymptomatic vulnerable plaques that rupture are believed to prompt such fatal heart attacks and strokes. The role of microcalcifications in the vulnerable plaque rupture mechanics is still debated. Recent studies suggest the microcalcifications increase the plaque vulnerability. In this manuscript we present a numerical study of the role of microcalcifications in plaque vulnerability in an eccentric stenosis model using a transient fluid-structure interaction (FSI) analysis. Two cases are being compared (i) in the absence of a microcalcification (ii) with a microcalcification spot fully embedded in the fibrous cap. Critical plaque stress/strain conditions were affected considerably by the presence of a calcified spot, and were dependent on the timing (phase) during the flow cycle. The vulnerable plaque with the embedded calcification spot presented higher wall stress concentration region in the fibrous cap a bit upstream to the calcified spot, with stress propagating to the deformable parts of the structure around the calcified spot. Following previous studies, this finding supports the hypothesis that microcalcifications increase the plaque vulnerability. Further studies in which the effect of additional microcalcifications and parametric studies of critical plaque cap thickness based on plaque properties and thickness, will help to establish the mechanism by which microcalcifications weaken the plaque and may lead to its rupture.

A decoupled coil detector array for fast image acquisition in magnetic resonance imaging

A method for magnetic resonance imaging (MRI) is investigated here, whereby an object is put under a homogeneous magnetic field, and the image is obtained by applying inverse source procedures to the data collected in an array of coil detectors surrounding the object. The induced current in each coil due to the precession of the magnetic dipole in each voxel depends on the characteristics of both the magnetic dipole frequency and strength, together with its distance from the coil, the coil direction in space, and the electrical properties of the coils. By calculating the induced current signals over an array of coil detectors, a relationship is established between the set of signals and the structure of the body under investigation. The linear relation can then be represented in matrix notation, and inversion of this matrix will produce an image of the body. Important problems which must be considered in the proposed method are signal-to-noise ratio (SNR) and coupling between adjacent coils. Solutions to these problems will provide a new method for obtaining an instantaneous image by NMR, with no need for gradient switching for encoding. A general algorithm for decoupling of the coils is presented and fast sampling of the signal, instead of filtering, is used in order to reduce both noise and numerical roundoff errors at the same time. Sensitivity considerations are made with respect to the number of coils that is required and its connection with coil radius and SNR. A computer simulation demonstrates the feasibility of this new modality. Based on the solutions presented here for the problems involved in the use of a large number of coils for a simultaneous recording of the signal, an improved method of multicoil recording is suggested, whereby it is combined with the conventional zeugmatographic method with read and phase gradients, to result in a novel method of magnetic resonance imaging. In the combined method, there are no phase-encoding gradients. Only a single slice-selecting gradient, to be followed by a single read-gradient. Instead of phase-encoding gradients, use is made of an equivalent number of coils. The number of coils now is reduced significantly. This method suggests a single slice image taken within a single echo time, and where a 128 x 128 resolution is possible with only 128 coils. The applicability of the method is based on a successful decoupling procedure for the detectors (coils and cables) and the availability of highly accurate, high-gain, low-noise amplifiers with a broad dynamic range.

Analysis of air flow patterns in the human nose

The nasal cavity is the main passage for air flow between the ambient atmosphere and the lungs. A preliminary requisite for any investigation of the mechanisms of each of its main physiological functions, such as filtration, air-conditioning and olfaction, is a basic knowledge of the air-flow pattern in this cavity. However, its complex three-dimensional structure and inaccessibility has traditionally prevented a detailed examination of internalin vivo orin vitro airflow patterns. To gain more insight into the flow pattern in inaccessible regions of the nasal cavity we have conducted a mathematical simulation of asymmetric airflow patterns through the nose. Development of a nose-like model, which resembles the complex structure of the nasal cavity, has allowed for a detailed analysis of various boundary conditions and structural parameters. The coronal and sagittal cross-sections of the cavity were modeled as trapezoids. The inferior and middle turbinates were represented by curved plates that emerge from the lateral walls. The airflow was considered to be incompressible, steady and laminar. Numerical computations show that the main air flux is along the cavity floor, while the turbinate structures direct the flow in an anterior-posterior direction. The presence of the turbinates and the trapezoidal shape of the cavity force more air flux towards the olfactory organs at the top of the cavity.

Abdominal aortic aneurysm risk of rupture: patient-specific FSI simulations using anisotropic model.

Abdominal aortic aneurysm (AAA) rupture represents a major cardiovascular risk, combining complex vascular mechanisms weakening the abdominal artery wall coupled with hemodynamic forces exerted on the arterial wall. At present, a reliable method to predict AAA rupture is not available. Recent studies have introduced fluid structure interaction (FSI) simulations using isotropic wall properties to map regions of stress concentrations developing in the aneurismal wall as a much better alternative to the current clinical criterion, which is based on the AAA diameter alone. A new anisotropic material model of AAA that closely matches observed biomechanical AAA material properties was applied to FSI simulations of patient-specific AAA geometries in order to develop a more reliable predictor for its risk of rupture. Each patient-specific geometry was studied with and without an intraluminal thrombus (ILT) using two material models-the more commonly used isotropic material model and an anisotropic material model-to delineate the ILT contribution and the dependence of stress distribution developing within the aneurismal wall on the material model employed. Our results clearly indicate larger stress values for the anisotropic material model and a broader range of stress values as compared to the isotropic material, indicating that the latter may underestimate the risk of rupture. While the locations of high and low stresses are consistent in both material models, the differences between the anisotropic and isotropic models become pronounced at large values of strain-a range that becomes critical when the AAA risk of rupture is imminent. As the anisotropic model more closely matches the biomechanical behavior of the AAA wall and resolves directional strength ambiguities, we conclude that it offers a more reliable predictor of AAA risk of rupture.

Dynamics of Blood Flow and Platelet Transport in Pathological Vessels

Abstract: Arterial disease, characterized by arterial occlusion (stenosis), is a leading cause of cardiovascular diseases and a major healthcare problem in the Western world. One of the main mechanisms leading to vessel occlusion is thrombus formation, which may be initiated by platelet activation. Shear rates and flow patterns (fluid dynamics factors) and concentration of coagulation factors and platelet agonists (biological factors) modulate platelet function and may lead to platelet activation and aggregation. Here, we examine the flow‐induced mechanisms leading to platelet activation in models of stenosed coronary vessels. Experimental and numerical methods were used to investigate and characterize the influence of the flow field on platelet activation. As it passes through pathological geometries characteristic of arterial stenosis, a platelet is exposed to varying levels of shear stress. The cumulative effect of the shear stress level and the duration of the platelets exposure to it determine whether the platelet is brought beyond its activation threshold. Stress histories of individual platelets can be tracked within the flow field to locate the regions where activated platelets might be found and subsequently aggregate and/or adhere to the wall.

Patient-Based Abdominal Aortic Aneurysm Rupture Risk Prediction with Fluid Structure Interaction Modeling

Elective repair of abdominal aortic aneurysm (AAA) is warranted when the risk of rupture exceeds that of surgery, and is mostly based on the AAA size as a crude rupture predictor. A methodology based on biomechanical considerations for a reliable patient-specific prediction of AAA risk of rupture is presented. Fluid–structure interaction (FSI) simulations conducted in models reconstructed from CT scans of patients who had contained ruptured AAA (rAAA) predicted the rupture location based on mapping of the stresses developing within the aneurysmal wall, additionally showing that a smaller rAAA presented a higher rupture risk. By providing refined means to estimate the risk of rupture, the methodology may have a major impact on diagnostics and treatment of AAA patients.

Evaluation of center-line extraction algorithms in quantitative coronary angiography

Objective testing of centerline extraction accuracy in quantitative coronary angiography (QCA) algorithms is a very difficult task. Standard tools for this task are not yet available. The authors present a simulation tool that generates synthetic angiographic images of a single coronary artery with predetermined centerline and diameter function. This simulation tool was used creating a library of images for the objective comparison and evaluation of QCA algorithms. This technique also provides the means for understanding the relationship between the algorithms performance and limitations and the vessels geometrical parameters. In this paper, two algorithms are evaluated and the results are presented.

Device Thrombogenicity Emulation: A Novel Method for Optimizing Mechanical Circulatory Support Device Thromboresistance

Mechanical circulatory support (MCS) devices provide both short and long term hemodynamic support for advanced heart failure patients. Unfortunately these devices remain plagued by thromboembolic complications associated with chronic platelet activation – mandating complex, lifelong anticoagulation therapy. To address the unmet need for enhancing the thromboresistance of these devices to extend their long term use, we developed a universal predictive methodology entitled Device Thrombogenicity Emulation (DTE) that facilitates optimizing the thrombogenic performance of any MCS device – ideally to a level that may obviate the need for mandatory anticoagulation. DTE combines in silico numerical simulations with in vitro measurements by correlating device hemodynamics with platelet activity coagulation markers – before and after iterative design modifications aimed at achieving optimized thrombogenic performance. DTE proof-of-concept is demonstrated by comparing two rotary Left Ventricular Assist Devices (LVADs) (DeBakey vs HeartAssist 5, Micromed Houston, TX), the latter a version of the former following optimization of geometrical features implicated in device thrombogenicity. Cumulative stresses that may drive platelets beyond their activation threshold were calculated along multiple flow trajectories and collapsed into probability density functions (PDFs) representing the device ‘thrombogenic footprint’, indicating significantly reduced thrombogenicity for the optimized design. Platelet activity measurements performed in the actual pump prototypes operating under clinical conditions in circulation flow loops – before and after the optimization with the DTE methodology, show an order of magnitude lower platelet activity rate for the optimized device. The robust capability of this predictive technology – demonstrated here for attaining safe and cost-effective pre-clinical MCS thrombo-optimization – indicates its potential for reducing device thrombogenicity to a level that may significantly limit the extent of concomitant antithrombotic pharmacotherapy needed for safe clinical device use.