Fluid flow and dispersion predictions in porous media by utilizing advanced numerical methods

Diplomityön tavoitteena oli tarkastella numeerisen virtauslaskennan avulla virtaukseen liittyviä ilmiöitä ja kaasun dispersiota. Diplomityön sisältö on jaettu viiteen osaan; johdantoon, teoriaan, katsaukseen virtauksen mallinnukseen huokoisessa materiaalissa liittyviin tutkimusselvityksiin, numeeriseen mallinnukseen sekä tulosten esittämiseen ja johtopäätöksiin. Diplomityön alussa kiinnitettiin huomiota erilaisiin kokeellisiin, numeerisiin ja teoreettisiin mallinnusmenetelmiin, joilla voidaan mallintaa virtausta huokoisessa materiaalissa. Kirjallisuusosassa tehtiin katsaus aikaisemmin julkaistuihin puoliempiirisiin ja empiirisiin tutkimusselvityksiin, jotka liittyvät huokoisen materiaalin aiheuttamaan painehäviöön. Numeerisessa virtauslaskenta osassa rakennettiin ja esitettiin huokoista materiaalia kuvaavat numeeriset mallit käyttäen kaupallista FLUENT -ohjelmistoa. Työn lopussa arvioitiin teorian, numeerisen virtauslaskennan ja kokeellisten tutkimusselvitysten tuloksia. Kolmiulotteisen huokoisen materiaalinnumeerisessa mallinnuksesta saadut tulokset vaikuttivat lupaavilta. Näiden tulosten perusteella tehtiin suosituksia ajatellen tulevaa virtauksen mallinnusta huokoisessa materiaalissa. Osa tässä diplomityössä esitetyistä tuloksista tullaan esittämään 55. Kanadan Kemiantekniikan konferenssissa Torontossa 1619 Lokakuussa 2005. ASME :n kansainvälisessä tekniikan alan julkaisussa. Työ on hyväksytty esitettäväksi esitettäväksi laskennallisen virtausmekaniikan (CFD) aihealueessa 'Peruskäsitteet'. Lisäksi työn yksityiskohtaiset tulokset tullaan lähettämään myös CES:n julkaisuun.

Calculations of Boiling Two-Phase Flow Using a Porous Media Model

Boiling two-phase flow and the equations governing the motion of fluid in two-phase flows are discussed in this thesis. Disposition of the governing equations in three-dimensional complex geometries is considered from the perspective of the porous medium concept. The equations governing motion in two-phase flows were formulated, discretized and implemented in a subroutine for pressure-velocity solution utilizing the SIMPLE algorithm modified for two-phase flow. The subroutine was included in PORFLO, which is a three-dimensional 5-equation porous media model developed at VTT by Jaakko Miettinen. The development of two-phase flow and the resulting void fraction distribution was predicted in a geometry resembling a section of BWR fuel bundle in a couple of test cases using PORFLO.

Modeling the transport phenomena within arterial wall: porous media approach

The transport of macromolecules, such as low-density lipoprotein (LDL), and their accumulation in the layers of the arterial wall play a critical role in the creation and development of atherosclerosis. Atherosclerosis is a disease of large arteries e.g., the aorta, coronary, carotid, and other proximal arteries that involves a distinctive accumulation of LDL and other lipid-bearing materials in the arterial wall. Over time, plaque hardens and narrows the arteries. The flow of oxygen-rich blood to organs and other parts of the body is reduced. This can lead to serious problems, including heart attack, stroke, or even death. It has been proven that the accumulation of macromolecules in the arterial wall depends not only on the ease with which materials enter the wall, but also on the hindrance to the passage of materials out of the wall posed by underlying layers. Therefore, attention was drawn to the fact that the wall structure of large arteries is different than other vessels which are disease-resistant. Atherosclerosis tends to be localized in regions of curvature and branching in arteries where fluid shear stress (shear rate) and other fluid mechanical characteristics deviate from their normal spatial and temporal distribution patterns in straight vessels. On the other hand, the smooth muscle cells (SMCs) residing in the media layer of the arterial wall respond to mechanical stimuli, such as shear stress. Shear stress may affect SMC proliferation and migration from the media layer to intima. This occurs in atherosclerosis and intimal hyperplasia. The study of blood flow and other body fluids and of heat transport through the arterial wall is one of the advanced applications of porous media in recent years. The arterial wall may be modeled in both macroscopic (as a continuous porous medium) and microscopic scales (as a heterogeneous porous medium). In the present study, the governing equations of mass, heat and momentum transport have been solved for different species and interstitial fluid within the arterial wall by means of computational fluid dynamics (CFD). Simulation models are based on the finite element (FE) and finite volume (FV) methods. The wall structure has been modeled by assuming the wall layers as porous media with different properties. In order to study the heat transport through human tissues, the simulations have been carried out for a non-homogeneous model of porous media. The tissue is composed of blood vessels, cells, and an interstitium. The interstitium consists of interstitial fluid and extracellular fibers. Numerical simulations are performed in a two-dimensional (2D) model to realize the effect of the shape and configuration of the discrete phase on the convective and conductive features of heat transfer, e.g. the interstitium of biological tissues. On the other hand, the governing equations of momentum and mass transport have been solved in the heterogeneous porous media model of the media layer, which has a major role in the transport and accumulation of solutes across the arterial wall. The transport of Adenosine 5´-triphosphate (ATP) is simulated across the media layer as a benchmark to observe how SMCs affect on the species mass transport. In addition, the transport of interstitial fluid has been simulated while the deformation of the media layer (due to high blood pressure) and its constituents such as SMCs are also involved in the model. In this context, the effect of pressure variation on shear stress is investigated over SMCs induced by the interstitial flow both in 2D and three-dimensional (3D) geometries for the media layer. The influence of hypertension (high pressure) on the transport of lowdensity lipoprotein (LDL) through deformable arterial wall layers is also studied. This is due to the pressure-driven convective flow across the arterial wall. The intima and media layers are assumed as homogeneous porous media. The results of the present study reveal that ATP concentration over the surface of SMCs and within the bulk of the media layer is significantly dependent on the distribution of cells. Moreover, the shear stress magnitude and distribution over the SMC surface are affected by transmural pressure and the deformation of the media layer of the aorta wall. This work reflects the fact that the second or even subsequent layers of SMCs may bear shear stresses of the same order of magnitude as the first layer does if cells are arranged in an arbitrary manner. This study has brought new insights into the simulation of the arterial wall, as the previous simplifications have been ignored. The configurations of SMCs used here with elliptic cross sections of SMCs closely resemble the physiological conditions of cells. Moreover, the deformation of SMCs with high transmural pressure which follows the media layer compaction has been studied for the first time. On the other hand, results demonstrate that LDL concentration through the intima and media layers changes significantly as wall layers compress with transmural pressure. It was also noticed that the fraction of leaky junctions across the endothelial cells and the area fraction of fenestral pores over the internal elastic lamina affect the LDL distribution dramatically through the thoracic aorta wall. The simulation techniques introduced in this work can also trigger new ideas for simulating porous media involved in any biomedical, biomechanical, chemical, and environmental engineering applications.

Non-resorbable glass fibre-reinforced composite with porous surface as bone substitute material: Experimental studies in vitro and in vivo focused on bone-implant interface

The development of load-bearing osseous implant with desired mechanical and surface properties in order to promote incorporation with bone and to eliminate risk of bone resorption and implant failure is a very challenging task. Bone formation and resoption processes depend on the mechanical environment. Certain stress/strain conditions are required to promote new bone growth and to prevent bone mass loss. Conventional metallic implants with high stiffness carry most of the load and the surrounding bone becomes virtually unloaded and inactive. Fibre-reinforced composites offer an interesting alternative to metallic implants, because their mechanical properties can be tailored to be equal to those of bone, by the careful selection of matrix polymer, type of fibres, fibre volume fraction, orientation and length. Successful load transfer at bone-implant interface requires proper fixation between the bone and implant. One promising method to promote fixation is to prepare implants with porous surface. Bone ingrowth into porous surface structure stabilises the system and improves clinical success of the implant. The experimental part of this work was focused on polymethyl methacrylate (PMMA) -based composites with dense load-bearing core and porous surface. Three-dimensionally randomly orientated chopped glass fibres were used to reinforce the composite. A method to fabricate those composites was developed by a solvent treatment technique and some characterisations concerning the functionality of the surface structure were made in vitro and in vivo. Scanning electron microscope observations revealed that the pore size and interconnective porous architecture of the surface layer of the fibre-reinforced composite (FRC) could be optimal for bone ingrowth. Microhardness measurements showed that the solvent treatment did not have an effect on the mechanical properties of the load-bearing core. A push-out test, using dental stone as a bone model material, revealed that short glass fibre-reinforced porous surface layer is strong enough to carry load. Unreacted monomers can cause the chemical necrosis of the tissue, but the levels of leachable resisidual monomers were considerably lower than those found in chemically cured fibre-reinforced dentures and in modified acrylic bone cements. Animal experiments proved that surface porous FRC implant can enhance fixation between bone and FRC. New bone ingrowth into the pores was detected and strong interlocking between bone and the implant was achieved.

Porous silicon optical filters in gas sensing applications

In this thesis, the gas sensing properties of porous silicon-based thin-film optical filters are explored. The effects of surface chemistry on the adsorption and desorption of various gases are studied in detail. Special emphasis is placed on investigating thermal carbonization as a stabilization method for optical sensing applications. Moreover, the possibility of utilizing the increased electrical conductivity of thermally carbonized porous silicon for implementing a multiparametric gas sensor, which would enable simultaneous monitoring of electrical and optical parameters, is investigated. In addition, different porous silicon-based optical filter-structures are prepared, and their properties in sensing applications are evaluated and compared. First and foremost, thermal carbonization is established as a viable method to stabilize porous silicon optical filters for chemical sensing applications. Furthermore, a multiparametric sensor, which can be used for increasing selectivity in gas sensing, is also demonstrated. Methods to improve spectral quality in multistopband mesoporous silicon rugate filters are studied, and structural effects to gas sorption kinetics are evaluated. Finally, the stability of thermally carbonized optical filters in basic environments is found to be superior in comparison to other surface chemistries currently available for porous silicon. The results presented in this thesis are of particular interest for developing novel reliable sensing systems based on porous silicon, e.g., label-free optical biosensors.

Understanding the in vitro dissolution rate of glasses with respect to future clinical applications

Glass is a unique material with a long history. Several glass products are used daily in our everyday life, often unnoticed. Glass can be found not only in obvious applications such as tableware, windows, and light bulbs, but also in tennis rackets, windmill turbine blades, optical devices, and medical implants. The glasses used at present as implants are inorganic silica-based melt-derived compositions mainly for hard-tissue repair as bone graft substitute in dentistry and orthopedics. The degree of glass reactivity desired varies according to implantation situation and it is vital that the ion release from any glasses used in medical applications is controlled. Understanding the in vitro dissolution rate of glasses provides a first approximation of their behavior in vivo. Specific studies concerning dissolution properties of bioactive glasses have been relatively scarce and mostly concentrated to static condition studies. The motivation behind this work was to develop a simple and accurate method for quantifying the in vitro dissolution rate of highly different types of glass compositions with interest for future clinical applications. By combining information from various experimental conditions, a better knowledge of glass dissolution and the suitability of different glasses for different medical applications can be obtained. Thus, two traditional and one novel approach were utilized in this thesis to study glass dissolution. The chemical durability of silicate glasses was tested in water and TRIS-buffered solution at static and dynamic conditions. The traditional in vitro testing with a TRISbuffered solution under static conditions works well with bioactive or with readily dissolving glasses, and it is easy to follow the ion dissolution reactions. However, in the buffered solution no marked differences between the more durable glasses were observed. The hydrolytic resistance of the glasses was studied using the standard procedure ISO 719. The relative scale given by the standard failed to provide any relevant information when bioactive glasses were studied. However, the clear differences in the hydrolytic resistance values imply that the method could be used as a rapid test to get an overall idea of the biodegradability of glasses. The standard method combined with the ion concentration and pH measurements gives a better estimate of the hydrolytic resistance because of the high silicon amount released from a glass. A sensitive on-line analysis method utilizing inductively coupled plasma optical emission spectrometer and a flow-through micro-volume pH electrode was developed to study the initial dissolution of biocompatible glasses. This approach was found suitable for compositions within a large range of chemical durability. With this approach, the initial dissolution of all ions could be measured simultaneously and quantitatively, which gave a good overall idea of the initial dissolution rates for the individual ions and the dissolution mechanism. These types of results with glass dissolution were presented for the first time during the course of writing this thesis. Based on the initial dissolution patterns obtained with the novel approach using TRIS, the experimental glasses could be divided into four distinct categories. The initial dissolution patterns of glasses correlated well with the anticipated bioactivity. Moreover, the normalized surface-specific mass loss rates and the different in vivo models and the actual in vivo data correlated well. The results suggest that this type of approach can be used for prescreening the suitability of novel glass compositions for future clinical applications. Furthermore, the results shed light on the possible bioactivity of glasses. An additional goal in this thesis was to gain insight into the phase changes occurring during various heat treatments of glasses with three selected compositions. Engineering-type T-T-T curves for glasses 1-98 and 13-93 were stablished. The information gained is essential in manufacturing amorphous porous implants or for drawing of continuous fibers of the glasses. Although both glasses can be hot worked to amorphous products at carefully controlled conditions, 1-98 showed one magnitude greater nucleation and crystal growth rate than 13-93. Thus, 13-93 is better suited than 1-98 for working processes which require long residence times at high temperatures. It was also shown that amorphous and partially crystalline porous implants can be sintered from bioactive glass S53P4. Surface crystallization of S53P4, forming Na2O∙CaO∙2SiO2, was observed to start at 650°C. The secondary crystals of Na2Ca4(PO4)2SiO4, reported for the first time in this thesis, were detected at higher temperatures, from 850°C to 1000°C. The crystal phases formed affected the dissolution behavior of the implants in simulated body fluid. This study opens up new possibilities for using S53P4 to manufacture various structures, while tailoring their bioactivity by controlling the proportions of the different phases. The results obtained in this thesis give valuable additional information and tools to the state of the art for designing glasses with respect to future clinical applications. With the knowledge gained we can identify different dissolution patters and use this information to improve the tuning of glass compositions. In addition, the novel online analysis approach provides an excellent opportunity to further enhance our knowledge of glass behavior in simulated body conditions.

Microstructural and environmental effects on the mechanical properties of porous particle-polymer composites

Papper kan anses vara ett av de mest använda materialen i det dagliga livet. Tidskrifter, tidningar, böcker och diverse förpackningar är några exempel på pappersbaserade produkter. Papperets egenskaper måste anpassas till användningsändamålet. En tidskrift kräver t.ex. hög ljushet, opacitet och en slät yta hos papperet, medan dessa egenskaper är mindre viktiga för en dagstidning. Allt tryckpapper behöver vissa mekaniska egenskaper för att tåla vidarebearbetning såsom kalandrering, tryckning och vikning. Man kan bestryka papper för att förbättra dess optiska egenskaper och tryckbarhetsegenskaper. Vid bestrykning appliceras en dispersion av mineralpigment och polymerbindemedel som ett tunt lager på papperets yta. Bestrykningsskiktet kan ses som ett komplext, poröst kompositmaterial som även bidrar till papperets mekaniska egenskaper och dess processerbarhet i diverse konverteringsoperationer. Kravet på framställning av förmånligt papper med tillräckliga styrkeegenskaper ställer allt högre krav på optimeringen av pappersbestrykningsskiktets egenskaper och produktionskostnader. Målet med detta arbete var att förstå sambandet mellan pigmentbestrykningsskiktets mikrostruktur och dess makroskopiska, mekaniska egenskaper. Resultaten visar att adhesionen i gränsytan mellan pigment och bindemedel är kritisk för bestrykningsskiktets förmåga att bära mekanisk belastning. Polära vätskor är vanliga i tryckfärger och kan, eftersom de påverkar syra/bas-interaktionerna mellan pigment och latexbindemedel, försvaga denna adhesion. Resultaten tyder på att ytstyrkan hos bestruket papper kan höjas genom användning av bifunktionella dispergeringsmedel för mineralpigment. Detta medför inbesparingar i pappersproduktionen eftersom mängden bindemedel, den dyraste komponenten i bestrykningsskiktet, kan minskas.

Development of porous glass-fiber reinforced composite for bone implants. Evaluation of antimicrobial effect and implant fixation

Cranial bone reconstructions are necessary for correcting large skull bone defects due to trauma, tumors, infections and craniotomies. Traditional synthetic implant materials include solid or mesh titanium, various plastics and ceramics. Recently, biostable glass-fiber reinforced composites (FRC), which are based on bifunctional methacrylate resin, were introduced as novel implant solution. FRCs were originally developed and clinically used in dental applications. As a result of further in vitro and in vivo testing, these composites were also approved for clinical use in cranial surgery. To date, reconstructions of large bone defects were performed in 35 patients. This thesis is dedicated to the development of a novel FRC-based implant for cranial reconstructions. The proposed multi-component implant consists of three main parts: (i) porous FRC structure; (ii) bioactive glass granules embedded between FRC layers and (iii) a silver-polysaccharide nanocomposite coating. The porosity of the FRC structure should allow bone ingrowth. Bioactive glass as an osteopromotive material is expected to stimulate the formation of new bone. The polysaccharide coating is expected to prevent bacterial colonization of the implant. The FRC implants developed in this study are based on the porous network of randomly-oriented E-glass fibers bound together by non-resorbable photopolymerizable methacrylate resin. These structures had a total porosity of 10–70 volume %, of which > 70% were open pores. The pore sizes > 100 μm were in the biologically-relevant range (50-400 μm), which is essential for vascularization and bone ingrowth. Bone ingrowth into these structures was simulated by imbedding of porous FRC specimens in gypsum. Results of push-out tests indicated the increase in the shear strength and fracture toughness of the interface with the increase in the total porosity of FRC specimens. The osteopromotive effect of bioactive glass is based on its dissolution in the physiological environment. Here, calcium and phosphate ions, released from the glass, precipitated on the glass surface and its proximity (the FRC) and formed bone-like apatite. The biomineralization of the FRC structure, due to the bioactive glass reactions, was studied in Simulated Body Fluid (SBF) in static and dynamic conditions. An antimicrobial, non-cytotoxic polysaccharide coating, containing silver nanoparticles, was obtained through strong electrostatic interactions with the surface of FRC. In in vitro conditions the lactose-modified chitosan (chitlac) coating showed no signs of degradation within seven days of exposure to lysozyme or one day to hydrogen peroxide (H2O2). The antimicrobial efficacy of the coating was tested against Staphylococcus aureus and Pseudomonas aeruginosa. The contact-active coating had an excellent short time antimicrobial effect. The coating neither affected the initial adhesion of microorganisms to the implant surface nor the biofilm formation after 24 h and 72 h of incubation. Silver ions released to the aqueous environment led to a reduction of bacterial growth in the culture medium.

Analyzing hydraulic head of fluid flow in underground porous medium

Hydraulic head is distributed through a medium with porous aspect. The analysis of hydraulic head from one point to another is used by the Richard's equation. This equation is equivalent to the groundwater ow equation that predicts the volumetric water contents. COMSOL 3.5 is used for computation applying Richard's equation. A rectangle of 100 meters of length and 10 meters of large (depth) with 0,1 m/s fl ux of inlet as source of our fl uid is simulated. The domain have Richards' equation model in two dimension (2D). Hydraulic head increases proportional with moisture content.

Porous Bodies, Porous Minds. Emotions and the Supernatural in the Íslendingasögur (ca. 1200–1400)

In my PhD Thesis, I study the conceptions and representation of emotions in medieval 13th and 14th-century Iceland. I have used Icelandic saga literature as my source material and Icelandic Family sagas (Íslendingasögur) as my main sources. Firstly, I wished to explore in my study the medieval Icelandic folk theory of emotions: what emotions were thought to be, from what they originated and how they operated? Secondly, in earlier research it has been shown that emotions were seldom described in Íslendingasögur. They were mostly represented in dialogue, poetry or in somatic changes (e.g. turning pale). Consequently, I examined whether medieval Icelanders had alternative emotion discourses in literature, in addition to the usual manner of representation. My study consists of qualitative case studies, and I have analysed the sources intertextually. I suggest that medieval Icelanders regarded emotions as movements of the mind. The mind existed in the heart. As a consequence, emotions were considered physical in nature. The human body and therefore also the human mind was considered porous: if the mind of the person was not strong enough, supernatural agents and forces could penetrate theboundaries of his/her body as winds or sharp projectiles. Correspondingly, minds of strong-willed people could penetrate the minds of others. As a result, illness and emotions could upspring. People did not always distinguish between emotions and physical illnesses. Excessive emotions could cause illness, even death. Especially fear, grief and emotions of moral responsibility (e.g. guilt) made people vulnerable to the supernatural influence. Guilt was considered part of the emotional experience of misfortune (ógæfa), and in literature guilt could also be represented as eye pain that was inflicted upon the sufferer by a supernatural agent in a dream. Consequently, supernatural forces and beings were part of the upspring of emotions, but also part of the representation of emotions in literature: They caused the emotion but their presence also represented the emotional turmoil in the lives of the people that the supernatural agents harassed; emotions that had followed from norm transgressions, betrayal and other forms of social disequilibrium. Medieval readers and listeners of the Íslendingasögur were used to interpreting such different layers of meaning in texts.

Superkondensaattoreiden valmistus selluloosageeleistä ja hiilimustasta

Tässä työssä tavoitteena oli selvittää selluloosageelien soveltuvuutta superkondensaattorien elektrodimateriaaliksi. Superkondensaattoreiden elektrodit valmistetaan yleensä hyvin huokoisesta hiilimateriaalista. Näin ollen selluloosageeliin joukkoon sekoitettiin sähkönjohtavuutta parantamaan hiilimustaa. Kirjallisessa osassa keskityttiin superkondensaattorin toimintaperiaatteeseen, rakenteeseen ja yleisimmin käytössä oleviin elektrodimateriaaleihin. Yksityiskohtaisemmin paneuduttiin selluloosan käyttöön superkondensaattorien rakennemateriaalina. Kokeellisessa osassa valmistettiin erilaisia selluloosageelejä ja tämän jälkeen elektrodeja hiilimusta-selluloosageeli-seoksesta. Eri sekoitustapojen vaikutusta elektrodin sähköisiin ominaisuuksiin tarkasteltiin sekoittamalla hiilimustaa selluloosageeliin kahdella eri tavalla; korkeapainehajottajalla ja sekoittimella. Tämän lisäksi selvitettiin myös eri kuivaustapojen ja kuivauslämpötilojen vaikutusta elektrodiarkkien ominaisuuksiin. Arkkien sähköiset ominaisuudet määritettiin syklistä voltammetriaa (SV) käyttäen. Saavutetut tulokset osoittivat, että selluloosageelin valmistus on yksi kriittisimmistä kohdista elektrodimateriaalin valmistuksessa. Riittävä entsyymikäsittely vaaditaan, jotta saadaan muodostettua sopiva geeli, johon hiilimusta-partikkelit saadaan takertumaan. Hiilimustan sekoitustavalla sekä elektrodiarkin kuivauslämpötilalla ja – tavalla on myös suuresti vaikutusta elektrodin sähköisiin ominaisuuksiin.

Atomic Force Microscopy Investigation of NI Particles Deposited on Porous Silicon

In this thesis properties and influence of modification techniques of porous silicon were studied by Atomic Force Microscope (AFM). This device permits to visualize the surface topography and to study properties of the samples on atomic scale, which was necessary for recent investigation. Samples of porous silicon were obtained by electrochemical etching. Nickel particles were deposited by two methods: electrochemical deposition and extracting from NiCl2 ethanol solution. Sample growth was conducted in Saint-Petersburg State Electrotechnical University, LETI. Kelvin probe force microscopy (KPFM) and Magnetic force microscopy (MFM) were utilized for detailed information about surface properties of the samples. Measurements showed the difference in morphology correlating with initial growth conditions. Submicron size particles were clearly visible on surfaces of the treated samples. Although their nature was not clarified due to limitations of AFM technique. It is expected that surfaces were covered by nanometer scale Ni particles, which can be verified by implication of RAMAN device.

Flow in porous media of oil-water systems

Fluid flow behaviour in porous media is a conundrum. Therefore, this research is focused on filtration-volumetric characterisation of fractured-carbonate sediments, coupled with their proper simulation. For this reason, at laboratory rock properties such as pore volume, permeability and porosity are measured, later phase permeabilities and oil recovery in function of flow rate are assessed. Furthermore, the rheological properties of three oils are measured and analysed. Finally based on rock and fluid properties, a model using COMSOL Multiphysics is built in order to compare the experimental and simulated results. The rock analyses show linear relation between flow rate and differential pressure, from which phase permeabilities and pressure gradient are determined, eventually the oil recovery under low and high flow rate is established. In addition, the oils reveal thixotropic properties as well as non-Newtonian behaviour described by Bingham model, consequently Carreau viscosity model for the used oil is given. Given these points, the model for oil and water is built in COMSOL Multiphysics, whereupon successfully the reciprocity between experimental and simulated results is analysed and compared. Finally, a two-phase displacement model is elaborated.

Coolability of porous core debris beds : Effects of bed geometry and multi-dimensional flooding

This thesis addresses the coolability of porous debris beds in the context of severe accident management of nuclear power reactors. In a hypothetical severe accident at a Nordic-type boiling water reactor, the lower drywell of the containment is flooded, for the purpose of cooling the core melt discharged from the reactor pressure vessel in a water pool. The melt is fragmented and solidified in the pool, ultimately forming a porous debris bed that generates decay heat. The properties of the bed determine the limiting value for the heat flux that can be removed from the debris to the surrounding water without the risk of re-melting. The coolability of porous debris beds has been investigated experimentally by measuring the dryout power in electrically heated test beds that have different geometries. The geometries represent the debris bed shapes that may form in an accident scenario. The focus is especially on heap-like, realistic geometries which facilitate the multi-dimensional infiltration (flooding) of coolant into the bed. Spherical and irregular particles have been used to simulate the debris. The experiments have been modeled using 2D and 3D simulation codes applicable to fluid flow and heat transfer in porous media. Based on the experimental and simulation results, an interpretation of the dryout behavior in complex debris bed geometries is presented, and the validity of the codes and models for dryout predictions is evaluated. According to the experimental and simulation results, the coolability of the debris bed depends on both the flooding mode and the height of the bed. In the experiments, it was found that multi-dimensional flooding increases the dryout heat flux and coolability in a heap-shaped debris bed by 47–58% compared to the dryout heat flux of a classical, top-flooded bed of the same height. However, heap-like beds are higher than flat, top-flooded beds, which results in the formation of larger steam flux at the top of the bed. This counteracts the effect of the multi-dimensional flooding. Based on the measured dryout heat fluxes, the maximum height of a heap-like bed can only be about 1.5 times the height of a top-flooded, cylindrical bed in order to preserve the direct benefit from the multi-dimensional flooding. In addition, studies were conducted to evaluate the hydrodynamically representative effective particle diameter, which is applied in simulation models to describe debris beds that consist of irregular particles with considerable size variation. The results suggest that the effective diameter is small, closest to the mean diameter based on the number or length of particles.