Kuumalanka-anemometri virtausmittauksissa

Tässä työssä selostetaan kuumalanka-anemometrin käyttö virtausmittauksissa. Kuumalanka-anemometrilla saadaan mitattua virtausnopeuden ja -suunnan lisäksi nopeusheilahteluja. Mittaustaajuus on tyypillisesti useita kymmeniä tuhansia mittauksia sekunnissa ja signaali on jatkuva. Nykytekniikalla pystytään helposti tallentamaan mittauslaitteistolta saatu viesti tietokoneelle ja muuntamaan se nopeudeksi. Hetkellisten nopeuksien avulla voidaan laskea turbulenttisen virtauksen ominaisuuksia, kuten turbulenssin intensiteetti ja spektri. Kuumalanka-anemometrissa lämmitetään sähköisesti ohutta lankaa, joka on mitattavassa virtauksessa. Langan sähköteho on suunnilleen yhtäsuuri kuin langasta konvektiolla siirtyvä lämpöteho. Tällöin on teoreettisesti mahdollista laskea virtausnopeus lämpötehosta lämmönsiirtokorrelaatioilla. Käytännössä laitteisto joudutaan kuitenkin erikseen kalibroimaan, mutta sähkötehon teoreettista riippuvuutta konvektiosta käytetään hyväksi. Kuumalangan lämmitettävä osuus on tyypillisesti halkaisijaltaan 5 µm ja pituudeltaan noin 1 mm. Sitä käytetään pääasiassa kaasuvirtausten mittaamiseen ja valtaosassa mittauksissa virtausaineena on ilma. Kuumalanka voi olla toteutettu kuumakalvotekniikalla, jossa halkaisijaltaan noin 50 - 70 µm paksuinen kuitu on päällystetty ohuella sähköä johtavalla kalvolla. Kuumakalvoanturin ei tarvitse olla muodoltaan sylinterimäinen, se voi olla mm. kartiomainen tai kiilamainen. Erikoispäällystetyllä kuumakalvoanturilla on mahdollista mitata myös nestevirtauksia. Mitattaessa kaasuvirtauksia kuumakalvon etuna on selvästi parempi kestävyys verrattuna kuumalankaan. Nimitystä kuumalanka-anemometri käytetään yleisesti molemmista anturityypeistä Tämän työn alussa käsitellään sylinterin yli tapahtuvaan virtaukseen liittyvää virtausmekaniikkaa ja lämmönsiirtoa. Anemometrin sähköinen osa, laitteisto ja sen kalibrointi käydään läpi. Langan suuntariippuvuuden laskentaan esitetään tarvittavat yhtälöt. Työssä esitellään kolme laitteistolla tehtyä perusmittausta: anturin kohtauskulman muuttaminen, pyörähdyssymmerisen suihkun nopeuskenttä ja tuulitunnelin rajakerros. Lisäksi esitellään yksi käytännöllinen ja vaativampi mittaus, jossa on mitattu nopeusprofiili radiaalikompressorin diffuusorin loppuosassa.

Komponenttien lämpöteknisten parametrien karakterisointilaitteiston suunnittelu ja toteutus

Elektroniikan pienentyessä ja tehotiheyden kasvaessa tarvitaan paremmin suunniteltuja jäähdytysjärjestelyjä. Suunnittelun apuna voidaan käyttää perinteisiä lämmönsiirtoyhtälöitä, virtauslaskentaa ja mittauksia. Tässä työssä on suunniteltu ja rakennettu mittauslaitteisto nimenomaan elektroniikan komponenttien ja jäähdytyselementtien lämpötilojen ja painehäviöiden mittaamiseen. Koeajona laitteistolla on mitattu jäähdytyselementin lämpötilajakauma kahdella eri virtausnopeudella ja lämpötilalla sekä verrattu niitä aikaisemmin mitattuihin tapauksiin sekä Fluent- virtauslaskentaohjelmalla laskettuihin tuloksiin.

ADVANCED NUMERICAL METHODS FOR SIMULATING MICROBUBBLE IN AN AERATED MIXING TANK

Fluid mixing in mechanically agitated tanks is one of the major unit operations in many industries. Bubbly flows have been of interest among researchers in physics, medicine, chemistry and technology over the centuries. The aim of this thesis is to use advanced numerical methods for simulating microbubble in an aerated mixing tank. Main components of the mixing tank are a cylindrical vessel, a rotating Rushton turbine and the air nozzle. The objective of Computational Fluid Dynamics (CFD) is to predict fluid flow, heat transfer, mass transfer and chemical reactions. The CFD simulations of a turbulent bubbly flow are carried out in a cylindrical mixing tank using large eddy simulation (LES) and volume of fluid (VOF) method. The Rushton turbine induced flow is modeled by using a sliding mesh method. Numerical results are used to describe the bubbly flows in highly complex liquid flow. Some of the experimental works related to turbulent bubbly flow in a mixing tank are briefly reported. Numerical simulations are needed to complete and interpret the results of the experimental work. Information given by numerical simulations has a major role in designing and scaling-up mixing tanks. The results of this work have been reported in the following scientific articles: ·Honkanen M., Koohestany A., Hatunen T., Saarenrinne P., Zamankhan P., Large eddy simulations and PIV experiments of a two-phase air-water mixer, in Proceedings of ASME Fluids Engineering Summer Conference (2005). ·Honkanen M., Koohestany A., Hatunen T., Saarenrinne P., Zamankhan P., Dynamical States of Bubbling in an Aerated Stirring Tank, submitted to J. Computational Physics.

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.

Virtausjakauman kokeellinen määritys levylämmönsiirtimissä

Diplomityön tavoitteena oli kokeellisen tutkimuksen keinoin selvittää juotettujen levylämmönsiirtimien levypakkarakenteessa virtausten käyttäytyminen ja jakautuminen sekä löytää ideoita ja kehitysehdotuksia levylämmönsiirtimen levypakan ja levyprofiilin kehittämiseksi. Kokeellinen tutkimus suoritettiin Oy Danfoss Ab LPM:n levylämmönsiirtimien tutkimuslaboratoriossa. Virtausjakauman tutkimusta varten suunniteltiin ja valittiin tutkimuslaitteisto, joka koostui termoelementtiantureista, tiedonkeruulaitteistosta sekä ohjelmistosta. Lämmönsiirtimistä mitattiin ensiö- ja toisiopuolen tilavuusvirrat ja painehäviöt sekä lämpötilat ennen ja jälkeen lämmönsiirtimen. Tutkimuslaitteiston avulla mitattiin lämpötiloja lämmönsiirtimen sisältä levyväleistä. Mittaukset suoritettiin neljällä levypakkarakenteella useilla massavirran arvoilla. Mittaustuloksista määritettiin levylämmönsiirtimien lämpö- ja virtaustekniset ominaisuudet nesteen Reynoldsin luvun funktiona sekä selvitettiin nesteen virtausjakaumat. Mittaustuloksien perusteella laskettuja virtausjakauman arvoja verrattiin teorian mukaan laskettuihin jakaumiin. Mitatuista siirtimistä lasketut massavirrat viittaavat siihen, että suurin osa nesteestä virtaa siirtimien keskeltä tai lähempää loppupäätä kuin alkupäästä. Teorian mukaan suurin nestemäärä virtaisi siirtimen alkupäästä vähentyen tasaisesti kohti levypakan loppupäätä. Teorian mukaiselle virtausjakaumalle ja lasketuille jakaumille ei löydetty yhteyttä. Tutkimuksessa havaittiin suuria, jopa yli 20 asteen, lämpötilaeroja levyväleistä ulostulevissa virtauksissa. Levyvälien virtauksen käyttäytymisen ja jakautumisen tutkiminen nähdäänkin levypakan pitkittäistä kehittämistä suurempana mielenkiinnon ja kehittämisen kohteena.

Effect of hydrodynamics on modelling, monitoring and control of crystallization

Crystallization is a purification method used to obtain crystalline product of a certain crystal size. It is one of the oldest industrial unit processes and commonly used in modern industry due to its good purification capability from rather impure solutions with reasonably low energy consumption. However, the process is extremely challenging to model and control because it involves inhomogeneous mixing and many simultaneous phenomena such as nucleation, crystal growth and agglomeration. All these phenomena are dependent on supersaturation, i.e. the difference between actual liquid phase concentration and solubility. Homogeneous mass and heat transfer in the crystallizer would greatly simplify modelling and control of crystallization processes, such conditions are, however, not the reality, especially in industrial scale processes. Consequently, the hydrodynamics of crystallizers, i.e. the combination of mixing, feed and product removal flows, and recycling of the suspension, needs to be thoroughly investigated. Understanding of hydrodynamics is important in crystallization, especially inlargerscale equipment where uniform flow conditions are difficult to attain. It is also important to understand different size scales of mixing; micro-, meso- and macromixing. Fast processes, like nucleation and chemical reactions, are typically highly dependent on micro- and mesomixing but macromixing, which equalizes the concentrations of all the species within the entire crystallizer, cannot be disregarded. This study investigates the influence of hydrodynamics on crystallization processes. Modelling of crystallizers with the mixed suspension mixed product removal (MSMPR) theory (ideal mixing), computational fluid dynamics (CFD), and a compartmental multiblock model is compared. The importance of proper verification of CFD and multiblock models is demonstrated. In addition, the influence of different hydrodynamic conditions on reactive crystallization process control is studied. Finally, the effect of extreme local supersaturation is studied using power ultrasound to initiate nucleation. The present work shows that mixing and chemical feeding conditions clearly affect induction time and cluster formation, nucleation, growth kinetics, and agglomeration. Consequently, the properties of crystalline end products, e.g. crystal size and crystal habit, can be influenced by management of mixing and feeding conditions. Impurities may have varying impacts on crystallization processes. As an example, manganese ions were shown to replace magnesium ions in the crystal lattice of magnesium sulphate heptahydrate, increasing the crystal growth rate significantly, whereas sodium ions showed no interaction at all. Modelling of continuous crystallization based on MSMPR theory showed that the model is feasible in a small laboratoryscale crystallizer, whereas in larger pilot- and industrial-scale crystallizers hydrodynamic effects should be taken into account. For that reason, CFD and multiblock modelling are shown to be effective tools for modelling crystallization with inhomogeneous mixing. The present work shows also that selection of the measurement point, or points in the case of multiprobe systems, is crucial when process analytical technology (PAT) is used to control larger scale crystallization. The thesis concludes by describing how control of local supersaturation by highly localized ultrasound was successfully applied to induce nucleation and to control polymorphism in reactive crystallization of L-glutamic acid.

Determination of slurry’s viscosity using case based reasoning approach

Case-based reasoning (CBR) is a recent approach to problem solving and learning that has got a lot of attention over the last years. In this work, the CBR methodology is used to reduce the time and amount of resources spent on carry out experiments to determine the viscosity of the new slurry. The aim of this work is: to develop a CBR system to support the decision making process about the type of slurries behavior, to collect a sufficient volume of qualitative data for case base, and to calculate the viscosity of the Newtonian slurries. Firstly in this paper, the literature review about the types of fluid flow, Newtonian and non-Newtonian slurries is presented. Some physical properties of the suspensions are also considered. The second part of the literature review provides an overview of the case-based reasoning field. Different models and stages of CBR cycles, benefits and disadvantages of this methodology are considered subsequently. Brief review of the CBS tools is also given in this work. Finally, some results of work and opportunities for system modernization are presented. To develop a decision support system for slurry viscosity determination, software application MS Office Excel was used. Designed system consists of three parts: workspace, the case base, and section for calculating the viscosity of Newtonian slurries. First and second sections are supposed to work with Newtonian and Bingham fluids. In the last section, apparent viscosity can be calculated for Newtonian slurries.

Utilization of Particle Image velocimetry in PPOOLEX condensation experiments

Particle Image Velocimetry, PIV, is an optical measuring technique to obtain velocity information of a flow in interest. With PIV it is possible to achieve two or three dimensional velocity vector fields from a measurement area instead of a single point in a flow. Measured flow can be either in liquid or in gas form. PIV is nowadays widely applied to flow field studies. The need for PIV is to obtain validation data for Computational Fluid Dynamics calculation programs that has been used to model blow down experiments in PPOOLEX test facility in the Lappeenranta University of Technology. In this thesis PIV and its theoretical background are presented. All the subsystems that can be considered to be part of a PIV system are presented as well with detail. Emphasis is also put to the mathematics behind the image evaluation. The work also included selection and successful testing of a PIV system, as well as the planning of the installation to the PPOOLEX facility. Already in the preliminary testing PIV was found to be good addition to the measuring equipment for Nuclear Safety Research Unit of LUT. The installation to PPOOLEX facility was successful even though there were many restrictions considering it. All parts of the PIV system worked and they were found out to be appropriate for the planned use. Results and observations presented in this thesis are a good background to further PIV use.

Energy-efficient control strategies for variable speed driven parallel pumping systems based on pump operation point monitoring with frequency converters

The pumping processes requiring wide range of flow are often equipped with parallelconnected centrifugal pumps. In parallel pumping systems, the use of variable speed control allows that the required output for the process can be delivered with a varying number of operated pump units and selected rotational speed references. However, the optimization of the parallel-connected rotational speed controlled pump units often requires adaptive modelling of both parallel pump characteristics and the surrounding system in varying operation conditions. The available information required for the system modelling in typical parallel pumping applications such as waste water treatment and various cooling and water delivery pumping tasks can be limited, and the lack of real-time operation point monitoring often sets limits for accurate energy efficiency optimization. Hence, alternatives for easily implementable control strategies which can be adopted with minimum system data are necessary. This doctoral thesis concentrates on the methods that allow the energy efficient use of variable speed controlled parallel pumps in system scenarios in which the parallel pump units consist of a centrifugal pump, an electric motor, and a frequency converter. Firstly, the suitable operation conditions for variable speed controlled parallel pumps are studied. Secondly, methods for determining the output of each parallel pump unit using characteristic curve-based operation point estimation with frequency converter are discussed. Thirdly, the implementation of the control strategy based on real-time pump operation point estimation and sub-optimization of each parallel pump unit is studied. The findings of the thesis support the idea that the energy efficiency of the pumping can be increased without the installation of new, more efficient components in the systems by simply adopting suitable control strategies. An easily implementable and adaptive control strategy for variable speed controlled parallel pumping systems can be created by utilizing the pump operation point estimation available in modern frequency converters. Hence, additional real-time flow metering, start-up measurements, and detailed system model are unnecessary, and the pumping task can be fulfilled by determining a speed reference for each parallel-pump unit which suggests the energy efficient operation of the pumping system.

Modelling of ultrasound assisted mixing in Newtonian and non-Newtonian liquids

The purpose of this work is to obtain a better understanding of behaviour of possible ultrasound appliance on fluid media mixing. The research is done in the regard to Newtonian and non-Newtonian fluids. The process of ultrasound appliance on liquids is modelled in COMSOL Multiphysics software. The influence of ultrasound using is introduced as waveform equation. Turbulence modelling is fulfilled by the k-ε model in Newtonian fluid. The modeling of ultrasound assisted mixing in non-Newtonian fluids is based on the power law. To verify modelling results two practical methods are used: Particle Image Velocimetry and measurements of mixing time. Particle Image Velocimetry allows capturing of velocity flow field continuously and presents detailed depiction of liquid dynamics. The second way of verification is the comparison of mixing time of homogeneity. Experimentally achievement of mixing time is done by conductivity measurements. In modelling part mixing time is achieved by special module of COMSOL Multiphysics – the transport of diluted species. Both practical and modelling parts show similar radial mechanism of fluid flow under ultrasound appliance – from the horn tip fluid moves to the bottom and along the walls goes back. Velocity profiles are similar in modelling and experimental part in the case of Newtonian fluid. In the case of non-Newtonian fluid velocity profiles do not agree. The development track of ultrasound-assisted mixing modelling is presented in the thesis.

Korkeapainevesisumusuuttimen virtaustekninen laadunhallinta

Tämä työ vastaa tarpeeseen hallita korkeapainevesisumusuuttimen laatua virtausmekaniikan työkalujen avulla. Työssä tutkitaan suutinten testidatan lisäksi virtauksen käyttäytymistä suuttimen sisällä CFD-laskennan avulla. Virtausmallinnus tehdään Navier-Stokes –pohjaisella laskentamenetelmällä. Työn teoriaosassa käsitellään virtaustekniikkaa ja sen kehitystä yleisesti. Lisäksi esitetään suuttimen laskennassa käytettävää perusteoriaa sekä teknisiä ratkaisuja. Teoriaosassa käydään myös läpi laskennalliseen virtausmekaniikkaan (CFD-laskenta) liittyvää perusteoriaa. Tutkimusosiossa esitetään käsitellyt suutintestitulokset sekä mallinnetaan suutinvirtausta ajasta riippumattomaan virtauslaskentaan perustuvalla laskentamenetelmällä. Virtauslaskennassa käytetään OpenFOAM-laskentaohjelmiston SIMPLE-virtausratkaisijaa sekä k-omega SST –turbulenssimallia. Tehtiin virtausmallinnus kaikilla paineilla, joita suuttimen testauksessa myös todellisuudessa käytetään. Lisäksi selvitettiin mahdolliset kavitaatiokohdat suuttimessa ja suunniteltiin kavitaatiota ehkäisevä suutingeometria. Todettiin myös lämpötilan ja epäpuhtauksien vaikuttavan kavitaatioon sekä mallinnettiin lämpötilan vaikutusta. Luotiin malli, jolla suuttimen suunnitteluun liittyviin haasteisiin voidaan vastata numeerisella laskennalla.

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.

Static waves in corporate space : characterizing oscillating trading patterns in New York stock exchange

Various researches in the field of econophysics has shown that fluid flow have analogous phenomena in financial market behavior, the typical parallelism being delivered between energy in fluids and information on markets. However, the geometry of the manifold on which market dynamics act out their dynamics (corporate space) is not yet known. In this thesis, utilizing a Seven year time series of prices of stocks used to compute S&P500 index on the New York Stock Exchange, we have created local chart to the corporate space with the goal of finding standing waves and other soliton like patterns in the behavior of stock price deviations from the S&P500 index. By first calculating the correlation matrix of normalized stock price deviations from the S&P500 index, we have performed a local singular value decomposition over a set of four different time windows as guides to the nature of patterns that may emerge. I turns out that in almost all cases, each singular vector is essentially determined by relatively small set of companies with big positive or negative weights on that singular vector. Over particular time windows, sometimes these weights are strongly correlated with at least one industrial sector and certain sectors are more prone to fast dynamics whereas others have longer standing waves.

Pain-evoked alterations on gingival blood flow and gingival crevicular fluid (GCF) neuropeptide SP and collagenase-2 (MMP-8) levels

Previous studies have demonstrated that clinical pulpal pain can induce the expression of pro-inflammatory neuropeptides in the adjacent gingival crevice fluid (GCF). Vasoactive agents such as substance P (SP) are known to contribute to the inflammatory type of pain and are associated with increased blood flow. More recent animal studies have shown that application of capsaicin on alveolar mucosa provokes pain and neurogenic vasodilatation in the adjacent gingiva. Pain-associated inflammatory reactions may initiate expression of several pro- and anti-inflammatory mediators. Collagenase-2 (MMP-8) has been considered to be the major destructive protease, especially in the periodontitis-affected gingival crevice fluid (GCF). MMP-8 originates mostly from neutrophil leukocytes, the first line of defence cells that exist abundantly in GCF, especially in inflammation. With this background, we wished to clarify the spatial extensions and differences between tooth-pain stimulation and capsaicin-induced neurogenic vasodilatation in human gingiva. Experiments were carried out to study whether tooth stimulation and capsaicin stimulation of alveolar mucosa would induce changes in GCF MMP-8 levels and whether tooth stimulation would release neuropeptide SP in GCF. The experiments were carried out on healthy human volunteers. During the experiments, moderate and high intensity painful tooth stimulation was performed by a constant current tooth stimulator. Moderate tooth stimulation activates A-delta fibres, while high stimulation also activates C-fibres. Painful stimulation of the gingiva was achieved by topical application of capsaicin-moistened filter paper on the mucosal surface. Capsaicin is known to activate selectively nociceptive C-fibres of stimulated tissue. Pain-evoked vasoactive changes in gingivomucosal tissues were mapped by laser Doppler imaging (LDI), which is a sophisticated and non-invasive method for studying e.g. spatial and temporal characteristics of pain- and inflammation-evoked blood flow changes in gingivomucosal tissues. Pain-evoked release of MMP-8 in GCF samples was studied by immunofluorometric assay (IFMA) and Western immunoblotting. The SP levels in GCF were analysed by Enzyme immunoassay (EIA). During the experiments, subjective stimulus-evoked pain responses were determined by a visual analogue pain scale. Unilateral stimulation of alveolar mucosa and attached gingiva by capsaicin evoked a distinct neurogenic vasodilatation in the ipsilateral gingiva, which attenuated rapidly at the midline. Capsaicin stimulation of alveolar mucosa provoked clear inflammatory reactions. In contrast to capsaicin stimuli, tooth stimulation produced symmetrical vasodilatations bilaterally in the gingiva. The ipsilateral responses were significantly smaller during tooth stimulation than during capsaicin stimuli. The current finding – that tooth stimulation evokes bilateral vasodilatation while capsaicin stimulation of the gingiva mainly produces unilateral vasodilatation – emphasises the usefulness of LDI in clarifying spatial features of neurogenic vasoactive changes in the intra-oral tissues. Capsaicin stimulation of the alveolar mucosa induced significant elevations in MMP-8 levels and activation in GCF of the adjacent teeth. During the experiments, no marked changes occurred in MMP-8 levels in the GCF of distantly located teeth. Painful stimulation of the upper incisor provoked elevations in GCF MMP-8 and SP levels of the stimulated tooth. The GCF MMP-8 and SP levels of the non-stimulated teeth were not changed. These results suggest that capsaicin-induced inflammatory reactions in gingivomucosal tissues do not cross the midline in the anterior maxilla. The enhanced reaction found during stimulation of alveolar mucosa indicates that alveolar mucosa is more sensitive to chemical irritants than the attached gingiva. Analysis of these data suggests that capsaicin-evoked neurogenic inflammation in the gingiva can trigger the expression and activation of MMP-8 in GCF of the adjacent teeth. In this study, it is concluded that experimental tooth pain at C-fibre intensity can induce local elevations in MMP-8 and SP levels in GCF. Depending on the role of MMP-8 in inflammation, in addition to surrogated tissue destruction, the elevated MMP-8 in GCF may also reflect accelerated local defensive and anti-inflammatory reactions.

Application of Computational Fluid Dynamics in Modeling Blood Flow in Human Thoracic Aorta

The aim of this study was to simulate blood flow in thoracic human aorta and understand the role of flow dynamics in the initialization and localization of atherosclerotic plaque in human thoracic aorta. The blood flow dynamics in idealized and realistic models of human thoracic aorta were numerically simulated in three idealized and two realistic thoracic aorta models. The idealized models of thoracic aorta were reconstructed with measurements available from literature, and the realistic models of thoracic aorta were constructed by image processing Computed Tomographic (CT) images. The CT images were made available by South Karelia Central Hospital in Lappeenranta. The reconstruction of thoracic aorta consisted of operations, such as contrast adjustment, image segmentations, and 3D surface rendering. Additional design operations were performed to make the aorta model compatible for the numerical method based computer code. The image processing and design operations were performed with specialized medical image processing software. Pulsatile pressure and velocity boundary conditions were deployed as inlet boundary conditions. The blood flow was assumed homogeneous and incompressible. The blood was assumed to be a Newtonian fluid. The simulations with idealized models of thoracic aorta were carried out with Finite Element Method based computer code, while the simulations with realistic models of thoracic aorta were carried out with Finite Volume Method based computer code. Simulations were carried out for four cardiac cycles. The distribution of flow, pressure and Wall Shear Stress (WSS) observed during the fourth cardiac cycle were extensively analyzed. The aim of carrying out the simulations with idealized model was to get an estimate of flow dynamics in a realistic aorta model. The motive behind the choice of three aorta models with distinct features was to understand the dependence of flow dynamics on aorta anatomy. Highly disturbed and nonuniform distribution of velocity and WSS was observed in aortic arch, near brachiocephalic, left common artery, and left subclavian artery. On the other hand, the WSS profiles at the roots of branches show significant differences with geometry variation of aorta and branches. The comparison of instantaneous WSS profiles revealed that the model with straight branching arteries had relatively lower WSS compared to that in the aorta model with curved branches. In addition to this, significant differences were observed in the spatial and temporal profiles of WSS, flow, and pressure. The study with idealized model was extended to study blood flow in thoracic aorta under the effects of hypertension and hypotension. One of the idealized aorta models was modified along with the boundary conditions to mimic the thoracic aorta under the effects of hypertension and hypotension. The results of simulations with realistic models extracted from CT scans demonstrated more realistic flow dynamics than that in the idealized models. During systole, the velocity in ascending aorta was skewed towards the outer wall of aortic arch. The flow develops secondary flow patterns as it moves downstream towards aortic arch. Unlike idealized models, the distribution of flow was nonplanar and heavily guided by the artery anatomy. Flow cavitation was observed in the aorta model which was imaged giving longer branches. This could not be properly observed in the model with imaging containing a shorter length for aortic branches. The flow circulation was also observed in the inner wall of the aortic arch. However, during the diastole, the flow profiles were almost flat and regular due the acceleration of flow at the inlet. The flow profiles were weakly turbulent during the flow reversal. The complex flow patterns caused a non-uniform distribution of WSS. High WSS was distributed at the junction of branches and aortic arch. Low WSS was distributed at the proximal part of the junction, while intermedium WSS was distributed in the distal part of the junction. The pulsatile nature of the inflow caused oscillating WSS at the branch entry region and inner curvature of aortic arch. Based on the WSS distribution in the realistic model, one of the aorta models was altered to induce artificial atherosclerotic plaque at the branch entry region and inner curvature of aortic arch. Atherosclerotic plaque causing 50% blockage of lumen was introduced in brachiocephalic artery, common carotid artery, left subclavian artery, and aortic arch. The aim of this part of the study was first to study the effect of stenosis on flow and WSS distribution, understand the effect of shape of atherosclerotic plaque on flow and WSS distribution, and finally to investigate the effect of lumen blockage severity on flow and WSS distributions. The results revealed that the distribution of WSS is significantly affected by plaque with mere 50% stenosis. The asymmetric shape of stenosis causes higher WSS in branching arteries than in the cases with symmetric plaque. The flow dynamics within thoracic aorta models has been extensively studied and reported here. The effects of pressure and arterial anatomy on the flow dynamic were investigated. The distribution of complex flow and WSS is correlated with the localization of atherosclerosis. With the available results we can conclude that the thoracic aorta, with complex anatomy is the most vulnerable artery for the localization and development of atherosclerosis. The flow dynamics and arterial anatomy play a role in the localization of atherosclerosis. The patient specific image based models can be used to diagnose the locations in the aorta vulnerable to the development of arterial diseases such as atherosclerosis.