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.

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.

Optimizing the Performance of a Heat Exchanger Network

This research is the continuation and a joint work with a master thesis that has been done in this department recently by Hemamali Chathurangani Yashika Jayathunga. The mathematical system of the equations in the designed Heat Exchanger Network synthesis has been extended by adding a number of equipment; such as heat exchangers, mixers and dividers. The solutions of the system is obtained and the optimal setting of the valves (Each divider contains a valve) is calculated by introducing grid-based optimization. Finding the best position of the valves will lead to maximization of the transferred heat in the hot stream and minimization of the pressure drop in the cold stream. The aim of the following thesis will be achieved by practicing the cost optimization to model an optimized network.

Product Standardization of Shell and Tube Heat Exchangers

Diplomityön tavoitteena oli selvittää, kuinka kohdeyrityksessä nykyhetkellä pitkältiräätälöitynä tuotetut putkilämmönvaihtimet voitaisiin standardisoida, jotta pystyttäisiin kohdistamaan resursseja oikein ja parantamaan siten yrityksen taloudellista suorituskykyä. Tutkimuksen edetessä määräävimmäksi tekijäksi työn tuloksia ajatellen nousi taloudellinen standardisointiaste, jonka perusteeksi kohdeyrityksen putkilämmönvaihtimissa valittiin standardisoitu suunnitteluohjeisto. Suunnitteluohjeiston lisäksi tuotestandardiin liitettiin ehdotus toimintapojen automatisoimisesta parametrisen valintataulukon sekä sähköisen manuaalin avulla. Työssä tuotestandardisointia tarkastellaan rationalisointi-investointina ja sen kannattavuutta on tarkasteltu investointilaskelmien sekä herkkyys analyysin avulla. Alkuarvot investointilaskelmissa perustuvat työssä yhdelle putkilämmönvaihtimelletoteutettuun tuotestandardisointitapaan, jos tätä tapaa käytettäisiin jatkossa kaikkiin kyseisen yksikön putkilämmönvaihtimiin.

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.

Modeling a Network of Heat Exchangers

This research work addresses the problem of building a mathematical model for the given system of heat exchangers and to determine the temperatures, pressures and velocities at the intermediate positions. Such model could be used in nding an optimal design for such a superstructure. To limit the size and computing time a reduced network model was used. The method can be generalized to larger network structures. A mathematical model which includes a system of non-linear equations has been built and solved according to the Newton-Raphson algorithm. The results obtained by the proposed mathematical model were compared with the results obtained by the Paterson approximation and Chen's Approximation. Results of this research work in collaboration with a current ongoing research at the department will optimize the valve positions and hence, minimize the pumping cost and maximize the heat transfer of the system of heat exchangers.

Analysis of heat transfer coefficient in silica nanofluids with water as base fluid under transitional flow

The effect of Reynolds number variation in a vertical double pipe counterflow heat exchanger due to the changes in viscosity can cause the change in flow regime, for instance, when heats up and cools down, it can convert from turbulent to laminar or inversely, that can have significant effect on heat transfer coefficient and pressure drop. Mainly, the range of transition phase has been studied in this study with the investigation of silica nanofluid dispersed in water in three different concentrations. The results have been compared with distilled water sample and showed a remarkable raise in heat transfer coefficient while pressure drop has been increased respectively, as well. Although pumping power has to go up at the same time and it is a drawback, heat transfer efficiency grows for diluted samples. On the other hand, for the most concentrated sample, effect of pressure drop dominates which leads to decline in the overall efficiency.

Vaneri- ja kertopuutehtaiden kuivatuskoneiden lämmöntalteenottojen järjestelmäselvitys

Vanerin tai kertopuun valmistusprosessissaviilun kuivaukseen käytetään suurin osa koko valmistusprosessin primäärienergiasta. Viilunkuivauskoneessa viilun sisältämä vesi siirretään tyypillisesti prosessihöyryllä lämmitettyyn viilunkuivaajan kiertoilmaan höyrystämällä ja poistetaanviilunkuivaajasta poistoilman mukana. Viilunkuivaajan poistoilma on lämmintä jaerittäin suuren kosteuspitoisuutensa takia sisältää runsaasti energiaa. Tyypillisellä viilunkuivaajalla poistoilmaan sitoutunut lämpöteho vaihtelee prosessiolosuhteista riippuen välillä 2,7-5,7 MW. Diplomityössä tutkittiin viilunkuivaajan poistoilman sisältämän lämmön talteenottoa laitteistolla, johon kuuluu lämmöntalteenottopesuri, jossa poistoilmalla lämmitetään tuotantolaitoksen tukkipuun hautomon kiertovettä sekä ilma-ilma-lämmönsiirrin, jolla lämmitetään pesurista poistuvan ilman jäännöslämmöllä ulkoilmaa tehdassalin tuloilmakäyttöön. Työn tavoitteena oli kehittää lämmöntalteenottojärjestelmän suunnittelua, mitoitusta ja ajotapoja. Työssä analysoitiin teoreettisesti pesuria ja ilmalämmönsiirrintä, kehitettiin lämmöntalteenottopesurin simulointimenetelmä ja mitattiin toiminnassa olevia talteenottolaitteistoja. Tutkimuksessa todettiin lämmöntalteenottohyötysuhteen vaihtelevan lämmityskaudella välillä 50-70 %. Lämmöntalteenottolaitteiston pesurin veteen saatava teho riippuu ensisijaisesti viilunkuivaajan poistoilman lämpösisällöstä, joka on enimmäkseen kosteusriippuvainen ja ilmanvaihtoilmaan saatava teho ulkolämpö-tilan määräämästä tehontarpeesta. Pesurin vesijärjestelmän vaikutusmekanismit pesurin suorituskykyyn tunnistettiin ja niiden pohjalta annetaan suositukset mitoitukseen ja ajotapaan. Lämmöntalteenottolaitteiston lämpötehon tasapainottamiseen pesurin ja ilma-ilma-lämmönsiirtimen välillä mitoituksen avulla esitellään työkalut.

Huoneistokohtaisen ilmanvaihtokoneen lämmöntalteenoton ja äänenvaimentimien kehitys

Tämän työn tavoitteena oli selvittää erään ilmanvaihtokoneen lämmöntalteenottokennon vuosihyötysuhde japohtia voitaisiinko ilmastointikoneen energiatehokkuutta parantaa käyttämällä nykyisen ristivirtalämmönsiirtimen tilalla vastavirtalämmönsiirrintä tai poistoilmalämpöpumppua. Lisäksi tavoitteena oli mitata kahden ilmanvaihtokoneeseen liitettävän äänenvaimentimen vaimennuskyky ja pyrkiä parantamaan näiden vaimennusta. Työn teoriaosassa käsitellään yleisellä tasolla eri ilmanvaihtotapojen ja lämmöntalteenottomenetelmien toimintaperiaatteita ja soveltuvuuksiaeri tilanteisiin. Työn kokeellisessa osassa mitataan ilmanvaihtokoneen lämmöntalteenoton toiminta eri sääolosuhteissa ja lasketaan tulosten perusteella laitteelle lämpötila- ja entalpiavuosihyötysuhteet. Äänenvaimentimesta rakennetaan kaksi kehitysversiota, joiden mittaustuloksia verrataan alkuperäisen vaimentimen vaimennuskykyyn. Lämmöntalteenottoja vertailtaessa havaittiin, että suositeltavin kehitysvaihtoehto olisi vastavirtalämmönsiirtimeen siirtyminen sen yksinkertaisen rakenteen ja hyvän hyötysuhteen ansiosta.Äänenvaimentimien rakenteessa kannattaa mittausten ja pohdinnan perusteella siirtyä seinämillä jaettuun vaimentimeen.

Lämmöntalteenotto kosteasta poistoilmasta prosessiteollisuudessa

Tämän diplomityön päämääränä on tehdä prosessiteollisuuden tarpeisiin Excel-taulukkolaskentaohjelmassa toimiva putkilämmönsiirtimen mitoitusohjelma. Prosessiteollisuudessa lämmönvaihtimien toimintaympäristöt ja olosuhteet vaihtelevat merkittävästi ja niinpä jokaisen vaihtimen suunnittelu ja mitoitus on toteutettava tapauskohtaisesti. Työssä käsitellään rekuperatiivisen ristivirtaputkilämmönvaihtimen yleinen lämpötekninen mitoitus sisältäen putken sisäpinnalle tapahtuvan mahdollisen lauhtumisen. Mitoitettava vaihdinkoostuu pystysuorista putkista, joissa lämmin ja kostea ilma virtaa putkien sisäpuolella ja kylmä kuiva ilma vaippapuolella vaakasuoraan. Vaihdinmateriaalina käytetään ruostumatonta AISI 304 -tai haponkestävää AISI 316 terästä. Kuuman ilman tila vaihtelee tarkasteltavan kohteen mukaan. Paperiteollisuuden kuivausyksiköiltä poistuva ilma on usein lämmintä ja kosteaa, ja infrakuivaimilta poistuva ilma on kuumaa. Mitoitettavalle lämmönvaihtimelle tulevan kuuman ilman lämpötila tapauksesta riippuen voi vaihdella 30°C, maksimissaan +300°C:een saakka, vesisisällön ollessa välillä 0,010...0,200 kg/kg ki tai jopa tämän ylikin. Vaihtimen mitoitus perustuu energiataseyhtälöiden käyttöön. Laskennassa määritetään vaihtimen pintalämpötila sekä mahdollinen kostean ilman lauhtuminen putken sisäpinnalle. Lisäksi teoria käsittää molempien virtausten tilanmuutosten laskennan. Työssä on esitetty esimerkkilaskelma, jossa on laskettu ilma- kostea ilma lämmönsiirrinkonstruktio. Esimerkissä on tarkasteltu vaihtimen hyötysuhdetta, virtausten lämpö- ja kosteuskäyttäytymistä ulkoilman lämpötilan funktiona. Ohjelmasta saadaan tulostettua mitoitettavanvaihtimen dimensiot; putkien lukumäärät syvyys- ja pituussuunnassa sekä kokonaisputkilukumäärä, putkien väliset etäisyydet toisiinsa nähden sekä syvyys, että pituussuunnassa, putkipituus ja putken sisä- ja ulkohalkaisijat. Nämä tiedot suunnittelija itse syöttää ohjelmalle alkuarvoina. Laskettuna tietona ohjelma antaa molempien virtausten poistolämpötilat, kuuman ilman poistuvan absoluuttisen kosteuden, kondenssivesimäärän, vaihtimen tehon ja painehäviöt vaippa- ja putkipuolelle. Lisäksi ohjelma laskee kuuman ilman ominaisentalpiat vaihtimen sisään- ja ulostulossa. Tämä mahdollistaa ilman tilapisteiden piirtämisen Mollier-piirrokseen.

Putkilämmönsiirtimien mitoitustyökalujen kehittäminen

Diplomityön tarkoituksena on kehittää tietokoneohjelma putkilämmönsiirtimen vaippapuolen painehäviön laskemiseksi. Ohjelmalla voidaan varmistaa lämmönsiirtimen mitoitusvaiheessa, että vaippapuolen painehäviö ei ylitä sallittuja rajoja. Ohjelmatäydentää olemassa olevia mitoitusohjelmia. Tässä diplomityössä käsitellään ainoastaan höyryvoimalaitosprosesseissa käytettäviä putkilämmönsiirtimiä. Työn kirjallisessa osassa on selvitetty periaate höyryvoimalaitosprosessista ja siinä käytettävistä putkilämmönsiirtimistä sekä esitetty putkilämmönsiirtimien rakenne, yleinen suunnittelu ja lämpö- ja virtaustekninen mitoitus. Painehäviön laskennassa käytetyt ja lämpö- ja virtausteknistä mitoitusta käsittelevässä kappaleessa esitetyt yhtälöt perustuvat Bell-Delawaren menetelmään. Painehäviönlaskentaohjelma on toteutettu hyväksikäyttäen Microsoft Excel taulukkolaskentaa ja Visual Basic -ohjelmointikieltä. Painehäviön laskenta perustuu segmenttivälilevyillä varustetun putkilämmönsiirtimen vaippapuolen yksifaasivirtaukseen. Lämmönsiirtimen lauhdutinosan painehäviö oletetaan merkityksettömäksi, joten kokonaispainehäviö muodostuu höyryn- ja lauhteenjäähdyttimessä. Kehitetty ohjelma on suunniteltu erityisesti lauhteenjäähdyttimessä muodostuvan painehäviön laskentaan. Ohjelmalla laskettuja painehäviön arvoja on verrattu todellisesta lämmönsiirtimestä mitattuihin arvoihin. Lasketut arvotvastaavat hyvin mittaamalla saatuja, eikä tuloksissa ilmene mitään systemaattista virhettä. Ohjelma on valmis käytettäväksi putkilämmönsiirtimien mitoitustyökaluna. Diplomityön pohjalta on tehty ehdotukset ohjelman edelleen kehittämiseksi.

Levylämmönsiirtimen mitoitus ja optimointi erityisesti lämpöpumpuissa

Työn tavoitteena on laatia käsikirjamainen läpileikkaus levylämmönsiirtimen rakenteesta ja käytöstä lämpöpumppulaitteistoissa. Kylmäprosessin tarkastelun lisäksi on eri lähteistä haettu yhtälöitä lämmönsiirron ja painehäviönlaskentaan. Lähdeaineistona on käytetty lämmönsiirron oppikirjoja, joiden lisäksi on käyty läpi suuri joukko tieteellisen tutkimuksen julkaisemia tutkimusraportteja levylämmönsiirtimen mitoituksesta erilaisissa käyttökohteissa. Oppikirjoissa ei ole esitetty varsinaisesti levylämmönsiirtimen laskentamenetelmiä, vaan niissä esitetään lämmönsiirron perusyhtälöt. Varsinainen lämmönsiirtolevyprofiilin laskentaan perustuva lähdeaineisto on löytynyt lämmönsiirtoon erikoistuneista julkaisuista. Lämmönsiirto tapahtuu aina kuumemmasta kylmempään tilaan. Lämmönsiirto eri virtausaineiden välillä toteutetaan lämmönsiirtimien avulla. Lämmönsiirrintyyppejä on olemassa lukuisia, joista yksi yleisesti käytetty tyyppi on levylämmönsiirrin. Tässä konstruktiossa on mahdutettu paljon lämmönsiirtopintaa ulkomitoiltaan pieneen tilaan. Tämä siirrintyyppi on eduksi silloin, kun virtaavat aineet ovat puhtaita ja niillä ei ole likaavaa vaikutusta lämmönsiirtopinnoille. Lämpöpumpulla tarkoitetaan laitetta, jolla voidaan käyttää hyödyksi lämmönlähteen matalaa lämpötilatasoa nostamalla lämpötilatasoa kompressorin puristustyön avulla korkeampaan lämpötilatasoon. Lämpöpumppulaitteiston toiminta perustuu kylmäprosessin toimintaan. Kylmäprosessin läpikäynti auttaa lukijaa hahmottamaan, millainen prosessi on kysymyksessä ja mitä komponentteja liittyy kylmäprosessiin. Tässä diplomityössä esitetyt yhtälöt antavat suuntaa, millä tavalla levylämmönsiirtimien ominaisuuksia voidaan laskea ja mitkä tekijät vaikuttavat siirtimien mitoittamiseen. Tarkemmat vaihdinkohtaiset laskentakorrelaatiot muotoutuvat vasta sitten, kun valmis tuote on testattu laboratorio-olosuhteissa ja siitä on saatu lämpötila-, virtaus- ja painesuhteet selville. Tämän jälkeen voidaan mittaustuloksiin perustuen rakentaa matemaattinen malli, jolla laskennallisesti määritelläänvaihtimien ominaisuudet. Lisäksi on esitetty yleisiä tapoja, joilla voidaan määritellä lämmönsiirtimien lämpöpintoja.

Hitsausmenetelmät ja niiden kehittämismahdollisuudet painelaitteiden valmistuksessa

Työn tavoitteena oli painelaitteiden hitsaustuotannon kehittäminen. Keinoina olivat käytössä olevien hitsausprosessienkehittäminen tai niiden käytön tehostaminen ja mahdollisesti muiden hitsausprosessien, kuten plasma-, laser- tai hybridihitsauksen, käyttöönotto. Työn teoriaosassa käsiteltiin painelaitteiden valmistuksessa käytettävät hitsausprosessit ja niiden kehittämiskeinot, hitsauskustannuksen muodostuminen sekä miten näihin kustannuksiin voidaan vaikuttaa. Lisäksi tarkasteltiin painelaitteiden valmistuksessa käytettävien materiaalien hitsattavuutta ja laadunvarmistusta. Tutkimusosuudessa käytiin läpi valmistuksen nykytila ja tehtiin tuotannon analysointi kehittämiskohteeksi valitulle lämmönsiirtimelle. Kehitysehdotuksenatuotannon tehostamiseksi ehdotettiin MAG-täytelanka- ja jauhekaarihitsauksen käytön lisäämistä. Lisäksi ehdotettiin juurituen käyttöä ja hitsauksen mekanisoinnin lisäämistä. Plasma- tai hybridihitsauksen käytöllä todettiin pystyttävän lyhentämään läpimenoaikoja. Ennen laitteistohankintojen tekoa täytyy tehdä lisäselvityksiä soveltuvuuden ja kannattavuuden varmistamiseksi.

Selective catalytic reduction of nitrogen oxides with ammonia in forced unsteady state reactors - Case based reasoning and mathematical model simulation reasoning

The application of forced unsteady-state reactors in case of selective catalytic reduction of nitrogen oxides (NOx) with ammonia (NH3) is sustained by the fact that favorable temperature and composition distributions which cannot be achieved in any steady-state regime can be obtained by means of unsteady-state operations. In a normal way of operation the low exothermicity of the selective catalytic reduction (SCR) reaction (usually carried out in the range of 280-350°C) is not enough to maintain by itself the chemical reaction. A normal mode of operation usually requires supply of supplementary heat increasing in this way the overall process operation cost. Through forced unsteady-state operation, the main advantage that can be obtained when exothermic reactions take place is the possibility of trapping, beside the ammonia, the moving heat wave inside the catalytic bed. The unsteady state-operation enables the exploitation of the thermal storage capacity of the catalyticbed. The catalytic bed acts as a regenerative heat exchanger allowing auto-thermal behaviour when the adiabatic temperature rise is low. Finding the optimum reactor configuration, employing the most suitable operation model and identifying the reactor behavior are highly important steps in order to configure a proper device for industrial applications. The Reverse Flow Reactor (RFR) - a forced unsteady state reactor - corresponds to the above mentioned characteristics and may be employed as an efficient device for the treatment of dilute pollutant mixtures. As a main disadvantage, beside its advantages, the RFR presents the 'wash out' phenomena. This phenomenon represents emissions of unconverted reactants at every switch of the flow direction. As a consequence our attention was focused on finding an alternative reactor configuration for RFR which is not affected by the incontrollable emissions of unconverted reactants. In this respect the Reactor Network (RN) was investigated. Its configuration consists of several reactors connected in a closed sequence, simulating a moving bed by changing the reactants feeding position. In the RN the flow direction is maintained in the same way ensuring uniformcatalyst exploitation and in the same time the 'wash out' phenomena is annulated. The simulated moving bed (SMB) can operate in transient mode giving practically constant exit concentration and high conversion levels. The main advantage of the reactor network operation is emphasizedby the possibility to obtain auto-thermal behavior with nearly uniformcatalyst utilization. However, the reactor network presents only a small range of switching times which allow to reach and to maintain an ignited state. Even so a proper study of the complex behavior of the RN may give the necessary information to overcome all the difficulties that can appear in the RN operation. The unsteady-state reactors complexity arises from the fact that these reactor types are characterized by short contact times and complex interaction between heat and mass transportphenomena. Such complex interactions can give rise to a remarkable complex dynamic behavior characterized by a set of spatial-temporal patterns, chaotic changes in concentration and traveling waves of heat or chemical reactivity. The main efforts of the current research studies concern the improvement of contact modalities between reactants, the possibility of thermal wave storage inside the reactor and the improvement of the kinetic activity of the catalyst used. Paying attention to the above mentioned aspects is important when higher activity even at low feeding temperatures and low emissions of unconverted reactants are the main operation concerns. Also, the prediction of the reactor pseudo or steady-state performance (regarding the conversion, selectivity and thermal behavior) and the dynamicreactor response during exploitation are important aspects in finding the optimal control strategy for the forced unsteady state catalytic tubular reactors. The design of an adapted reactor requires knowledge about the influence of its operating conditions on the overall process performance and a precise evaluation of the operating parameters rage for which a sustained dynamic behavior is obtained. An apriori estimation of the system parameters result in diminution of the computational efforts. Usually the convergence of unsteady state reactor systems requires integration over hundreds of cycles depending on the initial guess of the parameter values. The investigation of various operation models and thermal transfer strategies give reliable means to obtain recuperative and regenerative devices which are capable to maintain an auto-thermal behavior in case of low exothermic reactions. In the present research work a gradual analysis of the SCR of NOx with ammonia process in forced unsteady-state reactors was realized. The investigation covers the presentationof the general problematic related to the effect of noxious emissions in the environment, the analysis of the suitable catalysts types for the process, the mathematical analysis approach for modeling and finding the system solutions and the experimental investigation of the device found to be more suitable for the present process. In order to gain information about the forced unsteady state reactor design, operation, important system parameters and their values, mathematical description, mathematicalmethod for solving systems of partial differential equations and other specific aspects, in a fast and easy way, and a case based reasoning (CBR) approach has been used. This approach, using the experience of past similarproblems and their adapted solutions, may provide a method for gaining informations and solutions for new problems related to the forced unsteady state reactors technology. As a consequence a CBR system was implemented and a corresponding tool was developed. Further on, grooving up the hypothesis of isothermal operation, the investigation by means of numerical simulation of the feasibility of the SCR of NOx with ammonia in the RFRand in the RN with variable feeding position was realized. The hypothesis of non-isothermal operation was taken into account because in our opinion ifa commercial catalyst is considered, is not possible to modify the chemical activity and its adsorptive capacity to improve the operation butis possible to change the operation regime. In order to identify the most suitable device for the unsteady state reduction of NOx with ammonia, considering the perspective of recuperative and regenerative devices, a comparative analysis of the above mentioned two devices performance was realized. The assumption of isothermal conditions in the beginningof the forced unsteadystate investigation allowed the simplification of the analysis enabling to focus on the impact of the conditions and mode of operation on the dynamic features caused by the trapping of one reactant in the reactor, without considering the impact of thermal effect on overall reactor performance. The non-isothermal system approach has been investigated in order to point out the important influence of the thermal effect on overall reactor performance, studying the possibility of RFR and RN utilization as recuperative and regenerative devices and the possibility of achieving a sustained auto-thermal behavior in case of lowexothermic reaction of SCR of NOx with ammonia and low temperature gasfeeding. Beside the influence of the thermal effect, the influence of the principal operating parameters, as switching time, inlet flow rate and initial catalyst temperature have been stressed. This analysis is important not only because it allows a comparison between the two devices and optimisation of the operation, but also the switching time is the main operating parameter. An appropriate choice of this parameter enables the fulfilment of the process constraints. The level of the conversions achieved, the more uniform temperature profiles, the uniformity ofcatalyst exploitation and the much simpler mode of operation imposed the RN as a much more suitable device for SCR of NOx with ammonia, in usual operation and also in the perspective of control strategy implementation. Theoretical simplified models have also been proposed in order to describe the forced unsteady state reactors performance and to estimate their internal temperature and concentration profiles. The general idea was to extend the study of catalytic reactor dynamics taking into account the perspectives that haven't been analyzed yet. The experimental investigation ofRN revealed a good agreement between the data obtained by model simulation and the ones obtained experimentally.