An extended computational model of the circulatory system for designing ventricular assist devices.

An extended computational model of the circulatory system has been developed to predict blood flow in the presence of ventricular assist devices (VADs). A novel VAD, placed in the descending aorta, intended to offload the left ventricle (LV) and augment renal perfusion is being studied. For this application, a better understanding of the global hemodynamic response of the VAD, in essence an electrically driven pump, and the cardiovascular system is necessary. To meet this need, a model has been established as a nonlinear, lumped-parameter electrical analog, and simulated results under different states [healthy, congestive heart failure (CHF), and postinsertion of VAD] are presented. The systemic circulation is separated into five compartments and the descending aorta is composed of three components to accurately yield the system response of each section before and after the insertion of the VAD. Delays in valve closing time and blood inertia in the aorta were introduced to deliver a more realistic model. Pump governing equations and optimization are based on fundamental theories of turbomachines and can serve as a practical initial design point for rotary blood pumps. The model's results closely mimic established parameters for the circulatory system and confirm the feasibility of the intra-aortic VAD concept. This computational model can be linked with models of the pump motor to provide a valuable tool for innovative VAD design.

A controller for a miniature intra-aortic ventricular assist device.

Two control algorithms have been developed for a minimally invasive axial-flow ventricular assist device (VAD) for placement in the descending aorta. The purpose of the device is to offload the left ventricle and to augment lower body perfusion in patients with moderate congestive heart failure. The VAD consists of an intra-aortic impeller with a built-in permanent magnet rotor and an extra-aortic stator. The control algorithms, which use pressure readings upstream and downstream of the VAD to determine the pump status, have been tested in a mock circulatory system under two conditions, namely with or without afterload sensitivity. The results give an insight into controller design for an intra-aortic blood pump working in series with the heart.

Design method of an innovative motor for an intra-aortic ventricular assist device

A permanent-magnet motor has been designed for an innovative axial-flow ventricular assist device (VAD), to be placed in the descending aorta, intended to offload the left ventricle and augment renal perfusion in patients with congestive heart failure (CHF). For this application, an intra-aortic impeller with a built-in permanent magnet rotor is driven by an extraaortic stator working in synchronism with the natural heart. To meet this need, a two-dimensional analytical model has been developed in the MATLAB environment to estimate machine parameters; finite element analysis (FEA) has been used to refine the results. A prototype blood pump equipped with an innovative motor designed from the procedure above has been tested in a mock loop representing the human circulatory system. The performance of VAD incorporating the motor is presented. © 2009 IEEE.

Effects of turbulent reynolds number on the displacement speed statistics in the thin reaction zones regime of turbulent premixed combustion

The effects of turbulent Reynolds number on the statistical behaviour of the displacement speed have been studied using three-dimensional Direct Numerical Simulation of statistically planar turbulent premixed flames. The probability of finding negative values of the displacement speed is found to increase with increasing turbulent Reynolds number when the Damkhler number is held constant. It has been shown that the statistical behaviour of the Surface Density Function, and its strain rate and curvature dependence, plays a key role in determining the response of the different components of displacement speed. Increasing the turbulent Reynolds number is shown to reduce the strength of the correlations between tangential strain rate and dilatation rate with curvature, although the qualitative nature of the correlations remains unaffected. The dependence of displacement speed on strain rate and curvature is found to weaken with increasing turbulent Reynolds number when either Damkhler or Karlovitz number is held constant, but the qualitative nature of the correlation remains unaltered. The implications of turbulent Reynolds number effects in the context of Flame Surface Density (FSD) modelling have also been addressed, with emphasis on the influence of displacement speed on the curvature and propagation terms in the FSD balance equation. © 2011 Nilanjan Chakraborty et al.

Statistics and modelling of turbulent kinetic energy transport in different regimes of premixed combustion

The statistical behaviour of turbulent kinetic energy transport in turbulent premixed flames is analysed using data from three-dimensional Direct Numerical Simulation (DNS) of freely propagating turbulent premixed flames under decaying turbulence. For flames within the corrugated flamelets regime, it is observed that turbulent kinetic energy is generated within the flame brush. By contrast, for flames within the thin reaction zones regime it has been found that the turbulent kinetic energy decays monotonically through the flame brush. Similar trends are observed also for the dissipation rate of turbulent kinetic energy. Within the corrugated flamelets regime, it is demonstrated that the effects of the mean pressure gradient and pressure dilatation within the flame are sufficient to overcome the effects of viscous dissipation and are responsible for the observed augmentation of turbulent kinetic energy in the flame brush. In the thin reaction zones regime, the effects of the mean pressure gradient and pressure dilatation terms are relatively much weaker than those of viscous dissipation, resulting in a monotonic decay of turbulent kinetic energy across the flame brush. The modelling of the various unclosed terms of the turbulent kinetic energy transport equation has been analysed in detail. The predictions of existing models are compared with corresponding quantities extracted from DNS data. Based on this a-priori DNS assessment, either appropriate models are identified or new models are proposed where necessary. It is shown that the turbulent flux of turbulent kinetic energy exhibits counter-gradient (gradient) transport wherever the turbulent scalar flux is counter-gradient (gradient) in nature. A new model has been proposed for the turbulent flux of turbulent kinetic energy, and is found to capture the qualitative and quantitative behaviour obtained from DNS data for both the corrugated flamelets and thin reaction zones regimes without the need to adjust any of the model constants. © 2010 Springer Science+Business Media B.V.

Scalar dissipation rate transport in the corrugated flamelets and thin reaction zones regime for turbulent premixed flames

The characteristics of the scalar dissipation rate transport in the corrugated flamelets and the thin reaction zones regimes are studied based on two three-dimensional Direct Numerical Simulation (DNS) databases for freely propagating statistically planar turbulent premixed flames. The turbulent flame parameters are so chosen that the database which represents the corrugated flamelets regime has a global Damköhler number Da>1 whereas the database representing the thin reaction zones regime has Da <1. It is demonstrated that the terms originating from the correlation between fluctuating velocity and scalar gradient T1 shows strong Da dependence. The terms originating from dilatation T2, the scalar inner product of gradients of velocity and scalar fields T3 and the correlation between reaction rate and scalar gradients T4 and the dissipation term D2 remain important for both the flames. However, T3 dissipates scalar dissipation rate in the Da > 1 flame while it produces scalar dissipation rate in the Da < 1 flame. This difference is because of the change in the alignment between scalar and velocity gradients

Effects of flow geometry and Damkohler number on turbulent premixed flame

DNS of planar turbulent flame and turbulent V-flame has been conducted to investigate turbulence-scalar interaction in relatively practical turbulent combustion. Several turbulence quantities are examined for the understandings of fundamental characteristics of flow field in V-flame. Due to the additional turbulence production by the hot-rod, turbulence does not simply decay in V-flame. Turbulence-scalar interaction, scalar alignments with the principal strain rate in other words, is then clarified. The competition of turbulence and dilatation can be found in the conditional PDF of flame normal alignment. The results suggests that the alignment characteristics in high Da flames are applicable to low Da flames in the region of intense heat release.