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.

Radiaalikompressorin spiraalivirtauksen numeerinen mallinnus

Diplomityössä mallinnetaan numeerisesti radiaalikompressorin spiraalin virtaus. Spiraalin tarkoituksena radiaalikompressorissa on kerätä tasaisesti virtaus diffuusorin kehältä. Spiraaliin on viime aikoina kiinnitetty enemmän huomiota, koska on havaittu, että kompressorin hyötysuhdetta voidaan parantaa spiraalia optimoimalla. Spiraalin toimintaa tarkastellaan kolmella eri massavirralla. Työn alussa käsitellään spiraalin toimintaperiaatteita. Numeerisena ratkaisijana käytetään Teknillisessä korkeakoulussa kehitettyä FIN-FLO -koodia. FINFLO -laskentaohjelmassa ratkaistaan Navier-Stokes yhtälöt kolmeulotteiselle laskenta-alueelle. Diskretointi perustuu kontrollitilavuus menetelmään. Työssä käsitellään laskentakoodin toimintaperiaatteitta. Turbulenssia mallinnetaan algebrallisella Baldwin-Lomaxin ja kahden yhtälön Chienin k-e turbulenssimalleilla. Laskentatuloksia verrataan Lappeenrannan teknillisessä korkeakoulussa tehtyihin mittauksiin kyseessä olevasta kompressorin spiraalista. Myös eri turbulenssimalleilla ja hilatasoilla saatuja tuloksia verrataan keskenään. Laskentatuloksien jälkikäsittelyä varten ohjelmoitiin neljä eri tietokoneohjelmaa. Laskennalla pyritään saamaan lisäselvyyttä virtauksen käyttäytymiseen spiraalissa ja erityisesti ns. kielen alueella. Myös kahden eri turbulenssimallin toimivuutta kompressorin numeerisessa mallinnuksessa tutkitaan.

Ilmanjaon päätelaitteen virtauskentän mallintaminen

Tässä työssä esitetään tuloilmalaitteen virtausten laskennallinen tarkastelu huonetilassa virtauslaskentaohjelmiston avulla. Tarkastelu suoritetaan kolmiulotteisena sekä isotermisille että termisille tapauksille, jolloin huonetilan joihinkin pintoihin asetetaan lämpövuo ja epäpuhtauslähde. Työn tarkoituksena on arvioida CFX -virtauslaskentaohjelmiston soveltuvuutta huonetilavirtausten mallintamiseen ja tarkastella huonetilan ilmanvaihdon tehokkuutta eri tapauksissa. Arviointi suoritetaan vertaamalla laskentatuloksia mittaustuloksiin. Turbulenssimallien vaikutusta on arvioitu käyttämällä nollan yhtälön ja standardi k-e turbulenssimalleja. CFX:n k-e turbulenssimallilla lasketut nopeuden alenemat eri tilavuusvirroilla vastasivat paremmin mittaustuloksia kuin nollan yhtälön turbulenssimallilla, jolla erot olivat huomattavia. Huonetilan ilmanvaihdon tehokkuus vaihteli suuresti tuloilmalaitteen sijoituspaikan ja lämpövuon mukaan.

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.

Sähkönmyyntisopimusten markkinalähtöinen hinnoittelu

Vuonna 1995 alkanut sähkömarkkinoiden vapautuminen on muuttanut sähköyhtiöiden myyntisopimusten hinnoittelua merkittävästi. Ennen markkinauudistusta voitiin sähkön myyntisopimukset hinnoitella perustuen oman sähkön tuotannon kustannuksiin ja haluttuun katteeseen. Nykyään sähköpörssissä noteerattava sähkön hinta muodostaa perustan kaikkien myyntisopimusten hinnoittelulle. Sähkön markkinahinnan lisäksi on myyntisopimusten hinnoittelussa otettava huomioon sähkömarkkinoiden ominaispiirteistä aiheutuvat riskit sähkön myyjälle. Tässä työssä mallinnetaan Lappeenrannan Energia Oy:n markkinalähtöiset hinnoittelumenetelmät kahdelle sähkönmyyntisopimustyypille. Lisäksi tutkitaan markkinalähtöisen hinnoittelun tärkeimpien riskikomponenttien, aluehintaeron sekä profiililisän, merkitystä markkinalähtöisten myyntisopimusten hinnoittelussa. Aluehintaeron hyväksikäyttöä myyntisopimusten hinnoittelussa on tutkittu selvittämällä Suomen hinta-alueen CfD-johdannaisten riskipreemiot. Profiililisän merkitystä myyntisopimusten hinnoittelussa on tutkittu havainnoimalla profiililisän muutoksia hinta- ja kulutusaikasarjoissa sekä suojaushinta ja tehotasossa tapahtuneiden muutosten suhteen. Ennustetun ja toteutuneen profiililisän eroja on tutkittu laskemalla ne seitsemälle Lappeenrannan Energia Oy:n merkittävälle asiakkaalle. Lisäksi on tarkasteltu profiililisän laskentaan tarvittavan hinta-aikasarjan mallintamiseen käytettyjen termiinituotteiden hintojen merkitystä lasketun profiililisän suuruuteen. Työn lopuksi esitetään kaksi vaihtoehtoista sähkösopimusten hinnoittelumenetelmää ja vertaillaan esitettyjä hinnoittelumenetelmiä keskenään. Työssä havaitaan, että aluehintaerolta suojautumiseen käytettävien CfD-johdannaisten avulla olisi ollut mahdollista lisätä markkinalähtöisten myyntisopimusten tuottoa viimeisen kolmen vuoden aikana. Suoritettujen herkkyysanalyysien perusteella voidaan todeta, että toteutuneen ja ennustetun profiililisän erot johtuvat laskentaan käytettävien hinta- sekä kulutusaikasarjojen epätarkkuudesta. Lappeenrannan Energia Oy:n käyttämät profiililisät osoittautuivat ex-post -tarkastelussa liian suuriksi yhtä asiakasta lukuun ottamatta. Lisäksi tarjousaikana laskettujen profiililisien voidaan katsoa muuttuvan täysin samassa suhteessa hinta-aikasarjan mallintamiseen käytettyjen termiinituotteiden volatiliteettien kanssa. Esiteltyjen vaihtoehtoisten myyntisopimusten hinnoittelumenetelmien voidaan katsoa antavan varsin samanlaisia tuloksia kuin Lappeenrannan Energia Oy:n nykyinen hinnoittelumenetelmä. Saatuihin tuloksiin vaikuttavat kuitenkin painokertoimien estimointiin käytetyn vuoden volatiliteetti sekä profiililisän laskentamenetelmä

Pellettien laatu ja ominaisuudet eri raaka-aineilla ja niiden seoksilla

Diplomityö tehtiin osana Vapon toteuttamaa monivuotista pelletin kehitysohjelmaa. Kehitysohjelma koostuu useista pienemmistä osaprojekteista, jotka täydentävät toinen toisiaan. Pellettien raaka-ainepohjan laajentaminen on eräs näistä osaprojekteista. Tutkimustyön tavoitteena oli selvittää erilaisten potentiaalisten bioraaka-aineiden soveltuvuutta pelletointiin joko sellaisenaan tai erilaisina seoksina. Raaka-aineiden pelletoitavuutta tutkittiin kenttäolosuhteissa mobiilipelletointilaitoksella. Laitoksen pääkomponentit muodostivat Kahl C 38–780 tasomatriisipuristin, jäähdytin ja täryseula. Pelletointikokeissa tutkittuja raaka-aineita olivat mäntysahanpuru, männynkuori, harvennusranka, haapa, koivu, jyrsinturve ja ruokohelpi. Raaka-aineiden irtotiheys käyttökosteudessa vaihteli välillä 73–244 kg/m3 ja keskimääräinen kosteuspitoisuus 6,5–15 %. Useissa tapauksissa säkitettyjä raaka-aineita säkkikostutettiin haluttuun kosteuspitoisuuteen ennen pelletointia. Säkkikostutettujen raaka-aineiden kosteuspitoisuudet vaihtelivat tällöin välillä 12–14 m- %. Valtaosa tutkituista raaka-aineista ja niiden seoksista pystyttiin pelletoimaan puristimen matriisilla 8/40 mm, jossa puristuskanavan halkaisija oli 8 mm ja kanavan suoran osan pituus 40 mm. Vaikeuksia tuotti ainoastaan pelkän koivupurun ja ruokohelven pelletointi. Käytetty matriisi oli kanavapituudeltaan liian pitkä koivupurun pelletointiin nostaen puristusvastuksen suureksi. Ruokohelven pelletoinnin vaikeudet johtuivat pääasiassa pelletointiin liian karkeasta raaka-aineesta. Myös matriisia 8/55 mm kokeiltiin, mutta se osoittautui liian ”tiukaksi” valtaosalle puuraaka-aineista. Ainoastaan jyrsinturpeen pelletointi onnistui tällä matriisilla. Männynkuoren pelletointia ei matriisilla 8/55 mm yritetty. Kenttäkokeissa valmistetuista pelleteistä määritettiin erilaisia ominaisuuksia, kuten keskipituus, kosteuspitoisuus, irtotiheys, hienoaineksen määrä ja käsittelykestävyys. Lujuus mitattiin sekä Ligno-testillä että CEN-rummutuslujuuden määrityksellä. Lisäksi pelleteille määritettiin alkuaineanalyysi, tuhkapitoisuus ja lämpöarvo ENAS Oy:n laboratoriossa Jyväskylässä. Ligno-testauksessa parhaimman luokan pelletin tulee yltää arvoon 97,5 %. Pelletoitaessa raaka-aineita ja niiden seoksia tasomatriisikoneella sopivalla matriisilla yllettiin usein näihin tai parempiin tuloksiin. Puumateriaaleilla raaka-aineen optimaalinen lähtökosteus oli välillä 12–14 m- % ja turpeella sekä ruokohelvellä 14–16 m- %. Pelletointi onnistui tällöin vaivattomasti, kunhan sopivat puristimen ajoparametrit oli löydetty. Pellettiä alkoi muodostua matriisin puristuskanavien lämpötilan kohotessa noin 70 ºC. Pellettien lämpötila stabiilitilanteessa heti pelletoinnin jälkeen oli useissa tapauksissa 80–90 ºC. Pelletoinnin aikainen tehontarve vaihteli välillä 90–150 kWh/t, ollen suurimmillaan irtotiheydeltään keveillä materiaaleilla. Raaka-aineen suuri partikkelikoko kasvatti puristimen tehontarvetta. Tämä havaittiin selvästi lisättäessä karkeaa ruokohelpisilppua eri raaka-aineiden joukkoon. Kestävyydeltään erinomaisia pellettejä saatiin, kun raaka-aineena oli jyrsinturve, harvennusranka tai mäntypuru. Varsinkin jyrsinturpeen ja harvennusrangan seoksesta valmistetut pelletit osoittautuivat erittäin kestäviksi. Myös jyrsinturpeen ja ruokohelven sekä mäntypurun ja ruokohelven seoksien pelleteille määritettiin hyviä kestävyysarvoja. Männynkuoresta valmistetut pelletit jäivät Ligno-testauksessa kestävyydeltään alle 97,5 % rajan. Pääsyynä tähän oli kuoren pelletointiin käytetyn matriisin 8/40 mm liian lyhyet puristuskanavat.

Eri raaka-ainemuotojen hyödyntäminen puumuovikomposiitin valmistuslinjastossa

Diplomityössä tutkittiin erilaisten raaka-ainemuotojen hyödyntämistä puumuovikomposiitin valmistuslinjastossa ja raaka-ainemuotojen vaikutusta valmistuslinjaston kannattavuuteen. Työssä etsittiin kannattavinta yhdistelmää raaka-ainemuotojen ja valmistusvaihtoehtojen välille. Työn johdanto osuudessa käsitellään aiheeseen liittyvää kirjallisuutta ja tutkimuksia. Siinä esitellään valmistusprosessin vaiheet ja raaka-aineet, niihin vaikuttavat tekijät, valmistuksen taloudellisuus ja investointikäsitteet. Käytännön osuus koostuu valmistusprosessin taloudellisesta analyysista ja laboratoriossa suoritetuista koeajoista. Taloudellisessa analyysissa selvitettiin raaka-ainemuotojen kannattavuudet eri valmistusvaihtoehdoille ja koeajoissa tutkittiin raaka-ainemuotojen todellista käyttäytymistä ekstruusiossa kahdella erilaisella profiililla. Taloudellisen analyysin perusteella paras kannattavuus valmistusprosessissa saavutetaan puupelleteillä ja yksivaiheinen valmistus on kannattavampaa kuin kaksivaiheinen valmistus kaikilla raaka-ainemuodoilla. Kaksivaiheinen valmistus voisi kuitenkin teoriassa olla kannattavampaa kuin yksivaiheinen. Koeajoissa suurin maksimituotos saavutettiin puumuovipuristeilla, mutta ajettavuuden kannalta sopivin raaka-ainemuoto linjalle olisi jauhemainen sahanpuru. Puumuovipuristeita käytettäessä oli erittäin vaikeaa löytää ekstruuderille sopivia muuttujien arvoja ja muuttujien säätäminen muuttui vaikeammaksi, kun käytössä oli monimutkaisempi terassilautaprofiili.

Large Break Blowdown Test Facility Study

Työn tarkoituksena on kerätä yhteen tiedot kaikista maailmalta löytyvistä ison LOCA:n ulospuhallusvaiheen tutkimiseen käytetyistä koelaitteistoista. Työn tarkoituksena on myös antaa pohjaa päätökselle, onko tarpeellista rakentaa uusi koelaitteisto nesterakenne-vuorovaikutuskoodien laskennan validoimista varten. Ennen varsinaisen koelaitteiston rakentamista olisi tarkoituksenmukaista myös rakentaa pienempi pilottikoelaitteisto, jolla voitaisiin testata käytettäviä mittausmenetelmiä. Sopivaa mittausdataa tarvitaan uusien CFD-koodien ja rakenneanalyysikoodien kytketyn laskennan validoimisessa. Näitä koodeja voidaan käyttää esimerkiksi arvioitaessa reaktorin sisäosien rakenteellista kestävyyttä ison LOCA:n ulospuhallusvaiheen aikana. Raportti keskittyy maailmalta löytyviin koelaitteistoihin, uuden koelaitteiston suunnitteluperusteisiin sekä aiheeseen liittyviin yleisiin asioihin. Raportti ei korvaa olemassa olevia validointimatriiseja, mutta sitä voi käyttää apuna etsittäessä validointitarkoituksiin sopivaa ison LOCA:n ulospuhallusvaiheen koelaitteistoa.

Fluid flow in T-junction of pipes

The aim of this work is to study flow properties at T-junction of pipe, pressure loss suffered by the flow after passing through T-junction and to study reliability of the classical engineering formulas used to find head loss for T-junction of pipes. In this we have compared our results with CFD software packages with classical formula and made an attempt to determine accuracy of the classical formulas. In this work we have studies head loss in T-junction of pipes with various inlet velocities, head loss in T-junction of pipes when the angle of the junction is slightly different from 90 degrees and T-junction with different area of cross-section of the main pipe and branch pipe. In this work we have simulated the flow at T-junction of pipe with FLUENT and Comsol Multiphysics and observed flow properties inside the T-junction and studied the head loss suffered by fluid flow after passing through the junction. We have also compared pressure (head) losses obtained by classical formulas by A. Vazsonyi and Andrew Gardel and formulas obtained by assuming T-junction as combination of other pipe components and observations obtained from software experiments. One of the purposes of this study is also to study change in pressure loss with change in angle of T-junction. Using software we can have better view of flow inside the junction and study turbulence, kinetic energy, pressure loss etc. Such simulations save a lot of time and can be performed without actually doing the experiment. There were no real life experiments made, the results obtained completely rely on accuracy of software and numerical methods used.

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.

La Vendetta : from "The Marriage of Figaro"

Äänitetty: 7.-8.5.1953, Los Angeles.

Per questa bella mano

Äänitetty: 7.-8.5.1953, Los Angeles.

Hai già vinta - Vedrò mentr'io : from "The Marriage of Figaro"

Äänitetty: 7.-8.5.1953, Los Angeles.

Se vuol ballare : from "The Marriage of Figaro"

Äänitetty: 7.-8.5.1953, Los Angeles.