Fluid dynamic analysis of the 50 cc Penn State artificial heart under physiological operating conditions using particle image velocimetry

In order to bridge the gap of existing artificial heart technology to the diverse needs of the patient population, we have been investigating the viability of a scaled-down design of the current 70 cc Penn State artificial heart. The issues of clot formation and hemolysis may become magnified within a 50 cc chamber compared to the existing 70 cc one. Particle image velocimetry (PIV) was employed to map the entire 50 cc Penn State artificial heart chamber. Flow fields constructed from PIV data indicate a rotational flow pattern that provides washout during diastole. In addition, shear rate maps were constructed for the inner walls of the heart chamber The lateral walls of the mitral and aortic ports experience high shear rates while the upper and bottom walls undergo low shear rates, with sufficiently long exposure times to potentially induce platelet activation or thrombus formation. In this study, we have demonstrated that PIV may adequately map the flow fields accurately in a reasonable amount of time. Therefore, the potential exists of employing PIV as a design tool.     

Flow field analysis in a compliant acinus replica model using particle image velocimetry (PIV)

Inhaled particles reaching the alveolar walls have the potential to cross the blood–gas barrier and enter the blood stream. Experimental evidence of pulmonary dosimetry, however, cannot be explained by current whole lung dosimetry models. Numerical and experimental studies shed some light on the mechanisms of particle transport, but realistic geometries have not been investigated. In this study, a three dimensional expanding model including two generations of respiratory bronchioles and five terminal alveolar sacs was created from a replica human lung cast. Flow visualization techniques were employed to quantify the fluid flow while utilizing streamlines to evaluate recirculation. Pathlines were plotted to track the fluid motion and estimate penetration depth of inhaled air. This study provides evidence that the two generations immediately proximal to the terminal alveolar sacs do not have recirculating eddies, even for intense breathing. Results of Peclet number calculations indicate that substantial convective motion is present in vivo for the case of deep breathing, which significantly increases particle penetration into the alveoli. However, particle diffusion remains the dominant mechanism of particle transport over convection, even for intense breathing because inhaled particles do not reach the alveolar wall in a single breath by convection alone. Examination of the velocity fields revealed significant uneven ventilation of the alveoli during a single breath, likely due to variations in size and location. This flow field data, obtained from replica model geometry with realistic breathing conditions, provides information to better understand fluid and particle behavior in the acinus region of the lung.     

Visualization of the intracavitary blood flow in systemic ventricles of Fontan patients by contrast echocardiography using particle image velocimetry

Objectives

Two-dimensional strain echocardiography (2DSE) technique has enabled accurate quantification of regional myocardial function. This experimental study was aimed to investigate the value of 2DSE in detection of segmental regional myocardial dysfunction induced by fibrosis following myocardial infarction in a small animal (rat) model.

Methods

A rat model of myocardial infarction was established by ligation of the proximal left anterior descending coronary artery in 17 SD rats. Regional myocardial function was detected by 2DSE at baseline and 4-weeks post-infarction, including end-systolic radial strain and strain rate (SR and SrR) and end-systolic circumferential strain and strain rate (SC and SrC) of each of six segments at papillary level. According to the size of scar found by histologic Masson staining, the optimal cutoff points of parameters for detecting scar area were analyzed and the sensitivity and specificity of every parameter to detect myocardial scar were obtained using ROC.

Results

(1) Comparing with parameters measured at baseline, there were significant decreases in SR, SrR, SC and SrC of each segment at 4 weeks post-infarction, with the worst in the infarct area (32.90 ± 8.79 vs 11.18 ± 3.89, 6.28 ± 1.35 vs 3.18 ± 0.47, -14.46 ± 2.21 vs -6.30 ± 2.17 and 4.93 ± 0.95 vs 2.59 ± 1.16, respectively) (all P < 0.05). (2)By 4 weeks, the myocardium of infarct area (anteroseptum, anterior and anterolateral) had fibrosis (31.33 ± 9.89, 73.42 ± 13.21 and 13.99 ± 3.24%, respectively) with minimal fibrosis in inferoseptal segment (0.32 ± 0.19%), no fibrosis was found in the inferior and inferolateral segments. (3)Significant negative correlations were found between the size of segmental scar and 2DSE parameters (r-value -0.61 ~ -0.80, all P < 0.01) with the strongest correlation in SR. SR less than 10% has 84% sensitivity and 98% specificity for detecting segments of scar area greater than 30% with AUC = 0.97.

Conclusions

2DSE is able to assess regional myocardial dysfunction in a rat model of myocardial infarction and has high accuracy in detecting infarct segments with scar area greater than 30%.     

Three-dimensional heart motion estimation using endoscopic monocular vision system: From artificial landmarks to texture analysis

In robot-assisted beating heart surgery, motion of the heart surface might be virtually stabilized to let the surgeon work as in on-pump cardiac surgery. Virtual stabilization means to compensate physically the relative motion between the instrument tool tip and the region of interest on the heart surface, and to offer surgeon a stable visual display of the scene. To this end, motion of the heart must be estimated. This article focusses on motion estimation of the heart surface. Two approaches are considered in the paper. The first one is based on landmark tracking allowing 3D pose estimation. The second is based on texture tracking. Classical computer vision methods, as well as a new texture-based tracking scheme has been applied to track the heart motion, and, when possible, reconstruct 3D distance to the heart surface. Experimental results obtained on in vivo images show the estimated motion of heart surface points.     

Turbulence intensity measurements using particle image velocimetry in diseased carotid artery models: Effect of stenosis severity,plaque eccentricity,and ulceration

Clinical decision-making for the treatment of patients with diseased carotid artery is mainly based on the severity of the stenosis. However, stenosis severity alone is not a sensitive indicator, and other local factors for the assessment of stroke risk are required. Flow disturbance is of particular interest due to its proven association with increased thromboembolic activities. The objective of this study was to investigate the level of turbulence intensity (TI) with regards to certain geometrical features of the plaque – namely stenosis severity, eccentricity, and ulceration. A family of eight carotid-artery bifurcation models was examined using particle image velocimetry. Results showed a marked difference in turbulence intensity among these models; increasing degree of stenosis severity resulted in increased turbulence intensity, going from 0.12 m/s for mild stenosis to 0.37 m/s for severe stenosis (with concentric geometry). Moreover, independent of stenosis severity, eccentricity led to further elevations in turbulence intensity, increasing TI by 0.05–0.10 m/s over the counterpart concentric plaque. The presence of ulceration (in a 50% eccentric plaque) produced a larger portion of moderate turbulence intensity (~0.10 m/s) compared to the non-ulcerated model, more proximal to the bifurcation apex in the post-stenotic recirculation zone. The effect of plaque eccentricity and ulceration in enhancing the downstream turbulence has potential clinical implications for a more sensitive assessment of stroke risk beyond stenosis severity alone.     

Multilaboratory particle image velocimetry analysis of the FDA benchmark nozzle model to support validation of computational fluid dynamics simulations

This study is part of a FDA-sponsored project to evaluate the use and limitations of computational fluid dynamics (CFD) in assessing blood flow parameters related to medical device safety. In an interlaboratory study, fluid velocities and pressures were measured in a nozzle model to provide experimental validation for a companion round-robin CFD study. The simple benchmark nozzle model, which mimicked the flow fields in several medical devices, consisted of a gradual flow constriction, a narrow throat region, and a sudden expansion region where a fluid jet exited the center of the nozzle with recirculation zones near the model walls. Measurements of mean velocity and turbulent flow quantities were made in the benchmark device at three independent laboratories using particle image velocimetry (PIV). Flow measurements were performed over a range of nozzle throat Reynolds numbers (Re(throat)) from 500 to 6500, covering the laminar, transitional, and turbulent flow regimes. A standard operating procedure was developed for performing experiments under controlled temperature and flow conditions and for minimizing systematic errors during PIV image acquisition and processing. For laminar (Re(throat)=500) and turbulent flow conditions (Re(throat)≥3500), the velocities measured by the three laboratories were similar with an interlaboratory uncertainty of ～10% at most of the locations. However, for the transitional flow case (Re(throat)=2000), the uncertainty in the size and the velocity of the jet at the nozzle exit increased to ～60% and was very sensitive to the flow conditions. An error analysis showed that by minimizing the variability in the experimental parameters such as flow rate and fluid viscosity to less than 5% and by matching the inlet turbulence level between the laboratories, the uncertainties in the velocities of the transitional flow case could be reduced to ～15%. The experimental procedure and flow results from this interlaboratory study (available at http://fdacfd.nci.nih.gov) will be useful in validating CFD simulations of the benchmark nozzle model and in performing PIV studies on other medical device models.     

Flow pattern and shear stress distribution of distal end-to-side anastomoses. A comparison of the instantaneous velocity fields obtained by particle image velocimetry

OBJECTIVE: To describe the local hemodynamics and pressure losses of crural bypass anastomoses using instantaneous velocity fields acquired by particle image velocimetry (PIV). METHODS: Silastic models of a Taylor patch, a Miller cuff and a femoro-crural patch prosthesis (FCPP) were attached to a circuit driven by a Berlin Heart, providing a pulsatile flow with an amplitude of 450 to 25 ml/min (mean 200 ml/min). An outflow resistance of 0.5 mmHg/ml/min (peripheral resistance units, PRU) was modeled using small silastic tubes providing a phase shift of -12 degrees between flow and pressure curves. The working fluid consisted of a glycerine/water mixture with a viscosity of 4 mPas. Hollow glass spheres with a mean size of 9-13 microm were used as tracer particles. Instantaneous velocity fields were obtained by means of PIV and shear rates as well as shear stresses were calculated. Triggered by the flowmeter signal, 10 measurements at 100 ms intervals per cardiac cycle were obtained. The pressures were measured on the inflow and at both distal outflows. The resulting mean pressure losses due to flow separation and distal fluid acceleration were calculated. RESULTS: Inside the Taylor patch anastomosis a large flow separation at the hood containing a clockwise rotating vortex was found. Additionally a smaller flow separation at the heel and a flow stagnation zone on the floor of the recipient artery were observed. Conversely, inside the Miller cuff a counterclockwise rotating vortex was seen inside a large heel flow separation. The FCPP also showed typical separation areas at the hood and heel of the anastomosis, although these were smaller compared to the other anastomoses. Inside the FCPP anastomosis no vortex creation was observed throughout the cardiac cycle. The mainstream velocities at the inlet levels were comparable for the three anastomoses. A significant fluid acceleration was present at the antegrade as well as the retrograde outlets of the Taylor and Miller cuff, while the fluid acceleration at the antegrade outflow of the FCPP was small, which was attributed to the end-to-end configuration of the antegrade FCPP leg. The calculated normalized antegrade and retrograde pressure losses for the Taylor form were 0.90 and 0.88, for the Miller cuff 0.89 and 0.86 and for the FCPP 0.94 and 0.86, respectively. The shear stresses inside the flow separations of the three anastomoses were significantly lower than normal wall shear stresses. High shear stress levels were found inside the transition zones between flow separation and high velocity mainstream. CONCLUSIONS: The flow pattern inside cuffed or funnel shaped anastomoses consists of large flow separation zones, which are thought to be associated with intimal hyperplasia development. In addition, fluid accelerations at the distal outlets result in pressure losses, which may contribute to impaired crural perfusion.     

Intelligent image analysis (IIA) using artificial neural network (ANN) for non-invasive estimation of chlorophyll content in micropropagated plants of potato

An artificial neural network (ANN) was used to analyze photometric features extracted from the digitized images of leaves from in vitro-regenerated potato plants for non-invasive estimation of chlorophyll content. A MATLAB®-based, feed-forward, backpropagation-type network was developed for an input layer (three input elements), with one hidden layer (one node) and one output layer representing the predicted chlorophyll content. A significant influence of training function during optimization of ANN modeling was observed. Among the 11 training functions tested, “trainlm” was found to be the best on the basis of comparative analysis of root-mean-square error (RMSE) at zero epoch. A significant correlation between the model-predicted and Soil-Plant Analysis Development (SPAD) meter-measured relative chlorophyll contents was obtained when the mean brightness ratio (rgb) parameters were used. Compared to a red (R), green (G), and blue (B) color space model, the rgb model exhibited better performance with a significant correlation (R ² = 0.85). Incorporation of photometric features, such as luminosity (L), blue (B)/L, and green (G)/L, with rgb failed to improve the performance of the network. The developed Intelligent image analysis (IIA) system was able to estimate in real time the chlorophyll content of in vitro-regenerated leaves for assessment of plant nutrient status during micropropagation.     

A novel method for the estimation of soybean chlorophyll content using a smartphone and image analysis

The development of smartphones, specifically their cameras, and imaging technologies has enabled their use as sensors/measurement tools. Here we aimed to evaluate the applicability of a fast and noninvasive method for the estimation of total chlorophyll (Chl), Chl a, Chl b, and carotenoids (Car) content of soybean plants using a smartphone camera. Single leaf disc images were obtained using a smartphone camera. Subsequently, for the same leaf discs, a Chl meter was used to obtain the relative index of Chl and the photosynthetic pigments were then determined using a classic method. The RGB, HSB and CIELab color models were extracted from the smartphone images and correlated to Chl values obtained using a Chl meter and by a standard laboratory protocol. The smartphone camera was sensitive enough to capture successfully a broad range of Chl and Car contents seen in soybean leaves. Although there was a variation between color models, some of the proposed regressions (e.g., the S and b index from HSB and Lab color models and NRI [RGB model]) were very close to the Chl meter values. Based on our findings, smartphones can be used for rapid and accurate estimation of soybean and Car contents in soybean leaves.     

Placement of the alpha-sarcin loop within the 50S subunit: evidence derived using a photolabile oligodeoxynucleotide probe.

We report the synthesis of a radioactive, photolabile oligodeoxyribonucleotide probe and its exploitation in identifying 50S ribosomal subunit components neighboring the alpha-sarcin loop. The probe is complementary to 23S rRNA nt 2653-2674. Photolysis of the complex formed between the probe and 50S subunits leads to site-specific probe photoincorporation into proteins L2, the most highly labeled protein, L1, L15, L16 and L27, labeled to intermediate extents, and L5, L9, L17 and L24, each labeled to a minor extent. Portions of each of these proteins thus lie within 23 A of nt U2653. These results lead us to conclude that the alpha-sarcin loop is located at the base of the L1 projection within the 50S subunit. Such placement, near the peptidyl transferase center, provides a rationale for the extreme sensitivity of ribosomal function to cleavage of the alpha-sarcin loop.     

Reconstitution of biologically active 50S ribosomal subunits with artificial 5S RNA molecules carrying disturbances in the base pairing within the molecular stalk.

Bacillus stearothermophilus 50S ribosomal subunits were reconstituted in vitro using artificial 5S RNA molecules constructed by combining parts of major and minor type (Raué et al. (1976) Europ. J. Biochem. 68, 169-176) B. licheniformis 5S RNA. The artificial 5S RNA molecules carry defined disturbances (A.C juxtapositions and extra G.U pairs) in the base pairing between the 5'- and 3'-terminal sequences of the molecule (the molecular stalk region). The biological activity of the reconstituted subunits was determined in an E. coli cell-free system programmed with poly-U. The results show that conservation of the base pairing within the molecular stalk is not required for biological activity of 5S RNA. Disturbances of the base pairing within this region do reduce the rate of reconstitution, however. Normal base pairing in the molecular stalk is thus required to ensure efficient ribosome assembly.     

Correlation of atrial and ventricular contractility. The possibility of using the atria for controlling the artificial heart]

A possibility of the artificial ventricle regulation using information on the atrial contraction is discussed. Correlations between atrial wall tension and ventricle tension as well as the pressure under radom heart rhythm variations were studied. The wall tension of the appropriate heart portion was recorded by the arched tensiometers, and intraventricular pressure under cavity catheterization by the electric manometer. The cardiac rhythm varied between 2.0 c-1 and 4.0 c-1. A correlation between the tension and interpulse interval, (correlation coefficient being 0.62 +/- 0.05) was established. The close relationships "atrial wall tension-intramuscular pressure" were observed. For such relationships the correlation coefficient changed from 0.713 +/- 0,09 to 0,874 +/- 0.02 depending on the average value of interpulse interval. It is concluded that the information on the atrial contraction can be used for artificial heart regulation.     

A quantitative 'in-situ' hybridization method using computer-assisted image analysis. Validation and measurement of atrial-natriuretic-factor mRNA in the rat heart.

We have validated a quantitative 'in-situ' hybridization method and computer-assisted image analysis of autoradiographs as a technique for measuring atrial-natriuretic-factor (ANF) mRNA in tissue sections of rat heart by: (i) producing radioactive standards to calibrate the autoradiograms and (ii) assessing: (a) specificity (through RNAase A background subtraction, comparison of ANF mRNA and non-ANF mRNA probe binding to sections, Northern analysis and section-thickness titration curves); (b) sensitivity (by calculating the limit of detection for ventricular levels of ANF mRNA); (c) precision [inter-assay CV (coefficient of variation) less than 10%; intra-assay CV 6-7%]; and (d) accuracy. We have found with this technique that deoxycortone and saline treatment of rats elevates ANF mRNA to a larger extent in the ventricles than in the atria and that, in neonatal-rat hearts, ANF mRNA is elevated in all cardiac chambers relative to adult levels.     

In situ variability of carrageenan content and biomass in the cultivated red macroalga <Emphasis Type="Italic">Kappaphycus alvarezii</Emphasis> with an estimation of its carrageenan stock at the scale of the Malasoro Bay (Indonesia) using satellite image processing

The combination of in situ measurements, hydrodynamical modeling, and satellite imagery processing presents a complete tool to improve seaweed culture management. This study measured the evolution of carrageenan content during 1-year period and estimated the biomass of Kappaphycus alvarezii at a bay scale in Malasoro Bay, Indonesia. It allowed the determination of the carrageenan stock at the scale of the bay. The carrageenan content was assessed from different parts, i.e., basal and apical, of the thallus. The biomass from T ₀ (beginning of the cultivation) to T ₄₅ (harvesting time) was determined at two seasons. Satellite image processing was performed to estimate the biomass at bay scale using parcels resulting from a semi-automatic delineation process. As no long-term in situ environmental data existed in the study area, a 3D hydrodynamical model (Model for Application at Regional Scale 3D) was implemented to obtain high-frequency salinity, water temperature, and currents. The high carrageenan yield was obtained from April to September 2015 (54–63%) when temperature ranged from 25 to 30 °C, salinity range from 33.8 to 34.8 psu, and precipitation below 0.5 mm. The biomass at bay scale was estimated at 2590 t with the highest carrageenan stock at 1.8 t ha^?1 in May 2015. The carrageenan yields and quality obtained in this study fulfilled the specification recommended for industry, and the cultivation of K. alvarezii can be carried out at optimal environmental conditions on April–September. The comprehensive approaches used in this study provide information for carrageenan stock and seaweed culture management as an important economic activity to support Indonesian coastal communities.