Aineominaisuuksien mallinnustarkkuuden vaikutukset lämmönsiirtimien lämpöteknisessä simuloinnissa

Työssä on käsitelty fluidien aineominaisuuksien vaikutuksia paperikoneiden kuivatusosissa käytettävien lämmönsiirtimien lämpöteknisessä simuloinnissa. Pääkohteena selvitettiin kostean ilman ja veden fysikaalisien aineominaisuuksien mallinnustarkkuuden vaikutuksia lämpövirtaan lauhduttamattomissa ja lauhduttavissa tapauksissa. Asiaa tutkittiin tekemällä herkkyysanalyysi työssä kehitetyille termodynaamisille malleille. Perinteisen herkkyysanalyysin lisäksi herkkyyksiä tutkittiin myös Bayesiläisellä tilastoanalyysillä. Työssä käsiteltiin myös aineominaisuuksien käyttäytymistä ja mallintamista lämmönsiirtimissä. Kirjallisuudesta etsittiin aineominaisuusmallit, joilla kostean ilman ja veden fysikaalisia aineominaisuuksia voidaan kuvata riittävän tarkasti. Työssä havaittiin, että yksittäisistä aineominaisuuksista selkeästi suurimmat vaikutukset on ominaisentalpioiden mallinnuksen epätarkkuuksilla. Myös kaikkien aineominaisuuksien epätarkkuuksilla havaittiin olevan huomattavan suuret yhteisvaikutukset lämpövirran laskentatarkkuuteen. Viiden prosentin epätarkkuus kaikkien aineominaisuuksien mallinnuksessa johtaa 3 - 7 %:n epätarkkuuteen lämpövirran laskennassa. Näin ollen kaikkien aineominaisuuksien mallintamiseen tulee kiinnittää huomiota.

Meesauunin polttovyöhykkeen lämmönsiirron ja palamisen mallinnus

Tälle diplomityölle on antanut alkusysäyksen tarve kehittää käytössä olevaa meesauunin simulointiohjelmaa. Simulointiohjelma mallintaa meesauunin stationaaritilan toimintaa. Sillä voidaan tutkia uunin konstruktion muutosten vaikutuksia uunin toimintaan. Ohjelma on tehty 1980-luvun alkupuolella. Sen aikaisten tietokoneiden laskentatehojen vuoksi ohjelman käyttämään laskentamalliin on jouduttu tekemään joukko erilaisia yksinkertaistuksia laskenta-ajan lyhentämiseksi. Tässä diplomityössä keskityttiin tutkimaan meesauunin polttovyöhykkeen lämmönsiirron ja palamisen mallinnusta. Työssä luotiin aluksi tarvittavat massa- ja energiataseet sekä esitettiin tarvittavat lämmönsiirtoyhtälöt. Sen jälkeen kehitettiin uusi polttoaineen 1D-palamismalli. Palamismalli tehtiin VTT:n tekemien 3D-mallinnusten perusteella. Polttoaineina käytettiin maakaasua ja polttoöljyä. Uusi palamismalli lisättiin simulointiohjelmaan. Lisäksi simulointiohjelmasta muutettiin savukaasun emissiviteetin laskenta ja lämmönsiirto uunin ulkopinnasta ympäristöön. Tuloksena saatiin aikaan uusi tarkempi kuvaus lämmönsiirrosta ja palamisesta meesauunissa.

Passiivisten turvallisuusjärjestelmien luontaiset virhetoiminnot ydinvoimalaitoksessa

Uusissa ydinvoimalaitostyypeissä aiotaan käyttää aiempaa enemmän passiivisia turvallisuusjärjestelmiä. Näistä järjestelmistä on vielä vähän käyttökokemusta aktiivisiin turvallisuusjärjestelmiin verrattuna. Työssä tarkastellaan passiivisten turvallisuusjärjestelmien toimintaa sekä etsitään niiden mahdollisia luontaisia vikatilanteita. Luontaisten vikatilanteiden seurauksia järjestelmän suorituskykyyn arvioitiin yksinkertaisilla laskuilla ja mallintamalla RELAP5/MOD3.2.2 beta -termohydrauliikkaohjelmalla. Tarkastelu rajattiin kahden erityyppisen ydinvoimalaitoksen passiivisiin turvallisuusjärjestelmiin. Turvallisuusjärjestelmien suuntaa antavat mitat ja käyttötilanteiden parametrit saatiin laitosvalmistajien laitoskuvauksista. Osoittautui, että vikatilanteissa passiivisissa turvallisuusjärjestelmissä geometrialla on merkittävä vaikutus järjestelmän kapasiteettiin. Tarkasteluissa saatiin myös selville, että laitosmittakaavassa painovoimaisen hätälisävesijärjestelmän turvallisuustoiminto voi toteutua vaikka esiintyisi lyhytaikaisia toimintahäiriöitä, kuten lauhtumista hätälisävesisäiliössä. Sen sijaan lämmönsiirtopiirin virtausreittien tukkeutuminen voi olla fysikaalisesti merkittävä toimintaa haittaava tekijä.

Voimalaitoskattilan pyörivän regeneratiivisen ilman esilämmittimen suorituskyky ja mitoitus

Diplomityön tarkoituksena oli kehittää laskentaohjelma pyörivälle regeneratiiviselle lämmönsiirtimelle. Työ tehtiin Foster Wheeler Energia Oy:n Varkauden toimipisteessä. Työn ensimmäisessä osuudessa tutkittiin kirjallisuuden avulla regeneraattoreihin liittyvää teoriaa. Tämä osuus sisältää regeneratiivisen lämmönsiirron perusteita, regeneraattoreiden luokittelua, sektorijakoja, vuotoja, lämpöpintojen geometriaa ja likaantumista. Soveltavassa osassa tehtiin laskentaohjelma, jonka avulla voidaan laskea pyörivän regeneraattorin mitoitus- ja suorituskykylaskuja. Lisäksi ohjelman avulla voidaan laskea regeneraattorin lämpötilaprofiili.

Development of Furnace Model for Bubbling Fluidized Bed Boiler

The aim of this thesis was to develop a model, which can predict heat transfer, heat release distribution and vertical temperature profile of gas phase in the furnace of a bubbling fluidized bed (BFB) boiler. The model is based on three separate model components that take care of heat transfer, heat release distribution and mass and energy balance calculations taking into account the boiler design and operating conditions. The model was successfully validated by solving the model parameters on the basis of commercial size BFB boiler test run information and by performing parametric studies with the model. Implementation of the developed model for the Foster Wheeler BFB design procedures will require model validation with existing BFB database and possibly more detailed measurements at the commercial size BFB boilers.

Characterization Method for Comminution Behavior of Limestone and Fuel Ash in a Circulating Fluidized Bed Boiler

Työn tavoite oli kehittää karakterisointimenetelmät kalkkikiven ja polttoaineen tuhkan jauhautumisen ennustamiselle kiertoleijukattilan tulipesässä. Kiintoainekäyttäytymisen karakterisoinnilla ja mallintamisella voidaan tarkentaa tulipesän lämmönsiirron ja tuhkajaon ennustamista. Osittain kokeelliset karakterisointimenetelmät perustuvat kalkkikiven jauhautumiseen laboratoriokokoluokan leijutetussa kvartsiputkireaktorissa ja tuhkan jauhatumiseen rotaatiomyllyssä. Karakterisointimenetelmät ottavat huomioon eri-laiset toimintaolosuhteet kaupallisen kokoluokan kiertoleijukattiloissa. Menetelmät kelpoistettiin kaupallisen kokoluokan kiertoleijukattiloista mitattujen ja fraktioittaisella kiintoainemallilla mallinnettujen taseiden avulla. Kelpoistamistaseiden vähäisyydestä huolimatta karakterisointimenetelmät arvioitiin virhetarkastelujen perusteella järkeviksi. Karakterisointimenetelmien kehittämistä ja tarkentamista tullaan jatkamaan.

Turbulenttisen putkivirtauksen numeerinen mallinnus.

Numeerisella mallinnuksella on tavoitteena täydentää ja korvata kokeellista tutkimusta. Tässä tutkimuksessa on mallinnettu CFX 4.1- ja CFX 4.2-ohjelmien avulla lämmönsiirtoa putken sisäpinnalla. Virtausaineena putkessa on käytetty vettä ja vesi-monopropyleeniglykoliliuosta. Tarkasteltujen virtaustapausten Reynoldsin luku vaihtelee 200 - 30000. Kun glykolipitoisuus on suuri ja liuoksen lämpötila on pieni virtaus on laminaarista ja tällöin lämmönsiirtymiskerroin on pieni. Lämmönsiirron tehostamiseksi putkeen on asennettu turbulaattorilanka. Työssä on selvitetty edellytyksiä mallintaa hydraulisesti sileässä putkessa tapahtuvaa virtausta. Reynoldsin luvun ollessa alle 2300 mallinnuksessa on käytetty laminaarimallia. Reynoldsin luvuilla 2300-30000, turbulenttisella alueella, on käytetty pienten Reynoldsin luvun k-ɛ-mallia. Malli vaatii toimiakseen tiheän laskentaverkon putken seinämän läheisyydessä. Tarkastellulla alueella virtauksen ja lämmönsiirron mallinnuksen tulokset ovat vastaavat kuin teorian perusteella lasketut ja kokeellisista mittauksista saatavat tulokset. Lämmönsiirron tehostamiseksi putkeen on asennettu turbulaattorilanka. Tässä työssä on numeerisin menetelmin (pienten Reynoldsin luvun k-ɛ-malli ja k-ɛ-malli) suoritettu laskentaa yhdellä turbulaattorilankarakenteella. Laskennan vertailuaineistona on käytetty aikaisemmasta kokeellisesta tutkimuksesta saatua mittausdataa. Kokeellisessa tutkimuksessa turbulaattorirakenteena on käytetty putken seinämällä kiertyvää turbulaattorilankaa. Todellinen kolmiulotteinen geometria osoittautui vaikeaksi mallintaa. Toimivaa mallia ei ollut mahdollista toteuttaa aikataulun puitteissa ja mallin laskentakapasiteetin tarve kasvoi liian suureksi. Lankarakenne yksinkertaistettiin tasavälein toistuvaksi riparakenteeksi, joka on helpompi mallintaa aksisymmetriaa käyttäen kaksiulotteisena. Mallin tuloksista painehäviö asettuu kirjallisuudesta saatavan vertailuaineiston kanssa samalle tasolle, mutta lämmönsiirtymiskerroin on vertailuaineistoa huomattavasti suurempi.

Polttoaineseoksen ja prosessiolosuhteiden vaikutus kerrosleijukattilan tulistinalueen likaantumiseen

Työn tavoitteena oli löytää yhteys kerrosleijukattilaan syötetyn polttoaineseoksen sekä kattilan likaantumisen välille. Tulistinalueen likaantumista tutkittiin kahden kerrostumasondin avulla 28 päivää kestäneellä mittausjaksolla. Mittausten aikana otettiin näytteitä polttoaineseoksesta, lentotuhkasta sekä nuohouksen aikaisesta tuhkavirrasta. Mittausjakson jälkeen myös sondien tuhkakerrostumien pitoisuudet määritettiin. Mittausjakson aikana sondien jättöpinnoille muodostui tuhkakerrostumaa, joka voitiin poistaa nuohouksella, kun taas tulopinnalle syntyi pysyvää kerrostumaa. Nuohouksella irtoavan tuhkakerrostuman havaittiin sisältävän suuremmat pitoisuudet piitä, alumiinia, natriumia ja kaliumia kuin kattilan läpi jatkuvasti kulkevan lentotuhkan. Tuhkakerrostuma, jota ei nuohouksella saatu poistettua, sisälsi enemmän natriumia, rikkiä, kaliumia ja lyijyä kuin muut tuhkanäytteet. Polttoaineista turpeella oli merkittävin vaikutus likaantumiseen. Turpeen osuuden ollessa suurimmillaan jäivät jättöpinnan tuhkakerrostuman lämpövastukset pienemmiksi kuin muulloin eli lyhytaikaista kerrostumaa syntyi tällöin vähemmän. Pysyvän kerrostuman kasvu hidastui, kun turpeen osuus oli suuri, ja jopa pysähtyi, kun turpeen osuus oli 42-51 %. Prosessiolosuhteista tutkittiin kattilan kuorman vaikutusta likaantumiseen. Havaittiin, että ajettaessa kattilaa isommilla kuormilla, syntyi lyhytaikaista kerrostumaa vähemmän kuin muulloin.

Kipsin kiteytymiseen vaikuttavat tekijät ammoniumsulfaattikiteyttämön lämmönvaihtimissa

Työssä pyrittiin selvittämään syitä kipsin saostumiseen ammoniumsulfaattikiteyttämön putkilämmönvaihtimien pinnalle ja miten epätoivottua saostumista voitaisiin estää. Lämmönvaihtimissa virtaa ammoniumsulfaattia, jossa on epäpuhtautena kalsiumia, joka saostuu pinnoille kalsiumsulfaattina. Kirjallisuusosassa tarkasteltiin kiteytymisen mekanismia ja kipsin kiteytymiseen vaikuttavia tekijöitä. Saostumien estoaineita ja niiden vaikutusta kipsin kiteytymiseen sekä kipsin liukoisuutta ammoniumsulfaattiliuoksessa käsiteltiin myös kirjallisuusosassa. Kipsin kiteytymiseen vaikuttavia tekijöitä selvitettiin laboratoriokokeilla, joissa pyrittiin simuloimaan lämmönvaihdinta lämpövastuksella. Laboratoriokokeissa kokeiltiin erilaisia saostuman estoaineita ja pyrittiin löytämään prosessiin mahdollisimman tehokas kipsin kiteytymisen estoaine. Lämmönvaihtimien toiminnan tehokkuutta eli muodostuneen saostuman vaikutusta lämmönsiirtymiseen tutkittiin veden luovuttaman lämpövirran avulla. Lämmönvaihtimien tukkeutumista selvitettiin putkien vaihdon tarpeen perusteella. Kalsiumpitoisuuden vaihteluja prosessivirroissa selvitettiin kalsiumtaseen avulla. Saostumiseen vaikuttavien tekijöiden lisäksi selvitettiin mistä ja kuinka paljon kalsiumia kulkeutuu prosessiin ja poistuu sieltä. Työn tarkoituksena oli löytää ratkaisu, jolla epätoivottua saostumista lämmönvaihdin-ten pinnoille pystyttäisiin vähentämään joko kemiallisesti tai muuttamalla prosessi-muuttujia. Kalsiumia havaittiin olevan eniten pelkistämön sisäisissä ammoniumsul-faattiliuoskierroissa. Kalsiumtaseen perusteella kalsiumia poistuu pelkistämöltä eniten kipsinä ammoniumsulfaattituotteen mukana. Laboratoriokokeissa havaittiin polykar-boksylaattien estävän kipsin kasvua parhaiten, joskin estoaineen oikealla annostuksel-la havaittiin olevan suuri vaikutus. Lämmönvaihtimien saostuman havaittiin olevan kipsin ja glauberiitin seos. Vaippapuolen luovuttamien lämpövirran arvojen perusteel-la pystyttiin seuraamaan putkien tukkeutumista.

Cooling system optimization in medium-speed engines

Tässä työssä optimoidaan keskinopean Wärtsilä 32 -dieselmoottorin jäähdytysjärjestelmää ja tutkitaan taajuusmuuttajien käyttömahdollisuutta kiertopumppujen yhteydessä niin, että järjestelmässä saataisiin kiertämään vain kulloinkin tarvittava määrä vettä. Tutkimuksen mallinnus on toteutettu laatimalla aiemmin käytössä olleista yksinkertaisista simulointimalleista yksi malli, johon on sisällytetty sekä virtauksen että lämmönsiirron laskenta, jotka on aiemmin mallinnettu erillisillä ohjelmilla. Diplomityö on osa projektia, joka on tehty Sähkötekniikan osaston tutkijan Mikko Pääkkösen kanssa yhteistyössä. Tämän diplomityö keskittyy lähinnä virtausteknisiin ja lämmönsiirtoon liittyviin asioihin, kun taas sähkötekniikan osuus on esitetty Mikko Pääkkösen raportissa. Tulosten perustella voidaan sanoa, että taajuusmuuttajakäyttö kannattaa kiertopumppujen yhteydessä. Käyttämällä pumppujen virtaussäätöä voidaan jäähdytysjärjestelmästä jättää monia komponentteja, kuten termostaattiventtiilejä pois. Mallinnetut yksinkertaiset piiriratkaisut näyttävät toimivan ainakin yleisellä tasolla. Tutkimusta pumppujen säädöstä ja tässä projektissa luoduista jäähdytysjärjestelmäkonfiguraatioista kannattaa jatkaa.

Fischer-Tropsch Slurry Bubble Column Reactor Modeling

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

Ydinreaktorin painesäiliön ulkoisen jäähdytyksen lämmönsiirtokokeita

The purpose of this work was to design and carry out thermal-hydraulic experiments dealing with overcooling transients of a VVER-440-type nuclear reactor pressure vessel. Sudden overcooling accident could have negative effect on the mechanical strength of the pressure vessel. If part of the pressure vessel is compromised, the intense pressure inside a pressurized water reactor could cause the wall to fracture. Information on the heat transfer along the outside of the pressure vessel wall is necessary for stress analysis. Basic knowledge of the overcooling accident and heat transfer types on the outside of the pressure vessel is presented as background information. Test facility was designed and built based to study and measure heat transfer during specific overcooling scenarios. Two test series were conducted with the first one concentrating on the very beginning of the transient and the second one concentrating on steady state heat transfer. Heat transfer coefficients are calculated from the test data using an inverse method, which yields better results in fast transients than direct calculation from the measurement results. The results show that heat transfer rate varies considerably during the transient, being very high in the beginning and dropping to steady state in a few minutes. The test results show that appropriate correlations can be used in future analysis.

Modeling the transport phenomena within arterial wall: porous media approach

The transport of macromolecules, such as low-density lipoprotein (LDL), and their accumulation in the layers of the arterial wall play a critical role in the creation and development of atherosclerosis. Atherosclerosis is a disease of large arteries e.g., the aorta, coronary, carotid, and other proximal arteries that involves a distinctive accumulation of LDL and other lipid-bearing materials in the arterial wall. Over time, plaque hardens and narrows the arteries. The flow of oxygen-rich blood to organs and other parts of the body is reduced. This can lead to serious problems, including heart attack, stroke, or even death. It has been proven that the accumulation of macromolecules in the arterial wall depends not only on the ease with which materials enter the wall, but also on the hindrance to the passage of materials out of the wall posed by underlying layers. Therefore, attention was drawn to the fact that the wall structure of large arteries is different than other vessels which are disease-resistant. Atherosclerosis tends to be localized in regions of curvature and branching in arteries where fluid shear stress (shear rate) and other fluid mechanical characteristics deviate from their normal spatial and temporal distribution patterns in straight vessels. On the other hand, the smooth muscle cells (SMCs) residing in the media layer of the arterial wall respond to mechanical stimuli, such as shear stress. Shear stress may affect SMC proliferation and migration from the media layer to intima. This occurs in atherosclerosis and intimal hyperplasia. The study of blood flow and other body fluids and of heat transport through the arterial wall is one of the advanced applications of porous media in recent years. The arterial wall may be modeled in both macroscopic (as a continuous porous medium) and microscopic scales (as a heterogeneous porous medium). In the present study, the governing equations of mass, heat and momentum transport have been solved for different species and interstitial fluid within the arterial wall by means of computational fluid dynamics (CFD). Simulation models are based on the finite element (FE) and finite volume (FV) methods. The wall structure has been modeled by assuming the wall layers as porous media with different properties. In order to study the heat transport through human tissues, the simulations have been carried out for a non-homogeneous model of porous media. The tissue is composed of blood vessels, cells, and an interstitium. The interstitium consists of interstitial fluid and extracellular fibers. Numerical simulations are performed in a two-dimensional (2D) model to realize the effect of the shape and configuration of the discrete phase on the convective and conductive features of heat transfer, e.g. the interstitium of biological tissues. On the other hand, the governing equations of momentum and mass transport have been solved in the heterogeneous porous media model of the media layer, which has a major role in the transport and accumulation of solutes across the arterial wall. The transport of Adenosine 5´-triphosphate (ATP) is simulated across the media layer as a benchmark to observe how SMCs affect on the species mass transport. In addition, the transport of interstitial fluid has been simulated while the deformation of the media layer (due to high blood pressure) and its constituents such as SMCs are also involved in the model. In this context, the effect of pressure variation on shear stress is investigated over SMCs induced by the interstitial flow both in 2D and three-dimensional (3D) geometries for the media layer. The influence of hypertension (high pressure) on the transport of lowdensity lipoprotein (LDL) through deformable arterial wall layers is also studied. This is due to the pressure-driven convective flow across the arterial wall. The intima and media layers are assumed as homogeneous porous media. The results of the present study reveal that ATP concentration over the surface of SMCs and within the bulk of the media layer is significantly dependent on the distribution of cells. Moreover, the shear stress magnitude and distribution over the SMC surface are affected by transmural pressure and the deformation of the media layer of the aorta wall. This work reflects the fact that the second or even subsequent layers of SMCs may bear shear stresses of the same order of magnitude as the first layer does if cells are arranged in an arbitrary manner. This study has brought new insights into the simulation of the arterial wall, as the previous simplifications have been ignored. The configurations of SMCs used here with elliptic cross sections of SMCs closely resemble the physiological conditions of cells. Moreover, the deformation of SMCs with high transmural pressure which follows the media layer compaction has been studied for the first time. On the other hand, results demonstrate that LDL concentration through the intima and media layers changes significantly as wall layers compress with transmural pressure. It was also noticed that the fraction of leaky junctions across the endothelial cells and the area fraction of fenestral pores over the internal elastic lamina affect the LDL distribution dramatically through the thoracic aorta wall. The simulation techniques introduced in this work can also trigger new ideas for simulating porous media involved in any biomedical, biomechanical, chemical, and environmental engineering applications.

Fouling tendency in biomass-fired boilers

Although it is well known that firing biomass fuels leads to increased deposition buildup on heat transfer surfaces in boiler compared with firing coal, existing empirical knowledge about combustion of different types of biofuels is limited. The aim of this study is to give greater awareness and understanding of the circumstances which are able to decrease considerably deposition build up on heat transfer surfaces when firing different types of biofuels. The crucial part of this thesis is experimental investigation of fouling tendency while firing biomass fuels, such as straw, bark, and peat having different chemical composition. In order to give comprehensive overview of ash deposition phenomena the number of not less important issues such as mechanisms of ash deposition, effect of fouling on heat transfer, and design of boilers subjected to ash buildup were examined as well. The answers obtained in this study may be a step towards a better knowledge of firing biofuels as separately as in mixtures, and may provide solutions for successful combustion technique of biomass fuels.

CFD simulation of complex phenomena containing suspensions and flow throughporous media

There is an increasing reliance on computers to solve complex engineering problems. This is because computers, in addition to supporting the development and implementation of adequate and clear models, can especially minimize the financial support required. The ability of computers to perform complex calculations at high speed has enabled the creation of highly complex systems to model real-world phenomena. The complexity of the fluid dynamics problem makes it difficult or impossible to solve equations of an object in a flow exactly. Approximate solutions can be obtained by construction and measurement of prototypes placed in a flow, or by use of a numerical simulation. Since usage of prototypes can be prohibitively time-consuming and expensive, many have turned to simulations to provide insight during the engineering process. In this case the simulation setup and parameters can be altered much more easily than one could with a real-world experiment. The objective of this research work is to develop numerical models for different suspensions (fiber suspensions, blood flow through microvessels and branching geometries, and magnetic fluids), and also fluid flow through porous media. The models will have merit as a scientific tool and will also have practical application in industries. Most of the numerical simulations were done by the commercial software, Fluent, and user defined functions were added to apply a multiscale method and magnetic field. The results from simulation of fiber suspension can elucidate the physics behind the break up of a fiber floc, opening the possibility for developing a meaningful numerical model of the fiber flow. The simulation of blood movement from an arteriole through a venule via a capillary showed that the model based on VOF can successfully predict the deformation and flow of RBCs in an arteriole. Furthermore, the result corresponds to the experimental observation illustrates that the RBC is deformed during the movement. The concluding remarks presented, provide a correct methodology and a mathematical and numerical framework for the simulation of blood flows in branching. Analysis of ferrofluids simulations indicate that the magnetic Soret effect can be even higher than the conventional one and its strength depends on the strength of magnetic field, confirmed experimentally by Völker and Odenbach. It was also shown that when a magnetic field is perpendicular to the temperature gradient, there will be additional increase in the heat transfer compared to the cases where the magnetic field is parallel to the temperature gradient. In addition, the statistical evaluation (Taguchi technique) on magnetic fluids showed that the temperature and initial concentration of the magnetic phase exert the maximum and minimum contribution to the thermodiffusion, respectively. In the simulation of flow through porous media, dimensionless pressure drop was studied at different Reynolds numbers, based on pore permeability and interstitial fluid velocity. The obtained results agreed well with the correlation of Macdonald et al. (1979) for the range of actual flow Reynolds studied. Furthermore, calculated results for the dispersion coefficients in the cylinder geometry were found to be in agreement with those of Seymour and Callaghan.