Methane combustion kinetic rate constants determination: an ill-posed inverse problem analysis

Methane combustion was studied by the Westbrook and Dryer model. This well-established simplified mechanism is very useful in combustion science, for computational effort can be notably reduced. In the inversion procedure to be studied, rate constants are obtained from [CO] concentration data. However, when inherent experimental errors in chemical concentrations are considered, an ill-conditioned inverse problem must be solved for which appropriate mathematical algorithms are needed. A recurrent neural network was chosen due to its numerical stability and robustness. The proposed methodology was compared against Simplex and Levenberg-Marquardt, the most used methods for optimization problems.

Kinetic modeling of the biosorption of Cd2+ ions from aqueous solutions onto Eichhornia crassipes roots using potentiometry: low-cost alternative to conventional methods

This work presents the use of potentiometric measurements for kinetic studies of biosorption of Cd2+ ions from aqueous solutions on Eichhornia crassipes roots. The open circuit potential of the Cd/Cd2+ electrode of the first kind was measured during the bioadsorption process. The amount of Cd2+ ions accumulated was determined in real time. The data were fit to different models, with the pseudo-second-order model proving to be the best in describing the data. The advantages and limitations of the methodology proposed relative to the traditional method are discussed.

A new approach to evaluate immobilization of β-galactosidase on Eupergit® C: structural, kinetic, and thermal characterization

This study aimed to evaluate β-galactosidase immobilization. For this purpose, the ionic strength of the buffer, reaction time, amount of the immobilization support, and pH were evaluated by a central composite design. Assay 8, which consisted of 1.5 mol L-1 phosphate buffer (pH 7.5) and a reaction time of 2 h, produced the maximum yield. Eupergit® C (400 mg) was subsequently used as an immobilization support. Immobilization kinetics wereinvestigated, and a significant increase in the yield was obtained after immobilization compared with that obtained from assay 8 (22.0 U mL-1 vs. 15.6 U mL-1). The enzyme efficiency of actuation was evaluated using o-nitrophenyl-β-D-galactopyranoside and lactose, with lactose providing better results. The reuse of β-galactosidase was evaluated, and more than 50% of the initial enzyme activity was maintained after five cycles of use. Enzyme characterization revealed that immobilization improved some aspects of the thermostability of β-galactosidase.

KINETIC STUDY OF SELECTIVE GAS-PHASE OXIDATION OF ISOPROPANOL TO ACETONE USING MONOCLINIC ZRO2 AS A CATALYST

Zirconia was prepared by a precipitation method and calcined at 723 K, 1023 K, and 1253 K in order to obtain monoclinic zirconia. The prepared zirconia was characterized by XRD, SEM, EDX, surface area and pore size analyzer, and particle size analyzer. Monoclinic ZrO2 as a catalyst was used for the gas-phase oxidation of isopropanol to acetone in a Pyrex-glass-flow-type reactor with a temperature range of 443 K - 473 K. It was found that monoclinic ZrO2 shows remarkable catalytic activity (68%) and selectivity (100%) for the oxidation of isopropanol to acetone. This kinetic study reveals that the oxidation of isopropanol to acetone follows the L-H mechanism.

Thermal modelling of semiconductor switches in three-level Neutral Point Clamped Medium Voltage DTC-inverter

The aim of this thesis is to investigate the thermal loading of medium voltage three-level NPC inverter’s semiconductor IGCT switches in different operation points. The objective is to reach both a fairly accurate off-line simulation program and also so simple a simulation model that its implementation into an embedded system could be reasonable in practice and a real time use should become feasible. Active loading limitation of the inverter can be realized with a thermal model which is practical in a real time use. Determining of the component heating has been divided into two parts; defining of component losses and establishing the structure of a thermal network. Basics of both parts are clarified. The simulation environment is Matlab-Simulink. Two different models are constructed – a more accurate one and a simplified one. Potential simplifications are clarified with the help of the first one. Simplifications are included in the latter model and the functionalities of both models are compared. When increasing the calculation time step a decreased number of considered components and time constants of the thermal network can be used in the simplified model. Heating of a switching component is dependent on its topological position and inverter’s operation point. The output frequency of the converter defines mainly which one of the switching components is – because of its losses and heating – the performance limiting component of the converter. Comparison of results given by different thermal models demonstrates that with larger time steps, describing of fast occurring switching losses becomes difficult. Generally articles and papers dealing with this subject are written for two-level inverters. Also inverters which apply direct torque control (DTC) are investigated rarely from the heating point of view. Hence, this thesis completes the former material.

Demographic Modelling of Human Population Growth

This work presents models and methods that have been used in producing forecasts of population growth. The work is intended to emphasize the reliability bounds of the model forecasts. Leslie model and various versions of logistic population models are presented. References to literature and several studies are given. A lot of relevant methodology has been developed in biological sciences. The Leslie modelling approach involves the use of current trends in mortality,fertility, migration and emigration. The model treats population divided in age groups and the model is given as a recursive system. Other group of models is based on straightforward extrapolation of census data. Trajectories of simple exponential growth function and logistic models are used to produce the forecast. The work presents the basics of Leslie type modelling and the logistic models, including multi- parameter logistic functions. The latter model is also analysed from model reliability point of view. Bayesian approach and MCMC method are used to create error bounds of the model predictions.

Dehydration and volatilization non isothermic kinetic of the solid state aluminium 8-hydroxyquinolinate

Al(C9H6ON)3.2.5H2O was precipitated from the mixture of an aqueous solution of aluminium ion and an acid solution of 8-hydroxyquinoline, by increasing the pH value to 9.5 with ammonia aqueous solution. The TG curves in nitrogen atmosphere present mass losses due to dehydration, partial volatilisation (sublimation plus vaporisation) of the anhydrous compound followed by thermal decomposition with the formation of a mixture of carbonaceous and residues. The relation between sublimation and vaporisation depends on the heating rate used. The non isothermic integral isoconventional methods as linear equations of Ozawa-Flynn-Wall and Kissinger-Akahira-Sunose (KAS) were used to obtain the kinetic parameters from TG and DTA curves, respectively. Despite the fact that both dehydration and volatilisation reactions follow the linearity by using both methods, only for the volatilisation reaction the validity condition, 20<= E/RT<= 50, was verified.

Kinetic evaluation of the dehydration of Yb(III) Lu(III) and Y(III) 4-chlorobenzylidenepyruvate by thermogravimetry (TG)

The dehydration kinetic of Yb, Lu and Y 4-chlorobenzylidenepyruvate was studied by using thermogravimetry and the kinetics parameters obtained by Flynn and Wall method suggest that the dehydration step follows a first order mechanism. The activation energies calculated were 103.6, 96.6 and 97.2 kJ/mol and the lifetime considering the temperature of 31 and 101 º C for the dehydration of these compounds were 23, 26, 31 minutes and 0.6, 1.3 and 1.4 seconds, respectively. The results have similar values and suggest that the water is attached in the same way.

Numerical modelling of small supersonic axial flow turbines

Supersonic axial turbine stages typically exhibit lower efficiencies than subsonic axial turbine stages. One reason for the lower efficiency is the occurrence of shock waves. With higher pressure ratios the flow inside the turbine becomes relatively easily supersonic if there is only one turbine stage. Supersonic axial turbines can be designed in smaller physical size compared to subsonic axial turbines of same power. This makes them good candidates for turbochargers in large diesel engines, where space can be a limiting factor. Also the production costs are lower for a supersonic axial turbine stage than for two subsonic stages. Since supersonic axial turbines are typically low reaction turbines, they also create lower axial forces to be compensated with bearings compared to high reaction turbines. The effect of changing the stator-rotor axial gap in a small high (rotational) speed supersonic axial flow turbine is studied in design and off-design conditions. Also the effect of using pulsatile mass flow at the supersonic stator inlet is studied. Five axial gaps (axial space between stator and rotor) are modeled using threedimensional computational fluid dynamics at the design and three axial gaps at the off-design conditions. Numerical reliability is studied in three independent studies. An additional measurement is made with the design turbine geometry at intermediate off-design conditions and is used to increase the reliability of the modelling. All numerical modelling is made with the Navier-Stokes solver Finflo employing Chien’s k ¡ ² turbulence model. The modelling of the turbine at the design and off-design conditions shows that the total-to-static efficiency of the turbine decreases when the axial gap is increased in both design and off-design conditions. The efficiency drops almost linearily at the off-design conditions, whereas the efficiency drop accelerates with increasing axial gap at the design conditions. The modelling of the turbine stator with pulsatile inlet flow reveals that the mass flow pulsation amplitude is decreased at the stator throat. The stator efficiency and pressure ratio have sinusoidal shapes as a function of time. A hysteresis-like behaviour is detected for stator efficiency and pressure ratio as a function of inlet mass flow, over one pulse period. This behaviour arises from the pulsatile inlet flow. It is important to have the smallest possible axial gap in the studied turbine type in order to maximize the efficiency. The results for the whole turbine can also be applied to some extent in similar turbines operating for example in space rocket engines. The use of a supersonic stator in a pulsatile inlet flow is shown to be possible.

Spectroscopic analysis of C80 doping with group III and group V elements using semiempirical PM3 molecular modelling technique

Semiempirical calculations at the level of PM3 of theory were carried out to study the structural and electronic properties of C80 and some of its doped derivatives with the elements of group III and V at the level of PM3 of theory. We have selected these elements to be substituted in the fullerene-C80 cage in order to show the effect of such structural change on the electronic properties of the molecules studied. The theoretical IR spectra, some of physical and chemical properties of the molecules studied are obtained and discussed.

A novel indicator reaction for the catalytic determination of V(V) at ppb levels by the kinetic spectrophotometric method

A novel sensitive and relatively selective kinetic method is presented for the determination of V(V), based on its catalytic effect on the oxidation reaction of Ponceau Xylydine by potassium bromate in presence of 5-sulfosalicylic acid (SSA) as activator. The reaction was monitored spectrophotometrically by measuring the decrease in absorbance of Ponceau Xylydine at 640 nm between 0.5 to 7 min (the fixed time method) in H3PO4 medium at 25ºC. The effect of various parameters such as concentrations of H3PO4, SSA, bromate and Ponceau Xylydine, temperature and ionic strength on the rate of net reaction were studied. The method is free from most interferences, especially from large amounts of V(IV). The decrease in absorbance is proportional to the concentration of V(V) over the entire concentration range tested (1-15 ng mL−1) with a detection limit of 0.46 ng mL-1 (according to statistical 3Sblank/k criterion) and a coefficient of variation (CV) of 1.8% (for ten replicate measurement at 95% confidence level). The proposed method suffers few interferences such as Cr(VI) and Hg(II) ions. The method was successfully applied to the determination of V(V) in tap water, drinking water, bottled mineral water samples and a certified standard reference material such as SRM-1640 with satisfactory results. The vanadium contents of water samples were also determined by FAAS for a comparison. The recovery of spiked vanadium(V) was found to be quantitative and the reproducibility was satisfactory. It was observed that the results of the SRM 1640 were in good agreement with the certified value.

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.

Effect of hydrodynamics on modelling, monitoring and control of crystallization

Crystallization is a purification method used to obtain crystalline product of a certain crystal size. It is one of the oldest industrial unit processes and commonly used in modern industry due to its good purification capability from rather impure solutions with reasonably low energy consumption. However, the process is extremely challenging to model and control because it involves inhomogeneous mixing and many simultaneous phenomena such as nucleation, crystal growth and agglomeration. All these phenomena are dependent on supersaturation, i.e. the difference between actual liquid phase concentration and solubility. Homogeneous mass and heat transfer in the crystallizer would greatly simplify modelling and control of crystallization processes, such conditions are, however, not the reality, especially in industrial scale processes. Consequently, the hydrodynamics of crystallizers, i.e. the combination of mixing, feed and product removal flows, and recycling of the suspension, needs to be thoroughly investigated. Understanding of hydrodynamics is important in crystallization, especially inlargerscale equipment where uniform flow conditions are difficult to attain. It is also important to understand different size scales of mixing; micro-, meso- and macromixing. Fast processes, like nucleation and chemical reactions, are typically highly dependent on micro- and mesomixing but macromixing, which equalizes the concentrations of all the species within the entire crystallizer, cannot be disregarded. This study investigates the influence of hydrodynamics on crystallization processes. Modelling of crystallizers with the mixed suspension mixed product removal (MSMPR) theory (ideal mixing), computational fluid dynamics (CFD), and a compartmental multiblock model is compared. The importance of proper verification of CFD and multiblock models is demonstrated. In addition, the influence of different hydrodynamic conditions on reactive crystallization process control is studied. Finally, the effect of extreme local supersaturation is studied using power ultrasound to initiate nucleation. The present work shows that mixing and chemical feeding conditions clearly affect induction time and cluster formation, nucleation, growth kinetics, and agglomeration. Consequently, the properties of crystalline end products, e.g. crystal size and crystal habit, can be influenced by management of mixing and feeding conditions. Impurities may have varying impacts on crystallization processes. As an example, manganese ions were shown to replace magnesium ions in the crystal lattice of magnesium sulphate heptahydrate, increasing the crystal growth rate significantly, whereas sodium ions showed no interaction at all. Modelling of continuous crystallization based on MSMPR theory showed that the model is feasible in a small laboratoryscale crystallizer, whereas in larger pilot- and industrial-scale crystallizers hydrodynamic effects should be taken into account. For that reason, CFD and multiblock modelling are shown to be effective tools for modelling crystallization with inhomogeneous mixing. The present work shows also that selection of the measurement point, or points in the case of multiprobe systems, is crucial when process analytical technology (PAT) is used to control larger scale crystallization. The thesis concludes by describing how control of local supersaturation by highly localized ultrasound was successfully applied to induce nucleation and to control polymorphism in reactive crystallization of L-glutamic acid.