Theoretical and numerical study of thermomagnetic convection in magnetic fluids

In this thesis, the magnetic field control of convection instabilities and heat and mass transfer processesin magnetic fluids have been investigated by numerical simulations and theoretical considerations. Simulation models based on finite element and finite volume methods have been developed. In addition to standard conservation equations, themagnetic field inside the simulation domain is calculated from Maxwell equations and the necessary terms to take into account for the magnetic body force and magnetic dissipation have been added to the equations governing the fluid motion.Numerical simulations of magnetic fluid convection near the threshold supportedexperimental observations qualitatively. Near the onset of convection the competitive action of thermal and concentration density gradients leads to mostly spatiotemporally chaotic convection with oscillatory and travelling wave regimes, previously observed in binary mixtures and nematic liquid crystals. In many applications of magnetic fluids, the heat and mass transfer processes including the effects of external magnetic fields are of great importance. In addition to magnetic fluids, the concepts and the simulation models used in this study may be applied also to the studies of convective instabilities in ordinary fluids as well as in other binary mixtures and complex fluids.

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.

Development of a PC Tool for Estimating Pressure Drops in Forced-convection Boiler Tubes

Kaksifaasivirtauksen kuvaamiseen käytettävät mallit, ja menetelmät kaksifaasivirtauksen painehäviön määrittämiseksi kehittyvät yhä monimutkaisimmiksi. Höyrystinputkissa tapahtuvien painehäviöiden arvioinnin vaatiman laskennan suorittamiseksi tietokoneohjelman kehittäminen on välttämätöntä. Tässä työssä on kehitetty itsenäinen PC-ohjelma painehäviöiden arvioimiseksi pakotetulle konvektiovirtaukselle pystysuorissa höyrykattilan höyrystinputkissa. Veden ja vesihöyryn aineominaisuuksien laskentaan käytetään IAPWS-IF97 –yhtälökokoelmaa sekä muita tarvittavia IAPWS:n suosittelemia yhtälöitä. Höyrystinputkessa kulloinkin vallitsevan virtausmuodon määrittämiseen käytetään sovelluskelpoisia virtausmuotojen välisiä rajoja kuvaavia yhtälöitä. Ohjelmassa käytetään painehäviön määritykseen kirjallisuudessa julkaistuja yhtälöitä, virtausmuodosta riippuen, alijäähtyneelle virtaukselle, kupla-, tulppa- ja rengasvirtaukselle sekä tulistetun höyryn virtaukselle. Ohjelman laskemia painehäviöarvioita verrattiin kirjallisuudesta valittuihin mittaustuloksiin. Laskettujen painehäviöiden virhe vaihteli välillä –19.5 ja +23.9 %. Virheiden itseisarvojen keskiarvo oli 12.8 %.

A generalized 1D particle transport method for convection diffusion reaction model

This work is devoted to the development of numerical method to deal with convection diffusion dominated problem with reaction term, non - stiff chemical reaction and stiff chemical reaction. The technique is based on the unifying Eulerian - Lagrangian schemes (particle transport method) under the framework of operator splitting method. In the computational domain, the particle set is assigned to solve the convection reaction subproblem along the characteristic curves created by convective velocity. At each time step, convection, diffusion and reaction terms are solved separately by assuming that, each phenomenon occurs separately in a sequential fashion. Moreover, adaptivities and projection techniques are used to add particles in the regions of high gradients (steep fronts) and discontinuities and transfer a solution from particle set onto grid point respectively. The numerical results show that, the particle transport method has improved the solutions of CDR problems. Nevertheless, the method is time consumer when compared with other classical technique e.g., method of lines. Apart from this advantage, the particle transport method can be used to simulate problems that involve movingsteep/smooth fronts such as separation of two or more elements in the system.

Lämmönlähteiden aiheuttamien konvektiovirtausten vaikutus vastakkaiseen tuloilmasuihkuun

Huoneen sisäisten lämpölähteiden aiheuttamilla konvektiovirtauksilla eli pluumeilla on huomattava vaikutus ilmastointipalkkijärjestelmäntoimintaan. Tätä vaikutusta ei ole tutkittu kovin paljon. Tässä työssä tutkittiin huoneen päätyseinässä olevan sähkölämmityspatterin vaikutusta aktiivisen ilmastointipalkin toimintaan. Tutkimus kohdistui vastakkaisten suihkujen törmäykseen. Suihkuina olivat pluumit ja tuloilmavirtaus. Tutkimukseen kuului sekä teoreettinen että kokeellinen osio. Teoreettisessa tarkastelussa käytiin läpi ilmastointipalkin teoriaa sekä konvektiovirtausten muodostumisen teoriaa. Lisäksi esitettiin, kuinka suihkuista lasketaan liikemäärävirta. Kokeellisessa osiossa mitattiin huoneen sisäisten lämpölähteiden aiheuttamia pluumeja. Lämpölähteet oli sijoitettuina tutkimushuoneen päätyseinän viereen. Mittauksissa keskityttiin tarkemmin sähkölämmityspatteriin ja sitä mallinnettiin eri tehoilla. Lisäksi mitattiin ilmastointipalkin tuloilmavirtausta eri ilmamäärillä. Pluumien ja tuloilmavirtauksen törmäyskäyttäytymistä selvitettiin savukuvauksella. Analyysiosiossa sähkölämmityspatterin aiheuttaman pluumin ja palkin tuloilmavirtauksen törmäyskäyttäytymiselle kehitettiin malli. Mallissa huomioitiin vastaikkaisten törmäävien suihkujen liikemäärät ja esitettiin kohdat tutkimushuoneen pituussuunnassa, jossa tuloilmavirtaus kääntyy oleskeluvyöhykkeelle. Ilmastointipalkkijärjestelmässä tuloilmavirtauksen tulisi kääntyä oleskeluvyöhykkeelle vasta huoneen seinustalla. Tutkimuksen perusteella näin käy, jos pluumin liikemäärä on noin 10% tai vähemmän tuloilman liikemäärästä. Muulloin palkin tuloilmavirtaus kääntyy oleskeluvyöhykkeelle liian aikaisin ja aiheuttaa mahdollisestivedontunnetta.

Uusien energianormien vaikutus pientalojen rakenteisiin ympäristön kannalta

Tehokkaimpia keinoja vähentää rakennusten lämmitysenergian kulutusta ja lämmityksen aiheuttavia hiilidioksidi- ja happamoitavia päästöjä on tiukentaa rakentamismääräysten lämmöneristysvaatimuksia. Hyvin lämmöneristetyissä, tiiveissä ja ilmanvaihdoltaan optimoiduissa taloissa on pienet lämpöhäviöt. Näin ympäristöä kuormittava vaikutus saadaan paljon vähemmäksi kuin nykynormien mukaisissa asuinrakennuksissa. Johtumislämpöhäviö pienenee suoraan eristekerroksia paksuntamalla ja siihen on helpointa vaikuttaa. Mitä suurempiin eristepaksuuksiin mennään sen suuremmaksi tulee konvektion osuus kokonaislämpöhäviöstä. Tulevaisuudessa parempia ratkaisuja haetaan erityisesti konvektiosta ja säteilystä aiheutuvien lämpöhäviöiden pienentämiseksi. Eristeen osastointi ilmanpitävillä, vesihöyryä diffuusisesti läpäisevillä pystysuuntaisilla konvektiokatkoilla vähentää tehokkaasti paksun seinäeristeen kuljettumis-ilmavirtauksia. Katkoina käytetään erilaisia kalvoja ja rakennuspapereita, joilla on pieni emissiviteetti. Katkojen merkitys kasvaa, kun mennään uusien normien mukaisiin eristepaksuuksiin. Lämmöneriste voidaan toteuttaa myös kokoamalla ohuita kalvoja paketiksi, jotka jakavat ilmatilan ja siis eristeelle varatun paksuuden suljettuihin ilmaväleihin. Kun kalvoiksi valitaan pieniemissiviteettisiä pintoja, saadaan säteilylämmönsiirto lähes eliminoiduksi. Tällaisen ilmatilan lämmönjohtumisluku lähestyy paikallaan pysyvän ilman lämmönjohtumislukua, l = 0,025 W/Km, eli tällä rakennesysteemillä on mahdollista toteuttaa ohuempia rakenteita kuin perinteisillä eristeillä. Hygroskooppisen massan käyttö sisäilman kosteutta tasaavana rakenteena voi olla tulevaisuutta. Kehitystyö tuottaa uusia, kosteusteknisesti toimivia sovelluksia. Toisaalta palomääräykset tulevat kehitystyötä vastaan. Hygroskooppinen pintamateriaali on kevyt (pieni tiheys) ja paloteknisesti arka. Suoraa sähkölämmitystä ei voida pitää ympäristöystävällisenä. Sen jalostusketju on pitkä ja monivaiheinen. Millä peruspolttoaineella sähköä tuotetaan, vaikuttaa asiaan luonnollisestikin. Suoraa sähkölämmitystä voidaan suositella vain yksinäisen ihmisen taloudessa lämmitysmuotona taloudellisista syistä. Halvan polttoaineen säästöllä ei voida maksaa suuria laiteinvestointeja. Aurinkoenergian hyvä hyödyntäminen edellyttää hyvää säätöä, joka kytkee lämmityksen pois päältä silloin, kun aurinko lämmittää. Auringon hetkelliset säteilytehot ovat suuria verrattuna rakenteen lämpöhäviöihin ja huonetilojen lämmöntarpeeseen. Ratkaisu aurinkoenergian hetkellisyyteen ja paikallisuuteen on energian siirtäminen lämmöntarpeen mukaan rakennuksen eri osiin ja sen varastoiminen päivätasolla. Kun varastoivasta massasta ei ole suoraa yhteyttä ulos, voidaan kerääjäeristeeltä saatu lämpö käyttää häviöttömästi huonetilojen lämmittämiseen. Vaikka lämmitysenergian käytössä päästään 30 % vähennyksiin uudisrakennusten osalta, ei kokonaisenergian käyttö merkittävästi pienene, jos taloussähkön kulutus pysyy vakiona. Sama pätee myös CO2 -päästöihin. Saavutettava etu lämmitys-energian kulutuksessa voidaan hukata yhä suurenevaksi taloussähkön käytöksi, mikä olisi erityisen huono asia ympäristön kannalta.

Lämmönsiirron numeerinen laskenta profiloidulle kanavalle.

Työn tavoitteena oli tutkia aaltomaisen profiloinnin vaikutusta suorien jäähdytyskanavien lämmönsiirtoon ja painehäviöön. Erilaisia profiileja oli kymmenen kappaletta ja ne olivat 5mm leveitä ja 30cm pitkiä kukin. Ne laskettiin kolmeulotteisina tapauksina FINFLO-virtausratkaisijalla kolmella eri Reynoldsin luvulla, jotka vastasivat laminaarista, osittain turbulenttista ja lähes kokonaan turbulenttista virtausta. Lämmönsiirtoaine oli kuiva +30°C ilma ja profiloinnin toteutustapa oli toisiaan sivuavat ympyräkaaret kolmella erilaisella säteen arvolla ja kolmella erilaisella aallonpituuden arvolla. Lisäksi laskettiin saman levyisen tasokanavan arvot jokaisella Reynoldsin luvulla kaksiulotteisina tapauksina. Näitä profiloimattomia kanavia pidettiin referenssitapauksina. Tuloksena havaittiin että profiloimalla saadaan yksiselitteisesti suurempi lämpöteho ulos samasta tilavuudesta. Lämmönsiirtokerroin kasvaa profiloinnin avulla parhaimmillaan n. 20% käytetystä turbulenssimallista tai lämmönsiirtokertoimen määritelmästä riippumatta. Painehäviö kasvaa myös aina, mutta kitkakerroin voi hieman pienentyä. Profiilin varsinaisena hyvyyskriteerinä pidettiin lämmönsiirtokertoimen ja kitkakertoimen suhdetta h/f. Se osoittautui riippuvaksi Reynoldsin luvusta ja turbulenssimallista; ASM ja Chien k-έ -mallit ennustavat transitioetäisyyden eri tavalla. Laminaarisilla virtauksilla h/f :n vaihtelu oli vähäistä; suhde vaihteli vain ±5% eri profiilien kesken. ASM-mallilla havaittiin sekundääripyörteilyä, ehkä siksi että se mallintaa anisotrooppisen turbulenssin. Chien k-έ malli ennusti suuremman ja aikaisemmin alkavan turbulenttisuuden kuin ASM. Lisäksi havaittiin mm. että tietyillä profiileilla muodostuu kanavan kapeimpaan kohtaan selvä virtausnopeuden paikallinen minimi seinämän läheisyyden takia.

Kupariliekkisulatusuunin lampohavioiden maarittaminen

Työn teoriaosassa käydään läpi kuparisulaton prosessi pääpiirteittäin, jonka jälkeen siirrytään lämmönsiirron teoriaan. Lämpö siirtyy kolmella penis mekanismilla: johtumalla, säteilemällä ja konvektiolla. Naiden ilmiöiden vaikutusta kupariliekkisulatusuunin eri osissa tapahtuviin lämpöhäviöihin on käyty läpi. Lisaksi työstä löytyy kupariliekkisulatusuunin lämpö- ja ainetase sekä erilaisten lämpötilanmittauslaitteiden vertailua sekä mitä tulee ottaa huomioon lämpötiloja mitatessa. Työn kokeellisessa osassa esitellään lämpötilanmittausten perusteella lasketut kupariliekkisulatusuunin eri osien lämpöhäviöiden suuruudet mittauskampanjan ajalta. Saatuja tuloksia on vertailtu kirjallisuudesta löydettyihin aikaisempiin arvoihin. Uunin keskimääräiseksi kokonaislämpöhäviöksi saatiin 9,8 MW mittauskampanjan ajalle. Lisäksi on vertailtu uunin erilaisten ajoarvojen, kuten rikaste-, lentopöly-, murska-, ja kuonarikastesyötön suuruuden sekä öljyn syötön ja rikasteseoksen vaihtumisen vaikutusta lämpöhäviöihin sekä mietitty muita mahdollisia syitä lämpöhäviöiden vaihtelulle.

Water circulations for effective bleaching of high-brightness mechanical pulps

The amount of water available is usually restricted, which leads to a situation where a complete understanding of the process, including water circulations and the influence of water components, is essential. The main aim of this thesis was to clarify the possibilities for the efficient use of residual peroxide by means of water circulation rearrangements. Rearranging water circulations and the reduction of water usage may cause new problems, such as metal induced peroxide decomposition that needs to be addressed. This thesis introduces theoretical methods of water circulations to combine two variables; effective utilization of residual peroxide and avoiding manganese in the alkaline peroxide bleaching stage. Results are mainly based on laboratory and mill site experiments concerning the utilization of residual peroxide. A simulation model (BALAS) was used to evaluate the manganese contents and residual peroxide doses. It was shown that with optimum recirculation of residual peroxide the brightness can be improved or chemical costs can be decreased. From the scientific perspective, it was also very important to discover that recycled peroxide was more effective pre-bleaching agent compared to fresh peroxide. This can be due to the organic acids i.e. per acetic acid in wash press filtrate that have been formed in alkaline bleaching stage. Even short retention time was adequate and the activation of residual peroxide using sodium hydroxide was not necessary. There are several possibilities for using residual peroxide in practice regarding bleaching. A typical modern mechanical pulping process line consist of defibering, screening, a disc filter, a bleach press, high consistency (HC) peroxide bleaching and a wash press. Furthermore there usually is not a particular medium consistency (MC) pre-bleaching stage that includes additional thickening equipment. The most advisable way to utilize residual peroxide in this kind of process is to recycle the wash press filtrate to the dilution of disc filter pulp (low MC pre-bleaching stage). An arrangement such as this would be beneficial in terms of the reduced convection of manganese to the alkaline bleaching stage. Manganese originates from wood material and will be removed to the water phase already in the early stages of the process. Recycling residual peroxide prior to the disc filter is not recommended because of low consistencies. Regarding water circulations, the novel point of view is that, it would be beneficial to divide water circulations into two sections and the critical location for the division is the disc filter. Both of these two sections have their own priority. Section one before the disc filter: manganese removal. Section two after the disc filter: brightening of pulp. This division can be carried out if the disc filter pulp is diluted only by wash press filtrate before the MC storage tower. The situation is even better if there is an additional press after the disc filter, which will improve the consistency of the pulp. This has a significant effect on the peroxide concentration in the MC pre-bleaching stage. In terms of manganese content, it is essential to avoid the use of disc filter filtrate in the bleach press and wash press showers. An additional cut-off press would also be beneficial for manganese removal. As a combination of higher initial brightness and lower manganese content, the typical brightness increase varies between approximately 0.5 and 1% ISO units after the alkaline peroxide bleaching stage. This improvement does not seem to be remarkable, but as it is generally known, the final brightness unit is the most expensive and difficult to achieve. The estimation of cost savings is not unambiguous. For example in GW/TMP mill case 0.6% ISO units higher final brightness gave 10% savings in the costs of bleaching chemicals. With an hypothetical 200 000 ton annual production, this means that the mill could save in the costs of bleaching chemicals more than 400 000 euros per year. In general, it can be said that there were no differences between the behavior of different types of processes (GW, PGW, TMP and BCTMP). The enhancement of recycling gave a similar response in all cases. However, we have to remember that the utilization of residual peroxide in older mills depends a great deal on the process equipment, the amount of water available and existing pipeline connections. In summary, it can be said that processes are individual and the same solutions cannot be applied to all cases.

Kestomagneettigeneraattorin lämmönsiirron mallinnus tuulivoimakäytössä

Pyörivien sähkökoneiden suunnittelussa terminen suunnittelu on yhtä tärkeää kuin sähköinen ja mekaaninen suunnittelukin. Tässä diplomityössä tarkoituksena on kehittää ilmajäähdytteisten kestomagneettigeneraattorien laskentaan soveltuva lämmönsiirtymismalli, jolla staattorin lämpötilajakauma voitaisiin selvittää. Kehitetty lämmönsiirtymismalli perustuu kolmiulotteiseen äärellisen erotuksen (finite difference) menetelmään. Malli ottaa huomioon lämmönjohtumisen staattorin aktiiviosissa ja konvektion jäähdytysilmavirtaan. Mallissa on myös yksinkertainen painehäviölaskenta jäähdytysjärjestelmän komponenttien mitoittamista varten. Laskentamallilla lasketaan esimerkkitapauksena 4,3 MW:n kestomagneettigeneraattorin jäähdytystä eri toimintapisteissä. Tuloksia verrataan CFD-mallinnuksen antamiin tuloksiin sekä kokeellisten mittausten antamiin tuloksiin.

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.

Fundamentals of heat transfer

Following over 170+ pages and additional appendixes are formed based on content of Course: Fundamentals of Heat Transfer. Mainly this summarizes relevant parts on Book of Fundamentals of Heat and Mass Transfer (Incropera), but also other references introducing the same concepts are included. Student’s point of view has been consideredwith following highlights: (1) Relevant topics are presented in a nutshell to provide fast digestion of principles of heat transfer. (2) Appendixes include terminology dictionary. (3) Totally 22 illustrating examples are connecting theory to practical applications and quantifying heat transfer to understandable forms as: temperatures, heat transfer rates, heat fluxes, resistances and etc. (4) Most important Learning outcomes are presented for each topic separately. The Book, Fundamentals of Heat and Mass Transfer (Incropera), is certainly recommended for those going beyond basic knowledge of heat transfer. Lecture Notes consists of four primary content-wise objectives: (1) Give understanding to physical mechanisms of heat transfer, (2)Present basic concepts and terminology relevant for conduction, convection and radiation (3) Introduce thermal performance analysis methods for steady state and transient conduction systems. (4) Provide fast-to-digest phenomenological understanding required for basic design of thermal models

The bifurcational behaviour of the spatially distributed Rayleigh equation

At the present work the bifurcational behaviour of the solutions of Rayleigh equation and corresponding spatially distributed system is being analysed. The conditions of oscillatory and monotonic loss of stability are obtained. In the case of oscillatory loss of stability, the analysis of linear spectral problem is being performed. For nonlinear problem, recurrent formulas for the general term of the asymptotic approximation of the self-oscillations are found, the stability of the periodic mode is analysed. Lyapunov-Schmidt method is being used for asymptotic approximation. The correlation between periodic solutions of ODE and PDE is being investigated. The influence of the diffusion on the frequency of self-oscillations is being analysed. Several numerical experiments are being performed in order to support theoretical findings.

Modelling of hydrogen stratification and gas mixing in containment using the APROS code

Hydrogen stratification and atmosphere mixing is a very important phenomenon in nuclear reactor containments when severe accidents are studied and simulated. Hydrogen generation, distribution and accumulation in certain parts of containment may pose a great risk to pressure increase induced by hydrogen combustion, and thus, challenge the integrity of NPP containment. The accurate prediction of hydrogen distribution is important with respect to the safety design of a NPP. Modelling methods typically used for containment analyses include both lumped parameter and field codes. The lumped parameter method is universally used in the containment codes, because its versatility, flexibility and simplicity. The lumped parameter method allows fast, full-scale simulations, where different containment geometries with relevant engineering safety features can be modelled. Lumped parameter gas stratification and mixing modelling methods are presented and discussed in this master’s thesis. Experimental research is widely used in containment analyses. The HM-2 experiment related to hydrogen stratification and mixing conducted at the THAI facility in Germany is calculated with the APROS lump parameter containment package and the APROS 6-equation thermal hydraulic model. The main purpose was to study, whether the convection term included in the momentum conservation equation of the 6-equation modelling gives some remarkable advantages compared to the simplified lumped parameter approach. Finally, a simple containment test case (high steam release to a narrow steam generator room inside a large dry containment) was calculated with both APROS models. In this case, the aim was to determine the extreme containment conditions, where the effect of convection term was supposed to be possibly high. Calculation results showed that both the APROS containment and the 6-equation model could model the hydrogen stratification in the THAI test well, if the vertical nodalisation was dense enough. However, in more complicated cases, the numerical diffusion may distort the results. Calculation of light gas stratification could be probably improved by applying the second order discretisation scheme for the modelling of gas flows. If the gas flows are relatively high, the convection term of the momentum equation is necessary to model the pressure differences between the adjacent nodes reasonably.

Rheological characterization of polysaccharide gels

In this Master’s Thesis work the rheological properties of different polysaccharide gels have been studied. The results of this study are used as a starting point for further investigations of potential applications. In order to understand rheological behavior of studied materials, the commercial hydrocolloids such as sodium carboxymethyl cellulose, xanthan gum and guar gum were used as reference and comparison material for rheological studies. As a part the rheological research the development and implementation of proper measurement methods for studied materials were carried out. In the literature review, short introductions of studied materials and application areas of rheological modifiers are summarized. In addition, basic rheological concepts and key fundamentals are explained. In the experimental part the focus was on the rheological characterization of aqueous suspensions of studied materials. Especially, gel strength and solution stability were investigated. The rheological measurements included both rotational and oscillatory measurements in different conditions, where several chemical and physical properties were measured with Anton Paar MCR302 dynamic rotational rheometer. Studied polysaccharide gels can be clearly defined to be shear thinning and thixotropic materials. They have strong gel forming properties even at low concentrations, which explains the superior thickening behavior for some of the samples. Along with rheological characterization of selected materials the factors behind different phenomena were investigated. To reveal value and potential use of polysaccharide gels the influence of various factors such as concentration, temperature and ionic strength were determined. The measurements showed a clear difference between studied materials under investigated external parameters.