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.

An alternative finite element formulation for determination of streamlines in two-dimensional problems

It is well known that the numerical solutions of incompressible viscous flows are of great importance in Fluid Dynamics. The graphics output capabilities of their computational codes have revolutionized the communication of ideas to the non-specialist public. In general those codes include, in their hydrodynamic features, the visualization of flow streamlines - essentially a form of contour plot showing the line patterns of the flow - and the magnitudes and orientations of their velocity vectors. However, the standard finite element formulation to compute streamlines suffers from the disadvantage of requiring the determination of boundary integrals, leading to cumbersome implementations at the construction of the finite element code. In this article, we introduce an efficient way - via an alternative variational formulation - to determine the streamlines for fluid flows, which does not need the computation of contour integrals. In order to illustrate the good performance of the alternative formulation proposed, we capture the streamlines of three viscous models: Stokes, Navier-Stokes and Viscoelastic flows.

Experimental analysis of pressure drop and flow redistribution in axial flows in rod bundles

Fuel elements of PWR type nuclear reactors consist of rod bundles, arranged in a square array, and held by spacer grids. The coolant flows, mainly, axially along the rods. Although such elements are laterally open, experiments are performed in closed type test sections, originating the appearance of subchannels with different geometries. In the present work, utilizing a test section of two bundles of 4x4 pins each, experiments were performed to determine the friction and the grid drag coefficients for the different subchannels and to observe the effect of the grids in the crossflow, in cases of inlet flow maldistribution.

Flows of Bingham Materials Through Ideal Porous Media: an Experimental and Theoretical Study

The flow of Bingham liquids through porous media has been studied. Experiments have been performed to determine the flow rate / pressure drop relationship for the flow of a grease of Binghamian rheological behavior through an array of rods of circular cross section. The yield stress and plastic viscosity of the grease have been determined with the aid of a controlled stress rotational rheometer. To investigate a wider range of the flow parameters, the mass and momentum conservation equations have been solved numerically, in conjunction with the generalized Newtonian constitutive law and the bi-viscosity model. The finite volume method has been employed to obtain the numerical solution. These numerical results also yielded a flow rate / pressure drop relationship, which is in very good agreement with the experimental results. A capillaric theory has been developed to determine an analytical relationship between the flow rate and pressure drop for flows of Bingham liquids through porous media. It is shown that the predictions of this theory are in good agreement with the experimental and numerical results.

Wind tunnel simulation of atmospheric boundary layer flows

The present work shows how thick boundary layers can be produced in a short wind tunnel with a view to simulate atmospheric flows. Several types of thickening devices are analysed. The experimental assessment of the devices was conducted by considering integral properties of the flow and the spectra: skin-friction, mean velocity profiles in inner and outer co-ordinates and longitudinal turbulence. Designs based on screens, elliptic wedge generators, and cylindrical rod generators are analysed. The paper describes in detail the experimental arrangement, including the features of the wind tunnel and of the instrumentation. The results are compared with experimental data published by other authors and with naturally developed flows.

Roles of leadership in complex environments : Enhancing Knowledge Flows in Organisational Constellations through Practice-Based Innovation Processes

Global challenges, complexity and continuous uncertainty demand development of leadership approaches, employees and multi-organisation constellations. Current leadership theories do not sufficiently address the needs of complex business environments. First of all, before successful leadership models can be applied in practice, leadership needs to shift from the industrial age to the knowledge era. Many leadership models still view leadership solely through the perspective of linear process thinking. In addition, there is not enough knowledge or experience in applying these newer models in practice. Leadership theories continue to be based on the assumption that leaders possess or have access to all the relevant knowledge and capabilities to decide future directions without external advice. In many companies, however, the workforce consists of skilled professionals whose work and related interfaces are so challenging that the leaders cannot grasp all the linked viewpoints and cross-impacts alone. One of the main objectives of this study is to understand how to support participants in organisations and their stakeholders to, through practice-based innovation processes, confront various environments. Another aim is to find effective ways of recognising and reacting to diverse contexts, so companies and other stakeholders are better able to link to knowledge flows and shared value creation processes in advancing joint value to their customers. The main research question of this dissertation is, then, to seek understanding of how to enhance leadership in complex environments. The dissertation can, on the whole, be characterised as a qualitative multiple-case study. The research questions and objectives were investigated through six studies published in international scientific journals. The main methods applied were interviews, action research and a survey. The empirical focus was on Finnish companies, and the research questions were examined in various organisations at the top levels (leaders and managers) and bottom levels (employees) in the context of collaboration between organisations and cooperation between case companies and their client organisations. However, the emphasis of the analysis is the internal and external aspects of organisations, which are conducted in practice-based innovation processes. The results of this study suggest that the Cynefin framework, complexity leadership theory and transformational leadership represent theoretical models applicable to developing leadership through practice-based innovation. In and of themselves, they all support confronting contemporary challenges, but an implementable method for organisations may be constructed by assimilating them into practice-based innovation processes. Recognition of diverse environments, their various contexts and roles in the activities and collaboration of organisations and their interest groups is ever-more important to achieving better interaction in which a strategic or formal status may be bypassed. In innovation processes, it is not necessarily the leader who is in possession of the essential knowledge; thus, it is the role of leadership to offer methods and arenas where different actors may generate advances. Enabling and supporting continuous interaction and integrated knowledge flows is of crucial importance, to achieve emergence of innovations in the activities of organisations and various forms of collaboration. The main contribution of this dissertation relates to applying these new conceptual models in practice. Empirical evidence on the relevance of different leadership roles in practice-based innovation processes in Finnish companies is another valuable contribution. Finally, the dissertation sheds light on the significance of combining complexity science with leadership and innovation theories in research.

Data Flows and Engaging Visualisations for Your Repositories

Poster at Open Repositories 2014, Helsinki, Finland, June 9-13, 2014

Slurry flows in progressive cavity and in hose pumps

The aim of this thesis was to research how slurry’s viscosity and rheology affect to pumping in peristaltic hose pump and in eccentric progressive cavity pump. In addition, it was researched the formed pressure pulsation in hose pump. Pressure pulsation was studied by pumping different slurries and by using different pipe materials. Pressure and power curves were determined for both used pumps. It was also determined NPSHR curve for the progressive cavity pump. Literature part of the thesis considered to distribute fluids to different rheology types, as well as theories and models to identify different rheology types. Special attention was paid to non-Newtonian fluids, which were also used in experimental part of this thesis. In addition, the literature part discusses about pumps, parameters for pump sizing, and pressure pulsation in hose pump. Starch, bentonite, and carboxymethyl cellulose slurries were used in the experimental part of this thesis. The slurries were pumped with Flowrox peristaltic hose pump (LPP-T32) and eccentric progressive cavity pump (C10/10). From the each slurry was taken a sample, and the samples were analyzed for concentration, viscosity and rheology type. The used pipe materials in pressure pulsation experiments were steel and elastic, and it was also used a prototype of pulsation dampener. The pulsation experiments indicated that the elastic pipe and the prototype of pulsation dampener attenuated pressure pulsation better than the steel pipe at low pressure levels. The differences between different materials disappeared when pressure level and pump rotation speed increased. In slurry experiments, pulsation was different depending on rheology and viscosity of the slurry. According to experiments, the rheology did not significantly affect to pump power consumption or efficiency.

Large-eddy simulation of wind flows over complex terrains for wind energy applications

Wind energy has obtained outstanding expectations due to risks of global warming and nuclear energy production plant accidents. Nowadays, wind farms are often constructed in areas of complex terrain. A potential wind farm location must have the site thoroughly surveyed and the wind climatology analyzed before installing any hardware. Therefore, modeling of Atmospheric Boundary Layer (ABL) flows over complex terrains containing, e.g. hills, forest, and lakes is of great interest in wind energy applications, as it can help in locating and optimizing the wind farms. Numerical modeling of wind flows using Computational Fluid Dynamics (CFD) has become a popular technique during the last few decades. Due to the inherent flow variability and large-scale unsteadiness typical in ABL flows in general and especially over complex terrains, the flow can be difficult to be predicted accurately enough by using the Reynolds-Averaged Navier-Stokes equations (RANS). Large- Eddy Simulation (LES) resolves the largest and thus most important turbulent eddies and models only the small-scale motions which are more universal than the large eddies and thus easier to model. Therefore, LES is expected to be more suitable for this kind of simulations although it is computationally more expensive than the RANS approach. With the fast development of computers and open-source CFD software during the recent years, the application of LES toward atmospheric flow is becoming increasingly common nowadays. The aim of the work is to simulate atmospheric flows over realistic and complex terrains by means of LES. Evaluation of potential in-land wind park locations will be the main application for these simulations. Development of the LES methodology to simulate the atmospheric flows over realistic terrains is reported in the thesis. The work also aims at validating the LES methodology at a real scale. In the thesis, LES are carried out for flow problems ranging from basic channel flows to real atmospheric flows over one of the most recent real-life complex terrain problems, the Bolund hill. All the simulations reported in the thesis are carried out using a new OpenFOAM® -based LES solver. The solver uses the 4th order time-accurate Runge-Kutta scheme and a fractional step method. Moreover, development of the LES methodology includes special attention to two boundary conditions: the upstream (inflow) and wall boundary conditions. The upstream boundary condition is generated by using the so-called recycling technique, in which the instantaneous flow properties are sampled on aplane downstream of the inlet and mapped back to the inlet at each time step. This technique develops the upstream boundary-layer flow together with the inflow turbulence without using any precursor simulation and thus within a single computational domain. The roughness of the terrain surface is modeled by implementing a new wall function into OpenFOAM® during the thesis work. Both, the recycling method and the newly implemented wall function, are validated for the channel flows at relatively high Reynolds number before applying them to the atmospheric flow applications. After validating the LES model over simple flows, the simulations are carried out for atmospheric boundary-layer flows over two types of hills: first, two-dimensional wind-tunnel hill profiles and second, the Bolund hill located in Roskilde Fjord, Denmark. For the twodimensional wind-tunnel hills, the study focuses on the overall flow behavior as a function of the hill slope. Moreover, the simulations are repeated using another wall function suitable for smooth surfaces, which already existed in OpenFOAM® , in order to study the sensitivity of the flow to the surface roughness in ABL flows. The simulated results obtained using the two wall functions are compared against the wind-tunnel measurements. It is shown that LES using the implemented wall function produces overall satisfactory results on the turbulent flow over the two-dimensional hills. The prediction of the flow separation and reattachment-length for the steeper hill is closer to the measurements than the other numerical studies reported in the past for the same hill geometry. The field measurement campaign performed over the Bolund hill provides the most recent field-experiment dataset for the mean flow and the turbulence properties. A number of research groups have simulated the wind flows over the Bolund hill. Due to the challenging features of the hill such as the almost vertical hill slope, it is considered as an ideal experimental test case for validating micro-scale CFD models for wind energy applications. In this work, the simulated results obtained for two wind directions are compared against the field measurements. It is shown that the present LES can reproduce the complex turbulent wind flow structures over a complicated terrain such as the Bolund hill. Especially, the present LES results show the best prediction of the turbulent kinetic energy with an average error of 24.1%, which is a 43% smaller than any other model results reported in the past for the Bolund case. Finally, the validated LES methodology is demonstrated to simulate the wind flow over the existing Muukko wind farm located in South-Eastern Finland. The simulation is carried out only for one wind direction and the results on the instantaneous and time-averaged wind speeds are briefly reported. The demonstration case is followed by discussions on the practical aspects of LES for the wind resource assessment over a realistic inland wind farm.

Helsinki hub – new traffic flows and business networks

This thesis is the Logistics Development Forum's assignment and the work dealing with the development of the Port of Helsinki as part of Helsinki hub. The Forum aims to develop logistics efficiency through public-private co-operation and development of the port is clearly dependent on both factors. Freight volumes in the Port of Helsinki are the biggest single factor in hub and, therefore, the role of the port of the entire hub development is strong. The aim is to look at how the port will develop as a result of changes in the foreign trade of Finland and the Northern European logistics trends in 25 years time period. Work includes the current state analysis and scenario work. The analyses are intended to find out, which trends are the most important in the port volume development. The change and effect of trends is examined through scenarios based on current state. Based on the work, the structure of Finnish export industry and international demand are in the key role in the port volume development. There is significant difference between demands of Finnish exporting products in different export markets and the development between the markets has different impacts on the port volumes by mass and cargo type. On the other hand, the Finnish economy is stuck in a prolonged recession and competition between ports has become a significant factor in the individual port's volume development. Ecological valuesand regulations have changed the competitive landscape and maritime transport emissions reductions has become an important competitive factor for short routes in the Baltic Sea, such as in the link between Helsinki and Tallinn.

Trade Creation, Trade Diversion and Investment Flows: Evidence from the Canada-Chile Free Trade Agreement.

Rapport de recherche

La méthode IIM pour une membrane immergée dans un fluide incompressible

La méthode IIM (Immersed Interface Method) permet d'étendre certaines méthodes numériques à des problèmes présentant des discontinuités. Elle est utilisée ici pour étudier un fluide incompressible régi par les équations de Navier-Stokes, dans lequel est immergée une membrane exerçant une force singulière. Nous utilisons une méthode de projection dans une grille de différences finies de type MAC. Une dérivation très complète des conditions de saut dans le cas où la viscosité est continue est présentée en annexe. Deux exemples numériques sont présentés : l'un sans membrane, et l'un où la membrane est immobile. Le cas général d'une membrane mobile est aussi étudié en profondeur.

Une nouvelle méthode smoothed particle hydrodynamics : simulation des interfaces immergées et de la dynamique Brownienne des molécules avec des interactions hydrodynamiques

Dans cette thèse, nous présentons une nouvelle méthode smoothed particle hydrodynamics (SPH) pour la résolution des équations de Navier-Stokes incompressibles, même en présence des forces singulières. Les termes de sources singulières sont traités d'une manière similaire à celle que l'on retrouve dans la méthode Immersed Boundary (IB) de Peskin (2002) ou de la méthode régularisée de Stokeslets (Cortez, 2001). Dans notre schéma numérique, nous mettons en oeuvre une méthode de projection sans pression de second ordre inspirée de Kim et Moin (1985). Ce schéma évite complètement les difficultés qui peuvent être rencontrées avec la prescription des conditions aux frontières de Neumann sur la pression. Nous présentons deux variantes de cette approche: l'une, Lagrangienne, qui est communément utilisée et l'autre, Eulerienne, car nous considérons simplement que les particules SPH sont des points de quadrature où les propriétés du fluide sont calculées, donc, ces points peuvent être laissés fixes dans le temps. Notre méthode SPH est d'abord testée à la résolution du problème de Poiseuille bidimensionnel entre deux plaques infinies et nous effectuons une analyse détaillée de l'erreur des calculs. Pour ce problème, les résultats sont similaires autant lorsque les particules SPH sont libres de se déplacer que lorsqu'elles sont fixes. Nous traitons, par ailleurs, du problème de la dynamique d'une membrane immergée dans un fluide visqueux et incompressible avec notre méthode SPH. La membrane est représentée par une spline cubique le long de laquelle la tension présente dans la membrane est calculée et transmise au fluide environnant. Les équations de Navier-Stokes, avec une force singulière issue de la membrane sont ensuite résolues pour déterminer la vitesse du fluide dans lequel est immergée la membrane. La vitesse du fluide, ainsi obtenue, est interpolée sur l'interface, afin de déterminer son déplacement. Nous discutons des avantages à maintenir les particules SPH fixes au lieu de les laisser libres de se déplacer. Nous appliquons ensuite notre méthode SPH à la simulation des écoulements confinés des solutions de polymères non dilués avec une interaction hydrodynamique et des forces d'exclusion de volume. Le point de départ de l'algorithme est le système couplé des équations de Langevin pour les polymères et le solvant (CLEPS) (voir par exemple Oono et Freed (1981) et Öttinger et Rabin (1989)) décrivant, dans le cas présent, les dynamiques microscopiques d'une solution de polymère en écoulement avec une représentation bille-ressort des macromolécules. Des tests numériques de certains écoulements dans des canaux bidimensionnels révèlent que l'utilisation de la méthode de projection d'ordre deux couplée à des points de quadrature SPH fixes conduit à un ordre de convergence de la vitesse qui est de deux et à une convergence d'ordre sensiblement égale à deux pour la pression, pourvu que la solution soit suffisamment lisse. Dans le cas des calculs à grandes échelles pour les altères et pour les chaînes de bille-ressort, un choix approprié du nombre de particules SPH en fonction du nombre des billes N permet, en l'absence des forces d'exclusion de volume, de montrer que le coût de notre algorithme est d'ordre O(N). Enfin, nous amorçons des calculs tridimensionnels avec notre modèle SPH. Dans cette optique, nous résolvons le problème de l'écoulement de Poiseuille tridimensionnel entre deux plaques parallèles infinies et le problème de l'écoulement de Poiseuille dans une conduite rectangulaire infiniment longue. De plus, nous simulons en dimension trois des écoulements confinés entre deux plaques infinies des solutions de polymères non diluées avec une interaction hydrodynamique et des forces d'exclusion de volume.