Hydrodynamic phenomena in the downcomer during flow rate transients in the primary circuit of a PWR

Mixing phenomena observed when the flow rate in a single loop of the primary circuit is changed can influence the operation of pressurized water reactor (PWR) by inducing local gradients of boron concentration or coolant temperature. Analysis of one-dimensional Laser Doppler Anemometry (LDA) measurements during the start-up and shutdown of pump on a single loop of the ROCOM test facility has been performed. The effect of a step change and a ramped change in the flow rate on the axial and azimuthal velocities was examined. Numerical simulations were also performed for the step change in the flow rate that gave quantitative agreement with the axial velocities. Phenomenological agreement was made on the turbulent kinetic energy; however, observed values were a factor of 2.5 less than the turbulent kinetic energy derived from the measurements. © 2007.

Modelling and measurements on a condensing counter-flow heat exchanger

Experimental investigations and computer modelling studies have been made on the refrigerant-water counterflow condenser section of a small air to water heat pump. The main object of the investigation was a comparative study between the computer modelling predictions and the experimental observations for a range of operating conditions but other characteristics of a counterflow heat exchanger are also discussed. The counterflow condenser consisted of 15 metres of a thermally coupled pair of copper pipes, one containing the R12 working fluid and the other water flowing in the opposite direction. This condenser was mounted horizontally and folded into 0.5 metre straight sections. Thermocouples were inserted in both pipes at one metre intervals and transducers for pressure and flow measurement were also included. Data acquisition, storage and analysis was carried out by a micro-computer suitably interfaced with the transducers and thermocouples. Many sets of readings were taken under a variety of conditions, with air temperature ranging from 18 to 26 degrees Celsius, water inlet from 13.5 to 21.7 degrees, R12 inlet temperature from 61.2 to 81.7 degrees and water mass flow rate from 6.7 to 32.9 grammes per second. A Fortran computer model of the condenser (originally prepared by Carrington[1]) has been modified to match the information available from experimental work. This program uses iterative segmental integration over the desuperheating, mixed phase and subcooled regions for the R12 working fluid, the water always being in the liquid phase. Methods of estimating the inlet and exit fluid conditions from the available experimental data have been developed for application to the model. Temperature profiles and other parameters have been predicted and compared with experimental values for the condenser for a range of evaporator conditions and have shown that the model gives a satisfactory prediction of the physical behaviour of a simple counterflow heat exchanger in both single phase and two phase regions.

Computer simulation of liquid flow patterns on distillation trays

This thesis describes work carried out to improve the fundamental modelling of liquid flows on distillation trays. A mathematical model is presented based on the principles of computerised fluid dynamics. It models the liquid flow in the horizontal directions allowing for the effects of the vapour through the use of an increased liquid turbulence, modelled by an eddy viscosity, and a resistance to liquid flow caused by the vapour being accelerated horizontally by the liquid. The resultant equations are similar to the Navier-Stokes equations with the addition of a resistance term.A mass-transfer model is used to calculate liquid concentration profiles and tray efficiencies. A heat and mass transfer analogy is used to compare theoretical concentration profiles to experimental water-cooling data obtained from a 2.44 metre diameter air-water distillation simulation rig. The ratios of air to water flow rates are varied in order to simulate three pressures: vacuum, atmospheric pressure and moderate pressure.For simulated atmospheric and moderate pressure distillation, the fluid mechanical model constantly over-predicts tray efficiencies with an accuracy of between +1.7% and +11.3%. This compares to -1.8% to -10.9% for the stagnant regions model (Porter et al. 1972) and +12.8% to +34.7% for the plug flow plus back-mixing model (Gerster et al. 1958). The model fails to predict the flow patterns and tray efficiencies for vacuum simulation due to the change in the mechanism of liquid transport, from a liquid continuous layer to a spray as the liquid flow-rate is reduced. This spray is not taken into account in the development of the fluid mechanical model. A sensitivity analysis carried out has shown that the fluid mechanical model is relatively insensitive to the prediction of the average height of clear liquid, and a reduction in the resistance term results in a slight loss of tray efficiency. But these effects are not great. The model is quite sensitive to the prediction of the eddy viscosity term. Variations can produce up to a 15% decrease in tray efficiency. The fluid mechanical model has been incorporated into a column model so that statistical optimisation techniques can be employed to fit a theoretical column concentration profile to experimental data. Through the use of this work mass-transfer data can be obtained.

The effect of innate compliance on the performance of a counterpulsation device

Cardiovascular disease (CVD) continues to be one of the top causes of mortality in the world. World Heart Organization (WHO) reported that in 2004, CVD contributed to almost 30% of death from estimated worldwide death figures of 58 million[1]. Heart failure treatment varies from lifestyle adjustment to heart transplantation; its aims are to reduce HF symptoms, prolong patient survival and minimize risk [2]. One alternative available in the market for HF treatment is Left Ventricular Assist Device (LVAD). Chronic Intermittent Mechanical Support (CIMS) device is a novel (LVAD) heart failure treatment using counterpulsation similar to Intra Aortic Balloon Pumps (IABP). However, the implantation site of the CIMS balloon is in the ascending aorta just distal to aortic valve contrasted with IABP in the descending aorta. Counterpulsation coupled with implantation close to the aortic valve enables comparable flow augmentation with reduced balloon volume. Two prototypes of the CIMS balloon were constructed using rapid prototyping: the straight-body model is a cylindrical tube with a silicone membrane lining with zero expansive compliance. The compliant-body model had a bulging structure that allowed the membrane to expand under native systolic pressure increasing the device’s static compliance to 1.5 mL/mmHg. This study examined the effect of device compliance and vascular compliance on counterpulsating flow augmentation. Both prototypes were tested on a two-element Windkessel model human mock circulatory loop (MCL). The devices were placed just distal to aortic valve and left coronary artery. The MCL mimicked HF with cardiac output of 3 L/min, left ventricular pressure of 85/15 mmHg, aortic pressure of 70/50 mmHg and left coronary artery flow rate of 66 mL/min. The mean arterial pressure (MAP) was calculated to be 57 mmHg. Arterial compliance was set to be1.25 mL/mmHg and 2.5 mL/mmHg. Inflation of the balloon was triggered at the dicrotic notch while deflation was at minimum aortic pressure prior to systole. Important haemodynamics parameters such as left ventricular pressure (LVP), aortic pressure (AoP), cardiac output (CO), left coronary artery flowrate (QcorMean), and dP (Peak aortic diastolic augmentation pressure – AoPmax ) were simultaneously recorded for both non-assisted mode and assisted mode. ANOVA was used to analyse the effect of both factors (balloon and arterial compliance) to flow augmentation. The results showed that for cardiac output and left coronary artery flowrate, there were significant difference between balloon and arterial compliance at p < 0.001. Cardiac output recorded maximum output at 18% for compliant body and stiff arterial compliance. Left coronary artery flowrate also recorded around 20% increase due to compliant body and stiffer arterial compliance. Resistance to blood ejection recorded highest difference for combination of straight body and stiffer arterial compliance. From these results it is clear that both balloon and arterial compliance are statistically significant factors for flow augmentation on peripheral artery and reduction of resistance. Although the result for resistance reduction was different from flow augmentation, these results serves as an important aspect which will influence the future design of the CIMS balloon and its control strategy. References: 1. Mathers C, Boerma T, Fat DM. The Global Burden of disease:2004 update. Geneva: World Heatlh Organization; 2008. 2. Jessup M, Brozena S. Heart Failure. N Engl J Med 2003;348:2007-18.

Flow across microvessel walls through the endothelial surface glycocalyx and the interendothelial cleft

A mathematical model is presented for steady fluid flow across microvessel walls through a serial pathway consisting of the endothelial surface glycocalyx and the intercellular cleft between adjacent endothelial cells, with junction strands and their discontinuous gaps. The three-dimensional flow through the pathway from the vessel lumen to the tissue space has been computed numerically based on a Brinkman equation with appropriate values of the Darcy permeability. The predicted values of the hydraulic conductivity Lp, defined as the ratio of the flow rate per unit surface area of the vessel wall to the pressure drop across it, are close to experimental measurements for rat mesentery microvessels. If the values of the Darcy permeability for the surface glycocalyx are determined based on the regular arrangements of fibres with 6nm radius and 8nm spacing proposed recently from the detailed structural measurements, then the present study suggests that the surface glycocalyx could be much less resistant to flow compared to previous estimates by the one-dimensional flow analyses, and the intercellular cleft could be a major determinant of the hydraulic conductivity of the microvessel wall.