Simulation of surface waves in porous media

We present a novel numerical algorithm for the simulation of seismic wave propagation in porous media, which is particularly suitable for the accurate modelling of surface wave-type phenomena. The differential equations of motion are based on Biot's theory of poro-elasticity and solved with a pseudospectral approach using Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. The time solver is a splitting algorithm that accounts for the stiffness of the differential equations. Due to the Chebyshev operator the grid spacing in the vertical direction is non-uniform and characterized by a denser spatial sampling in the vicinity of interfaces, which allows for a numerically stable and accurate evaluation of higher order surface wave modes. We stretch the grid in the vertical direction to increase the minimum grid spacing and reduce the computational cost. The free-surface boundary conditions are implemented with a characteristics approach, where the characteristic variables are evaluated at zero viscosity. The same procedure is used to model seismic wave propagation at the interface between a fluid and porous medium. In this case, each medium is represented by a different grid and the two grids are combined through a domain-decomposition method. This wavefield decomposition method accounts for the discontinuity of variables and is crucial for an accurate interface treatment. We simulate seismic wave propagation with open-pore and sealed-pore boundary conditions and verify the validity and accuracy of the algorithm by comparing the numerical simulations to analytical solutions based on zero viscosity obtained with the Cagniard-de Hoop method. Finally, we illustrate the suitability of our algorithm for more complex models of porous media involving viscous pore fluids and strongly heterogeneous distributions of the elastic and hydraulic material properties.

Acurácia e calibração de sonda de capacitância em Latossolo Vermelho cultivado com cafeeiro

O objetivo deste trabalho foi determinar a acurácia da sonda de multisensores de capacitância "Delta-T Profile probe PR2/6", na avaliação do conteúdo de água do solo com uso de calibrações padrão do fabricante, realizar a calibração para condições específicas de locais e profundidades de amostragem do solo e obter coeficientes de calibração para medições acuradas em tempo real. Em janeiro de 2010, foram coletadas amostras de solo com estrutura preservada a diferentes profundidades, nas linhas de plantio do cafeeiro e nas entrelinhas. As análises foram realizadas em laboratório, com o sensor ML2x Theta probe. Após a obtenção das leituras do sensor, o teor de água foi determinado por meio do método gravimétrico. Foram utilizadas amostras de Latossolo Vermelho distrófico muito argiloso. As calibrações padrão do fabricante (mineral e orgânica) não não se mostraram adequadas para emprego nas condições de manejo (locais e profundidades de amostragem) avaliadas. Na impossibilidade de averiguar a acurácia obtida pelo método recomendado pelo fabricante, o uso de ajustes de regressão linear ou da ferramenta Solver mostrou-se útil no processo de calibração. São necessárias apenas duas equações de calibração para avaliação do teor de água das situações contrastantes de manejo.

Computational and Experimental Analysis of Flow Field in the Diffusers of Centrifugal Compressors

Centrifugal compressors are widely used for example in process industry, oil and gas industry, in small gas turbines and turbochargers. In order to achieve lower consumption of energy and operation costs the efficiency of the compressor needs to be improve. In the present work different pinches and low solidity vaned diffusers were utilized in order to improve the efficiency of a medium size centrifugal compressor. In this study, pinch means the decrement of the diffuser flow passage height. First different geometries were analyzed using computational fluid dynamics. The flow solver Finflo was used to solve the flow field. Finflo is a Navier-Stokes solver. The solver is capable to solve compressible, incompressible, steady and unsteady flow fields. Chien's k-e turbulence model was used. One of the numerically investigated pinched diffuser and one low solidity vaned diffuser were studied experimentally. The overall performance of the compressor and the static pressure distribution before and after the diffuser were measured. The flow entering and leaving the diffuser was measured using a three-hole Cobra-probe and Kiel-probes. The pinch and the low solidity vaned diffuser increased the efficiency of the compressor. Highest isentropic efficiency increment obtained was 3\% of the design isentropic efficiency of the original geometry. It was noticed in the numerical results that the pinch made to the hub and the shroud wall was most beneficial to the operation of the compressor. Also the pinch made to the hub was better than the pinchmade to the shroud. The pinch did not affect the operation range of the compressor, but the low solidity vaned diffuser slightly decreased the operation range.The unsteady phenomena in the vaneless diffuser were studied experimentally andnumerically. The unsteady static pressure was measured at the diffuser inlet and outlet, and time-accurate numerical simulation was conducted. The unsteady static pressure showed that most of the pressure variations lay at the passing frequency of every second blade. The pressure variations did not vanish in the diffuser and were visible at the diffuser outlet. However, the amplitude of the pressure variations decreased in the diffuser. The time-accurate calculations showed quite a good agreement with the measured data. Agreement was very good at the design operation point, even though the computational grid was not dense enough inthe volute and in the exit cone. The time-accurate calculation over-predicted the amplitude of the pressure variations at high flow.

Description of a robotics-oriented relational positioning methodology

This paper presents a relational positioning methodology for flexibly and intuitively specifying offline programmed robot tasks, as well as for assisting the execution of teleoperated tasks demanding precise movements.In relational positioning, the movements of an object can be restricted totally or partially by specifying its allowed positions in terms of a set of geometric constraints. These allowed positions are found by means of a 3D sequential geometric constraint solver called PMF – Positioning Mobile with respect to Fixed. PMF exploits the fact that in a set of geometric constraints, the rotational component can often be separated from the translational one and solved independently.

Short-term Hydropower Planning System

The Thesis gives a decision support framework that has significant impact on the economic performance and viability of a hydropower company. The studyaddresses the short-term hydropower planning problem in the Nordic deregulated electricity market. The basics of the Nordic electricity market, trading mechanisms, hydropower system characteristics and production planning are presented in the Thesis. The related modelling theory and optimization methods are covered aswell. The Thesis provides a mixed integer linear programming model applied in asuccessive linearization method for optimal bidding and scheduling decisions inthe hydropower system operation within short-term horizon. A scenario based deterministic approach is exploited for modelling uncertainty in market price and inflow. The Thesis proposes a calibration framework to examine the physical accuracy and economic optimality of the decisions suggested by the model. A calibration example is provided with data from a real hydropower system using a commercial modelling application with the mixed integer linear programming solver CPLEX.

Sähkön omatuotannon kannattavuustarkastelu

Tässä diplomityössä tutkitaan sähkön omatuotannon kannattavuutta M-realin Simpeleen tehtaiden voimalaitoksella. Erityisesti työssä arvioidaan lauhdesähköntuotannon kannattavuutta polttoaine- ja päästökustannuksista muodostuvien marginaalikustannusten osalta. Koska voimalaitoksen rakennusaste on varsin alhainen,sähköntuotannon kannattavuutta on tarkasteltu arvioimalla sähkön ja lämmön yhteistuotannon kustannuksia ja jakamalla syntyneet kustannukset suhdemenetelmän avulla eri tuotteille. Diplomityössä etsitään kustannustehokkain seospolttosuhde annettujen reunaehtojen puitteissa muodostamalla polttoaineista aiheutuvista kustannuksista laskentamalli, jota optimoidaan Microsoft Excelin Solver-toiminnolla. Lauhdesähköntuotannon marginaalikustannuksia verrataan Nord Poolin SPOT-tuntihintaan. Lauhdesähköntuotanto voimalaitoksella on kannattavaa, mikäli SPOT-tuntihinnan vuorokautinen keskiarvo ylittää lauhdesähköntuotannon marginaalikustannukset.

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.

Experimental and Numerical Analysis of Different Volutes in a Centrifugal Compressor

In a centrifugal compressor the flow around the diffuser is collected and led to the pipe system by a spiral-shaped volute. In this study a single-stage centrifugal compressor with three different volutes is investigated. The compressorwas first equipped with the original volute, the cross-section of which was a combination of a rectangle and semi-circle. Next a new volute with a fully circular cross-section was designed and manufactured. Finally, the circular volute wasmodified by rounding the tongue and smoothing the tongue area. The overall performance of the compressor as well as the static pressure distribution after the impeller and on the volute surface were measured. The flow entering the volute was measured using a three-hole Cobra-probe, and flow visualisations were carriedout in the exit cone of the volute. In addition, the radial force acting on theimpeller was measured using magnetic bearings. The complete compressor with thecircular volute (inlet pipe, full impeller, diffuser, volute and outlet pipe) was also modelled using computational fluid dynamics (CFD). A fully 3-D viscous flow was solved using a Navier-Stokes solver, Finflo, developed at Helsinki University of Technology. Chien's k-e model was used to take account of the turbulence. The differences observed in the performance of the different volutes were quite small. The biggest differences were at low speeds and high volume flows,i.e. when the flow entered the volute most radially. In this operating regime the efficiency of the compressor with the modified circular volute was about two percentage points higher than with the other volutes. Also, according to the Cobra-probe measurements and flow visualisations, the modified circular volute performed better than the other volutes in this operating area. The circumferential static pressure distribution in the volute showed increases at low flow, constant distribution at the design flow and decrease at high flow. The non-uniform static pressure distribution of the volute was transmitted backwards across the vaneless diffuser and observed at the impeller exit. At low volume flow a strong two-wave pattern developed into the static pressure distribution at the impeller exit due to the response of the impeller to the non-uniformity of pressure. The radial force of the impeller was the greatest at the choke limit, the smallest atthe design flow, and moderate at low flow. At low flow the force increase was quite mild, whereas the increase at high flow was rapid. Thus, the non-uniformityof pressure and the force related to it are strong especially at high flow. Theforce caused by the modified circular volute was weaker at choke and more symmetric as a function of the volume flow than the force caused by the other volutes.

Kiinnirullaimen dynamiikan ja värähtelyjen simulointi

Diplomityössä tutkittiin kaupallisen monikappaledynamiikkaohjelmiston soveltuvuutta kiinnirullaimen dynamiikan ja värähtelyjen tutkimiseen. Erityisen kiinnostuneita oltiin nipin kuvauksesta sekä nipissä tapahtuvista värähtelyistä. Tässä diplomityössä mallinnettiin kiinnirullaimen ensiö- ja toisiokäytöt sekä tampuuritela. Malli yhdistettiin myöhemmin Metso Paper Järvenpäässä rinnakkaisena diplomityönä tehtyyn malliin, joista muodostui kahteen ratkaisijaan perustuva simulointimalli. Simulointimalli rakennettiin käyttämään kahta erillistä ratkaisijaa, joista toinen on mekaniikkamallin rakentamisessa käytetty ADAMS-ohjelmisto ja toinen säätöjärjestelmää ja hydraulipiirejä kuvaava Simulink-malli. Nipin mallintamiseksi tampuuritela ja rullaussylinteri mallinnettiin joustaviksi käyttäen keskitettyjen massojen menetelmää. Siirtolaitteissa sekä runkorakenteissa tapahtuvat joustot kuvattiin yhden vapausasteen jousi-vaimennin voimilla kuvattuina järjestelminä. Tässä diplomityössä on myös keskitytty esittelemään ADAMS-ohjelmiston toimintaa ohjeistavasti sekä käsittelemään parametrisen mallintamisen etuja. Työssä havaittiin monikappaledynamiikan soveltuvuus kiinnirullaimen dynamiikan sekä dynaamisten voimien aiheuttamien värähtelyjen tutkimiseen. Suoritetuista värähtelymittauksista voitiin tehdä vain arvioita. Mallin havaittiin vaativan lisätutkimusta ja kehitystyötä

Talonrakennushankkeen logististen kustannusten mallintaminen

Logistiikan kehittämisen merkitys on suuri rakennusalan yritysten kilpailukyvylle. Tutkimuksen kohdeyrityksessä ei tunnettu logistiikkakustannuksia ja tämän takia eri toimitusketjuvaihtoehtojen kustannusvaikutuksia ei voitu vertailla. Ilman kustannustietoja toiminnan kehittäminen ja kilpailukyky kärsivät. Rakentamisen kehittämisen nopein tehostamis- ja kustannussäästöpotentiaali onkin logistiikassa. Tutkimuksen tavoite oli luoda helppokäyttöinen, nopea, ymmärrettävä ja luotettavaa informaatiota tarjoava toimitusketjuratkaisija, jonka avulla voidaan löytää kustannustehokkain toimitusketju projektitietojen perusteella. Alustava mallinnus ja testaus suoritettiin Excel-ohjelmalla. Malli laski logistiset kustannukset kustannustekijäkohtaisesti yhteen huomioiden tutkimuksen rajaukset. Osa rakennusalan tuotteista aiheuttaa huomattavan määrän kokonais- sekä logistisista kustannuksista, kun taas tiettyjen tuotteiden kustannukset ovat vain pieniä. Eri talotyypeille voidaan löytää yleisiä ohjesääntöjä toimitusketjujen käytöstä, mutta työmailla on aina myös tuotteita, jotka kannattaa tuoda joko terminaali- tai suorana toimituksena. Työmaiden olosuhteet ja toimitusmäärät ovatkin varsin ratkaisevassa asemassa toimitustapaa valittaessa. Edellä mainittuja tutkimustuloksia tukivat alan kirjallisuus sekä mallin testaustulokset.

Lämmönsiirron numeerinen laskenta profiloidulle kanavalle.

Työn tavoitteena oli tutkia aaltomaisen profiloinnin vaikutusta suorien jäähdytyskanavien lämmönsiirtoon ja painehäviöön. Erilaisia profiileja oli kymmenen kappaletta ja ne olivat 5mm leveitä ja 30cm pitkiä kukin. Ne laskettiin kolmeulotteisina tapauksina FINFLO-virtausratkaisijalla kolmella eri Reynoldsin luvulla, jotka vastasivat laminaarista, osittain turbulenttista ja lähes kokonaan turbulenttista virtausta. Lämmönsiirtoaine oli kuiva +30°C ilma ja profiloinnin toteutustapa oli toisiaan sivuavat ympyräkaaret kolmella erilaisella säteen arvolla ja kolmella erilaisella aallonpituuden arvolla. Lisäksi laskettiin saman levyisen tasokanavan arvot jokaisella Reynoldsin luvulla kaksiulotteisina tapauksina. Näitä profiloimattomia kanavia pidettiin referenssitapauksina. Tuloksena havaittiin että profiloimalla saadaan yksiselitteisesti suurempi lämpöteho ulos samasta tilavuudesta. Lämmönsiirtokerroin kasvaa profiloinnin avulla parhaimmillaan n. 20% käytetystä turbulenssimallista tai lämmönsiirtokertoimen määritelmästä riippumatta. Painehäviö kasvaa myös aina, mutta kitkakerroin voi hieman pienentyä. Profiilin varsinaisena hyvyyskriteerinä pidettiin lämmönsiirtokertoimen ja kitkakertoimen suhdetta h/f. Se osoittautui riippuvaksi Reynoldsin luvusta ja turbulenssimallista; ASM ja Chien k-έ -mallit ennustavat transitioetäisyyden eri tavalla. Laminaarisilla virtauksilla h/f :n vaihtelu oli vähäistä; suhde vaihteli vain ±5% eri profiilien kesken. ASM-mallilla havaittiin sekundääripyörteilyä, ehkä siksi että se mallintaa anisotrooppisen turbulenssin. Chien k-έ malli ennusti suuremman ja aikaisemmin alkavan turbulenttisuuden kuin ASM. Lisäksi havaittiin mm. että tietyillä profiileilla muodostuu kanavan kapeimpaan kohtaan selvä virtausnopeuden paikallinen minimi seinämän läheisyyden takia.

Matrix Analytic Analysis of Delays in Communication Systems

Tässä päättötyössä annetaan kuvaus kehitetystä sovelluksesta Quasi Birth Death processien ratkaisuun. Tämä ohjelma on tähän mennessä ainutlaatuinen ja sen avulla voi ratkaista sarjan tehtäviä ja sitä tarvitaan kommunikaatio systeemien analyysiin. Mainittuun sovellukseen on annettu kuvaus ja määritelmä. Lyhyt kuvaus toisesta sovelluksesta Quasi Birth Death prosessien tehtävien ratkaisuun on myös annettu

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.

Oferta de trabajo y análisis impositivo con el ordenador: tutorial y actividad práctica

En el presente documento proponemos un ejercicio tutorial y una práctica docente que permita al estudiante a asimilar y entender mejor el concepto de exceso de gravamen o pérdida irrecuperable de eficiencia asociada a la imposición. Para ello planteamos, en primer lugar, un modelo de elección individual entre ocio y trabajo/consumo. Centrándonos en el individuo representativo, estudiamos su elección de horas de trabajo en ausencia de impuestos. A continuación, introducimos un impuesto proporcional sobre la renta del trabajo y, alternativamente, un impuesto de suma fija capaz de generar una recaudación igual. Verificamos que el impuesto proporcional siempre deja al individuo con un nivel de bienestar inferior respecto al escenario del impuesto de suma fija, es decir, que el impuesto proporcional genera una pérdida irrecuperable de eficiencia. Aunque la resolución manual del problema es factible, proponemos utilizar la herramienta Solver de Excel. Esto nos permite jugar de manera más ágil con los valores de los diferentes parámetros y verificar que hay pérdida de eficiencia independientemente de que el individuo decida trabajar más o menos horas en respuesta a la introducción del impuesto sobre la renta del trabajo