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.

Metsähakkeen jalostusarvon nostaminen eri kuivausmenetelmillä

The thesis was made for Hyötypaperi Oy. To the business activity of the company belongs the recycling of materials and carry out for re-use and the manufacture of solid biofuels and solid recovered fuels. Hyötypaperi Oy delivers forest chips to its partner incineration plants by day and night though the year. The value of the forest chips is based on its percentage of dry material. It is important to dry forest chips well before its storage in piles and delivering to incineration plants. In the thesis was examed the increasing of the degree of refinement of forest chips by different drying methods. In the thesis was examined four different drying methods of the forest chips. The methods were field-, plate-, platform- and channel drying. In the channel drying was used a mechanical blower and the other drying methods were based on the weather conditions. By all drying methods were made test dryings during the summer 2007. In the thesis was examined also the economical profitableness of the field- and channel drying. The last examination in the thesis was measuring of the forest chips’s humidity by humidity measuring equipment of sawn timber during November 2007. The field drying on the surface of asphalt is the only method of drying, which is used by Hyötypaperi Oy in its own production. There do not exist earlier properly examined facts of any drying methods of forest chips, because the drying of forest chips is a new branch. By field- and platform drying achieved lower humidity of forest chips than by plate drying. The object by using the humidity measuring equipment was to be informed of the humidity of forest chips. At present the humidity will find out after 24 hours when the sample of humidity quantity has been dried in the oven. The Lappeenranta University of Technology had the humidity measuring equipment of sawn timber. The values of humidity measured by the equipment from the sample of forest chips were 2 – 9 percent lower than the real values of humidity specified by drying in oven.

Sähkön tuotantokustannusvertailu

Työssä vertaillaan eri sähköntuotantovaihtoehtojen taloudellista kannattavuutta. Kannattavuusvertailu suoritetaan pelkkää sähköä tuottaville voimalaitoksille. Sähkön ja lämmön yhteistuotannon lisärakentaminen tulee kattamaan tietyn osuuden lähitulevaisuuden sähkön hankinnan vajeesta, mutta sen lisäksi tarvitaan myös uutta lauhdetuotantokapasiteettia. Tutkittavat voimalaitostyypit ovat: ydinvoimalaitos, maakaasukombilauhdevoimalaitos, kivihiililauhdevoimalaitos, turvelauhdevoimalaitos, puulauhdevoimalaitos ja tuulivoimala. Kannattavuustarkastelu suoritetaan annuiteettimenetelmällä käyttäen 5 % reaalikorkoa ja tammikuun 2008 hintatasoa. Laskelmien perusteella 8000 tunnin huipunkäyttöajalla ydinsähkön tuotantokustannus olisi 35,0 € /MWh, kaasusähkön 59,2 €/MWh ja hiilisähkön 64,4 €/MWh, kun hiilidioksidipäästöoikeuden hintana käytetään 23 €/t. Ilman päästökauppaa kaasusähkön hinta on 51,2 €/MWh ja hiilisähkön 45,7 €/MWh ydinsähkön hinnan pysyessä ennallaan. Herkkyystarkastelun tulosten perusteella ydinvoiman kilpailukyky korostuu muihin tarkasteltuihin tuotantomuotoihin verrattuna. Ydinpolttoaineen suurellakaan hinnan muutoksella ei ole merkittävää vaikutusta ydinsähkön tuotantokustannukseen, kun taas maakaasusähkö on erittäin riippuvainen polttoaineen hinnasta. Myös päästöoikeuden hinnan kasvu lisää merkittävästi ydinvoiman kilpailukykyä kaasu- ja hiilisähköön verrattuna. Ydinvoimainvestoinnin kannattavuutta ja takaisinmaksua tarkastellaan myös yksinään siten, että investoinnilla saavutettavien tuottojen laskennassa käytetään useita eri sähkön markkinahintoja. Investoinnin kannattavuus on erittäin hyvä, kun sähkön markkinahinta on 50 €/MWh tai suurempi.

Comparison of electricity generation costs

The economical competitiveness of various power plant alternatives is compared. The comparison comprises merely electricity producing power plants. Combined heat and power (CHP) producing power will cover part of the future power deficit in Finland, but also condensing power plants for base load production will be needed. The following types of power plants are studied: nuclear power plant, combined cycle gas turbine plant, coal-fired condensing power plant, peat-fired condensing power plant, wood-fired condensing power plant and wind power plant. The calculations are carried out by using the annuity method with a real interest rate of 5 % per annum and with a fixed price level as of January 2008. With the annual peak load utilization time of 8000 hours (corresponding to a load factor of 91,3 %) the production costs would be for nuclear electricity 35,0 €/MWh, for gas based electricity 59,2 €/MWh and for coal based electricity 64,4 €/MWh, when using a price of 23 €/tonCO2 for the carbon dioxide emission trading. Without emission trading the production cost of gas electricity is 51,2 €/MWh and that of coal electricity 45,7 €/MWh and nuclear remains the same (35,0 €/MWh) In order to study the impact of changes in the input data, a sensitivity analysis has been carried out. It reveals that the advantage of the nuclear power is quite clear. E.g. the nuclear electricity is rather insensitive to the changes of nuclear fuel price, whereas for natural gas alternative the rising trend of gas price causes the greatest risk. Furthermore, increase of emission trading price improves the competitiveness of the nuclear alternative. The competitiveness and payback of the nuclear power investment is studied also as such by using various electricity market prices for determining the revenues generated by the investment. The profitability of the investment is excellent, if the market price of electricity is 50 €/MWh or more.

Radiative Heat Transfer in Boiler Furnaces

Better models are needed for radiative heat transfer in boiler furnaces. If the process is known better, combustion in the furnace can be optimized to produce low emissions. It makes the process to be environmental friendly. Furthermore, if there is a better model of the furnace it can more fully explain what is happening inside the furnace. Using of the model one can quickly and easily analyze how it operates with bio fuels, moist fuels or difficult fuels and improve the operation. Models helps with better estimation of furnace dimensions and result in more accurate understanding of operation. Key component lacking in these models is radiative heat transfer in particle laden gases. If there are no particles than radiative heat transfer can be calculated approximately. There are two problems with current models when used with flow modeling. The first one is a need to account for a particle laden gas and the second one is an absence of a fast algorithm. Fast calculation is needed if radiative heat transfer calculation is done for a large CDF model. Computations slow down if time is required for calculating radiative properties over and over again. This thesis presents a band model for radiative heat transfer in boiler furnaces. Advantage is a quickness of calculation and account of particles in the process.

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.

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.

Heat transfer in superheaters of the bubbling fluidized bed boiler conversions

It is often reasonable to convert old boiler to bubbling fluidized bed boiler instead of building a new one. Converted boiler consists of old and new heat surfaces which must be fitted to operate together. Prediction of heat transfer in not so ideal conditions sets challenges for designers. Two converted boilers situated in Poland were studied on the grounds of acceptance tests and further studies. Calculation of boiler process was performed with boiler design program. Main interest was heat transfer in superheaters and factors affecting it. Theory for heat transfer is presented according to information found from literature. Results obtained from experimental studies and calculations have been compared. With correct definitions calculated parameters corresponded well to measured data at boiler maximum design load. However overload situations revealed to be difficult to model at least without considering changes in the combustion process which requires readjustments to the design program input values.

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.

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.