Corrosion behaviour of boiler tube materials during combustion of fuels containing Zn and Pb

Många förbränningsanläggningar som bränner utmanande bränslen såsom restfraktioner och avfall råkar ut för problem med ökad korrosion på överhettare och/eller vattenväggar pga. komponenter i bränslena som är korrosiva. För att minimera problemen i avfallseldade pannor hålls ångparametrarna på en relativt låg nivå, vilket drastiskt minskar energiproduktionen. Beläggningarna i avfallseldade pannor består till största delen av element som är förknippade med högtemperaturkorrosion: Cl, S, alkalimetaller, främst K och Na, och tungmetaller som Pb och Zn, och det finns också indikationer av Br-förekomst. Det låga ångtrycket i avfallseldade pannor påverkar också stålrörens temperatur i pannväggarna i eldstaden. I dagens läge hålls temperaturen normalt vid 300-400 °C. Alkalikloridorsakad (KCl, NaCl) högtemperaturkorrosion har inte rapporterats vara relevant vid såpass låga temperaturer, men närvaro av Zn- och Pb-komponenter i beläggningarna har påvisats förorsaka ökad korrosion redan vid 300-400 °C. Vid förbränning kan Zn och Pb reagera med S och Cl och bilda klorider och sulfater i rökgaserna. Dessa tungmetallföreningar är speciellt problematiska pga. de bildar lågsmältande saltblandningar. Dessa lågsmältande gasformiga eller fasta föreningar följer rökgasen och kan sedan fastna eller kondensera på kallare ytor på pannväggar eller överhettare för att sedan bilda aggressiva beläggningar. Tungmetallrika (Pb, Zn) klorider och sulfater ökar risken för korrosion, och effekten förstärks ytterligare vid närvaro av smälta. Motivet med den här studien var att få en bättre insikt i högtemperaturkorrosion förorsakad av Zn och Pb, samt att undersöka och prediktera beteendet och motståndskraften hos några stålkvaliteter som används i överhettare och pannväggar i tungmetallrika förhållanden och höga materialtemperaturer. Omfattande laboratorie-, småskale- och fullskaletest utfördes. Resultaten kan direkt utnyttjas i praktiska applikationer, t.ex. vid materialval, eller vid utveckling av korrosionsmotverkande verktyg för att hitta initierande faktorer och förstå deras effekt på högtemperaturkorrosion.

Analysis and validation of space averaged drag model for numerical simulations of gas-solid flows in fluidized beds

This thesis presents an approach for formulating and validating a space averaged drag model for coarse mesh simulations of gas-solid flows in fluidized beds using the two-fluid model. Proper modeling for fluid dynamics is central in understanding any industrial multiphase flow. The gas-solid flows in fluidized beds are heterogeneous and usually simulated with the Eulerian description of phases. Such a description requires the usage of fine meshes and small time steps for the proper prediction of its hydrodynamics. Such constraint on the mesh and time step size results in a large number of control volumes and long computational times which are unaffordable for simulations of large scale fluidized beds. If proper closure models are not included, coarse mesh simulations for fluidized beds do not give reasonable results. The coarse mesh simulation fails to resolve the mesoscale structures and results in uniform solids concentration profiles. For a circulating fluidized bed riser, such predicted profiles result in a higher drag force between the gas and solid phase and also overestimated solids mass flux at the outlet. Thus, there is a need to formulate the closure correlations which can accurately predict the hydrodynamics using coarse meshes. This thesis uses the space averaging modeling approach in the formulation of closure models for coarse mesh simulations of the gas-solid flow in fluidized beds using Geldart group B particles. In the analysis of formulating the closure correlation for space averaged drag model, the main parameters for the modeling were found to be the averaging size, solid volume fraction, and distance from the wall. The closure model for the gas-solid drag force was formulated and validated for coarse mesh simulations of the riser, which showed the verification of this modeling approach. Coarse mesh simulations using the corrected drag model resulted in lowered values of solids mass flux. Such an approach is a promising tool in the formulation of appropriate closure models which can be used in coarse mesh simulations of large scale fluidized beds.

Limestone Reactions in a Fluidized Bed Heat Exchanger of an Oxy-fuel Circulating Fluidized Bed Boiler

Oxy-fuel combustion in a circulating fluidized bed (CFB) boiler appears to be a promising option for capturing CO2 in power plants. Oxy-fuel combustion is based on burning of fuel in the mixture of oxygen and re-circulated flue gas instead of air. Limestone (CaCO3) is typically used for capturing of SO2 in CFB boilers where limestone calcines to calcium oxide (CaO). Because of high CO2 concentration in oxy-fuel combustion, calcination reaction may be hindered or carbonation, the reverse reaction of calcination, may occur. Carbonation of CaO particles can cause problems especially in the circulation loop of a CFB boiler where temperature level is lower than in the furnace. The aim of the thesis was to examine carbonation of CaO in a fluidized bed heat exchanger of a CFB boiler featuring oxy-fuel combustion. The calculations and analyzing were based on measurement data from an oxy-fuel pilot plant and on 0-dimensional (0D) gas balance of a fluidized bed heat exchanger. Additionally, the objective was to develop a 1-dimensional (1D) model of a fluidized bed heat exchanger by searching a suitable pre-exponential factor for a carbonation rate constant. On the basis of gas measurement data and the 0D gas balance, it was found that the amount of fluidization gas decreased as it flew through the fluidized bed heat exchanger. Most likely the reason for this was carbonation of CaO. It was discovered that temperature has a promoting effect on the reaction rate of carbonation. With the 1D model, a suitable pre-exponential factor for the equation of carbonation rate constant was found. However, during measurements there were several uncertainties, and in the calculations plenty of assumptions were made. Besides, the temperature level in the fluidized bed heat exchanger was relatively low during the measurements. Carbonation should be considered when fluidized bed heat exchangers and the capacity of related fans are designed for a CFB boiler with oxy-fuel combustion.

Mathematical modeling of laboratory scale three-phase fixed bed reactors

Detta arbete fokuserar på modellering av katalytiska gas-vätskereaktioner som genomförs i kontinuerliga packade bäddar. Katalyserade gas-vätskereaktioner hör till de mest typiska reaktionerna i kemisk industri; därför behandlas här packade bäddreaktorer som ett av de populäraste alternativen, då kontinuerlig drift eftersträvas. Tack vare en stor katalysatormängd per volym har de en kompakt struktur, separering av katalysatorn behövs inte och genom en professionell design kan den mest fördelaktiga strömningsbilden upprätthållas i reaktorn. Packade bäddreaktorer är attraktiva p.g.a. lägre investerings- och driftskostnader. Även om packade bäddar används intensivt i industri, är det mycket utmanande att modellera. Detta beror på att tre faser samexisterar och systemets geometri är komplicerad. Existensen av flera reaktioner gör den matematiska modelleringen även mera krävande. Många förenklingar blir därmed nödvändiga. Modellerna involverar typiskt flera parametrar som skall justeras på basis av experimentella data. I detta arbete studerades fem olika reaktionssystem. Systemen hade studerats experimentellt i vårt laboratorium med målet att nå en hög produktivitet och selektivitet genom ett optimalt val av katalysatorer och driftsbetingelser. Hydrering av citral, dekarboxylering av fettsyror, direkt syntes av väteperoxid samt hydrering av sockermonomererna glukos och arabinos användes som exempelsystem. Även om dessa system hade mycket gemensamt, hade de också unika egenskaper och krävde därför en skräddarsydd matematisk behandling. Citralhydrering var ett system med en dominerande huvudreaktion som producerar citronellal och citronellol som huvudprodukter. Produkterna används som en citrondoftande komponent i parfymer, tvålar och tvättmedel samt som plattform-kemikalier. Dekarboxylering av stearinsyra var ett specialfall, för vilket en reaktionsväg för produktion av långkedjade kolväten utgående från fettsyror söktes. En synnerligen hög produktselektivitet var karakteristisk för detta system. Även processuppskalning modellerades för dekarboxylerings-reaktionen. Direkt syntes av väteperoxid hade som målsättning att framta en förenklad process att producera väteperoxid genom att låta upplöst väte och syre reagera direkt i ett lämpligt lösningsmedel på en aktiv fast katalysator. I detta system förekommer tre bireaktioner, vilka ger vatten som oönskad produkt. Alla dessa tre reaktioner modellerades matematiskt med hjälp av dynamiska massbalanser. Målet med hydrering av glukos och arabinos är att framställa produkter med en hög förädlingsgrad, nämligen sockeralkoholer, genom katalytisk hydrering. För dessa två system löstes ämnesmängd- och energibalanserna simultant för att evaluera effekter inne i porösa katalysatorpartiklar. Impulsbalanser som bestämmer strömningsbetingelser inne i en kemisk reaktor, ersattes i alla modelleringsstudier med semi-empiriska korrelationsuttryck för vätskans volymandel och tryckförlust och med axiell dispersionsmodell för beskrivning av omblandningseffekter. Genom att justera modellens parametrar kunde reaktorns beteende beskrivas väl. Alla experiment var genomförda i laboratorieskala. En stor mängd av kopplade effekter samexisterade: reaktionskinetik inklusive adsorption, katalysatordeaktivering, mass- och värmeöverföring samt strömningsrelaterade effekter. En del av dessa effekter kunde studeras separat (t.ex. dispersionseffekter och bireaktioner). Inverkan av vissa fenomen kunde ibland minimeras genom en noggrann planering av experimenten. På detta sätt kunde förenklingar i modellerna bättre motiveras. Alla system som studerades var industriellt relevanta. Utveckling av nya, förenklade produktionsteknologier för existerande kemiska komponenter eller nya komponenter är ett gigantiskt uppdrag. Studierna som presenterades här fokuserade på en av den teknisk-vetenskapliga utfärdens första etapper.

Fast pyrolysis in circulating fluidized bed: feasibility study

Today the limitedness of fossil fuel resources is clearly realized. For this reason there is a strong focus throughout the world on shifting from fossil fuel based energy system to biofuel based energy system. In this respect Finland with its proven excellent forestry capabilities has a great potential to accomplish this goal. It is regarded that one of the most efficient ways of wood biomass utilization is to use it as a feedstock for fast pyrolysis process. By means of this process solid biomass is converted into liquid fuel called bio-oil which can be burnt at power plants, used for hydrogen generation through a catalytic steam reforming process and as a source of valuable chemical compounds. Nowadays different configurations of this process have found their applications in several pilot plants worldwide. However the circulating fluidized bed configuration is regarded as the one with the highest potential to be commercialized. In the current Master’s Thesis a feasibility study of circulating fluidized bed fast pyrolysis process utilizing Scots pine logs as a raw material was conducted. The production capacity of the process is 100 000 tonne/year of bio-oil. The feasibility study is divided into two phases: a process design phase and economic feasibility analysis phase. The process design phase consists of mass and heat balance calculations, equipment sizing, estimation of pressure drops in the pipelines and development of plant layout. This phase resulted in creation of process flow diagrams, equipment list and Microsoft Excel spreadsheet that calculates the process mass and heat balances depending on the bio-oil production capacity which can be set by a user. These documents are presented in the current report as appendices. In the economic feasibility analysis phase there were at first calculated investment and operating costs of the process. Then using these costs there was calculated the price of bio-oil which is required to reach the values of internal rate of return of 5%, 10%, 20%, 30%, 40%, and 50%.

Modeling of temperature regimes in a furnace of bubbling fluidized bed boiler

The work is mainly focused on the technology of bubbling fluidized bed combustion. Heat transfer and hydrodynamics of the process were examined in the work in detail. Special emphasis was placed on the process of heat exchange in a freeboard zone of bubbling fluidized bed boiler. Operating mode of bubbling fluidized bed boiler depends on many parameters. To assess the influence of some parameters on a temperature regime inside the furnace a simplified method of zonal modeling was used in the work. Thus, effects of bed material fineness, excess air ratio and changes in boiler load were studied. Besides the technology of combustion in bubbling fluidized bed, other common technologies of solid fuels combustion were reviewed. In addition, brief survey of most widely used types of solid fuel was performed in the work.

Bio-ethanol valorization towards C4s including 1-butanol over metal modified alumina and zeolite catalysts

Bio-ethanol has been used as a fuel additive in modern society aimed at reducing CO2-emissions and dependence on oil. However, ethanol is unsuitable as fuel supplement in higher proportions due to its physico-chemical properties. One option to counteract the negative effects is to upgrade ethanol in a continuous fixed bed reactor to more valuable C4 products such as 1-butanol providing chemical similarity with traditional gasoline components. Bio-ethanol based valorization products also have other end-uses than just fuel additives. E.g. 1-butanol and ethyl acetate are well characterised industrial solvents and platform chemicals providing greener alternatives. The modern approach is to apply heterogeneous catalysts in the investigated reactions. The research was concentrated on aluminium oxide (Al2O3) and zeolites that were used as catalysts and catalyst supports. The metals supported (Cu, Ni, Co) gave very different product profiles and, thus, a profound view of different catalyst preparation methods and characterisation techniques was necessary. Additionally, acidity and basicity of the catalyst surface have an important role in determining the product profile. It was observed that ordinary determination of acid strength was not enough to explain all the phenomena e.g. the reaction mechanism. One of the main findings of the thesis is based on the catalytically active site which originates from crystallite structure. As a consequence, the overall evaluation of different by-products and intermediates was carried out by combining the information. Further kinetic analysis was carried out on metal (Cu, Ni, Co) supported self-prepared alumina catalysts. The thesis gives information for further catalyst developments aimed to scale-up towards industrially feasible operations.

Particle model for simulating limestone reactions in novel fluidised bed energy applications

The development of carbon capture and storage (CCS) has raised interest towards novel fluidised bed (FB) energy applications. In these applications, limestone can be utilized for S02 and/or CO2 capture. The conditions in the new applications differ from the traditional atmospheric and pressurised circulating fluidised bed (CFB) combustion conditions in which the limestone is successfully used for SO2 capture. In this work, a detailed physical single particle model with a description of the mass and energy transfer inside the particle for limestone was developed. The novelty of this model was to take into account the simultaneous reactions, changing conditions, and the effect of advection. Especially, the capability to study the cyclic behaviour of limestone on both sides of the calcination-carbonation equilibrium curve is important in the novel conditions. The significances of including advection or assuming diffusion control were studied in calcination. Especially, the effect of advection in calcination reaction in the novel combustion atmosphere was shown. The model was tested against experimental data; sulphur capture was studied in a laboratory reactor in different fluidised bed conditions. Different Conversion levels and sulphation patterns were examined in different atmospheres for one limestone type. The Conversion curves were well predicted with the model, and the mechanisms leading to the Conversion patterns were explained with the model simulations. In this work, it was also evaluated whether the transient environment has an effect on the limestone behaviour compared to the averaged conditions and in which conditions the effect is the largest. The difference between the averaged and transient conditions was notable only in the conditions which were close to the calcination-carbonation equilibrium curve. The results of this study suggest that the development of a simplified particle model requires a proper understanding of physical and chemical processes taking place in the particle during the reactions. The results of the study will be required when analysing complex limestone reaction phenomena or when developing the description of limestone behaviour in comprehensive 3D process models. In order to transfer the experimental observations to furnace conditions, the relevant mechanisms that take place need to be understood before the important ones can be selected for 3D process model. This study revealed the sulphur capture behaviour under transient oxy-fuel conditions, which is important when the oxy-fuel CFB process and process model are developed.

Development of remote monitoring process tools for biomass fired bubbling fluidized bed boilers

Remote monitoring of a power boiler allows the supplying company to make sure that equipment is used as supposed to and gives a good chance for process optimization. This improves co-operation between the supplier and the customer and creates an aura of trust that helps securing future contracts. Remote monitoring is already in use with recovery boilers but the goal is to expand especially to biomass-fired BFB-boilers. To make remote monitoring possible, data has to be measured reliably on site and the link between the power plant and supplying company’s server has to work reliably. Data can be gathered either with the supplier’s sensors or with measurements originally installed in the power plant if the plant in question is not originally built by the supplying company. Main goal in remote monitoring is process optimization and avoiding unnecessary accidents. This can be achieved for instance by following the efficiency curves and fouling in different parts of the process and comparing them to past values. The final amount of calculations depends on the amount of data gathered. Sudden changes in efficiency or fouling require further notice and in such a case it’s important that dialogue toward the power plant in question also works.

Treatment of Fe-containing Acid Mine Waters in Fixed Bed Adsorption Column with Calcite Side-Stone

The acid mining drainage is considered the most significant environmental pollution problem around the world for the extensive formation acidic leachates containing heavy metals. Adsorption is widely used methods in water treatment due to it easy operation and the availability of a wide variety of commercial adsorbent (low cost). The primary goal of this thesis was to investigate the efficiency of neutralizing agents, CaCO3 and CaSiO3, and metal adsorption materials with unmodified limestone from Company Nordkalk Oy. In addition to this, the side materials of limestone mining were tested for iron adsorption from acidic model solution. This study was executed at Lappeenranta University of Technology, Finland. The work utilised fixed-bed adsorption column as the main equipment and large fluidized column. Atomic absorption spectroscopy (AAS) and x-ray diffraction (XRD) was used to determine ferric removal and the composition of material respectively. The results suggest a high potential for the studied materials to be used a low cost adsorbents in acid mine drainage treatment. From the two studied adsorbents, the FS material was more suitable than the Gotland material. Based on the findings, it is recommended that further studies might include detailed analysis of Gotland materials.

Methyl and ethyl chloride synthesis in microreactors

Methyl chloride is an important chemical intermediate with a variety of applications. It is produced today in large units and shipped to the endusers. Most of the derived products are harmless, as silicones, butyl rubber and methyl cellulose. However, methyl chloride is highly toxic and flammable. On-site production in the required quantities is desirable to reduce the risks involved in transportation and storage. Ethyl chloride is a smaller-scale chemical intermediate that is mainly used in the production of cellulose derivatives. Thus, the combination of onsite production of methyl and ethyl chloride is attractive for the cellulose processing industry, e.g. current and future biorefineries. Both alkyl chlorides can be produced by hydrochlorination of the corresponding alcohol, ethanol or methanol. Microreactors are attractive for the on-site production as the reactions are very fast and involve toxic chemicals. In microreactors, the diffusion limitations can be suppressed and the process safety can be improved. The modular setup of microreactors is flexible to adjust the production capacity as needed. Although methyl and ethyl chloride are important chemical intermediates, the literature available on potential catalysts and reaction kinetics is limited. Thus the thesis includes an extensive catalyst screening and characterization, along with kinetic studies and engineering the hydrochlorination process in microreactors. A range of zeolite and alumina based catalysts, neat and impregnated with ZnCl2, were screened for the methanol hydrochlorination. The influence of zinc loading, support, zinc precursor and pH was investigated. The catalysts were characterized with FTIR, TEM, XPS, nitrogen physisorption, XRD and EDX to identify the relationship between the catalyst characteristics and the activity and selectivity in the methyl chloride synthesis. The acidic properties of the catalyst were strongly influenced upon the ZnCl2 modification. In both cases, alumina and zeolite supports, zinc reacted to a certain amount with specific surface sites, which resulted in a decrease of strong and medium Brønsted and Lewis acid sites and the formation of zinc-based weak Lewis acid sites. The latter are highly active and selective in methanol hydrochlorination. Along with the molecular zinc sites, bulk zinc species are present on the support material. Zinc modified zeolite catalysts exhibited the highest activity also at low temperatures (ca 200 °C), however, showing deactivation with time-onstream. Zn/H-ZSM-5 zeolite catalysts had a higher stability than ZnCl2 modified H-Beta and they could be regenerated by burning the coke in air at 400 °C. Neat alumina and zinc modified alumina catalysts were active and selective at 300 °C and higher temperatures. However, zeolite catalysts can be suitable for methyl chloride synthesis at lower temperatures, i.e. 200 °C. Neat γ-alumina was found to be the most stable catalyst when coated in a microreactor channel and it was thus used as the catalyst for systematic kinetic studies in the microreactor. A binder-free and reproducible catalyst coating technique was developed. The uniformity, thickness and stability of the coatings were extensively characterized by SEM, confocal microscopy and EDX analysis. A stable coating could be obtained by thermally pretreating the microreactor platelets and ball milling the alumina to obtain a small particle size. Slurry aging and slow drying improved the coating uniformity. Methyl chloride synthesis from methanol and hydrochloric acid was performed in an alumina-coated microreactor. Conversions from 4% to 83% were achieved in the investigated temperature range of 280-340 °C. This demonstrated that the reaction is fast enough to be successfully performed in a microreactor system. The performance of the microreactor was compared with a tubular fixed bed reactor. The results obtained with both reactors were comparable, but the microreactor allows a rapid catalytic screening with low consumption of chemicals. As a complete conversion of methanol could not be reached in a single microreactor, a second microreactor was coupled in series. A maximum conversion of 97.6 % and a selectivity of 98.8 % were reached at 340°C, which is close to the calculated values at a thermodynamic equilibrium. A kinetic model based on kinetic experiments and thermodynamic calculations was developed. The model was based on a Langmuir Hinshelwood-type mechanism and a plug flow model for the microreactor. The influence of the reactant adsorption on the catalyst surface was investigated by performing transient experiments and comparing different kinetic models. The obtained activation energy for methyl chloride was ca. two fold higher than the previously published, indicating diffusion limitations in the previous studies. A detailed modeling of the diffusion in the porous catalyst layer revealed that severe diffusion limitations occur starting from catalyst coating thicknesses of 50 μm. At a catalyst coating thickness of ca 15 μm as in the microreactor, the conditions of intrinsic kinetics prevail. Ethanol hydrochlorination was performed successfully in the microreactor system. The reaction temperature was 240-340°C. An almost complete conversion of ethanol was achieved at 340°C. The product distribution was broader than for methanol hydrochlorination. Ethylene, diethyl ether and acetaldehyde were detected as by-products, ethylene being the most dominant by-product. A kinetic model including a thorough thermodynamic analysis was developed and the influence of adsorbed HCl on the reaction rate of ethanol dehydration reactions was demonstrated. The separation of methyl chloride using condensers was investigated. The proposed microreactor-condenser concept enables the production of methyl chloride with a high purity of 99%.

Deoxygenation of fatty acids for production of fuels and chemicals

The decreasing fossil fuel resources combined with an increasing world energy demand has raised an interest in renewable energy sources. The alternatives can be solar, wind and geothermal energies, but only biomass can be a substitute for the carbon–based feedstock, which is suitable for the production of transportation fuels and chemicals. However, a high oxygen content of the biomass creates challenges for the future chemical industry, forcing the development of new processes which allow a complete or selective oxygen removal without any significant carbon loss. Therefore, understanding and optimization of biomass deoxygenation processes are crucial for the future bio–based chemical industry. In this work, deoxygenation of fatty acids and their derivatives was studied over Pd/C and TiO2 supported noble metal catalysts (Pt, Pt–Re, Re and Ru) to obtain future fuel components. The 5 % Pd/C catalyst was investigated in semibatch and fixed bed reactors at 300 °C and 1.7–2 MPa of inert and hydrogen–containing atmospheres. Based on extensive kinetic studies, plausible reaction mechanisms and pathways were proposed. The influence of the unsaturation in the deoxygenation of model compounds and industrial feedstock – tall oil fatty acids – over a Pd/C catalyst was demonstrated. The optimization of the reaction conditions suppressed the formation of by–products, hence high yields and selectivities towards linear hydrocarbons and catalyst stability were achieved. Experiments in a fixed bed reactor filled with a 2 % Pd/C catalyst were performed with stearic acid as a model compound at different hydrogen–containing gas atmospheres to understand the catalyst stability under various conditions. Moreover, prolonged experiments were carried out with concentrated model compounds to reveal the catalyst deactivation. New materials were proposed for the selective deoxygenation process at lower temperatures (~200 °C) with a tunable selectivity to hydrodeoxygenation by using 4 % Pt/TiO2 or decarboxylation/decarbonylation over 4 % Ru/TiO2 catalysts. A new method for selective hydrogenation of fatty acids to fatty alcohols was demonstrated with a 4 % Re/TiO2 catalyst. A reaction pathway and mechanism for TiO2 supported metal catalysts was proposed and an optimization of the process conditions led to an increase in the formation of the desired products.

Application and development of numerical methods for the modelling of innovative gas cooled fission reactors

Innovative gas cooled reactors, such as the pebble bed reactor (PBR) and the gas cooled fast reactor (GFR) offer higher efficiency and new application areas for nuclear energy. Numerical methods were applied and developed to analyse the specific features of these reactor types with fully three dimensional calculation models. In the first part of this thesis, discrete element method (DEM) was used for a physically realistic modelling of the packing of fuel pebbles in PBR geometries and methods were developed for utilising the DEM results in subsequent reactor physics and thermal-hydraulics calculations. In the second part, the flow and heat transfer for a single gas cooled fuel rod of a GFR were investigated with computational fluid dynamics (CFD) methods. An in-house DEM implementation was validated and used for packing simulations, in which the effect of several parameters on the resulting average packing density was investigated. The restitution coefficient was found out to have the most significant effect. The results can be utilised in further work to obtain a pebble bed with a specific packing density. The packing structures of selected pebble beds were also analysed in detail and local variations in the packing density were observed, which should be taken into account especially in the reactor core thermal-hydraulic analyses. Two open source DEM codes were used to produce stochastic pebble bed configurations to add realism and improve the accuracy of criticality calculations performed with the Monte Carlo reactor physics code Serpent. Russian ASTRA criticality experiments were calculated. Pebble beds corresponding to the experimental specifications within measurement uncertainties were produced in DEM simulations and successfully exported into the subsequent reactor physics analysis. With the developed approach, two typical issues in Monte Carlo reactor physics calculations of pebble bed geometries were avoided. A novel method was developed and implemented as a MATLAB code to calculate porosities in the cells of a CFD calculation mesh constructed over a pebble bed obtained from DEM simulations. The code was further developed to distribute power and temperature data accurately between discrete based reactor physics and continuum based thermal-hydraulics models to enable coupled reactor core calculations. The developed method was also found useful for analysing sphere packings in general. CFD calculations were performed to investigate the pressure losses and heat transfer in three dimensional air cooled smooth and rib roughened rod geometries, housed inside a hexagonal flow channel representing a sub-channel of a single fuel rod of a GFR. The CFD geometry represented the test section of the L-STAR experimental facility at Karlsruhe Institute of Technology and the calculation results were compared to the corresponding experimental results. Knowledge was gained of the adequacy of various turbulence models and of the modelling requirements and issues related to the specific application. The obtained pressure loss results were in a relatively good agreement with the experimental data. Heat transfer in the smooth rod geometry was somewhat under predicted, which can partly be explained by unaccounted heat losses and uncertainties. In the rib roughened geometry heat transfer was severely under predicted by the used realisable k − epsilon turbulence model. An additional calculation with a v2 − f turbulence model showed significant improvement in the heat transfer results, which is most likely due to the better performance of the model in separated flow problems. Further investigations are suggested before using CFD to make conclusions of the heat transfer performance of rib roughened GFR fuel rod geometries. It is suggested that the viewpoints of numerical modelling are included in the planning of experiments to ease the challenging model construction and simulations and to avoid introducing additional sources of uncertainties. To facilitate the use of advanced calculation approaches, multi-physical aspects in experiments should also be considered and documented in a reasonable detail.

Simulation of the Interaction Mean Free Path between Neutrons and TRISO Particles

The interaction mean free path between neutrons and TRISO particles is simulated using scripts written in MATLAB to solve the increasing error present with an increase in the packing factor in the reactor physics code Serpent. Their movement is tracked both in an unbounded and in a bounded space. Their track is calculated, depending on the program, linearly directly using the position vectors of the neutrons and the surface equations of all the fuel particles; by dividing the space in multiple subspaces, each of which contain a fraction of the total number of particles, and choosing the particles from those subspaces through which the neutron passes through; or by choosing the particles that lie within an infinite cylinder formed on the movement axis of the neutron. The estimate from the current analytical model, based on an exponential distribution, for the mean free path, utilized by Serpent, is used as a reference result. The results from the implicit model in Serpent imply a too long mean free path with high packing factors. The received results support this observation by producing, with a packing factor of 17 %, approximately 2.46 % shorter mean free path compared to the reference model. This is supported by the packing factor experienced by the neutron, the simulation of which resulted in a 17.29 % packing factor. It was also observed that the neutrons leaving from the surfaces of the fuel particles, in contrast to those starting inside the moderator, do not follow the exponential distribution. The current model, as it is, is thus not valid in the determination of the free path lengths of the neutrons.

Analysis and modelling of chemical looping combustion process with and without oxygen uncoupling

Effective control and limiting of carbon dioxide (CO₂) emissions in energy production are major challenges of science today. Current research activities include the development of new low-cost carbon capture technologies, and among the proposed concepts, chemical combustion (CLC) and chemical looping with oxygen uncoupling (CLOU) have attracted significant attention allowing intrinsic separation of pure CO₂ from a hydrocarbon fuel combustion process with a comparatively small energy penalty. Both CLC and CLOU utilize the well-established fluidized bed technology, but several technical challenges need to be overcome in order to commercialize the processes. Therefore, development of proper modelling and simulation tools is essential for the design, optimization, and scale-up of chemical looping-based combustion systems. The main objective of this work was to analyze the technological feasibility of CLC and CLOU processes at different scales using a computational modelling approach. A onedimensional fluidized bed model frame was constructed and applied for simulations of CLC and CLOU systems consisting of interconnected fluidized bed reactors. The model is based on the conservation of mass and energy, and semi-empirical correlations are used to describe the hydrodynamics, chemical reactions, and transfer of heat in the reactors. Another objective was to evaluate the viability of chemical looping-based energy production, and a flow sheet model representing a CLC-integrated steam power plant was developed. The 1D model frame was succesfully validated based on the operation of a 150 kWth laboratory-sized CLC unit fed by methane. By following certain scale-up criteria, a conceptual design for a CLC reactor system at a pre-commercial scale of 100 MWth was created, after which the validated model was used to predict the performance of the system. As a result, further understanding of the parameters affecting the operation of a large-scale CLC process was acquired, which will be useful for the practical design work in the future. The integration of the reactor system and steam turbine cycle for power production was studied resulting in a suggested plant layout including a CLC boiler system, a simple heat recovery setup, and an integrated steam cycle with a three pressure level steam turbine. Possible operational regions of a CLOU reactor system fed by bituminous coal were determined via mass, energy, and exergy balance analysis. Finally, the 1D fluidized bed model was modified suitable for CLOU, and the performance of a hypothetical 500 MWth CLOU fuel reactor was evaluated by extensive case simulations.