Working fluid selection and design of small-scale waste heat recovery systems based on organic Rankine cycles

Demand for the use of energy systems, entailing high efficiency as well as availability to harness renewable energy sources, is a key issue in order to tackling the threat of global warming and saving natural resources. Organic Rankine cycle (ORC) technology has been identified as one of the most promising technologies in recovering low-grade heat sources and in harnessing renewable energy sources that cannot be efficiently utilized by means of more conventional power systems. The ORC is based on the working principle of Rankine process, but an organic working fluid is adopted in the cycle instead of steam. This thesis presents numerical and experimental results of the study on the design of small-scale ORCs. Two main applications were selected for the thesis: waste heat re- covery from small-scale diesel engines concentrating on the utilization of the exhaust gas heat and waste heat recovery in large industrial-scale engine power plants considering the utilization of both the high and low temperature heat sources. The main objective of this work was to identify suitable working fluid candidates and to study the process and turbine design methods that can be applied when power plants based on the use of non-conventional working fluids are considered. The computational work included the use of thermodynamic analysis methods and turbine design methods that were based on the use of highly accurate fluid properties. In addition, the design and loss mechanisms in supersonic ORC turbines were studied by means of computational fluid dynamics. The results indicated that the design of ORC is highly influenced by the selection of the working fluid and cycle operational conditions. The results for the turbine designs in- dicated that the working fluid selection should not be based only on the thermodynamic analysis, but requires also considerations on the turbine design. The turbines tend to be fast rotating, entailing small blade heights at the turbine rotor inlet and highly supersonic flow in the turbine flow passages, especially when power systems with low power outputs are designed. The results indicated that the ORC is a potential solution in utilizing waste heat streams both at high and low temperatures and both in micro and larger scale appli- cations.

Benefits of gas driven ground source heat pump systems

This thesis studies energy efficiencies and technical properties of gas driven ground source heat pumps and pump systems. The research focuses on two technologies: gas engine driven compressor heat pump and thermally driven gas absorption heat pump. System consist of a gas driven compressor or absorption ground source heat pump and a gas condensing boiler, which covers peak load. The reference system is a standard electrically powered compressor heat pump with electric heating elements for peak load. The systems are compared through primary energy ratios. Coefficient of performances of different heat pump technologies are also compared. At heat pump level, gas driven heat pumps are having lower coefficient of performances as compared with corresponding electric driven heat pump. However, gas heat pumps are competitive when primary energy ratios, where electricity production losses are counted in, are compared. Technically, gas heat pumps can potentially achieve a slightly higher temperatures with greater total energy efficiency as compared to the electric driven heat pump. The primary energy ratios of gas heat pump systems in relation to EHP-system improves when the share of peak load increases. Electric heat pump system's overall energy efficiency is heavily dependent on the electricity production efficiency. Economy as well as CO2-emissions were not examined in this thesis, which however, would be good topics for further study.

Failure modes of passive decay heat removing safety systems of modern nuclear power plants

The purpose of this master’s thesis is to gain an understanding of passive safety systems’ role in modern nuclear reactors projects and to research the failure modes of passive decay heat removal safety systems which use phenomenon of natural circulation. Another purpose is to identify the main physical principles and phenomena which are used to establish passive safety tools in nuclear power plants. The work describes passive decay heat removal systems used in AES-2006 project and focuses on the behavior of SPOT PG system. The descriptions of the main large-scale research facilities of the passive safety systems of the AES-2006 power plant are also included. The work contains the calculations of the SPOT PG system, which was modeled with thermal-hydraulic system code TRACE. The dimensions of the calculation model are set according to the dimensions of the real SPOT PG system. In these calculations three parameters are investigated as a function of decay heat power: the pressure of the system, the natural circulation mass flow rate around the closed loop, and the level of liquid in the downcomer. The purpose of the calculations is to test the ability of the SPOT PG system to remove the decay heat from the primary side of the nuclear reactor in case of failure of one, two, or three loops out of four. The calculations show that three loops of the SPOT PG system have adequate capacity to provide the necessary level of safety. In conclusion, the work supports the view that passive systems could be widely spread in modern nuclear projects.

Structural and operational optimisation : applications in energy systems

In this work mathematical programming models for structural and operational optimisation of energy systems are developed and applied to a selection of energy technology problems. The studied cases are taken from industrial processes and from large regional energy distribution systems. The models are based on Mixed Integer Linear Programming (MILP), Mixed Integer Non-Linear Programming (MINLP) and on a hybrid approach of a combination of Non-Linear Programming (NLP) and Genetic Algorithms (GA). The optimisation of the structure and operation of energy systems in urban regions is treated in the work. Firstly, distributed energy systems (DES) with different energy conversion units and annual variations of consumer heating and electricity demands are considered. Secondly, district cooling systems (DCS) with cooling demands for a large number of consumers are studied, with respect to a long term planning perspective regarding to given predictions of the consumer cooling demand development in a region. The work comprises also the development of applications for heat recovery systems (HRS), where paper machine dryer section HRS is taken as an illustrative example. The heat sources in these systems are moist air streams. Models are developed for different types of equipment price functions. The approach is based on partitioning of the overall temperature range of the system into a number of temperature intervals in order to take into account the strong nonlinearities due to condensation in the heat recovery exchangers. The influence of parameter variations on the solutions of heat recovery systems is analysed firstly by varying cost factors and secondly by varying process parameters. Point-optimal solutions by a fixed parameter approach are compared to robust solutions with given parameter variation ranges. In the work enhanced utilisation of excess heat in heat recovery systems with impingement drying, electricity generation with low grade excess heat and the use of absorption heat transformers to elevate a stream temperature above the excess heat temperature are also studied.

Effects of Fin Parameters on Gas-Side Heat Transfer Performance of H-Type Finned Tube Bundle in a Waste Heat Recovery Boiler

Alfa Laval Aalborg Oy designs and manufactures waste heat recovery systems utilizing extended surfaces. The waste heat recovery boiler considered in this thesis is a water-tube boiler where exhaust gas is used as the convective heat transfer medium and water or steam flowing inside the tubes is subject to cross-flow. This thesis aims to contribute to the design of waste heat recovery boiler unit by developing a numerical model of the H-type finned tube bundle currently used by Alfa Laval Aalborg Oy to evaluate the gas-side heat transfer performance. The main objective is to identify weaknesses and potential areas of development in the current H-type finned tube design. In addition, numerical simulations for a total of 15 cases with varying geometric parameters are conducted to investigate the heat transfer and pressure drop performance dependent on H-type fin geometry. The investigated geometric parameters include fin width and height, fin spacing, and fin thickness. Comparison between single and double tube type configuration is also conducted. Based on the simulation results, the local heat transfer and flow behaviour of the H-type finned tube is presented including boundary layer development between the fins, the formation of recirculation zone behind the tubes, and the local variations of flow velocity and temperature within the tube bundle and on the fin surface. Moreover, an evaluation of the effects of various fin parameters on heat transfer and pressure drop performance of H-type finned tube bundle has been provided. It was concluded that from the studied parameters fin spacing and fin width had the most significant effect on tube bundle performance and the effect of fin thickness was the least important. Furthermore, the results suggested that the heat transfer performance would increase due to enhanced turbulence if the current double tube configuration is replaced with single tube configuration, but further investigation and experimental measurements are required in order to validate the results.

Vaneri- ja kertopuutehtaiden kuivatuskoneiden lämmöntalteenottojen järjestelmäselvitys

Vanerin tai kertopuun valmistusprosessissaviilun kuivaukseen käytetään suurin osa koko valmistusprosessin primäärienergiasta. Viilunkuivauskoneessa viilun sisältämä vesi siirretään tyypillisesti prosessihöyryllä lämmitettyyn viilunkuivaajan kiertoilmaan höyrystämällä ja poistetaanviilunkuivaajasta poistoilman mukana. Viilunkuivaajan poistoilma on lämmintä jaerittäin suuren kosteuspitoisuutensa takia sisältää runsaasti energiaa. Tyypillisellä viilunkuivaajalla poistoilmaan sitoutunut lämpöteho vaihtelee prosessiolosuhteista riippuen välillä 2,7-5,7 MW. Diplomityössä tutkittiin viilunkuivaajan poistoilman sisältämän lämmön talteenottoa laitteistolla, johon kuuluu lämmöntalteenottopesuri, jossa poistoilmalla lämmitetään tuotantolaitoksen tukkipuun hautomon kiertovettä sekä ilma-ilma-lämmönsiirrin, jolla lämmitetään pesurista poistuvan ilman jäännöslämmöllä ulkoilmaa tehdassalin tuloilmakäyttöön. Työn tavoitteena oli kehittää lämmöntalteenottojärjestelmän suunnittelua, mitoitusta ja ajotapoja. Työssä analysoitiin teoreettisesti pesuria ja ilmalämmönsiirrintä, kehitettiin lämmöntalteenottopesurin simulointimenetelmä ja mitattiin toiminnassa olevia talteenottolaitteistoja. Tutkimuksessa todettiin lämmöntalteenottohyötysuhteen vaihtelevan lämmityskaudella välillä 50-70 %. Lämmöntalteenottolaitteiston pesurin veteen saatava teho riippuu ensisijaisesti viilunkuivaajan poistoilman lämpösisällöstä, joka on enimmäkseen kosteusriippuvainen ja ilmanvaihtoilmaan saatava teho ulkolämpö-tilan määräämästä tehontarpeesta. Pesurin vesijärjestelmän vaikutusmekanismit pesurin suorituskykyyn tunnistettiin ja niiden pohjalta annetaan suositukset mitoitukseen ja ajotapaan. Lämmöntalteenottolaitteiston lämpötehon tasapainottamiseen pesurin ja ilma-ilma-lämmönsiirtimen välillä mitoituksen avulla esitellään työkalut.

Maalämpöpumpun mitoitusjärjestelmän kehittäminen

Työn tavoitteena oli kehittää maalämpöpumppujärjestelmään kuuluvien komponenttien mitoitusta. Työ tehtiin Alufer Oy nimiselle yritykselle, joka on työskennellyt jo kolme vuotta maalämpöpumppujärjestelmän tuotekehityksen parissa. Maalämpöpumppujärjestelmä tullaan suunnittelemaan mahdollisimman suorituskykyiseksi ja joustavaksi. Suunnittelun lähtökohtana on, että lämmitysjärjestelmä on ns. matalalämpö-järjestelmä, joka käytännössä usein toteutetaan lattialämmityksenä. Ensimmäiseksi työssä on selvitetty mitä maalämpö on ja mitkä ovat yleisimmät maalämpöpumppujärjestelmän lämmönkeruuputkistojen asennustavat. Tällä hetkellä käytössä on joko vaakaan (maa, vesi) tai pystyyn asennettava lämmönkeruuputkisto (porakaivo). Seuraavaksi työssä on lähdetty selvittämään maalämpöpumppumarkkinoita Suomessa sekä selvitetty kolmen suurimman valmistajan Geopro Systemsin, Suomen Lämpöpumpputekniikan ja Ekowellin tuotteita sekä tekniikkaa. Työssä selvitetään myös muutamien Eurooppalaisten maiden markkinat. Mitoitusjärjestelmässä tarkastelu on aloitettu uudisrakennuksen lämmitystehon tarpeesta ja käyttöveden lämmityksen tarvitsemasta tehosta. Tarvittavan lämmitysenergian perusteella määriteltiin lämpöpumppujärjestelmään kuuluvat komponentit. Maalämpöpumppujärjestelmä koostuu seuraavista pääkomponenteista: höyrystin, kompressori, lauhdutin ja paisuntaventtiili. Höyrystimen tehon mitoituksessa on huomioitu lämmönkeruuputkistossa kulkevan nesteen aineominaisuudet, massavirta ja lämpötilaero höyrystimen nesteen ulostulon sekä sisään menon välillä. Kompressorin teho on määritetty valitun kylmäaineen (R407C) lg p-h piirroksesta tai määritetty teoreettisesti kompressorivalmistajien omista valintaohjelmista. Lauhduttimen teho on määritelty höyrystimen sekä kompressorin tehon summasta. Samalla määräytyy myös uudisrakennuksen lämmitystehontarve. Lopuksi työssä on käsitelty maalämpöpumppujärjestelmän kehitysmahdollisuuksia. Vaihtoehtoina on huomioitu tulistin, alijäähdytin ja varaaja, joilla voidaan huomattavasti parantaa maalämpöpumpun lämpökerrointa.

Lämpöpumppujärjestelmät integroidussa kylmä- ja lämpötehon tuotossa

Tässä työssä tutkittiin lämpöpumppujärjestelmiä, joilla tuotetaan samaan aikaan kylmä- ja lämpötehoa. Höyryn puristus lämpöpumppu on yleisimmin käytetty lämpöpumpputyyppi ja sen pääkomponentit ovat kompressori, lauhdutin, paisuntaventtiili ja höyrystin. Lämpöpumppu tuottaa samaan aikaan kylmätehoa höyrystimellä ja lämpötehoa lauhduttimella. Lämpöpumpun toiminta-arvoihin vaikuttaa valittujen lämpötilatasojen lisäksi voimakkaasti valitun kiertoaineen termodynaamiset ominaisuudet sekä kompressorin painesuhteeseen verrannollinen isentrooppihyötysuhde. Uusissa lämpöpumpuissa käytetään HFC yhdisteitä sekä sekoituksia kiertoaineina, mutta myös luonnolliset aineet, kuten ammoniakki, ovat lupaavia korvikkeita CFC yhdisteille. Sopivia sovelluskohteita kylmä- ja lämpötehon yhteistuotannolle ovat kauppa- ja asuinrakennukset, hotellit, toimistot, elintarviketeollisuus ja -myymälät sekä vierekkäiset jää- ja uimahallit ja hiihtoputket. Kylmä- ja lämpötehon yhteistuotannolla voitaisiin saavuttaa merkittäviä säästöjä ja päästövähennyksiä. Esimerkiksi jäähallien kylmäkoneiden lauhdelämmön hyödyntämisessä olisi Suomessa potentiaalia 6-10 miljoonan euron vuotuisiin säästöihin. Kylmä- ja lämpötehon yhteistuotanto voidaan toteuttaa hyödyntämällä kylmäkoneen lauhdelämpöä toisella lämpöpumpulla. Toinen vaihtoehto on käyttää eri tilojen samanaikaiseen lämmittämiseen ja jäähdyttämiseen HPS lämpöpumppua tai moniyksikköistä lämpöpumppua.

The bifurcational behaviour of the spatially distributed Rayleigh equation

At the present work the bifurcational behaviour of the solutions of Rayleigh equation and corresponding spatially distributed system is being analysed. The conditions of oscillatory and monotonic loss of stability are obtained. In the case of oscillatory loss of stability, the analysis of linear spectral problem is being performed. For nonlinear problem, recurrent formulas for the general term of the asymptotic approximation of the self-oscillations are found, the stability of the periodic mode is analysed. Lyapunov-Schmidt method is being used for asymptotic approximation. The correlation between periodic solutions of ODE and PDE is being investigated. The influence of the diffusion on the frequency of self-oscillations is being analysed. Several numerical experiments are being performed in order to support theoretical findings.

Uusiutuvan energian käyttömahdollisuudet Liikuntakeskus Pajulahdessa

Tämä diplomityö on tehty Lappeenrannan teknillisen yliopiston, Liikuntakeskus Pajulahden ja Nastolan Energiasäätiön yhteisprojektina. Työn tarkoituksena on tutkia aurinkosähköjärjestelmien taloudellista kannattavuutta ostosähkön korvaamisessa Liikuntakeskus Pajulahden kiinteistössä. Työssä tehdään laskelmat aurinkopaneelien tuotoista, tuotetun sähkön hinnasta, investointien takaisinmaksuajoista sekä tuotetun aurinkosähkön synnyttämistä hiilidioksidipäästövähenemistä. Lisäksi tarkastellaan energiaveron vaikutusta pientuotantolaitoksen kannattavuuteen. Tämän hetkisillä aurinkosähköjärjestelmien hinnoilla pystytään tuottamaan puhdasta, ympäristöystävällistä sähköä kilpailukykyiseen hintaan. Työssä tehdään yleiskatsaus myös muiden uusiutuviin energialähteisiin perustuvien tuotantomenetelmien käyttömahdollisuuksista. Yhtenä erityistarkastelunkohteena on Pajulahden monitoimihalliin vuonna 2008 suunniteltu vesistöä lämmönlähteenä käyttävä lämpöpumppujärjestelmä, jonka toteutukseen ei saatu tarvittavia lupia. Tarkastelussa tutkitaan myös aiemmin Suomessa toteutettuja vastaavia hankkeita ja niiden lupien saantia. Lisäksi on haettu tosiasioihin perustuvia argumentteja keskusteluun vesistölämmön hyödyntämisestä. Tarkastelun kohteena oleva vesistölämpöhanke voidaan todeta olevan erittäin kannattava, taloudellisesti sekä ympäristöllisesti. Tutkimuksen pohjalta Nastolan Energiasäätiölle tuotetaan havainnollinen ja selkeä aineisto aurinkoenergian hyödyntämismahdollisuuksista Nastolan kunnassa.

Asuinkerrostalon lämmöntuottojärjestelmän optimointi elinkaarikustannusten perusteella hyödynnettäessä lämpöpumpulla useita lämmönlähteitä

The basis of this thesis was to optimize heat pump that uses multiple heat sources to get competitive heating system for residential building when life cycle costs are considered. The objectives were to compile necessary information to calculate life cycle costs for heating system of residential building and start to compose of designing program for heat pump based heating systems. Examinations were made for the purchase energy need of residential building. Features of heat pump, considered refrigerant and potential heat sources were examined to find out heat production potential of heat pumps. Necessary information for life cycle cost calculation was also examined. Collected data was used in two case analyses to design selected heat production systems and calculate their life cycle costs. On the basis of case analyses heat pump based hybrid heat production systems are very competitive on life cycle cost comparison against district heating when residential building uses a lot of energy. New buildings use considerably less energy and achieved energy cost savings with heat pump systems may not be enough to cover the relatively high investment cost in reasonable time period compared to district heating system. The calculation method was found to require further development to at least include the cooling energy need of the building. Cooling demand will continue to grow in the future, which improves the heat pump based heat production systems competitiveness compared to other systems.

Particle transport method for convection-diffusion-reaction problems

Convective transport, both pure and combined with diffusion and reaction, can be observed in a wide range of physical and industrial applications, such as heat and mass transfer, crystal growth or biomechanics. The numerical approximation of this class of problemscan present substantial difficulties clue to regions of high gradients (steep fronts) of the solution, where generation of spurious oscillations or smearing should be precluded. This work is devoted to the development of an efficient numerical technique to deal with pure linear convection and convection-dominated problems in the frame-work of convection-diffusion-reaction systems. The particle transport method, developed in this study, is based on using rneshless numerical particles which carry out the solution along the characteristics defining the convective transport. The resolution of steep fronts of the solution is controlled by a special spacial adaptivity procedure. The serni-Lagrangian particle transport method uses an Eulerian fixed grid to represent the solution. In the case of convection-diffusion-reaction problems, the method is combined with diffusion and reaction solvers within an operator splitting approach. To transfer the solution from the particle set onto the grid, a fast monotone projection technique is designed. Our numerical results confirm that the method has a spacial accuracy of the second order and can be faster than typical grid-based methods of the same order; for pure linear convection problems the method demonstrates optimal linear complexity. The method works on structured and unstructured meshes, demonstrating a high-resolution property in the regions of steep fronts of the solution. Moreover, the particle transport method can be successfully used for the numerical simulation of the real-life problems in, for example, chemical engineering.

Methods for construction and analysis of computational models in systems biology : applications to the modelling of the heat shock response and the self-assembly of intermediate filaments

Systems biology is a new, emerging and rapidly developing, multidisciplinary research field that aims to study biochemical and biological systems from a holistic perspective, with the goal of providing a comprehensive, system- level understanding of cellular behaviour. In this way, it addresses one of the greatest challenges faced by contemporary biology, which is to compre- hend the function of complex biological systems. Systems biology combines various methods that originate from scientific disciplines such as molecu- lar biology, chemistry, engineering sciences, mathematics, computer science and systems theory. Systems biology, unlike “traditional” biology, focuses on high-level concepts such as: network, component, robustness, efficiency, control, regulation, hierarchical design, synchronization, concurrency, and many others. The very terminology of systems biology is “foreign” to “tra- ditional” biology, marks its drastic shift in the research paradigm and it indicates close linkage of systems biology to computer science. One of the basic tools utilized in systems biology is the mathematical modelling of life processes tightly linked to experimental practice. The stud- ies contained in this thesis revolve around a number of challenges commonly encountered in the computational modelling in systems biology. The re- search comprises of the development and application of a broad range of methods originating in the fields of computer science and mathematics for construction and analysis of computational models in systems biology. In particular, the performed research is setup in the context of two biolog- ical phenomena chosen as modelling case studies: 1) the eukaryotic heat shock response and 2) the in vitro self-assembly of intermediate filaments, one of the main constituents of the cytoskeleton. The range of presented approaches spans from heuristic, through numerical and statistical to ana- lytical methods applied in the effort to formally describe and analyse the two biological processes. We notice however, that although applied to cer- tain case studies, the presented methods are not limited to them and can be utilized in the analysis of other biological mechanisms as well as com- plex systems in general. The full range of developed and applied modelling techniques as well as model analysis methodologies constitutes a rich mod- elling framework. Moreover, the presentation of the developed methods, their application to the two case studies and the discussions concerning their potentials and limitations point to the difficulties and challenges one encounters in computational modelling of biological systems. The problems of model identifiability, model comparison, model refinement, model inte- gration and extension, choice of the proper modelling framework and level of abstraction, or the choice of the proper scope of the model run through this thesis.

Development of Craniofacial Bone Defect Reconstruction Implant Based on Fibre-Reinforced Composite with Photopolymerisable Resin Systems. Experimental studies in vitro and in vivo.

Reconstruction of defects in the craniomaxillofacial (CMF) area has mainly been based on bone grafts or metallic fixing plates and screws. Particularly in the case of large calvarial and/or craniofacial defects caused by trauma, tumours or congenital malformations, there is a need for reliable reconstruction biomaterials, because bone grafts or metallic fixing systems do not completely fulfill the criteria for the best possible reconstruction methods in these complicated cases. In this series of studies, the usability of fibre-reinforced composite (FRC) was studied as a biostable, nonmetallic alternative material for reconstructing artificially created bone defects in frontal and calvarial areas of rabbits. The experimental part of this work describes the different stages of the product development process from the first in vitro tests with resin-impregnated fibrereinforced composites to the in vivo animal studies, in which this FRC was tested as an implant material for reconstructing different size bone defects in rabbit frontal and calvarial areas. In the first in vitro study, the FRC was polymerised in contact with bone or blood in the laboratory. The polymerised FRC samples were then incubated in water, which was analysed for residual monomer content by using high performance liquid chromatography (HPLC). It was found that this in vitro polymerisation in contact with bone and blood did not markedly increase the residual monomer leaching from the FRC. In the second in vitro study, different adhesive systems were tested in fixing the implant to bone surface. This was done to find an alternative implant fixing system to screws and pins. On the basis of this study, it was found that the surface of the calvarial bone needed both mechanical and chemical treatments before the resinimpregnated FRC could be properly fixed onto it. In three animal studies performed with rabbit frontal bone defects and critical size calvarial bone defect models, biological responses to the FRC implants were evaluated. On the basis of theseevaluations, it can be concluded that the FRC, based on E-glass (electrical glass) fibres forming a porous fibre veil enables the ingrowth of connective tissues to the inner structures of the material, as well as the bone formation and mineralization inside the fibre veil. Bone formation could be enhanced by using bioactive glass granules fixed to the FRC implants. FRC-implanted bone defects healed partly; no total healing of defects was achieved. Biological responses during the follow-up time, at a maximum of 12 weeks, to resin-impregnated composite implant seemed to depend on the polymerization time of the resin matrix of the FRC. Both of the studied resin systems used in the FRC were photopolymerised and the heat-induced postpolymerisation was used additionally.

Mass transfer and particleinteractions in granular systems and suspension flows

The purpose of this study was to investigate some important features of granular flows and suspension flows by computational simulation methods. Granular materials have been considered as an independent state ofmatter because of their complex behaviors. They sometimes behave like a solid, sometimes like a fluid, and sometimes can contain both phases in equilibrium. The computer simulation of dense shear granular flows of monodisperse, spherical particles shows that the collisional model of contacts yields the coexistence of solid and fluid phases while the frictional model represents a uniform flow of fluid phase. However, a comparison between the stress signals from the simulations and experiments revealed that the collisional model would result a proper match with the experimental evidences. Although the effect of gravity is found to beimportant in sedimentation of solid part, the stick-slip behavior associated with the collisional model looks more similar to that of experiments. The mathematical formulations based on the kinetic theory have been derived for the moderatesolid volume fractions with the assumption of the homogeneity of flow. In orderto make some simulations which can provide such an ideal flow, the simulation of unbounded granular shear flows was performed. Therefore, the homogeneous flow properties could be achieved in the moderate solid volume fractions. A new algorithm, namely the nonequilibrium approach was introduced to show the features of self-diffusion in the granular flows. Using this algorithm a one way flow can beextracted from the entire flow, which not only provides a straightforward calculation of self-diffusion coefficient but also can qualitatively determine the deviation of self-diffusion from the linear law at some regions nearby the wall inbounded flows. Anyhow, the average lateral self-diffusion coefficient, which was calculated by the aforementioned method, showed a desirable agreement with thepredictions of kinetic theory formulation. In the continuation of computer simulation of shear granular flows, some numerical and theoretical investigations were carried out on mass transfer and particle interactions in particulate flows. In this context, the boundary element method and its combination with the spectral method using the special capabilities of wavelets have been introduced as theefficient numerical methods to solve the governing equations of mass transfer in particulate flows. A theoretical formulation of fluid dispersivity in suspension flows revealed that the fluid dispersivity depends upon the fluid properties and particle parameters as well as the fluid-particle and particle-particle interactions.