Fluid flow and dispersion predictions in porous media by utilizing advanced numerical methods

Diplomityön tavoitteena oli tarkastella numeerisen virtauslaskennan avulla virtaukseen liittyviä ilmiöitä ja kaasun dispersiota. Diplomityön sisältö on jaettu viiteen osaan; johdantoon, teoriaan, katsaukseen virtauksen mallinnukseen huokoisessa materiaalissa liittyviin tutkimusselvityksiin, numeeriseen mallinnukseen sekä tulosten esittämiseen ja johtopäätöksiin. Diplomityön alussa kiinnitettiin huomiota erilaisiin kokeellisiin, numeerisiin ja teoreettisiin mallinnusmenetelmiin, joilla voidaan mallintaa virtausta huokoisessa materiaalissa. Kirjallisuusosassa tehtiin katsaus aikaisemmin julkaistuihin puoliempiirisiin ja empiirisiin tutkimusselvityksiin, jotka liittyvät huokoisen materiaalin aiheuttamaan painehäviöön. Numeerisessa virtauslaskenta osassa rakennettiin ja esitettiin huokoista materiaalia kuvaavat numeeriset mallit käyttäen kaupallista FLUENT -ohjelmistoa. Työn lopussa arvioitiin teorian, numeerisen virtauslaskennan ja kokeellisten tutkimusselvitysten tuloksia. Kolmiulotteisen huokoisen materiaalinnumeerisessa mallinnuksesta saadut tulokset vaikuttivat lupaavilta. Näiden tulosten perusteella tehtiin suosituksia ajatellen tulevaa virtauksen mallinnusta huokoisessa materiaalissa. Osa tässä diplomityössä esitetyistä tuloksista tullaan esittämään 55. Kanadan Kemiantekniikan konferenssissa Torontossa 1619 Lokakuussa 2005. ASME :n kansainvälisessä tekniikan alan julkaisussa. Työ on hyväksytty esitettäväksi esitettäväksi laskennallisen virtausmekaniikan (CFD) aihealueessa 'Peruskäsitteet'. Lisäksi työn yksityiskohtaiset tulokset tullaan lähettämään myös CES:n julkaisuun.

Calculations of Boiling Two-Phase Flow Using a Porous Media Model

Boiling two-phase flow and the equations governing the motion of fluid in two-phase flows are discussed in this thesis. Disposition of the governing equations in three-dimensional complex geometries is considered from the perspective of the porous medium concept. The equations governing motion in two-phase flows were formulated, discretized and implemented in a subroutine for pressure-velocity solution utilizing the SIMPLE algorithm modified for two-phase flow. The subroutine was included in PORFLO, which is a three-dimensional 5-equation porous media model developed at VTT by Jaakko Miettinen. The development of two-phase flow and the resulting void fraction distribution was predicted in a geometry resembling a section of BWR fuel bundle in a couple of test cases using PORFLO.

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.

Flow in porous media of oil-water systems

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

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.

Coolability of porous core debris beds : Effects of bed geometry and multi-dimensional flooding

This thesis addresses the coolability of porous debris beds in the context of severe accident management of nuclear power reactors. In a hypothetical severe accident at a Nordic-type boiling water reactor, the lower drywell of the containment is flooded, for the purpose of cooling the core melt discharged from the reactor pressure vessel in a water pool. The melt is fragmented and solidified in the pool, ultimately forming a porous debris bed that generates decay heat. The properties of the bed determine the limiting value for the heat flux that can be removed from the debris to the surrounding water without the risk of re-melting. The coolability of porous debris beds has been investigated experimentally by measuring the dryout power in electrically heated test beds that have different geometries. The geometries represent the debris bed shapes that may form in an accident scenario. The focus is especially on heap-like, realistic geometries which facilitate the multi-dimensional infiltration (flooding) of coolant into the bed. Spherical and irregular particles have been used to simulate the debris. The experiments have been modeled using 2D and 3D simulation codes applicable to fluid flow and heat transfer in porous media. Based on the experimental and simulation results, an interpretation of the dryout behavior in complex debris bed geometries is presented, and the validity of the codes and models for dryout predictions is evaluated. According to the experimental and simulation results, the coolability of the debris bed depends on both the flooding mode and the height of the bed. In the experiments, it was found that multi-dimensional flooding increases the dryout heat flux and coolability in a heap-shaped debris bed by 47–58% compared to the dryout heat flux of a classical, top-flooded bed of the same height. However, heap-like beds are higher than flat, top-flooded beds, which results in the formation of larger steam flux at the top of the bed. This counteracts the effect of the multi-dimensional flooding. Based on the measured dryout heat fluxes, the maximum height of a heap-like bed can only be about 1.5 times the height of a top-flooded, cylindrical bed in order to preserve the direct benefit from the multi-dimensional flooding. In addition, studies were conducted to evaluate the hydrodynamically representative effective particle diameter, which is applied in simulation models to describe debris beds that consist of irregular particles with considerable size variation. The results suggest that the effective diameter is small, closest to the mean diameter based on the number or length of particles.

Liquid transportation and distribution during (re)wetting : impact on coated papers in printing

In many industrial applications, such as the printing and coatings industry, wetting of porous materials by liquids includes not only imbibition and permeation into the bulk but also surface spreading and evaporation. By understanding these phenomena, valuable information can be obtained for improved process control, runnability and printability, in which liquid penetration and subsequent drying play important quality and economic roles. Knowledge of the position of the wetting front and the distribution/degree of pore filling within the structure is crucial in describing the transport phenomena involved. Although exemplifying paper as a porous medium in this work, the generalisation to dynamic liquid transfer onto a surface, including permeation and imbibition into porous media, is of importance to many industrial and naturally occurring environmental processes. This thesis explains the phenomena in the field of heatset web offset printing but the content and the analyses are applicable in many other printing methods and also other technologies where water/moisture monitoring is crucial in order to have a stable process and achieve high quality end products. The use of near-infrared technology to study the water and moisture response of porous pigmented structures is presented. The use of sensitive surface chemical and structural analysis, as well as the internal structure investigation of a porous structure, to inspect liquid wetting and distribution, complements the information obtained by spectroscopic techniques. Strong emphasis has been put on the scale of measurement, to filter irrelevant information and to understand the relationship between interactions involved. The near-infrared spectroscopic technique, presented here, samples directly the changes in signal absorbance and its variation in the process at multiple locations in a print production line. The in-line non-contact measurements are facilitated by using several diffuse reflectance probes, giving the absolute water/moisture content from a defined position in the dynamic process in real-time. The nearinfrared measurement data illustrate the changes in moisture content as the paper is passing through the printing nips and dryer, respectively, and the analysis of the mechanisms involved highlight the roles of the contacting surfaces and the relative liquid carrier properties of both non-image and printed image areas. The thesis includes laboratory studies on wetting of porous media in the form of coated paper and compressed pigment tablets by mono-, dual-, and multi-component liquids, and paper water/moisture content analysis in both offline and online conditions, thus also enabling direct sampling of temporal water/moisture profiles from multiple locations. One main focus in this thesis was to establish a measurement system which is able to monitor rapid changes in moisture content of paper. The study suggests that near-infrared diffuse reflectance spectroscopy can be used as a moisture sensitive system and to provide accurate online qualitative indicators, but, also, when accurately calibrated, can provide quantification of water/moisture levels, its distribution and dynamic liquid transfer. Due to the high sensitivity, samples can be measured with excellent reproducibility and good signal to noise ratio. Another focus of this thesis was on the evolution of the moisture content, i.e. changes in moisture content referred to (re)wetting, and liquid distribution during printing of coated paper. The study confirmed different wetting phases together with the factors affecting each phase both for a single droplet and a liquid film applied on a porous substrate. For a single droplet, initial capillary driven imbibition is followed by equilibrium pore filling and liquid retreat by evaporation. In the case of a liquid film applied on paper, the controlling factors defining the transportation were concluded to be the applied liquid volume in relation to surface roughness, capillarity and permeability of the coating giving the liquid uptake capacity. The printing trials confirmed moisture gradients in the printed sheet depending on process parameters such as speed, fountain solution dosage and drying conditions as well as the printed layout itself. Uneven moisture distribution in the printed sheet was identified to be one of the sources for waving appearance and the magnitude of waving was influenced by the drying conditions.

Näyttöpäätetyön kuormituksen tasaaminen media-assistenttiopiskelijoilla

Opinnäytetyömme tarkoituksena oli antaa Omnia ammattikoulun media-assistenttiopiskelijoille tietoa näyttöpäätetyöhön liittyvistä kuormitustekijöistä ja ergonomiasta. Opinnäytetyömme tavoitteena on vastata kysymykseen " Mitä media-assistenttiopiskelijoiden tulisi tietää, jotta he osaisivat kuormittaa itseään sopivasti näyttöpäätetyössä?". Koostimme opinnäytetyömme kirjallisuudesta ja tutkimuksista saatujen tietojen pohjalta. Näistä muodostui opinnäytetyömme teoreettinen viitekehys sekä perustelut opinnäytetyömme tuotoksen sisällölle. Tutkimusten mukaan tuki- ja liikuntaelimistön kuormituksen tasaamiseen voi vaikuttaa vähentämällä haitallisesti kuormittavia työjaksoja, vaihtamalla työasentoja, säätämällä työvälineitä ja tauottamalla työntekoa. Opinnäytetyömme kirjallisessa osiossa perehdyimme näyttöpäätetyön ergonomiaan ja kuormittavuuteen media-assistenttiopiskelijoilla. Opinnäytetyömme tuotoksena oli "Kevennä kuormitusta näyttöpäätteen äärellä"- CD-Rom, jossa on 39 diaa. CD-Rom sisältää tietoa niska-hartiaseudun ja selän rakenteesta sekä toiminnasta, lihastyötavoista, niskaa ja selkää kuormittavista tekijöistä näyttöpäätteen äärellä työskennellessä, työpisteen säätelystä ja optimaalisesta istuma-asennosta, sekä vinkkejä kuormituksen keventämiseen mm. tauottamisen avulla. Testasimme CD-Rom:n sisältöä yhdellä Omnia ammattikoulun media-assistenttiopiskelijaryhmällä järjestämällä ergonomiaohjauksen. Palautteen pohjalta kehittelimme työn lopulliseen muotoonsa. Opinnäytetyömme tuotosta käytetään Omnia ammattikoulun media-assistenttiopiskelijoiden ergonomiaohjauksen välineenä. Työmme soveltuu myös muiden näyttöpäätteen äärellä työskentelevien henkilöiden ja ryhmien ergonomiaohjauksen välineeksi. "Kevennä kuormitusta näyttöpäätteen äärellä" - CD-Romia voivat hyödyntää fysioterapeutit, muut ergonomiaohjauksia antavat toimihenkilöt sekä Kunto-Stadia.

Compaction of growth media based on Sphagnum peat during one-year culturing of container seedlings

Tiivistelmä: Turvepohjaisten kasvualustojen tiivistyminen yksivuotisessa paakkutaimikasvatuksessa