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ä.

Eetteröintiprosessien simulointivalmiuksien kehittäminen

Diplomityön tarkoituksena oli tutkia ETBE-prosessien simulointimallin kehittämistä. Simulointia varten valittiin eetteröinnin syöttövirroiksi etanoli ja tyypillinen FCC-kolonnin hiilivetysyöttö. Kirjallisuusosassa paneuduttiin tutkimaan syitä ETBE:n käytölle bensiinissä, valmistukseen tarvittavien raaka-aineiden lähteitä ja mahdollisia korvaavia raaka-aineita sekä mietittiin etanolissa olevien epäpuhtauksien vaikutusta prosessin tuotespesifikaatioihin. Kirjallisuusosassa tarkasteltiin lisäksi eri valmistajien markkinoimia eetteröintiteknologioita. Eetteröintiteknlogiat jaettiin perinteiseen tislaukseen pohjautuviin prosesseihin, reaktiiviseen tislaukseen pohjautuviin prosesseihin sekä Neste Engineering:in markkinoimaan NExETHERS-teknologiaan. Työn kokeellisessa osassa tutkittiin sekä etanolin epäpuhtauksien kulkeutumista prosessissa että kokeellisesti saadusta höyry-nestetasapainotiedosta määritettyjen Wilsonin yhtälön binääri-interaktioparametrien vaikutusta eetteröintiprosessin simulointituloksiin. Lopuksi simulointiin prosessia, jossa oli prosessin kannalta kriittiseksi havaitut etanolin epäpuhtaudet sekä hyviksi todetut binääri-interaktioparametrit etanolin ja FCC-syöttövirran C4-hiilivetyjen välillä. Uusilla binääriparametreilla saatuja simulointituloksia vertailtiin aikaisemmin samasta mallista vanhoilla binääriparametreilla saatuihin tuloksiin. Lopuksi tehtiin yhteenveto työn tuloksista ja annettiin ehdotukset jatkotutkimuksia varten.

Tapahtumapohjainen simulointi kuljetusketjun automatisoinnin kehittämisessä

Työn tavoitteena oli tutkia kuljetusketjun automatisoinnin kehittämistä tapahtumapohjaisen simuloinnin avulla. Työn teoria osassa käsitellään yleisellä tasolla simuloinnin teoriaa, siihen liittyviä käsitteitä ja erityisesti simulointiprojektin läpiviennin vaiheita. Tässä osassa saadaan vastaukset seuraaviin kysymyksiin: - Mitä tapahtumapohjaisella simuloinnilla tarkoitetaan? - Mitkä ovat simuloinnin hyödyt ja rajoitteet? - Mitä käyttökohteita simuloinnilla on? - Minkälaisia ohjelmia simulointiin käytetään? - Mitä vaiheita simulointiprojekti sisältää? Soveltavassa osassa tutkitaan kuinka tapahtumapohjaista simulointia hyödynnettiin Veto-Ketju -projektissa ja millaisia tuloksia rakennetulla simulointimallilla saavutettiin. Veto-Ketju -projekti käsittelee kuljetusketjun sekä varasto- ja satamakäsittelyn tehostamista automatisoinnin avulla. Projektin vetäjänä ja automaattisen tavarankäsittelyjärjestelmän toimittajana on Pesmel Oy ja satamatoimintojen ja logistiikan asiantuntijana SysOpen-konserniin kuuluva suunnittelu- ja konsulttitoimisto EP-Logistics Oy. Paperiteollisuudesta mukana ovat UPM-Kymmene ja M-real, satamaoperaattorina toimii Rauma Stevedoring ja kuljetusoperaattorina VR Cargo. Veto-Ketju –projektin simulointitutkimuksessa varmistettiin suunnitellun automaattisen junavaunujen purkausjärjestelmän toiminta ennen sen käyttöönottoa ja tutkittiin erilaisten toimintatapojen vaikutuksia satamatoimintoihin, satamassa tarvittavien resurssien määrään ja määritettiin tarvittavien varastojen koko. Simulointimallin avulla pystyttiin osoittamaan selkeästi erilaisten toimintavaihtoehtojen erot. Tällä tavoin saatiin tuotetuksi lisää tietoa päätöksenteon tueksi muun muassa järjestelmän sijoituspaikan ja mahdollisen uuden varaston rakentamisen suhteen.

Mekatronisten koneiden reaaliaikainen simulointi Linux-ympäristössä

Työn tavoitteena oli selvittää mahdollisuuksia käyttää Linux-ympäristöä mekatronisten koneiden reaaliaikaisessa simuloinnissa. Työssä tutkittiin C-kielellä tehdyn reaaliaikaisen simulointimallin ratkaisua Linux-käyttöjärjestelmässä RTLinux-reaaliaikalaajennuksen avulla. Reaaliaikainen simulointi onnistui RTLinuxin avulla tehokkaasti ja mallinnusmenetelmien rajoissa tarkasti. I/O-toimintojen lisäämistä erillisten I/O-korttien avulla ei tarkasteltu tässä työssä. Reaaliaikaista Linuxia ei ole aikaisemmin käytetty mekatronisten koneiden simulointiin. Tämän vuoksi valmiita työkaluja ei ole olemassa. Linux-ympäristö ei näin ollen sovellu kovin hyvin yleiseen koneensuunnitteluun mallintamisen työläyden vuoksi.

Tietoverkon varmistusten simulointi ja optimointi

Tässä diplomityössä määritellään varmistusjärjestelmän simulointimalli eli varmistusmalli. Varmistusjärjestelmän toiminta optimoidaan kyseisen varmistusmallin avulla. Optimoinnin tavoitteena on parantaa varmistusjärjestelmän tehokkuutta. Parannusta etsitään olemassa olevien varmistusjärjestelmän resurssien maksimaalisella hyödyntämisellä. Varmistusmalli optimoidaan evoluutioalgoritmin avulla. Optimoinnissa on useita tavoitteita, jotka ovat ristiriidassa keskenään. Monitavoiteoptimointiongelma muunnetaan yhden tavoitteen optimointiongelmaksi muodostamalla tavoitefunktio painotetun summan menetelmän avulla. Rinnakkain edellisen menetelmän kanssa käytetään myös Pareto-optimointia. Pareto-optimaalisen rintaman pisteiden etsintä ohjataan lähelle painotetun summan menetelmän optimipistettä. Evoluutioalgoritmin toteutuksessa käytetään hyväksi varmistusjärjestelmiin liittyvää ongelmakohtaista tietoa. Työn tuloksena saadaan varmistusjärjestelmän simulointi- sekä optimointityökalu. Simulointityökalua käytetään kartoittamaan nykyisen varmistusjärjestelmän toimivuutta. Optimoinnin avulla tehostetaan varmistusjärjestelmän toimintaa. Työkalua voidaan käyttää myös uusien varmistusjärjestelmien suunnittelussa sekä nykyisten varmistusjärjestelmien laajentamisessa.

Teollisuusverkon stabiilisuus saarekekäyttöön siirryttäessä

Työn aluksi on käsitelty sähkövoimajärjestelmien stabiilisuutta yleisesti. Myös teollisuuden sähköverkkojen erityispiirteitä on tarkasteltu. Teollisuusverkkojen komponentteja koskevassa selvityksessä eri laitteet on käsitelty siten, että niiden tärkeimmät ominaisuudet verkon erilaisissa muutostilanteissa tulevat esille. Tahtigeneraattorit ja epätahtikoneet on käsitelty yksityiskohtaisesti. Tutkimuksessa on esitelty teollisuusverkon saarekekäyttöön johtavia syitä, saarekkeeseen siirtymisen kriteereitä ja saarekekäytön toteutustapoja. Tehonvajaussuojaukseen liittyen on tutustuttu kuormanpudotusjärjestelyihin. Myös höyryverkon kuormien tasausten mahdollisuuksiin on luotu katsaus. Verkostolaskentaa käsittelevässä luvussa on keskitytty sopivan laskentamallin luomiseen liittyviin menetelmiin ja rajoituksiin. Lisäksi on kerrottu työn käytännön osan laskentamallin luomisen keskeiset periaatteet. Työn lopuksi on suoritettu simulointeja, joissa on tutkittu erään teollisuusverkon stabiilisuutta saarekekäyttöön siirryttäessä Calpos-verkostolaskentaohjelmistolla.

Virtuaaliprototyypin käyttö roottorin dynamiikan analysoinnissa

Työn tavoitteena oli selvittää kaupallisen dynamiikansimulointiohjelmiston soveltuvuus roottoridynamiikan analysointiin. Työssä keskityttiin erityisesti roottorin dynamiikkaan vaikuttavien epäideaalisuuksien mallintamiseen. Simulointitulosten tarkkuutta selvitettiin mittauksilla. Lisäksi vertailtiin yleiskäyttöisen dynamiikan simulointiohjelmiston ja roottoridynamiikan erikoisohjelmiston teoriaa. Tutkittava roottori oli paperikoneen putkitela. Telan joustavuus kuvattiin elementtimenetelmällä ratkaistujen moodien avulla. Elementtimallissa huomioitiin telan vaipan seinämänpaksuusvaihtelu, joka vaikuttaa telan massa- ja jäykkyysjakaumaan. Dynamiikkaohjelmistossa mallinnettiin telan tuennasta tulevat herätteet. Dynamiikkaohjelmistona käytettiin ADAMS:ia ja FEM-ohjelmana ANSYS:stä. Tuloksista havaittiin käytetyn menetelmän soveltuvan roottoridynamiikan analysointiin ja roottorin epäideaalisuuksien mallintamiseen. Simulointimallilla saatiin esille murtolukukriittiset pyörimisnopeudet ja telan kriittinen pyörimisnopeus vastasi hyvin mittaustuloksia.

Virtauslaskentaa ja energia-analyysiä hyödyntävä huoneilman lämpöolosuhteiden simulointimalli

Huonetilojen lämpöolosuhteiden hallinta on tärkeä osa talotekniikan suunnittelua. Tavallisesti huonetilan lämpöolosuhteita mallinnetaan menetelmillä, joissa lämpödynamiikkaa lasketaan huoneilmassa yhdessä laskentapisteessä ja rakenteissa seinäkohtaisesti. Tarkastelun kohteena on yleensä vain huoneilman lämpötila. Tämän diplomityön tavoitteena oli kehittää huoneilman lämpöolosuhteiden simulointimalli, jossa rakenteiden lämpödynamiikka lasketaan epästationaarisesti energia-analyysilaskennalla ja huoneilman virtauskenttä mallinnetaan valittuna ajanhetkenä stationaarisesti virtauslaskennalla. Tällöin virtauskentälle saadaan jakaumat suunnittelun kannalta olennaisista suureista, joita tyypillisesti ovat esimerkiksi ilman lämpötila ja nopeus. Simulointimallin laskentatuloksia verrattiin testihuonetiloissa tehtyihin mittauksiin. Tulokset osoittautuivat riittävän tarkoiksi talotekniikan suunnitteluun. Mallilla simuloitiin kaksi huonetilaa, joissa tarvittiin tavallista tarkempaa mallinnusta. Vertailulaskelmia tehtiin eri turbulenssimalleilla, diskretointitarkkuuksilla ja hilatiheyksillä. Simulointitulosten havainnollistamiseksi suunniteltiin asiakastuloste, jossa on esitetty suunnittelun kannalta olennaiset asiat. Simulointimallilla saatiin lisätietoa varsinkin lämpötilakerrostumista, joita tyypillisesti on arvioitu kokemukseen perustuen. Simulointimallin kehityksen taustana käsiteltiin rakennusten sisäilmastoa, lämpöolosuhteita ja laskentamenetelmiä sekä mallinnukseen soveltuvia kaupallisia ohjelmia. Simulointimallilla saadaan entistä tarkempaa ja yksityiskohtaisempaa tietoa lämpöolosuhteiden hallinnan suunnitteluun. Mallin käytön ongelmia ovat vielä virtauslaskennan suuri laskenta-aika, turbulenssin mallinnus, tuloilmalaitteiden reunaehtojen tarkka määritys ja laskennan konvergointi. Kehitetty simulointimalli tarjoaa hyvän perustan virtauslaskenta- ja energia-analyysiohjelmien kehittämiseksi ja yhdistämiseksi käyttäjäystävälliseksi talotekniikan suunnittelutyökaluksi.

Hylky- ja vesijärjestelmän vaikutus kartongin tuotepuhtauteen

Kirjallisuusosassa käsitellään paperi- ja kartonkikoneiden eri vesijärjestelmiä, vedenkäytön vähentämistä ja sen vaikutusta prosessiin. Mikrobiologiaa, mikrobien aiheuttamia prosessiongelmia, mikrobien torjuntaa ja simulointia on myös käsitelty kirjallisuusosassa. Kokeellisessa osassa laadittiin Kartonkikone 1:n massa- ja vesijärjestelmän mikrobiologisen puhtauden simulointimalli Balas-prosessisimulointiohjelmalla. Kartongin mikrobiologinen puhtaus määritettiin kokonaispesäkelukuna. Kokeellinen osa sisälsi pilot-mittakaavan suodatinkoeajoja, joissa määritettiin tarvittavat erotusparametrit simulointimallille, selvitettiin Dynasand-hiekkasuodattimen soveltuvuus suihkuvesien puhdistukseen, Dynadisc-kiekkosuodattimen käyttökelpoisuus kiertoveden puhdistukseen kiertovesisuodattimena ja hiekkasuodattimien esisuodattimena sekä UV-reaktorin soveltuvuus suodatettujen vesien desifiointiin. Simulointimallilla tarkasteltiin eri prosessivaihtoehtojen vaikutus kartongin mikrobiologiseen puhtauteen. Suihkuvesien puhdistaminen hiekkasuodattimilla alentaa kokonaispesäkelukua vain vähän. Koneen ajettavuus paranee suihkuputkien roskaryöppyjen ja viiraosan limoittumisen vähenemisen takia. Kiertovesijärjestelmän pelkkä jakaminen hylkyprosessi- ja sakeudensäätövedeksi aiheuttaa kokonaispesäkeluvun kasvun prosessista poistettavien vesien puhdistumisen takia. Kokonaispesäkeluku minimoidaan poistamalla ylimääräinen vesi hylkysaostajan suodoksena. Suodatettua hylkysaostajan suodoksen ja viiraveden seosta käytetään sekä sakeudensäätö- että hylkyprosessivedeksi. Hylkymassalla on dominoiva vaikutus kokonaispesäkelukuun. Suurin vähennys kokonaispesäkeluvussa saadaan estämällä mikrobikasvu hylkymassassa. Tämä voidaan tehdä kuumentamalla hylkymassa 80 °C:n höyryllä. Suodatuskoeajojen perusteella Dynasand-hiekkasuodattimella voidaan puhdistaa pisaralämmönvaihtimelta saatavaa lämmintä suihkuvettä. Esisuodatettua kiertovettä voidaan puhdistaa Dynasand-hiekkasuodattimella viiran johtotelojen suihkuvedeksi. Dynadisc-kiekkosuodatinta voidaan käyttää hiekkasuodattimen esisuodattimena ja kiertovesisuodattimena. UV-säteily poistaa tehokkaasti mikrobeja suodatetuista vesistä.

Liikesimulaattorialustan kehittäminen virtuaaliprototyypin avulla

Työn tarkoituksena oli kehittää kuuden vapausasteen liikesimulaattorialusta virtuaaliprototyypin avulla siten, että alustan dynamiikka mallinnetaan Adams-ohjelmistolla ja ohjaus- ja säätöpiiri Matlab Simulink:llä. Tarkoituksena oli tutkia dynaamisen mallin ja säätöteknisen mallin yhdistämistä ja niiden yhteen toimimista. Tarkoituksena oli myös selvittää tulevaisuudessa rakennettavan Stewart:n alustan mekaniikan mitat ja hydraulikomponenttien koot. Työssä tutkittiin alustan käyttäytymistä halutulla liikealueella, nopeuksia joita saavutetaan ja mekaanisia rajoitteita. Työn tuloksesta on tarkoitus rakentaa fyysinen prototyyppi liittyen KONSI-projektiin, jossa kehitetään satamanosturisimulaattori nosturiohjaajan koulutuksen tueksi. Malli tullaan kytkemään Teppo Lehtisen diplomityönä tehtyyn satamanosturi kontin simulointimalliin ja koko järjestelmän on tarkoitus toimia reaaliaikaisena.

Lajinvaihtokustannusten vaikutus tuotevalikoiman kannattavuuteen

Asiakkaat vaativat yhä suurempaa tuotannon joustavuutta ja lyhyitä toimitusaikoja paperiteollisuudelta. Samanaikaisesti kiristyvä kilpailu ajaa yritykset tavoittelemaan parempaa kannattavuutta kasvattamalla tuotevalikoimaa. Edellä mainitut tekijät lisäävät lajinvaihtojen määrää ja aiheuttavat kustannuksia. Diplomityön pääasiallisena tavoitteena on selvittää kohdeyrityksen paperikoneen lajinvaihdoista aiheutuvat kustannukset ja niiden vaikutus tuotevalikoiman kannattavuuteen. Toisena merkittävänä tavoitteena voidaan mainita lajinvaihtokustannusten simulointimallin rakentaminen, jota voidaan käyttää apuna tuotannonsuunnittelussa. Lisäksi saatujen tutkimustulosten perusteella pyritään löytämään keinoja kannattavuuden parantamiseksi. Paperikoneella tehdyn seurannan avulla selvitettiin eri tyyppisten lajinvaihtojen aiheuttamat tuotannonmenetykset. Tämän perusteella laskettiin niiden kustannusvaikutukset. Tarkastelua laajennettiin kahden toteutuneen syklin avulla syklitasolle ja lopuksi tuotevalikoimatasolle. Lajinvaihtokustannukset lisääntyivät paperin laatusuureiden muutoksen kasvaessa. Kalanterointimuutokset osoittautuivat kalleimmiksi ja niiden määrä tulisi minimoida. Kustannukset eivät osoittautuneet syklitasolla kovin suuriksi, mutta laji- ja ajokohtaisesti ne ovat merkittäviä erityisesti pienivolyymisilla tuotteilla. Niiden ajofrekvenssin harventaminen vähentäisi kustannuksia ja parantaisi kannattavuutta. Lajinvaihto-kustannukset tulisi jatkossa huomioida tuotekustannuslaskennassa.

TMP- laitoksen mallinnus ja simulointi

Diplomityössä tehtiin Balas- simulointimalli Stora Enso Publication Papers Oy Ltd Varkauden tehtaiden kuumahierrelaitoksesta, johon sovitettiin tehtaan lämpö-, kuitu- ja kiertovesivirtaukset sekä kiertovesissä olevien liuenneiden ja kolloidisten aineiden virtaukset. Kirjallisuusosassa perehdyttiin simulaattorimallin luontiin ja käyttötarkoituksiin. Siinä käsitellään Balas- simulaattorin ominaisuuksia ja laiteparametrointia. Tarkastellaan kuumahierreprosessin eri vaiheita, olosuhteita ja laitteiden toimintaa. Perehdytään prosessin energian kulutukseen ja talteenottoon sekä kiertovesien liuenneiden ja kolloidisten aineiden mittaussuureisiin ja vaikutuksiin prosessissa. Simulaattorimallin tekemiseen kuului tehtaan virtauskaavion ja simulaattorimallin tekeminen, mittauskoesuunnittelu, tehdasmittaukset, simulaattorin parametrointi, mittaus- ja simulaattoritulosten yhteensovittaminen sekä mallin kelpoistaminen. Tehtaan virtauskaavion piirtämisessä käytettiin apuna tehtaan prosessi- ja instrumentointikaavioita, joiden pohjalta mittauskoesuunnitelma ja simulaattorimalli tehtiin. Koesuunnitelman mittaussuureiksi valittiin virtaukset, sakeudet, lämpötilat sekä kiertovesien kiintoaineen-, kokonaisorgaanisen hiilen- ja liuenneen orgaanisen aineen pitoisuuksien määritykset. Mittauksilla saatiin tietoja prosessivirtausten taseista, joita käytettiin simulaattorimallin keskeisempien laiteparametriarvojen määrityksessä. Exceliä hyödynnettiin mittaus- ja simulaattoritulosten taulukoinnissa, laiteparametrien laskemisessa sekä arvojen syötössä ja vastaanotossa Excelin ja simulaattorin välillä. Sitä käytettiin myös graafisessa mittaus- ja simulaattoritulosten yhteensovittamisessa, jolla pystyttiin havainnollisesti näkemään eri syöttöparametrien muutoksien vaikutukset simulaattorin laskemissa virtaustuloksissa. Mallin antamien arvojen ja todellisten prosessimittausarvojen yhteensovittamisen ja mallista varmistumisen tuloksista voidaan todeta mallin korreloivan todellista prosessia melko hyvin.

Sähkökäyttöisen konejärjestelmän vuorovaikutteinen simulointi

Työn tavoitteena oli toteuttaa simulointimalli, jolla pystytään tutkimaan kestomagnetoidun tahtikoneen aiheuttaman vääntömomenttivärähtelyn vaikutuksia sähkömoottoriin liitetyssä mekaniikassa. Tarkoitus oli lisäksi selvittää kuinka kyseinen simulointimalli voidaan toteuttaa nykyaikaisia simulointiohjelmia käyttäen. Saatujen simulointitulosten oikeellisuus varmistettiin tätä työtä varten rakennetulla verifiointilaitteistolla. Tutkittava rakenne koostui akselista, johon kiinnitettiin epäkeskotanko. Epäkeskotankoon kiinnitettiin massa, jonka sijaintia voitiin muunnella. Massan asemaa muuttamalla saatiin rakenteelle erilaisia ominaistaajuuksia. Epäkeskotanko mallinnettiin joustavana elementtimenetelmää apuna käyttäen. Mekaniikka mallinnettiin dynamiikan simulointiin tarkoitetussa ADAMS –ohjelmistossa, johon joustavana mallinnettu epäkeskotanko tuotiin ANSYS –elementtimenetelmäohjelmasta. Mekaniikan malli siirrettiin SIMULINK –ohjelmistoon, jossa mallinnettiin myös sähkökäyttö. SIMULINK –ohjelmassa mallinnettiin sähkökäyttö, joka kuvaa kestomagnetoitua tahtikonetta. Kestomagnetoidun tahtikoneen yhtälöt perustuvat lineaarisiin differentiaaliyhtälöihin, joihin hammasvääntömomentin vaikutus on lisätty häiriösignaalina. Sähkökäytön malli tuottaa vääntömomenttia, joka syötetään ADAMS –ohjelmistolla mallinnettuun mekaniikkaan. Mekaniikan mallista otetaan roottorin kulmakiihtyvyyden arvo takaisinkytkentänä sähkömoottorin malliin. Näin saadaan aikaiseksi yhdistetty simulointi, joka koostuu sähkötoimilaitekäytöstä ja mekaniikasta. Tulosten perusteella voidaan todeta, että sähkökäyttöjen ja mekaniikan yhdistetty simulointi on mahdollista toteuttaa valituilla menetelmillä. Simuloimalla saadut tulokset vastaavat hyvin mitattuja tuloksia.

CITRICULTURE ECONOMIC AND FINANCIAL EVALUATION UNDER CONDITIONS OF UNCERTAINTY

ABSTRACT The citriculture consists in several environmental risks, as weather changes and pests, and also consists in considerable financial risk, mainly due to the period ofreturn on the initial investment. This study was motivated by the need to assess the risks of a business activity such as citriculture. Our objective was to build a stochastic simulation model to achieve the economic and financial analysis of an orange producer in the Midwest region of the state of Sao Paulo, under conditions of uncertainty. The parameters used were the Net Present Value (NPV), the Modified Internal Rate of Return(MIRR), and the Discounted Payback. To evaluate the risk conditions we built a probabilistic model of pseudorandom numbers generated with Monte Carlo method. The results showed that the activity analyzed provides a risk of 42.8% to reach a NPV negative; however, the yield assessed by MIRR was 7.7%, higher than the yield from the reapplication of the positive cash flows. The financial investment pays itself after the fourteenth year of activity.

Sinkkirikasteen liuotusprosessin kineettinen mallinnus

Työssä on tehty kineettinen simulointimalli sinkkirikasteen liuotusprosessista. Prosessi on pieni osa Kokkolan sinkinvalmistusprosessia, jonka muita osia ovat: pasutus, neutraaliliuotus, konversio, liuospuhdistus ja elektrolyysi. Rikasteen liuotukseen tulee konversioprosessin liuos ja liuotuksesta lähtevä neste menee takaisin neutraaliliuotukseen. Saostunut jarosiitti läjitetään. Kokkolan liuotusprosessi koostuu liettoreaktorista ja kahdesta neljän liuotusreaktorin sarjasta. Liuotukseen syötetään paluuhappoa liettoreaktoriin ja liuotuspiirien ensimmäisiin liuotusreaktoreihin. Happea syötetään kaikkiin liuotusreaktoreihin. Prosessin mallintamiseen käytettiin Aspen Plus-simulointiohjelmaa, johon pystyttiin syöttämään kineettisiä yhtälöitä. Reaktionopeusyhtälöitä käytettiin raudan hapetuksen, sulfidien liuotuksen ja jarosiitiin saostumisen mallintamiseen, eli kaikkiin liuotusreaktoreissa tapahtuviin reaktioihin. Kineettiset yhtälöt etsittiin kirjallisuudesta. Liettoreaktori puolestaan mallinnettiin syöttämällä ohjelmaan reaktioyhtälöt ja antamalla niille etenemisasteet. Jarosiitin liukenemisesta työssä on tehty laboratoriokokeita, koska aiheesta ei kirjallisuudesta löytynyt kineettistä tietoa. Liuotuskokeissa käytetyn kiintoaineen kuitenkin todettiin sisältävän liikaa götiittiä, että tuloksista olisi voitu laskea kinetiikkaa jarosiitin liukenemiselle. Simulointimallilla laskettiin yksi tapaus vertailukohdaksi, johon malliin tehtyjä muutoksia verrattiin. Mallilla tutkittiin konversiosta tulevan jarosiitin määrän vaikutusta, reaktorikoon merkitystä ja rikasteen liuotuksen sekä jarosiitin saostuksen reaktionopeuksien muutoksen vaikutuksia. Käytetyillä kineettisillä yhtälöillä reaktioiden todettiin tarvitsevan vain ¾ käytetystä reaktiotilavuudesta, rikasteen liuotusnopeuden kohtalaisen pienellä hidastamisen todettiin vähentävän sinkin saantoa ja jarosiitin saostuksen reaktionopeuden kasvulla todettiin myös olevan negatiivinen vaikutus sinkin saantoon. Simulointimallissa käytettyjen reaktionopeusyhtälöiden varmentaminen kokeilla todettiin tarpeelliseksi, sillä jo kohtalaisen pienillä muutoksilla havaittiin olevan merkitystä prosessin toimivuuteen. Lisäksi todettiin jarosiitin liukenemisen huomioimisen olevan tarpeen.