CFD study of penetration and mixing of fuel in a circulating fluidized bed furnace

The aim of this thesis is to study the mixing of fuel and, also to some extent, the mixing of air in a circulating fluidized bed boiler. In the literature survey part of this thesis, a review is made of the previous experimental studies related to the fuel and air mixing in the circulating fluidized beds. In the simulation part of it the commercial computational fluid dynamics software (FLUENT) is used with the Eulerian multiphase model for studying the fuel mixing in the two and three-dimensional furnace geometries. The results of the three-dimensional simulations are promising and, therefore suggestions are made for the future simulations. The two-dimensional studies give new information of the effects of the fluidization velocity, fuel particle size and fuel density on the fuel mixing. However, the present results show that three-dimensional models produce more realistic representation of the circulating fluidized bed behavior.

Modelling of ultrasound assisted mixing in Newtonian and non-Newtonian liquids

The purpose of this work is to obtain a better understanding of behaviour of possible ultrasound appliance on fluid media mixing. The research is done in the regard to Newtonian and non-Newtonian fluids. The process of ultrasound appliance on liquids is modelled in COMSOL Multiphysics software. The influence of ultrasound using is introduced as waveform equation. Turbulence modelling is fulfilled by the k-ε model in Newtonian fluid. The modeling of ultrasound assisted mixing in non-Newtonian fluids is based on the power law. To verify modelling results two practical methods are used: Particle Image Velocimetry and measurements of mixing time. Particle Image Velocimetry allows capturing of velocity flow field continuously and presents detailed depiction of liquid dynamics. The second way of verification is the comparison of mixing time of homogeneity. Experimentally achievement of mixing time is done by conductivity measurements. In modelling part mixing time is achieved by special module of COMSOL Multiphysics – the transport of diluted species. Both practical and modelling parts show similar radial mechanism of fluid flow under ultrasound appliance – from the horn tip fluid moves to the bottom and along the walls goes back. Velocity profiles are similar in modelling and experimental part in the case of Newtonian fluid. In the case of non-Newtonian fluid velocity profiles do not agree. The development track of ultrasound-assisted mixing modelling is presented in the thesis.

Phosphorus extractability in surface soil samples as affected by mixing with subsoil

Selostus: Muokkauskerrokseen sekoitetun pohjamaan vaikutus maasta uuttuvan fosforin määrään

Linjasekoittimensekoituskyvyn määrittäminen

Työn tavoitteena oli kehittää sekoituskyvyn empiirinen monimuuttujafunktio eräälle dynaamiselle linjasekoittimelle. Monimuuttujafunktio oli tarkoitus kehittää sekoittimen kierrosnopeuden, virtaavien materiaalien ja tilavuusvirtausten perus-teella. Työn kirjallisuusosassa tarkasteltiin pääosin dynaamisia linjasekoittimia ja niiden toimintaperiaatteita sekä sekoituskyvyn määrittämiseen soveltuvia mittausmene-telmiä. Sekoitettavien aineiden epähomogeenisuus määritettiin mittaamalla termoelemen-teillä sekoittimen jälkeisen virran lämpötilavaihtelut. Sekoitettavina aineina käy-tettiin kolmea erilaista seosta: vesi-vesi-, kuitususpensio-vesi- ja karboksyylime-tyyliselluloosa-vesiseoksia. Muita muuttujia kokeissa olivat neljä eri sekoittimen kierrosnopeutta, neljä eri päävirran tilavuusvirtausta ja kolme päävirran sekäsii-hen lisättävän sivuvirran suhdetta. Monimuuttujafunktiokehitettiin vain kuitususpensio-vesiajojen pohjalta muilla seoksilla tehtyjen koeajojen osittaisen epäonnistumisen vuoksi. Hallitsevaksi pa-rametriksi monimuuttujafunktiossa osoittautui sekoittimen kierrosnopeus. Osittai-nen epäonnistuminen kokeissa johtui osaksi termoelementtien käytön sopimatto-muudesta sekoituskyvyn määrittämiseen ja osaksi sekoitettavien aineiden valin-nasta. Jatkotutkimuksiaajatellen muista sekoituskyvyn mittaamiseen käytettävissä olevista menetelmistäkäyttökelpoisimmaksi arvioitiin sähkötomografia.

Particle transport method for convection-diffusion-reaction problems

Convective transport, both pure and combined with diffusion and reaction, can be observed in a wide range of physical and industrial applications, such as heat and mass transfer, crystal growth or biomechanics. The numerical approximation of this class of problemscan present substantial difficulties clue to regions of high gradients (steep fronts) of the solution, where generation of spurious oscillations or smearing should be precluded. This work is devoted to the development of an efficient numerical technique to deal with pure linear convection and convection-dominated problems in the frame-work of convection-diffusion-reaction systems. The particle transport method, developed in this study, is based on using rneshless numerical particles which carry out the solution along the characteristics defining the convective transport. The resolution of steep fronts of the solution is controlled by a special spacial adaptivity procedure. The serni-Lagrangian particle transport method uses an Eulerian fixed grid to represent the solution. In the case of convection-diffusion-reaction problems, the method is combined with diffusion and reaction solvers within an operator splitting approach. To transfer the solution from the particle set onto the grid, a fast monotone projection technique is designed. Our numerical results confirm that the method has a spacial accuracy of the second order and can be faster than typical grid-based methods of the same order; for pure linear convection problems the method demonstrates optimal linear complexity. The method works on structured and unstructured meshes, demonstrating a high-resolution property in the regions of steep fronts of the solution. Moreover, the particle transport method can be successfully used for the numerical simulation of the real-life problems in, for example, chemical engineering.

Ash Particle Erosion on Steam Boiler Convective Section

In this study, equations for the calculation of erosion wear caused by ash particles on convective heat exchanger tubes of steam boilers are presented. Anew, three-dimensional test arrangement was used in the testing of the erosion wear of convective heat exchanger tubes of steam boilers. When using the sleeve-method, three different tube materials and three tube constructions could be tested. New results were obtained from the analyses. The main mechanisms of erosionwear phenomena and erosion wear as a function of collision conditions and material properties have been studied. Properties of fossil fuels have also been presented. When burning solid fuels, such as pulverized coal and peat in steam boilers, most of the ash is entrained by the flue gas in the furnace. In bubbling andcirculating fluidized bed boilers, particle concentration in the flue gas is high because of bed material entrained in the flue gas. Hard particles, such as sharp edged quartz crystals, cause erosion wear when colliding on convective heat exchanger tubes and on the rear wall of the steam boiler. The most important ways to reduce erosion wear in steam boilers is to keep the velocity of the flue gas moderate and prevent channelling of the ash flow in a certain part of the cross section of the flue gas channel, especially near the back wall. One can do this by constructing the boiler with the following components. Screen plates can beused to make the velocity and ash flow distributions more even at the cross-section of the channel. Shield plates and plate type constructions in superheaters can also be used. Erosion testing was conducted with three types of tube constructions: a one tube row, an inline tube bank with six tube rows, and a staggered tube bank with six tube rows. Three flow velocities and two particle concentrations were used in the tests, which were carried out at room temperature. Three particle materials were used: quartz, coal ash and peat ash particles. Mass loss, diameter loss and wall thickness loss measurements of the test sleeves were taken. Erosion wear as a function of flow conditions, tube material and tube construction was analyzed by single-variable linear regression analysis. In developing the erosion wear calculation equations, multi-variable linear regression analysis was used. In the staggered tube bank, erosion wear had a maximum value in a tube row 2 and a local maximum in row 5. In rows 3, 4 and 6, the erosion rate was low. On the other hand, in the in-line tube bank the minimum erosion rate occurred in tube row 2 and in further rows the erosion had an increasing value, so that in a six row tube bank, the maximum value occurred in row 6.

Fabrication and characterization of silicon position sensitive particle detectors

Position sensitive particle detectors are needed in high energy physics research. This thesis describes the development of fabrication processes and characterization techniques of silicon microstrip detectors used in the work for searching elementary particles in the European center for nuclear research, CERN. The detectors give an electrical signal along the particles trajectory after a collision in the particle accelerator. The trajectories give information about the nature of the particle in the struggle to reveal the structure of the matter and the universe. Detectors made of semiconductors have a better position resolution than conventional wire chamber detectors. Silicon semiconductor is overwhelmingly used as a detector material because of its cheapness and standard usage in integrated circuit industry. After a short spread sheet analysis of the basic building block of radiation detectors, the pn junction, the operation of a silicon radiation detector is discussed in general. The microstrip detector is then introduced and the detailed structure of a double-sided ac-coupled strip detector revealed. The fabrication aspects of strip detectors are discussedstarting from the process development and general principles ending up to the description of the double-sided ac-coupled strip detector process. Recombination and generation lifetime measurements in radiation detectors are discussed shortly. The results of electrical tests, ie. measuring the leakage currents and bias resistors, are displayed. The beam test setups and the results, the signal to noise ratio and the position accuracy, are then described. It was found out in earlier research that a heavy irradiation changes the properties of radiation detectors dramatically. A scanning electron microscope method was developed to measure the electric potential and field inside irradiated detectorsto see how a high radiation fluence changes them. The method and the most important results are discussed shortly.

Propagation of the Ignition Front against Airflow in Packed Beds of Wood Particles

Combustion of wood is increasing because of the needs of decreasing the emissions of carbon dioxide and the amount of waste going to landfills. Wood based fuels are often scattered on a large area. The transport distances should be short enough to prevent too high costs, and so the size of heating and power plants using wood fuels is often rather small. Combustion technologies of small-size units have to be developed to reach efficient and environmentally friendly energy production. Furnaces that use different packed bed combustion or gasification techniques areoften most economic in small-scale energy production. Ignition front propagation rate affects the stability, heat release rate and emissions of packed bed combustion. Ignition front propagation against airflow in packed beds of wood fuels has been studied. The research has been carried out mainly experimentally. Theoretical aspects have been considered to draw conclusions about the experimental results. The effects of airflow rate, moisture content of the fuel, size, shape and density of particles, and porosity of the bed on the propagation rate of the ignition front have been studied. The experiments were carried out in a pot furnace. The fuels used in the experiments were mainly real wood fuels that are often burned in the production of energy. The fuel types were thin wood chips, saw dust, shavings, wood chips, and pellets with different sizes. Also a few mixturesof the above were tested. Increase in the moisture content of the fuel decreases the propagation rates of the ignition front and makes the range of possible airflow rates narrower because of the energy needed for the evaporation of water and the dilution of volatile gases due to evaporated steam. Increase in the airflow rate increases the ignition rate until a maximum rate of propagation is reached after which it decreases. The maximum flame propagation rate is not always reached in stoichiometric combustion conditions. Increase in particle size and density transfers the optimum airflow rate towards fuel lean conditions. Mixing of small and large particles is often advantageous, because small particles make itpossible to reach the maximum ignition rate in fuel rich conditions, and large particles widen the range of possible airflow rates. A correlation was found forthe maximum rate of ignition front propagation in different wood fuels. According to the correlation, the maximum ignition mass flux is increased when the sphericity of the particles and the porosity of the bed are increased and the moisture content of the fuel is decreased. Another fit was found between sphericity and porosity. Increase in sphericity decreases the porosity of the bed. The reasons of the observed results are discussed.

ADVANCED NUMERICAL METHODS FOR SIMULATING MICROBUBBLE IN AN AERATED MIXING TANK

Fluid mixing in mechanically agitated tanks is one of the major unit operations in many industries. Bubbly flows have been of interest among researchers in physics, medicine, chemistry and technology over the centuries. The aim of this thesis is to use advanced numerical methods for simulating microbubble in an aerated mixing tank. Main components of the mixing tank are a cylindrical vessel, a rotating Rushton turbine and the air nozzle. The objective of Computational Fluid Dynamics (CFD) is to predict fluid flow, heat transfer, mass transfer and chemical reactions. The CFD simulations of a turbulent bubbly flow are carried out in a cylindrical mixing tank using large eddy simulation (LES) and volume of fluid (VOF) method. The Rushton turbine induced flow is modeled by using a sliding mesh method. Numerical results are used to describe the bubbly flows in highly complex liquid flow. Some of the experimental works related to turbulent bubbly flow in a mixing tank are briefly reported. Numerical simulations are needed to complete and interpret the results of the experimental work. Information given by numerical simulations has a major role in designing and scaling-up mixing tanks. The results of this work have been reported in the following scientific articles: ·Honkanen M., Koohestany A., Hatunen T., Saarenrinne P., Zamankhan P., Large eddy simulations and PIV experiments of a two-phase air-water mixer, in Proceedings of ASME Fluids Engineering Summer Conference (2005). ·Honkanen M., Koohestany A., Hatunen T., Saarenrinne P., Zamankhan P., Dynamical States of Bubbling in an Aerated Stirring Tank, submitted to J. Computational Physics.

Particle Size Analysis of Pharmaceutical Powders by Laser Diffraction

Raaka-aineen hiukkaskoko on lääkekehityksessä keskeinen materiaaliparametri. Lääkeaineen partikkelikoko vaikuttaa moneen lääketuotteen tärkeään ominaisuuteen, esimerkiksi lääkkeen biologiseen hyväksikäytettävyyteen. Tässä diplomityössä keskityttiin jauhemaisten lääkeaineiden hiukkaskoon määrittämiseen laserdiffraktiomenetelmällä. Menetelmä perustuu siihen, että partikkeleista sironneen valon intensiteetin sirontakulmajakauma on riippuvainen partikkelien kokojakaumasta. Työn kirjallisuusosassa esiteltiin laserdiffraktiomenetelmän teoriaa. PIDS (Polarization Intensity Differential Scattering) tekniikka, jota voidaan käyttää laserdiffraktion yhteydessä, on myös kuvattu kirjallisuusosassa. Muihin menetelmiin perustuvista analyysimenetelmistä tutustuttiin mikroskopiaan sekä aerodynaamisen lentoajan määrittämiseen perustuvaan menetelmään. Kirjallisuusosassa esiteltiin myös partikkelikoon yleisimpiä esitystapoja. Työn kokeellisen osan tarkoituksena oli kehittää ja validoida laserdiffraktioon perustuva partikkelikoon määritysmenetelmä tietylle lääkeaineelle. Menetelmäkehitys tehtiin käyttäen Beckman Coulter LS 13 320 laserdiffraktoria. Laite mahdollistaa PIDS-tekniikan käytön laserdiffraktiotekniikan ohella. Menetelmäkehitys aloitettiin arvioimalla, että kyseinen lääkeaine soveltuu parhaiten määritettäväksi nesteeseen dispergoituna. Liukoisuuden perusteella väliaineeksi valittiin tällä lääkeaineella kyllästetty vesiliuos. Dispergointiaineen sekä ultraäänihauteen käyttö havaittiin tarpeelliseksi dispergoidessa kyseistä lääkeainetta kylläiseen vesiliuokseen. Lopuksi sekoitusnopeus näytteensyöttöyksikössä säädettiin sopivaksi. Validointivaiheessa kehitetyn menetelmän todettiin soveltuvan hyvin kyseiselle lääkeaineelle ja tulosten todettiin olevan oikeellisia sekä toistettavia. Menetelmä ei myöskään ollut herkkä pienille häiriöille.

Flokkulaation seuraaminen partikkelikokoanalysaattorilla

Tässä työssä tutkittiin kahden erilaisen partikkelikokoanalysaattorin, PSyA:n ja PIA:n soveltuvuutta flokkuloinnin online-seurantaan. Kummallekin menetelmälle määritettiin raja-arvot, kuten lietteen maksimisakeus. Lisäksi tutkittiin flokkulanttiannostuksen, sekoitusnopeuden, sekoitusajan ja lietteen kiintoainepitoisuuden vaikutusta flokkikokojakaumaan. Kirjallisuusosassa tarkasteltiin kolloidisen suspension ominaispiirteitä, koaguloinnin ja flokkuloinnin teoriaa, flokkulaation kokeellista tutkimista sekä prosessin jatkuvatoimiseen seurantaan soveltuvia laitteita. Lisäksi esitettiin taustaa hydrometallurgisesta prosessista, johon työ liittyy. Flokkauskokeissa käytettiin jätevettä, jonka koostumus vastasi metalliteollisuuden peittausjätevesien tyypillistä koostumusta. Tutkittava jätevesimäärä käsiteltiin ensin kalkkimaidolla, jonka jälkeen saostunut kiintoaine flokattiin synteettisellä polymeeriflokkulantilla. Lietteen keskimääräinen kiintoainepitoisuus oli n. 10 g/l. Esikokeiden perusteella PSyA:lla voitiin mitata ilman laimennusta, mutta PIA:lla tuloksia ei saatu ilman laimentamista kiintoainepitoisuuteen n. 2,5 g/l. Kokeiden aikana havaittiin, että flokit muodostuivat erittäin nopeasti. Flokkien hajoaminen alkoi pian sen jälkeen, kun flokkulantin annostelu lopetettiin. Sekoitusnopeudella 40 r/min tai alle flokit alkoivat laskeutua astian pohjalle sekoituksesta huolimatta ja ne pysyivät pitempään koossa kuin suuremmilla sekoitusnopeuksilla. 5 - 10 minuutin kuluttua flokkulantin lisäämisestä saavutettiin tasapaino, jolloin flokkien kokojakauma ei enää muuttunut. Sekoitusnopeuksilla 80 r/min ja 120 r/min tasapainotilanteen koko-jakauma oli selvästi kapeampi kuin pienimmällä sekoitusnopeudella. Alkuperäisessä lietteessä flokit olivat suurempia kuin laimennetussa lietteessä. PSyA:lla jännepituusjakaumien määrittäminen oli varsin hidasta prosessissa tapahtuviin muutoksiin verrattuna, ja tuloksissa oli suurta hajontaa. PIA:lla saadut partikkelikokojakaumat sitä vastoin olivat johdonmukaisempia, vaikka suurimpien flokkien määrittäminen osoittautuikin epämääräiseksi. Menetelmän suurimmaksi puutteeksi todettiin soveltumattomuus sakeiden lietteiden analysointiin. Kumpikaan menetelmä ei ilman modifiointia sovellu tutkitun lietteen kaltaisten prosessilietteiden flokkuloinnin seurantaan.

Effects of vacancy-type defects in silicon based particle detectors

The semiconductor particle detectors used at CERN experiments are exposed to radiation. Under radiation, the formation of lattice defects is unavoidable. The defects affect the depletion voltage and leakage current of the detectors, and hence affect on the signal-to-noise ratio of the detectors. This shortens the operational lifetime of the detectors. For this reason, the understanding of the formation and the effects of radiation induced defects is crucial for the development of radiation hard detectors. In this work, I have studied the effects of radiation induced defects-mostly vacancy related defects-with a simulation package, Silvaco. Thus, this work essentially concerns the effects of radiation induced defects, and native defects, on leakage currents in particle detectors. Impurity donor atom-vacancy complexes have been proved to cause insignificant increase of leakage current compared with the trivacancy and divacancy-oxygen centres. Native defects and divacancies have proven to cause some of the leakage current, which is relatively small compared with trivacancy and divacancy-oxygen.