Characteristics of urban heat island (UHI) in a high latitude coastal city - a case study of Turku, SW Finland

The term urban heat island (UHI) refers to the common situation in which the city is warmer than its rural surroundings. In this dissertation, the local climate, and especially the UHI, of the coastal city of Turku (182,000 inh.), SW Finland, was studied in different spatial and temporal scales. The crucial aim was to sort out the urban, topographical and water body impact on temperatures at different seasons and times of the day. In addition, the impact of weather on spatiotemporal temperature differences was studied. The relative importance of environmental factors was estimated with different modelling approaches and a large number of explanatory variables with various spatial scales. The city centre is the warmest place in the Turku area. Temperature excess relative to the coldest sites, i.e. rural areas about 10 kilometers to the NE from the centre, is on average 2 °C. Occasionally, the UHI intensity can be even 10 °C. The UHI does not prevail continuously in the Turku area, but occasionally the city centre can be colder than its surroundings. Then the term urban cool island or urban cold island (UCI) is used. The UCI is most common in daytime in spring and in summer, whereas during winter the UHI prevails throughout the day. On average, the spatial temperature differences are largest in summer, whereas the single extreme values are often observed in winter. The seasonally varying sea temperature causes the shift of relatively warm areas towards the coast in autumn and inland in spring. In the long term, urban land use was concluded to be the most important factor causing spatial temperature differences in the Turku area. The impact was mainly a warming one. The impact of water bodies was emphasised in spring and autumn, when the water temperature was relatively cold and warm, respectively. The impact of topography was on average the weakest, and was seen mainly in proneness of relatively low-lying places for cold air drainage during night-time. During inversions, however, the impact of topography was emphasised, occasionally outperforming those of urban land use and water bodies.

Effect of heat input on the mechanical properties of welded joints in high strength steels

Experiments were carried out to determine the properties of the welded joints in 8mm thick high-strength steels produced by quenching and tempering and thermomechanical rolling with accelerated cooling (tensile strength 821–835 MPa). The dependence of the strength, elongation, hardness, impact energy and crack opening displacement on the heat input in the range 1.0–0.7 kJ mm21 was determined. The results show that the dependence of the strength of the welded joints decreases and that of the elongation increases. The heat input has only a slight effect on the impact energy and crack opening displacement in the heat-affected zone.

Effect of Noncondensable Gases on Circulation of Primary Coolant in Nuclear Power Plants in Abnormal Situations

The present study focuses on two effects of the presence of a noncondensable gas on the thermal-hydraulic behavior of thecoolant of the primary circuit of a nuclear reactor in the VVER-440 geometry inabnormal situations. First, steam condensation with the presence of air was studied in the horizontal tubes of the steam generator (SG) of the PACTEL test facility. The French thermal-hydraulic CATHARE code was used to study the heat transfer between the primary and secondary side in conditions derived from preliminary experiments performed by VTT using PACTEL. In natural circulation and single-phase vapor conditions, the injection of a volume of air, equivalent to the totalvolume of the primary side of the SG at the entrance of the hot collector, did not stop the heat transfer from the primary to the secondary side. The calculated results indicate that air is located in the second half-length (from the mid-length of the tubes to the cold collector) in all the tubes of the steam generator The hot collector remained full of steam during the transient. Secondly, the potential release of the nitrogen gas dissolved in the water of the accumulators of the emergency core coolant system of the Loviisa nuclear power plant (NPP) was investigated. The author implemented a model of the dissolution and release ofnitrogen gas in the CATHARE code; the model created by the CATHARE developers. In collaboration with VTT, an analytical experiment was performed with some components of PACTEL to determine, in particular, the value of the release time constant of the nitrogen gas in the depressurization conditions representative of the small and intermediate break transients postulated for the Loviisa NPP. Such transients, with simplified operating procedures, were calculated using the modified CATHARE code for various values of the release time constant used in the dissolution and release model. For the small breaks, nitrogen gas is trapped in thecollectors of the SGs in rather large proportions. There, the levels oscillate until the actuation of the low-pressure injection pumps (LPIS) that refill the primary circuit. In the case of the intermediate breaks, most of the nitrogen gas is expelled at the break and almost no nitrogen gas is trapped in the SGs. In comparison with the cases calculated without taking into account the release of nitrogen gas, the start of the LPIS is delayed by between 1 and 1.75 h. Applicability of the obtained results to the real safety conditions must take into accountthe real operating procedures used in the nuclear power plant.

Numerical and ExperimentalModelling of Gas Flow and Heat Transfer in the Air Gap of An Electric Machine

This work deals with the cooling of high-speed electric machines, such as motors and generators, through an air gap. It consists of numerical and experimental modelling of gas flow and heat transfer in an annular channel. Velocity and temperature profiles are modelled in the air gap of a high-speed testmachine. Local and mean heat transfer coefficients and total friction coefficients are attained for a smooth rotor-stator combination at a large velocity range. The aim is to solve the heat transfer numerically and experimentally. The FINFLO software, developed at Helsinki University of Technology, has been used in the flow solution, and the commercial IGG and Field view programs for the grid generation and post processing. The annular channel is discretized as a sector mesh. Calculation is performed with constant mass flow rate on six rotational speeds. The effect of turbulence is calculated using three turbulence models. The friction coefficient and velocity factor are attained via total friction power. The first part of experimental section consists of finding the proper sensors and calibrating them in a straight pipe. After preliminary tests, a RdF-sensor is glued on the walls of stator and rotor surfaces. Telemetry is needed to be able to measure the heat transfer coefficients at the rotor. The mean heat transfer coefficients are measured in a test machine on four cooling air mass flow rates at a wide Couette Reynolds number range. The calculated values concerning the friction and heat transfer coefficients are compared with measured and semi-empirical data. Heat is transferred from the hotter stator and rotor surfaces to the coolerair flow in the air gap, not from the rotor to the stator via the air gap, althought the stator temperature is lower than the rotor temperature. The calculatedfriction coefficients fits well with the semi-empirical equations and precedingmeasurements. On constant mass flow rate the rotor heat transfer coefficient attains a saturation point at a higher rotational speed, while the heat transfer coefficient of the stator grows uniformly. The magnitudes of the heat transfer coefficients are almost constant with different turbulence models. The calibrationof sensors in a straight pipe is only an advisory step in the selection process. Telemetry is tested in the pipe conditions and compared to the same measurements with a plain sensor. The magnitudes of the measured data and the data from the semi-empirical equation are higher for the heat transfer coefficients than thenumerical data considered on the velocity range. Friction and heat transfer coefficients are presented in a large velocity range in the report. The goals are reached acceptably using numerical and experimental research. The next challenge is to achieve results for grooved stator-rotor combinations. The work contains also results for an air gap with a grooved stator with 36 slots. The velocity field by the numerical method does not match in every respect the estimated flow mode. The absence of secondary Taylor vortices is evident when using time averagednumerical simulation.

Prediction of heat transfer on kraft recovery boiler

Työssä tutkittiin Andritz-Ahlstrom toimittamien soodakattiloiden lämmönsiirtoa ANITA 2.20- suunnitteluohjelmalla feedback- laskentaa apuna käyttäen. Data laskentaan saatiin kattiloiden takuukokeissa mitatuista arvoista. Mittaukset on suoritettiin Andritz-Ahlstromin henkilökunnan toimesta tehdashenkilökunnan avustuksella. Feedback -laskenta tapahtui mittaustulosten perusteella, joten tiettyä virhettä luonnollisesti esiintyi. Aluksi laskettiin taseet molempien ekojen yli erikseen sekä molemmat yhdessä Excel-taulukkolaskentaohjelmalla. Täältä saatiin oletettu savukaasuvirtaus kattilassa. Tämän jälkeen lämpöpintoja muokattiin todellisuutta vastaaviksi yleislikaisuuskerrointa muuttamalla (overall fouling factor). Kertoimet ovat liikkuivat noin 0.4 ja 1.6 välillä riipuen kattilan tyypistä ja ANITAn oletuksesta lämpöpintojen likaisuudelle. Havaittin että yhtä varsinaista syytä lämpöpintojen eroavaisuuteen oletetusta ei saatu. Syitä toiminnan poikkeamiseen oli monia. Mm. etukammion koolla havaittiin olevan suurtakin vaikutusta tulistimien, etenkin savukaasuvirrassa ensimmäisen tulistimen toimintaan. Yleisesti todettiin muiden tulistimien vastaavan oletettua toimintaa. Keittopinnan ja ekonomiserien toimintaa tutkittiin hivenen suppeammin ja havaittiin niiden toimivan huomattavasti stabiilimmin kuin tulistimien. Likaisuus kertoimet oletetusta vaihtelivat noin ±20 %.

Hydrodynamics and Heat Transfer in an Airlift Reactor

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

Coolant Flow and Heat Transfer in Pebble-Bed Reactor Core

This thesis gathers knowledge about ongoing high-temperature reactor projects around the world. Methods for calculating coolant flow and heat transfer inside a pebble-bed reactor core are also developed. The thesis begins with the introduction of high-temperature reactors including the current state of the technology. Process heat applications that could use the heat from a high-temperature reactor are also introduced. A suitable reactor design with data available in literature is selected for the calculation part of the thesis. Commercial computational fluid dynamics software Fluent is used for the calculations. The pebble-bed is approximated as a packed-bed, which causes sink terms to the momentum equations of the gas flowing through it. A position dependent value is used for the packing fraction. Two different models are used to calculate heat transfer. First a local thermal equilibrium is assumed between the gas and solid phases and a single energy equation is used. In the second approach, separate energy equations are used for the phases. Information about steady state flow behavior, pressure loss, and temperature distribution in the core is obtained as results of the calculations. The effect of inlet mass flow rate to pressure loss is also investigated. Data found in literature and the results correspond each other quite well, considered the amount of simplifications in the calculations. The models developed in this thesis can be used to solve coolant flow and heat transfer in a pebble-bed reactor, although additional development and model validation is needed for better accuracy and reliability.

Statistical Optimum Design of Heat Exchangers

The optimal design of a heat exchanger system is based on given model parameters together with given standard ranges for machine design variables. The goals set for minimizing the Life Cycle Cost (LCC) function which represents the price of the saved energy, for maximizing the momentary heat recovery output with given constraints satisfied and taking into account the uncertainty in the models were successfully done. Nondominated Sorting Genetic Algorithm II (NSGA-II) for the design optimization of a system is presented and implemented inMatlab environment. Markov ChainMonte Carlo (MCMC) methods are also used to take into account the uncertainty in themodels. Results show that the price of saved energy can be optimized. A wet heat exchanger is found to be more efficient and beneficial than a dry heat exchanger even though its construction is expensive (160 EUR/m2) compared to the construction of a dry heat exchanger (50 EUR/m2). It has been found that the longer lifetime weights higher CAPEX and lower OPEX and vice versa, and the effect of the uncertainty in the models has been identified in a simplified case of minimizing the area of a dry heat exchanger.

Heat Transfer in Plate Heat Exchangers

This study illustrates the different types of plate heat exchangers that are commonly used in various domestic and industrial applications. The main purpose of this paper was to devise a methodology that is capable of calculating optimum number of plates in the design of a plate heat exchanger. To obtain the appropriate number of plates, typically several iterations must be made before a final acceptable design is completed, since plate amount depends on many factors such as, flow velocities, physical properties of the streams, flow channel geometry, allowable pressure drop, plate dimensions, and the gap between the plates. The methodology presented here can be used as a general guide for designing a plate heat exchanger. To investigate the effects of relevant parameters on the thermal-hydraulic design of a plate heat exchanger, several experiments were carried out for single-phase and counter flow arrangement with two brazed plate heat exchangers by varying the flow rates and the inlet temperatures of the fluid streams. The actual heat transfer coefficients obtained based on the experiment were nearly close to the calculated values and to improve the design, a correction factor was introduced. Besides, the effect of flow channel velocity on the pressure drop inside the unit is presented.

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.

Effect of forest-based biofuels production on carbon footprint, Case: Integrated LWC paper mill

International energy and climate strategies also set Finland’s commitments to increasing the use of renewable energy sources and reducing greenhouse gas emissions. The target can be achieved by, for example, increasing the use of energy wood. Finland’s forest biomass potential is significant compared with current use. Increased use will change forest management and wood harvesting methods however. The thesis examined the potential for integrated pulp and paper mills to increase bioenergy production. The effects of two bioenergy production technologies on the carbon footprint of an integrated LWC mill were studied at mill level and from the cradle-to-customer approach. The LignoBoost process and FT diesel production were chosen as bioenergy cases. The data for the LignoBoost process were obtained from Metso and for the FT diesel process from Neste Oil. The rest of the information is based on the literature and databases of the KCL-ECO life-cycle computer program and Ecoinvent. In both case studies, the carbon footprint was reduced. From the results, it can be concluded that it is possible to achieve a fossil-fuel-free pulp mill with the LignoBoost process. By using steam from the FT diesel process, the amount of auxiliary fuel can be reduced considerably and the bark boiler can be replaced. With a choice of auxiliary fuels for use in heat production in the paper mill and the production methods for purchased electricity, it is possible to affect the carbon footprints even more in both cases.

Effect of hydrodynamics on modelling, monitoring and control of crystallization

Crystallization is a purification method used to obtain crystalline product of a certain crystal size. It is one of the oldest industrial unit processes and commonly used in modern industry due to its good purification capability from rather impure solutions with reasonably low energy consumption. However, the process is extremely challenging to model and control because it involves inhomogeneous mixing and many simultaneous phenomena such as nucleation, crystal growth and agglomeration. All these phenomena are dependent on supersaturation, i.e. the difference between actual liquid phase concentration and solubility. Homogeneous mass and heat transfer in the crystallizer would greatly simplify modelling and control of crystallization processes, such conditions are, however, not the reality, especially in industrial scale processes. Consequently, the hydrodynamics of crystallizers, i.e. the combination of mixing, feed and product removal flows, and recycling of the suspension, needs to be thoroughly investigated. Understanding of hydrodynamics is important in crystallization, especially inlargerscale equipment where uniform flow conditions are difficult to attain. It is also important to understand different size scales of mixing; micro-, meso- and macromixing. Fast processes, like nucleation and chemical reactions, are typically highly dependent on micro- and mesomixing but macromixing, which equalizes the concentrations of all the species within the entire crystallizer, cannot be disregarded. This study investigates the influence of hydrodynamics on crystallization processes. Modelling of crystallizers with the mixed suspension mixed product removal (MSMPR) theory (ideal mixing), computational fluid dynamics (CFD), and a compartmental multiblock model is compared. The importance of proper verification of CFD and multiblock models is demonstrated. In addition, the influence of different hydrodynamic conditions on reactive crystallization process control is studied. Finally, the effect of extreme local supersaturation is studied using power ultrasound to initiate nucleation. The present work shows that mixing and chemical feeding conditions clearly affect induction time and cluster formation, nucleation, growth kinetics, and agglomeration. Consequently, the properties of crystalline end products, e.g. crystal size and crystal habit, can be influenced by management of mixing and feeding conditions. Impurities may have varying impacts on crystallization processes. As an example, manganese ions were shown to replace magnesium ions in the crystal lattice of magnesium sulphate heptahydrate, increasing the crystal growth rate significantly, whereas sodium ions showed no interaction at all. Modelling of continuous crystallization based on MSMPR theory showed that the model is feasible in a small laboratoryscale crystallizer, whereas in larger pilot- and industrial-scale crystallizers hydrodynamic effects should be taken into account. For that reason, CFD and multiblock modelling are shown to be effective tools for modelling crystallization with inhomogeneous mixing. The present work shows also that selection of the measurement point, or points in the case of multiprobe systems, is crucial when process analytical technology (PAT) is used to control larger scale crystallization. The thesis concludes by describing how control of local supersaturation by highly localized ultrasound was successfully applied to induce nucleation and to control polymorphism in reactive crystallization of L-glutamic acid.

Limestone Reactions in a Fluidized Bed Heat Exchanger of an Oxy-fuel Circulating Fluidized Bed Boiler

Oxy-fuel combustion in a circulating fluidized bed (CFB) boiler appears to be a promising option for capturing CO2 in power plants. Oxy-fuel combustion is based on burning of fuel in the mixture of oxygen and re-circulated flue gas instead of air. Limestone (CaCO3) is typically used for capturing of SO2 in CFB boilers where limestone calcines to calcium oxide (CaO). Because of high CO2 concentration in oxy-fuel combustion, calcination reaction may be hindered or carbonation, the reverse reaction of calcination, may occur. Carbonation of CaO particles can cause problems especially in the circulation loop of a CFB boiler where temperature level is lower than in the furnace. The aim of the thesis was to examine carbonation of CaO in a fluidized bed heat exchanger of a CFB boiler featuring oxy-fuel combustion. The calculations and analyzing were based on measurement data from an oxy-fuel pilot plant and on 0-dimensional (0D) gas balance of a fluidized bed heat exchanger. Additionally, the objective was to develop a 1-dimensional (1D) model of a fluidized bed heat exchanger by searching a suitable pre-exponential factor for a carbonation rate constant. On the basis of gas measurement data and the 0D gas balance, it was found that the amount of fluidization gas decreased as it flew through the fluidized bed heat exchanger. Most likely the reason for this was carbonation of CaO. It was discovered that temperature has a promoting effect on the reaction rate of carbonation. With the 1D model, a suitable pre-exponential factor for the equation of carbonation rate constant was found. However, during measurements there were several uncertainties, and in the calculations plenty of assumptions were made. Besides, the temperature level in the fluidized bed heat exchanger was relatively low during the measurements. Carbonation should be considered when fluidized bed heat exchangers and the capacity of related fans are designed for a CFB boiler with oxy-fuel combustion.

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.

Communication Solutions for LVDC Island Networks

Production and generation of electrical power is evolving to more environmental friendly technologies and schemes. Pushed by the increasing cost of fossil fuels, the operational costs of producing electrical power with fossil fuels and the effect in the environment, like pollution and global warming, renewable energy sources gain con-stant impulse into the global energy economy. In consequence, the introduction of distributed energy sources has brought a new complexity to the electrical networks. In the new concept of smart grids and decen-tralized power generation; control, protection and measurement are also distributed and requiring, among other things, a new scheme of communication to operate with each other in balance and improve performance. In this research, an analysis of different communication technologies (power line communication, Ethernet over unshielded twisted pair (UTP), optic fiber, Wi-Fi, Wi-MAX, and Long Term Evolution) and their respective characteristics will be carried out. With the objective of pointing out strengths and weaknesses from different points of view (technical, economical, deployment, etc.) to establish a richer context on which a decision for communication approach can be done depending on the specific application scenario of a new smart grid deployment. As a result, a description of possible optimal deployment solutions for communication will be shown considering different options for technologies, and a mention of different important considerations to be taken into account will be made for some of the possible network implementation scenarios.