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.

Towards desired crystalline product properties: In-situ monitoring of batch crystallization

The objective of industrial crystallization is to obtain a crystalline product which has the desired crystal size distribution, mean crystal size, crystal shape, purity, polymorphic and pseudopolymorphic form. Effective control of the product quality requires an understanding of the thermodynamics of the crystallizing system and the effects of operation parameters on the crystalline product properties. Therefore, obtaining reliable in-line information about crystal properties and supersaturation, which is the driving force of crystallization, would be very advantageous. Advanced techniques, such asRaman spectroscopy, attenuated total reflection Fourier transform infrared (ATR FTIR) spectroscopy, and in-line imaging techniques, offer great potential for obtaining reliable information during crystallization, and thus giving a better understanding of the fundamental mechanisms (nucleation and crystal growth) involved. In the present work, the relative stability of anhydrate and dihydrate carbamazepine in mixed solvents containing water and ethanol were investigated. The kinetics of the solvent mediated phase transformation of the anhydrate to hydrate in the mixed solvents was studied using an in-line Raman immersion probe. The effects of the operation parameters in terms of solvent composition, temperature and the use of certain additives on the phase transformation kineticswere explored. Comparison of the off-line measured solute concentration and the solid-phase composition measured by in-line Raman spectroscopy allowedthe identification of the fundamental processes during the phase transformation. The effects of thermodynamic and kinetic factors on the anhydrate/hydrate phase of carbamazepine crystals during cooling crystallization were also investigated. The effect of certain additives on the batch cooling crystallization of potassium dihydrogen phosphate (KDP) wasinvestigated. The crystal growth rate of a certain crystal face was determined from images taken with an in-line video microscope. An in-line image processing method was developed to characterize the size and shape of thecrystals. An ATR FTIR and a laser reflection particle size analyzer were used to study the effects of cooling modes and seeding parameters onthe final crystal size distribution of an organic compound C15. Based on the obtained results, an operation condition was proposed which gives improved product property in terms of increased mean crystal size and narrowersize distribution.

Experimental insights into deformation dynamics and intermittency in rapid granular shear flows

This work concerns the experimental study of rapid granular shear flows in annular Couette geometry. The flow is induced by continuous driving of the horizontal plate at the top of the granular bed in an annulus. The compressive pressure, driving torque, instantaneous bed height and rotational speed of the shearing plate are measured. Moreover, local stress fluctuations are measured in a medium made of steel spheres 2 and 3 mm in diameter. Both monodisperse packing and bidisperse packing are investigated to reveal the influence of size diversity in intermittent features of granular materials. Experiments are conducted in an annulus that can contain up to 15 kg of spherical steel balls. The shearing granular medium takes place via the rotation of the upper plate which compresses the material loaded inside the annulus. Fluctuations of compressive force are locally measured at the bottom of the annulus using a piezoelectric sensor. Rapid shear flow experiments are pursued at different compressive forces and shear rates and the sensitivity of fluctuations are then investigated by different means through monodisperse and bidisperse packings. Another important feature of rapid granular shear flows is the formation of ordered structures upon shearing. It requires a certain range for the amount of granular material (uniform size distribution) loaded in the system in order to obtain stable flows. This is studied more deeply in this thesis. The results of the current work bring some new insights into deformation dynamics and intermittency in rapid granular shear flows. The experimental apparatus is modified in comparison to earlier investigations. The measurements produce data for various quantities continuously sampled from the start of shearing to the end. Static failure and dynamic shearing ofa granular medium is investigated. The results of this work revealed some important features of failure dynamics and structure formation in the system. Furthermore, some computer simulations are performed in a 2D annulus to examine the nature of kinetic energy dissipation. It is found that turbulent flow models can statistically represent rapid granular flows with high accuracy. In addition to academic outcomes and scientific publications our results have a number of technological applications associated with grinding, mining and massive grain storages.

Mass transfer and particleinteractions in granular systems and suspension flows

The purpose of this study was to investigate some important features of granular flows and suspension flows by computational simulation methods. Granular materials have been considered as an independent state ofmatter because of their complex behaviors. They sometimes behave like a solid, sometimes like a fluid, and sometimes can contain both phases in equilibrium. The computer simulation of dense shear granular flows of monodisperse, spherical particles shows that the collisional model of contacts yields the coexistence of solid and fluid phases while the frictional model represents a uniform flow of fluid phase. However, a comparison between the stress signals from the simulations and experiments revealed that the collisional model would result a proper match with the experimental evidences. Although the effect of gravity is found to beimportant in sedimentation of solid part, the stick-slip behavior associated with the collisional model looks more similar to that of experiments. The mathematical formulations based on the kinetic theory have been derived for the moderatesolid volume fractions with the assumption of the homogeneity of flow. In orderto make some simulations which can provide such an ideal flow, the simulation of unbounded granular shear flows was performed. Therefore, the homogeneous flow properties could be achieved in the moderate solid volume fractions. A new algorithm, namely the nonequilibrium approach was introduced to show the features of self-diffusion in the granular flows. Using this algorithm a one way flow can beextracted from the entire flow, which not only provides a straightforward calculation of self-diffusion coefficient but also can qualitatively determine the deviation of self-diffusion from the linear law at some regions nearby the wall inbounded flows. Anyhow, the average lateral self-diffusion coefficient, which was calculated by the aforementioned method, showed a desirable agreement with thepredictions of kinetic theory formulation. In the continuation of computer simulation of shear granular flows, some numerical and theoretical investigations were carried out on mass transfer and particle interactions in particulate flows. In this context, the boundary element method and its combination with the spectral method using the special capabilities of wavelets have been introduced as theefficient numerical methods to solve the governing equations of mass transfer in particulate flows. A theoretical formulation of fluid dispersivity in suspension flows revealed that the fluid dispersivity depends upon the fluid properties and particle parameters as well as the fluid-particle and particle-particle interactions.

Steady-state global optimization of metabolic non-linear dynamic models through recasting into power-law canonical models

Background: Design of newly engineered microbial strains for biotechnological purposes would greatly benefit from the development of realistic mathematical models for the processes to be optimized. Such models can then be analyzed and, with the development and application of appropriate optimization techniques, one could identify the modifications that need to be made to the organism in order to achieve the desired biotechnological goal. As appropriate models to perform such an analysis are necessarily non-linear and typically non-convex, finding their global optimum is a challenging task. Canonical modeling techniques, such as Generalized Mass Action (GMA) models based on the power-law formalism, offer a possible solution to this problem because they have a mathematical structure that enables the development of specific algorithms for global optimization. Results: Based on the GMA canonical representation, we have developed in previous works a highly efficient optimization algorithm and a set of related strategies for understanding the evolution of adaptive responses in cellular metabolism. Here, we explore the possibility of recasting kinetic non-linear models into an equivalent GMA model, so that global optimization on the recast GMA model can be performed. With this technique, optimization is greatly facilitated and the results are transposable to the original non-linear problem. This procedure is straightforward for a particular class of non-linear models known as Saturable and Cooperative (SC) models that extend the power-law formalism to deal with saturation and cooperativity. Conclusions: Our results show that recasting non-linear kinetic models into GMA models is indeed an appropriate strategy that helps overcoming some of the numerical difficulties that arise during the global optimization task.

Insertion of a malE B-Galactosidase fusion protein into the envelope of Escherichia coli disrupts biogenesis of outer membrane proteins and processing of inner membrane proteins

The synthesis of a membrane-bound MalE ,B-galactosidase hybrid protein, when induced by growth of Escherichia coli on maltose, leads to inhibition of cell division and eventually a reduced rate of mass increase. In addition, the relative rate of synthesis of outer membrane proteins, but not that of inner membrane proteins, was reduced by about 50%o. Kinetic experiments demonstrated that this reduction coincided with the period of maximum synthesis of the hybrid protein (and another maltose-inducible protein, LamB). The accumulation of this abnormal protein in the envelope therefore appeared specifically to inhibit the synthesis, the assembly of outer membrane proteins, or both, indicating that the hybrid protein blocks some export site or causes the sequestration of some limiting factor(s) involved in the export process. Since the MalE protein is normally located in the periplasm, the results also suggest that the synthesis of periplasmic and outer membrane proteins may involve some steps in common. The reduced rate of synthesis of outer membrane proteins was also accompanied by the accumulation in the envelope of at least one outer membrane protein and at least two inner membrane proteins as higher-molecular-weight forms, indicating that processing (removal of the N-terminal signal sequence) was also disrupted by the presence of the hybrid protein. These results may indicate that the assembly of these membrane proteins is blocked at a relatively late step rather than at the level of primary recognition of some site by the signal sequence. In addition, the results suggest that some step common to the biogenesis of quite different kinds of envelope protein is blocked by the presence of the hybrid protein.

Application of Information Statistical Theory to the Description of the Effect of Heat Conduction on the Chemical Reaction Rate in Gases

The effect of the heat flux on the rate of chemical reaction in dilute gases is shown to be important for reactions characterized by high activation energies and in the presence of very large temperature gradients. This effect, obtained from the second-order terms in the distribution function (similar to those obtained in the Burnett approximation to the solution of the Boltzmann equation), is derived on the basis of information theory. It is shown that the analytical results describing the effect are simpler if the kinetic definition for the nonequilibrium temperature is introduced than if the thermodynamic definition is introduced. The numerical results are nearly the same for both definitions

The Dynamics of Force Fluctuations in Shear Granular Flows

Tämän työn tavoitteena oli tutkia rakeisen materiaalin kinematiikkaa ja rakentaa koelaitteisto rakeisen materiaalin leikkausjännitysvirtauksien tutkimiseen. Kokeellisessa osassa on keskitytty sisäisiin voimaheilahteluihin ja niiden ymmärtämiseen. Teoriaosassa on käyty läpi rakeisen materiaalin yleisiä ominaisuuksia ja lisäksi on esitetty kaksi eri tapaa mallintaa fysikaalisien ominaisuuksien heilahteluja rakeisessa materiaalissa. Nämä kaksi esitettyä mallinnusmenetelmää ovat skalaarinen q-malli ja simulointi. Skalaarinen q-malli määrittelee jokaiseen yksittäiseen rakeeseen kohdistuvan jännityksen, rakeen ollessa osa 2- tai 3-dimensionaalista asetelmaa. Tämän mallin perusidea on kuvata jännityksien epähomogeenisuutta, joka johtuu rakeiden satunnaisasettelusta. Simulointimallinnus perustuu event-driven algoritmiin, missä systeemin dynamiikkaa kuvataan yksittäisillä partikkelien törmäyksillä. Törmäyksien vaiheet ratkaistiin käyttämällä liikemääräyhtälöitä ja restituution määritelmää. Teoriaosuudessa käytiin vielä pieniltä osin läpi syitä jännitysheilahteluihin ja rakeisen materiaalin lukkiintumiseen. Tutkimuslaitteistolla tutkittiin rakeisen materiaalin käyttäytymistä rengasmaisessa leikkausjännitysvirtauksessa. Tutkimusosuuden päätavoitteena oli mitata partikkelien kosketuksista ja törmäyksistä johtuvia hetkellisiä voimaheilahteluja rengastilavuuden pohjalta. Rakeisena materiaalina tutkimuksessa käytettiin teräskuulia. Jännityssignaali ajan funktiona osoittaa suurta heilahtelua, joka voi olla jopa kertalukua keskiarvosta suurempaa. Tällainen suuren amplitudin omaava heilahtelu on merkittävä haittapuoli yleisesti rakeisissa materiaaleissa käytettyjen jatkuvuusmallien kanssa. Tällainen heilahtelu tekee käytetyt jatkuvuusmallit epäpäteviksi. Yleisellä tasolla jännityksien todennäköisyysjakauma on yhtäpitävä skalaarisen q-mallin tuloksien kanssa. Molemmissa tapauksissa todennäköisyysjakaumalla on eksponentiaalinen muoto.

Arvometallien selektiivinen erottaminen kloridiliuoksista

Työssä tutkittiin jalometallien selektiivistä erottamista kloridiliuoksista synteettisten polymeerihartsien avulla. Laboratoriokokeissa keskityttiin tutkimaan kullan erottamista hydrofiilisen polymetakrylaattipohjaisen adsorbentin avulla. Lähtökohtana oli platinarikaste, joka sisälsi kullan lisäksi platinaa, palladiumia, hopeaa, kuparia, rautaa, vismuttia, seleeniä ja telluuria. Mittauksissa tutkittiin eri metallien ja puolimetallien adsorptiota hartsiin tasapaino-, kinetiikka- ja kolonnikokeilla. Työssä käytettiin myös adsorption simulointiin monikomponenttierotuksen dynaamiseen mallintamiseen tarkoitettua tietokoneohjelmaa, johon tarvittavat parametrit estimoitiin kokeellisen datan avulla. Tasapainokokeet yhtä metallia sisältäneistä liuoksista osoittivat, että hartsi adsorboi tehokkaasti kultaa kaikissa tutkituissa suolahappopitoisuuksissa (1-6 M). Kulta muodostaa hartsiin hyvin adsorboituvia tetrakloroauraatti(III)ioneja, [AuCl4]-, jotka ovat erittäin stabiileja pieniin kloridipitoisuuksiin saakka. Suolahappopitoisuudella oli merkitystä ainoastaan raudan adsorptioon, joka kasvoi huomattavasti suolahappopitoisuuden noustessa johtuen raudan taipumuksesta muodostaa hyvin adsorboituvia [FeCl4]--ioneja väkevissä suolahappopitoisuuksissa. Muiden tutkittujen alkuaineiden adsorptiot jäivät alhaisiksi kaikilla suolahappopitoisuuksilla. Rikasteliuoksella tehdyt tasapainokokeet osoittivat, että adsorptiokapasiteetti kullalle riippuu voimakkaasti muista läsnäolevista komponenteista. Kilpaileva adsorptio kuvattiin Langmuir-Freundlich-isotermillä. Kolonnikokeet osoittivat, että hartsi adsorboi kullan lisäksi hieman myös rautaa ja telluuria, jotka saatiin kuitenkin eluoitua hartsista täysin 5 M suolahappopesulla ja sitä seuraavalla 1 M suolahappopesulla. Tehokkaaksi liuokseksi kullan desorboimiseen osoittautui asetonin ja 1 M suolahapon seos. Kolonnierotuksen eri vaiheet pystyttiin tyydyttävästi kuvaamaan simulointimallilla.

Kineettinen tarkastelu keskisakeusalueella toimivalle happivaiheelle

Työn tavoitteena oli tarkastella keskisakeudessa toimivan happidelignifioinnin kinetiikkaa ja vertailla sitä olemassa olevien kineettisten tutkimuksien kanssa. Lisäksi tutkittiin kuinka happivaiheen suodos vaikuttaa delignifiointiasteeseen sekä selektiivisyyteen. Työn kirjallisuusosassa perehdyttiin happidelignifioinnin kineettisiin malleihin, jotka on kehitetty ennustamaan alkalin, paineen, lämpötilan ja ajan vaikutusta ligniinin poistumiselle sekä hiilihydraattiketjujen katkeilulle. Delignifioituminen on nopeaa ensimmäisen kymmenen minuutin aikana ja hidastuu jäännösdelignifiointivaiheessa. Laboratoriotutkimuksissa selvitettiin kuinka prosessimuuttujat vaikuttavat kappaluvun ja viskositeetin muutokseen. Tuloksia tarkasteltiin lineaarisella regressioanalyysillä, jonka avulla muuttujien vaikutukset saadaan selkeästi esille. Happivaiheen ensimmäisessä vaiheessa suurin vaikutus delignifiointiasteeseen sekä viskositeetin alenemiseen on alkaliannoksella. Pienillä hapen annostuksella saavutettiin hyvä selektiivisyys, suurillakin alkaliannoksilla. Lämpötilan vaikutus on lähes olematon, kun paine on pieni. Kun paine kasvaa, kasvaa myös lämpötilan vaikutus delignifiointiasteeseen. Hiilihydraattien depolymerisoituminen lisääntyy lämpötilan kasvaessa lähes lineaarisesti, riippumatta paineesta. Hyvä loppuselektiivisyys saavutetaan kun ensimmäisen vaiheen selektiivisyys on hyvä. Ensimmäisessä vaiheessa saavutettiin hyvä selektiivisyys, kun alkali ei pääse kulumaan loppuun eikä delignifioituminen jatku liian pitkälle. Toisessa vaiheessa selektiivisyys säilyy parhaiten, kun lisättävä alkaliannos on pieni. Lämpötila ja paine vaikuttavat alkaliannokseen verrattuna hyvin vähän loppuselektiivisyyteen.

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.

A simple model to describe the thixotropic behavior of paints.

We propose a simple rheological model to describe the thixotropic behavior of paints, since the classical hysteresis area, which is usually used, is not enough to evaluate thixotropy. The model is based on the assumption that viscosity is a direct measure of the structural level of the paint. The model depends on two equations: the Cross-Carreau equation to describe the equilibrium viscosity and a second order kinetic equation to express the time dependence of viscosity. Two characteristic thixotropic times are differentiated: one for the net structure breakdown, which is defined as a power law function of shear rate, and an other for the net structure buildup, which is not dependent on the shear rate. The knowledge of both kinetic processes can be used to improve the quality and applicability of paints. Five representative commercial protective marine paints are tested. They are based on chlorinated rubber, acrylic, alkyd, vinyl, and epoxy resins. The temperature dependence of the rheological behavior is also studied with the temperature ranging from 5 ºC to 35 ºC. It is found that the paints exhibit both shear thinning and thixotropic behavior. The model fits satisfactorily the thixotropy of the studied paints. It is also able to predict the thixotropy dependence on temperature. Both viscosity and the degree of thixotropy increase as the temperature decreases.

Thomas-Fermi theory for atomic nuclei revisited

The recently developed semiclassical variational Wigner-Kirkwood (VWK) approach is applied to finite nuclei using external potentials and self-consistent mean fields derived from Skyrme inter-actions and from relativistic mean field theory. VWK consist s of the Thomas-Fermi part plus a pure, perturbative h 2 correction. In external potentials, VWK passes through the average of the quantal values of the accumulated level density and total en energy as a function of the Fermi energy. However, there is a problem of overbinding when the energy per particle is displayed as a function of the particle number. The situation is analyzed comparing spherical and deformed harmonic oscillator potentials. In the self-consistent case, we show for Skyrme forces that VWK binding energies are very close to those obtained from extended Thomas-Fermi functionals of h 4 order, pointing to the rapid convergence of the VWK theory. This satisfying result, however, does not cure the overbinding problem, i.e., the semiclassical energies show more binding than they should. This feature is more pronounced in the case of Skyrme forces than with the relativistic mean field approach. However, even in the latter case the shell correction energy for e.g.208 Pb turns out to be only ∼ −6 MeV what is about a factor two or three off the generally accepted value. As an adhoc remedy, increasing the kinetic energy by 2.5%, leads to shell correction energies well acceptable throughout the periodic table. The general importance of the present studies for other finite Fermi systems, self-bound or in external potentials, is pointed out.

Implisiittisen kalibroinnin soveltamismahdollisuudet kemiantekniikassa

Työn taustalla oli tavoite parantaa erään teollisuusprosessin toimintaa ja sen tuottoa mallintamalla reaktiovaiheen alussa tapahtuvan välituotteen muodostumisen reaktiokinetiikka sekä perinteisellä tavalla että implisiittisellä kalibroinnilla. Toisena tavoitteena oli selvittää, kuinka implisiittistä kalibrointia voidaan yleisemmin hyödyntää kemiantekniikassa. Implisiittinen kalibrointi on menetelmä, jolla voidaan ratkaista jonkin teoreettisen mallin parametrit suoraan epäsuorasta mittausdatasta (esimerkiksi spektreistä) lähes kokonaan ilman off-line analyysejä. Tämän työn kirjallisuusosassa on esitetty implisiittisen kalibroinnin toimintaperiaate sekä lyhyesti FTIR-spektrometrian perusteita. Työn kokeellisessa osassa on estimoitu tutkitun välituotteen muodostumisen kineettiset parametrit sekä tavanomaisella parametriestimoinnilla että implisiittisellä kalibroinnilla. Lisäksi kokeellisessa osassa on selvitetty lyhyesti tutkitun prosessin FTIR-spektrien lämpötilariippuvuuksia ja esitetty neljä mahdollista uutta sovelluskohdetta implisiittiselle kalibroinnille. Tavanomaisella parametriestimoinnilla saatiin estimoitua varsin yksiselitteiset arvot kineettisille parametreille. Myös mallin sovitus koedataan on hyvä kolmessa kokeessa viidestä. Parametriestimointi implisiittisellä kalibroinnilla onnistui lupaavasti vaikka tulokset eivät ole aivan niin hyviä kuin tavanomaisessa parametriestimoinnissa. Parhaat tulokset implisiittisessä kalibroinnissa saavutettiin suoralla kalibrointitavalla GRR (Generalized Ridge Regression)-kalibrointimenetelmää käyttämällä.