Comparison of Model-Based Flow Rate Estimation Methods in Frequency-Converter-Driven Pumps and Fans

Fluid handling systems such as pump and fan systems are found to have a significant potential for energy efficiency improvements. To deliver the energy saving potential, there is a need for easily implementable methods to monitor the system output. This is because information is needed to identify inefficient operation of the fluid handling system and to control the output of the pumping system according to process needs. Model-based pump or fan monitoring methods implemented in variable speed drives have proven to be able to give information on the system output without additional metering; however, the current model-based methods may not be usable or sufficiently accurate in the whole operation range of the fluid handling device. To apply model-based system monitoring in a wider selection of systems and to improve the accuracy of the monitoring, this paper proposes a new method for pump and fan output monitoring with variable-speed drives. The method uses a combination of already known operating point estimation methods. Laboratory measurements are used to verify the benefits and applicability of the improved estimation method, and the new method is compared with five previously introduced model-based estimation methods. According to the laboratory measurements, the new estimation method is the most accurate and reliable of the model-based estimation methods.

Computational Analysis and Optimization of Real Gas Flow In Small Centrifugal Compressors

Small centrifugal compressors are more and more widely used in many industrialsystems because of their higher efficiency and better off-design performance comparing to piston and scroll compressors as while as higher work coefficient perstage than in axial compressors. Higher efficiency is always the aim of the designer of compressors. In the present work, the influence of four partsof a small centrifugal compressor that compresses heavy molecular weight real gas has been investigated in order to achieve higher efficiency. Two parts concern the impeller: tip clearance and the circumferential position of the splitter blade. The other two parts concern the diffuser: the pinch shape and vane shape. Computational fluid dynamics is applied in this study. The Reynolds averaged Navier-Stokes flow solver Finflo is used. The quasi-steady approach is utilized. Chien's k-e turbulence model is used to model the turbulence. A new practical real gas model is presented in this study. The real gas model is easily generated, accuracy controllable and fairly fast. The numerical results and measurements show good agreement. The influence of tip clearance on the performance of a small compressor is obvious. The pressure ratio and efficiency are decreased as the size of tip clearance is increased, while the total enthalpy rise keeps almost constant. The decrement of the pressure ratio and efficiency is larger at higher mass flow rates and smaller at lower mass flow rates. The flow angles at the inlet and outlet of the impeller are increased as the size of tip clearance is increased. The results of the detailed flow field show that leakingflow is the main reason for the performance drop. The secondary flow region becomes larger as the size of tip clearance is increased and the area of the main flow is compressed. The flow uniformity is then decreased. A detailed study shows that the leaking flow rate is higher near the exit of the impeller than that near the inlet of the impeller. Based on this phenomenon, a new partiallyshrouded impeller is used. The impeller is shrouded near the exit of the impeller. The results show that the flow field near the exit of the impeller is greatly changed by the partially shrouded impeller, and better performance is achievedthan with the unshrouded impeller. The loading distribution on the impeller blade and the flow fields in the impeller is changed by moving the splitter of the impeller in circumferential direction. Moving the splitter slightly to the suction side of the long blade can improve the performance of the compressor. The total enthalpy rise is reduced if only the leading edge of the splitter ismoved to the suction side of the long blade. The performance of the compressor is decreased if the blade is bended from the radius direction at the leading edge of the splitter. The total pressure rise and the enthalpy rise of thecompressor are increased if pinch is used at the diffuser inlet. Among the fivedifferent pinch shape configurations, at design and lower mass flow rates the efficiency of a straight line pinch is the highest, while at higher mass flow rate, the efficiency of a concave pinch is the highest. The sharp corner of the pinch is the main reason for the decrease of efficiency and should be avoided. The variation of the flow angles entering the diffuser in spanwise direction is decreased if pinch is applied. A three-dimensional low solidity twisted vaned diffuser is designed to match the flow angles entering the diffuser. The numerical results show that the pressure recovery in the twisted diffuser is higher than in a conventional low solidity vaned diffuser, which also leads to higher efficiency of the twisted diffuser. Investigation of the detailed flow fields shows that the separation at lower mass flow rate in the twisted diffuser is later than in the conventional low solidity vaned diffuser, which leads to a possible wider flow range of the twisted diffuser.

Numerical and ExperimentalModelling of Gas Flow and Heat Transfer in the Air Gap of An Electric Machine

This work deals with the cooling of high-speed electric machines, such as motors and generators, through an air gap. It consists of numerical and experimental modelling of gas flow and heat transfer in an annular channel. Velocity and temperature profiles are modelled in the air gap of a high-speed testmachine. Local and mean heat transfer coefficients and total friction coefficients are attained for a smooth rotor-stator combination at a large velocity range. The aim is to solve the heat transfer numerically and experimentally. The FINFLO software, developed at Helsinki University of Technology, has been used in the flow solution, and the commercial IGG and Field view programs for the grid generation and post processing. The annular channel is discretized as a sector mesh. Calculation is performed with constant mass flow rate on six rotational speeds. The effect of turbulence is calculated using three turbulence models. The friction coefficient and velocity factor are attained via total friction power. The first part of experimental section consists of finding the proper sensors and calibrating them in a straight pipe. After preliminary tests, a RdF-sensor is glued on the walls of stator and rotor surfaces. Telemetry is needed to be able to measure the heat transfer coefficients at the rotor. The mean heat transfer coefficients are measured in a test machine on four cooling air mass flow rates at a wide Couette Reynolds number range. The calculated values concerning the friction and heat transfer coefficients are compared with measured and semi-empirical data. Heat is transferred from the hotter stator and rotor surfaces to the coolerair flow in the air gap, not from the rotor to the stator via the air gap, althought the stator temperature is lower than the rotor temperature. The calculatedfriction coefficients fits well with the semi-empirical equations and precedingmeasurements. On constant mass flow rate the rotor heat transfer coefficient attains a saturation point at a higher rotational speed, while the heat transfer coefficient of the stator grows uniformly. The magnitudes of the heat transfer coefficients are almost constant with different turbulence models. The calibrationof sensors in a straight pipe is only an advisory step in the selection process. Telemetry is tested in the pipe conditions and compared to the same measurements with a plain sensor. The magnitudes of the measured data and the data from the semi-empirical equation are higher for the heat transfer coefficients than thenumerical data considered on the velocity range. Friction and heat transfer coefficients are presented in a large velocity range in the report. The goals are reached acceptably using numerical and experimental research. The next challenge is to achieve results for grooved stator-rotor combinations. The work contains also results for an air gap with a grooved stator with 36 slots. The velocity field by the numerical method does not match in every respect the estimated flow mode. The absence of secondary Taylor vortices is evident when using time averagednumerical simulation.

Process intensification: From optimised flow patterns to microprocess technology

This thesis is focused on process intensification. Several significant problems and applications of this theme are covered. Process intensification is nowadays one of the most popular trends in chemical engineering and attempts have been made to develop a general, systematic methodology for intensification. This seems, however, to be very difficult, because intensified processes are often based on creativity and novel ideas. Monolith reactors and microreactors are successful examples of process intensification. They are usually multichannel devices in which a proper feed technique is important for creating even fluid distribution into the channels. Two different feed techniques were tested for monoliths. In the first technique a shower method was implemented by means of perforated plates. The second technique was a dispersion method using static mixers. Both techniques offered stable operation and uniform fluid distribution. The dispersion method enabled a wider operational range in terms of liquid superficial velocity. Using dispersion method, a volumetric gas-liquid mass transfer coefficient of 2 s-1 was reached. Flow patterns play a significant role in terms of the mixing performance of micromixers. Although the geometry of a T-mixer is simple, channel configurations and dimensions had a clear effect on mixing efficiency. The flow in the microchannel was laminar, but the formation of vortices promoted mixing in micro T-mixers. The generation of vortices was dependent on the channel dimensions, configurations and flow rate. Microreactors offer a high ratio of surface area to volume. Surface forces and interactions between fluids and surfaces are, therefore, often dominant factors. In certain cases, the interactions can be effectively utilised. Different wetting properties of solid materials (PTFE and stainless steel) were applied in the separation of immiscible liquid phases. A micro-scale plate coalescer with hydrophilic and hydrophobic surfaces was used for the continuous separation of organic and aqueous phases. Complete phase separation occurred in less than 20 seconds, whereas the separation time by settling exceeded 30 min. Fluid flows can be also intensified in suitable conditions. By adding certain additives into turbulent fluid flow, it was possible to reduce friction (drag) by 40 %. Drag reduction decreases frictional pressure drop in pipelines which leads to remarkable energy savings and decreases the size or number of pumping facilities required, e.g., in oil transport pipes. Process intensification enables operation often under more optimal conditions. The consequent cost savings from reduced use of raw materials and reduced waste lead to greater economic benefits in processing.

Flow phenomena, heat and mass transfer in microchannel reactors

The studies of flow phenomena, heat and mass transfer in microchannel reactors are beneficial to estimate and evaluate the ability of microchannel reactors to be operated for a given process reaction such as Fischer-Tropsch synthesis. The flow phenomena, for example, the flow regimes and flow patterns in microchannel reactors for both single phase and multiphase flow are affected by the configuration of the flow channel. The reviews of the previous works about the analysis of related parameters that affect the flow phenomena are shown in this report. In order to predict the phenomena of Fischer-Tropsch synthesis in microchannel reactors, the 3-dimensional computational fluid dynamic simulation with commercial software package FLUENT was done to study the flow phenomena and heat transfer for gas phase Fischer-Tropsch products flow in rectangular microchannel with hydraulic diameter 500 ¿m and length 15 cm. Numerical solution with slip boundary condition was used in the simulation and the flowphenomena and heat transfer were determined.

Synteesikaasun liukoisuus ja aineensiirto fermentointiliuokseen

Työssä tutkittiin synteesikaasun komponenttien: hiilimonoksidin, vedyn ja hiilidioksidin liukoisuutta ja aineensiirtonopeutta fermentointiliuokseen. Kirjallisuusosassa käsitellään kaasujen liukoisuuksiin ja kaasu-nesteaineensiirtoon vaikuttavia tekijöitä ja esitellään korrelaatioita, jotka on kehitetty volumetrisen aineensiirtokertoimen ennustamiseen sekoitetussa fermentorissa. Kirjallisuus-osassa esitetään myös synteesikaasun komponenttien liukoisuudet veteen, etanoliin ja etikkahappoon 37 ºC lämpötilassa ja esitellään Flowbat-simulointiohjelman MHV2-mallin käyttöä kaasu-neste tasapainojen mallin-nuksessa. Työn kokeellisessa osassa tutkittiin synteesikaasun komponenttien liukoisuuksia veteen ja kasvatusalustaan sekä kehitettiin mittausmenetelmä kaasu-neste tasapainojen mittaukseen. Tasapainomittauksissa tutkittiin etanolin ja etikkahapon konsentraatioiden vaikutusta synteesikaasun liukoisuuteen. Lisäksi mallinnettiin kaasu-neste tasapainoja monikomponenttisysteemeissä MHV2-mallin avulla. Kokeellisen osan aineensiirtomittauksissa tutkittiin sekoitusnopeuden ja kaasun volumetrisen syöttönopeuden vaikutusta hiilimonoksidin ja vedyn volumetriseen aineensiirtokertoimeen kLa kahden litran tilavuuksisessa laboratoriofermentorissa. Mittaustulosten perusteella kasvatusalustan komponentit vaikuttavat merkittävästi hiilidioksidin liukoisuuteen. Lisäys etanolin ja etikkahapon konsentraatioissa parantaa hiilimonoksidin liukoisuutta kasvatusalustaan. Kaasun volumetrinen syöttönopeus ja sekoitusnopeus vaikuttavat voimakkaasti volumetrisen aineensiirtokertoimen arvoon. Tutkitussa systeemissä korkein teoreettinen solutiheys, joka voitiin saavuttaa suurimmalla hiilimonoksidin aineensiirto-nopeudella, oli 3 g/L. Tämä on kaksinkertainen verrattuna aiemmissa VTT:n kokeissa saavutettuihin solutiheyksiin.

Stock flow through a drum screen

Tässä diplomityössä haluttiin mallintaa kuidutusrummun toimintaa Fluent virtausmallinusohjelman avulla. Aikaisempi tieto ja kehitystyö on perustunut kokemukseen ja käytännön kokeisiin. Kehitystyön alkuaikoina on suoritettu muutamia laskelmia koskien rummun tuottoa mutta sen jälkeen ei toimintaa ole laskennallisesti kuvattu. Työn ensimmäinen osa käsittelee yleisesti keräyspaperin käsittelyyn liittyviä laitteita ja menetelmiä. Toimintaperiaatteita on kuvattu yleisellä tasolla ja FibreFlow® rumpu on sitten käsitelty muita laitteita tarkemmin. Työn toinen osa sisältää sitten laboratoriotestit paikallisilta tahteilta hankittujen näytteiden viskositeettien ja tiheyksien määrittämiseksi. Kokeet suoritettiin Kemiantekniikan osastolla Brookfield viskoosimetrillä. Joitain alustavia laskentoja tuotosta suoritettiin aikaisempien tietojen perusteella. Rumpua kun on valmistettu vuodesta 1976, on tietoa kertynyt runsaasti vuosien mittaan. Laskelmia varten valittiin mallinnettavaksi alueeksi vain yksittäinen reikä sihdistä jolle laskettiin massavirta. Käytetyt laadut olivat OCC ja DIP. Myös eri rumpukoot otettiin jossain määrin huomioon.

Coolant Flow and Heat Transfer in Pebble-Bed Reactor Core

This thesis gathers knowledge about ongoing high-temperature reactor projects around the world. Methods for calculating coolant flow and heat transfer inside a pebble-bed reactor core are also developed. The thesis begins with the introduction of high-temperature reactors including the current state of the technology. Process heat applications that could use the heat from a high-temperature reactor are also introduced. A suitable reactor design with data available in literature is selected for the calculation part of the thesis. Commercial computational fluid dynamics software Fluent is used for the calculations. The pebble-bed is approximated as a packed-bed, which causes sink terms to the momentum equations of the gas flowing through it. A position dependent value is used for the packing fraction. Two different models are used to calculate heat transfer. First a local thermal equilibrium is assumed between the gas and solid phases and a single energy equation is used. In the second approach, separate energy equations are used for the phases. Information about steady state flow behavior, pressure loss, and temperature distribution in the core is obtained as results of the calculations. The effect of inlet mass flow rate to pressure loss is also investigated. Data found in literature and the results correspond each other quite well, considered the amount of simplifications in the calculations. The models developed in this thesis can be used to solve coolant flow and heat transfer in a pebble-bed reactor, although additional development and model validation is needed for better accuracy and reliability.

Kromatografinen erotus bioetanolin valmistuksessa lignoselluloosasta

Bioetanolin tuotanto kiinnostaa monissa maissa johtuen kansainvälisissä sopimuksissa määritellyistä ilmastotavoitteista. Työssä tutkittiin laboratorio-oloissa ioninvaihtohartsien ominaisuuksien ja erotuksen olosuhteiden vaikutusta rikkihapon ja glukoosin kromatografiseen erotukseen. Tehokkaimmaksi hartsiksi osoittautui polysulfonoitu mesohuokoinen vahva kationinvaihtohartsi Finex CS100C. CS100C:lla voitiin erottaa rikkihappoa ja glukoosia tehokkaimmin korkeissa 25 p-% ja 36 p-% glukoosi- ja rikkihappo-pitoisuuksissa. Lisäksi sillä havaittiin suurin tuotto simuloidussa liikkuvassa pedissä. Yhdessä kolonnissa suoritetuissa erotuskokeissa tutkittiin hartsien erotuskykyä rikkihapolle ja glukoosille sekä virtausnopeuden vaikutusta erotukseen lämpötilassa 22 °C. Saatujen tulosten pohjalta valittiin CS11GC, CS16GC ja CS100C tarkempaan isotermin määritykseen ja simulointiin hyvän erotuskyvyn sekä keskinäisten erojen takia. Adsorptioisotermit määritettiin kolonnikokein sekä 22 °C:n että 50 °C:n lämpötilassa. Isotermeistä havaittiin, että tasapaino kiinto- ja liuosfaasien välille saavutetaan rikkihapolla alhaisella 1 cm3/min virtausnopeudella varmemmin kuin suuremmalla 2,5 cm3/min virtausnopeudella. 50 °C:n lämpötilassa hapon ja glukoosin isotermit olivat jyrkempiä kuin 22 °C:n lämpötilassa. Määritettyihin hapon ja sokerin isotermeihin sovitettiin mallit, joiden parametreja käytettiin yksittäisen kolonnin simulointiin. Simuloinnissa oli estimoitavia parametreja yhdellä kolonnilla aineensiirtokertoimet sekä läpäisykäyristä määritetyt isotermiparametrit glukoosille sekä rikkihapolle ja SMB–erotuksessa vyöhykkeiden 2 ja 3 suhteelliset virtausnopeudet. Siirryttäessä lämpötilojen 22 °C ja 50 °C välillä hartsien parametrit muuttuivat sokerille täysin ja hapolle vain aineensiirtokertoimen osalta. CS100C oli tehokkain SMB–erotuksessa korkeimmalla 0,11 cm3/min tuottavuudella 95 %:n saannon saavuttamiseksi 95 % tuotepuhtaudella raffinaatissa ja ekstraktissa.

Influence of reactant absorption and dissolution in heterogeneous precipitation processes

In this research work, the aim was to investigate the volumetric mass transfer coefficient [kLa] of oxygen in stirred tank in the presence of solid particle experimentally. The kLa correlations as a function of propeller rotation speed and flow rate of gas feed were studied. The O2 and CO2 absorption in water and in solid-liquid suspensions and heterogeneous precipitation of MgCO3 were thoroughly examined. The absorption experiments of oxygen were conducted in various systems like pure water and in aqueous suspensions of quartz and calcium carbonate particles. Secondly, the precipitation kinetics of magnesium carbonate was also investigated. The experiments were performed to study the reactive crystallization with magnesium hydroxide slurry and carbon dioxide gas by varying the feed rates of carbon dioxide and rotation speeds of mixer. The results of absorption and precipitation are evaluated by titration, total carbon (TC analysis), and ionic chromatrography (IC). For calcium carbonate, the particle concentration was varied from 17.4 g to 2382 g with two size fractions: 5 µm and 45-63 µm sieves. The kLa and P/V values of 17.4 g CaCO3 with particle size of 5µm and 45-63 µm were 0.016 s-1 and 2400 W/m3. At 69.9 g concentration of CaCO3, the achieved kLa is 0.014 s-1 with particle size of 5 µm and 0.017 s-1 with particle size of 45 to 63 µm. Further increase in concentration of calcium carbonate, i.e. 870g and 2382g , does not affect volumetric mass transfer coeffienct of oxygen. It could be concluded from absorption results that maximum value of kLa is 0.016 s-1. Also particle size and concentration does affect the transfer rate to some extend. For precipitation experiments, the constant concentration of Mg(OH)2 was 100 g and the rotation speed varied from 560 to 750 rpm, whereas the used feed rates of CO2 were 1 and 9 L/min. At 560 rpm and feed rate of CO2 is 1 L/min, the maximum value of Mg ion and TC were 0.25 mol/litre and 0.12 mol/litre with the residence time of 40 min. When flow rate of CO2 increased to 9 L/min with same 560 rpm, the achieved value of Mg and TC were 0.3 mol/litre and 0.12 mol/L with shorter residence time of 30 min. It is concluded that feed rate of CO2 is dominant in precipitation experiments and it has a key role in dissociation and reaction of magnesium hydroxide in precipitation of magnesium carbonate.

Flow in porous media of oil-water systems

Fluid flow behaviour in porous media is a conundrum. Therefore, this research is focused on filtration-volumetric characterisation of fractured-carbonate sediments, coupled with their proper simulation. For this reason, at laboratory rock properties such as pore volume, permeability and porosity are measured, later phase permeabilities and oil recovery in function of flow rate are assessed. Furthermore, the rheological properties of three oils are measured and analysed. Finally based on rock and fluid properties, a model using COMSOL Multiphysics is built in order to compare the experimental and simulated results. The rock analyses show linear relation between flow rate and differential pressure, from which phase permeabilities and pressure gradient are determined, eventually the oil recovery under low and high flow rate is established. In addition, the oils reveal thixotropic properties as well as non-Newtonian behaviour described by Bingham model, consequently Carreau viscosity model for the used oil is given. Given these points, the model for oil and water is built in COMSOL Multiphysics, whereupon successfully the reciprocity between experimental and simulated results is analysed and compared. Finally, a two-phase displacement model is elaborated.

Öljynjalostamon tyhjöyksikön esilämmitysuunin optimointi

Työn tarkoituksena oli vertailla kahta kaupallista prosessiuunisimulointiohjelmaa sekä selvittää onko toinen ohjelma soveltuva Neste Jacobs Oy:n käyttöön. Tavoitteena oli selvittää ohjelmien ominaisuuksien eroja ja rajoitteita sekä tuloksien yhdenmukaisuutta. Kirjallisuusosassa perehdyttiin prosessiuunien toimintaperiaatteisiin, yleisiin mitoitusperiaatteisiin ja säteilylämmönsiirron perusteisiin. Kirjallisuusosassa käsiteltiin myös prosessiuunien polttimien toimintaa ja NOx-päästöjä vähentäviä polttotekniikoita. Molemmilla laskentaohjelmilla laadittiin kuumaöljy- ja tyhjöuuni, joiden tuloksia vertailtiin tärkeimpien prosessiarvojen perusteella ja toisen ohjelman osalta myös laitevalmistajan ilmoittamiin arvoihin. Tällä ohjelmalla tutkittiin tyhjöuunin höyrynruiskutuksen ja koksaantumisen vaikutuksia prosessiin sekä säteilyosan leveyden muutoksen vaikutusta säteilylämmönsiirtoon. Ohjelmien ominaisuuksien keskeisimmät erot liittyivät kaksifaasivirtauksen laskentaan, prosessiaineen massavirran rajoituksiin, ilman esilämmittimen käyttöön ja laskentamallien paineen kontrollointiin. Tulokset erosivat savukaasun ja tuubien lämpötilojen ja painehäviön laskennassa. Tuloksien ja ominaisuuksien perusteella toisen ohjelman ei todettu olevan täysin soveltuva Neste Jacobs Oy:n käyttöön.

Aerodynamical instability of web

Diplomityön tavoitteena oli tutkia miten ilman turbulenttisuus vaikuttaa tasaisesti liikkuvan rainan tilaan. Yhtenä sovelluskohteena teollisuudessa voidaan mainita esimerkiksi leiju-kuivain. Tiedetään, että konenopeuksien kasvu ja siitä johtuva ilmavirran nopeuden kasvu aiheuttaa voimavaikutuksia rainaan ja voi aiheuttaa lepatusta. Lepatus johtaa dynaamiseen epästabiilisuuteen, joka voidaan havaita, kun lineaarinen systeemi tulee epävakaaksi ja joh-taa epälineaariseen, rajoitettuun värähtelyyn. Lepatus huonontaa tuotteiden laatua ja voi johtaa ratakatkoihin. Työssä on esitetty tietoa ilman ja rainan vuorovaikutuksesta, jota hyödyntämällä voidaan kehittää yksinkertaistettu malli, jonka avulla liikkuvaa rainaa voidaan simuloida kuivaimes-sa. Kaasufaasin virtausyhtälöt on ratkaistu eri turbulenttimalleja käyttäen. Myös viskoelas-tisen rainan muodonmuutosta on tarkasteltu. Koska rainalle ei ole kirjallisuudesta saatavilla tarkkoja fysikaalisia ja mekaanisia arvoja, näitä ominaisuuksia testattiin eri arvoilla, jotta rainan käyttäytymistä jännityksen alaisena voidaan tarkastella. Näiden ominaisuuksien tun-teminen on ensiarvoisen tärkeää määritettäessä rainan aeroviskoelastista käyttäytymistä. Virtaussimulointi on kallista ja aikaa vievää. Tämä tarkoittaa uusien tutkimusmenetelmien omaksumista. Tässä työssä vaihtoehtoisena lähestymistapana on esitetty yksinkertaistettu malli, joka sisältää ilman ja rainan vuorovaikutusta kuvaavat ominaisuudet. Mallin avulla saadaan tietoa epälineaarisuuden ja turbulenssin vaikutuksesta sekä monimutkaisesta yh-teydestä stabiilisuuden ja ulkoisesti aikaansaadun värähtelyn sekä itse aiheutetun värähtelyn välillä. Työn lopussa on esitetty havainnollinen esimerkki, jolla voidaan kuvata olosuhteita, jossa rainan tasainen liike muuttuu epävakaaksi. Kun turbulenttisuudesta johtuva painevaih-telu ylittää tietyn rajan, rainan värähtely kasvaa muuttuen satunnaisesta järjestäytyneeksi. Saaduttulokset osoittavat, että turbulenttisuudella on suuri vaikutus eikä sitä voi jättää huomioimatta. Myös rainan viskoelastiset ominaisuudet tulee huomioida, jotta rainan käyt-täytymistä voidaan kuvata tarkasti.

Ilmavirtausten laskennallinen mallinnus

Diplomityössä tutkitaan kolmea erilaista virtausongelmaa CFD-mallinnuksella. Yhteistä näille ongelmille on virtaavana aineena oleva ilma. Lisäksi tapausten perinteinen mittaus on erittäin vaikeaa tai mahdotonta. Ensimmäinen tutkimusongelma on tarrapaperirainan kuivain, jonka tuotantomäärä halutaan nostaa kaksinkertaiseksi. Tämä vaatii kuivatustehon kaksinkertaistamista, koska rainan viipymäaika kuivausalueella puolittuu. Laskentayhtälöillä ja CFD-mallinnuksella tutkitaan puhallussuihkun nopeuden ja lämpötilan muutoksien vaikutusta rainan pinnan lämmön- ja massansiirtokertoimiin. Tuloksena saadaan varioitujen suureiden sekä massan- ja lämmönsiirtokertoimien välille riippuvuuskäyrät, joiden perusteella kuivain voidaan säätää parhaallamahdollisella tavalla. Toinen ongelma käsittelee suunnitteilla olevan kuparikonvertterin sekundaarihuuvan sieppausasteen optimointia. Ilman parannustoimenpiteitä käännetyn konvertterin päästöistä suurin osa karkaa ohi sekundaarihuuvan. Tilannetta tutkitaan konvertterissa syntyvän konvektiivisen nostevirtauksen eli päästöpluumin sekä erilaisten puhallussuihkuratkaisujen CFD-mallinnuksella. Tuloksena saadaan puhallussuihkuilla päästöpluumia poikkeuttava ilmaverho. Suurin osa nousevasta päästöpluumista indusoituu ilmaverhoon ja kulkeutuu poistokanavaan. Kolmas tutkittava kohde on suunnitteilla oleva kuparielektrolyysihalli, jossa ilmanvaihtoperiaatteena on luonnollinen ilmanvaihto ja mekaaninen happosumun keräysjärjestelmä. Ilmanvaihtosysteemin tehokkuus ja sisäilman virtaukset halutaan selvittää ennen hallin rakentamista. CFD-mallinnuksella ja laskentayhtälöillä tutkitaan lämpötila- ja virtauskentät sekä hallin läpi virtaava ilmamäärä ja ilmanvaihtoaste. Tulo- ja poistoilma-aukkojen mitoitukseen ja sijoitukseen liittyvät suunnitteluarvot varmennetaan sekä löydetään ilmanvaihdon ongelmakohdat. Ongelmakohtia tutkitaan ja niille esitetään parannusehdotukset.

Stability Evaluation of a Paper Machine Wet End

Työssä analysoidaanprosessin vaikutusta paperikoneen stabiiliuteen. Kaksi modernia sanomalehtipaperikonetta analysoitiin ja sen perusteella molemmista prosesseista rakennettiin fysiikan lakeihin perustuvat simulointimallit APROS Paper simulointiohjelmistolla. Työn tavoitteena on selvittää, miten kyseisten koneiden prosessit eroavat toisistaan ja arvioida, miten havaitut erot vaikuttavat prosessien stabiiliuteen. Työssä tarkastellaan periodisten häiriöiden vaimenemista prosessissa. Simuloinnissa herätteenä käytettiin puhdasta valkoista kohinaa, jonka avulla eri taajuistenperiodisten häiriöiden vaimenemista analysoitiin. Prosessien häiriövasteet esitetään taajuuskoordinaatistossa. Suurimmat erot prosessien välillä löytyivät viirakaivosta ja sen sekoitusdynamiikasta. Perinteisen viirakaivon todettiin muistuttavan käyttäytymiseltään sarjaan kytkettyjä ideaalisekoittimia, kun taas pienempitilavuuksisen fluumin todettiin käyttäytyvän lähes kuin putkiviive. Vaikka erotprosessitilavuudessa sekä viirakaivon sekoitusdynamiikassa olivat hyvin selkeät, havaittiin vain marginaalinen ero prosessin välillä periodisten häiriöiden vaimenemisessa, koska erot viiraretentiotasoissa vaikuttivat eniten simulointituloksia. Matalammalla viiraretentiolla operoivan paperikoneen todettiin vaimentavan tehokkaammin prosessihäiriöitä. Samalla retentiotasolla pienempitilavuuksisen prosessin todettiin vaimentavan hitaita prosessihäiriöitä marginaalisesti paremmin. Tutkituista paperikoneista toisella simuloitiin viiraosan vedenpoistomuutoksenvaikutusta viiraretentioon ja paperin koostumukseen. Lisäksi arvioitiin viiraretention säädön toimivuutta. Viiraosan listakengän vedenpoiston todettiin aiheuttavan merkittäviä sakeus- ja retentiohäiriöitä, mikäli sen avulla poistettavan kiintoaineen virtaus tuplaantuisi. Viiraretention säädön todettiin estävän häiriöiden kierron prosessissa, mutta siirtävän ne suoraan rainaan. Retention säädön eikuitenkaan todettu olevan suoranainen häiriön lähde.