Modeling of non-isothermalvapor membrane separation with thermodynamic models and generalized mass transfer equations

A rigorous unit operation model is developed for vapor membrane separation. The new model is able to describe temperature, pressure, and concentration dependent permeation as wellreal fluid effects in vapor and gas separation with hydrocarbon selective rubbery polymeric membranes. The permeation through the membrane is described by a separate treatment of sorption and diffusion within the membrane. The chemical engineering thermodynamics is used to describe the equilibrium sorption of vapors and gases in rubbery membranes with equation of state models for polymeric systems. Also a new modification of the UNIFAC model is proposed for this purpose. Various thermodynamic models are extensively compared in order to verify the models' ability to predict and correlate experimental vapor-liquid equilibrium data. The penetrant transport through the selective layer of the membrane is described with the generalized Maxwell-Stefan equations, which are able to account for thebulk flux contribution as well as the diffusive coupling effect. A method is described to compute and correlate binary penetrant¿membrane diffusion coefficients from the experimental permeability coefficients at different temperatures and pressures. A fluid flow model for spiral-wound modules is derived from the conservation equation of mass, momentum, and energy. The conservation equations are presented in a discretized form by using the control volume approach. A combination of the permeation model and the fluid flow model yields the desired rigorous model for vapor membrane separation. The model is implemented into an inhouse process simulator and so vapor membrane separation may be evaluated as an integralpart of a process flowsheet.

Ion-exchange resins as stationary phase in reactive chromatography

Dynamic behavior of bothisothermal and non-isothermal single-column chromatographic reactors with an ion-exchange resin as the stationary phase was investigated. The reactor performance was interpreted by using results obtained when studying the effect of the resin properties on the equilibrium and kinetic phenomena occurring simultaneously in the reactor. Mathematical models were derived for each phenomenon and combined to simulate the chromatographic reactor. The phenomena studied includes phase equilibria in multicomponent liquid mixture¿ion-exchange resin systems, chemicalequilibrium in the presence of a resin catalyst, diffusion of liquids in gel-type and macroporous resins, and chemical reaction kinetics. Above all, attention was paid to the swelling behavior of the resins and how it affects the kinetic phenomena. Several poly(styrene-co-divinylbenzene) resins with different cross-link densities and internal porosities were used. Esterification of acetic acid with ethanol to produce ethyl acetate and water was used as a model reaction system. Choosing an ion-exchange resin with a low cross-link density is beneficial inthe case of the present reaction system: the amount of ethyl acetate as well the ethyl acetate to water mole ratio in the effluent stream increase with decreasing cross-link density. The enhanced performance of the reactor is mainly attributed to increasing reaction rate, which in turn originates from the phase equilibrium behavior of the system. Also mass transfer considerations favor the use ofresins with low cross-link density. The diffusion coefficients of liquids in the gel-type ion-exchange resins were found to fall rapidly when the extent of swelling became low. Glass transition of the polymer was not found to significantlyretard the diffusion in sulfonated PS¿DVB ion-exchange resins. It was also shown that non-isothermal operation of a chromatographic reactor could be used to significantly enhance the reactor performance. In the case of the exothermic modelreaction system and a near-adiabatic column, a positive thermal wave (higher temperature than in the initial state) was found to travel together with the reactive front. This further increased the conversion of the reactants. Diffusion-induced volume changes of the ion-exchange resins were studied in a flow-through cell. It was shown that describing the swelling and shrinking kinetics of the particles calls for a mass transfer model that explicitly includes the limited expansibility of the polymer network. A good description of the process was obtained by combining the generalized Maxwell-Stefan approach and an activity model that was derived from the thermodynamics of polymer solutions and gels. The swelling pressure in the resin phase was evaluated by using a non-Gaussian expression forthe polymer chain length distribution. Dimensional changes of the resin particles necessitate the use of non-standard mathematical tools for dynamic simulations. A transformed coordinate system, where the mass of the polymer was used as a spatial variable, was applied when simulating the chromatographic reactor columns as well as the swelling and shrinking kinetics of the resin particles. Shrinking of the particles in a column leads to formation of dead volume on top of the resin bed. In ordinary Eulerian coordinates, this results in a moving discontinuity that in turn causes numerical difficulties in the solution of the PDE system. The motion of the discontinuity was eliminated by spanning two calculation grids in the column that overlapped at the top of the resin bed. The reactive and non-reactive phase equilibrium data were correlated with a model derived from thethermodynamics of polymer solution and gels. The thermodynamic approach used inthis work is best suited at high degrees of swelling because the polymer matrixmay be in the glassy state when the extent of swelling is low.

Comprehensive Study of Crystal Growth from Solution

Crystal growth is an essential phase in crystallization kinetics. The rate of crystal growth provides significant information for the design and control of crystallization processes; nevertheless, obtaining accurate growth rate data is still challenging due to a number of factors that prevail in crystal growth. In industrial crystallization, crystals are generally grown from multi-componentand multi-particle solutions under complicated hydrodynamic conditions; thus, it is crucial to increase the general understanding of the growth kinetics in these systems. The aim of this work is to develop a model of the crystal growth rate from solution. An extensive literature review of crystal growth focuses on themodelling of growth kinetics and thermodynamics, and new measuring techniques that have been introduced in the field of crystallization. The growth of a singlecrystal is investigated in binary and ternary systems. The binary system consists of potassium dihydrogen phosphate (KDP, crystallizing solute) and water (solvent), and the ternary system includes KDP, water and an organic admixture. The studied admixtures, urea, ethanol and 1-propanol, are employed at relatively highconcentrations (of up to 5.0 molal). The influence of the admixtures on the solution thermodynamics is studied using the Pitzer activity coefficient model. Theprediction method of the ternary solubility in the studied systems is introduced and verified. The growth rate of the KDP (101) face in the studied systems aremeasured in the growth cell as a function of supersaturation, the admixture concentration, the solution velocity over a crystal and temperature. In addition, the surface morphology of the KDP (101) face is studied using ex situ atomic force microscopy (AFM). The crystal growth rate in the ternary systems is modelled on the basis of the two-step growth model that contains the Maxwell-Stefan (MS) equations and a surface-reaction model. This model is used together with measuredcrystal growth rate data to develop a new method for the evaluation of the model parameters. The validation of the model is justified with experiments. The crystal growth rate in an imperfectly mixed suspension crystallizer is investigatedusing computational fluid dynamics (CFD). A solid-liquid suspension flow that includes multi-sized particles is described by the multi-fluid model as well as by a standard k-epsilon turbulence model and an interface momentum transfer model. The local crystal growth rate is determined from calculated flow information in a diffusion-controlled crystal growth regime. The calculated results are evaluated 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.

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.

Kuplaleijupetikattilan dynaaminen malli

Työn tavoitteena oli rakentaa dynaaminen malli kuplaleijupetikattilasta APROS- ohjelmistoa käyttäen. Tarkoituksena oli selvittää kyseisen ohjelmiston soveltuvuutta nykyaikaisen voimalaitoskattilan mallintamiseen. Mallin rakentamisen perustana oli toiminnassa oleva kuplaleijupetillä varustettu voimalaitoskattila. Näin oli käytettävissä riittävä määrä aineistoa mallin rakenteen luomiseen ja valmiin mallin sovittamiseen. Työ on luonteeltaan kaksiosainen. Ensimmäinen osa on kirjallisuusosa, jossa esitellään mallinnuksen kohteena olevaa tekniikkaa. Tekniikasta annetaan kuva esittelemällä perusteoria ja käytännön sovellukset. Lisäksi esitellään kattilassa käytettävät polttoaineet. Kirjallisuusosassa esitellään myös käytettävä APROS-mallinnusohjelmisto. Ohjelmiston laskennan perusteita ei erikseen esitellä. Ne pohjautuvat yleiseen termodynamiikan ja lämmönsiirron teoriaan. Ohjelmiston käytöstä ja sen toiminnasta yleensä annetaan yleisluontoinen selostus. Toisessa osassa mallin rakentaminen esitellään vaiheittain ja siinä järjestyksessä kuin se mallia rakennettaessa tehtiin. Kattilamallin toimintaa testattiin vertaamalla kattilan mitoitustilaan viritettyä mallia takuukokeiden mittaustuloksiin. Lisäksi testattiin mallin toimintaa osakuormalla koeajojakson soveltuvasta osakuormatilasta saatuihin mittausarvoihin. Mallin jatkokehitys pitää sisällään laajamittaisen automaation luomisen ja erilaisten muutostilojen testaamista mallilla.

Lämmöntalteenottokattilan dynaaminen simulointimalli

Työssä mallinnettiin kombivoimalaitoksen lämmöntalteenottokattila Apros-simulointiohjelmalla. Simulointimalli valmistettiin vastaamaan Helsingin Energian Vuosaari B:n voimalaitoksen lämmöntalteenottokattilaa, joka toimii kahdella painetasolla. Kattila on Foster Wheelerin valmistama. Ennen mallinnuksen aloittamista tutustuttiin laitoksen termodynamiikkaan, jolloin saatiin riittävä teoreettinen tieto koko laitoksen toiminnasta. Kattilan reunaehtoina ovat kaasuturbiiniprosessi ja laitoksen höyrykierto. Kaasuturbiini korvattiin laskentayhtälöillä, jotka antavat alkuarvot mm. savukaasun massavirralle ja lämpötilalle ennen kattilaa kaasuturbiinin tehon funktiona. Kattila liitetään höyrykiertoon tuorehöyry- ja syöttövesilinjasta, jolloin reunaehtoina annetaan lämpötilat ja paineet massavirroille. Valmistettua mallia testattiin ylösajo- ja kuormanmuutostilanteessa. Ylösajotilanteessa saatuja laskentatuloksia verrattiin todellisen laitoksen mittaustuloksiin, jolloin varmistuttiin simulointimallin oikeasta fysikaalisesta toiminnasta. Kuormanmuutostilanteissa kaasuturbiinin tehoa muutettiin ja samalla seurattiin kattilan reagointia muutostilanteessa. Kuormanmuutosmittauksessa varmistettiin vielä, että kattila reagoi kuormanmuutokseen oikealla tavalla, eikä muutos aiheuta kattilan toiminnalle haitallista värähtelyä.

Selvityksiä reaktiokalorimetrin soveltuvuudesta hydrometallurgiaan

Tämän työn tavoitteena oli reaktiokalorimetrin käyttöönotto sekä sen käyttökelpoisuuden selvittäminen hydrometallurgisten sovellusten ja erityisesti sinkkisulfidin liuotuksen tutkimiseen. Työn kirjallisuusosassa on käsitelty yleisellä tasolla kalorimetrian ja reaktiokalorimetrian teoriaa, termodynamiikkaa sekä sinkkirikasteen liuotuksen kemiaa. Lisäksi työssä esitellään erilaisten kalorimetrien ja termoanalyyttisten mittauslaitteiden toimintaperiaatteita. Työn kokeellisessa osassa selvitettiin reaktiokalorimetrin mittaustulosten tarkkuutta vesi- kokeiden avulla. Laitteistolla määritettiin myös reaktiolämmöt sinkkisulfidin liukenemisreaktiolle sekä elohopeasuolan saostusreaktiolle. Lisäksi tutkittiin sekoitusnopeuden, lämpötilan ja kiintoainepitoisuuden vaikutusta mittaustuloksiin. Reaktiokalorimetrillä suoritettujen kokeiden perusteella havaittiin, että reaktiolämpöjen absoluuttisten arvojen määrittäminen laitteistolla on käytännössä vaikeaa. Koska reaktiokalorimetrillä pystytään määrittämään vain mittauksen aikana tapahtunut kokonaislämpömuutos, vaikuttavat mahdolliset faasimuutokset ja reaktorin lämpöhäviöt mittaustuloksiin. Näiden tekijöiden vaikutus on pyrittävä eliminoimaan tai niiden vaikutus on tunnettava tarkkaan, jos laitteella halutaan saada luotettavia reaktiolämpömittaustuloksia. Laitteiston mittaustarkkuus huononee huomattavasti, kun reaktorin lämpötila nousee yli 60 °C:een. Laitteistolla mitatut reaktiolämmöt poikkeavat huomattavasti vastaavista kirjallisuuden arvoista. Vedelle määritetyt ominaislämpökapasiteetit poikkeavat kirjallisuuden arvoista enintään 5 alle 90 °C:een lämpötilassa.

Antibioottien poisto vedestä adsorptiolla

Antibiootit ovat yleisessä käytössä olevia lääkeaineita, joilla on kyky hidastaa mikrobien kasvua. Osa lääkeaineesta poistuu elimistöstä muuntumattomana. Koska tavanomainen jäteveden käsittelyprosessi ei riitä poistamaan antibiootteja jätevedestä, ne päätyvät vesistöihin, joissa ne häiritsevät ekosysteemiä ja voivat aiheuttaa vastustuskykyisten bakteerikantojen muodostumisen. Antibiootit voitaisiin poistaa vedestä adsorptiolla. Työssä on esitelty antibioottien adsorptiomekanismeja, vedenkäsittelyssä käytettäviä adsorptioprosesseja ja adsorbentteja, jotka soveltuvat antibioottien erottamiseen vedestä. Lisäksi tarkastellaan yleisesti adsorption kinetiikkaa ja termodynamiikkaa sekä mallinnetaan panoskokeiden tulosten perusteella jatkuvatoimisen adsorptiokolonnin toimintaa.

Economics of combined power plant and pyrolysis bio-oil production

There are reasons of necessity in bio-fuel use and bio-energy fast development. It includes the material about bio-energy technologies, applications and methods. There are basic thermodynamics and economic theories. The economic calculation presents the comparison between two combinations. There are boiler plant below 20 MW in combination with ablative pyrolysis plant for bio-oil production and CHP plant below 100 MW in combination with the RTP pyrolysis bio-oil production technology. It provides a material about wood chips and bio-oil characteristics and explains it nature, presents the situation around the bio-fuel market or bio-fuel trade. There is a description of pyrolysis technologies such as ablative and RTP. The liquid product of the pyrolysis processes is bio-oil. The bio-oil could be different even of the same production process, because of the raw material nature and characteristics. The calculation shows advantages and weaknesses of combinations and obtained a proof of suppositions. The next thing, proven by this work is the fact that to get more efficiency from energy project it is good possibility to built plants in combinations.

Motores de reacción modernos por turbinas de gas

This Bachelor’s thesis investigates the different types of jet engines used nowadays, their performance and applications. The thesis includes a general study of dynamics of fly, engine thermodynamics and contamination.

Review of water electrolysis technologies and design of renewable hydrogen production systems

An electric system based on renewable energy faces challenges concerning the storage and utilization of energy due to the intermittent and seasonal nature of renewable energy sources. Wind and solar photovoltaic power productions are variable and difficult to predict, and thus electricity storage will be needed in the case of basic power production. Hydrogen’s energetic potential lies in its ability and versatility to store chemical energy, to serve as an energy carrier and as feedstock for various industries. Hydrogen is also used e.g. in the production of biofuels. The amount of energy produced during hydrogen combustion is higher than any other fuel’s on a mass basis with a higher-heating-value of 39.4 kWh/kg. However, even though hydrogen is the most abundant element in the universe, on Earth most hydrogen exists in molecular forms such as water. Therefore, hydrogen must be produced and there are various methods to do so. Today, the majority hydrogen comes from fossil fuels, mainly from steam methane reforming, and only about 4 % of global hydrogen comes from water electrolysis. Combination of electrolytic production of hydrogen from water and supply of renewable energy is attracting more interest due to the sustainability and the increased flexibility of the resulting energy system. The preferred option for intermittent hydrogen storage is pressurization in tanks since at ambient conditions the volumetric energy density of hydrogen is low, and pressurized tanks are efficient and affordable when the cycling rate is high. Pressurized hydrogen enables energy storage in larger capacities compared to battery technologies and additionally the energy can be stored for longer periods of time, on a time scale of months. In this thesis, the thermodynamics and electrochemistry associated with water electrolysis are described. The main water electrolysis technologies are presented with state-of-the-art specifications. Finally, a Power-to-Hydrogen infrastructure design for Lappeenranta University of Technology is presented. Laboratory setup for water electrolysis is specified and factors affecting its commissioning in Finland are presented.

Kiertoleijupetikattilateknologia ja sen matemaattinen mallinnus energiantuotantoprosessissa

Työn teoreettisessa osuudessa tehdään katsaus kiertoleijupetiteknologian eri osa-alueisiin: leijupedin virtausdynamiikkaan, hiukkaserottimeen ja kiintoaineen palautusmekanismiin. Myös teknologian historiaa ja muita käyttötarkoituksia energiantuotannon ohella käydään läpi. Termodynamiikkaa sekä lämmönsiirron ja voimalaitosprosessien teoriaa käsitellään mallinnuksessa tarvittavilta osin. Mallinnusosiossa käydään läpi kiertoleijupetihöyrykattilan matemaattisen mallin tekoprosessia. Malli perustuu yleisesti saatavilla oleviin yhtälöihin ja korrelaatioihin. Mallintaminen koostuu höyrykattilan jakamisesta lämpöpintoihin ja niiden mitoittamisesta. Mallissa esitetään myös näkemys siitä, miten lämpö siirtyy savukaasuun ja miten petimateriaalin kierto tapahtuu tulipesässä.

A high efficiency microturbine concept

There is a growing trend towards decentralized electricity and heat production throughout the world. Reciprocating engines and gas turbines have an essential role in the global decentralized energy markets and any improvement in their electrical efficiency has a significant impact from the environmental and economic viewpoints. This paper introduces an inter-cooled and recuperated two-shaft microturbine at 500 kW electric output range. The microturbine is optimized for a realistic combination of the turbine inlet temperature, the recuperation rate and the pressure ratio. The new microturbine design aims to achieve significantly increased performance within the range of microturbines and even competing with the efficiencies achieved in large industrial gas turbines. The simulated electrical efficiency is 45%. Improving the efficiency of combined heat and power (CHP) systems will significantly decrease the emissions and operating costs of decentralized heat and electricity production. Cost-effective, compact and environmentally friendly micro-and small-scale CHP turbine systems with high electrical efficiency will have an opportunity to successfully compete against reciprocating engines, which today are used in heat and power generation all over the world and manufactured in large production series. This paper presents a small-scale gas turbine process, capable of competing with reciprocating engine in terms of electrical efficiency.