Heat flux measurement inside internal combustion engine with gradient heat flux sensor.

This master’s thesis is devoted to study different heat flux measurement techniques such as differential temperature sensors, semi-infinite surface temperature methods, calorimetric sensors and gradient heat flux sensors. The possibility to use Gradient Heat Flux Sensors (GHFS) to measure heat flux in the combustion chamber of compression ignited reciprocating internal combustion engines was considered in more detail. A. Mityakov conducted an experiment, where Gradient Heat Flux Sensor was placed in four stroke diesel engine Indenor XL4D to measure heat flux in the combustion chamber. The results which were obtained from the experiment were compared with model’s numerical output. This model (a one – dimensional single zone model) was implemented with help of MathCAD and the result of this implementation is graph of heat flux in combustion chamber in relation to the crank angle. The values of heat flux throughout the cycle obtained with aid of heat flux sensor and theoretically were sufficiently similar, but not identical. Such deviation is rather common for this type of experiment.

Design of an amplifier and an instrumentation setup for a harsh environment heat flux measurement application – case axial flux electric machine

Measurement is a tool for researching. Therefore, it is important that the measuring process is carried out correctly, without distorting the signal or the measured event. Researches of thermoelectric phenomena have been focused more on transverse thermoelectric phenomena during recent decades. Transverse Seebeck effect enables to produce thinner and faster heat flux sensor than before. Studies about transverse Seebeck effect have so far focused on materials, so in this Master’s Thesis instrumentation of transverse Seebeck effect based heat flux sensor is studied, This Master’s Thesis examines an equivalent circuit of transverse Seebeck effect heat flux sensors, their connectivity to electronics and choosing and design a right type amplifier. The research is carried out with a case study which is Gradient Heat Flux Sensors and an electrical motor. In this work, a general equivalent circuit was presented for the transverse Seebeck effect-based heat flux sensor. An amplifier was designed for the sensor of the case study, and the solution was produced for the measurement of the local heat flux of the electric motor to improve the electromagnetic compatibility.

Automatization of feeding the heat flux data to water circulation simulations of recovery boilers

Tässä työssä perehdytään soodakattiloiden vesikiertomallin rakentamiseen. Työn päätavoitteena on kehittää simulointimallia varten taulukkolaskentapohja, jonka avulla soodakattilan lämpövuotietoja on yksinkertaista ja nopeaa käsitellä ja siirtää Apros 6 -simulointiohjelmaan. Lisäksi tarkoituksena on pyrkiä automatisoimaan työvaiheet mahdollisimman pitkälle, jolloin vesikiertolaskennan tekeminen yksinkertaistuisi, yhtenäistyisi ja tarkentuisi. Tämä on mahdollista Excel- makrojen ja Apros 6:n uusien toimintojen avulla. Apros 6:ssa on nyt mahdollista hyödyntää SCL- komentotiedostoja, joiden avulla sujuva tiedonsiirto Aproksen ja Excelin välillä vodaan toteuttaa. Vesikiertolaskentaan käytettävän datan käsittely on aikaisemmin ollut työlästä ja sen tarkkuus on pitkälti riippunut mallintajasta. Tässä diplomityössä päästään hyödyntämään uusimpia ja realistisempia soodakattiloiden CFD- malleja, joiden avulla pystytään luomaan aikaisempaa tarkemmat lämpövuojakaumat soodakattilan lämpöpinnoille. Tämä muutos parantaa vesikiertolaskennan tarkkuutta. Työn kokeellisessa osassa uutta Excel laskentatyökalua ja uusia lämpövuoarvoja testataan käytännössä. Eräs vanha Apros- vesikiertomalli päivitetään uusilla lämpövuoarvoilla ja sen rakenteeseen tehdään muutoksia tarkkuuden parantamiseksi. Uuden mallin toimivuutta testataan myös 115 %:n kapasiteetilla ja tutkitaan kuinka kyseinen vesikiertopiiri reagoi suurempaan lämpötehoon. Näitä kolmea eri tilannetta vertaillaan toisiinsa ja tarkastellaan eroavaisuuksia niiden vesi-höyrypiireissä.

A three dimensional study of heat transfer in a supercritical circulating fluidized bed combustor

From the boiler design point of view, it is imperative to know and understand the operation of the boiler. Since comprehensive measurement of a large furnace is impossible, the furnace can be modeled in order to study its behavior and phenomena. This requires the used model to be validated to correspond with the physical furnace behavior. In this thesis, a three dimensional furnace model is validated to match a bituminous coal utilizing, supercritical once-through circulating fluidized bed combustor based on measurement data. The validated model is used for analyzing the furnace heat transfer. Other heat transfer analysis methods are energy balance method based on tube surface temperature measurements and a method based on measured temperature difference between the tube crest and the fin. The latter method was developed in the thesis using Fluent-software. In the theory part, literature is reviewed and the fundamental aspects of circulating fluidized bed are discussed. These aspects are solid particle behavior in fluidization known as hydrodynamics, behavior of fuel and combustion and heat transfer. Fundamental aspects of modeling are also presented.

Heat transfer inside internal combustion engine: modelling and comparison with experimental data

The paper is devoted to study specific aspects of heat transfer in the combustion chamber of compression ignited reciprocating internal combustion engines and possibility to directly measure the heat flux by means of Gradient Heat Flux Sensors (GHFS). A one – dimensional single zone model proposed by Kyung Tae Yun et al. and implemented with the aid of Matlab, was used to obtain approximate picture of heat flux behavior in the combustion chamber with relation to the crank angle. The model’s numerical output was compared to the experimental results. The experiment was accomplished by A. Mityakov at four stroke diesel engine Indenor XL4D. Local heat fluxes on the surface of cylinder head were measured with fast – response, high – sensitive GHFS. The comparison of numerical data with experimental results has revealed a small deviation in obtained heat flux values throughout the cycle and different behavior of heat flux curve after Top Dead Center.

Determination of the thermal parameters of high-power batteries by local heat

A new approach to the determination of the thermal parameters of high-power batteries is introduced here. Application of local heat flux measurement with a gradient heat flux sensor (GHFS) allows determination of the cell thermal parameters in di_erent surface points of the cell. The suggested methodology is not cell destructive as it does not require deep discharge of the cell or application of any charge/discharge cycles during measurements of the thermal parameters of the cell. The complete procedure is demonstrated on a high-power Li-ion pouch cell, and it is verified on a sample with well-known thermal parameters. A comparison of the experimental results with conventional thermal characterization methods shows an acceptably low error. The dependence of the cell thermal parameters on state of charge (SoC) and measurement points on the surface was studied by the proposed measurement approach.

Kerrostumien muodostumisen määritys sondin lämmönsiirron perusteella

Diplomityön tavoitteena oli tutkia kerrostumien muodostumista voimalaitoskattiloiden lämpöpinnoille ja niiden vaikutusta lämmönsiirtoon. Kerrostumien vaikutusta lämpöpintojen lämpötiloihin tutkittiin kerrostumasondin avulla kierto- ja kerrosleijukattiloissa. Saadun mittausaineiston muodostettiin mittaustilannetta vastaava laskentamalli, jonka avulla pystytään selvittämään voimalaitoksen lämpöpintojen likaantumista käyttöperiodin aikana. Työn alussa tarkasteltiin kerrostumien muodostumismekanismeja lämmönsiirtopinnoille, kerrostumien ominaisuuksia ja virtausympäristön vaikutusta lämmönsiirtoon. Poltossa syntyvä tuhka kiinnittyy eri mekanismeilla kattilan seinille ja lämpöpinnoille riippuen polttoaineesta ja ympäröivistä olosuhteista. Likakerrosten muodostuminen lämmönsiirtimen pinnalle pienentää putkeen siirtyvää lämpövirtaa ja alentaa kattilan hyötysuhdetta. Kerrostumien lämmönjohtumisesta ja mikrorakenteesta on hyvin vähän kokeellisen tietoa, mikä vaikeuttaa reaalisten mallien muodostamista kattilan käyttäytymisestä sen likaantuessa.Työn kokeellisessa osassa tarkasteltiin kahdella eri kattilalla tehtyjä mittauksia. Mittaukset tehtiin kerrostumasondilla, jonka lämmönsiirrosta luotiin laskentamalli SIMULINK-simulointiohjelmalla. Mittaustuloksina saatiin kolmen eri pinnan lämpötilat, jotka muuttuivat kerrostuman ja jäähdytyksen vaikutuksesta. Laskentamallista muokattiin mittaustilanteita vastaava, jolloin lämpötilamuutoksista nähdään likakerroksen ominaisuuksien vaikutus lämmönsiirtoon. Sondista muodostettiin myös FLUENT-malli, jolla tarkasteltiin yksittäisen putken virtauskenttää sekä kahden lähekkäin olevan putken vaikutusta virtaukseen.

Muutosilmiöt soodakattilakeossa

Soodakattilan sulakeon epästationaarinen käyttäytyminen sekä keon pitkä jäähtymisaika alasajon jälkeen ovat aiheuttaneet ongelmia kattilan taloudellisessa käytettävyydessä. Keon käyttäytymisestä on luotu CFD-malleja, joiden tavoitteena on havainnollistaa keon lämpötilajakaumaa ja rakennetta. Mallien ongelmana on se, että niissä huomioidaan vain keon aktiivinen pintakerros. Keon sisäosan rakennetta ja siinä tapahtuvia prosesseja ei toistaiseksi tunneta kunnolla luotettavan, koko keon kattavan mallin luomiseen. Tässä työssä tutkittiin sulakeon käytön aikana havaittuja muutosilmiöitä, jotka vaikuttavat keon rakenteeseen ja ominaisuuksiin sekä tutkittiin ilmiöiden taustalla olevia tekijöitä. Näitä tekijöitä ovat keon sisässä tapahtuvat kemialliset ja fyysiset prosessit, jotka aiheuttavat muutoksia niin lämpöteknisesti kuin fyysisesti sekä ulkoapäin tulevat tekijät, jotka aiheutuvat ajotilanteiden seurauksena tapahtuvista muutoksista. Työn kokeellisena osana luotiin sulakeon jäähtymismalli käyttäen 1-dimensionaalista ADL-mallia. Mallin pohjana käytettiin StoraEnso Oy:ltä Oulun soodakattilan sulakeosta saatua mittausraporttia. ADL-mallin avulla luotiin keon jäähtymiskäyrät lämpötilan ja syvyyden funktiona. Saadut käyrät täsmäsivät hyvin mittausraportin tuloksiin. Mallin avulla keolle saatiin muodostettua energiatase, jonka tuloksena keosta 12 tunnin aikana poistuva lämpövirta pinnalla oli noin 9.8kW/m2 ja pinnan lämmönsiirtokerroin 58.3W/m²°C. Pohjan poistuvaksi lämpövirraksi saatiin 14.1kW/m2 ja lämmönsiirtokertoimeksi 75.4W/m²°C. Termiseksi diffuusiokertoimeksi saatiin 3.9•10-7m²/s.

Statistical Analysis and Numerics of Heat Exchanger Models

The identifiability of the parameters of a heat exchanger model without phase change was studied in this Master’s thesis using synthetically made data. A fast, two-step Markov chain Monte Carlo method (MCMC) was tested with a couple of case studies and a heat exchanger model. The two-step MCMC-method worked well and decreased the computation time compared to the traditional MCMC-method. The effect of measurement accuracy of certain control variables to the identifiability of parameters was also studied. The accuracy used did not seem to have a remarkable effect to the identifiability of parameters. The use of the posterior distribution of parameters in different heat exchanger geometries was studied. It would be computationally most efficient to use the same posterior distribution among different geometries in the optimisation of heat exchanger networks. According to the results, this was possible in the case when the frontal surface areas were the same among different geometries. In the other cases the same posterior distribution can be used for optimisation too, but that will give a wider predictive distribution as a result. For condensing surface heat exchangers the numerical stability of the simulation model was studied. As a result, a stable algorithm was developed.

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.

Arinakattilan lämmönsiirron laskentamallin kehittäminen

Kattilalaitosten polttoaineen syötössä ilmenevät häiriöt ja biopolttoaineiden laatuvaihtelut aiheuttavat epävakaata palamista ja tekevät prosessin hallinnasta vaikeampaa. Polttoaineen laatuvaihtelut vaikuttavat koko prosessiin ja näkyvät lopulta myös höyryntuotannossa. Kompensoinnilla pyritään estämään häiriöiden suuret vaikutukset höyryn tuotantoon. Tarkoituksena on saada kattilan toiminta ja tehon tuotanto tasaisemmaksi ja helpommin hallittavaksi. Tässä diplomityössä tarkastellaan polttoaineen ominaisuuksien, säteilylämmönsiirron sekä säätöjen vaikutusta toisiinsa ja merkitystä kattilan toiminnan kannalta. Työssä muodostetaan säteilylämmönsiirron laskentamalli arinakattilan tulipesälle käyttäen hyväksi hyvin sekoittuneen tulipesän menetelmää. Menetelmällä voidaan määrittääsavukaasujen keskimääräinen lämpötila tulipesässä, lämpövirta tulipesän seiniin tai poltossa vapautuva lämpöteho. Mallin avulla voidaan paremmin ymmärtää prosessin käyttäytymistä polttoaineen laadun muuttuessa sekä helpottaa ja nopeuttaa kattilan käyttäytymisen ennustamista. Laskentamalli tehtiin Excel –laskentaohjelmaan, jossa se testataan. Verifioinnin jälkeen malli on tarkoitus siirtää toimimaan apros –simulointiympäristöön.

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.

Modeling of residual stresses and distortion due to welding in fillet welds

Welding has a growing role in modern world manufacturing. Welding joints are extensively used from pipes to aerospace industries. Prediction of welding residual stresses and distortions is necessary for accurate evaluation of fillet welds in relation to design and safety conditions. Residual stresses may be beneficial or detrimental, depending whether they are tensile or compressive and the loading. They directly affect the fatigue life of the weld by impacting crack growth rate. Beside theoretical background of residual stresses this study calculates residual stresses and deformations due to localized heating by welding process and subsequent rapid cooling in fillet welds. Validated methods are required for this purpose due to complexity of process, localized heating, temperature dependence of material properties and heat source. In this research both empirical and simulation methods were used for the analysis of welded joints. Finite element simulation has become a popular tool of prediction of welding residual stresses and distortion. Three different cases with and without preload have been modeled during this study. Thermal heat load set is used by calculating heat flux from the given heat input energy. First the linear and then nonlinear material behavior model is modeled for calculation of residual stresses. Experimental work is done to calculate the stresses empirically. The results from both the methods are compared to check their reliability. Residual stresses can have a significant effect on fatigue performance of the welded joints made of high strength steel. Both initial residual stress state and subsequent residual stress relaxation need to be considered for accurate description of fatigue behavior. Tensile residual stresses are detrimental and will reduce the fatigue life and compressive residual stresses will increase it. The residual stresses follow the yield strength of base or filler material and the components made of high strength steel are typically thin, where the role of distortion is emphasizing.

Two-phase flow measurements with spacer grid in a rod bundle geometry

Heat transfer effectiveness in nuclear rod bundles is of great importance to nuclear reactor safety and economics. An important design parameter is the Critical Heat Flux (CHF), which limits the transferred heat from the fuel to the coolant. The CHF is determined by flow behaviour, especially the turbulence created inside the fuel rod bundle. Adiabatic experiments can be used to characterize the flow behaviour separately from the heat transfer phenomena in diabatic flow. To enhance the turbulence, mixing vanes are attached to spacer grids, which hold the rods in place. The vanes either make the flow swirl around a single sub-channel or induce cross-mixing between adjacent sub-channels. In adiabatic two-phase conditions an important phenomenon that can be investigated is the effect of the spacer on canceling the lift force, which collects the small bubbles to the rod surfaces leading to decreased CHF in diabatic conditions and thus limits the reactor power. Computational Fluid Dynamics (CFD) can be used to simulate the flow numerically and to test how different spacer configurations affect the flow. Experimental data is needed to validate and verify the used CFD models. Especially the modeling of turbulence is challenging even for single-phase flow inside the complex sub-channel geometry. In two-phase flow other factors such as bubble dynamics further complicate the modeling. To investigate the spacer grid effect on two-phase flow, and to provide further experimental data for CFD validation, a series of experiments was run on an adiabatic sub-channel flow loop using a duct-type spacer grid with different configurations. Utilizing the wire-mesh sensor technology, the facility gives high resolution experimental data in both time and space. The experimental results indicate that the duct-type spacer grid is less effective in canceling the lift force effect than the egg-crate type spacer tested earlier.