Kipsin kiteytymiseen vaikuttavat tekijät ammoniumsulfaattikiteyttämön lämmönvaihtimissa

Työssä pyrittiin selvittämään syitä kipsin saostumiseen ammoniumsulfaattikiteyttämön putkilämmönvaihtimien pinnalle ja miten epätoivottua saostumista voitaisiin estää. Lämmönvaihtimissa virtaa ammoniumsulfaattia, jossa on epäpuhtautena kalsiumia, joka saostuu pinnoille kalsiumsulfaattina. Kirjallisuusosassa tarkasteltiin kiteytymisen mekanismia ja kipsin kiteytymiseen vaikuttavia tekijöitä. Saostumien estoaineita ja niiden vaikutusta kipsin kiteytymiseen sekä kipsin liukoisuutta ammoniumsulfaattiliuoksessa käsiteltiin myös kirjallisuusosassa. Kipsin kiteytymiseen vaikuttavia tekijöitä selvitettiin laboratoriokokeilla, joissa pyrittiin simuloimaan lämmönvaihdinta lämpövastuksella. Laboratoriokokeissa kokeiltiin erilaisia saostuman estoaineita ja pyrittiin löytämään prosessiin mahdollisimman tehokas kipsin kiteytymisen estoaine. Lämmönvaihtimien toiminnan tehokkuutta eli muodostuneen saostuman vaikutusta lämmönsiirtymiseen tutkittiin veden luovuttaman lämpövirran avulla. Lämmönvaihtimien tukkeutumista selvitettiin putkien vaihdon tarpeen perusteella. Kalsiumpitoisuuden vaihteluja prosessivirroissa selvitettiin kalsiumtaseen avulla. Saostumiseen vaikuttavien tekijöiden lisäksi selvitettiin mistä ja kuinka paljon kalsiumia kulkeutuu prosessiin ja poistuu sieltä. Työn tarkoituksena oli löytää ratkaisu, jolla epätoivottua saostumista lämmönvaihdin-ten pinnoille pystyttäisiin vähentämään joko kemiallisesti tai muuttamalla prosessi-muuttujia. Kalsiumia havaittiin olevan eniten pelkistämön sisäisissä ammoniumsul-faattiliuoskierroissa. Kalsiumtaseen perusteella kalsiumia poistuu pelkistämöltä eniten kipsinä ammoniumsulfaattituotteen mukana. Laboratoriokokeissa havaittiin polykar-boksylaattien estävän kipsin kasvua parhaiten, joskin estoaineen oikealla annostuksel-la havaittiin olevan suuri vaikutus. Lämmönvaihtimien saostuman havaittiin olevan kipsin ja glauberiitin seos. Vaippapuolen luovuttamien lämpövirran arvojen perusteel-la pystyttiin seuraamaan putkien tukkeutumista.

Virtausjakauman kokeellinen määritys levylämmönsiirtimissä

Diplomityön tavoitteena oli kokeellisen tutkimuksen keinoin selvittää juotettujen levylämmönsiirtimien levypakkarakenteessa virtausten käyttäytyminen ja jakautuminen sekä löytää ideoita ja kehitysehdotuksia levylämmönsiirtimen levypakan ja levyprofiilin kehittämiseksi. Kokeellinen tutkimus suoritettiin Oy Danfoss Ab LPM:n levylämmönsiirtimien tutkimuslaboratoriossa. Virtausjakauman tutkimusta varten suunniteltiin ja valittiin tutkimuslaitteisto, joka koostui termoelementtiantureista, tiedonkeruulaitteistosta sekä ohjelmistosta. Lämmönsiirtimistä mitattiin ensiö- ja toisiopuolen tilavuusvirrat ja painehäviöt sekä lämpötilat ennen ja jälkeen lämmönsiirtimen. Tutkimuslaitteiston avulla mitattiin lämpötiloja lämmönsiirtimen sisältä levyväleistä. Mittaukset suoritettiin neljällä levypakkarakenteella useilla massavirran arvoilla. Mittaustuloksista määritettiin levylämmönsiirtimien lämpö- ja virtaustekniset ominaisuudet nesteen Reynoldsin luvun funktiona sekä selvitettiin nesteen virtausjakaumat. Mittaustuloksien perusteella laskettuja virtausjakauman arvoja verrattiin teorian mukaan laskettuihin jakaumiin. Mitatuista siirtimistä lasketut massavirrat viittaavat siihen, että suurin osa nesteestä virtaa siirtimien keskeltä tai lähempää loppupäätä kuin alkupäästä. Teorian mukaan suurin nestemäärä virtaisi siirtimen alkupäästä vähentyen tasaisesti kohti levypakan loppupäätä. Teorian mukaiselle virtausjakaumalle ja lasketuille jakaumille ei löydetty yhteyttä. Tutkimuksessa havaittiin suuria, jopa yli 20 asteen, lämpötilaeroja levyväleistä ulostulevissa virtauksissa. Levyvälien virtauksen käyttäytymisen ja jakautumisen tutkiminen nähdäänkin levypakan pitkittäistä kehittämistä suurempana mielenkiinnon ja kehittämisen kohteena.

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.

Numerical Investigation of Flow Past a Circular Cylinder and in a Staggered Tube Bundle Using Various Turbulence Models

Transitional flow past a three-dimensional circular cylinder is a widely studied phenomenon since this problem is of interest with respect to many technical applications. In the present work, the numerical simulation of flow past a circular cylinder, performed by using a commercial CFD code (ANSYS Fluent 12.1) with large eddy simulation (LES) and RANS (κ - ε and Shear-Stress Transport (SST) κ - ω! model) approaches. The turbulent flow for ReD = 1000 & 3900 is simulated to investigate the force coefficient, Strouhal number, flow separation angle, pressure distribution on cylinder and the complex three dimensional vortex shedding of the cylinder wake region. The numerical results extracted from these simulations have good agreement with the experimental data (Zdravkovich, 1997). Moreover, grid refinement and time-step influence have been examined. Numerical calculations of turbulent cross-flow in a staggered tube bundle continues to attract interest due to its importance in the engineering application as well as the fact that this complex flow represents a challenging problem for CFD. In the present work a time dependent simulation using κ – ε, κ - ω! and SST models are performed in two dimensional for a subcritical flow through a staggered tube bundle. The predicted turbulence statistics (mean and r.m.s velocities) have good agreement with the experimental data (S. Balabani, 1996). Turbulent quantities such as turbulent kinetic energy and dissipation rate are predicted using RANS models and compared with each other. The sensitivity of grid and time-step size have been analyzed. Model constants sensitivity study have been carried out by adopting κ – ε model. It has been observed that model constants are very sensitive to turbulence statistics and turbulent quantities.

Laitosmittaukset turbiinilaitoksen suorituskyvyn seurantaan ydinvoimalaitoksen vuosihuoltojen yhteydessä

Työssä kehitettiin höyryvoimalaitosprosessin tilan seurantaan laitosmittausrutiini, jolla tarkastellaan erityisesti turbiinilaitoksen suorituskykyä ennen ja jälkeen vuosihuoltojen. Laitosmittaukset on kehitetty Teollisuuden Voima Oyj:n Olkiluodon ydinvoimalaitosyksikköjen OL1 ja OL2 tarpeisiin. Mittauksilla saadaan tietoa vuosihuollon kokonaisvaikutuksista prosessiin eli muuttuiko mikään vuosihuollon aikana ja mitkä syyt johtivat muutokseen. Lisäksi säännöllisillä mittauksilla tavoitellaan pitkän aikavälin referenssitietoa prosessin tilasta. Esimerkkimittaukset tehtiin kevään 2008 vuosihuolloissa, R108 ja R208. Työssä on esitetty mittausten suorittaminen, tulosten laskenta ja tarkastelu sekä raportointi. Luotu mittausrutiini pohjautuu höyryturbiinien ja lämmönsiirtimien vastaanottokokeita käsitteleviin standardeihin sekä laitosten viimeisimpiin prosessimuutosten yhteydessä toteutettuihin suorituskykymittauksiin. Laitosmittauksista tehtiin vastaanottokokeita kevyemmät ja yksinkertaisemmat mittaustapahtumat, joilla saadaan kuitenkin riittävän luotettavaa tietoa. Tunnuslukuja, joita ei käytettävissä olevien mittatietojen avulla saatu suoraan selville tai laskettua massa- ja energiataseilla, määritettiin stationaarisen tilan mallinnusohjelmalla TEMPO:lla. Esimerkkimittaukset sujuivat hyvin ja saadut tulokset antoivat kaivattua lisätietoa vuosihuollon vaikutuksesta prosessiin. Generaattorin bruttoteho aleni OL1:llä ja pysyi ennallaan OL2:lla. OL1:n tehon aleneminen selittyi syöttöveden virtausmäärän vähenemisellä, johon vaikutti virtausmittalaitteen lähettimen vaihto. OL2:lla ei havaittu turbiiniprosessissa muutoksia, mutta lauhduttimien suorituskyky parani. Lauhduttimien suorituskyvyn paraneminen ilmeni asteisuuksien pienentymisenä, johon vaikuttivat puhdistaminen ja ilmavuotojen väheneminen.

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.

Variable-speed-drive-based detection methods for problems in fluid handling systems

Fluid handling systems account for a significant share of the global consumption of electrical energy. They also suffer from problems, which reduce their energy efficiency and increase life-cycle costs. Detecting or predicting these problems in time can make fluid handling systems more environmentally and economically sustainable to operate. In this Master’s Thesis, significant problems in fluid systems were studied and possibilities to develop variable-speed-drive-based detection methods for them was discussed. A literature review was conducted to find significant problems occurring in fluid handling systems containing pumps, fans and compressors. To find case examples for evaluating the feasibility of variable-speed-drive-based methods, queries were sent to industrial companies. As a result of this, the possibility to detect heat exchanger fouling with a variable-speed drive was analysed with data from three industrial cases. It was found that a mass flow rate estimate, which can be generated with a variable speed drive, can be used together with temperature measurements to monitor a heat exchanger’s thermal performance. Secondly, it was found that the fouling-related increase in the pressure drop of a heat exchanger can be monitored with a variable speed drive. Lastly, for systems where the flow device is speed controlled with by a pressure measurement, it was concluded that increasing rotational speed can be interpreted as progressing fouling in the heat exchanger.

Interactions between the physical soil environment and a horizontal ground coupled heat pump, for a domestic site in the UK

There is currently an increased interest of Government and Industry in the UK, as well as at the European Community level and International Agencies (i.e. Department of Energy, American International Energy Agency), to improve the performance and uptake of Ground Coupled Heat Pumps (GCHP), in order to meet the 2020 renewable energy target. A sound knowledge base is required to help inform the Government Agencies and advisory bodies; detailed site studies providing reliable data for model verification have an important role to play in this. In this study we summarise the effect of heat extraction by a horizontal ground heat exchanger (installed at 1 m depth) on the soil physical environment (between 0 and 1 m depth) for a site in the south of the UK. Our results show that the slinky influences the surrounding soil by significantly decreasing soil temperatures. Furthermore, soil moisture contents were lower for the GCHP soil profile, most likely due to temperature-gradient related soil moisture migration effects and a decreased hydraulic conductivity, the latter as a result of increased viscosity (caused by the lower temperatures for the GCHP soil profile). The effects also caused considerable differences in soil thermal properties. This is the first detailed mechanistic study conducted in the UK with the aim to understand the interactions between the soil, horizontal heat exchangers and the aboveground environment. An increased understanding of these interactions will help to achieve an optimum and sustainable use of the soil heat resources in the future. The results of this study will help to calibrate and verify a simulation model that will provide UK-wide recommendations to improve future GCHP uptake and performance, while safeguarding the soil physical resources.

Simulation of a domestic ground source heat pump system using a three-dimensional numerical borehole heat exchanger model

Common approaches to the simulation of borehole heat exchangers (BHEs) assume heat transfer in circulating fluid and grout to be in a quasi-steady state and ignore fluctuations in fluid temperature due to transport of the fluid around the loop. However, in domestic ground source heat pump (GSHP) systems, the heat pump and circulating pumps switch on and off during a given hour; therefore, the effect of the thermal mass of the circulating fluid and the dynamics of fluid transport through the loop has important implications for system design. This may also be important in commercial systems that are used intermittently. This article presents transient simulation of a domestic GSHP system with a single BHE using a dynamic three-dimensional (3D) numerical BHE model. The results show that delayed response associated with the transit of fluid along the pipe loop is of some significance in moderating swings in temperature during heat pump operation. In addition, when 3D effects are considered, a lower heat transfer rate is predicted during steady operations. These effects could be important when considering heat exchanger design and system control. The results will be used to develop refined two-dimensional models.

Photovoltaic and solar-assisted ground-source heat pump systems

This review investigates the performance of photovoltaic and solar-assisted ground-source heat pumps in which solar heat is transferred to the ground to improve the coefficient of performance. A number of studies indicate that, for systems with adequately sized ground heat exchangers, the effect on system efficiency is small: about 1% improvement if the heat source is photovoltaic, a 1–2% decline if the source is solar thermal. With possible exceptions for systems in which the ground heat exchanger is undersized, or natural recharge from ground water is insufficient, solar thermal energy is better used for domestic hot water than to recharge ground heat. This appears particularly true outside the heating season, as although much of the heat extracted from the ground can be replaced, it seems to have little effect on the coefficient of performance. Any savings in electrical consumption that do result from an improved coefficient can easily be outweighed by an inefficient control system for the circulation pumps.

Identifiering av parametrar till ackumulatortanken Solus1050 från consolar - en metodbeskrivning

This report describes the work done creating a computer model of a kombi tank from Consolar. The model was created with Presim/Trnsys and Fittrn and DF were used to identify the parameters. Measurements were carried out and were used to identify the values of the parameters in the model. The identifications were first done for every circuit separately. After that, all parameters are normally identified together using all the measurements. Finally the model should be compared with other measurements, preferable realistic ones. The two last steps have not yet been carried out, because of problems finding a good model for the domestic hot water circuit.The model of the domestic hot water circuit give relatively good results for low flows at 5 l/min, but is not good for higher flows. In the report suggestions for improving the model are given. However, there was not enough time to test this within the project as much time was spent trying to solve problems with the model crashing. Suggestions for improving the model for the domestic circuit are given in chapter 4.4. The improved equations that are to be used in the improved model are given by equation 4.18, 4.19 and 4.22.Also for the boiler circuit and the solar circuit there are improvements that can be done. The model presented here has a few shortcomings, but with some extra work, an improved model can be created. In the attachment (Bilaga 1) is a description of the used model and all the identified parameters.A qualitative assessment of the store was also performed based on the measurements and the modelling carried out. The following summary of this can be given: Hot Water PreparationThe principle for controlling the flow on the primary side seems to work well in order to achieve good stratification. Temperatures in the bottom of the store after a short use of hot water, at a coldwater temperature of 12°C, was around 28-30°C. This was almost independent of the temperature in the store and the DHW-flow.The measured UA-values of the heat exchangers are not very reliable, but indicates that the heat transfer rates are much better than for the Conus 500, and in the same range as for other stores tested at SERC.The function of the mixing valve is not perfect (see diagram 4.3, where Tout1 is the outlet hot water temperature, and Tdhwo and Tdhw1 is the inlet temperature to the hot and cold side of the valve respectively). The outlet temperature varies a lot with different temperatures in the storage and is going down from 61°C to 47°C before the cold port is fully closed. This gives a problem to find a suitable temperature setting and gives also a risk that the auxiliary heating is increased instead of the set temperature of the valve, when the hot water temperature is to low.Collector circuitThe UA-value of the collector heat exchanger is much higher than the value for Conus 500, and in the same range as the heat exchangers in other stores tested at SERC.Boiler circuitThe valve in the boiler circuit is used to supply water from the boiler at two different heights, depending on the temperature of the water. At temperatures from the boiler above 58.2°C, all the water is injected to the upper inlet. At temperatures below 53.9°C all the water is injected to the lower inlet. At 56°C the water flow is equally divided between the two inlets. Detailed studies of the behaviour at the upper inlet shows that better accuracy of the model would have been achieved using three double ports in the model instead of two. The shape of the upper inlet makes turbulence, that could be modelled using two different inlets. Heat lossesThe heat losses per m3 are much smaller for the Solus 1050, than for the Conus 500 Storage. However, they are higher than those for some good stores tested at SERC. The pipes that are penetrating the insulation give air leakage and cold bridges, which could be a major part of the losses from the storage. The identified losses from the bottom of the storage are exceptionally high, but have less importance for the heat losses, due to the lower temperatures in the bottom. High losses from the bottom can be caused by air leakage through the insulation at the pipe connections of the storage.

Rapport från/Report from IEA Task-14 mötet/meeting on solar heating system, Hameln Tyskland/Germany August 1992

Participation as observer at the meeting of Task 14 of IEA's Solar Heating and Cooling Projects held in Hameln, Germany has led to greater understanding of interesting developments underway in several countries. This will be of use during the development of small scale systems suitable for Swedish conditions. A summary of the work carried out by the working groups within Task 14 is given, with emphasis on the Domestic Hot Water group. Experiences of low-flow systems from several countries are related, and the conclusion is drawn that the maximum theoretical possible increase in performance of 20% has not been achieved due to poor heat exchangers and poor stratification in the storage tanks. Positive developments in connecting tubes and pumps is noted. Further participation as observer in Task 14 meetings is desired, and is looked on favourably by the members of the group. Another conclusion is that SERC should carry on with work on Swedish storage tanks, with emphasis on better stratification and heat exchangers, and possible modelling of system components. Finally a German Do-it-Vourself kit is described and judged in comparison with prefabricated models and Swedish Do-it-Yourself kits.

Influence of hydraulics and control of thermal storage in solar assisted heat pump combisystems

This paper studies the influence of hydraulics and control of thermal storage in systems combined with solar thermal and heat pump for the production of warm water and space heating in dwellings. A reference air source heat pump system with flat plate collectors connected to a combistore was defined and modeled together with the IEA SHC Task 44 / HPP Annex 38 (T44A38) “Solar and Heat Pump Systems” boundary conditions of Strasbourg climate and SFH45 building. Three and four pipe connections as well as use of internal and external heat exchangers for DHW preparation were investigated as well as sensor height for charging of the DHW zone in the store. The temperature in this zone was varied to ensure the same DHW comfort was achieved in all cases. The results show that the four pipe connection results in 9% improvement in SPF compared to three pipe and that the external heat exchanger for DHW preparation leads to a 2% improvement compared to the reference case. Additionally the sensor height for charging the DHW zone of the store should not be too low, otherwise system performance is adversely affected

Non-newtonian heat transfer on a plate heat exchanger with generalized configurations

For the configuration optimization of plate heat exchangers (PHEs), the mathematical models for heat transfer and pressure drop must be valid for a wide range of operational conditions of all configurations of the exchanger or the design results may be compromised. In this investigation, the thermal model of a PHE is adjusted to fit experimental data obtained from non-Newtonian heat transfer for eight different configurations, using carboxymethylcellulose solutions (CMC) as test fluid. Although it is possible to successfully adjust the model parameters, Newtonian and non-Newtonian heat transfer cannot be represented by a single generalized correlation. In addition, the specific heat, thermal conductivity and power-law rheological parameters of CMC solutions were correlated with temperature, over a range compatible with a continuous pasteurization process.

Transient non-Darcy forced convective heat transfer from a flat plate embedded in a fluid-saturated porous medium

Transient non-Darcy forced convection on a flat plate embedded in a porous medium is investigated using the Forchheimer-extended Darcy law. A sudden uniform pressure gradient is applied along the flat plate, and at the same time, its wall temperature is suddenly raised to a high temperature. Both the momentum and energy equations are solved by retaining the unsteady terms. An exact velocity solution is obtained and substituted into the energy equation, which then is solved by means of a quasi-similarity transformation. The temperature field can be divided into the one-dimensional transient (downstream) region and the quasi-steady-state (upstream) region. Thus the transient local heat transfer coefficient can be described by connecting the quasi-steady-state solution and the one-dimensional transient solution. The non-Darcy porous inertia works to decrease the velocity level and the time required for reaching the steady-state velocity level. The porous-medium inertia delays covering of the plate by the steady-state thermal boundary layer. © 1990.