APROS code benchmarking inboiling water reactor accident analyses

Diplomityö käsittelee kiehutusvesilaitosten transienttien ja onnettomuuksien analysointia APROS-ohjelmiston avulla. Työ on tehty Teollisuuden Voima Oy:n (TVO) Olkiluoto 1 ja 2 laitosyksiköiden mallin pohjalta. Raportissa esitetään ohjelmiston käyttämiä yhtälöitäja laskentamalleja yleisellä tasolla. Työssä esitellään laitoksen yleispiirteet turvallisuustoimintoineen ja kuvataan ohjelmaan suureksi osaksi aiemmin luotua laskentamallia. Työssä on luetteloitu voimassa olevatlisensiointianalyysit, joiden joukosta on valittu laskentatapauksia ohjelmiston suorituskyvyn arviointia varten. Lisäksi työhön on valittu laskentatapauksia muilla kuin lisensointiin käytetyillä ohjelmilla lasketuista analyyseistä. Lisäksi on suoritettu vertailulaskuja konservatiivisen ja realistisen mallin erojen esille saamiseksi. Laskentatapauksia ovat mm. ylipainetransientti, jäähdytteen menetysonnettomuus ja oletettavissa oleva käyttöhäiriö, jossa pikasulku ei toimi (ATWS). Diplomityön edetessä laitosmallia on kehitetty edelleen lisäämällä joitakin järjestelmiä ja tarkentamalla joidenkin komponenttien kuvausta. Työssä ilmeni, että APROS soveltuu jäähdytteenmenetysonnettomuuden ja suojarakennuksen yhtäaikaiseen analyysiin. APROS.n vaste nopeisiin transientteihin jäi kuitenkin vertailutasosta. Tämän työn perusteella APROS-mallia kehitys jatkuu edelleen siten, että se soveltuisi entistä paremmin myös nopeiden transienttien ja ATWS-tilanteiden kuvaamiseen. Työssä olevaa lisensointianalyysien kuvausta tullaan käyttämään hyväksi selvitettäessä laitoksen turvallisuuden väliarviossa tarvittavien analyysien määrää ja laatua. Nyt saatuja kokemuksia voidaan hyödyntää myös mahdollisen kolmiulotteisen sydänmallin hankinnassa APROS-ohjelmistoon. Tässä diplomityössä esitettyjä parannuksia voidaan käyttää hyväksi SAFIRtutkimusohjelman hankkeiden suunnittelussa.

Gas-liquid mass transfer in bubbly flow: Estimation of mass transfer, bubble size and reactor performance in various applications

Gas-liquid mass transfer is an important issue in the design and operation of many chemical unit operations. Despite its importance, the evaluation of gas-liquid mass transfer is not straightforward due to the complex nature of the phenomena involved. In this thesis gas-liquid mass transfer was evaluated in three different gas-liquid reactors in a traditional way by measuring the volumetric mass transfer coefficient (kLa). The studied reactors were a bubble column with a T-junction two-phase nozzle for gas dispersion, an industrial scale bubble column reactor for the oxidation of tetrahydroanthrahydroquinone and a concurrent downflow structured bed.The main drawback of this approach is that the obtained correlations give only the average volumetric mass transfer coefficient, which is dependent on average conditions. Moreover, the obtained correlations are valid only for the studied geometry and for the chemical system used in the measurements. In principle, a more fundamental approach is to estimate the interfacial area available for mass transfer from bubble size distributions obtained by solution of population balance equations. This approach has been used in this thesis by developing a population balance model for a bubble column together with phenomenological models for bubble breakage and coalescence. The parameters of the bubble breakage rate and coalescence rate models were estimated by comparing the measured and calculated bubble sizes. The coalescence models always have at least one experimental parameter. This is because the bubble coalescence depends on liquid composition in a way which is difficult to evaluate using known physical properties. The coalescence properties of some model solutions were evaluated by measuring the time that a bubble rests at the free liquid-gas interface before coalescing (the so-calledpersistence time or rest time). The measured persistence times range from 10 msup to 15 s depending on the solution. The coalescence was never found to be instantaneous. The bubble oscillates up and down at the interface at least a coupleof times before coalescence takes place. The measured persistence times were compared to coalescence times obtained by parameter fitting using measured bubble size distributions in a bubble column and a bubble column population balance model. For short persistence times, the persistence and coalescence times are in good agreement. For longer persistence times, however, the persistence times are at least an order of magnitude longer than the corresponding coalescence times from parameter fitting. This discrepancy may be attributed to the uncertainties concerning the estimation of energy dissipation rates, collision rates and mechanisms and contact times of the bubbles.

Hydrodynamics and Heat Transfer in an Airlift Reactor

Työn teoriaosassa esitetään kirjallisuudessa esiintyviä teoreettisia ja kokeellisia yhtälöitä nesteen nopeuden, kaasun tilavuusosuuden, painehäviön ja lämmönsiirron laskemiseksi. Lisäksi käsitellään airlift-reaktoreiden toimintaa, rakennetta ja teollisia sovelluksia, sekä sekoitusta ja geometrian vaikutusta airlift-reaktoreiden hydrodynaamisiin ominaisuuksiin. Kokeellisessa osassa kuvataan käytetty koelaitteisto ja mittausmenetelmät sekä esitetään saadut koetulokset. Koelaitteisto on viidellä nousuputkella varustettu ulkoisen kierron airlift-reaktori. Kokeellisessa osassa pyritään ratkaisemaan tällaisessa reaktorissa mahdollisesti esiintyviä ongelmia, kuten "slug flown" muodostuminen nousuputkissa sekä fluidien epätasainen jakautuminen nousuputkiin. Lisäksi tutkitaan erilaisten muuttujien, kuten kaasun tilavuusvirran, nesteen viskositeetin, suutinkoon ja nesteen jakoputken rakenteen, vaikutusta kaasun tilavuusosuuteen ja nesteen nopeuteen nousuputkissa. Nesteen nopeudet mitataan merkkiainemenetelmällä ja kaasun tilavuusosuudet manometrimenetelmällä. Lämmönsiirtoa tutkitaan mittaamalla lämpötilaeroja nousuputkissa NiCr-Ni –termoelementeillä. Mittaustulosten perusteella muokataan korrelaatiot kaasun tilavuusosuudelle ja nesteen tyhjäputkinopeudelle. Korrelaatioista lasketut tulokset sopivat kohtuullisen hyvin yhteen mitattujen tulosten kanssa. "Slug flown" ei todettu muodostuvan ongelmaksi 2.5 mPa s pienemmillä viskositeetin arvoilla 2 metriä pitkissä ja 19 mm halkaisijaltaan olevissa putkissa. Lisäksi todettiin, että kaasu- ja nestefaasien jakautumisongelmat voidaan ratkaista rakenteellisesti.

Fischer-Tropsch Slurry Bubble Column Reactor Modeling

The literature part of the work reviews overall Fischer-Tropsch process, Fischer-Tropsch reactors and catalysts. Fundamentals of Fischer-Tropsch modeling are also presented. The emphasis is on the reactor unit. Comparison of the reactors and the catalysts is carried out to choose the suitable reactor setup for the modeling work. The effects of the operation conditions are also investigated. Slurry bubble column reactor model operating with cobalt catalyst is developed by taking into account the mass transfer of the reacting components (CO and H2) and the consumption of the reactants in the liquid phase. The effect of hydrostatic pressure and the change in total mole flow rate in gas phase are taken into account in calculation of the solubilities. The hydrodynamics, reaction kinetics and product composition are determined according to literature. The cooling system and furthermore the required heat transfer area and number of cooling tubes are also determined. The model is implemented in Matlab software. Commercial scale reactor setup is modeled and the behavior of the model is investigated. The possible inaccuraries are evaluated and the suggestions for the future work are presented. The model is also integrated to Aspen Plus process simulation software, which enables the usage of the model in more extensive Fischer-Tropsch process simulations. Commercial scale reactor of diameter of 7 m and height of 30 m was modeled. The capacity of the reactor was calculated to be about 9 800 barrels/day with CO conversion of 75 %. The behavior of the model was realistic and results were in the right range. The highest uncertainty to model was estimated to be caused by the determination of the kinetic rate.

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.

Electrochemical Characterisation of Filter Media Properties and Their Exploitation in Enhanced Filtration

Streaming potential measurements for the surface charge characterisation of different filter media types and materials were used. The equipment was developed further so that measurements could be taken along the surfaces, and so that tubular membranes could also be measured. The streaming potential proved to be a very useful tool in the charge analysis of both clean and fouled filter media. Adsorption and fouling could be studied, as could flux, as functions of time. A module to determine the membrane potential was also constructed. The results collected from the experiments conducted with these devices were used in the study of the theory of streaming potential as an electrokinetic phenomenon. Several correction factors, which are derived to take into account the surface conductance and the electrokinetic flow in very narrow capillaries, were tested in practice. The surface materials were studied using FTIR and the results compared with those from the streaming potentials. FTIR analysis was also found to be a useful tool in the characterisation of filters, as well as in the fouling studies. Upon examination of the recorded spectra from different depths in a sample it was possible to determine the adsorption sites. The influence of an external electric field on the cross flow microflltration of a binary protein system was investigated using a membrane electroflltration apparatus. The results showed that a significant improvement could be achieved in membrane filtration by using the measured electrochemical properties to help adjust the process conditions.

Program hosting in communications environment with Reactor and thread pool approach

In this thesis concurrent communication event handling is implemented using thread pool approach. Concurrent events are handled with a Reactor design pattern and multithreading is implemented using a Leader/Followers design pattern. Main focus is to evaluate behaviour of implemented model by different numbers of concurrent connections and amount of used threads. Furthermore, model feasibility in a PeerHood middleware is evaluated. Implemented model is evaluated with created test environment which enables concurrent message sending from multiple connections to the system under test. Messages round trip times are measured in the tester application. In the evaluation processing delay into system is simulated and influence of delay to the average round trip time is analysed.

Design and implementation of a multipurpose reactor for scale-up studies

In the theoretical part, the different polymerisation catalysts are introduced and the phenomena related to mixing in the stirred tank reactor are presented. Also the advantages and challenges related to scale-up are discussed. The aim of the applied part was to design and implement an intermediate-sized reactor useful for scale-up studies. The reactor setting was tested making one batch of Ziegler–Natta polypropylene catalyst. The catalyst preparation with a designed equipment setting succeeded and the catalyst was analysed. The analyses of the catalyst were done, because the properties of the catalyst were compared to the normal properties of Ziegler–Natta polypropylene catalyst. The total titanium content of the catalyst was slightly higher than in normal Ziegler–Natta polypropylene catalyst, but the magnesium and aluminium content of the catalyst were in the normal level. By adjusting the siphonation tube and adding one washing step the titanium content of the catalyst could be decreased. The particle size of the catalyst was small, but the activity was in a normal range. The size of the catalyst particles could be increased by decreasing the stirring speed. During the test run, it was noticed that some improvements for the designed equipment setting could be done. For example more valves for the chemical feed line need to be added to ensure inert conditions during the catalyst preparation. Also nitrogen for the reactor needs to separate from other nitrogen line. With this change the pressure in the reactor can be kept as desired during the catalyst preparation. The proposals for improvements are presented in the applied part. After these improvements are done, the equipment setting is ready for start-up. The computational fluid dynamics model for the designed reactor was provided by cooperation with Lappeenranta University of Technology. The experiments showed that for adequate mixing with one impeller, stirring speed of 600 rpm is needed. The computational fluid dynamics model with two impellers showed that there was no difference in the mixing efficiency if the upper impeller were pumping downwards or upwards.

APROS model validation in case of fast transients in a boiling water reactor

This Master´s thesis investigates the performance of the Olkiluoto 1 and 2 APROS model in case of fast transients. The thesis includes a general description of the Olkiluoto 1 and 2 nuclear power plants and of the most important safety systems. The theoretical background of the APROS code as well as the scope and the content of the Olkiluoto 1 and 2 APROS model are also described. The event sequences of the anticipated operation transients considered in the thesis are presented in detail as they will form the basis for the analysis of the APROS calculation results. The calculated fast operational transient situations comprise loss-of-load cases and two cases related to a inadvertent closure of one main steam isolation valve. As part of the thesis work, the inaccurate initial data values found in the original 1-D reactor core model were corrected. The input data needed for the creation of a more accurate 3-D core model were defined. The analysis of the APROS calculation results showed that while the main results were in good accordance with the measured plant data, also differences were detected. These differences were found to be caused by deficiencies and uncertainties related to the calculation model. According to the results the reactor core and the feedwater systems cause most of the differences between the calculated and measured values. Based on these findings, it will be possible to develop the APROS model further to make it a reliable and accurate tool for the analysis of the operational transients and possible plant modifications.

Electrochemical and spectroelectrochemical characterization of alkylated and thiolated polyanilines

Vi omges i vardagen av alla sorters plaster, som kemiskt kallas för polymerer. Vi anknyter dem oftast till vardagliga föremål såsom muggar, leksaker eller platskassar. Det finns dock en särklass av polymerer som fås elektriskt ledande genom en så kallad dopningsprocess. Dopning innebär i detta fall oxidation eller reduktion av konjugerade dubbelbindningar i polymerstrukturen. Detta har lett till utveckling av elektriska apparater där dyra, och i vissa fall även sällsynta, metall och halvledarmaterial ersätts genom av elektriskt ledande polymerer (plaster). Utöver elektronisk ledningsförmåga uppvisar dessa polymerer också jonisk ledningsförmåga. Denna kombination av unika egenskaper möjliggör skapandet av t.ex. nya sensormaterial som kan överföra kemisk information till en mätbar elektronisk signal. Detta öppnar i sin tur möjligheter att göra snabba, billiga och känsliga sensorer för bl.a. mediciniska analyser. I denna avhandling karakteriserades elektrokemiskt och spektroelektrokemiskt N- och ring-substituerade polyanilinfilmer. Polyanilin (PANI) hör till de mest studerade elektriskt ledande polymererna. Den är stabil och lätt att framställa. Substituerade polyaniliner har ändå studerats måttligt, mest p.g.a. att substituerade PANIs ledningsförmåga är lägre än PANIs och deras framställning kan vara svår. De nya grupperna i PANI-kedjan ger dock en möjlighet att binda ytterligare molekyler av intresse till PANI-kedjan, som t.ex. jonselektiva grupper. Kovalent bundna selektiva molekyler ger upphov till stabila, känsliga och selektiva sensormaterial. Karakteriseringen av de studerade polymerer är viktig för den fundamentala förståelsen av deras unika egenskaper och för utvecklingen av framtidens sensormaterial. -------------------------------------- Käytämme joka päivä monenlaisia muoveja, joita kutsutaan kemiassa myös polymeereiksi. Olemme tottuneet yhdistämään muovit arkisiin esineisiin kuten mukeihin, leluihin tai muovikasseihin. On kuitenkin olemassa erityisiä polymeerejä, jotka voidaan saada sähköä johtaviksi hapetus- ja pelkistysreaktioiden avulla. Tästä johtuen on kehitteillä sähköisiä laitteita, joissa kalliit ja jossain tapauksissa myös harvinaisia metalleja sisältävät osat ja puolijohteet voidaan korvata johdepolymeereillä (eli muoveilla). Sähkönjohtavuuden lisäksi johdepolymeereillä on myös ionijohtavuutta. Näiden erityislaatuisten ominaisuuksien yhdistelmä on mahdollistanut mm. sensorimateriaalin kehittämisen, sillä kemiallinen tieto voidaan kääntää mitattavaksi sähköiseksi signaaliksi. Tämä taas omalta osaltaan mahdollistaisi nopeiden, halpojen ja herkkien sensorien valmistuksen, mm. diagnostiikkaan. Tässä väitöksessä on tutkittu sähkökemiallisesti valmistettuja N- ja rengassubstituoituja polyaniliinikalvoja. Polyaniliini (PANI) on yksi eniten tutkituista johdepolymeereistä. Se on stabiili ja helppo valmistaa. Substituoidut polyaniliinit ovat herättäneet vain kohtalaista tieteellistä kiinnostusta, lähinnä, koska niiden sähköinen johdekyky on alhaisempi kuin PANIn. Myös niiden valmistus voi olla vaikeaa. Substituoidut molekyylit PANI-ketjussa mahdollistavat kuitenkin, että niihin voi liittää uusia molekyylejä, esim. ioniherkkiä ryhmiä. Kovalentisti sitoutuneilla selektiivisillä molekyyleillä saadaan tehtyä stabiileja, herkkiä ja selektiivisiä sensorimateriaaleja. Väitöksessä käytettyjen polymeerien karakterisointi on tärkeää, jotta niiden erityisominaisuuksia pystyttäisiin ymmärtämään paremmin ja myös kehittämään sopivia tulevaisuuden sensorimateriaaleja.

Electrochemical Treatment of Wastewaters

Electrocoagulation is a process in which wastewater is treated under electrical current. Coagulant is formed during the process through the metal anode dissolution to respective ions which react with hydroxyl ions released in cathode. These metal hydroxides form complexes with pollutant ions. Pollutants are removed among metal hydroxide precipitates. This study was concentrated on describing chemistry and device structures in which electrochemical treatment operations are based on. Studied pollutants were nitrogen compounds, sulphate, trivalent and pentavalent arsenic, heavy metals, phosphate, fluoride, chloride, and bromide. In experimental part, removal of ammonium, nitrate, and sulphate during electrochemical treatment was studied separately. Main objective of this study was to find suitable metal plate material for ammonium, nitrate, and sulphate removal, respectively. Also other parameters such as pH of solution, concentration of pollutant and sodium chloride, and current density were optimized. According to this study the most suitable material for ammonium and sulphate removal by electrochemical treatment was stainless steel. Respectively, iron was the optimum material for nitrate removal. Rise in the pH of solution at the final stage of electrochemical treatment of ammonium, nitrate, and sulphate was detected. Conductivities of solutions decreased during ammonium removal in electrochemical processes. When nitrate and sulphate were removed electrochemically conductivities of solutions increased. Concentrations of residual metals in electrochemically treated solutions were not significant. Based on this study electrochemical treatment processes are recommended to be used in treatment of industrial wastewaters. Treatment conditions should be optimized for each wastewater matrix.

Development of a process for the direct synthesis of hydrogen peroxide in a novel microstructured reactor

Microreactors have proven to be versatile tools for process intensification. Over recent decades, they have increasingly been used for product and process development in chemical industries. Enhanced heat and mass transfer in the reactors due to the extremely high surfacearea- to-volume ratio and interfacial area allow chemical processes to be operated at extreme conditions. Safety is improved by the small holdup volume of the reactors and effective control of pressure and temperature. Hydrogen peroxide is a powerful green oxidant that is used in a wide range of industries. Reduction and auto-oxidation of anthraquinones is currently the main process for hydrogen peroxide production. Direct synthesis is a green alternative and has potential for on-site production. However, there are two limitations: safety concerns because of the explosive gas mixture produced and low selectivity of the process. The aim of this thesis was to develop a process for direct synthesis of hydrogen peroxide utilizing microreactor technology. Experimental and numerical approaches were applied for development of the microreactor. Development of a novel microreactor was commenced by studying the hydrodynamics and mass transfer in prototype microreactor plates. The prototypes were designed and fabricated with the assistance of CFD modeling to optimize the shape and size of the microstructure. Empirical correlations for the mass transfer coefficient were derived. The pressure drop in micro T-mixers was investigated experimentally and numerically. Correlations describing the friction factor for different flow regimes were developed and predicted values were in good agreement with experimental results. Experimental studies were conducted to develop a highly active and selective catalyst with a proper form for the microreactor. Pd catalysts supported on activated carbon cloths were prepared by different treatments during the catalyst preparation. A variety of characterization methods were used for catalyst investigation. The surface chemistry of the support and the oxidation state of the metallic phase in the catalyst play important roles in catalyst activity and selectivity for the direct synthesis. The direct synthesis of hydrogen peroxide was investigated in a bench-scale continuous process using the novel microreactor developed. The microreactor was fabricated based on the hydrodynamic and mass transfer studies and provided a high interfacial area and high mass transfer coefficient. The catalysts were prepared under optimum treatment conditions. The direct synthesis was conducted at various conditions. The thesis represents a step towards a commercially viable direct synthesis. The focus is on the two main challenges: mitigating the safety problem by utilization of microprocess technology and improving the selectivity by catalyst development.

Catalytic direct synthesis of hydrogen peroxide in a novel microstructured reactor

The direct synthesis from hydrogen and oxygen is a green alternative for production of hydrogen peroxide. However, this process suffers from two challenges. Firstly, mixtures of hydrogen and oxygen are explosive over a wide range of concentrations (4-94% H2 in O2). Secondly, the catalytic reaction of hydrogen and oxygen involves several reaction pathways, many of them resulting in water production and therfore decreasing selectivity. The present work deals with these two challenges. The safety problem was dealed by employing a novel microstructured reactor. Selectivity of the reaction was highly improved by development a set of new catalysts. The final goal was to develop an effective and safe continuous process for direct synthesis of hydrogen peroxide from H2 and O2. Activated carbon cloth and Sibunit were examined as the catalysts’ supports. Palladium and gold monometallic and palladium-gold bimetallic catalysts were thoroughly investigated by numerous kinetic experiments performed in a tailored batch reactor and several catalyst charachterization methods. A complete set of data for direct synthesis of H2O2 and its catalytic decomposition and hydrogenation was obtained. These data were used to assess factors influencing selectivity and activity of the catalysts in direct synthesis of H2O2 as well as its decomposition and hydrogenation. A novel microstructured reactor was developed based on hydrodynamics and mass transfer studies in prototype microstractural plates. The shape and the size of the structural elements in the microreactor plate were optimized in a way to get high gas-liquid interfacial area and gas-liquid mass transfer. Finally, empirical correlations for the volumetric mass transfer coefficient were derived. A bench-scale continuous process was developed by using the novel microstructral plate reactor. A series of kinetic experiments were performed to investigate the effects of the gas and the liquid feed rates and their ratio, the amount of the catalyst, the gas feed composition and pressure on the final rate of H2O2 production and selectivity.

Development of an HTR-10 model in the Serpent reactor physics code

The use of exact coordinates of pebbles and fuel particles of pebble bed reactor modelling becoming possible in Monte Carlo reactor physics calculations is an important development step. This allows exact modelling of pebble bed reactors with realistic pebble beds without the placing of pebbles in regular lattices. In this study the multiplication coefficient of the HTR-10 pebble bed reactor is calculated with the Serpent reactor physics code and, using this multiplication coefficient, the amount of pebbles required for the critical load of the reactor. The multiplication coefficient is calculated using pebble beds produced with the discrete element method and three different material libraries in order to compare the results. The received results are lower than those from measured at the experimental reactor and somewhat lower than those gained with other codes in earlier studies.