Coolant Flow and Heat Transfer in Pebble-Bed Reactor Core

This thesis gathers knowledge about ongoing high-temperature reactor projects around the world. Methods for calculating coolant flow and heat transfer inside a pebble-bed reactor core are also developed. The thesis begins with the introduction of high-temperature reactors including the current state of the technology. Process heat applications that could use the heat from a high-temperature reactor are also introduced. A suitable reactor design with data available in literature is selected for the calculation part of the thesis. Commercial computational fluid dynamics software Fluent is used for the calculations. The pebble-bed is approximated as a packed-bed, which causes sink terms to the momentum equations of the gas flowing through it. A position dependent value is used for the packing fraction. Two different models are used to calculate heat transfer. First a local thermal equilibrium is assumed between the gas and solid phases and a single energy equation is used. In the second approach, separate energy equations are used for the phases. Information about steady state flow behavior, pressure loss, and temperature distribution in the core is obtained as results of the calculations. The effect of inlet mass flow rate to pressure loss is also investigated. Data found in literature and the results correspond each other quite well, considered the amount of simplifications in the calculations. The models developed in this thesis can be used to solve coolant flow and heat transfer in a pebble-bed reactor, although additional development and model validation is needed for better accuracy and reliability.

Process intensification: From optimised flow patterns to microprocess technology

This thesis is focused on process intensification. Several significant problems and applications of this theme are covered. Process intensification is nowadays one of the most popular trends in chemical engineering and attempts have been made to develop a general, systematic methodology for intensification. This seems, however, to be very difficult, because intensified processes are often based on creativity and novel ideas. Monolith reactors and microreactors are successful examples of process intensification. They are usually multichannel devices in which a proper feed technique is important for creating even fluid distribution into the channels. Two different feed techniques were tested for monoliths. In the first technique a shower method was implemented by means of perforated plates. The second technique was a dispersion method using static mixers. Both techniques offered stable operation and uniform fluid distribution. The dispersion method enabled a wider operational range in terms of liquid superficial velocity. Using dispersion method, a volumetric gas-liquid mass transfer coefficient of 2 s-1 was reached. Flow patterns play a significant role in terms of the mixing performance of micromixers. Although the geometry of a T-mixer is simple, channel configurations and dimensions had a clear effect on mixing efficiency. The flow in the microchannel was laminar, but the formation of vortices promoted mixing in micro T-mixers. The generation of vortices was dependent on the channel dimensions, configurations and flow rate. Microreactors offer a high ratio of surface area to volume. Surface forces and interactions between fluids and surfaces are, therefore, often dominant factors. In certain cases, the interactions can be effectively utilised. Different wetting properties of solid materials (PTFE and stainless steel) were applied in the separation of immiscible liquid phases. A micro-scale plate coalescer with hydrophilic and hydrophobic surfaces was used for the continuous separation of organic and aqueous phases. Complete phase separation occurred in less than 20 seconds, whereas the separation time by settling exceeded 30 min. Fluid flows can be also intensified in suitable conditions. By adding certain additives into turbulent fluid flow, it was possible to reduce friction (drag) by 40 %. Drag reduction decreases frictional pressure drop in pipelines which leads to remarkable energy savings and decreases the size or number of pumping facilities required, e.g., in oil transport pipes. Process intensification enables operation often under more optimal conditions. The consequent cost savings from reduced use of raw materials and reduced waste lead to greater economic benefits in processing.

CFD simulation of complex phenomena containing suspensions and flow throughporous media

There is an increasing reliance on computers to solve complex engineering problems. This is because computers, in addition to supporting the development and implementation of adequate and clear models, can especially minimize the financial support required. The ability of computers to perform complex calculations at high speed has enabled the creation of highly complex systems to model real-world phenomena. The complexity of the fluid dynamics problem makes it difficult or impossible to solve equations of an object in a flow exactly. Approximate solutions can be obtained by construction and measurement of prototypes placed in a flow, or by use of a numerical simulation. Since usage of prototypes can be prohibitively time-consuming and expensive, many have turned to simulations to provide insight during the engineering process. In this case the simulation setup and parameters can be altered much more easily than one could with a real-world experiment. The objective of this research work is to develop numerical models for different suspensions (fiber suspensions, blood flow through microvessels and branching geometries, and magnetic fluids), and also fluid flow through porous media. The models will have merit as a scientific tool and will also have practical application in industries. Most of the numerical simulations were done by the commercial software, Fluent, and user defined functions were added to apply a multiscale method and magnetic field. The results from simulation of fiber suspension can elucidate the physics behind the break up of a fiber floc, opening the possibility for developing a meaningful numerical model of the fiber flow. The simulation of blood movement from an arteriole through a venule via a capillary showed that the model based on VOF can successfully predict the deformation and flow of RBCs in an arteriole. Furthermore, the result corresponds to the experimental observation illustrates that the RBC is deformed during the movement. The concluding remarks presented, provide a correct methodology and a mathematical and numerical framework for the simulation of blood flows in branching. Analysis of ferrofluids simulations indicate that the magnetic Soret effect can be even higher than the conventional one and its strength depends on the strength of magnetic field, confirmed experimentally by Völker and Odenbach. It was also shown that when a magnetic field is perpendicular to the temperature gradient, there will be additional increase in the heat transfer compared to the cases where the magnetic field is parallel to the temperature gradient. In addition, the statistical evaluation (Taguchi technique) on magnetic fluids showed that the temperature and initial concentration of the magnetic phase exert the maximum and minimum contribution to the thermodiffusion, respectively. In the simulation of flow through porous media, dimensionless pressure drop was studied at different Reynolds numbers, based on pore permeability and interstitial fluid velocity. The obtained results agreed well with the correlation of Macdonald et al. (1979) for the range of actual flow Reynolds studied. Furthermore, calculated results for the dispersion coefficients in the cylinder geometry were found to be in agreement with those of Seymour and Callaghan.

Effect of meal ingestion on liver stiffness in patients with cirrhosis and portal hypertension.

BACKGROUND AND AIMS: Liver stiffness is increasingly used in the non-invasive evaluation of chronic liver diseases. Liver stiffness correlates with hepatic venous pressure gradient (HVPG) in patients with cirrhosis and holds prognostic value in this population. Hence, accuracy in its measurement is needed. Several factors independent of fibrosis influence liver stiffness, but there is insufficient information on whether meal ingestion modifies liver stiffness in cirrhosis. We investigated the changes in liver stiffness occurring after the ingestion of a liquid standard test meal in this population. METHODS: In 19 patients with cirrhosis and esophageal varices (9 alcoholic, 9 HCV-related, 1 NASH; Child score 6.9±1.8), liver stiffness (transient elastography), portal blood flow (PBF) and hepatic artery blood flow (HABF) (Doppler-Ultrasound) were measured before and 30 minutes after receiving a standard mixed liquid meal. In 10 the HVPG changes were also measured. RESULTS: Post-prandial hyperemia was accompanied by a marked increase in liver stiffness (+27±33%; p<0.0001). Changes in liver stiffness did not correlate with PBF changes, but directly correlated with HABF changes (r = 0.658; p = 0.002). After the meal, those patients showing a decrease in HABF (n = 13) had a less marked increase of liver stiffness as compared to patients in whom HABF increased (n = 6; +12±21% vs. +62±29%,p<0.0001). As expected, post-prandial hyperemia was associated with an increase in HVPG (n = 10; +26±13%, p = 0.003), but changes in liver stiffness did not correlate with HVPG changes. CONCLUSIONS: Liver stiffness increases markedly after a liquid test meal in patients with cirrhosis, suggesting that its measurement should be performed in standardized fasting conditions. The hepatic artery buffer response appears an important factor modulating postprandial changes of liver stiffness. The post-prandial increase in HVPG cannot be predicted by changes in liver stiffness.

Determination of sulfate in the wet-process of phosphoric acid by reverse flow injection

An improved method based on reverse flow injection is proposed for determining sulfate concentration in the wet-process of phosphoric acid (WPA). The effect of reagent composition, flow rate, temperature, acid concentration, length of the reaction coil, and linear response range on the flow system is discussed in detail. Optimal conditions are established for determining sulfate in the WPA samples. Baseline drift is avoided by a periodic washing step with EDTA in an alkaline medium. A linear response is observed within a range of 20 - 360 mg L-1, given by the equation A = 0.0020C (mg L-1) + 0.0300, R² = 0.9991. The detection limit of the proposed method for sulfate analysis is 3 mg L-1, and the relative standard deviation (n = 12) of sulfate absorbance peak is less than 1.60%. This method has a rate of up to 29 samples per hour, and the results compare well with those obtained with gravimetric method.

Numerical modelling of small supersonic axial flow turbines

Supersonic axial turbine stages typically exhibit lower efficiencies than subsonic axial turbine stages. One reason for the lower efficiency is the occurrence of shock waves. With higher pressure ratios the flow inside the turbine becomes relatively easily supersonic if there is only one turbine stage. Supersonic axial turbines can be designed in smaller physical size compared to subsonic axial turbines of same power. This makes them good candidates for turbochargers in large diesel engines, where space can be a limiting factor. Also the production costs are lower for a supersonic axial turbine stage than for two subsonic stages. Since supersonic axial turbines are typically low reaction turbines, they also create lower axial forces to be compensated with bearings compared to high reaction turbines. The effect of changing the stator-rotor axial gap in a small high (rotational) speed supersonic axial flow turbine is studied in design and off-design conditions. Also the effect of using pulsatile mass flow at the supersonic stator inlet is studied. Five axial gaps (axial space between stator and rotor) are modeled using threedimensional computational fluid dynamics at the design and three axial gaps at the off-design conditions. Numerical reliability is studied in three independent studies. An additional measurement is made with the design turbine geometry at intermediate off-design conditions and is used to increase the reliability of the modelling. All numerical modelling is made with the Navier-Stokes solver Finflo employing Chien’s k ¡ ² turbulence model. The modelling of the turbine at the design and off-design conditions shows that the total-to-static efficiency of the turbine decreases when the axial gap is increased in both design and off-design conditions. The efficiency drops almost linearily at the off-design conditions, whereas the efficiency drop accelerates with increasing axial gap at the design conditions. The modelling of the turbine stator with pulsatile inlet flow reveals that the mass flow pulsation amplitude is decreased at the stator throat. The stator efficiency and pressure ratio have sinusoidal shapes as a function of time. A hysteresis-like behaviour is detected for stator efficiency and pressure ratio as a function of inlet mass flow, over one pulse period. This behaviour arises from the pulsatile inlet flow. It is important to have the smallest possible axial gap in the studied turbine type in order to maximize the efficiency. The results for the whole turbine can also be applied to some extent in similar turbines operating for example in space rocket engines. The use of a supersonic stator in a pulsatile inlet flow is shown to be possible.

Pulsatile Blood Flow Simulations in Computed Tomography(CT) Scan-Based and Idealized Geometries of Human Aorta

Blood flow in human aorta is an unsteady and complex phenomenon. The complex patterns are related to the geometrical features like curvature, bends, and branching and pulsatile nature of flow from left ventricle of heart. The aim of this work was to understand the effect of aorta geometry on the flow dynamics. To achieve this, 3D realistic and idealized models of descending aorta were reconstructed from Computed Tomography (CT) images of a female patient. The geometries were reconstructed using medical image processing code. The blood flow in aorta was assumed to be laminar and incompressible and the blood was assumed to be Newtonian fluid. A time dependent pulsatile and parabolic boundary condition was deployed at inlet. Steady and unsteady blood flow simulations were performed in real and idealized geometries of descending aorta using a Finite Volume Method (FVM) code. Analysis of Wall Shear Stress (WSS) distribution, pressure distribution, and axial velocity profiles were carried out in both geometries at steady and unsteady state conditions. The results obtained in thesis work reveal that the idealization of geometry underestimates the values of WSS especially near the region with sudden change of diameter. However, the resultant pressure and velocity in idealized geometry are close to those in real geometry

Simulating Fluid Flow in a Compressed Valve

The objective of the work is to study fluid flow behavior through a pinch valve and to estimate the flow coefficient (KV ) at different opening positions of the valve. The flow inside a compressed valve is more complex than in a straight pipe, and it is one of main topics of interest for engineers in process industry. In the present work, we have numerically simulated compressed valve flow at different opening positions. In order to simulate the flow through pinch valve, several models of the elastomeric valve tube (pinch valve tube) at different opening positions were constructed in 2D-axisymmetric and 3D geometries. The numerical simulations were performed with the CFD packages; ANSYS FLUENT and ANSYS CFX by using parallel computing. The distributions of static pressure, velocity and turbulent kinetic energy have been studied at different opening positions of the valve in both 2D-axisymmetric and 3D experiments. The flow coefficient (KV ) values have been measured at different valve openings and are compared between 2D-axisymmetric and 3D simulation results.

Numerical Investigation of Flow Past a Circular Cylinder and in a Staggered Tube Bundle Using Various Turbulence Models

Transitional flow past a three-dimensional circular cylinder is a widely studied phenomenon since this problem is of interest with respect to many technical applications. In the present work, the numerical simulation of flow past a circular cylinder, performed by using a commercial CFD code (ANSYS Fluent 12.1) with large eddy simulation (LES) and RANS (κ - ε and Shear-Stress Transport (SST) κ - ω! model) approaches. The turbulent flow for ReD = 1000 & 3900 is simulated to investigate the force coefficient, Strouhal number, flow separation angle, pressure distribution on cylinder and the complex three dimensional vortex shedding of the cylinder wake region. The numerical results extracted from these simulations have good agreement with the experimental data (Zdravkovich, 1997). Moreover, grid refinement and time-step influence have been examined. Numerical calculations of turbulent cross-flow in a staggered tube bundle continues to attract interest due to its importance in the engineering application as well as the fact that this complex flow represents a challenging problem for CFD. In the present work a time dependent simulation using κ – ε, κ - ω! and SST models are performed in two dimensional for a subcritical flow through a staggered tube bundle. The predicted turbulence statistics (mean and r.m.s velocities) have good agreement with the experimental data (S. Balabani, 1996). Turbulent quantities such as turbulent kinetic energy and dissipation rate are predicted using RANS models and compared with each other. The sensitivity of grid and time-step size have been analyzed. Model constants sensitivity study have been carried out by adopting κ – ε model. It has been observed that model constants are very sensitive to turbulence statistics and turbulent quantities.

A pH optode based on thymol blue: application to determination of CO2 using flow injection analysis system

An optode based on thymol blue (TB), an acid-based indicator, has been constructed and evaluated as a detector in FIA system for CO2 determination. The dye was chemically immobilised on the surface of a bifurcated glass optical fibre bundle, using silanisation in organic media. In FIA system, hydrogen carbonate or carbonate samples are injected in a buffer carrier solution, and then are mixed with phosphoric acid solution to generate CO2, which diffuses through a PTFE membrane, in order to be collected in an acceptor carrier fluid, pumped towards to detection cell, in which the optode was adapted. The proposed system presents two linear response ranges, from 1.0 x 10-3 to 1.0 x 10-2 mol l-1, and from 2.0 x 10-2 to 0.10 mol l-1. The sampling frequency was 11 sample h-1, with good repeatability (R.S.D < 4 %, n = 10). In flow conditions the optode lifetime was 170 h. The system was applied in the analysis of commercial mineral water and the results obtained in the hydrogen carbonate determination did not differ significantly from those obtained by potentiometry, at a confidence level of 95 %.

CFD Modelling of Direct Contact Condensation in Suppression Pools by Applying Condensation Models of Separated Flow

The condensation rate has to be high in the safety pressure suppression pool systems of Boiling Water Reactors (BWR) in order to fulfill their safety function. The phenomena due to such a high direct contact condensation (DCC) rate turn out to be very challenging to be analysed either with experiments or numerical simulations. In this thesis, the suppression pool experiments carried out in the POOLEX facility of Lappeenranta University of Technology were simulated. Two different condensation modes were modelled by using the 2-phase CFD codes NEPTUNE CFD and TransAT. The DCC models applied were the typical ones to be used for separated flows in channels, and their applicability to the rapidly condensing flow in the condensation pool context had not been tested earlier. A low Reynolds number case was the first to be simulated. The POOLEX experiment STB-31 was operated near the conditions between the ’quasi-steady oscillatory interface condensation’ mode and the ’condensation within the blowdown pipe’ mode. The condensation models of Lakehal et al. and Coste & Lavi´eville predicted the condensation rate quite accurately, while the other tested ones overestimated it. It was possible to get the direct phase change solution to settle near to the measured values, but a very high resolution of calculation grid was needed. Secondly, a high Reynolds number case corresponding to the ’chugging’ mode was simulated. The POOLEX experiment STB-28 was chosen, because various standard and highspeed video samples of bubbles were recorded during it. In order to extract numerical information from the video material, a pattern recognition procedure was programmed. The bubble size distributions and the frequencies of chugging were calculated with this procedure. With the statistical data of the bubble sizes and temporal data of the bubble/jet appearance, it was possible to compare the condensation rates between the experiment and the CFD simulations. In the chugging simulations, a spherically curvilinear calculation grid at the blowdown pipe exit improved the convergence and decreased the required cell count. The compressible flow solver with complete steam-tables was beneficial for the numerical success of the simulations. The Hughes-Duffey model and, to some extent, the Coste & Lavi´eville model produced realistic chugging behavior. The initial level of the steam/water interface was an important factor to determine the initiation of the chugging. If the interface was initialized with a water level high enough inside the blowdown pipe, the vigorous penetration of a water plug into the pool created a turbulent wake which invoked the chugging that was self-sustaining. A 3D simulation with a suitable DCC model produced qualitatively very realistic shapes of the chugging bubbles and jets. The comparative FFT analysis of the bubble size data and the pool bottom pressure data gave useful information to distinguish the eigenmodes of chugging, bubbling, and pool structure oscillations.

Application of response surface methodology to study the biological removal of nitrogen from effluent of cattle slaughterhouse in a sequencing batch reactor

The aim of this study was to evaluate the efficiency of a sequencing batch reactor (SBR) on biological removal of nitrogen from cattle slaughterhouse wastewater by nitrification/denitrification processes. The effects of initial concentration of ammoniacal nitrogen were investigated at 100; 150 and 200 mg L-1 and air flow rate at 0.125; 0.375 and 0.625 L min¹ Lreactor-1 on the nitrogen compounds removal, by a Central Composite Rotational Design (CCRD) configuration. There were variations from 9.2 to 94.9%, 4.0 to 19.6% and 20.8 to 92.0% in the conversion of ammoniacal nitrogen to nitrate and nitrite concentration and removal of total nitrogen, respectively. The increase of air flow rate and decrease of the initial concentration of ammoniacal nitrogen resulted in higher efficiencies of total nitrogen removal, as well as the conversion of ammoniacal nitrogen to nitrate. During the pre-established intervals of this study, the removal and conversion efficiencies of nitrogen compounds above 85% were achieved in air flow rate variations from 0.375 to 0.725 L min-1 Lreactor-1 and initial concentration of ammoniacal nitrogen from 80 to 200 mg L-1. On denitrification process, we obtained efficiencies from 91.5 to 96.9% on the removal of nitrite/nitrate and from 78.3 to 87.9% on the removal of organic matter.

Portable flow board for storage of fruits and vegetables in mini-chambers with controlled atmosphere

ABSTRACT A portable flow board system was developed in the present study with the aim to facilitate lab-scale experiments of controlled atmosphere (CA) with fruits and vegetables. This sturdy flow board combines ease fabrication, low cost and gas economy. Its functionality is provided by manifolds and gas mixers. Each gaseous component is supplied by a gas cylinder through a differential valve of adjusted pressure control, generally at 6 kPa, and forced through 13 standardized restrictors coupled to each manifold output. Controlled atmospheres are then formed with one, two or three gases in 13 gas mixers affixed to the flow board base, which are further conducted through flexible tubes to storage mini-chambers that can also be used to study metabolic consumption and production of gaseous components. The restrictors used in the flow gaseous components were manufactured from microhematocrit test-type capillary glass tubes following the hot forming method under continuous air flow. The portable flow board showed to be low cost and simple post-harvest equipment that allows preparing controlled atmospheres in open systems with stable composition and flow, in a manner similar to traditional flow boards with control of gas escape by barostats.

Software for design of experiments and response modelling of cake ltration applications

Filtration is a widely used unit operation in chemical engineering. The huge variation in the properties of materials to be ltered makes the study of ltration a challenging task. One of the objectives of this thesis was to show that conventional ltration theories are di cult to use when the system to be modelled contains all of the stages and features that are present in a complete solid/liquid separation process. Furthermore, most of the ltration theories require experimental work to be performed in order to obtain critical parameters required by the theoretical models. Creating a good overall understanding of how the variables a ect the nal product in ltration is somewhat impossible on a purely theoretical basis. The complexity of solid/liquid separation processes require experimental work and when tests are needed, it is advisable to use experimental design techniques so that the goals can be achieved. The statistical design of experiments provides the necessary tools for recognising the e ects of variables. It also helps to perform experimental work more economically. Design of experiments is a prerequisite for creating empirical models that can describe how the measured response is related to the changes in the values of the variable. A software package was developed that provides a ltration practitioner with experimental designs and calculates the parameters for linear regression models, along with the graphical representation of the responses. The developed software consists of two software modules. These modules are LTDoE and LTRead. The LTDoE module is used to create experimental designs for di erent lter types. The lter types considered in the software are automatic vertical pressure lter, double-sided vertical pressure lter, horizontal membrane lter press, vacuum belt lter and ceramic capillary action disc lter. It is also possible to create experimental designs for those cases where the variables are totally user de ned, say for a customized ltration cycle or di erent piece of equipment. The LTRead-module is used to read the experimental data gathered from the experiments, to analyse the data and to create models for each of the measured responses. Introducing the structure of the software more in detail and showing some of the practical applications is the main part of this thesis. This approach to the study of cake ltration processes, as presented in this thesis, has been shown to have good practical value when making ltration tests.

Alpha2-Adrenoceptor-mediated vascular responses induced by dexmedetomidine

Alpha2-Adrenoceptors are cell-surface G protein coupled receptors that mediate many of the effects of the catecholamines noradrenaline and adrenaline. The three human α2-adrenoceptor subtypes are widely expressed in different tissues and organs, and they mediate many different physiological and pharmacological effects in the central and peripheral nervous system and as postsynaptic receptors in target organs. Previous studies have demonstrated that α2-adrenoceptors mediate both vascular constriction and dilatation in humans. Large inter-individual variation has been observed in the vascular responses to α2-adrenoceptor activation in clinical studies. All three receptor subtypes are potential drug targets. It was therefore considered important to further elucidate the details of adrenergic vascular regulation and its genetic variation, since such knowledge may help to improve the development of future cardiovascular drugs and intensive care therapies. Dexmedetomidine is the most selective and potent α2-adrenoceptor agonist currently available for clinical use. When given systemically, dexmedetomidine induces nearly complete sympatholysis already at low concentrations, and postsynaptic effects, such vasoconstriction, can be observed with increasing concentrations. Thus, local infusions of small doses of dexmedetomidine into dorsal hand veins and the application of pharmacological sympathectomy with brachial plexus block provide a means to assess drug-induced peripheral vascular responses without interference from systemic pharmacological effects and autonomic nervous system regulation. Dexmedetomidine was observed to have biphasic effects on haemodynamics, with an initial decrease in blood pressure at low concentrations followed by substantial increases in blood pressure and coronary vascular resistance at high concentrations. Plasma concentrations of dexmedetomidine that significantly exceeded the recommended therapeutic level did not reduce myocardial blood flow below the level that is observed with the usual therapeutic concentrations and did not induce any evident myocardial ischaemia in healthy subjects. Further, it was demonstrated that dexmedetomidine also had significant vasodilatory effects through activation of endothelial nitric oxide synthesis, and thus when the endothelial component of the blood vessel response to dexmedetomidine was inhibited, peripheral vasoconstriction was augmented. Hand vein constriction responses to α2-adrenoceptor activation by dexmedetomidine were only weakly associated with the constriction responses to α1-adrenoceptor activation, pointing to independent cellular regulation by these two adrenoceptor classes. Substantial inter-individual variation was noted in the venous constriction elicited by activation of α2-adrenoceptors by dexmedetomidine. In two study populations from two different continents, a single nucleotide polymorphism in the PRKCB gene was found to be associated with the dorsal hand vein constriction response to dexmedetomidine, suggesting that protein kinase C beta may have an important role in the vascular α2-adrenoceptor signalling pathways activated by dexmedetomidine.