A generalized 1D particle transport method for convection diffusion reaction model

This work is devoted to the development of numerical method to deal with convection diffusion dominated problem with reaction term, non - stiff chemical reaction and stiff chemical reaction. The technique is based on the unifying Eulerian - Lagrangian schemes (particle transport method) under the framework of operator splitting method. In the computational domain, the particle set is assigned to solve the convection reaction subproblem along the characteristic curves created by convective velocity. At each time step, convection, diffusion and reaction terms are solved separately by assuming that, each phenomenon occurs separately in a sequential fashion. Moreover, adaptivities and projection techniques are used to add particles in the regions of high gradients (steep fronts) and discontinuities and transfer a solution from particle set onto grid point respectively. The numerical results show that, the particle transport method has improved the solutions of CDR problems. Nevertheless, the method is time consumer when compared with other classical technique e.g., method of lines. Apart from this advantage, the particle transport method can be used to simulate problems that involve movingsteep/smooth fronts such as separation of two or more elements in the system.

Quantification of Proteins and Cells: Luminometric Nonspecific Particle-Based Methods

New luminometric particle-based methods were developed to quantify protein and to count cells. The developed methods rely on the interaction of the sample with nano- or microparticles and different principles of detection. In fluorescence quenching, timeresolved luminescence resonance energy transfer (TR-LRET), and two-photon excitation fluorescence (TPX) methods, the sample prevents the adsorption of labeled protein to the particles. Depending on the system, the addition of the analyte increases or decreases the luminescence. In the dissociation method, the adsorbed protein protects the Eu(III) chelate on the surface of the particles from dissociation at a low pH. The experimental setups are user-friendly and rapid and do not require hazardous test compounds and elevated temperatures. The sensitivity of the quantification of protein (from 40 to 500 pg bovine serum albumin in a sample) was 20-500-fold better than in most sensitive commercial methods. The quenching method exhibited low protein-to-protein variability and the dissociation method insensitivity to the assay contaminants commonly found in biological samples. Less than ten eukaryotic cells were detected and quantified with all the developed methods under optimized assay conditions. Furthermore, two applications, the method for detection of the aggregation of protein and the cell viability test, were developed by utilizing the TR-LRET method. The detection of the aggregation of protein was allowed at a more than 10,000 times lower concentration, 30 μg/L, compared to the known methods of UV240 absorbance and dynamic light scattering. The TR-LRET method was combined with a nucleic acid assay with cell-impermeable dye to measure the percentage of dead cells in a single tube test with cell counts below 1000 cells/tube.

Utilization of Particle Image velocimetry in PPOOLEX condensation experiments

Particle Image Velocimetry, PIV, is an optical measuring technique to obtain velocity information of a flow in interest. With PIV it is possible to achieve two or three dimensional velocity vector fields from a measurement area instead of a single point in a flow. Measured flow can be either in liquid or in gas form. PIV is nowadays widely applied to flow field studies. The need for PIV is to obtain validation data for Computational Fluid Dynamics calculation programs that has been used to model blow down experiments in PPOOLEX test facility in the Lappeenranta University of Technology. In this thesis PIV and its theoretical background are presented. All the subsystems that can be considered to be part of a PIV system are presented as well with detail. Emphasis is also put to the mathematics behind the image evaluation. The work also included selection and successful testing of a PIV system, as well as the planning of the installation to the PPOOLEX facility. Already in the preliminary testing PIV was found to be good addition to the measuring equipment for Nuclear Safety Research Unit of LUT. The installation to PPOOLEX facility was successful even though there were many restrictions considering it. All parts of the PIV system worked and they were found out to be appropriate for the planned use. Results and observations presented in this thesis are a good background to further PIV use.

Microstructural and environmental effects on the mechanical properties of porous particle-polymer composites

Papper kan anses vara ett av de mest använda materialen i det dagliga livet. Tidskrifter, tidningar, böcker och diverse förpackningar är några exempel på pappersbaserade produkter. Papperets egenskaper måste anpassas till användningsändamålet. En tidskrift kräver t.ex. hög ljushet, opacitet och en slät yta hos papperet, medan dessa egenskaper är mindre viktiga för en dagstidning. Allt tryckpapper behöver vissa mekaniska egenskaper för att tåla vidarebearbetning såsom kalandrering, tryckning och vikning. Man kan bestryka papper för att förbättra dess optiska egenskaper och tryckbarhetsegenskaper. Vid bestrykning appliceras en dispersion av mineralpigment och polymerbindemedel som ett tunt lager på papperets yta. Bestrykningsskiktet kan ses som ett komplext, poröst kompositmaterial som även bidrar till papperets mekaniska egenskaper och dess processerbarhet i diverse konverteringsoperationer. Kravet på framställning av förmånligt papper med tillräckliga styrkeegenskaper ställer allt högre krav på optimeringen av pappersbestrykningsskiktets egenskaper och produktionskostnader. Målet med detta arbete var att förstå sambandet mellan pigmentbestrykningsskiktets mikrostruktur och dess makroskopiska, mekaniska egenskaper. Resultaten visar att adhesionen i gränsytan mellan pigment och bindemedel är kritisk för bestrykningsskiktets förmåga att bära mekanisk belastning. Polära vätskor är vanliga i tryckfärger och kan, eftersom de påverkar syra/bas-interaktionerna mellan pigment och latexbindemedel, försvaga denna adhesion. Resultaten tyder på att ytstyrkan hos bestruket papper kan höjas genom användning av bifunktionella dispergeringsmedel för mineralpigment. Detta medför inbesparingar i pappersproduktionen eftersom mängden bindemedel, den dyraste komponenten i bestrykningsskiktet, kan minskas.

Effect of additives, precursors and process conditions on particle morphology

The objective of the thesis was to develop methods to manufacture and control calcium carbonate crystal nucleation and growth in precipitation process. The work consists of experimental part and literature part that addresses theory of nucleation, crystallization and precipitation. In the experimental part calcium carbonate was precipitated using carbonization reaction. Precipitation was carried out in presence of known morphology controlling agents (anionic polymers and sodium silicate) and by using different operation conditions. Formed material was characterized using SEM images, and its thermal stability was assessed. This work demonstrates that carbon dioxide feeding rate and concentrations of calcium hydroxide and additives can be used to control size, shape and amount of precipitating calcium carbonate.

Particle model for simulating limestone reactions in novel fluidised bed energy applications

The development of carbon capture and storage (CCS) has raised interest towards novel fluidised bed (FB) energy applications. In these applications, limestone can be utilized for S02 and/or CO2 capture. The conditions in the new applications differ from the traditional atmospheric and pressurised circulating fluidised bed (CFB) combustion conditions in which the limestone is successfully used for SO2 capture. In this work, a detailed physical single particle model with a description of the mass and energy transfer inside the particle for limestone was developed. The novelty of this model was to take into account the simultaneous reactions, changing conditions, and the effect of advection. Especially, the capability to study the cyclic behaviour of limestone on both sides of the calcination-carbonation equilibrium curve is important in the novel conditions. The significances of including advection or assuming diffusion control were studied in calcination. Especially, the effect of advection in calcination reaction in the novel combustion atmosphere was shown. The model was tested against experimental data; sulphur capture was studied in a laboratory reactor in different fluidised bed conditions. Different Conversion levels and sulphation patterns were examined in different atmospheres for one limestone type. The Conversion curves were well predicted with the model, and the mechanisms leading to the Conversion patterns were explained with the model simulations. In this work, it was also evaluated whether the transient environment has an effect on the limestone behaviour compared to the averaged conditions and in which conditions the effect is the largest. The difference between the averaged and transient conditions was notable only in the conditions which were close to the calcination-carbonation equilibrium curve. The results of this study suggest that the development of a simplified particle model requires a proper understanding of physical and chemical processes taking place in the particle during the reactions. The results of the study will be required when analysing complex limestone reaction phenomena or when developing the description of limestone behaviour in comprehensive 3D process models. In order to transfer the experimental observations to furnace conditions, the relevant mechanisms that take place need to be understood before the important ones can be selected for 3D process model. This study revealed the sulphur capture behaviour under transient oxy-fuel conditions, which is important when the oxy-fuel CFB process and process model are developed.

Imbuing airborne electronic warfare self-protection : possibilities to disguise

Tutkielman nimi käännettynä suomeksi on Lentokoneiden omasuojajärjestelmien kyllästäminen – mahdollisuuksia kätkeytymiseen. Siinä keskitytään etsimään aktiivisia, ilmauhkaa kohti säteileviä keinoja ilmahyökkääjän, lentokoneen tai helikopterin, omasuojajärjestelmän harhauttamiseksi. Järjestelmät pyrkivät varoittamaan lentäjää laser- ja tutkasäteilyhavainnoista sekä ohjuksen laukaisusta ja lähestymisestä. Omasuojajärjestelmiin kuuluviin vastakeinoihin tutkimus ei keskity. Tavoitteena harhauttamisella on ilmapuolustuksen todellisen määrän ja sijainnin suojaaminen. Epätietoisuus voi saada lentäjän tekemään vääriä johtopäätöksiä. Työ ei fokusoi harhauttamisen taktiseen kehykseen työn kannalta tarpeellista määrää enempää. Tietosuojasyistä tutkimus ei käsittele käytössä olevien omasuojajärjestelmien operatiivisia ominaisuuksia vaan käsittelee teoriaa niiden taustalla. Tutkimus tehtiin käyttäen taustatutkimukselle enemmän tyypillistä kirjallisuustutkimusta. Jo tutki-muksen alkuvaiheessa oli oletettavaa, että elektronisiin sensoreihin perustuvaa omasuojajärjestelmää pystytään harhauttamaan. Systeemianalysoinnilla pyrittiin löytämään vastauksia tähän olettamukseen. Omasuojajärjestelmästä muodostettiin jo tutkimuksen varhaisessa vaiheessa malli. Tieteellistä kirjallisuutta omasuojajärjestelmistä on olemassa jonkin verran, ja niistä löydettyä tietoa omasuojajärjestelmien ominaisuuksista ja sensoreista yhdistettiin malliin niin, että siitä saatiin mahdollisimman tarkka systeemin kuvaus todellisesta omasuojajärjestelmästä. Analyysin tavoitteena oli löytää niitä kriteereitä, joilla omasuojajärjestelmä saataisiin kohtuullisen tehokkaasti uskomaan harhautusta oikeaksi hälytykseksi. Ohjuksen laukaisusta varoittava sensori perustuu ohjuksen moottorin muodostaman pilven lämpö-säteilyyn. Säteily kuitenkin muuttuu lennon eri vaiheissa, mikä tuottaa haasteita järjestelmälle. Millimetrialueen tutkan käyttö varoittimen sensorina on myös yksi vaihtoehto. Laser-varoittimet toimivat koko sillä taajuusalueella, mitä sotilaskäytössä tulenjohtamiseen ja ohjusten ohjaamiseen käytetään. Tutkavaroittimen tutkimus on vielä kesken. Löydettyjä tuloksia analysoimalla tulen tässä vaiheessa hieman ristiriitaisiin tuloksiin. Lämpösäteilyn käytön suurin haaste on sen eteneminen ilmakehässä. Varoittimen tavoitekaan ei ole toimia kymmentä kilometriä pidemmälle. Yksi mahdollinen ratkaisu on suunnitella ja toteuttaa raketti, jonka tuottaa lämpösäteilyä kuten tietty puolustavan joukon käytössä oleva ohjus. Jos koneessa on kuitenkin myös millimetrialueen tutka tukemassa varoitusjärjestelmää, hankaloituu rakettiharhautus merkittävästi, koska sen pitäisi oletettavasti olla lentokoneen kanssa melko tarkasti kohtaavalla reitillä. Laser-varoitin on ilmeisesti herkin järjestelmistä, koska se voi tietyissä olosuhteissa ja varsinkin matalalla lentokorkeudella aiheuttaa paljon vääriä hälytyksiä ilman tarkoituksellista harhauttamista. Laserin käyttö yhdistettynä raketin laukaisuun saattaisi tuottaa halutun tuloksen. Tutkavaroittimen harhautus onnistuu, jos valelaitteen signaali on uskottavan tarkka.

Open publishing: lessons from particle physics

Syksy Räsänen's presentation at Kirjastoverkkopäivät, Helsinki 21.10.2015.

Use of airborne laser scanner data in demanding forest conditions

Most of the applications of airborne laser scanner data to forestry require that the point cloud be normalized, i.e., each point represents height from the ground instead of elevation. To normalize the point cloud, a digital terrain model (DTM), which is derived from the ground returns in the point cloud, is employed. Unfortunately, extracting accurate DTMs from airborne laser scanner data is a challenging task, especially in tropical forests where the canopy is normally very thick (partially closed), leading to a situation in which only a limited number of laser pulses reach the ground. Therefore, robust algorithms for extracting accurate DTMs in low-ground-point-densitysituations are needed in order to realize the full potential of airborne laser scanner data to forestry. The objective of this thesis is to develop algorithms for processing airborne laser scanner data in order to: (1) extract DTMs in demanding forest conditions (complex terrain and low number of ground points) for applications in forestry; (2) estimate canopy base height (CBH) for forest fire behavior modeling; and (3) assess the robustness of LiDAR-based high-resolution biomass estimation models against different field plot designs. Here, the aim is to find out if field plot data gathered by professional foresters can be combined with field plot data gathered by professionally trained community foresters and used in LiDAR-based high-resolution biomass estimation modeling without affecting prediction performance. The question of interest in this case is whether or not the local forest communities can achieve the level technical proficiency required for accurate forest monitoring. The algorithms for extracting DTMs from LiDAR point clouds presented in this thesis address the challenges of extracting DTMs in low-ground-point situations and in complex terrain while the algorithm for CBH estimation addresses the challenge of variations in the distribution of points in the LiDAR point cloud caused by things like variations in tree species and season of data acquisition. These algorithms are adaptive (with respect to point cloud characteristics) and exhibit a high degree of tolerance to variations in the density and distribution of points in the LiDAR point cloud. Results of comparison with existing DTM extraction algorithms showed that DTM extraction algorithms proposed in this thesis performed better with respect to accuracy of estimating tree heights from airborne laser scanner data. On the other hand, the proposed DTM extraction algorithms, being mostly based on trend surface interpolation, can not retain small artifacts in the terrain (e.g., bumps, small hills and depressions). Therefore, the DTMs generated by these algorithms are only suitable for forestry applications where the primary objective is to estimate tree heights from normalized airborne laser scanner data. On the other hand, the algorithm for estimating CBH proposed in this thesis is based on the idea of moving voxel in which gaps (openings in the canopy) which act as fuel breaks are located and their height is estimated. Test results showed a slight improvement in CBH estimation accuracy over existing CBH estimation methods which are based on height percentiles in the airborne laser scanner data. However, being based on the idea of moving voxel, this algorithm has one main advantage over existing CBH estimation methods in the context of forest fire modeling: it has great potential in providing information about vertical fuel continuity. This information can be used to create vertical fuel continuity maps which can provide more realistic information on the risk of crown fires compared to CBH.

Stochastic particle models: mean reversion and burgers dynamics. An application to commodity markets

The aim of this study is to propose a stochastic model for commodity markets linked with the Burgers equation from fluid dynamics. We construct a stochastic particles method for commodity markets, in which particles represent market participants. A discontinuity in the model is included through an interacting kernel equal to the Heaviside function and its link with the Burgers equation is given. The Burgers equation and the connection of this model with stochastic differential equations are also studied. Further, based on the law of large numbers, we prove the convergence, for large N, of a system of stochastic differential equations describing the evolution of the prices of N traders to a deterministic partial differential equation of Burgers type. Numerical experiments highlight the success of the new proposal in modeling some commodity markets, and this is confirmed by the ability of the model to reproduce price spikes when their effects occur in a sufficiently long period of time.