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.

The impact of bioavailability limitations on microbial stable isotope fractionation in a multiphase system

Stable isotope fractionation analysis of contaminants is a promising method for assessing biodegradation of contaminants in natural systems. However, standard procedures to determine stable isotope fractionation factors, so far, neglect the influence of pollutant bioavailability on stable isotope fractionation. On a microscale, bioavailability may vary due to the spatio-temporal variability of local contaminant concentrations, limited effective diffusivities of the contaminants and cell densities, and thus, the pollutant supply might not meet the intrinsic degradation capacity of the microorganisms. The aim of this study was to demonstrate the effect of bioavailability on the apparent stable isotope fractionation, using a multiphase laboratory setup. The data gained show that the apparent isotope fractionation factors observed during biodegradation processes depend on the amount of biomass and/or the rate of toluene mass transfer from a second to the aqueous phase. They indicate that physico-chemical processes need to be taken into account when stable isotope fractionation analysis is used for the quantification of environmental contaminant degradation.

Reaction Kinetics and Viscosity Modelling in the Fusion Syntheses of Ca- and Ca/Mg-resinates

Rosin is a natural product from pine forests and it is used as a raw material in resinate syntheses. Resinates are polyvalent metal salts of rosin acids and especially Ca- and Ca/Mg- resinates find wide application in the printing ink industry. In this thesis, analytical methods were applied to increase general knowledge of resinate chemistry and the reaction kinetics was studied in order to model the non linear solution viscosity increase during resinate syntheses by the fusion method. Solution viscosity in toluene is an important quality factor for resinates to be used in printing inks. The concept of critical resinate concentration, c crit, was introduced to define an abrupt change in viscosity dependence on resinate concentration in the solution. The concept was then used to explain the non-inear solution viscosity increase during resinate syntheses. A semi empirical model with two estimated parameters was derived for the viscosity increase on the basis of apparent reaction kinetics. The model was used to control the viscosity and to predict the total reaction time of the resinate process. The kinetic data from the complex reaction media was obtained by acid value titration and by FTIR spectroscopic analyses using a conventional calibration method to measure the resinate concentration and the concentration of free rosin acids. A multivariate calibration method was successfully applied to make partial least square (PLS) models for monitoring acid value and solution viscosity in both mid-infrared (MIR) and near infrared (NIR) regions during the syntheses. The calibration models can be used for on line resinate process monitoring. In kinetic studies, two main reaction steps were observed during the syntheses. First a fast irreversible resination reaction occurs at 235 °C and then a slow thermal decarboxylation of rosin acids starts to take place at 265 °C. Rosin oil is formed during the decarboxylation reaction step causing significant mass loss as the rosin oil evaporates from the system while the viscosity increases to the target level. The mass balance of the syntheses was determined based on the resinate concentration increase during the decarboxylation reaction step. A mechanistic study of the decarboxylation reaction was based on the observation that resinate molecules are partly solvated by rosin acids during the syntheses. Different decarboxylation mechanisms were proposed for the free and solvating rosin acids. The deduced kinetic model supported the analytical data of the syntheses in a wide resinate concentration region, over a wide range of viscosity values and at different reaction temperatures. In addition, the application of the kinetic model to the modified resinate syntheses gave a good fit. A novel synthesis method with the addition of decarboxylated rosin (i.e. rosin oil) to the reaction mixture was introduced. The conversion of rosin acid to resinate was increased to the level necessary to obtain the target viscosity for the product at 235 °C. Due to a lower reaction temperature than in traditional fusion synthesis at 265 °C, thermal decarboxylation is avoided. As a consequence, the mass yield of the resinate syntheses can be increased from ca. 70% to almost 100% by recycling the added rosin oil.

Membrane emulsification in preparation of single site catalyst

In the theory part the membrane emulsification was studied. Emulsions are used in many industrial areas. Traditionally emulsions are prepared by using high shear in rotor-stator systems or in high pressure homogenizer systems. In membrane emulsification two immiscible liquids are mixed by pressuring one liquid through the membrane into the other liquid. With this technique energy could be saved, more homogeneous droplets could be formed and the amount of surfactant could be decreased. Ziegler-Natta and single-site catalysts are used in olefin polymerization processes. Nowadays, these catalysts are prepared according to traditional mixing emulsification. More homogeneous catalyst particles that have narrower particle size distribution might be prepared with membrane emulsification. The aim of the experimental part was to examine the possibility to prepare single site polypropylene catalyst using membrane emulsification technique. Different membrane materials and solidification techniques of the emulsion were examined. Also the toluene-PFC phase diagram was successfully measured during this thesis work. This phase diagram was used for process optimization. The polytetrafluoroethylene membranes had the largest contact angles with toluene and also the biggest difference between the contact angles measured with PFC and toluene. Despite of the contact angle measurement results no significant difference was noticed between particles prepared using PTFE membrane or metal sinter. The particle size distributions of catalyst prepared in these tests were quite wide. This would probably be fixed by using a membrane with a more homogeneous pore size distribution. It is also possible that the solidification rate has an effect on the particle sizes and particle morphology. When polymeric membranes are compared PTFE is probably still the best material for the process as it had the best chemical durability.

Kinetic modeling of mechanisms of industrially important organic reactions in gas and liquid phase

This dissertation is based on 5 articles which deal with reaction mechanisms of the following selected industrially important organic reactions: 1. dehydrocyclization of n-butylbenzene to produce naphthalene 2. dehydrocyclization of 1-(p-tolyl)-2-methylbutane (MB) to produce 2,6-dimethylnaphthalene 3. esterification of neopentyl glycol (NPG) with different carboxylic acids to produce monoesters 4. skeletal isomerization of 1-pentene to produce 2-methyl-1-butene and 2-methyl-2-butene The results of initial- and integral-rate experiments of n-butylbenzene dehydrocyclization over selfmade chromia/alumina catalyst were applied when investigating reaction 2. Reaction 2 was performed using commercial chromia/alumina of different acidity, platina on silica and vanadium/calcium/alumina as catalysts. On all catalysts used for the dehydrocyclization, major reactions were fragmentation of MB and 1-(p-tolyl)-2-methylbutenes (MBes), dehydrogenation of MB, double bond transfer, hydrogenation and 1,6-cyclization of MBes. Minor reactions were 1,5-cyclization of MBes and methyl group fragmentation of 1,6- cyclization products. Esterification reactions of NPG were performed using three different carboxylic acids: propionic, isobutyric and 2-ethylhexanoic acid. Commercial heterogeneous gellular (Dowex 50WX2), macroreticular (Amberlyst 15) type resins and homogeneous para-toluene sulfonic acid were used as catalysts. At first NPG reacted with carboxylic acids to form corresponding monoester and water. Then monoester esterified with carboxylic acid to form corresponding diester. In disproportionation reaction two monoester molecules formed NPG and corresponding diester. All these three reactions can attain equilibrium. Concerning esterification, water was removed from the reactor in order to prevent backward reaction. Skeletal isomerization experiments of 1-pentene were performed over HZSM-22 catalyst. Isomerization reactions of three different kind were detected: double bond, cis-trans and skeletal isomerization. Minor side reaction were dimerization and fragmentation. Monomolecular and bimolecular reaction mechanisms for skeletal isomerization explained experimental results almost equally well. Pseudohomogeneous kinetic parameters of reactions 1 and 2 were estimated by usual least squares fitting. Concerning reactions 3 and 4 kinetic parameters were estimated by the leastsquares method, but also the possible cross-correlation and identifiability of parameters were determined using Markov chain Monte Carlo (MCMC) method. Finally using MCMC method, the estimation of model parameters and predictions were performed according to the Bayesian paradigm. According to the fitting results suggested reaction mechanisms explained experimental results rather well. When the possible cross-correlation and identifiability of parameters (Reactions 3 and 4) were determined using MCMC method, the parameters identified well, and no pathological cross-correlation could be seen between any parameter pair.

Electroactive ion exchange membranes based on conducting polymers

Ion exchange membranes are indispensable for the separation of ionic species. They can discriminate between anions and cations depending on the type of fixed ionic group present in the membrane. These conventional ion exchange membranes (CIX) have exceptional ionic conductivity, which is advantageous in various electromembrane separation processes such as electrodialysis, electrodeionisation and electrochemical ion exchange. The main disadvantage of CIX membranes is their high electrical resistance owing to the fact that the membranes are electronically non conductive. An alternative can be electroactive ion exchange membranes, which are ionically and electronically conducting. Polypyrrole (PPy) is a type of electroactive ion exchange material as well as a commonly known conducting polymer. When PPy membranes are repeatedly reduced and oxidised, ions are pumped through the membrane. The main aim of this thesis was to develop electroactive cation transport membranes based on PPy for the selective transport of divalent cations. Membranes developed composed of PPy films deposited on commercially available support materials. To carry out this study, cation exchange membranes based on PPy doped with immobile anions were prepared. Two types of dopant anions known to interact with divalent metal ions were considered, namely 4-sulphonic calix[6]arene (C6S) and carboxylated multiwalled carbon nanotubes (CNT). The transport of ions across membranes containing PPy doped with polystyrene sulphonate (PSS) and PPy doped with para-toluene sulphonate (pTS) was also studied in order to understand the nature of ion transport and permeability across PPy(CNT) and PPy(C6S) membranes. In the course of these studies, membrane characterisation was performed using electrochemical quartz crystal microbalance (EQCM) and scanning electron microscopy (SEM). Permeability of the membranes towards divalent cations was explored using a two compartment transport cell. EQCM results demonstrated that the ion exchange behaviour of polypyrrole is dependent on a number of factors including the type of dopant anion present, the type of ions present in the surrounding medium, the scan rate used during the experiment and the previous history of the polymer film. The morphology of PPy films was found to change when the dopant anion was varied and even when the thickness of the film was altered in some cases. In nearly all cases the permeability of the membranes towards metal ions followed the order K+ > Ca2+ > Mn2+. The one exception was PPy(C6S), for which the permeability followed the order Ca2+ ≥ K+ > Mn2+ > Co2+ > Cr3+. The above permeability sequences show a strong dependence on the size of the metal ions with metal ions having the smallest hydrated radii exhibiting the highest flux. Another factor that affected the permeability towards metal ions was the thickness of the PPy films. Films with the least thickness showed higher metal ion fluxes. Electrochemical control over ion transport across PPy(CNT) membrane was obtained when films composed of the latter were deposited on track-etched Nucleopore® membranes as support material. In contrast, the flux of ions across the same film was concentration gradient dependent when the polymer was deposited on polyvinylidene difluoride membranes as support material. However, electrochemical control over metal ion transport was achieved with a bilayer type of PPy film consisting of PPy(pTS)/PPy(CNT), irrespective of the type of support material. In the course of studying macroscopic charge balance during transport experiments performed using a two compartment transport cell, it was observed that PPy films were non-permselective. A clear correlation between the change in pH in the receiving solution and the ions transported across the membrane was observed. A decrease in solution pH was detected when the polymer membrane acted primarily as an anion exchanger, while an increase in pH occurred when it functioned as a cation exchanger. When there was an approximately equal flux of anions and cations across the polymer membrane, the pH in the receiving solution was in the range 6 - 8. These observations suggest that macroscopic charge balance during the transport of cations and anions across polypyrrole membranes was maintained by introduction of anions (OH-) and cations (H+) produced via electrolysis of water.

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.

α-Pinene oxide and verbenol oxide isomerizations over heterogeneous catalysts

Terpenes are a valuable natural resource for the production of fine chemicals. Turpentine, obtained from biomass and also as a side product of softwood industry, is rich in monoterpenes such as α-pinene and β-pinene, which are widely used as raw materials in the synthesis of flavors, fragrances and pharmaceutical compounds. The rearrangement of their epoxides has been thoroughly studied in recent years, as a method to obtain compounds which are further used in the fine chemical industry. The industrially most desired products of α-pinene oxide isomerization are campholenic aldehyde and trans-carveol. Campholenic aldehyde is an intermediate for the manufacture of sandalwood-like fragrances such as santalol. Trans-carveol is an expensive constituent of the Valencia orange essence oil used in perfume bases and food flavor composition. Furthermore it has been found to exhibit chemoprevention of mammary carcinogenesis. A wide range of iron and ceria supported catalysts were prepared, characterized and tested for α-pinene oxide isomerization in order to selective synthesis of above mentioned products. The highest catalytic activity in the preparation of campholenic aldehyde over iron modified catalysts using toluene as a solvent at 70 °C (total conversion of α-pinene oxide with a selectivity of 66 % to the desired aldehyde) was achieved in the presence of Fe-MCM-41. Furthermore, Fe-MCM-41 catalyst was successfully regenerated without deterioration of catalytic activity and selectivity. The most active catalysts in the synthesis of trans-carveol from α-pinene oxide over iron and ceria modified catalysts in N,N-dimethylacetamide as a solvent at 140 °C (total conversion of α-pinene oxide with selectivity 43 % to trans-carveol) were Fe-Beta-300 and Ce-Si-MCM-41. These catalysts were further tested for an analogous reaction, namely verbenol oxide isomerization. Verbenone is another natural organic compound which can be found in a variety of plants or synthesized by allylic oxidation of α-pinene. An interesting product which is synthesized from verbenone is (1R,2R,6S)-3-methyl-6-(prop-1-en-2-yl)cyclohex-3-ene-1,2-diol. It has been discovered that this diol possesses potent anti-Parkinson activity. The most effective way leading to desired diol starts from verbenone and includes three stages: epoxidation of verbenone to verbenone oxide, reduction of verbenone oxide and subsequent isomerization of obtained verbenol oxide, which is analogous to isomerization of α-pinene oxide. In the research focused on the last step of these synthesis, high selectivity (82 %) to desired diol was achieved in the isomerization of verbenol oxide at a conversion level of 96 % in N,N-dimethylacetamide at 140 °C using iron modified zeolite, Fe-Beta-300. This reaction displayed surprisingly high selectivity, which has not been achieved yet. The possibility of the reuse of heterogeneous catalysts without activity loss was demonstrated.

Neopentyyliglykolin ja propionihapon välisen esteröintireaktion kinetiikka

Glykolien esterit ovat haluttuja pintareaktiivisia aineita. Niitä voidaan valmistaa esteröintireaktiolla karboksyylihappojen kanssa katalyytin läsnä ollessa, jolloin toivottu reaktiotuote on yleensä muodostuva monoesteri. Monoesterin saannon lisäämiseksi reaktiossa muodostuvaa vettä voidaan poistaa jatkuvasti reaktiosta. Reaktion tasapainotilan tutkiminen on kuitenkin tärkeää, jotta reaktion kinetiikka tunnettaisiin mahdollisimman hyvin. Tällöin reaktiotuotteita ei poisteta reaktioseoksesta reaktion aikana. Glykolit esteröityvät happojen kanssa kahdessa vaiheessa. Ensimmäisessä vaiheessa muodostuu monoesteriä ja vettä ja toisessa vaiheessa diesteriä ja vettä. Kokeiden perusteella ensimmäinen vaihe on selvästi toista vaihetta nopeampi reaktio. Kirjallisuudessa on esitetty myös kaksi sivureaktiota, transesteröityminen ja disproportionaatio. Reaktion kinetiikka voidaan kuvata ilman näitä pieniä sivureaktiota, mutta täydellisen kuvaamisen vuoksi on ne myös otettava huomioon. Reaktion kinetiikan tutkimiseksi suoritettiin viisi laboratoriokoetta eri lämpötiloissa neopentyyliglykolilla ja propionihapolla homogeenisen para-tolueenisulfonihapon toimiessa katalyyttina. Lähtöaineiden ja tuotteiden konsentraatioita seurattiin ajan funktiona ja saatujen tulosten perusteella sovitettiin reaktiomekanismin differentiaaliyhtälöiden reaktionopeusvakiot. Nopeusvakioiden lämpötilariippuvuutta tutkittiin Arrheniuksen yhtälön avulla. Lisäksi määritettiin tasapainovakiot kullekin osareaktiolle.