Sulakaasutustutkimuslaitteiston kehittäminen ja testaus

Työ liittyy Lappeenrannan teknillisessä yliopistossa voimalaitostekniikan laboratoriossa vuonna 2004 aloitettuun sulapoltto- ja sulakaasutustutkimukseen. Päätavoitteena on suunnitella, rakentaa ja testata pienen kokoluokan sulakaasutuslaitetta. Kirjallisuusosassa luodaan lyhyt katsaus viime vuosina muualla tehtyyn sulapoltto- ja sulakaasutustutkimukseen. Laskentaosassa mallinnetaan sulapoltto- ja sulakaasutusprosessia eri polttoainekoostumuksilla. Tutkimusosassa suunnitellaan ja rakennetaan sulakaasutuskoelaitteisto kaasutuskokeita varten. Lopussa on kehittyneemmän tutkimuslaitteiston alustavaa suunnittelua.

Contribution of isotope ratio mass spectrometry to the investigation of improvised explosives: isotopic study of black powders ans ammonium nitrates

The establishment of legislative rules about explosives in the eighties has reduced the illicit use of military and civilian explosives. However, bomb-makers have rapidly taken advantage of substances easily accessible and intended for licit uses to produce their own explosives. This change in strategy has given rise to an increase of improvised explosive charges, which is moreover assisted by the ease of implementation of the recipes, widely available through open sources. While the nature of the explosive charges has evolved, instrumental methods currently used in routine, although more sensitive than before, have a limited power of discrimination and allow mostly the determination of the chemical nature of the substance. Isotope ratio mass spectrometry (IRMS) has been applied to a wide range of forensic materials. Conclusions drawn from the majority of the studies stress its high power of discrimination. Preliminary studies conducted so far on the isotopic analysis of intact explosives (pre-blast) have shown that samples with the same chemical composition and coming from different sources could be differentiated. The measurement of stable isotope ratios appears therefore as a new and remarkable analytical tool for the discrimination or the identification of a substance with a definite source. However, much research is still needed to assess the validity of the results in order to use them either in an operational prospect or in court. Through the isotopic study of black powders and ammonium nitrates, this research aims at evaluating the contribution of isotope ratio mass spectrometry to the investigation of explosives, both from a pre-blast and from a post-blast approach. More specifically, the goal of the research is to provide additional elements necessary to a valid interpretation of the results, when used in explosives investigation. This work includes a fundamental study on the variability of the isotopic profile of black powder and ammonium nitrate in both space and time. On one hand, the inter-variability between manufacturers and, particularly, the intra-variability within a manufacturer has been studied. On the other hand, the stability of the isotopic profile over time has been evaluated through the aging of these substances exposed to different environmental conditions. The second part of this project considers the applicability of this high-precision technology to traces and residues of explosives, taking account of the characteristics specific to the field, including their sampling, a probable isotopic fractionation during the explosion, and the interferences with the matrix of the site.

Using digital RNA counting and flow cytometry to compare mRNA with protein expression in acute leukemias.

BACKGROUND: The diagnosis of malignant hematologic diseases has become increasingly complex during the last decade. It is based on the interpretation of results from different laboratory analyses, which range from microscopy to gene expression profiling. Recently, a method for the analysis of RNA phenotypes has been developed, the nCounter technology (Nanostring® Technologies), which allows for simultaneous quantification of hundreds of RNA molecules in biological samples. We evaluated this technique in a Swiss multi-center study on eighty-six samples from acute leukemia patients. METHODS: mRNA and protein profiles were established for normal peripheral blood and bone marrow samples. Signal intensities of the various tested antigens with surface expression were similar to those found in previously performed Affymetrix microarray analyses. Acute leukemia samples were analyzed for a set of twenty-two validated antigens and the Pearson Correlation Coefficient for nCounter and flow cytometry results was calculated. RESULTS: Highly significant values between 0.40 and 0.97 were found for the twenty-two antigens tested. A second correlation analysis performed on a per sample basis resulted in concordant results between flow cytometry and nCounter in 44-100% of the antigens tested (mean = 76%), depending on the number of blasts present in a sample, the homogeneity of the blast population, and the type of leukemia (AML or ALL). CONCLUSIONS: The nCounter technology allows for fast and easy depiction of a mRNA profile from hematologic samples. This technology has the potential to become a valuable tool for the diagnosis of acute leukemias, in addition to multi-color flow cytometry.

Dysmegakaryopoiesis predicting response to therapy in acute myeloid leukaemia. A histologic and clinical study

With standard induction therapy between 50 to 85% of patients with Acute Myeloid Leukaemia (AML) achieve Complete Remission (CR). We investigated whether any morphological feature of bone marrow (BM) plastic embedded biopsies could predict failure of therapy. We reviewed BM plastic embedded biopsies from 54 adult patients presenting with untreated AML. The main histologic parameters analysed were cellularity, dysmegakaryopoiesis (DysM), percentage of marrow blasts and fibrosis. CR was obtained in 34 of 49 treated patients (69%). The rate of CR was significantly lower in the group of patients presenting with DysM: CR was achieved in 54% of the 28 treated patients with DysM and in 90% of the 21 treated patients without DysM (p less than 0.02). Patients with DysM had a significantly lower blood count and bone marrow blasts at presentation. Median age was not significantly different in the 2 groups. Cellularity and fibrosis were not predictive. DysM may be the hallmark of an AML subgroup with distinct clinical behaviour and lower rate of CR with conventional therapy. DysM should be carefully looked for on BM marrow biopsies and aspirate from AML patients at diagnosis.

Dynamic Model of a Small Once-Through Boiler

In this study, a model for the unsteady dynamic behaviour of a once-through counter flow boiler that uses an organic working fluid is presented. The boiler is a compact waste-heat boiler without a furnace and it has a preheater, a vaporiser and a superheater. The relative lengths of the boiler parts vary with the operating conditions since they are all parts of a single tube. The present research is a part of a study on the unsteady dynamics of an organic Rankine cycle power plant and it will be a part of a dynamic process model. The boiler model is presented using a selected example case that uses toluene as the process fluid and flue gas from natural gas combustion as the heat source. The dynamic behaviour of the boiler means transition from the steady initial state towards another steady state that corresponds to the changed process conditions. The solution method chosen was to find such a pressure of the process fluid that the mass of the process fluid in the boiler equals the mass calculated using the mass flows into and out of the boiler during a time step, using the finite difference method. A special method of fast calculation of the thermal properties has been used, because most of the calculation time is spent in calculating the fluid properties. The boiler was divided into elements. The values of the thermodynamic properties and mass flows were calculated in the nodes that connect the elements. Dynamic behaviour was limited to the process fluid and tube wall, and the heat source was regarded as to be steady. The elements that connect the preheater to thevaporiser and the vaporiser to the superheater were treated in a special way that takes into account a flexible change from one part to the other. The model consists of the calculation of the steady state initial distribution of the variables in the nodes, and the calculation of these nodal values in a dynamic state. The initial state of the boiler was received from a steady process model that isnot a part of the boiler model. The known boundary values that may vary during the dynamic calculation were the inlet temperature and mass flow rates of both the heat source and the process fluid. A brief examination of the oscillation around a steady state, the so-called Ledinegg instability, was done. This examination showed that the pressure drop in the boiler is a third degree polynomial of the mass flow rate, and the stability criterion is a second degree polynomial of the enthalpy change in the preheater. The numerical examination showed that oscillations did not exist in the example case. The dynamic boiler model was analysed for linear and step changes of the entering fluid temperatures and flow rates.The problem for verifying the correctness of the achieved results was that there was no possibility o compare them with measurements. This is why the only way was to determine whether the obtained results were intuitively reasonable and the results changed logically when the boundary conditions were changed. The numerical stability was checked in a test run in which there was no change in input values. The differences compared with the initial values were so small that the effects of numerical oscillations were negligible. The heat source side tests showed that the model gives results that are logical in the directions of the changes, and the order of magnitude of the timescale of changes is also as expected. The results of the tests on the process fluid side showed that the model gives reasonable results both on the temperature changes that cause small alterations in the process state and on mass flow rate changes causing very great alterations. The test runs showed that the dynamic model has no problems in calculating cases in which temperature of the entering heat source suddenly goes below that of the tube wall or the process fluid.

Ash Particle Erosion on Steam Boiler Convective Section

In this study, equations for the calculation of erosion wear caused by ash particles on convective heat exchanger tubes of steam boilers are presented. Anew, three-dimensional test arrangement was used in the testing of the erosion wear of convective heat exchanger tubes of steam boilers. When using the sleeve-method, three different tube materials and three tube constructions could be tested. New results were obtained from the analyses. The main mechanisms of erosionwear phenomena and erosion wear as a function of collision conditions and material properties have been studied. Properties of fossil fuels have also been presented. When burning solid fuels, such as pulverized coal and peat in steam boilers, most of the ash is entrained by the flue gas in the furnace. In bubbling andcirculating fluidized bed boilers, particle concentration in the flue gas is high because of bed material entrained in the flue gas. Hard particles, such as sharp edged quartz crystals, cause erosion wear when colliding on convective heat exchanger tubes and on the rear wall of the steam boiler. The most important ways to reduce erosion wear in steam boilers is to keep the velocity of the flue gas moderate and prevent channelling of the ash flow in a certain part of the cross section of the flue gas channel, especially near the back wall. One can do this by constructing the boiler with the following components. Screen plates can beused to make the velocity and ash flow distributions more even at the cross-section of the channel. Shield plates and plate type constructions in superheaters can also be used. Erosion testing was conducted with three types of tube constructions: a one tube row, an inline tube bank with six tube rows, and a staggered tube bank with six tube rows. Three flow velocities and two particle concentrations were used in the tests, which were carried out at room temperature. Three particle materials were used: quartz, coal ash and peat ash particles. Mass loss, diameter loss and wall thickness loss measurements of the test sleeves were taken. Erosion wear as a function of flow conditions, tube material and tube construction was analyzed by single-variable linear regression analysis. In developing the erosion wear calculation equations, multi-variable linear regression analysis was used. In the staggered tube bank, erosion wear had a maximum value in a tube row 2 and a local maximum in row 5. In rows 3, 4 and 6, the erosion rate was low. On the other hand, in the in-line tube bank the minimum erosion rate occurred in tube row 2 and in further rows the erosion had an increasing value, so that in a six row tube bank, the maximum value occurred in row 6.

Propagation of the Ignition Front against Airflow in Packed Beds of Wood Particles

Combustion of wood is increasing because of the needs of decreasing the emissions of carbon dioxide and the amount of waste going to landfills. Wood based fuels are often scattered on a large area. The transport distances should be short enough to prevent too high costs, and so the size of heating and power plants using wood fuels is often rather small. Combustion technologies of small-size units have to be developed to reach efficient and environmentally friendly energy production. Furnaces that use different packed bed combustion or gasification techniques areoften most economic in small-scale energy production. Ignition front propagation rate affects the stability, heat release rate and emissions of packed bed combustion. Ignition front propagation against airflow in packed beds of wood fuels has been studied. The research has been carried out mainly experimentally. Theoretical aspects have been considered to draw conclusions about the experimental results. The effects of airflow rate, moisture content of the fuel, size, shape and density of particles, and porosity of the bed on the propagation rate of the ignition front have been studied. The experiments were carried out in a pot furnace. The fuels used in the experiments were mainly real wood fuels that are often burned in the production of energy. The fuel types were thin wood chips, saw dust, shavings, wood chips, and pellets with different sizes. Also a few mixturesof the above were tested. Increase in the moisture content of the fuel decreases the propagation rates of the ignition front and makes the range of possible airflow rates narrower because of the energy needed for the evaporation of water and the dilution of volatile gases due to evaporated steam. Increase in the airflow rate increases the ignition rate until a maximum rate of propagation is reached after which it decreases. The maximum flame propagation rate is not always reached in stoichiometric combustion conditions. Increase in particle size and density transfers the optimum airflow rate towards fuel lean conditions. Mixing of small and large particles is often advantageous, because small particles make itpossible to reach the maximum ignition rate in fuel rich conditions, and large particles widen the range of possible airflow rates. A correlation was found forthe maximum rate of ignition front propagation in different wood fuels. According to the correlation, the maximum ignition mass flux is increased when the sphericity of the particles and the porosity of the bed are increased and the moisture content of the fuel is decreased. Another fit was found between sphericity and porosity. Increase in sphericity decreases the porosity of the bed. The reasons of the observed results are discussed.

Polttoaineen reaktiivisuuden määrittäminen kiertoleijupoltossa

Kehitettäessä leijupolttotekniikkaa entistä ympäristöystävällisemmäksi ja tehokkaammaksi tarvitaan lisää tietoa polttoaineen käyttäytymisestä tulipesässä. Polttoaineen palamisprofiili ja reaktiivisuus vaikuttavat oleellisesti esimerkiksi voimalaitoskattilan lämmönsiirtopintojen sijoitteluun ja suunnitteluun sekä ohjausjärjestelmän toteutukseen. Varsinkin monipolttoainekattiloilla ohjausjärjestelmän toimivuus joutuu koetukselle esimerkiksi kuormanmuutostilanteissa ja äkillisissä polttoaineen syöttöhäiriöissä. Tässä työssä on aluksi tutustuttu kiertoleijupolton ilmiöihin ja niiden matemaattiseen mallintamiseen. Lisäksi esitetään katsaus eri prosessiolosuhteiden vaikutuksesta palamiseen kiertoleijuolosuhteissa. Työn tutkimusosassa kehitettiin menetelmä polttoaineen reaktiivisuuden määrittämiseksi kiertoleijupoltossa. Kehitetty menetelmä koostuu koesarjasta ja matemaattisesta simulointimallista. Koetoiminta suoritettiin VTT Energian laboratoriokokoluokan kiertoleijukoelaitteella. Koesarja polttoaineen reaktiivisuuden määrittämiseksi sisältää eri kaasukomponenttien profiilimittauksia ja dynaamisia muutoskokeita. Menetelmän avulla voidaan tutkia eri polttoaineiden reaktiivisuuksia sekä polttoaineen reaktiivisuuden ja tietyn prosessiolosuhteen välistä riippuvuutta. Suorittamalla koeajomatriisin mukaiset kokeet tietyissä prosessiolosuhteissa voidaan polttoaineen reaktiivisuus selvittää koetulosten ja simulointimallin perusteella.

Pienten biopolttolaitosten modulointi

Biopolttoaineiden käyttö energiantuotannossa lisääntyy ja syrjäyttää samalla fossiilisia polttoaineita. Etenkin puupohjaisten polttoaineiden käyttö pienenergiantuotannossa on lisääntynyt nopeasti. Tätä kehitystä puoltavat tiukentuneet päästörajat ja yleinen mielipide fossiilisia polttoaineita vastaan sekä biopolttoaineiden verottomuus. Tästä johtuen biopolttoainetta käyttävien laitosten kysyntä kasvaa jatkuvasti. Kysynnän kasvaessa biopolttolaitosten valmistajien on pystyttävä nopeaan, tehokkaaseen ja laadukkaaseen toimintaan, sekä täyttämään asiakkaan tarpeet pystyäkseen vastaamaan alan kasvavaan kilpailuun. Tässä diplomityössä kehitetään standardoitua moduulirakennetta Sermet BioGrate ja BioGrate compact kattilalaitoksille. Moduloinnin tarkoituksena on helpottaa laitoskonseptin kokoonpanon suunnittelua, tarjousten ja tarjouslaskelmien tekemistä sekä yhtenäistää toimintatapoja, mikä vähentää työmäärää ja nopeuttaa toimituksia.

Kerrosleijukattilan sekundääri-ilmansyöttöjärjestelmän kehittäminen

Työn tarkoituksena oli kehittää kerrosleijukattilan sekundääri-ilmansyöttöä. Työssä tutkittiin tulipesässä ilmasuihkujen käyttäytymiseen vaikuttavia tekijöitä sekä eri ilmansyöttömallien vaikutuksia tulipesäolosuhteisiin. Sekundääri- ja tertiääri-ilmasuihkujen tehtävänä on sekoittaa tulipesän kaasuja sekä tuoda happi leijupetin yläpuolelle haihtuneiden aineiden palamisen loppuun saattamiseksi. Kaasujen sekoittumiseen vaikuttavat ilmasuihkun ja ristivirtauksen liikemäärien suhde, suutinten koko ja sijoittelu toisiinsa nähden. Ilmasuihku saavuttaa paremman tunkeutuvuuden ja siten tehokkaamman sekoittumisen suuttimen koon kasvaessa. Lisäksi tunkeutuvuus paranee suutinten välisen etäisyyden sekä ilmasuihkun ja ristivirtauksen liikemäärien suhteen kasvaessa. Optimaalisen suutinten välisen etäisyyden osoitettiin riippuvan suuttimen ja tulipesän koosta sekä ilmasuihkun ja ristivirtauksen liikemäärien suhteesta. Saatujen tulosten mukaan sekundääri- ja tertiääri-ilmatasoilla tulee suosia suuria harvaan ja lomittain sijoitettuja suuttimia. Ilmansyöttömalleja tutkittiin numeerisesti mallintamalla, ja saatujen tulosten perusteella tehokkain sekoittuminen saavutettiin sijoittamalla sekundääri-ilmasuuttimet tulipesän kahdelle vastakkaiselle seinälle.

Kerrostumien muodostumisen määritys sondin lämmönsiirron perusteella

Diplomityön tavoitteena oli tutkia kerrostumien muodostumista voimalaitoskattiloiden lämpöpinnoille ja niiden vaikutusta lämmönsiirtoon. Kerrostumien vaikutusta lämpöpintojen lämpötiloihin tutkittiin kerrostumasondin avulla kierto- ja kerrosleijukattiloissa. Saadun mittausaineiston muodostettiin mittaustilannetta vastaava laskentamalli, jonka avulla pystytään selvittämään voimalaitoksen lämpöpintojen likaantumista käyttöperiodin aikana. Työn alussa tarkasteltiin kerrostumien muodostumismekanismeja lämmönsiirtopinnoille, kerrostumien ominaisuuksia ja virtausympäristön vaikutusta lämmönsiirtoon. Poltossa syntyvä tuhka kiinnittyy eri mekanismeilla kattilan seinille ja lämpöpinnoille riippuen polttoaineesta ja ympäröivistä olosuhteista. Likakerrosten muodostuminen lämmönsiirtimen pinnalle pienentää putkeen siirtyvää lämpövirtaa ja alentaa kattilan hyötysuhdetta. Kerrostumien lämmönjohtumisesta ja mikrorakenteesta on hyvin vähän kokeellisen tietoa, mikä vaikeuttaa reaalisten mallien muodostamista kattilan käyttäytymisestä sen likaantuessa.Työn kokeellisessa osassa tarkasteltiin kahdella eri kattilalla tehtyjä mittauksia. Mittaukset tehtiin kerrostumasondilla, jonka lämmönsiirrosta luotiin laskentamalli SIMULINK-simulointiohjelmalla. Mittaustuloksina saatiin kolmen eri pinnan lämpötilat, jotka muuttuivat kerrostuman ja jäähdytyksen vaikutuksesta. Laskentamallista muokattiin mittaustilanteita vastaava, jolloin lämpötilamuutoksista nähdään likakerroksen ominaisuuksien vaikutus lämmönsiirtoon. Sondista muodostettiin myös FLUENT-malli, jolla tarkasteltiin yksittäisen putken virtauskenttää sekä kahden lähekkäin olevan putken vaikutusta virtaukseen.

Rikinpoiston tehostaminen leijukerroskattiloissa vähärikkisillä polttoaineilla

Leijukerroskattiloissa rikinpoisto tapahtuu useimmiten käyttämällä primäärisiä rikinpoistomenetelmiä. Tämä tarkoittaa jonkin rikkiä sitovan lisäaineen, yleensä kalkkikiven, syöttämistä tulipesään. Tämä menetelmä soveltuu hyvin leijukerrospolttoon sen edullisuuden ja tehokkuuden takia.Tulevaisuudessa tullaan kuitenkin tiukentamaan voimalaitosten rikkidioksidipäästöjen rajoja. Tällöin myös vähärikkisiä polttoaineita käyttävien voimalaitosten tulee alentaa rikkipäästöjään. Nykyisiä menetelmiä käyttäen nousisivat kustannukset korkeaksi saavutettavaan hyötyyn nähden. Tämän takia tullaan tulevaisuudessa tarvitsemaan erilaisia rikinpoiston tehostamiskeinoja.Tämän työn tarkoituksena on selvittää rikinpoistoon ja sen tehokkuuteen vaikuttavia tekijöitä sekä tutkia erilaisia keinoja rikinpoiston tehostamiseksi. Kirjallisuusosassa keskitytään selvittämään nykyiset käytössä tai tutkimusasteella olevat keinot rikinpoiston tehostamiseksi. Kokeellisessa osuudessa keskitytään lentotuhkan kierrättämisen hyödyntämiseen rikinpoistossa. Tehtävissä kokeissa tavoitteena oli rikinpoistoasteen parantaminen sekä puhtaan kalkkikiven kulutuksen vähentäminen käyttämällä hyväksi lentotuhkan sisältämää reagoimatonta kalsiumoksidia.

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.

Leijukerroslämmönsiirtimen leijutusjärjestelyjen säätö ja optimointi

Energiantuotantoa koskeva jatkuvasti tiukentuva lainsäädäntö ja yleinen tarve polttaa yhä vaativampia polttoaineita asettaa leijukerroskattiloille suuria haasteita. Eräs käyttökelpoinen työkalu tulevaisuuden haasteisiin vastaamisessa on säädettävä leijukerroslämmönsiirrin. Tiheässä leijukerroksessa toimiva lämmönsiirrin on vähemmän alttiina korroosiolle kuin savukaasukanavistossa sijaitsevat rakenteet. Rakenne tuo tarpeellista säädettävyyttä koko kattilalle, mikäli lämmönsiirtimen antamaa tehoa voidaan säätää tehokkaasti. Erityisen hyvä tilanne saavutetaan, jos lämmönsiirrin kykenee käyttämään hyväkseen tulipesän sisäistä kiintoainekiertoa. Tässä työssä tutkitaan kiertoleijukattilaan sijoitettavan leijukerroslämmönsiirtimen säädön toteuttamista leijutusjärjestelyjen avulla. Teoreettinen tarkastelu keskittyy leijukerroslämmönsiirtimen toimintaolosuhteisiin ja kiintoainevirtauksen säätömahdollisuuksiin. Kokeellisessa osuudessa etsitään käytännössä toimivaa säätöperiaatetta tutkimuksen kohteena olevalle leijukerroslämmönsiirrinrakenteelle.

Jäännöshiilen ja kalkkikiven reaktiokinetiikan mallinnus kiertoleijukattilan tulipesässä

Työn tavoitteena oli kehittää mallit jäännöshiilen ja kalkkikiven reaktiokinetiikan sekä jäännöshiilen jauhautumisen ennustamiseksi kiertoleijukattilan tulipesässä. Kehitetyt mallit toimivat tulipesämallin osamalleina. Perustuen mallinnettuihin reaktionopeuksiin ja jauhautumiskäyttäytymiseen tulipesämalli ennustaa kalkkikivihiukkasten rikinsidonnan ja jäännöshiilen jakautumisen erikokoisiksi hiukkasiksi tulipesässä ja tuhkissa. Työssä kehitetyt mallit perustuvat olemassa oleviin kalkkikiven ja polttoaineen reaktiivisuustesteihin laboratorio-kokoluokan leijukerrosreaktorissa. Mallit huomioivat myös tulipesän olosuhteet. Menetelmät kelpoistettiin onnistuneesti kaupallisen kokoluokan kiertoleijukattiloista mitattujen ja tulipesämallilla laskettujen taseiden avulla. Mallien kehittämistä tullaan jatkamaan.