Wet oxidation of recalcitrant lignin waters: Experimental and kinetic studies

The dissertation is based on four articles dealing with recalcitrant lignin water purification. Lignin, a complicated substance and recalcitrant to most treatment technologies, inhibits seriously pulp and paper industry waste management. Therefore, lignin is studied, using WO as a process method for its degradation. A special attention is paid to the improvement in biodegradability and the reduction of lignin content, since they have special importance for any following biological treatment. In most cases wet oxidation is not used as a complete ' mineralization method but as a pre treatment in order to eliminate toxic components and to reduce the high level of organics produced. The combination of wet oxidation with a biological treatment can be a good option due to its effectiveness and its relatively low technology cost. The literature part gives an overview of Advanced Oxidation Processes (AOPs). A hot oxidation process, wet oxidation (WO), is investigated in detail and is the AOP process used in the research. The background and main principles of wet oxidation, its industrial applications, the combination of wet oxidation with other water treatment technologies, principal reactions in WO, and key aspects of modelling and reaction kinetics are presented. There is also given a wood composition and lignin characterization (chemical composition, structure and origin), lignin containing waters, lignin degradation and reuse possibilities, and purification practices for lignin containing waters. The aim of the research was to investigate the effect of the operating conditions of WO, such as temperature, partial pressure of oxygen, pH and initial concentration of wastewater, on the efficiency, and to enhance the process and estimate optimal conditions for WO of recalcitrant lignin waters. Two different waters are studied (a lignin water model solution and debarking water from paper industry) to give as appropriate conditions as possible. Due to the great importance of re using and minimizing the residues of industries, further research is carried out using residual ash of an Estonian power plant as a catalyst in wet oxidation of lignin-containing water. Developing a kinetic model that includes in the prediction such parameters as TOC gives the opportunity to estimate the amount of emerging inorganic substances (degradation rate of waste) and not only the decrease of COD and BOD. The degradation target compound, lignin is included into the model through its COD value (CODligning). Such a kinetic model can be valuable in developing WO treatment processes for lignin containing waters, or other wastewaters containing one or more target compounds. In the first article, wet oxidation of "pure" lignin water was investigated as a model case with the aim of degrading lignin and enhancing water biodegradability. The experiments were performed at various temperatures (110 -190°C), partial oxygen pressures (0.5 -1.5 MPa) and pH (5, 9 and 12). The experiments showed that increasing the temperature notably improved the processes efficiency. 75% lignin reduction was detected at the lowest temperature tested and lignin removal improved to 100% at 190°C. The effect of temperature on the COD removal rate was lower, but clearly detectable. 53% of organics were oxidized at 190°C. The effect of pH occurred mostly on lignin removal. Increasing the pH enhanced the lignin removal efficiency from 60% to nearly 100%. A good biodegradability ratio (over 0.5) was generally achieved. The aim of the second article was to develop a mathematical model for "pure" lignin wet oxidation using lumped characteristics of water (COD, BOD, TOC) and lignin concentration. The model agreed well with the experimental data (R2 = 0.93 at pH 5 and 12) and concentration changes during wet oxidation followed adequately the experimental results. The model also showed correctly the trend of biodegradability (BOD/COD) changes. In the third article, the purpose of the research was to estimate optimal conditions for wet oxidation (WO) of debarking water from the paper industry. The WO experiments were' performed at various temperatures, partial oxygen pressures and pH. The experiments showed that lignin degradation and organics removal are affected remarkably by temperature and pH. 78-97% lignin reduction was detected at different WO conditions. Initial pH 12 caused faster removal of tannins/lignin content; but initial pH 5 was more effective for removal of total organics, represented by COD and TOC. Most of the decrease in organic substances concentrations occurred in the first 60 minutes. The aim of the fourth article was to compare the behaviour of two reaction kinetic models, based on experiments of wet oxidation of industrial debarking water under different conditions. The simpler model took into account only the changes in COD, BOD and TOC; the advanced model was similar to the model used in the second article. Comparing the results of the models, the second model was found to be more suitable for describing the kinetics of wet oxidation of debarking water. The significance of the reactions involved was compared on the basis of the model: for instance, lignin degraded first to other chemically oxidizable compounds rather than directly to biodegradable products. Catalytic wet oxidation of lignin containing waters is briefly presented at the end of the dissertation. Two completely different catalysts were used: a commercial Pt catalyst and waste power plant ash. CWO showed good performance using 1 g/L of residual ash gave lignin removal of 86% and COD removal of 39% at 150°C (a lower temperature and pressure than with WO). It was noted that the ash catalyst caused a remarkable removal rate for lignin degradation already during the pre heating for `zero' time, 58% of lignin was degraded. In general, wet oxidation is not recommended for use as a complete mineralization method, but as a pre treatment phase to eliminate toxic or difficultly biodegradable components and to reduce the high level of organics. Biological treatment is an appropriate post treatment method since easily biodegradable organic matter remains after the WO process. The combination of wet oxidation with subsequent biological treatment can be an effective option for the treatment of lignin containing waters.

Pitkäketjuisten keteenien valmistusreaktiot ja kinetiikka.

Tämän diplomityön tarkoituksena oli tutkia pitkäketjuisten rasvahappokloridien ja alkyyliketeenidimeerien (AKD) valmistusreaktioiden kinetiikkaa. Työn tavoitteena oli saada mittaustuloksia, joiden perusteella voitaisiin kehittää reaktioille kineettinen kineettinen malli ja suorittaa valmistusprosessien alustava optimointi. Teoreettisessa osassa on selvitetty rasvahappokloridien ja alkyyliketeenidimeerien (AKD) eri valmistustapoja. Lisäksi on perehdytty tarkemmin rasvahappokloridien dehydrohalogenointireaktioiden reaktiomekanismeihin ja valmistusprosessin problematiikkaan. Kokeellisessa osassa tutkittiin rasvahappokloridien valmistusta klooraamalla öljyhappoa 30 mol-% fosforitrikloridiylimäärällä lämpötiloissa 45, 50 ja 55OC. Alkyyliketeenidimeerien valmistusta tutkittiin dehydrohalogenoimalla palmitiini- ja öljyhappokloridia 0-30 mol-% trietyyliamiiniylimäärällä inertissä liuottimessa, lämpötiloissa 45 ja 50OC. Reaktiot toteutettiin puolipanosreaktiona, joissa reaktioastiassa olevaan reagenssin ja liuottimen seokseen lisättiin lähtöaine tasaisena massavirtana. Reaktioiden etenemistä seurattiin FT-IR- ja GC-analyysien avulla. Kalorimetrisilla kokeilla tutkittiin öljyhappokloridin dimeroitumisreaktion reaktiolämmön muodostumista ja UV-VIS-analyyseillä seurattiin öljyhappokloridin vanhenemista. Öljyhappokloridin valmistusreaktiolle saatiin hyvä kineettinen malli. Kineettisen mallin puutteena voidaan pitää sitä, ettei kokeissa saatu tietoa mahdollisten sivutuotteiden muodostumisesta. Sovitusohjelmalla saatiin sovitettua estimaatit reaktionopeusvakioille lämpötiloissa 45, 50 ja 55OC. Happokloridien dimeroitumisreaktioiden kineettisen mallin sopivuudelle saatiin suhteellisen hyvä kuva Kuten öljyhappokloridin tapauksessa, myöskään AKD:n valmistusreaktioissa syntynyt sivutuotteita, joten niiden osuutta oletettuun kineettiseen malliin ei tunneta. Sovitusohjelmalla saatiin sovitettua estimaatit reaktionopeusvakioille lämpötiloissa 45 ja 50OC.

Trimmerin kapasiteetin kasvattaminen dynamiikan mallinnuksen avulla

Työn tavoitteena oli kasvattaa sahan dimensiolaitoksella käytettävän trimmerin rakenteellista kapasiteettia. Tavoitteeseen pyrittiin modernisoimalla trimmerin teräyksikköä käyttävää toimilaite ja teräyksikön säätö dynamiikan mallinnuksen avulla. Trimmerin teräyksikön dynamiikka mallinnettiin MATLAB-matematiikkaohjelmistolla kaksiulotteisena kinematiikkamallina ja kolmeulotteisena kinetiikkamallina. Dynamiikkamallien tulosten perusteella valittin teräyksikköä käyttävä toimilaite komponentteineen. Kinetiikkamalliin mallinnettiin trimmeriä käyttävä hydraulipiiri valittuine komponentteineen keskittyneiden paineiden ja puoliempiirisen mallinnuksen periaatteita käyttäen. Teräyksikön työkiertoa säätämään mallinnettiin suljettu takaisinkytketty säätöpiiri. Tuloksien perusteella valittiin optimaalinen toimilaitteen asemointigeometria ja todettiin mallinnetun järjestelmän täyttävän asetetut vaatimukset. Järjestelmää testattiin muuttamalla jarjestelman parametreja ja tutkimalla muutosten vaikutuksia jarjestelman toimintaan. Lisaksi tutkittiin lyhyesti terayksikon rakenteen keventamisen vaikutuksia.

Modeling reaction kinetics and mass transfer in ozonation in water solutions

This dissertation is based on four articles dealing with modeling of ozonation. The literature part of this considers some models for hydrodynamics in bubble column simulation. A literature review of methods for obtaining mass transfer coefficients is presented. The methods presented to obtain mass transfer are general models and can be applied to any gas-liquid system. Ozonation reaction models and methods for obtaining stoichiometric coefficients and reaction rate coefficients for ozonation reactions are discussed in the final section of the literature part. In the first article, ozone gas-liquid mass transfer into water in a bubble column was investigated for different pH values. A more general method for estimation of mass transfer and Henry’s coefficient was developed from the Beltrán method. The ozone volumetric mass transfer coefficient and the Henry’s coefficient were determined simultaneously by parameter estimation using a nonlinear optimization method. A minor dependence of the Henry’s law constant on pH was detected at the pH range 4 - 9. In the second article, a new method using the axial dispersion model for estimation of ozone self-decomposition kinetics in a semi-batch bubble column reactor was developed. The reaction rate coefficients for literature equations of ozone decomposition and the gas phase dispersion coefficient were estimated and compared with the literature data. The reaction order in the pH range 7-10 with respect to ozone 1.12 and 0.51 the hydroxyl ion were obtained, which is in good agreement with literature. The model parameters were determined by parameter estimation using a nonlinear optimization method. Sensitivity analysis was conducted using object function method to obtain information about the reliability and identifiability of the estimated parameters. In the third article, the reaction rate coefficients and the stoichiometric coefficients in the reaction of ozone with the model component p-nitrophenol were estimated at low pH of water using nonlinear optimization. A novel method for estimation of multireaction model parameters in ozonation was developed. In this method the concentration of unknown intermediate compounds is presented as a residual COD (chemical oxygen demand) calculated from the measured COD and the theoretical COD for the known species. The decomposition rate of p-nitrophenol on the pathway producing hydroquinone was found to be about two times faster than the p-nitrophenol decomposition rate on the pathway producing 4- nitrocatechol. In the fourth article, the reaction kinetics of p-nitrophenol ozonation was studied in a bubble column at pH 2. Using the new reaction kinetic model presented in the previous article, the reaction kinetic parameters, rate coefficients, and stoichiometric coefficients as well as the mass transfer coefficient were estimated with nonlinear estimation. The decomposition rate of pnitrophenol was found to be equal both on the pathway producing hydroquinone and on the path way producing 4-nitrocathecol. Comparison of the rate coefficients with the case at initial pH 5 indicates that the p-nitrophenol degradation producing 4- nitrocathecol is more selective towards molecular ozone than the reaction producing hydroquinone. The identifiability and reliability of the estimated parameters were analyzed with the Marcov chain Monte Carlo (MCMC) method. @All rights reserved. No part of the publication may be reproduced, stored in a retrieval system, or transmitted, in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise, without the prior permission of the author.

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.

Acidic dissolution of iron oxides and regeneration of a ceramic filter medium

The dewatering of iron ore concentrates requires large capacity in addition to producing a cake with low moisture content. Such large processes are commonly energy intensive and means to lower the specific energy consumption are needed. Ceramic capillary action disc filters incorporate a novel filter medium enabling the harnessing of capillary action, which results in decreased energy consumption in comparison to traditional filtration technologies. As another benefit, the filter medium is mechanically and chemically more durable than, for example, filter cloths and can, thus, withstand harsh operating conditions and possible regeneration better than other types of filter media. In iron ore dewatering, the regeneration of the filter medium is done through a combination of several techniques: (1) backwashing, (2) ultrasonic cleaning, and (3) acid regeneration. Although it is commonly acknowledged that the filter medium is affected by slurry particles and extraneous compounds, published research, especially in the field of dewatering of mineral concentrates, is scarce. Whereas the regenerative effect of backwashing and ultrasound are more or less mechanical, regeneration with acids is based on chemistry. The chemistry behind the acid regeneration is, naturally, dissolution. The dissolution of iron oxide particles has been extensively studied over several decades but those studies may not necessarily be directly applicable in the regeneration of the filter medium which has undergone interactions with the slurry components. The aim of this thesis was to investigate if free particle dissolution indeed correlates with the regeneration of the filter medium. For this purpose, both free particle dissolution and dissolution of surface adhered particles were studied. The focus was on acidic dissolution of iron oxide particles and on the study of the ceramic filter medium used in the dewatering of iron ore concentrates. The free particle dissolution experiments show that the solubility of synthetic fine grained iron oxide particles in oxalic acid could be explained through linear models accounting for the effects of temperature and acid concentration, whereas the dissolution of a natural magnetite is not so easily explained by such models. In addition, the kinetic experiments performed both support and contradict the work of previous authors: the suitable kinetic model here supports previous research suggesting solid state reduction to be the reaction mechanism of hematite dissolution but the formation of a stable iron oxalate is not supported by the results of this research. Several other dissolution mechanisms have also been suggested for iron oxide dissolution in oxalic acid, indicating that the details of oxalate promoted reductive dissolution are not yet agreed and, in this respect, this research offers added value to the community. The results of the regeneration experiments with the ceramic filter media show that oxalic acid is highly effective in removing iron oxide particles from the surface of the filter medium. The dissolution of those particles did not, however, exhibit the expected behaviour, i.e. complete dissolution. The results of this thesis show that although the regeneration of the ceramic filter medium with acids incorporates the dissolution of slurry particles from the surface of the filter medium, the regeneration cannot be assessed purely based upon free particle dissolution. A steady state, dependent on temperature and on the acid concentration, was observed in the dissolution of particles from the surface even though the limit of solubility of free iron oxide particles had not been reached. Both the regeneration capacity and efficiency, with regards to the removal of iron oxide particles, was found to be temperature dependent, but was not affected by the acid concentration. This observation further suggests that the removal of the surface adhered particles does not follow the dissolution of free particles, which do exhibit a dependency on the acid concentration. In addition, changes in the permeability and in the pore structure of the filter medium were still observed after the bulk concentration of dissolved iron had reached a steady state. Consequently, the regeneration of the filter medium continued after the dissolution of particles from the surface had ceased. This observation suggests that internal changes take place at the final stages of regeneration. The regeneration process could, in theory, be divided into two, possibly overlapping, stages: (1) dissolution of surface-adhered particles, and (2) dissolution of extraneous compounds from within the pore structure. In addition to the fundamental knowledge generated during this thesis, tools to assess the effects of parameters on the regeneration of the ceramic filter medium are needed. It has become clear that the same tools used to estimate the dissolution of free particles cannot be used to estimate the regeneration of a filter medium unless only a robust characterisation of the order of regeneration efficiency is needed.

Extraction of polymeric galactoglucomannans from spruce wood by pressurized hot water

The major type of non-cellulosic polysaccharides (hemicelluloses) in softwoods, the partly acetylated galactoglucomannans (GGMs), which comprise about 15% of spruce wood, have attracted growing interest because of their potential to become high-value products with applications in many areas. The main objective of this work was to explore the possibilities to extract galactoglucomannans in native, polymeric form in high yield from spruce wood with pressurised hot-water, and to obtain a deeper understanding of the process chemistry involved. Spruce (Picea abies) chips and ground wood particles were extracted using an accelerated solvent extractor (ASE) in the temperature range 160 – 180°C. Detailed chemical analyses were done on both the water extracts and the wood residues. As much as 80 – 90% of the GGMs in spruce wood, i.e. about 13% based on the original wood, could be extracted from ground spruce wood with pure water at 170 – 180°C with an extraction time of 60 min. GGMs comprised about 75% of the extracted carbohydrates and about 60% of the total dissolved solids. Other substances in the water extracts were xylans, arabinogalactans, pectins, lignin and acetic acid. The yields from chips were only about 60% of that from ground wood. Both the GGMs and other non-cellulosic polysaccharides were extensively hydrolysed at severe extraction conditions when pH dropped to the level of 3.5. Addition of sodium bicarbonate increased the yields of polymeric GGMs at low additions, 2.5 – 5 mM, where the end pH remained around 3.9. However, at higher addition levels the yields decreased, mainly because the acetyl groups in GGMs were split off, leading to a low solubility of GGMs. Extraction with buffered water in the pH range 3.8 – 4.4 gave similar yields as with plain water, but gave a higher yield of polymeric GGMs. Moreover, at these pH levels the hydrolysis of acetyl groups in GGMs was significantly inhibited. It was concluded that hot-water extraction of polymeric GGMs in good yields (up to 8% of wood) demands appropriate control of pH, in a narrow range about 4. These results were supported by a study of hydrolysis of GGM at constant pH in the range of 3.8 – 4.2 where a kinetic model for degradation of GGM was developed. The influence of wood particle size on hot-water extraction was studied with particles in the range of 0.1 – 2 mm. The smallest particles (< 0.1 mm) gave 20 – 40% higher total yield than the coarsest particles (1.25 – 2 mm). The difference was greatest at short extraction times. The results indicated that extraction of GGMs and other polysaccharides is limited mainly by the mass transfer in the fibre wall, and for coarse wood particles also in the wood matrix. Spruce sapwood, heartwood and thermomechnical pulp were also compared, but only small differences in yields and composition of extracts were found. Two methods for isolation and purification of polymeric GGMs, i.e. membrane filtration and precipitation in ethanol-water, were compared. Filtration through a series of membranes with different pore sizes separated GGMs of different molar masses, from polymers to oligomers. Polysaccharides with molar mass higher than 4 kDa were precipitated in ethanol-water. GGMs comprised about 80% of the precipitated polysaccharides. Other polysaccharides were mainly arabinoglucuronoxylans and pectins. The ethanol-precipitated GGMs were by 13C NMR spectroscopy verified to be very similar to GGMs extracted from spruce wood in low yield at a much lower temperature, 90°C. The obtained large body of experimental data could be utilised for further kinetic and economic calculations to optimise technical hot-water extractionof softwoods.

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

Silicon-containing species in used lube oil re-refining

Nowadays, the re-refining of the used lube oils has gained worldwide a lot of attention due to the necessity for added environmental protection and increasingly stringent environmental legislation. One of the parameters determining the quality of the produced base oils is the composition of feedstock. Estimation of the chemical composition of the used oil collected from several European locations showed that the hydrocarbon structure of the motor oil is changed insignificantly during its operation and the major part of the changes is accounted for with depleted oil additives. In the lube oil re-refining industry silicon, coming mainly from antifoaming agents, is recognized to be a contaminant generating undesired solid deposits in various locations in the re-refining units. In this thesis, a particular attention was paid to the mechanism of solid product formation during the alkali treatment process of silicon-containing used lube oils. The transformations of a model siloxane, tetramethyldisiloxane (TMDS), were studied in a batch reactor at industrially relevant alkali treatment conditions (low temperature, short reaction time) using different alkali agents. The reaction mechanism involving solid alkali metal silanolates was proposed. The experimental data obtained demonstrated that the solids were dominant products at low temperature and short reaction time. The liquid products in the low temperature reactions were represented mainly by linear siloxanes. The prolongation of reaction time resulted in reduction of solids, whereas both temperature and time increase led to dominance of cyclic products in the reaction mixture. Experiments with the varied reaction time demonstrated that the concentration of cyclic trimer being the dominant in the beginning of the reaction diminished with time, whereas the cyclic tetramer tended to increase. Experiments with lower sodium hydroxide concentration showed the same effect. In addition, a decrease of alkali agent concentration in the initial reaction mixture accelerated TMDS transformation reactions resulting in solely liquid cyclic siloxanes yields. Comparison of sodium and potassium hydroxides applied as an alkali agent demonstrated that potassium hydroxide was more efficient, since the activation energy in KOH presence was almost 2-fold lower than that for sodium hydroxide containing reaction mixture. Application of potassium hydroxide for TMDS transformation at 100° C with 3 hours reaction time resulted in 20 % decrease of solid yields compared to NaOH-containing mixture. Moreover, TMDS transformations in the presence of sodium silanolate applied as an alkali agent led to formation of only liquid products without formation of the undesired solids. On the basis of experimental data and the proposed reaction mechanism, a kinetic model was developed, which provided a satisfactory description of the experimental results. Suitability of the selected siloxane as a relevant model of industrial silicon-containing compounds was verified by investigation of the commercially available antifoam agent in base-catalyzed conditions.

The role of metal ions in selective and sustainable processes

Sustainability and recycling are core values in today’s industrial operations. New materials, products and processes need to be designed in such a way as to consume fewer of the diminishing resources we have available and to put as little strain on the environment as possible. An integral part of this is cleaning and recycling. New processes are to be designed to improve the efficiency in this aspect. Wastewater, including municipal wastewaters, is treated in several steps including chemical and mechanical cleaning of waters. Well-cleaned water can be recycled and reused. Clean water for everyone is one of the greatest challenges we are facing today. Ferric sulphate, made by oxidation from ferrous sulphate, is used in water purification. The oxidation of ferrous sulphate, FeSO4, to ferric sulphate in acidic aqueous solutions of H2SO4 over finely dispersed active carbon particles was studied in a vigorously stirred batch reactor. Molecular oxygen was used as the oxidation agent and several catalysts were screened: active carbon, active carbon impregnated with Pt, Rh, Pd and Ru. Both active carbon and noble metal-active carbon catalysts enhanced the oxidation rate considerably. The order of the noble metals according to the effect was: Pt >> Rh > Pd, Ru. By the use of catalysts, the production capacities of existing oxidation units can be considerably increased. Good coagulants have a high charge on a long polymer chain effectively capturing dirty particles of the opposite charge. Analysis of the reaction product indicated that it is possible to obtain polymeric iron-based products with good coagulation properties. Systematic kinetic experiments were carried out at the temperature and pressure ranges of 60B100°C and 4B10 bar, respectively. The results revealed that both non-catalytic and catalytic oxidation of Fe2+ to Fe3+ take place simultaneously. The experimental data were fitted to rate equations, which were based on a plausible reaction mechanism: adsorption of dissolved oxygen on active carbon, electron transfer from Fe2+ ions to adsorbed oxygen and formation of surface hydroxyls. A comparison of the Fe2+ concentrations predicted by the kinetic model with the experimentally observed concentrations indicated that the mechanistic rate equations were able to describe the intrinsic oxidation kinetics of Fe2+ over active carbon and active carbon-noble metal catalysts. Engineering aspects were closely considered and effort was directed to utilizing existing equipment in the production of the new coagulant. Ferrous sulphate can be catalytically oxidized to produce a novel long-chained polymeric iron-based flocculent in an easy and affordable way in existing facilities. The results can be used for modelling the reactors and for scale-up. Ferric iron (Fe3+) was successfully applied for the dissolution of sphalerite. Sphalerite contains indium, gallium and germanium, among others, and the application can promote their recovery. The understanding of the reduction process of ferric to ferrous iron can be used to develop further the understanding of the dissolution mechanisms and oxidation of ferrous sulphate. Indium, gallium and germanium face an ever-increasing demand in the electronics industry, among others. The supply is, however, very limited. The fact that most part of the material is obtained through secondary production means that real production quota depends on the primary material production. This also sets the pricing. The primary production material is in most cases zinc and aluminium. Recycling of scrap material and the utilization of industrial waste, containing indium, gallium and geranium, is a necessity without real options. As a part of this study plausible methods for the recovery of indium, gallium and germanium have been studied. The results were encouraging and provided information about the precipitation of these valuables from highly acidic solutions. Indium and gallium were separated from acidic sulphuric acid solutions by precipitation with basic sulphates such as alunite or they were precipitated as basic sulphates of their own as galliunite and indiunite. Germanium may precipitate as a basic sulphate of a mixed composition. The precipitation is rapid and the selectivity is good. When the solutions contain both indium and gallium then the results show that gallium should be separated before indium to achieve a better selectivity. Germanium was separated from highly acidic sulphuric acid solutions containing other metals as well by precipitating with tannic acid. This is a highly selective method. According to the study other commonly found metals in the solution do not affect germanium precipitation. The reduction of ferric iron to ferrous, the precipitation of indium, gallium and germanium, and the dissolution of the raw materials are strongly depending on temperature and pH. The temperature and pH effect were studied and which contributed to the understanding and design of the different process steps. Increased temperature and reduced pH improve the reduction rate. Finally, the gained understanding in the studied areas can be employed to develop better industrial processes not only on a large scale but also increasingly on a smaller scale. The small amounts of indium, gallium and germanium may favour smaller and more locally bound recovery.

A reaction engineering approach to homogeneously catalyzed glycerol hydrochlorination

The production of biodiesel through transesterification has created a surplus of glycerol on the international market. In few years, glycerol has become an inexpensive and abundant raw material, subject to numerous plausible valorisation strategies. Glycerol hydrochlorination stands out as an economically attractive alternative to the production of biobased epichlorohydrin, an important raw material for the manufacturing of epoxy resins and plasticizers. Glycerol hydrochlorination using gaseous hydrogen chloride (HCl) was studied from a reaction engineering viewpoint. Firstly, a more general and rigorous kinetic model was derived based on a consistent reaction mechanism proposed in the literature. The model was validated with experimental data reported in the literature as well as with new data of our own. Semi-batch experiments were conducted in which the influence of the stirring speed, HCl partial pressure, catalyst concentration and temperature were thoroughly analysed and discussed. Acetic acid was used as a homogeneous catalyst for the experiments. For the first time, it was demonstrated that the liquid-phase volume undergoes a significant increase due to the accumulation of HCl in the liquid phase. Novel and relevant features concerning hydrochlorination kinetics, HCl solubility and mass transfer were investigated. An extended reaction mechanism was proposed and a new kinetic model was derived. The model was tested with the experimental data by means of regression analysis, in which kinetic and mass transfer parameters were successfully estimated. A dimensionless number, called Catalyst Modulus, was proposed as a tool for corroborating the kinetic model. Reactive flash distillation experiments were conducted to check the commonly accepted hypothesis that removal of water should enhance the glycerol hydrochlorination kinetics. The performance of the reactive flash distillation experiments were compared to the semi-batch data previously obtained. An unforeseen effect was observed once the water was let to be stripped out from the liquid phase, exposing a strong correlation between the HCl liquid uptake and the presence of water in the system. Water has revealed to play an important role also in the HCl dissociation: as water was removed, the dissociation of HCl was diminished, which had a retarding effect on the reaction kinetics. In order to obtain a further insight on the influence of water on the hydrochlorination reaction, extra semi-batch experiments were conducted in which initial amounts of water and the desired product were added. This study revealed the possibility to use the desired product as an ideal “solvent” for the glycerol hydrochlorination process. A co-current bubble column was used to investigate the glycerol hydrochlorination process under continuous operation. The influence of liquid flow rate, gas flow rate, temperature and catalyst concentration on the glycerol conversion and product distribution was studied. The fluid dynamics of the system showed a remarkable behaviour, which was carefully investigated and described. Highspeed camera images and residence time distribution experiments were conducted to collect relevant information about the flow conditions inside the tube. A model based on the axial dispersion concept was proposed and confronted with the experimental data. The kinetic and solubility parameters estimated from the semi-batch experiments were successfully used in the description of mass transfer and fluid dynamics of the bubble column reactor. In light of the results brought by the present work, the glycerol hydrochlorination reaction mechanism has been finally clarified. It has been demonstrated that the reactive distillation technology may cause drawbacks to the glycerol hydrochlorination reaction rate under certain conditions. Furthermore, continuous reactor technology showed a high selectivity towards monochlorohydrins, whilst semibatch technology was demonstrated to be more efficient towards the production of dichlorohydrins. Based on the novel and revealing discoveries brought by the present work, many insightful suggestions are made towards the improvement of the production of αγ-dichlorohydrin on an industrial scale.

Wet Oxidation of Concentrated Wastewaters: Process Combination and Reaction Kinetic Modeling

The accumulation of aqueous pollutants is becoming a global problem. The search for suitable methods and/or combinations of water treatment processes is a task that can slow down and stop the process of water pollution. In this work, the method of wet oxidation was considered as an appropriate technique for the elimination of the impurities present in paper mill process waters. It has been shown that, when combined with traditional wastewater treatment processes, wet oxidation offers many advantages. The combination of coagulation and wet oxidation offers a new opportunity for the improvement of the quality of wastewater designated for discharge or recycling. First of all, the utilization of coagulated sludge via wet oxidation provides a conditioning process for the sludge, i.e. dewatering, which is rather difficult to carry out with untreated waste. Secondly, Fe2(SO4)3, which is employed earlier as a coagulant, transforms the conventional wet oxidation process into a catalytic one. The use of coagulation as the post-treatment for wet oxidation can offer the possibility of the brown hue that usually accompanies the partial oxidation to be reduced. As a result, the supernatant is less colored and also contains a rather low amount of Fe ions to beconsidered for recycling inside mills. The thickened part that consists of metal ions is then recycled back to the wet oxidation system. It was also observed that wet oxidation is favorable for the degradation of pitch substances (LWEs) and lignin that are present in the process waters of paper mills. Rather low operating temperatures are needed for wet oxidation in order to destruct LWEs. The oxidation in the alkaline media provides not only the faster elimination of pitch and lignin but also significantly improves the biodegradable characteristics of wastewater that contains lignin and pitch substances. During the course of the kinetic studies, a model, which can predict the enhancements of the biodegradability of wastewater, was elaborated. The model includes lumped concentrations suchas the chemical oxygen demand and biochemical oxygen demand and reflects a generalized reaction network of oxidative transformations. Later developments incorporated a new lump, the immediately available biochemical oxygen demand, which increased the fidelity of the predictions made by the model. Since changes in biodegradability occur simultaneously with the destruction of LWEs, an attempt was made to combine these two facts for modeling purposes.

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.

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.