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.

Kuplakokojakauman hallinta happidelignifioinnissa

Työn tarkoituksena oli uutta kuvantamistekniikkaa hyödyntäen tutkia erilaisten tekijöiden vaikutusta kaasun dispergoitumiseen kemikaalisekoittimessa, kun kaasua sekoitetaan keskisakeaan massaan. Lisäksi työssä pyrittiin selvittämään, kuinka paljon kuitususpensioon tuotettu kaasufaasin kuplakokojakauma vaikuttaa happidelignifioinnin tulokseen. Kaasumaisten aineiden käyttäytymistä keskisakeissa kuitususpensiossa ei tarkkaan tunneta. Mikäli kaasumaisen hapen käyttäytymisestä saadaan uutta tietoa, tarjoaa tämä muun muassa uusia mahdollisuuksia kaasua sekoittavien laitteiden tuotekehityksessä. Työn kokeellinen osuus koostui kahdesta osasta, joista ensimmäisessä osassa selvitettiin sekoittimen roottorin pyörimisnopeuden, reaktorin kaasutilavuuden sekä suspension sakeuden vaikutusta muodostuvaan kaasun kuplakokojakaumaan. Työn jälkimmäisessä osassa arvioitiin yksivaiheisten keskisakeudessa tehtyjen happidelignifiointien perusteella suspensioon tuotetun kaasun kuplakokojakauman merkitystä happidelignifiointitulokseen. Kuplakokojakaumat määritettiin reaktoriin kiinnitetyllä kameralla kuvatuista valokuvista, joita otettiin sekoitustapahtuman aikana. Työn tuloksien perusteella sekoituksen voimakkuudella oli suurin vaikutus suspensioon muodostuvan kuplakokojakauman kannalta. Roottorin kierrosnopeuden kasvaessa kaasun keskimääräinen kuplakoko pieneni sekä havaittujen kuplien lukumäärää kasvoi huomattavasti. Myös suspension sakeuden kasvattamisen havaittiin vaikuttavan kuplakokoon pienentävästi. Happidelignifioinneissa saavutettiin paras kappareduktio, kun kaasun kuplakoko oli mahdollisimman pieni. Käytetty kuvantamistekniikka on tiettävästi ensimmäinen menetelmä, jolla saadaan reaaliaikaista tietoa vain muutamien kymmenien mikrometrien kokoisten kaasukuplien käyttäytymisestä oikeassa prosessitilanteessa.

Modeling of water/steam circulation in Circulating Fluidized Bed boiler

A distinctive design feature of steam boiler with natural circulation is the presence of the steam drum which plays a role of the separator of vapor from the flow of water-and-steam mixture coming into steam drum from the furnace tubes. Steam drum with unheated downcomer tubes, deducing from it, and riser (screen/furnace tubes) inside the furnace is a closed circulation loop in which movement of water (downcomer tubes) and water-and-steam mixture (riser tubes) is organized. The movement of the working fluid is appears due to occurrence of the natural pressure, determined by the difference in hydrostatic pressure and the mass of water and water-and-steam mixtures in downcomer and riser tubes and called the driving pressure of the natural circulation: S drive = H steam (ρ down + ρ mix) g where: ρ down - density of water in downcomer tubes; ρ mix - density of water in riser tubes; H steam - height of steam content section; g - acceleration of gravity. In steam boilers with natural circulation the circulation rate is usually between 10 and 30. Thus, consumption of water in the circulation circuit “circulation rate times” more than steam output of the boiler. There are two aspects of the design of natural water circulation loops. One is to ensure a sufficient mass flux of circulating water to avoid burnout of evaporator tubes. The other is to avoid tube wall temperature fluctuation and tube vibration due to oscillation of circulation velocity. The design criteria are therefore reduced, in principle, to those of critical heat flux, critical flow rate for burnout, and flow instability. In practical design, however, the circulation velocity and the void fraction at the evaporator tube outlet are used as the design criteria (Seikan I., et. al., 1999). This study has been made with assumption that the heat flux in the furnace of the boiler even all the time. The target of the study was to define the circulation rate of the boiler, thus average heat flux do not change it. I would like to acknowledge professionals from “Foster Wheeler” company for good and comfortable cooperation.

Aukko-osuuden mittausmenetelmien kartoitus PWR PACTEL -koelaitteiston pystyhöyrystimen sekundääripuolella

Diplomityössä kartoitetaan mahdollisuuksia aukko-osuuden mittaamiseksi ydinvoimalaitosta mallintavan PWR PACTEL -koelaitteiston pystyhöyrystimen sekundääripuolella. Työ on toteutettu osana kansallista SAFIR2014-ydinturvallisuustutkimusohjelmaa. Diplomityön tavoitteena on löytää kustannuksiltaan mahdollisimman järkevä ja toimiva menetelmä aukko-osuuden määrittämiseksi. Aukko-osuuden mittaaminen on tärkeää sekundääripuolen kaksifaasivirtauksen käyttäytymisen paremman tuntemuksen lisäämiseksi. Aukko-osuusmittausdataa tarvitaan muun muassa laskentakoodien validointiin. Diplomityössä perehdytään kaksifaasivirtauksen ja aukko-osuuden fysiikkaan sekä esitellään erityyppisiä aukko-osuuden mittausmenetelmiä. Kunkin mittausmenetelmän soveltuvuutta PWR PACTEL -koelaitteistoon arvioidaan erikseen. Aukko-osuuden mittaaminen höyrystimen sekundääripuolella osoittautuu käytännössä erittäin hankalaksi. Pääasiassa tämä johtuu höyrystimen rakenteesta sekä mittausmenetelmien korkeista kustannuksista. Tämän vuoksi työssä tarkastellaan myös edellytyksiä aukko-osuuden mittaamiselle erillisessä höyrystintä mallintavassa koelaitteistossa. Mikäli aukko-osuutta haluttaisiin mitata erilliskoelaitteistossa, tulisi höyrystinmallin rakennetta, materiaaleja tai kiertoainetta muuttaa PWR PACTELin höyrystimeen verrattuna.

A Practical Micro Mechanical Model Based Local Approach Methodology for the Analysis of Ductile Fracture of Welded T-Joints

This thesis concentrates on developing a practical local approach methodology based on micro mechanical models for the analysis of ductile fracture of welded joints. Two major problems involved in the local approach, namely the dilational constitutive relation reflecting the softening behaviour of material, and the failure criterion associated with the constitutive equation, have been studied in detail. Firstly, considerable efforts were made on the numerical integration and computer implementation for the non trivial dilational Gurson Tvergaard model. Considering the weaknesses of the widely used Euler forward integration algorithms, a family of generalized mid point algorithms is proposed for the Gurson Tvergaard model. Correspondingly, based on the decomposition of stresses into hydrostatic and deviatoric parts, an explicit seven parameter expression for the consistent tangent moduli of the algorithms is presented. This explicit formula avoids any matrix inversion during numerical iteration and thus greatly facilitates the computer implementation of the algorithms and increase the efficiency of the code. The accuracy of the proposed algorithms and other conventional algorithms has been assessed in a systematic manner in order to highlight the best algorithm for this study. The accurate and efficient performance of present finite element implementation of the proposed algorithms has been demonstrated by various numerical examples. It has been found that the true mid point algorithm (a = 0.5) is the most accurate one when the deviatoric strain increment is radial to the yield surface and it is very important to use the consistent tangent moduli in the Newton iteration procedure. Secondly, an assessment of the consistency of current local failure criteria for ductile fracture, the critical void growth criterion, the constant critical void volume fraction criterion and Thomason's plastic limit load failure criterion, has been made. Significant differences in the predictions of ductility by the three criteria were found. By assuming the void grows spherically and using the void volume fraction from the Gurson Tvergaard model to calculate the current void matrix geometry, Thomason's failure criterion has been modified and a new failure criterion for the Gurson Tvergaard model is presented. Comparison with Koplik and Needleman's finite element results shows that the new failure criterion is fairly accurate indeed. A novel feature of the new failure criterion is that a mechanism for void coalescence is incorporated into the constitutive model. Hence the material failure is a natural result of the development of macroscopic plastic flow and the microscopic internal necking mechanism. By the new failure criterion, the critical void volume fraction is not a material constant and the initial void volume fraction and/or void nucleation parameters essentially control the material failure. This feature is very desirable and makes the numerical calibration of void nucleation parameters(s) possible and physically sound. Thirdly, a local approach methodology based on the above two major contributions has been built up in ABAQUS via the user material subroutine UMAT and applied to welded T joints. By using the void nucleation parameters calibrated from simple smooth and notched specimens, it was found that the fracture behaviour of the welded T joints can be well predicted using present methodology. This application has shown how the damage parameters of both base material and heat affected zone (HAZ) material can be obtained in a step by step manner and how useful and capable the local approach methodology is in the analysis of fracture behaviour and crack development as well as structural integrity assessment of practical problems where non homogeneous materials are involved. Finally, a procedure for the possible engineering application of the present methodology is suggested and discussed.

Libyan breast cancer: Health services and biology. Diagnosis delay and prognostic value of DNA pploidy, S-phase fraction, and Ki-67 and Bcl-2 immunohistochemistry

The aim of this study was to investigate the diagnosis delay and its impact on the stage of disease. The study also evaluated a nuclear DNA content, immunohistochemical expression of Ki-67 and bcl-2, and the correlation of these biological features with the clinicopathological features and patient outcome. 200 Libyan women, diagnosed during 2008–2009 were interviewed about the period from the first symptoms to the final histological diagnosis of breast cancer. Also retrospective preclinical and clinical data were collected from medical records on a form (questionnaire) in association with the interview. Tumor material of the patients was collected and nuclear DNA content analysed using DNA image cytometry. The expression of Ki-67 and bcl-2 were assessed using immunohistochemistry (IHC). The studies described in this thesis show that the median of diagnosis time for women with breast cancer was 7.5 months and 56% of patients were diagnosed within a period longer than 6 months. Inappropriate reassurance that the lump was benign was an important reason for prolongation of the diagnosis time. Diagnosis delay was also associated with initial breast symptom(s) that did not include a lump, old age, illiteracy, and history of benign fibrocystic disease. The patients who showed diagnosis delay had bigger tumour size (p<0.0001), positive lymph nodes (p<0.0001), and high incidence of late clinical stages (p<0.0001). Biologically, 82.7% of tumors were aneuploid and 17.3% were diploid. The median SPF of tumors was 11% while the median positivity of Ki-67 was 27.5%. High Ki-67 expression was found in 76% of patients, and high SPF values in 56% of patients. Positive bcl-2 expression was found in 62.4% of tumors. 72.2% of the bcl-2 positive samples were ER-positive. Patients who had tumor with DNA aneuploidy, high proliferative activity and negative bcl-2 expression were associated with a high grade of malignancy and short survival. The SPF value is useful cell proliferation marker in assessing prognosis, and the decision cut point of 11% for SPF in the Libyan material was clearly significant (p<0.0001). Bcl-2 is a powerful prognosticator and an independent predictor of breast cancer outcome in the Libyan material (p<0.0001). Libyan breast cancer was investigated in these studies from two different aspects: health services and biology. The results show that diagnosis delay is a very serious problem in Libya and is associated with complex interactions between many factors leading to advanced stages, and potentially to high mortality. Cytometric DNA variables, proliferative markers (Ki-67 and SPF), and oncoprotein bcl-2 negativity reflect the aggressive behavior of Libyan breast cancer and could be used with traditional factors to predict the outcome of individual patients, and to select appropriate therapy.