Development of Self-Cleaning Box Culvert Designs, TR-545, 2009

The main function of a roadway culvert is to effectively convey drainage flow during normal and extreme hydrologic conditions. This function is often impaired due to the sedimentation blockage of the culvert. This research sought to understand the mechanics of sedimentation process at multi-box culverts, and develop self-cleaning systems that flush out sediment deposits using the power of drainage flows. The research entailed field observations, laboratory experiments, and numerical simulations. The specific role of each of these investigative tools is summarized below: a) The field observations were aimed at understanding typical sedimentation patterns and their dependence on culvert geometry and hydrodynamic conditions during normal and extreme hydrologic events. b) The laboratory experiments were used for modeling sedimentation process observed insitu and for testing alternative self-cleaning concepts applied to culverts. The major tasks for the initial laboratory model study were to accurately replicate the culvert performance curves and the dynamics of sedimentation process, and to provide benchmark data for numerical simulation validation. c) The numerical simulations enhanced the understanding of the sedimentation processes and aided in testing flow cases complementary to those conducted in the model reducing the number of (more expensive) tests to be conducted in the laboratory. Using the findings acquired from the laboratory and simulation works, self-cleaning culvert concepts were developed and tested for a range of flow conditions. The screening of the alternative concepts was made through experimental studies in a 1:20 scale model guided by numerical simulations. To ensure the designs are effective, performance studies were finally conducted in a 1:20 hydraulic model using the most promising design alternatives to make sure that the proposed systems operate satisfactory under closer to natural scale conditions.

Analysis of Industrial Granular Flow Applications by Using Advanced Collision Models

Granular flow phenomena are frequently encountered in the design of process and industrial plants in the traditional fields of the chemical, nuclear and oil industries as well as in other activities such as food and materials handling. Multi-phase flow is one important branch of the granular flow. Granular materials have unusual kinds of behavior compared to normal materials, either solids or fluids. Although some of the characteristics are still not well-known yet, one thing is confirmed: the particle-particle interaction plays a key role in the dynamics of granular materials, especially for dense granular materials. At the beginning of this thesis, detailed illustration of developing two models for describing the interaction based on the results of finite-element simulation, dimension analysis and numerical simulation is presented. The first model is used to describing the normal collision of viscoelastic particles. Based on some existent models, more parameters are added to this model, which make the model predict the experimental results more accurately. The second model is used for oblique collision, which include the effects from tangential velocity, angular velocity and surface friction based on Coulomb's law. The theoretical predictions of this model are in agreement with those by finite-element simulation. I n the latter chapters of this thesis, the models are used to predict industrial granular flow and the agreement between the simulations and experiments also shows the validation of the new model. The first case presents the simulation of granular flow passing over a circular obstacle. The simulations successfully predict the existence of a parabolic steady layer and show how the characteristics of the particles, such as coefficients of restitution and surface friction affect the separation results. The second case is a spinning container filled with granular material. Employing the previous models, the simulation could also reproduce experimentally observed phenomena, such as a depression in the center of a high frequency rotation. The third application is about gas-solid mixed flow in a vertically vibrated device. Gas phase motion is added to coherence with the particle motion. The governing equations of the gas phase are solved by using the Large eddy simulation (LES) and particle motion is predicted by using the Lagrangian method. The simulation predicted some pattern formation reported by experiment.

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.

Todennäköisyyspohjaisen paloanalyysimenetelmän kehittäminen ydinvoimalaitoksella

Diplomityössä on tutustuttu ydinvoimalaitosten paloriskejä käsittelevään todennäköisyyspohjaiseen turvallisuusanalyysiin. Tavoitteena on ollut Olkiluoto 1 ja 2 laitosyksiköiden paloanalyysimenetelmän kehittäminen. Työssä esitetään paloanalyysin pääpiirteet, kaksi erilaista palotaajuuksien estimointimenetelmää sekä palojen leviämisen arviointimenetelmiä. Palotaajuuksien estimointimenetelmistä keskitytään Berryn menetelmän sekä NUREG/CR-6850-palotaajuuslaskentamenetelmän tarkasteluun. Palon leviämisen arvioinnissa on esitetty kolmen erilaisen virtausteknisen laskentatyökalun perusteet sekä palon leviämistodennäköisyyksiä arvioivan Probabilistic Fire Simulator (PFS) -ohjelman käyttöä. Työn aikana on laskettu molemmilla palotaajuuden estimointimenetelmillä palotaajuuksia eri tyyppisille huonetiloille. Berryn menetelmän palotaajuudet olivat pääosin alhaisempia kuin NUREG/CR-6850-menetelmällä lasketut palotaajuudet. Palon leviämistarkastelussa on tutkittu ydinvoimalaitoksen relehuoneen tulipaloa. PFS:n avulla laskettujen leviämistodennäköisyyksien arvoja on vertailtu TVO:n paloanalyysissa käytettyihin kvalitatiivisiin peittokertoimiin. Palon leviämistodennäköisyys eri osajärjestelmien välillä todettiin suuresti riippuvan analyysissaoletetuista vaurioitumislämpötiloista. Tutkittuja menetelmiä hyödyntäen diplomityössä kehitettiin paloanalyysimenetelmäkuvaus. Menetelmäkuvauksessa huonetilojen paloriskit kartoitetaan aluksi Berryn menetelmällä. Näin kaikille laitoksen huonetiloille saadaan arvioitua palotaajuus sekä paloalkutapahtumaluokkien sydänvauriotaajuus. Seuraavaksi suoritetaan valintamenettely, jossa valitut kriteerit täyttäville huonetiloille tehdään tarkentava palotaajuuslaskenta. Tarkentava palotaajuuslaskenta perustuu NUREG/CR-6850-menetelmän mukaisesti huonetilojen realistisiin syttymislähteisiin. Kriittisimpien huonetilojen osalta palon leviämisen arviointiin on tarkoitus hyödyntää numeerista simulointia.

Modelado, simulación y control digital de autogiros

Este proyecto consiste en diseñar el algoritmo de control de un autogiro no tripulado. Su aplicación principal es llevar a cabo tareas rutinarias o peligrosas para el piloto como, por ejemplo, extinción de incendios, evaluación de riesgo químico o vigilancia de lugares de acceso restringido. Se realiza un estudio del movimiento del vehículo para obtener su modelo dinámico. A partir de las ecuaciones que describen su movimiento, se realiza una simulación numérica del vehículo. Se incorpora el controlador diseñado y se evalúa su funcionamiento. Finalmente, se implementa el sistema en un microcontrolador.

Motion artifact reduction and vessel enhancement for free-breathing navigator-gated coronary MRA using 3D k-space reordering.

Breathing-induced bulk motion of the myocardium during data acquisition may cause severe image artifacts in coronary magnetic resonance angiography (MRA). Current motion compensation strategies include breath-holding or free-breathing MR navigator gating and tracking techniques. Navigator-based techniques have been further refined by the applications of sophisticated 2D k-space reordering techniques. A further improvement in image quality and a reduction of relative scanning duration may be expected from a 3D k-space reordering scheme. Therefore, a 3D k-space reordered acquisition scheme including a 3D navigator gated and corrected segmented k-space gradient echo imaging sequence for coronary MRA was implemented. This new zonal motion-adapted acquisition and reordering technique (ZMART) was developed on the basis of a numerical simulation of the Bloch equations. The technique was implemented on a commercial 1.5T MR system, and first phantom and in vivo experiments were performed. Consistent with the results of the theoretical findings, the results obtained in the phantom studies demonstrate a significant reduction of motion artifacts when compared to conventional (non-k-space reordered) gating techniques. Preliminary in vivo findings also compare favorably with the phantom experiments and theoretical considerations. Magn Reson Med 45:645-652, 2001.

Particle transport method for convection-diffusion-reaction problems

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

Modeling dynamic behavior in tilting pad gas journal bearings

Thisresearch deals with the dynamic modeling of gas lubricated tilting pad journal bearings provided with spring supported pads, including experimental verification of the computation. On the basis of a mathematical model of a film bearing, a computer program has been developed, which can be used for the simulation of a special type of tilting pad gas journal bearing supported by a rotary spring under different loading conditions time dependently (transient running conditions due to geometry variations in time externally imposed). On the basis of literature, different transformations have been used in the model to achieve simpler calculation. The numerical simulation is used to solve a non-stationary case of a gasfilm. The simulation results were compared with literature results in a stationary case (steady running conditions) and they were found to be equal. In addition to this, comparisons were made with a number of stationary and non-stationary bearing tests, which were performed at Lappeenranta University of Technology using bearings designed with the simulation program. A study was also made using numerical simulation and literature to establish the influence of the different bearing parameters on the stability of the bearing. Comparison work was done with literature on tilting pad gas bearings. This bearing type is rarely used. One literature reference has studied the same bearing type as that used in LUT. A new design of tilting pad gas bearing is introduced. It is based on a stainless steel body and electron beam welding of the bearing parts. It has good operation characteristics and is easier to tune and faster to manufacture than traditional constructions. It is also suitable for large serial production.

Application of the N-alkane molecular alloys to thermally protected containers for catering

A thermally controlled transport device was designed and tested. As hot food needs to be transported at temperatures between 60 and 70ºC in order to avoid contamination by microorganisms, the use of Molecular Alloy Phase Change Materials (MAPCM) can lead to improvements in this field of application. A heat transfer numerical simulation of the box used for transporting the food was conducted. Despite obvious simplifications, a good agreement between numerical simulation and experimental results was obtained. Furthermore, we compared our experimental results with those from other experiments related to the transport of hot food. Here, pizza is taken as the example, and it is shown that delivering time can be increased three-fold.

Puumassapesurin suotautumisolojen kartoitus ja simulointi

Työn tavoitteena oli tutkia sellumassan suotautumiseen liittyviä tekijöitä Ahlstrom Machineryn Drum Displacer™ -puumassapesurissa (DD-pesuri). Teoriaosassa tarkasteltiin aluksi suotautumisen teoriaa kuitu-vesi-suspensiossa, minkä jälkeen esiteltiin suotautumisnopeuteen vaikuttavia fysikaalisia ja kemiallisia vaikutusmekanismeja. Seuraavaksi kuvattiin massan pesun yleisiä perusteita sekä teoriaa puumassapesureissa. Lopuksi tarkasteltiin pesurien kytkeytymistä muuhun kuitulinjaan sekä prosessista johtuvia pesun toiminnan ulkoisia häiriötekijöitä. Kokeellisen osan aluksi tarkasteltiin paine- ja lämpötilamittauksien avulla massapesurissa vallitsevia prosessioloja. Mittaustulosten perusteella pumppausolot pesurin suodoslinjoissa ovat vaikeahkot ja häiriötilanteita voi esiintyä, mutta käytäntö on osoittanut tästä olevan vain harvoin haittaa prosessin toiminnalle. Pesureissa toteutuneet syrjäytysnopeudet laskettiin ja niitä verrattiin syrjäytystestien antamiin tuloksiin. Kuitulinjasta riippuen testin vastaavuus tehdasprosessiin vaihteli suuresti. Syrjäytystesteillä kokeiltiin myös tehdasprosesseissa usein esiintyvien muuttujien vaikutusta sellukakun syrjäytettävyyteen. Kakun paksuus ja syrjäytyslämpötila vaikuttivat syrjäytysnopeuteen Darcyn lain mukaisesti. Alipaineen massakakun alapuolella havaittiin huonontavan syrjäytysnopeutta verrattuna tilanteeseen, jossa kakun alla vallitsi ilmanpaine. Tämä havainto on selvästi ristiriidassa suotautumisen teorian kanssa. Massakakun muodostumis-pH osoittautui ratkaisevaksi lopulliselle syrjäytysnopeudelle, sillä alkalisissa oloissa muodostetun kuitukakun syrjäytysnopeus ei enää parantunut happamalla syrjäytysnesteellä. Happamissa oloissa muodostetun kakun syrjäytysnopeus oli alkalisista parempi, mutta se alkoi hitaasti alentua, kun syrjäytysneste vaihtui alkaliseen. Massan laimentaminen ennen syrjäytystä alkalisella tehdassuodoksella puhtaan veden sijasta alensi ligniinipitoisella massalla lopullista syrjäytysnopeutta. Shirato-Tillerin mallilla ja Jönssonin staattisella mallilla simuloitiin numeerisesti syrjäytystestiä kahdessa eri pH:ssa, ja simulointituloksia verrattiin vastaavissa oloissa tehtyihin syrjäytystesteihin. Shirato-Tillerin mallin antamien syrjäytysnopeuksien havaittiin olevan lähellä syrjäytystestien nopeuksia, kun Jönssonin mallin antamat tulokset jäivät huomattavasti testituloksia alemmiksi. Herkkyystarkastelussa havaittiin mallien olevan varsin herkkiä parametrien virheille. Hajonta vaadittavien kuituparametrien määrityksissä ja menetelmien työläys rajoittavat numeerisen simuloinnin käytettävyyttä, sillä kuituparametrien määrityksen vaatima työmäärä on ainakin toistaiseksi syrjäytystestiä suurempi. Lopuksi todettiin, että oikeiden syrjäytysolosuhteiden käyttö on ensiarvoisen tärkeää oikeiden tulosten saamiseksi sekä kokeellisessa että numeerisessa simuloinnissa. Nykyinen syrjäytystestilaitteisto on pienin muutoksin käyttökelpoinen, kun massan testaus prosessioloissa tulee rutiininomaiseksi.

Synthesis of highly focused fields with circular polarization at any transverse plane

We develop a method for generating focused vector beams with circular polarization at any transverse plane. Based on the Richards-Wolf vector model, we derive analytical expressions to describe the propagation of these set of beams near the focal area. Since the polarization and the amplitude of the input beam are not uniform, an interferometric system capable of generating spatially-variant polarized beams has to be used. In particular, this wavefront is manipulated by means of spatial light modulators displaying computer generated holograms and subsequently focused using a high numerical aperture objective lens. Experimental results using a NA=0.85 system are provided: irradiance and Stokes images of the focused field at different planes near the focal plane are presented and compared with those obtained by numerical simulation.

Propagação da poluição atmosférica por flúor nas águas subterrâneas e solos de regiões próximas às indústrias de fertilizantes (Rio Grande, RS)

Fluoride concentration was determined in rainwater, ground water and soil in the zone of fertilizer industry in the city of Rio Grande. In contaminated rainwater fluorine concentration was registered up to the value of 4,4 mg.L-1. Fluorine concentration in the shallow ground water in general reflects its distribution in the atmosphere, but cannot be used as marker of atmosphere contamination in the urban area due to dissolving influence of residential effluents. The 0,01% HCl extracts from the set of surface soil samples demonstrates fluoride distribution in the zone of influence of industrial emissions, which coincides to the numerical simulation of fluorides dispersion in the air.

Chaotic dynamics of electric-field domains in periodically driven superlattices

Self-sustained time-dependent current oscillations under dc voltage bias have been observed in recent experiments on n-doped semiconductor superlattices with sequential resonant tunneling. The current oscillations are caused by the motion and recycling of the domain wall separating low- and high-electric-field regions of the superlattice, as the analysis of a discrete drift model shows and experimental evidence supports. Numerical simulation shows that different nonlinear dynamical regimes of the domain wall appear when an external microwave signal is superimposed on the dc bias and its driving frequency and driving amplitude vary. On the frequency-amplitude parameter plane, there are regions of entrainment and quasiperiodicity forming Arnold tongues. Chaos is demonstrated to appear at the boundaries of the tongues and in the regions where they overlap. Coexistence of up to four electric-field domains randomly nucleated in space is detected under ac+dc driving.

Application of a non isotropic turbulence model to stable atmospheric flows and dispersion over 3D topography

A non isotropic turbulence model is extended and applied to three dimensional stably stratified flows and dispersion calculations. The model is derived from the algebraic stress model (including wall proximity effects), but it retains the simplicity of the "eddy viscosity" concept of first order models. The "modified k-epsilon" is implemented in a three dimensional numerical code. Once the flow is resolved, the predicted velocity and turbulence fields are interpolated into a second grid and used to solve the concentration equation. To evaluate the model, various steady state numerical solutions are compared with small scale dispersion experiments which were conducted at the wind tunnel of Mitsubishi Heavy Industries, in Japan. Stably stratified flows and plume dispersion over three distinct idealized complex topographies (flat and hilly terrain) are studied. Vertical profiles of velocity and pollutant concentration are shown and discussed. Also, comparisons are made against the results obtained with the standard k-epsilon model.

Thermalhydraulic behavior of electrically heated rod during a critical heat flux transient

In nuclear reactors, the occurrence of critical heat flux leads to fuel rod overheating with clad fusion and radioactive products leakage. To predict the effects of such phenomenon, experiments are performed using electrically heated rods to simulate operational and accidental conditions of nuclear fuel rods. In the present work, it is performed a theoretical analysis of the drying and rewetting front propagation during a critical heat flux experiment, starting with the application of an electrical power step from steady state condition. After the occurrence of critical heat flux, the drying front propagation is predicted. After a few seconds, a power cut is considered and the rewetting front behavior is analytically observed. Studies performed with various values of coolant mass flow rate show that this variable has more influence on the drying front velocity than on the rewetting one.