Gas flow analysis in a Kraft recovery boiler

Demands for optimal boiler performance and increased concerns in lowering emission have always been the driving force in the reevaluation and evolution of the Kraft boiler: specifically the air distribution strategies that are directly related to achieving increased residence time of flue gas combustion inside the furnace which in turn lowers atmosphere emission levels and enhances boiler operation. This paper presents the results of a study that analyzes the interaction of the different multilevel air injections have on flue gas flow patterns including various quaternary air supply arrangements. Additionally, this study assesses the performance of the CFD (Computational Fluid Dynamics) model against data available in literature. Simulations were performed considering isothermal and incompressible flows, and did not take into account thermal phenomena or chemical reactions. The numerical solutions generated proved to be coherently related to the data available in literature, and provided proof of the efficiency of tertiary level air injection, as well as revealed that quaternary air injection ports arranged in a symmetrical configuration is most suitable for optimal equipment operation. (C) 2010 Elsevier B.V. All rights reserved.

Combustion of alternative biomass fuels in a domestic boiler and in a large-scale furnace

Programa Doutoral em Engenharia Mecânica.

Report on a review of the Boiler and Pressure Vessel Safety Program and the Elevator and Escalator Safety Program administered by Iowa Workforce Development for the period July 1, 2008 through June 30, 2010

Report on a review of the Boiler and Pressure Vessel Safety Program and the Elevator and Escalator Safety Program administered by Iowa Workforce Development for the period July 1, 2008 through June 30, 2010

Dynamic Model of a Small Once-Through Boiler

In this study, a model for the unsteady dynamic behaviour of a once-through counter flow boiler that uses an organic working fluid is presented. The boiler is a compact waste-heat boiler without a furnace and it has a preheater, a vaporiser and a superheater. The relative lengths of the boiler parts vary with the operating conditions since they are all parts of a single tube. The present research is a part of a study on the unsteady dynamics of an organic Rankine cycle power plant and it will be a part of a dynamic process model. The boiler model is presented using a selected example case that uses toluene as the process fluid and flue gas from natural gas combustion as the heat source. The dynamic behaviour of the boiler means transition from the steady initial state towards another steady state that corresponds to the changed process conditions. The solution method chosen was to find such a pressure of the process fluid that the mass of the process fluid in the boiler equals the mass calculated using the mass flows into and out of the boiler during a time step, using the finite difference method. A special method of fast calculation of the thermal properties has been used, because most of the calculation time is spent in calculating the fluid properties. The boiler was divided into elements. The values of the thermodynamic properties and mass flows were calculated in the nodes that connect the elements. Dynamic behaviour was limited to the process fluid and tube wall, and the heat source was regarded as to be steady. The elements that connect the preheater to thevaporiser and the vaporiser to the superheater were treated in a special way that takes into account a flexible change from one part to the other. The model consists of the calculation of the steady state initial distribution of the variables in the nodes, and the calculation of these nodal values in a dynamic state. The initial state of the boiler was received from a steady process model that isnot a part of the boiler model. The known boundary values that may vary during the dynamic calculation were the inlet temperature and mass flow rates of both the heat source and the process fluid. A brief examination of the oscillation around a steady state, the so-called Ledinegg instability, was done. This examination showed that the pressure drop in the boiler is a third degree polynomial of the mass flow rate, and the stability criterion is a second degree polynomial of the enthalpy change in the preheater. The numerical examination showed that oscillations did not exist in the example case. The dynamic boiler model was analysed for linear and step changes of the entering fluid temperatures and flow rates.The problem for verifying the correctness of the achieved results was that there was no possibility o compare them with measurements. This is why the only way was to determine whether the obtained results were intuitively reasonable and the results changed logically when the boundary conditions were changed. The numerical stability was checked in a test run in which there was no change in input values. The differences compared with the initial values were so small that the effects of numerical oscillations were negligible. The heat source side tests showed that the model gives results that are logical in the directions of the changes, and the order of magnitude of the timescale of changes is also as expected. The results of the tests on the process fluid side showed that the model gives reasonable results both on the temperature changes that cause small alterations in the process state and on mass flow rate changes causing very great alterations. The test runs showed that the dynamic model has no problems in calculating cases in which temperature of the entering heat source suddenly goes below that of the tube wall or the process fluid.

Ash Particle Erosion on Steam Boiler Convective Section

In this study, equations for the calculation of erosion wear caused by ash particles on convective heat exchanger tubes of steam boilers are presented. Anew, three-dimensional test arrangement was used in the testing of the erosion wear of convective heat exchanger tubes of steam boilers. When using the sleeve-method, three different tube materials and three tube constructions could be tested. New results were obtained from the analyses. The main mechanisms of erosionwear phenomena and erosion wear as a function of collision conditions and material properties have been studied. Properties of fossil fuels have also been presented. When burning solid fuels, such as pulverized coal and peat in steam boilers, most of the ash is entrained by the flue gas in the furnace. In bubbling andcirculating fluidized bed boilers, particle concentration in the flue gas is high because of bed material entrained in the flue gas. Hard particles, such as sharp edged quartz crystals, cause erosion wear when colliding on convective heat exchanger tubes and on the rear wall of the steam boiler. The most important ways to reduce erosion wear in steam boilers is to keep the velocity of the flue gas moderate and prevent channelling of the ash flow in a certain part of the cross section of the flue gas channel, especially near the back wall. One can do this by constructing the boiler with the following components. Screen plates can beused to make the velocity and ash flow distributions more even at the cross-section of the channel. Shield plates and plate type constructions in superheaters can also be used. Erosion testing was conducted with three types of tube constructions: a one tube row, an inline tube bank with six tube rows, and a staggered tube bank with six tube rows. Three flow velocities and two particle concentrations were used in the tests, which were carried out at room temperature. Three particle materials were used: quartz, coal ash and peat ash particles. Mass loss, diameter loss and wall thickness loss measurements of the test sleeves were taken. Erosion wear as a function of flow conditions, tube material and tube construction was analyzed by single-variable linear regression analysis. In developing the erosion wear calculation equations, multi-variable linear regression analysis was used. In the staggered tube bank, erosion wear had a maximum value in a tube row 2 and a local maximum in row 5. In rows 3, 4 and 6, the erosion rate was low. On the other hand, in the in-line tube bank the minimum erosion rate occurred in tube row 2 and in further rows the erosion had an increasing value, so that in a six row tube bank, the maximum value occurred in row 6.

Environmental requirements and environmental management system of boiler plant supplier

Työn tärkeimpänä tavoitteena oli edistää ympäristöjärjestelmän laatimista Kvaerner Pulping, Power divisioonan Kattilat-liiketoimintayksikölle. Lisäksi tavoitteena oli tarkastella yritykseen kohdistuvia ympäristövaatimuksia ja niiden vaikutusta yrityksen toimintaan. Aluksi työssä tarkasteltiin ympäristöjärjestelmästandardien sisältöjä ja niiden eroja. Työssä käsiteltiin myös erilaisia elinkaarijohtamisen malleja, joita voidaan hyödyntää yrityksen kokonaisvaltaisessa ympäristöjärjestelmässä. Työssä tarkasteltiin myös sidosryhmien vaikutusta yrityksen ympäristötoimintaan. Kattilalaitostoimittajan tärkeimpiä asiakkaita ovat sellu- ja paperiteollisuus. Näihin yrityksiin on kohdistunut runsaan kymmenen vuoden aikana paineita ympäristötoiminnan tehostamiseksi. Tämän kehityksen seurauksena vastaavat tehostamispaineet ovat siirtymässä myös alihankkijoille, kuten kattilaitostoimittajille. Tehokkaan ympäristöjohtamisen takaamiseksi työssä määriteltiin ympäristövastuut ja –valtuudet sekä ympäristöpäämäärät ja –tavoitteet. Lisäksi tunnistettiin yrityksen toimintaan liittyvät ympäristönäkökohdat ja –riskit. Työn yhteydessä laadittiin ympäristöjärjestelmän luonnos, ja se sisältää jatkuvan parantamisen periaatteen. Työn yhteydessä laadittiin lisäksi Kvaerner Pulping Oy:n koelaitokselle ympäristölupahakemus. Työssä on kuvattu koelaitosta ja sen ympäristölupahakemukseen liittyviä asioita esimerkkinä parantuneesta ympäristöasioiden hoidosta.

Development of a PC Tool for Estimating Pressure Drops in Forced-convection Boiler Tubes

Kaksifaasivirtauksen kuvaamiseen käytettävät mallit, ja menetelmät kaksifaasivirtauksen painehäviön määrittämiseksi kehittyvät yhä monimutkaisimmiksi. Höyrystinputkissa tapahtuvien painehäviöiden arvioinnin vaatiman laskennan suorittamiseksi tietokoneohjelman kehittäminen on välttämätöntä. Tässä työssä on kehitetty itsenäinen PC-ohjelma painehäviöiden arvioimiseksi pakotetulle konvektiovirtaukselle pystysuorissa höyrykattilan höyrystinputkissa. Veden ja vesihöyryn aineominaisuuksien laskentaan käytetään IAPWS-IF97 –yhtälökokoelmaa sekä muita tarvittavia IAPWS:n suosittelemia yhtälöitä. Höyrystinputkessa kulloinkin vallitsevan virtausmuodon määrittämiseen käytetään sovelluskelpoisia virtausmuotojen välisiä rajoja kuvaavia yhtälöitä. Ohjelmassa käytetään painehäviön määritykseen kirjallisuudessa julkaistuja yhtälöitä, virtausmuodosta riippuen, alijäähtyneelle virtaukselle, kupla-, tulppa- ja rengasvirtaukselle sekä tulistetun höyryn virtaukselle. Ohjelman laskemia painehäviöarvioita verrattiin kirjallisuudesta valittuihin mittaustuloksiin. Laskettujen painehäviöiden virhe vaihteli välillä –19.5 ja +23.9 %. Virheiden itseisarvojen keskiarvo oli 12.8 %.

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 %.

Advanced control methods for reducing nitrogen oxides in a fluidized bed boiler

Tämän työn tarkoituksena oli löytää keinoja erään leijukerroskattilan typenoksidipäästöjen vähentämiseksi. Koska päästöt olivat jo alunperin alhaiset leijukerrostekniikan ja hybridin SNCR/SCR –typenpoistolaitteiston ansiosta, päätettiin päästöjä lähteä vähentämään parantamalla ammoniakkiruiskutuksen säätöä. Alkuperäinen ammoniakkiruiskutuksen säätö oli liian hidas, jotta satunnaisten häiriöiden aiheuttamat typenoksidipiikit olisi pystytty poistamaan. Ammoniakkiruiskutusta parannettiin lisäämällä jokaiseen ammoniakkilinjaan mäntäpumput, joiden avulla ammoniakkia voidaan syöttää sinne, missä sitä eniten tarvitaan. Ammoniakkiruiskutuksen säätöön kehitettiin uusi sumeaan logiikkaan perustuva säätäjä. Myös muita kehittyneitä säätömenetelmiä kuten neuroverkkoa hyödynnettiin säätäjän kehityksessä. Ammoniakkiruiskutuksen säätäjää testattiin menestyksekkäästi Ruotsissa Brista Kraftin Märstassa sijaitsevalla voimalaitoksella

Development of Furnace Model for Bubbling Fluidized Bed Boiler

The aim of this thesis was to develop a model, which can predict heat transfer, heat release distribution and vertical temperature profile of gas phase in the furnace of a bubbling fluidized bed (BFB) boiler. The model is based on three separate model components that take care of heat transfer, heat release distribution and mass and energy balance calculations taking into account the boiler design and operating conditions. The model was successfully validated by solving the model parameters on the basis of commercial size BFB boiler test run information and by performing parametric studies with the model. Implementation of the developed model for the Foster Wheeler BFB design procedures will require model validation with existing BFB database and possibly more detailed measurements at the commercial size BFB boilers.

Characterization Method for Comminution Behavior of Limestone and Fuel Ash in a Circulating Fluidized Bed Boiler

Työn tavoite oli kehittää karakterisointimenetelmät kalkkikiven ja polttoaineen tuhkan jauhautumisen ennustamiselle kiertoleijukattilan tulipesässä. Kiintoainekäyttäytymisen karakterisoinnilla ja mallintamisella voidaan tarkentaa tulipesän lämmönsiirron ja tuhkajaon ennustamista. Osittain kokeelliset karakterisointimenetelmät perustuvat kalkkikiven jauhautumiseen laboratoriokokoluokan leijutetussa kvartsiputkireaktorissa ja tuhkan jauhatumiseen rotaatiomyllyssä. Karakterisointimenetelmät ottavat huomioon eri-laiset toimintaolosuhteet kaupallisen kokoluokan kiertoleijukattiloissa. Menetelmät kelpoistettiin kaupallisen kokoluokan kiertoleijukattiloista mitattujen ja fraktioittaisella kiintoainemallilla mallinnettujen taseiden avulla. Kelpoistamistaseiden vähäisyydestä huolimatta karakterisointimenetelmät arvioitiin virhetarkastelujen perusteella järkeviksi. Karakterisointimenetelmien kehittämistä ja tarkentamista tullaan jatkamaan.

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.

Offdesign operation of kraft recovery boiler

A model to solve heat and mass balances during the offdesign load calculations was created. These equations are complex and nonlinear. The main new ideas used in the created offdesign model of a kraft recovery boiler are the use of heat flows as torn iteration variables instead of the current practice of using the mass flows, vectorizing equation solving, thus speeding up the process, using non dimensional variables for solving the multiple heat transfer surface problem and using a new procedure for calculating pressure losses. Recovery boiler heat and mass balances are reduced to vector form. It is shown that these vectorized equations can be solved virtually without iteration. The iteration speed is enhanced by the use of the derived method of calculating multiple heat transfer surfaces simultaneously. To achieve this quick convergence the heat flows were used as the torn iteration parameters. A new method to handle pressure loss calculations with linearization was presented. This method enabled less time to be spent calculating pressure losses. The derived vector representation of the steam generator was used to calculate offdesign operation parameters for a 3000 tds/d example recovery boiler. The model was used to study recovery boiler part load operation and the effect of the black liquor dry solids increase on recovery boiler dimensioning. Heat flows to surface elements for part load calculations can be closely approximated with a previously defined exponent function. The exponential method can be used for the prediction of fouling in kraft recovery boilers. For similar furnaces the firing of 80 % dry solids liquor produces lower hearth heat release rate than the 65 % dry solids liquor if we fire at constant steam flow. The furnace outlet temperatures show that capacity increase with firing rate increase produces higher loadings than capacity increase with dry solids increase. The economizers, boiler banks and furnaces can be dimensioned smaller if we increase the black liquor dry solids content. The main problem with increased black liquor dry solids content is the decrease in the heat available to superheat. Whenever possible the furnace exit temperature should be increased by decreasing the furnace height. The increase in the furnace exit temperature is usually opposed because of fear of increased corrosion.

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.

Modern recovery boiler design criteria

Engineering and pricing of large recovery boiler were studied in this work. Engineering was carried out with Anita 4.2 which is an engineering program of Andritz. Key initial values were chosen with previous studies. Primary target of this work was to find out the consequences that furnace dimensions and furnace screen vertical part has to boiler pricing. Boilers that were engineered had different rate of furnace width and depth and different heat transfer plate count. Boiler balances were invariable. Boilers with different vertical screen construction were also calculated. First variation was boiler with vertical screen up to furnace roof. Other variation was to connect vertical screen to Pre-boiler generating bank inlet tubes. Total prices were calculated to engineered boilers. Pricing was sort out to heat transfers, high pressure pipes, steel structures, auxiliary equipments and civil/structural costs. This study did not notice parts of the boiler which costs do not vary with the construction of the boiler. Heat transfers had the largest share of costs. Boiler building had the most significant differences between the boilers. Furnace screen had also significant role especially to costs of the boiler building.