Oxidation rates of carbon and nitrogen in char residues from solid fuels

Computational fluid dynamics (CFD) modeling is an important tool in designing new combustion systems. By using CFD modeling, entire combustion systems can be modeled and the emissions and the performance can be predicted. CFD modeling can also be used to develop new and better combustion systems from an economical and environmental point of view. In CFD modeling of solid fuel combustion, the combustible fuel is generally treated as single fuel particles. One of the limitations with the CFD modeling concerns the sub-models describing the combustion of single fuel particles. Available models in the scientific literature are in many cases not suitable as submodels for CFD modeling since they depend on a large number of input parameters and are computationally heavy. In this thesis CFD-applicable models are developed for the combustion of single fuel particles. The single particle models can be used to improve the combustion performance in various combustion devices or develop completely new technologies. The investigated fields are oxidation of carbon (C) and nitrogen (N) in char residues from solid fuels. Modeled char-C oxidation rates are compared to experimental oxidation rates for a large number of pulverized solid fuel chars under relevant combustion conditions. The experiments have been performed in an isothermal plug flow reactor operating at 1123-1673 K and 3-15 vol.% O2. In the single particle model, the char oxidation is based on apparent kinetics and depends on three fuel specific parameters: apparent pre-exponential factor, apparent activation energy, and apparent reaction order. The single particle model can be incorporated as a sub-model into a CFD code. The results show that the modeled char oxidation rates are in good agreement with experimental char oxidation rates up to around 70% of burnout. Moreover, the results show that the activation energy and the reaction order can be assumed to be constant for a large number of bituminous coal chars under conditions limited by the combined effects of chemical kinetics and pore diffusion. Based on this, a new model based on only one fuel specific parameter is developed (Paper III). The results also show that reaction orders of bituminous coal chars and anthracite chars differ under similar conditions (Paper I and Paper II); reaction orders of bituminous coal chars were found to be one, while reaction orders of anthracite chars were determined to be zero. This difference in reaction orders has not previously been observed in the literature and should be considered in future char oxidation models. One of the most frequently used comprehensive char oxidation models could not explain the difference in the reaction orders. In the thesis (Paper II), a modification to the model is suggested in order to explain the difference in reaction orders between anthracite chars and bituminous coal chars. Two single particle models are also developed for the NO formation and reduction during the oxidation of single biomass char particles. In the models the char-N is assumed to be oxidized to NO and the NO is partly reduced inside the particle. The first model (Paper IV) is based on the concentration gradients of NO inside and outside the particle and the second model is simplified to such an extent that it is based on apparent kinetics and can be incorporated as a sub-model into a CFD code (Paper V). Modeled NO release rates from both models were in good agreement with experimental measurements from a single particle reactor of quartz glass operating at 1173-1323 K and 3-19 vol.% O2. In the future, the models can be used to reduce NO emissions in new combustion systems.

Häiriöanalyysejä tehonkorotetun kiehutusvesivoimalaitoksen käytettävyyden varmistamiseksi

TVO suunnittelee reaktoritehon 10 %:n korotusta Olkiluoto 1 ja 2 -voimalaitoksille. Reaktoriteho nostetaan 2500 MW:sta 2750 MW:iin polttoaineen rikastusastetta nostamalla ja pääkiertovirtausta kasvattamalla. Samalla syöttövesivirtaus reaktoriin ja tuorehöyryvirtaus turpiineille kasvaa. Lauhteenpuhdistusjärjestelmän kapasiteettia ei voida kuitenkaan kasvattaa, joten massavirran lisäys toteutetaan ottamalla käyttöön korkeapainesivulauhteen eteenpäinpumppaus. Lauhteen esilämmityslinjojen, lauhteenpuhdistuksen ja syöttövesipumppujen massavirta säilyy siten nykyisellään. Muita merkittäviä tehonkorotukseen liittyviä laitosmuutoksia ovat pääkiertopumppujen uusinta ja korkeapaineturpiinin muutokset. Tehonkorotetun prosessin käytettävyyden varmistamiseksi tehdään häiriöanalyysejä Apros-prosessisimulointiohjelmistoa käyttäen. OL1 ja OL2 -laitoksista on olemassa validoitu 2500 MW:n laitosmalli, josta muokatulla 2750 MW:n laitosmallilla simuloinnit tehdään. Häiriöanalyysien avulla selvitetään säätöjärjestelmien kyky pitää prosessin tila hallinnassa ilman suojausautomaation laukeamista. Simuloituihin tapauksiin kuuluu pumppujen ja venttiilien vikaantumistapauksia sekä turpiini- ja reaktoripuolen pikasulku- ja osittaispikasulkutapauksia. Myös meriveden lämpötilan vaikutusta häiriötilanteisiin tarkastellaan. Analyysien perusteella voimalaitosten ohjaus- ja suojausautomaatio toimivat hyvin myös korotetulla teholla. Tehonkorotuksen jälkeiset suuremmat massavirrat aiheuttavat kuitenkin voimakkaampia reaktoripaineen ja -tehon vaihteluita varsinkin venttiilien sulkeutumistapauksissa. Simuloinnit osoittivat, että tehonkorotus 2500 MW:sta 2750 MW:iin on mahdollinen, mutta aiheuttaa pieniä muutoksia laitoksen suojausjärjestelmien laukaisurajoihin.

Analyzing hydraulic head of fluid flow in underground porous medium

Hydraulic head is distributed through a medium with porous aspect. The analysis of hydraulic head from one point to another is used by the Richard's equation. This equation is equivalent to the groundwater ow equation that predicts the volumetric water contents. COMSOL 3.5 is used for computation applying Richard's equation. A rectangle of 100 meters of length and 10 meters of large (depth) with 0,1 m/s fl ux of inlet as source of our fl uid is simulated. The domain have Richards' equation model in two dimension (2D). Hydraulic head increases proportional with moisture content.

Computational studies for the design of process equipment with complex geometries

This study combines several projects related to the flows in vessels with complex shapes representing different chemical apparata. Three major cases were studied. The first one is a two-phase plate reactor with a complex structure of intersecting micro channels engraved on one plate which is covered by another plain plate. The second case is a tubular microreactor, consisting of two subcases. The first subcase is a multi-channel two-component commercial micromixer (slit interdigital) used to mix two liquid reagents before they enter the reactor. The second subcase is a micro-tube, where the distribution of the heat generated by the reaction was studied. The third case is a conventionally packed column. However, flow, reactions or mass transfer were not modeled. Instead, the research focused on how to describe mathematically the realistic geometry of the column packing, which is rather random and can not be created using conventional computeraided design or engineering (CAD/CAE) methods. Several modeling approaches were used to describe the performance of the processes in the considered vessels. Computational fluid dynamics (CFD) was used to describe the details of the flow in the plate microreactor and micromixer. A space-averaged mass transfer model based on Fick’s law was used to describe the exchange of the species through the gas-liquid interface in the microreactor. This model utilized data, namely the values of the interfacial area, obtained by the corresponding CFD model. A common heat transfer model was used to find the heat distribution in the micro-tube. To generate the column packing, an additional multibody dynamic model was implemented. Auxiliary simulation was carried out to determine the position and orientation of every packing element in the column. This data was then exported into a CAD system to generate desirable geometry, which could further be used for CFD simulations. The results demonstrated that the CFD model of the microreactor could predict the flow pattern well enough and agreed with experiments. The mass transfer model allowed to estimate the mass transfer coefficient. Modeling for the second case showed that the flow in the micromixer and the heat transfer in the tube could be excluded from the larger model which describes the chemical kinetics in the reactor. Results of the third case demonstrated that the auxiliary simulation could successfully generate complex random packing not only for the column but also for other similar cases.

Modeling of Liquid-Solid Flow in Industrial Scale

In the present work, liquid-solid flow in industrial scale is modeled using the commercial software of Computational Fluid Dynamics (CFD) ANSYS Fluent 14.5. In literature, there are few studies on liquid-solid flow in industrial scale, but any information about the particular case with modified geometry cannot be found. The aim of this thesis is to describe the strengths and weaknesses of the multiphase models, when a large-scale application is studied within liquid-solid flow, including the boundary-layer characteristics. The results indicate that the selection of the most appropriate multiphase model depends on the flow regime. Thus, careful estimations of the flow regime are recommended to be done before modeling. The computational tool is developed for this purpose during this thesis. The homogeneous multiphase model is valid only for homogeneous suspension, the discrete phase model (DPM) is recommended for homogeneous and heterogeneous suspension where pipe Froude number is greater than 1.0, while the mixture and Eulerian models are able to predict also flow regimes, where pipe Froude number is smaller than 1.0 and particles tend to settle. With increasing material density ratio and decreasing pipe Froude number, the Eulerian model gives the most accurate results, because it does not include simplifications in Navier-Stokes equations like the other models. In addition, the results indicate that the potential location of erosion in the pipe depends on material density ratio. Possible sedimentation of particles can cause erosion and increase pressure drop as well. In the pipe bend, especially secondary flows, perpendicular to the main flow, affect the location of erosion.

CTQ process flow development and the impact for cost of quality

Tämä työ tehtiin Kone Industrial Oy:lle Major Projects yksikköön, laatuosastolle. Kone Major Projects yksikkö keskittyy erikoisiin ja suuriin hissi- ja liukuporras projekteihin. Työn tavoitteena oli luoda harmonisoitu prosessi hissikomponenttien laaduntarkkailua varten sekä tarkastella ja vertailla kustannussäästöjä, jota tällä uudella prosessilla voidaan saavuttaa. Tavoitteena oli saavuttaa 80-prosentin kustannussäästöt laatukustannuksissa uuden laatuprosessin avulla. Työn taustana ja tutkimusongelmana ovat lisääntyneet erikoisprojektit ja niiden myötä lisääntynyt laaduntarkkailun tarve. Ongelmana laaduntarkkailussa voitiin pitää harmonisoidun ja selkeän prosessin puuttumista C-prosessikomponenttien valmistuksessa. Lisäksi kehitysprosessin aikana luotiin vanhojen työkalujen pohjalta keskeinen laaduntarkkailutyökalu, CTQ-työkalu. Työssä käsitellään ensin Konetta yhtiönä ja selvitetään Koneen keskeisimmät prosessit työn taustaksi. Teoria osuudessa käsitellään prosessin kehitykseen liittyviä teorioita sekä yleisiä laatukäsitteitä ja esitetään teorioita laadun asemasta nykypäivänä. Lopuksi käsitellään COQ eli laatukustannusten teoriaa ja esitellään teoria PAF-analyysille, jota käytetään työssä laatukustannusten vertailuun case esimerkin avulla. Työssä kuvataan CTQ prosessin luominen alusta loppuun ja case esimerkin avulla testataan uutta CTQ prosessia pilottihankkeessa. Tässä case esimerkissä projektin bracket eli johdekiinnitysklipsi tuotetaan uuden laatuprosessin avulla sekä tehdään kustannusvertailu saman projektin toisen bracketin kanssa, joka on tuotettu ennen uuden laatuprosessin implementoimista. Työn lopputuloksena CTQ prosessi saatiin luotua ja sitä pystyttiin testaamaan käytännössä case esimerkin avulla. Tulosten perusteella voidaan sanoa, että CTQ prosessin käyttö vähentää laatukustannuksia huomattavasti ja helpottaa laadunhallintaa C-prosessikomponenttien tuotannossa.

Microcracks and Vortices in Superconducting Thin Films

In this work, superconducting YBa₂ Cu₃O_6+x (YBCO) thin films have been studied with the experimental focus on the anisotropy of BaZrO₃ (BZO) doped YBCOthin films and the theoretical focus on modelling flux pinning by numerically solving Ginzburg- Landau equations. Also, the structural properties of undoped YBCO thin films grown on NdGaO₃ (NGO) and MgO substrates were investigated. The thin film samples were made by pulsed laser ablation on single crystal substrates. The structural properties of the thin films were characterized by X-ray diffraction and atomic force microscope measurements. The superconducting properties were investigated with a magnetometer and also with transport measurements in pulsed magnetic field up to 30 T. Flux pinning was modelled by restricting the value of the order parameter inside the columnar pinning sites and then solving the Ginzburg-Landau equations numerically with the restrictions in place. The computations were done with a parallel code on a supercomputer. The YBCO thin films were seen to develop microcracks when grown on NGO or MgO substrates. The microcrack formation was connected to the structure of the YBCO thin films in both cases. Additionally, the microcracks can be avoided by careful optimization of the deposition parameters and the film thickness. The BZO doping of the YBCO thin films was seen to decrease the effective electron mass anisotropy, which was seen by fitting the Blatter scaling to the angle dependence of the upper critical field. The Ginzburg-Landau simulations were able to reproduce the measured magnetic field dependence of the critical current density for BZO doped and undoped YBCO. The simulations showed that in addition to the large density also the large size of the BZO nanorods is a key factor behind the change in the power law behaviour between BZO doped and undoped YBCO. Additionally, the Ginzburg-Landau equations were solved for type I thin films where giant vortices were seen to appear depending on the film thickness. The simulations predicted that singly quantized vortices are stable in type I films up to quite large thicknesses and that the size of the vortices increases with decreasing film thickness, in a way that is similar to the behaviour of the interaction length of Pearl vortices.

Novel technology for preparation of optically active chemicals

Asymmetric synthesis using modified heterogeneous catalysts has gained lots of interest in the production of optically pure chemicals, such as pharmaceuticals, nutraceuticals, fragrances and agrochemicals. Heterogeneous modified catalysts capable of inducing high enantioselectivities are preferred in industrial scale due to their superior separation and handling properties. The topic has been intensively investigated both in industry and academia. The enantioselective hydrogenation of ethyl benzoylformate (EBF) to (R)-ethyl mandelate over (-)-cinchonidine (CD)-modified Pt/Al2O3 catalyst in a laboratory-scale semi-batch reactor was studied as a function of modifier concentration, reaction temperature, stirring rate and catalyst particle size. The main product was always (R)-ethyl mandelate while small amounts of (S)-ethyl mandelate were obtained as by product. The kinetic results showed higher enantioselectivity and lower initial rates approaching asymptotically to a constant value as the amount of modifier was increased. Additionally, catalyst deactivation due to presence of impurities in the feed was prominent in some cases; therefore activated carbon was used as a cleaning agent of the raw material to remove impurities prior to catalyst addition. Detailed characterizations methods (SEM, EDX, TPR, BET, chemisorption, particle size distribution) of the catalysts were carried out. Solvent effects were also studied in the semi-batch reactor. Solvents with dielectric constant (e) between 2 and 25 were applied. The enantiomeric excess (ee) increased with an increase of the dielectric coefficient up to a maximum followed by a nonlinear decrease. A kinetic model was proposed for the enantioselectivity dependence on the dielectric constant based on the Kirkwood treatment. The non-linear dependence of ee on (e) successfully described the variation of ee in different solvents. Systematic kinetic experiments were carried out in the semi-batch reactor. Toluene was used as a solvent. Based on these results, a kinetic model based on the assumption of different number of sites was developed. Density functional theory calculations were applied to study the energetics of the EBF adsorption on pure Pt(1 1 1). The hydrogenation rate constants were determined along with the adsorption parameters by non-linear regression analysis. A comparison between the model and the experimental data revealed a very good correspondence. Transient experiments in a fixed-bed reactor were also carried out in this work. The results demonstrated that continuous enantioselective hydrogenation of EBF in hexane/2-propanol 90/10 (v/v) is possible and that continuous feeding of (-)-cinchonidine is needed to maintain a high steady-state enantioselectivity. The catalyst showed a good stability and high enantioselectivity was achieved in the fixed-bed reactor. Chromatographic separation of (R)- and (S)-ethyl mandelate originating from the continuous reactor was investigated. A commercial column filled with a chiral resin was chosen as a perspective preparative-scale adsorbent. Since the adsorption equilibrium isotherms were linear within the entire investigated range of concentrations, they were determined by pulse experiments for the isomers present in a post-reaction mixture. Breakthrough curves were measured and described successfully by the dispersive plug flow model with a linear driving force approximation. The focus of this research project was the development of a new integrated production concept of optically active chemicals by combining heterogeneous catalysis and chromatographic separation technology. The proposed work is fundamental research in advanced process technology aiming to improve efficiency and enable clean and environmentally benign production of enantiomeric pure chemicals.

Methyl and ethyl chloride synthesis in microreactors

Methyl chloride is an important chemical intermediate with a variety of applications. It is produced today in large units and shipped to the endusers. Most of the derived products are harmless, as silicones, butyl rubber and methyl cellulose. However, methyl chloride is highly toxic and flammable. On-site production in the required quantities is desirable to reduce the risks involved in transportation and storage. Ethyl chloride is a smaller-scale chemical intermediate that is mainly used in the production of cellulose derivatives. Thus, the combination of onsite production of methyl and ethyl chloride is attractive for the cellulose processing industry, e.g. current and future biorefineries. Both alkyl chlorides can be produced by hydrochlorination of the corresponding alcohol, ethanol or methanol. Microreactors are attractive for the on-site production as the reactions are very fast and involve toxic chemicals. In microreactors, the diffusion limitations can be suppressed and the process safety can be improved. The modular setup of microreactors is flexible to adjust the production capacity as needed. Although methyl and ethyl chloride are important chemical intermediates, the literature available on potential catalysts and reaction kinetics is limited. Thus the thesis includes an extensive catalyst screening and characterization, along with kinetic studies and engineering the hydrochlorination process in microreactors. A range of zeolite and alumina based catalysts, neat and impregnated with ZnCl2, were screened for the methanol hydrochlorination. The influence of zinc loading, support, zinc precursor and pH was investigated. The catalysts were characterized with FTIR, TEM, XPS, nitrogen physisorption, XRD and EDX to identify the relationship between the catalyst characteristics and the activity and selectivity in the methyl chloride synthesis. The acidic properties of the catalyst were strongly influenced upon the ZnCl2 modification. In both cases, alumina and zeolite supports, zinc reacted to a certain amount with specific surface sites, which resulted in a decrease of strong and medium Brønsted and Lewis acid sites and the formation of zinc-based weak Lewis acid sites. The latter are highly active and selective in methanol hydrochlorination. Along with the molecular zinc sites, bulk zinc species are present on the support material. Zinc modified zeolite catalysts exhibited the highest activity also at low temperatures (ca 200 °C), however, showing deactivation with time-onstream. Zn/H-ZSM-5 zeolite catalysts had a higher stability than ZnCl2 modified H-Beta and they could be regenerated by burning the coke in air at 400 °C. Neat alumina and zinc modified alumina catalysts were active and selective at 300 °C and higher temperatures. However, zeolite catalysts can be suitable for methyl chloride synthesis at lower temperatures, i.e. 200 °C. Neat γ-alumina was found to be the most stable catalyst when coated in a microreactor channel and it was thus used as the catalyst for systematic kinetic studies in the microreactor. A binder-free and reproducible catalyst coating technique was developed. The uniformity, thickness and stability of the coatings were extensively characterized by SEM, confocal microscopy and EDX analysis. A stable coating could be obtained by thermally pretreating the microreactor platelets and ball milling the alumina to obtain a small particle size. Slurry aging and slow drying improved the coating uniformity. Methyl chloride synthesis from methanol and hydrochloric acid was performed in an alumina-coated microreactor. Conversions from 4% to 83% were achieved in the investigated temperature range of 280-340 °C. This demonstrated that the reaction is fast enough to be successfully performed in a microreactor system. The performance of the microreactor was compared with a tubular fixed bed reactor. The results obtained with both reactors were comparable, but the microreactor allows a rapid catalytic screening with low consumption of chemicals. As a complete conversion of methanol could not be reached in a single microreactor, a second microreactor was coupled in series. A maximum conversion of 97.6 % and a selectivity of 98.8 % were reached at 340°C, which is close to the calculated values at a thermodynamic equilibrium. A kinetic model based on kinetic experiments and thermodynamic calculations was developed. The model was based on a Langmuir Hinshelwood-type mechanism and a plug flow model for the microreactor. The influence of the reactant adsorption on the catalyst surface was investigated by performing transient experiments and comparing different kinetic models. The obtained activation energy for methyl chloride was ca. two fold higher than the previously published, indicating diffusion limitations in the previous studies. A detailed modeling of the diffusion in the porous catalyst layer revealed that severe diffusion limitations occur starting from catalyst coating thicknesses of 50 μm. At a catalyst coating thickness of ca 15 μm as in the microreactor, the conditions of intrinsic kinetics prevail. Ethanol hydrochlorination was performed successfully in the microreactor system. The reaction temperature was 240-340°C. An almost complete conversion of ethanol was achieved at 340°C. The product distribution was broader than for methanol hydrochlorination. Ethylene, diethyl ether and acetaldehyde were detected as by-products, ethylene being the most dominant by-product. A kinetic model including a thorough thermodynamic analysis was developed and the influence of adsorbed HCl on the reaction rate of ethanol dehydration reactions was demonstrated. The separation of methyl chloride using condensers was investigated. The proposed microreactor-condenser concept enables the production of methyl chloride with a high purity of 99%.

Loistehon säätö ja kompensointi

Kantaverkkoyhtiö Fingrid Oy on tarkastellut kantaverkosta otettua tai tuotettua loistehoa alueittain Fingridin määrittelemällä maantieteellisellä tavalla. Alue koostuu liittymispisteistä, jotka ovat yksittäisten verkkoyhtiöiden liittymispisteitä. Yhden alueen sisällä voi siis olla monia eri verkkoyhtiöiden liittymispisteitä. Vuodesta 2016 alkaen loistehon tarkastelu vaihtuu aluetarkastelusta liittymispistetarkasteluun. Liittymispistetarkasteluun siirtymisen myötä tulee verkkoyhtiöiden investoida erinäisiin kompensointiratkaisuihin, jotta yksittäisten liittymispisteiden loistehot pysyvät Fingridin asettaman loistehon tarkasteluikkunan sisällä. Työssä perehdytään erinäisiin kompensointiratkaisuihin ja niiden kytkentöihin sekä kytkentöjen vaikutuksiin olemassa olevaan sähköverkkoon. Lisäksi työssä käsitellään tuulivoimaloiden hyödyntämistä loistehon säädössä. Esimerkkitapauksessa tarkastellaan erästä Caruna Oy:n verkkoa ja pohditaan teknistä ja mahdollisimman taloudellista kompensointiratkaisua kyseiseen verkkoon. Työn lopputulos on, että maakaapeloinnin määrän kasvaessa loistehon kompensoinnin tarve kasvaa. Näin ollen Fingrid Oy:ltä tarvitaan selkeitä linjauksia moneen ongelmakohtaan, kuten tuotanto-kulutuspisteen loistehon tarkasteluun. Lisäksi Energiavirastolta tarvitaan kannusteita loistehon kompensointiin, sillä verkkoyhtiöiden liiketoiminta perustuu regulaatiomalliin ja kompensointilaitteiston yksikköhinnasto on hyvin suppea.