Weakly linked binary mixtures of F = 1 87Rb Bose-Einstein condensates

We present a study of binary mixtures of Bose-Einstein condensates confined in a double-well potential within the framework of the mean field Gross-Pitaevskii (GP) equation. We re-examine both the single component and the binary mixture cases for such a potential, and we investigate what are the situations in which a simpler two-mode approach leads to an accurate description of their dynamics. We also estimate the validity of the most usual dimensionality reductions used to solve the GP equations. To this end, we compare both the semi-analytical two-mode approaches and the numerical simulations of the one-dimensional (1D) reductions with the full 3D numerical solutions of the GP equation. Our analysis provides a guide to clarify the validity of several simplified models that describe mean-field nonlinear dynamics, using an experimentally feasible binary mixture of an F = 1 spinor condensate with two of its Zeeman manifolds populated, m = ±1.

Direct numerical simulations of interface dynamics to link capillary pressure and total surface energy

The flow of two immiscible fluids through a porous medium depends on the complex interplay between gravity, capillarity, and viscous forces. The interaction between these forces and the geometry of the medium gives rise to a variety of complex flow regimes that are difficult to describe using continuum models. Although a number of pore-scale models have been employed, a careful investigation of the macroscopic effects of pore-scale processes requires methods based on conservation principles in order to reduce the number of modeling assumptions. In this work we perform direct numerical simulations of drainage by solving Navier-Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and model the transition from stable flow to viscous fingering, we focus on the macroscopic capillary pressure and we compare different definitions of this quantity under quasi-static and dynamic conditions. We show that the difference between the intrinsic phase-average pressures, which is commonly used as definition of Darcy-scale capillary pressure, is subject to several limitations and it is not accurate in presence of viscous effects or trapping. In contrast, a definition based on the variation of the total surface energy provides an accurate estimate of the macroscopic capillary pressure. This definition, which links the capillary pressure to its physical origin, allows a better separation of viscous effects and does not depend on the presence of trapped fluid clusters.

Mesoscale numerical analysis of the historical November 1982 heavy precipitation event over Andorra (Eastern Pyrenees)

From 6 to 8 November 1982 one of the most catastrophic flash-flood events was recorded in the Eastern Pyrenees affecting Andorra and also France and Spain with rainfall accumulations exceeding 400 mm in 24 h, 44 fatalities and widespread damage. This paper aims to exhaustively document this heavy precipitation event and examines mesoscale simulations performed by the French Meso-NH non-hydrostatic atmospheric model. Large-scale simulations show the slow-evolving synoptic environment favourable for the development of a deep Atlantic cyclone which induced a strong southerly flow over the Eastern Pyrenees. From the evolution of the synoptic pattern four distinct phases have been identified during the event. The mesoscale analysis presents the second and the third phase as the most intense in terms of rainfall accumulations and highlights the interaction of the moist and conditionally unstable flows with the mountains. The presence of a SW low level jet (30 m s-1) around 1500 m also had a crucial role on focusing the precipitation over the exposed south slopes of the Eastern Pyrenees. Backward trajectories based on Eulerian on-line passive tracers indicate that the orographic uplift was the main forcing mechanism which triggered and maintained the precipitating systems more than 30 h over the Pyrenees. The moisture of the feeding flow mainly came from the Atlantic Ocean (7-9 g kg-1) and the role of the Mediterranean as a local moisture source was very limited (2-3 g kg-1) due to the high initial water vapour content of the parcels and the rapid passage over the basin along the Spanish Mediterranean coast (less than 12 h).

Numerical Modelling of Combustion and Emission Formation in a Diesel Engine with Different Fuel Injection Characteristics

Päästöjen vähentäminen on ollut viime vuosina tärkeässä osassa polttomoottoreita kehitettäessä.Monet viralliset tahot asettavat uusia tiukempia päästörajoituksia. Päästörajatovat tyypillisesti olleet tiukimmat autoteollisuuden valmistamille pienille nopeakäyntisille diesel-moottoreille, mutta viime aikoina paineita on kohdistunut myös suurempiin keskinopeisiin ja hidaskäyntisiin diesel-moottoreihin. Päästörajat ovat erilaisia riippuen moottorin tyypistä, käytetystä polttoaineesta ja paikasta missä moottoria käytetään johtuen erilaisista paikallisista laeista ja asetuksista. Eniten huomiota diesel-moottorin päästöissä täytyy kohdistaa typen oksideihin, savun muodostukseen sekä partikkeleihin. Laskennallisen virtausmekaniikan (CFD) avulla on hyvät mahdollisuudet tutkia diesel-moottorin sylinterissä tapahtuvia ilmiöitä palamisen aikana. CFD on hyödyllinen työkalu arvioitaessa moottorin suorituskykyä ja päästöjen muodostumista. CFD:llä on mahdollista testata erilaisten parametrien ja geometrioiden vaikutusta ilman kalliita moottorinkoeajoja. CFD:tä voidaan käyttää myös opetustarkoituksessa lisäämään paloprosessin tuntemusta. Tulevaisuudessa palamissimuloinnit CFD:llä tulevat epäilemättä olemaan tärkeä osa moottorin kehityksessä. Tässä diplomityössä on tehty palamissimuloinnit kahteen erilaisilla poittoaineenruiskutuslaitteistoilla varustettuun Wärtsilän keskinopeaan diesel-moottoriin. W46 moottorin ruiskutuslaitteisto on perinteinen mekaanisesti ohjattu pumppusuutin ja W46-CR moottorissa on elektronisesti ohjattu 'common rail' ruiskutuslaitteisto. Näiden moottorien ja käytössä olevien ruiskutusprofiilien lisäksi on simuloinneilla testattu erilaisia uusia ruiskutusprofiileja, jotta erityyppisten profiilien hyvät ja huonot ominaisuudet tulisivat selville. Matalalla kuormalla kiinnostuksen kohteena on nokipäästöjen muodostus ja täydellä kuormalla NOx-päästöjen muodostus ja polttoaineen kulutus. Simulointien tulokset osoittivat, että noen muodostusta matalalla kuormalla voidaan selvästi vähentää monivaiheisella ruiskutuksella, jossa yksi ruiskutusjakso jaetaan kahteen tai useampaan jaksoon. Erityisen tehokas noen vähentämisessä vaikuttaa olevan ns. jälkiruiskutus (post injection). Matalat NOx-päästöt ja hyvä polttoaineen kulutus täydellä kuormalla on mahdollista saavuttaaasteittain nostettavalla ruiskutusnopeudella.

Freezing Point Depressions of Dilute Solutions of Alkali Metal Chlorides and Bromides

Freezing point depressions (¿Tf) of dilute solutions of several alkali metal chlorides and bromides were calculated by means of the best activity coefficient equations. In the calculations, Hückel, Hamer and Pitzer equationswere used for activity coefficients. The experimental ¿Tf values available in the literature for dilute LiCl, NaCl and KBr solutions can be predicted within experimental error by the Hückel equations used. The experimental ¿Tf values for dilute LiCl and KBr solutions can also be accurately calculated by corresponding Pitzer equations and those for dilute NaCl solutions by the Hamer equation for this salt. Neither Hamer nor Pitzer equations predict accurately the freezing points reported in the literature for LiBr and NaBr solutions. The ¿Tf values available for dilute solutions of RbCl, CsCl or CsBr are not known at the moment accurately because the existing data for these solutions are not precise. The freezing point depressions are tabulated in the present study for LiCl, NaCl and KBr solutions at several rounded molalities. The ¿Tf values in this table can be highly recommended. The activity coefficient equations used in the calculation of these values have been tested with almost allhigh-precision electrochemical data measured at 298.15 K.

Numerical analysis of including phase change material in a domestic hot water cylinder

The domestic hot water cylinder incorporates encapsulated pcm placed in 57 vertical pipes. The use of PCM increases the thermal energy storage capacity of the cylinder and allows the use of low cost electricity during low peak periods. After experimental validation the numerical model developed in the project will be used to optimize the distribution of the pcm inside the water tank.

Experimental and Numerical Analysis of Different Volutes in a Centrifugal Compressor

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

Theoretical and numerical study of thermomagnetic convection in magnetic fluids

In this thesis, the magnetic field control of convection instabilities and heat and mass transfer processesin magnetic fluids have been investigated by numerical simulations and theoretical considerations. Simulation models based on finite element and finite volume methods have been developed. In addition to standard conservation equations, themagnetic field inside the simulation domain is calculated from Maxwell equations and the necessary terms to take into account for the magnetic body force and magnetic dissipation have been added to the equations governing the fluid motion.Numerical simulations of magnetic fluid convection near the threshold supportedexperimental observations qualitatively. Near the onset of convection the competitive action of thermal and concentration density gradients leads to mostly spatiotemporally chaotic convection with oscillatory and travelling wave regimes, previously observed in binary mixtures and nematic liquid crystals. In many applications of magnetic fluids, the heat and mass transfer processes including the effects of external magnetic fields are of great importance. In addition to magnetic fluids, the concepts and the simulation models used in this study may be applied also to the studies of convective instabilities in ordinary fluids as well as in other binary mixtures and complex fluids.

Impacts of upland open drains upon runoff generation: a numerical assessment of catchment-scale impacts

Shallow upland drains, grips, have been hypothesized as responsible for increased downstream flow magnitudes. Observations provide counterfactual evidence, often relating to the difficulty of inferring conclusions from statistical correlation and paired catchment comparisons, and the complexity of designing field experiments to test grip impacts at the catchment scale. Drainage should provide drier antecedent moisture conditions, providing more storage at the start of an event; however, grips have higher flow velocities than overland flow, thus potentially delivering flow more rapidly to the drainage network. We develop and apply a model for assessing the impacts of grips on flow hydrographs. The model was calibrated on the gripped case, and then the gripped case was compared with the intact case by removing all grips. This comparison showed that even given parameter uncertainty, the intact case had significantly higher flood peaks and lower baseflows, mirroring field observations of the hydrological response of intact peat. The simulations suggest that this is because delivery effects may not translate into catchment-scale impacts for three reasons. First, in our case, the proportions of flow path lengths that were hillslope were not changed significantly by gripping. Second, the structure of the grip network as compared with the structure of the drainage basin mitigated against grip-related increases in the concentration of runoff in the drainage network, although it did marginally reduce the mean timing of that concentration at the catchment outlet. Third, the effect of the latter upon downstream flow magnitudes can only be assessed by reference to the peak timing of other tributary basins, emphasizing that drain effects are both relative and scale dependent. However, given the importance of hillslope flow paths, we show that if upland drainage causes significant changes in surface roughness on hillslopes, then critical and important feedbacks may impact upon the speed of hydrological response. Copyright (c) 2012 John Wiley & Sons, Ltd.

PHOTOCATALYTIC OXIDATION OF AQUEOUS SOLUTIONS OF JET FUEL AND ICING INHIBITORS

The present paper is devoted to the results of experimental research undertaken into photocatalytical oxidation (PCO) of aqueous solutions of de-icing agents and aqueous extract of jet fuel. The report consists of introduction, literature review, description of materials and methods, discussion of results and conclusions. TiO2 was selected as a photocatalyst for the experiments with synthetic solutions of ethylene glycol, 2-ethoxyethanol and aqueous extract of jet fuel. To explain the PCO mechanisms affecting certain behaviour of de-icing agent under distinctive conditions, the following factors were studied: the impact of initial concentration of pollutant, the role of pH, the presence of tert-butanol as OH·-radicals scavenger and mineral admixtures. PCO under solar radiation performed in two ways: catalysed by irradiated TiO2 slurry or by TiO2 attached to buoyant hollow glass micro-spheres. Special attention was paid to the energy-saving PCO with reduced intensity mixing of the slurry. The effect of PCO was assessed by determination of residual chemical oxygen demand of solution (COD) and by measuring of concentration of glycols. The PCO process efficiency was assumed to be dependent on the TiO2 suspension fractional composition. Thus, the following effects of solutions’ media were viewed: presence of organic admixtures, pH influence, mixing mode during the PCO. The effects of mineral admixtures - Ca2+, Fe3+/2+, Mn2+, SO42- - that are often present in natural and wastewater systems or produced during the degradation of organic pollutants and which can affect the rate of PCO of de-icing agents, were also investigated.

Three-dimensional hydrodynamic numerical modelling of fold nappe formation, basement-cover deformation and slab detachment with applications to the helvetic nappe system (W Switzerland)

Many three-dimensional (3-D) structures in rock, which formed during the deformation of the Earth's crust and lithosphere, are controlled by a difference in mechanical strength between rock units and are often the result of a geometrical instability. Such structures are, for example, folds, pinch-and-swell structures (due to necking) or cuspate-lobate structures (mullions). These struc-tures occur from the centimeter to the kilometer scale and the related deformation processes con-trol the formation of, for example, fold-and-thrust belts and extensional sedimentary basins or the deformation of the basement-cover interface. The 2-D deformation processes causing these structures are relatively well studied, however, several processes during large-strain 3-D defor-mation are still incompletely understood. One of these 3-D processes is the lateral propagation of these structures, such as fold and cusp propagation in a direction orthogonal to the shortening direction or neck propagation in direction orthogonal to the extension direction. Especially, we are interested in fold nappes which are recumbent folds with amplitudes usually exceeding 10 km and they have been presumably formed by ductile shearing. They often exhibit a constant sense of shearing and a non-linear increase of shear strain towards their overturned limb. The fold axes of the Morcles fold nappe in western Switzerland plunges to the ENE whereas the fold axes in the more eastern Doldenhorn nappe plunges to the WSW. These opposite plunge direc-tions characterize the Rawil depression (Wildstrubel depression). The Morcles nappe is mainly the result of layer parallel contraction and shearing. During the compression the massive lime-stones were more competent than the surrounding marls and shales, which led to the buckling characteristics of the Morcles nappe, especially in the north-dipping normal limb. The Dolden-horn nappe exhibits only a minor overturned fold limb. There are still no 3-D numerical studies which investigate the fundamental dynamics of the formation of the large-scale 3-D structure including the Morcles and Doldenhorn nappes and the related Rawil depression. We study the 3-D evolution of geometrical instabilities and fold nappe formation with numerical simulations based on the finite element method (FEM). Simulating geometrical instabilities caused by sharp variations of mechanical strength between rock units requires a numerical algorithm that can accurately resolve material interfaces for large differences in material properties (e.g. between limestone and shale) and for large deformations. Therefore, our FE algorithm combines a nu-merical contour-line technique and a deformable Lagrangian mesh with re-meshing. With this combined method it is possible to accurately follow the initial material contours with the FE mesh and to accurately resolve the geometrical instabilities. The algorithm can simulate 3-D de-formation for a visco-elastic rheology. The viscous rheology is described by a power-law flow law. The code is used to study the 3-D fold nappe formation, the lateral propagation of folding and also the lateral propagation of cusps due to initial half graben geometry. Thereby, the small initial geometrical perturbations for folding and necking are exactly followed by the FE mesh, whereas the initial large perturbation describing a half graben is defined by a contour line inter-secting the finite elements. Further, the 3-D algorithm is applied to 3-D viscous nacking during slab detachment. The results from various simulations are compared with 2-D resulats and a 1-D analytical solution. -- On retrouve beaucoup de structures en 3 dimensions (3-D) dans les roches qui ont pour origines une déformation de la lithosphère terrestre. Ces structures sont par exemple des plis, des boudins (pinch-and-swell) ou des mullions (cuspate-lobate) et sont présentés de l'échelle centimétrique à kilométrique. Mécaniquement, ces structures peuvent être expliquées par une différence de résistance entre les différentes unités de roches et sont généralement le fruit d'une instabilité géométrique. Ces différences mécaniques entre les unités contrôlent non seulement les types de structures rencontrées, mais également le type de déformation (thick skin, thin skin) et le style tectonique (bassin d'avant pays, chaîne d'avant pays). Les processus de la déformation en deux dimensions (2-D) formant ces structures sont relativement bien compris. Cependant, lorsque l'on ajoute la troisiéme dimension, plusieurs processus ne sont pas complètement compris lors de la déformation à large échelle. L'un de ces processus est la propagation latérale des structures, par exemple la propagation de plis ou de mullions dans la direction perpendiculaire à l'axe de com-pression, ou la propagation des zones d'amincissement des boudins perpendiculairement à la direction d'extension. Nous sommes particulièrement intéressés les nappes de plis qui sont des nappes de charriage en forme de plis couché d'une amplitude plurikilométrique et étant formées par cisaillement ductile. La plupart du temps, elles exposent un sens de cisaillement constant et une augmentation non linéaire de la déformation vers la base du flanc inverse. Un exemple connu de nappes de plis est le domaine Helvétique dans les Alpes de l'ouest. Une de ces nap-pes est la Nappe de Morcles dont l'axe de pli plonge E-NE tandis que de l'autre côté de la dépression du Rawil (ou dépression du Wildstrubel), la nappe du Doldenhorn (équivalent de la nappe de Morcles) possède un axe de pli plongeant O-SO. La forme particulière de ces nappes est due à l'alternance de couches calcaires mécaniquement résistantes et de couches mécanique-ment faibles constituées de schistes et de marnes. Ces différences mécaniques dans les couches permettent d'expliquer les plissements internes à la nappe, particulièrement dans le flanc inver-se de la nappe de Morcles. Il faut également noter que le développement du flanc inverse des nappes n'est pas le même des deux côtés de la dépression de Rawil. Ainsi la nappe de Morcles possède un important flanc inverse alors que la nappe du Doldenhorn en est presque dépour-vue. A l'heure actuelle, aucune étude numérique en 3-D n'a été menée afin de comprendre la dynamique fondamentale de la formation des nappes de Morcles et du Doldenhorn ainsi que la formation de la dépression de Rawil. Ce travail propose la première analyse de l'évolution 3-D des instabilités géométriques et de la formation des nappes de plis en utilisant des simulations numériques. Notre modèle est basé sur la méthode des éléments finis (FEM) qui permet de ré-soudre avec précision les interfaces entre deux matériaux ayant des propriétés mécaniques très différentes (par exemple entre les couches calcaires et les couches marneuses). De plus nous utilisons un maillage lagrangien déformable avec une fonction de re-meshing (production d'un nouveau maillage). Grâce à cette méthode combinée il nous est possible de suivre avec précisi-on les interfaces matérielles et de résoudre avec précision les instabilités géométriques lors de la déformation de matériaux visco-élastiques décrit par une rhéologie non linéaire (n>1). Nous uti-lisons cet algorithme afin de comprendre la formation des nappes de plis, la propagation latérale du plissement ainsi que la propagation latérale des structures de type mullions causé par une va-riation latérale de la géométrie (p.ex graben). De plus l'algorithme est utilisé pour comprendre la dynamique 3-D de l'amincissement visqueux et de la rupture de la plaque descendante en zone de subduction. Les résultats obtenus sont comparés à des modèles 2-D et à la solution analytique 1-D. -- Viele drei dimensionale (3-D) Strukturen, die in Gesteinen vorkommen und durch die Verfor-mung der Erdkruste und Litosphäre entstanden sind werden von den unterschiedlichen mechani-schen Eigenschaften der Gesteinseinheiten kontrolliert und sind häufig das Resulat von geome-trischen Istabilitäten. Zu diesen strukturen zählen zum Beispiel Falten, Pich-and-swell Struktu-ren oder sogenannte Cusbate-Lobate Strukturen (auch Mullions). Diese Strukturen kommen in verschiedenen Grössenordungen vor und können Masse von einigen Zentimeter bis zu einigen Kilometer aufweisen. Die mit der Entstehung dieser Strukturen verbundenen Prozesse kontrol-lieren die Entstehung von Gerbirgen und Sediment-Becken sowie die Verformung des Kontaktes zwischen Grundgebirge und Stedimenten. Die zwei dimensionalen (2-D) Verformungs-Prozesse die zu den genannten Strukturen führen sind bereits sehr gut untersucht. Einige Prozesse wäh-rend starker 3-D Verformung sind hingegen noch unvollständig verstanden. Einer dieser 3-D Prozesse ist die seitliche Fortpflanzung der beschriebenen Strukturen, so wie die seitliche Fort-pflanzung von Falten und Cusbate-Lobate Strukturen senkrecht zur Verkürzungsrichtung und die seitliche Fortpflanzung von Pinch-and-Swell Strukturen othogonal zur Streckungsrichtung. Insbesondere interessieren wir uns für Faltendecken, liegende Falten mit Amplituden von mehr als 10 km. Faltendecken entstehen vermutlich durch duktile Verscherung. Sie zeigen oft einen konstanten Scherungssinn und eine nicht-lineare zunahme der Scherverformung am überkipp-ten Schenkel. Die Faltenachsen der Morcles Decke in der Westschweiz fallen Richtung ONO während die Faltenachsen der östicher gelegenen Doldenhorn Decke gegen WSW einfallen. Diese entgegengesetzten Einfallrichtungen charakterisieren die Rawil Depression (Wildstrubel Depression). Die Morcles Decke ist überwiegend das Resultat von Verkürzung und Scherung parallel zu den Sedimentlagen. Während der Verkürzung verhielt sich der massive Kalkstein kompetenter als der Umliegende Mergel und Schiefer, was zur Verfaltetung Morcles Decke führ-te, vorallem in gegen Norden eifallenden überkippten Schenkel. Die Doldenhorn Decke weist dagegen einen viel kleineren überkippten Schenkel und eine stärkere Lokalisierung der Verfor-mung auf. Bis heute gibt es keine 3-D numerischen Studien, die die fundamentale Dynamik der Entstehung von grossen stark verformten 3-D Strukturen wie den Morcles und Doldenhorn Decken sowie der damit verbudenen Rawil Depression untersuchen. Wir betrachten die 3-D Ent-wicklung von geometrischen Instabilitäten sowie die Entstehung fon Faltendecken mit Hilfe von numerischen Simulationen basiert auf der Finite Elemente Methode (FEM). Die Simulation von geometrischen Instabilitäten, die aufgrund von Änderungen der Materialeigenschaften zwischen verschiedenen Gesteinseinheiten entstehen, erfortert einen numerischen Algorithmus, der in der Lage ist die Materialgrenzen mit starkem Kontrast der Materialeigenschaften (zum Beispiel zwi-schen Kalksteineinheiten und Mergel) für starke Verfomung genau aufzulösen. Um dem gerecht zu werden kombiniert unser FE Algorithmus eine numerische Contour-Linien-Technik und ein deformierbares Lagranges Netz mit Re-meshing. Mit dieser kombinierten Methode ist es mög-lich den anfänglichen Materialgrenzen mit dem FE Netz genau zu folgen und die geometrischen Instabilitäten genügend aufzulösen. Der Algorithmus ist in der Lage visko-elastische 3-D Ver-formung zu rechnen, wobei die viskose Rheologie mit Hilfe eines power-law Fliessgesetzes beschrieben wird. Mit dem numerischen Algorithmus untersuchen wir die Entstehung von 3-D Faltendecken, die seitliche Fortpflanzung der Faltung sowie der Cusbate-Lobate Strukturen die sich durch die Verkürzung eines mit Sediment gefüllten Halbgraben bilden. Dabei werden die anfänglichen geometrischen Instabilitäten der Faltung exakt mit dem FE Netz aufgelöst wäh-rend die Materialgranzen des Halbgrabens die Finiten Elemente durchschneidet. Desweiteren wird der 3-D Algorithmus auf die Einschnürung während der 3-D viskosen Plattenablösung und Subduktion angewandt. Die 3-D Resultate werden mit 2-D Ergebnissen und einer 1-D analyti-schen Lösung verglichen.

Logistic Solutions of Imported Renovation Materials

Tämän diplomityön tavoitteena oli tutkia saneeraus- ja remontoimistuotteiden osalta maahantuontia Suomeen ja transitoliikennettä Suomen kautta Venäjälle. Lähtökohtana oli Suomen logistinen asema ja mahdolliset kilpailukykyyn vaikuttavat tekijät, joista on luotu potentiaalisia tulevaisuuden kuvia. Kaikki tulokset johdettiin LOADER/CLIENT –tutkimusprojektin yhteyksistä. Tutkimuksen tavoitteisiin päästiin empiiristä tutkimusta varten kootun teorian kautta, joka koostuu toimitusketjun hallinnasta ja tämänhetkisestä parhaimmasta käytännöstä eli ns. kolmen A:n menetelmästä. Lisäksi teoriaosuudessa tarkasteltiin toimitusketjun yritysten välisiä suhteita, kansainväliseen toimitusketjuun liittyviä tekijöitä sekä tulevaisuuden haasteita. Työn empiiristä tutkimustietoa kerättiin pääasiassa haastattelujen, mutta myös muiden tutkimuksen tuloksien kautta. Tutkimusmenetelmä oli laadullinen puolistrukturoitu teemahaastattelu, jota varten ryhmä remontoimiseen ja saneeraukseen liittyviä yrityksiä ja heidän käyttämiä logistisia operaattoreita sekä projektin johtoryhmää haastateltiin. Tutkimuksen tulosten mukaan suurin osa haastatelluista yrityksistä näkee Suomen logistisen aseman hyvänä tällä hetkellä ja lähivuosina. Etenkin arvotavaroiden transito- kuljetuksissa Suomen asema nähdään vahvana vuosien kokemuksen ansiosta. Tällä hetkellä arvotavaroiden toimitusaika ja kokonaiskustannukset varastoimisineen ja käsittelyineen ovat vastaavia tai jopa edullisempia verrattuna esimerkiksi reitteihin Baltian maiden kautta. Baltian maiden kehitys voi heikentää Suomen kilpailukykyä hetkellisesti, mutta todennäköisesti ei pidemmällä aikavälillä. Haastateltujen henkilöiden mielestä Suomen logistiseen kilpailukykyyn tulevaisuudessa vaikuttaa eniten kaksi kriittistä tekijää. Suomen tulisi ensinnäkin keskittyä ja kehittää arvotavaran transitokuljetuksia. Toiseksi, säilyttääkseen kilpailukykynsä, Suomen satamien palveluiden kuten tullin ja ahtaajien tulisi olla joustavasti käytettävissä 24 tuntia vuorokaudessa vuoden jokaisena päivänä.

The role of intramolecular barriers on the glass transition of polymers: computer simulations versus mode coupling theory

We present computer simulations of a simple bead-spring model for polymer melts with intramolecular barriers. By systematically tuning the strength of the barriers, we investigate their role on the glass transition. Dynamic observables are analyzed within the framework of the mode coupling theory (MCT). Critical nonergodicity parameters, critical temperatures, and dynamic exponents are obtained from consistent fits of simulation data to MCT asymptotic laws. The so-obtained MCT λ-exponent increases from standard values for fully flexible chains to values close to the upper limit for stiff chains. In analogy with systems exhibiting higher-order MCT transitions, we suggest that the observed large λ-values arise form the interplay between two distinct mechanisms for dynamic arrest: general packing effects and polymer-specific intramolecular barriers. We compare simulation results with numerical solutions of the MCT equations for polymer systems, within the polymer reference interaction site model (PRISM) for static correlations. We verify that the approximations introduced by the PRISM are fulfilled by simulations, with the same quality for all the range of investigated barrier strength. The numerical solutions reproduce the qualitative trends of simulations for the dependence of the nonergodicity parameters and critical temperatures on the barrier strength. In particular, the increase in the barrier strength at fixed density increases the localization length and the critical temperature. However the qualitative agreement between theory and simulation breaks in the limit of stiff chains. We discuss the possible origin of this feature.

Mesoscale numerical analysis of the historical November 1982 heavy precipitation event over Andorra (Eastern Pyrenees)

From 6 to 8 November 1982 one of the most catastrophic flash-flood events was recorded in the Eastern Pyrenees affecting Andorra and also France and Spain with rainfall accumulations exceeding 400 mm in 24 h, 44 fatalities and widespread damage. This paper aims to exhaustively document this heavy precipitation event and examines mesoscale simulations performed by the French Meso-NH non-hydrostatic atmospheric model. Large-scale simulations show the slow-evolving synoptic environment favourable for the development of a deep Atlantic cyclone which induced a strong southerly flow over the Eastern Pyrenees. From the evolution of the synoptic pattern four distinct phases have been identified during the event. The mesoscale analysis presents the second and the third phase as the most intense in terms of rainfall accumulations and highlights the interaction of the moist and conditionally unstable flows with the mountains. The presence of a SW low level jet (30 m s-1) around 1500 m also had a crucial role on focusing the precipitation over the exposed south slopes of the Eastern Pyrenees. Backward trajectories based on Eulerian on-line passive tracers indicate that the orographic uplift was the main forcing mechanism which triggered and maintained the precipitating systems more than 30 h over the Pyrenees. The moisture of the feeding flow mainly came from the Atlantic Ocean (7-9 g kg-1) and the role of the Mediterranean as a local moisture source was very limited (2-3 g kg-1) due to the high initial water vapour content of the parcels and the rapid passage over the basin along the Spanish Mediterranean coast (less than 12 h).