Hydrodynamics and Heat Transfer in an Airlift Reactor

Työn teoriaosassa esitetään kirjallisuudessa esiintyviä teoreettisia ja kokeellisia yhtälöitä nesteen nopeuden, kaasun tilavuusosuuden, painehäviön ja lämmönsiirron laskemiseksi. Lisäksi käsitellään airlift-reaktoreiden toimintaa, rakennetta ja teollisia sovelluksia, sekä sekoitusta ja geometrian vaikutusta airlift-reaktoreiden hydrodynaamisiin ominaisuuksiin. Kokeellisessa osassa kuvataan käytetty koelaitteisto ja mittausmenetelmät sekä esitetään saadut koetulokset. Koelaitteisto on viidellä nousuputkella varustettu ulkoisen kierron airlift-reaktori. Kokeellisessa osassa pyritään ratkaisemaan tällaisessa reaktorissa mahdollisesti esiintyviä ongelmia, kuten "slug flown" muodostuminen nousuputkissa sekä fluidien epätasainen jakautuminen nousuputkiin. Lisäksi tutkitaan erilaisten muuttujien, kuten kaasun tilavuusvirran, nesteen viskositeetin, suutinkoon ja nesteen jakoputken rakenteen, vaikutusta kaasun tilavuusosuuteen ja nesteen nopeuteen nousuputkissa. Nesteen nopeudet mitataan merkkiainemenetelmällä ja kaasun tilavuusosuudet manometrimenetelmällä. Lämmönsiirtoa tutkitaan mittaamalla lämpötilaeroja nousuputkissa NiCr-Ni –termoelementeillä. Mittaustulosten perusteella muokataan korrelaatiot kaasun tilavuusosuudelle ja nesteen tyhjäputkinopeudelle. Korrelaatioista lasketut tulokset sopivat kohtuullisen hyvin yhteen mitattujen tulosten kanssa. "Slug flown" ei todettu muodostuvan ongelmaksi 2.5 mPa s pienemmillä viskositeetin arvoilla 2 metriä pitkissä ja 19 mm halkaisijaltaan olevissa putkissa. Lisäksi todettiin, että kaasu- ja nestefaasien jakautumisongelmat voidaan ratkaista rakenteellisesti.

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.

Implementation of population balances in simulation of gas-liquid reactors using a CFD-PBM coupled approach

In bubbly flow simulations, bubble size distribution is an important factor in determination of hydrodynamics. Beside hydrodynamics, it is crucial in the prediction of interfacial area available for mass transfer and in the prediction of reaction rate in gas-liquid reactors such as bubble columns. Solution of population balance equations is a method which can help to model the size distribution by considering continuous bubble coalescence and breakage. Therefore, in Computational Fluid Dynamic simulations it is necessary to couple CFD and Population Balance Model (CFD-PBM) to get reliable distribution. In the current work a CFD-PBM coupled model is implemented as FORTRAN subroutines in ANSYS CFX 10 and it has been tested for bubbly flow. This model uses the idea of Multi Phase Multi Size Group approach which was previously presented by Sha et al. (2006) [18]. The current CFD-PBM coupled method considers inhomogeneous flow field for different bubble size groups in the Eulerian multi-dispersed phase systems. Considering different velocity field for bubbles can give the advantageof more accurate solution of hydrodynamics. It is also an improved method for prediction of bubble size distribution in multiphase flow compared to available commercial packages.

Kuplivan leijukerroksen poltto- ja lämpötekninen tarkastelu

Työn tavoitteena oli luoda menetelmä leijukerroskattilan petin lämpötilojen ja toiminta-alueen arvioimiseksi annetuilla lähtötiedoilla. Toimintatapana oli tutkia polttoaineen ominaisuuksien ja kattilan geometrian vaikutuksia petin energiataseeseen ja tätä kautta petin lämpötilaan kattilan ajotapa huomioiden. Työn teoriaosassa esitellään aluksi biopolttoaineiden ominaisuuksia ja leijukerroskattilan toimintaa petin toimintaympäristön selvittämiseksi. Teoriaosan lopussa käsitellään perusteellisesti petin hydrodynamiikkaa ja reaktiokinettikkaa sekä petin alueen mekaanista mitoitusta. Lisäksi polttoaineen syöttötapa ja palaminen petissä tuodaan esille. Tutkimusosassa selvitetään petin energiatasemallin määrittämistä teoriaan ja diplomityön aikana eri kattiloilla suoritettuihin mittauksiin perustuen. Mittauksista on erillinen selvitys seitsemännessä kappaleessa. Lopuksi kuvaillaan petin energiatasetta mallintavan menetelmän toimintaa sekä siihen syötettäviä lähtötietoja. Tuloksena on saatu suppealla alueella toimiva laajennettavissa oleva petin energiataseen malli. Menetelmän testaamiseen käytetään mittaustietoa kahdelta kattilalta, joita ei olla käytetty petin energiatasemallin virityksessä.

Rheology of high solids surface size starches

Tämän diplomityön tarkoituksena oli tutkia pintaliimatärkkelysten reologista käyttäytymistä korkeissa kuiva-ainepitoisuuksissa. Tarve työn suorittamiselle syntyi kun tutkittiin pintaliimausta filminsiirtopuristimella tavallista korkeammissa kuiva-ainepitoisuuksissa, sileän sauvan ollessa applikointilaitteena. Koska applikointi sileällä sauvalla tapahtuu hydrodynaamisten periaatteiden mukaisesti, sen käyttö edellyttää pintaliimojen reologisten ominaisuuksien tarkkaa tuntemusta ja hallintaa.Kiinnostuksen kohteena olevat ominaisuudet olivat tärkkelysten kuiva-ainepitoisuuden (8 – 30 %) vaikutus viskositeettiin eri lämpötiloissa (20, 30, 40 ja 50 ºC), leikkausnopeus alueella 1 s-1 - 700 000 s-1. Myös tärkkelysten myötörajat määritettiin tutkimuksessa. Viskositeetti eri leikkausnopeusalueilla mitattiin seuraavilla laitteilla: Bohlin VOR (matalat leikkausnopeudet ja myötöraja) ja Hercules HiShear (keskitason leikkausnopeudet) reometrit sekä Eklund kapillaariviskometri (korkeat leikkausno-peudet). Analysoidut tärkkelykset olivat kaksi anionista matalaviskoottista peruna (tärkkelys A) ja ohra (tärkkelys C) tärkkelystä, sekä yksi kationinen korkeaviskoottinen peruna tärkkelys (tärkkelys B). Tutkittujen tärkkelysten Brookfield viskositeetit (100 rpm) olivat (10 % liuos, 60 °C:ssa) tärkkelys A ja C: 25 ± 5 mPas ja tärkkelys B: 100 ± 20 mPas.Tärkkelysliuosten kuiva-ainepitoisuuden noustessa muuttui virtauskäyttäytyminen Newtoniaalisesta leikkausohenevaksi. Leikkausoheneva käyttäytyminen oli voimakkainta tärkkelys B:n kohdalla. Viskositeetti – lämpötila riippuvuus korkeissa leikkausnopeuksissa (esim. 500 000 s-1) oli vähäisempää, mitä oli oletettavissa Brookfield viskositeettiarvojen perusteella. Kaikki tarkkelykset osoittautuivat tiksotrooppisiksi, myös tiksotrooppisuus lisääntyi kuiva-ainepitoisuuden kasvaessa. Tärkkelysten myötörajat osoittautuivat odottamattoman alhaisiksi, kuitenkin varsinkin tärkkelys B:n myötörajat olivat selvästi riippuvaisia lämpötilasta ja kuiva-ainepitoisuudesta. Tutkittujen tärkkelysten virtauskäyttäytyminen oli kirjallisuudessa esitetyn kaltaista. Tärkkelysmolekyylien ketjun pituus oli tärkein tärkkelyksen reologisia ominaisuuksia määrittävä tekijä; mitä matalampi on tärkkelyksen molekyylimassa, sitä matalammat ovat viskositeetti ja myötöraja. Pintaliimauksessa tärkkelysmolekyylien ketjunpituudella on suuri vaikutus ajettavuuteen ja lopputuotteen ominaisuuksiin. Haasteellista pintaliimatärkkelyksen valinnassa on sellaisen yhdistelmän löytäminen, jossa sopivan reologisen käyttäytymisen omaava tärkkelys ja pintaliimatulle paperille tai kartongille asetetut vaatimukset kohtaavat.

Kiertoleijukattilan tulipesän lämmönsiirtoprofiilien määritys ja mallinnus

Tämä diplomityö on osa projektia, joka tähtää ylikriittisellä painealueella toimivan läpivirtaustekniikkaan perustuvan kiertoleijukattilan toteuttamiseen. Diplomityössä määritetään mittausten perusteella vaaka- ja pystysuuntaisia lämmönsiirtoprofiileja kaupallisen kiertoleijukattilan tulipesässä. Lisäksi samaisen kattilan tulipesä mallinnetaan ja mallin parametrejä säädetään niin, että mallin tulokset saadaan vastaamaan mitattuja tuloksia mahdollisimman hyvin. Työn tavoitteena on antaa lisätietoa tulipesässä vapautuvan lämpötehon jakautumisesta tulipesän eri osien kesken. Diplomityön kirjallisuuskatsauksessa käsitellään kiertoleijukattilan tulipesän perusilmiöitä kuten leijutusta, lämmönsiirtoa, hydrodynamiikkaa ja palamista sekä tutustutaan muutamiin tulipesämalleihin. Lisäksi tarkastellaan ylikriittisen läpivirtauskattilan tarjoamia etuja perinteisiin kattilaratkaisuihin nähden. Diplomityössä on myös esitelty työhön liittyvä mittausprojekti.

Mittausmenetelmien kehittaminen kuplakolonnille

Kuplakolonnireaktoreiden CFD-mallinnus on talla hetkella voimakkaasti kehittyva tutkimusalue. Kaksifaasivirtauksen luotettava simulointi ja mallintaminen on haastavaa kuplakolonnireaktorissa tapahtuvien ilmioiden monimutkaisuuden vuoksi. Reaktorin kayttaytymiseen vaikuttavat tekijat, kuten kolonnin hydrodynamiikka ja aineensiirto, tulee tuntea hyvin ennen mallien tekoa. Tassa tyossa on kokeellisesti tutkittu erilaisten mittausmenetelmien soveltuvuutta kuplakolonnin hydrodynamiikan tutkimiseen. Mittausmenetelmissa on keskitytty erityisesti CFD-mallien vaatimiin paikallisiin mittauksiin. Lisaksi tyossa on arvioitu mittausmenetelmien soveltuvuutta j a luotettavuutta CFD-mallien validointiin. Tyon kirjallisuusosassa on perehdytty kuplakolonnireaktorin hydrodynaamiseen kayttaytymiseen ja siihen vaikuttaviin tekijoihin. Naita ovat mm. reaktorityypit, kaasun dispergointi, virtaustyypit ja -alueet, kaasun tilavuusosuus, kaasukuplan koko ja kuplan nousunopeus. Mittauksia tehtiin kahdessa erikokoisessa kuplakolonnissa, joista pienemman halkaisija oli 0,078 m ja suuremman 0,182 m. Molempien kolonnien nestepinnan korkeus oli 4,62 m. Mittaukset tehtiin vesijohtovedella ja epaorgaanisella prosessiliuoksella. Hydrodynaamisista ominaisuuksista mitattiin kaasun tilavuusosuus, kaasukuplan koko seka kaasukuplan nousunopeus. Kaasun tilavuusosuusmittaukset tehtiin paaasiassa paine-eromittauksella ja joissakin tapauksissa pinnanmittausmenetelmalla. Kuplakoko- ja kuplan nousunopeusmittaukset tehtiin suumopeusvideokameralla ja laser Doppler-anemometrilla. Mittauksissa kaytettiin kahdeksaa erilaista kaasunjakolaitetta, joilla selvitettiin kaasunjakolaitteen ominaisuuksien vaikutusta kolonnin hydrodynamiikkaan. Tuloksista havaittiin, etta nestefaasin ominaisuuksilla oli suuri vaikutus kolonnin hydrodynaamiseen kayttaytymiseen. En kaasunjakolaitteilla vesijohtovedella mitatut hydrodynaamiset ominaisuudet eivat poikenneet paljoa toisistaan, kun taas prosessiliuoksella kaasunjakolaitteiden valille saatiin huomattavat erot. Mittausmenetelmista laser Doppler-anemometri ei kaytettavissa olleella optiikalla soveltunut kaasukuplien mittaamiseen. Kuplat olivat menetelmalle liian suuria. Suumopeusvideokamerallaja paine-eromittauksella paastiin hyviin tuloksiin.

Fischer-Tropsch Slurry Bubble Column Reactor Modeling

The literature part of the work reviews overall Fischer-Tropsch process, Fischer-Tropsch reactors and catalysts. Fundamentals of Fischer-Tropsch modeling are also presented. The emphasis is on the reactor unit. Comparison of the reactors and the catalysts is carried out to choose the suitable reactor setup for the modeling work. The effects of the operation conditions are also investigated. Slurry bubble column reactor model operating with cobalt catalyst is developed by taking into account the mass transfer of the reacting components (CO and H2) and the consumption of the reactants in the liquid phase. The effect of hydrostatic pressure and the change in total mole flow rate in gas phase are taken into account in calculation of the solubilities. The hydrodynamics, reaction kinetics and product composition are determined according to literature. The cooling system and furthermore the required heat transfer area and number of cooling tubes are also determined. The model is implemented in Matlab software. Commercial scale reactor setup is modeled and the behavior of the model is investigated. The possible inaccuraries are evaluated and the suggestions for the future work are presented. The model is also integrated to Aspen Plus process simulation software, which enables the usage of the model in more extensive Fischer-Tropsch process simulations. Commercial scale reactor of diameter of 7 m and height of 30 m was modeled. The capacity of the reactor was calculated to be about 9 800 barrels/day with CO conversion of 75 %. The behavior of the model was realistic and results were in the right range. The highest uncertainty to model was estimated to be caused by the determination of the kinetic rate.

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.

A three dimensional study of heat transfer in a supercritical circulating fluidized bed combustor

From the boiler design point of view, it is imperative to know and understand the operation of the boiler. Since comprehensive measurement of a large furnace is impossible, the furnace can be modeled in order to study its behavior and phenomena. This requires the used model to be validated to correspond with the physical furnace behavior. In this thesis, a three dimensional furnace model is validated to match a bituminous coal utilizing, supercritical once-through circulating fluidized bed combustor based on measurement data. The validated model is used for analyzing the furnace heat transfer. Other heat transfer analysis methods are energy balance method based on tube surface temperature measurements and a method based on measured temperature difference between the tube crest and the fin. The latter method was developed in the thesis using Fluent-software. In the theory part, literature is reviewed and the fundamental aspects of circulating fluidized bed are discussed. These aspects are solid particle behavior in fluidization known as hydrodynamics, behavior of fuel and combustion and heat transfer. Fundamental aspects of modeling are also presented.

An experimental and theoretical study of multiphase flow in a circulating fluidized bed

A systematic averaging procedure has been derived in order to obtain an integral form of conservation equations for dispersed multiphase flow, especially applicable to fluidized beds. A similar averaging method is applied further to formulate macroscopic integral equations, which can be used in one-dimensional and macroscopic multi dimensional models. Circulating fluid bed hydrodynamics has been studied experimentally and both macroscopic and microscopic flow profiles have been measured in a cold model. As an application of the theory, the one dimensional model has been used to study mass and momentum conservation of gas and solid in a circulating fluid bed. Axial solid mixing has also been modelled by the one dimensional model and mixing parameters have been evaluated.

Supersaturation-controlled crystallization

Identification of product requirements and quality, together with the management of production are key issues in chemical engineering. Quality control of crystalline products is part of the quality of many industrially manufactured products like paper, paintings, medicines and fertilizers. In most crystallization cases, quality is described with the size, polymorph, shape and purity of the crystal. The chemical composition, hydrodynamics and driving force, together with the operating temperature are in a key position when the properties of a crystalline product are controlled with the crystallization process. This study concentrates on managing the identified properties of a crystalline product with the control of a driving force. The controlling of the driving force can be based on the change of solubility or the change of concentration. Solubility can be changed with temperature, pressure and an antisolvent. The concentration of crystallizing compound, the solute can be changed with the evaporation of the solvent and with the addition of a reagent. The present study focuses on reagent addition and temperature change as methods of changing the level of the driving force. Three control structures for direct control of supersaturation are built, one for cooling crystallization and two for reactive crystallization. Closed loop feedback control structures are based on the measurement of the solute concentration with attenuated total reflection - Fourier transform infrared spectrometer. The details of the reagent feed are analyzed with experimental studies and with results of computational fluid dynamic simulations of the inert particle pulse in the premixer and inert particle injection to the mixing tank. Nucleation in conditions of controlled reactive crystallization is analyzed with Nielsen’s equation of homogeneous nucleation. The resulting control systems, based on regulation of supersaturation, can be used to produce the desired polymorph of an organic product. The polymorph composition of product crystals is controlled repeatably with the decision of a set value of supersaturation level.

Simulation of hydrogen peroxide direct synthesis in a microreactor

A set of models in Aspen plus was built to simulate the direct synthesis process of hydrogen peroxide in a micro-reactor system. This process model can be used to carry out material balance calculation under various experimental conditions. Three thermodynamic property methods were compared by calculating gas solubility and Uniquac-RK method was finally selected for process model. Two different operation modes with corresponding operation conditions were proposed as the starting point of future experiments. Simulations for these two modes were carried out to get the information of material streams. Moreover, some hydrodynamic parameters such as gas/liquid superficial velocity, gas holdup were also calculated with improved process model. These parameters proved the proposed experimental conditions reasonable to some extent. The influence of operation conditions including temperature, pressure and circulation ratio was analyzed for the first operation mode, where pure oxygen was fed into dissolving tank and hydrogen-carbon dioxide mixture was fed into microreactor directly. The preferred operation conditions for the system are low temperature (2°C) and high pressure (30 bar) in dissolving tank. High circulation ratio might be good in the sense that more oxygen could be dissolved and fed into reactor for reactions, but meanwhile hydrodynamics of microreactor should be considered. Furthermore, more operation conditions of reactor gas/liquid feeds in both of two operation modes were proposed to provide guidance for future experiment design and corresponding hydrodynamic parameters were also calculated. Finally, safety issue was considered from thermodynamic point of view and there is no explosion danger at given experimental plan since the released reaction heat will not cause solvent vaporization inside the microchannels. The improvement of process model still needs further study based on the future experimental results.

Microscopic simulation of pigment coating consolidation

The control of coating layer properties is becoming increasingly important as a result of an emerging demand for novel coated paper-based products and an increasing popularity of new coating application methods. The governing mechanisms of microstructure formation dynamics during consolidation and drying are nevertheless, still poorly understood. Some of the difficulties encountered by experimental methods can be overcome by the utilisation of numerical modelling and simulation-based studies of the consolidation process. The objective of this study was to improve the fundamental understanding of pigment coating consolidation and structure formation mechanisms taking place on the microscopic level. Furthermore, it is aimed to relate the impact of process and suspension properties to the microstructure of the coating layer. A mathematical model based on a modified Stokesian dynamics particle simulation technique was developed and applied in several studies of consolidation-related phenomena. The model includes particle-particle and particle-boundary hydrodynamics, colloidal interactions, Born repulsion, and a steric repulsion model. The Brownian motion and a free surface model were incorporated to enable the specific investigation of consolidation and drying. Filter cake stability was simulated in various particle systems, and subjected to a range of base substrate absorption rates and system temperatures. The stability of the filter cake was primarily affected by the absorption rate and size of particles. Temperature was also shown to have an influence. The consolidation of polydisperse systems, with varying wet coating thicknesses, was studied using imposed pilot trial and model-based drying conditions. The results show that drying methods have a clear influence on the microstructure development, on small particle distributions in the coating layer and also on the mobility of particles during consolidation. It is concluded that colloidal properties can significantly impact coating layer shrinkage as well as the internal solids concentration profile. Visualisations of particle system development in time and comparison of systems at different conditions are useful in illustrating coating layer structure formation mechanisms. The results aid in understanding the underlying mechanisms of pigment coating layer consolidation. Guidance is given regarding the relationship between coating process conditions and internal coating slurry properties and their effects on the microstructure of the coating.

Dispersal and metacommunity dynamics in a soft-sediment benthic system : how well is the seafloor connected?

Connectivity depends on rates of dispersal between communities. For marine soft-sediment communities continued small-scale dispersal as post-larvae and as adults can be equally important in maintaining community composition, as initial recruitment of substrate by pelagic larvae. In this thesis post-larval dispersal strategies of benthic invertebrates, as well as mechanisms by which communities are connected were investigated. Such knowledge on dispersal is scarce, due to the difficulties in actually measuring dispersal directly in nature, and dispersal has not previously been quantified in the Baltic Sea. Different trap-types were used underwater to capture dispersing invertebrates at different sites, while in parallel measuring waves and currents. Local community composition was found to change predictably under varying rates of dispersal and physical connectivity (waves and currents). This response was, however, dependent on dispersal-related traits of taxa. Actively dispersing taxa will be relatively better at maintaining their position, as they are not as dependent on hydrodynamic conditions for dispersal and will be less prone to be passively transported by currents. Taxa also dispersed in relative proportions that were distinctly different from resident community composition and a significant proportion (40 %) of taxa were found to lack a planktonic larval life-stage. Community assembly was re-started in a large-scale manipulative field experiment over one year across several sites, which revealed how patterns of community composition (α-, β- and λ-diversity) change depending on rates of dispersal. Results also demonstrated that in response to small-scale disturbance, initial recruitment was by nearby-dominant species after which other species arrived from successively further away. At later assembly time, the number of coexisting species increased beyond what was expected purely by local niche requirements (species sorting), transferring regional differences in community composition (β-diversity) to the local scale (α-diversity, mass effect). Findings of this thesis complement more theoretical studies in metacommunity ecology by demonstrating that understanding how and when individuals disperse relative to underlying environmental heterogeneity is key to interpreting how patterns of diversity change across different spatial scales. Such information from nature is critical when predicting responses to, for example, different types of disturbances or management actions in conservation.