A simple solution of the laminar dam break wave

The sudden release of a mass of fluid in a channel generates a highly unsteady flow motion, called dam break wave. While industrial fluids exhibit sometimes non-Newtonian behaviours, the viscous fluid flow assumption remains a useful approximation for simplified analyses. In this study, new solutions of laminar dam break wave are proposed for a semi-infinite reservoir based upon the method of characteristics. The solutions yield simple explicit expressions of the wave front location, wave front celerity and instantaneous free-surface profiles that compare favourably with experimental observations. Both horizontal and sloping channel configurations are treated. The simplicity of the equations may allow future extension to more complicated fluid flows.

Coletores de distribuição e sua otimização hidráulica, aplicados a uma instalação de produção de água arrefecida em sistemas de AVAC ( Fórum Picoas)

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia Mecânica

Analytical and numerical study of the electro-osmotic annular flow of viscoelastic fluids

In this work we present semi-analytical solutions for the electro-osmotic annular flow of viscoelastic fluids modeled by the Linear and Exponential PTT models. The viscoelastic fluid flows in the axial direction between two concentric cylinders under the combined influences of electrokinetic and pressure forcings. The analysis invokes the Debye-Hückel approximation and includes the limit case of pure electro-osmotic flow. The solution is valid for both no slip and slip velocity at the walls and the chosen slip boundary condition is the linear Navier slip velocity model. The combined effects of fluid rheology, electro-osmotic and pressure gradient forcings on the fluid velocity distribution are also discussed.

Measurements of the bulk and interfacial velocity profiles in oscillating Newtonian and Maxwellian fluids

We present the dynamic velocity profiles of a Newtonian fluid (glycerol) and a viscoelastic Maxwell fluid (CPyCl-NaSal in water) driven by an oscillating pressure gradient in a vertical cylindrical pipe. The frequency range explored has been chosen to include the first three resonance peaks of the dynamic permeability of the viscoelastic-fluid¿pipe system. Three different optical measurement techniques have been employed. Laser Doppler anemometry has been used to measure the magnitude of the velocity at the center of the liquid column. Particle image velocimetry and optical deflectometry are used to determine the velocity profiles at the bulk of the liquid column and at the liquid-air interface respectively. The velocity measurements in the bulk are in good agreement with the theoretical predictions of a linear theory. The results, however, show dramatic differences in the dynamic behavior of Newtonian and viscoelastic fluids, and demonstrate the importance of resonance phenomena in viscoelastic fluid flows, biofluids in particular, in confined geometries.

Modélisation par automate cellulaire des phénomènes diagénétiques des plateformes carbonatées : calibration et paramétrage à partir de deux cas d'études : l'Urgonien du Vercors (Crétacé inférieur, SE France) et les Calcaires Gris du Mont Compomolon (Lias, NE Italie)

Une fois déposé, un sédiment est affecté au cours de son enfouissement par un ensemble de processus, regroupé sous le terme diagenèse, le transformant parfois légèrement ou bien suffisamment pour le rendre méconnaissable. Ces modifications ont des conséquences sur les propriétés pétrophysiques qui peuvent être positives ou négatives, c'est-à-dire les améliorer ou bien les détériorer. Une voie alternative de représentation numérique des processus, affranchie de l'utilisation des réactions physico-chimiques, a été adoptée et développée en mimant le déplacement du ou des fluides diagénétiques. Cette méthode s'appuie sur le principe d'un automate cellulaire et permet de simplifier les phénomènes sans sacrifier le résultat et permet de représenter les phénomènes diagénétiques à une échelle fine. Les paramètres sont essentiellement numériques ou mathématiques et nécessitent d'être mieux compris et renseignés à partir de données réelles issues d'études d'affleurements et du travail analytique effectué. La représentation des phénomènes de dolomitisation de faible profondeur suivie d'une phase de dédolomitisation a été dans un premier temps effectuée. Le secteur concerne une portion de la série carbonatée de l'Urgonien (Barrémien-Aptien), localisée dans le massif du Vercors en France. Ce travail a été réalisé à l'échelle de la section afin de reproduire les géométries complexes associées aux phénomènes diagénétiques et de respecter les proportions mesurées en dolomite. De plus, la dolomitisation a été simulée selon trois modèles d'écoulement. En effet, la dédolomitisation étant omniprésente, plusieurs hypothèses sur le mécanisme de dolomitisation ont été énoncées et testées. Plusieurs phases de dolomitisation per ascensum ont été également simulées sur des séries du Lias appartenant aux formations du groupe des Calcaire Gris, localisées au nord-est de l'Italie. Ces fluides diagénétiques empruntent le réseau de fracturation comme vecteur et affectent préférentiellement les lithologies les plus micritisées. Cette étude a permis de mettre en évidence la propagation des phénomènes à l'échelle de l'affleurement. - Once deposited, sediment is affected by diagenetic processes during their burial history. These diagenetic processes are able to affect the petrophysical properties of the sedimentary rocks and also improve as such their reservoir capacity. The modelling of diagenetic processes in carbonate reservoirs is still a challenge as far as neither stochastic nor physicochemical simulations can correctly reproduce the complexity of features and the reservoir heterogeneity generated by these processes. An alternative way to reach this objective deals with process-like methods, which simplify the algorithms while preserving all geological concepts in the modelling process. The aim of the methodology is to conceive a consistent and realistic 3D model of diagenetic overprints on initial facies resulting in petrophysical properties at a reservoir scale. The principle of the method used here is related to a lattice gas automata used to mimic diagenetic fluid flows and to reproduce the diagenetic effects through the evolution of mineralogical composition and petrophysical properties. This method developed in a research group is well adapted to handle dolomite reservoirs through the propagation of dolomitising fluids and has been applied on two case studies. The first study concerns a mid-Cretaceous rudist and granular platform of carbonate succession (Urgonian Fm., Les Gorges du Nan, Vercors, SE France), in which several main diagenetic stages have been identified. The modelling in 2D is focused on dolomitisation followed by a dédolomitisation stage. For the second study, data collected from outcrops on the Venetian platform (Lias, Mont Compomolon NE Italy), in which several diagenetic stages have been identified. The main one is related to per ascensum dolomitisation along fractures. In both examples, the evolution of the effects of the mimetic diagenetic fluid on mineralogical composition can be followed through space and numerical time and help to understand the heterogeneity in reservoir properties. Carbonates, dolomitisation, dédolomitisation, process-like modelling, lattice gas automata, random walk, memory effect.

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.

Process intensification: From optimised flow patterns to microprocess technology

This thesis is focused on process intensification. Several significant problems and applications of this theme are covered. Process intensification is nowadays one of the most popular trends in chemical engineering and attempts have been made to develop a general, systematic methodology for intensification. This seems, however, to be very difficult, because intensified processes are often based on creativity and novel ideas. Monolith reactors and microreactors are successful examples of process intensification. They are usually multichannel devices in which a proper feed technique is important for creating even fluid distribution into the channels. Two different feed techniques were tested for monoliths. In the first technique a shower method was implemented by means of perforated plates. The second technique was a dispersion method using static mixers. Both techniques offered stable operation and uniform fluid distribution. The dispersion method enabled a wider operational range in terms of liquid superficial velocity. Using dispersion method, a volumetric gas-liquid mass transfer coefficient of 2 s-1 was reached. Flow patterns play a significant role in terms of the mixing performance of micromixers. Although the geometry of a T-mixer is simple, channel configurations and dimensions had a clear effect on mixing efficiency. The flow in the microchannel was laminar, but the formation of vortices promoted mixing in micro T-mixers. The generation of vortices was dependent on the channel dimensions, configurations and flow rate. Microreactors offer a high ratio of surface area to volume. Surface forces and interactions between fluids and surfaces are, therefore, often dominant factors. In certain cases, the interactions can be effectively utilised. Different wetting properties of solid materials (PTFE and stainless steel) were applied in the separation of immiscible liquid phases. A micro-scale plate coalescer with hydrophilic and hydrophobic surfaces was used for the continuous separation of organic and aqueous phases. Complete phase separation occurred in less than 20 seconds, whereas the separation time by settling exceeded 30 min. Fluid flows can be also intensified in suitable conditions. By adding certain additives into turbulent fluid flow, it was possible to reduce friction (drag) by 40 %. Drag reduction decreases frictional pressure drop in pipelines which leads to remarkable energy savings and decreases the size or number of pumping facilities required, e.g., in oil transport pipes. Process intensification enables operation often under more optimal conditions. The consequent cost savings from reduced use of raw materials and reduced waste lead to greater economic benefits in processing.

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.

Topological Invariants in Hydrodynamics and Hydromagnetics

In this thesis we are studying possible invariants in hydrodynamics and hydromagnetics. The concept of flux preservation and line preservation of vector fields, especially vorticity vector fields, have been studied from the very beginning of the study of fluid mechanics by Helmholtz and others. In ideal magnetohydrodynamic flows the magnetic fields satisfy the same conservation laws as that of vorticity field in ideal hydrodynamic flows. Apart from these there are many other fields also in ideal hydrodynamic and magnetohydrodynamic flows which preserves flux across a surface or whose vector lines are preserved. A general study using this analogy had not been made for a long time. Moreover there are other physical quantities which are also invariant under the flow, such as Ertel invariant. Using the calculus of differential forms Tur and Yanovsky classified the possible invariants in hydrodynamics. This mathematical abstraction of physical quantities to topological objects is needed for an elegant and complete analysis of invariants.Many authors used a four dimensional space-time manifold for analysing fluid flows. We have also used such a space-time manifold in obtaining invariants in the usual three dimensional flows.In chapter one we have discussed the invariants related to vorticity field using vorticity field two form w2 in E4. Corresponding to the invariance of four form w2 ^ w2 we have got the invariance of the quantity E. w. We have shown that in an isentropic flow this quantity is an invariant over an arbitrary volume.In chapter three we have extended this method to any divergence-free frozen-in field. In a four dimensional space-time manifold we have defined a closed differential two form and its potential one from corresponding to such a frozen-in field. Using this potential one form w1 , it is possible to define the forms dw1 , w1 ^ dw1 and dw1 ^ dw1 . Corresponding to the invariance of the four form we have got an additional invariant in the usual hydrodynamic flows, which can not be obtained by considering three dimensional space.In chapter four we have classified the possible integral invariants associated with the physical quantities which can be expressed using one form or two form in a three dimensional flow. After deriving some general results which hold for an arbitrary dimensional manifold we have illustrated them in the context of flows in three dimensional Euclidean space JR3. If the Lie derivative of a differential p-form w is not vanishing,then the surface integral of w over all p-surfaces need not be constant of flow. Even then there exist some special p-surfaces over which the integral is a constant of motion, if the Lie derivative of w satisfies certain conditions. Such surfaces can be utilised for investigating the qualitative properties of a flow in the absence of invariance over all p-surfaces. We have also discussed the conditions for line preservation and surface preservation of vector fields. We see that the surface preservation need not imply the line preservation. We have given some examples which illustrate the above results. The study given in this thesis is a continuation of that started by Vedan et.el. As mentioned earlier, they have used a four dimensional space-time manifold to obtain invariants of flow from variational formulation and application of Noether's theorem. This was from the point of view of hydrodynamic stability studies using Arnold's method. The use of a four dimensional manifold has great significance in the study of knots and links. In the context of hydrodynamics, helicity is a measure of knottedness of vortex lines. We are interested in the use of differential forms in E4 in the study of vortex knots and links. The knowledge of surface invariants given in chapter 4 may also be utilised for the analysis of vortex and magnetic reconnections.

Ontogenetic Changes in the Bell Morphology and Kinematics and Swimming Behavior of Rowing Medusae: the Special Case of the Limnomedusa Liriope tetraphylla

Swimming animals may experience significant changes in the Reynolds number (Re) of their surrounding fluid flows throughout ontogeny. Many medusae experience Re environments with significant viscous forces as small juveniles but inertially dominated Re environments as adults. These different environments may affect their propulsive strategies. In particular, rowing, a propulsive strategy with ecological advantages for large adults, may be constrained by viscosity for small juvenile medusae. We examined changes in the bell morphology and swimming kinematics of the limnomedusa Liriope tetraphylla at different stages of development. L. tetraphylla maintained an oblate bell (fineness ratio approximate to 0.5-0.6), large velar aperture ratio (R(v) approximate to 0.5-0.8), and rapid bell kinematics throughout development. These traits enabled it to use rowing propulsion at all stages except the very smallest sizes observed (diameter = 0.14 cm). During the juvenile stage, very rapid bell kinematics served to increase Re sufficiently for rowing propulsion. Other taxa that use rowing propulsion as adults, such as leptomedusae and scyphomedusae, typically utilize different propulsive strategies as small juveniles to function in low Re environments. We compared the performance values of the different propulsive modes observed among juvenile medusae.

Demonstration of Solar Heating and Cooling System using Sorption Integrated Solar Thermal Collectors

Producing cost-competitive small and medium-sized solar cooling systems is currently a significant challenge. Due to system complexity, extensive engineering, design and equipment costs; the installation costs of solar thermal cooling systems are prohibitively high. In efforts to overcome these limitations, a novel sorption heat pump module has been developed and directly integrated into a solar thermal collector. The module comprises a fully encapsulated sorption tube containing hygroscopic salt sorbent and water as a refrigerant, sealed under vacuum with no moving parts. A 5.6m2 aperture area outdoor laboratory-scale system of sorption module integrated solar collectors was installed in Stockholm, Sweden and evaluated under constant re-cooling and chilled fluid return temperatures in order to assess collector performance. Measured average solar cooling COP was 0.19 with average cooling powers between 120 and 200 Wm-2 collector aperture area. It was observed that average collector cooling power is constant at daily insolation levels above 3.6 kWhm-2 with the cooling energy produced being proportional to solar insolation. For full evaluation of an integrated sorption collector solar heating and cooling system, under the umbrella of a European Union project for technological innovation, a 180 m2 large-scale demonstration system has been installed in Karlstad, Sweden. Results from the installation commissioned in summer 2014 with non-optimised control strategies showed average electrical COP of 10.6 and average cooling powers between 140 and 250 Wm-2 collector aperture area. Optimisation of control strategies, heat transfer fluid flows through the collectors and electrical COP will be carried out in autumn 2014.

Simulação do escoamento monofásico em um estágio de uma bomba centrífuga utilizando técnicas de fluidodinâmica computacional

Oil production and exploration techniques have evolved in the last decades in order to increase fluid flows and optimize how the required equipment are used. The base functioning of Electric Submersible Pumping (ESP) lift method is the use of an electric downhole motor to move a centrifugal pump and transport the fluids to the surface. The Electric Submersible Pumping is an option that has been gaining ground among the methods of Artificial Lift due to the ability to handle a large flow of liquid in onshore and offshore environments. The performance of a well equipped with ESP systems is intrinsically related to the centrifugal pump operation. It is the pump that has the function to turn the motor power into Head. In this present work, a computer model to analyze the three-dimensional flow in a centrifugal pump used in Electric Submersible Pumping has been developed. Through the commercial program, ANSYS® CFX®, initially using water as fluid flow, the geometry and simulation parameters have been defined in order to obtain an approximation of what occurs inside the channels of the impeller and diffuser pump in terms of flow. Three different geometry conditions were initially tested to determine which is most suitable to solving the problem. After choosing the most appropriate geometry, three mesh conditions were analyzed and the obtained values were compared to the experimental characteristic curve of Head provided by the manufacturer. The results have approached the experimental curve, the simulation time and the model convergence were satisfactory if it is considered that the studied problem involves numerical analysis. After the tests with water, oil was used in the simulations. The results were compared to a methodology used in the petroleum industry to correct viscosity. In general, for models with water and oil, the results with single-phase fluids were coherent with the experimental curves and, through three-dimensional computer models, they are a preliminary evaluation for the analysis of the two-phase flow inside the channels of centrifugal pump used in ESP systems

Experimental analysis of heat recoveries with heat pipes operating in low inclination angle

Heat recovery devices are important in the optimization of thermal systems, since they can be used to reduce thermal losses to the environment. The use of heat pipes in these types of equipment can provide heat recoveries of higher efficiency, since both fluid flows are external and there are less contamination risks between the hot and cold fluids. The objective of this work is to study a heat recovery unit constructed with heat pipes and mainly, to analyze the influence of the inclination of the heat pipes on the performance of the equipment. For this analysis, a heat recovery unit was constructed which possesses 48 finned heat pipes in triangular geometry, the evaporator and condenser being of the same length. This unit was tested in an air-air system simulating a heat recovery process in which heat was supplied to the hot fluid by electrical resistances. The results have shown that there exists an inclination at which the system has a better performance, but for higher inclinations there is no significant increase of the efficiency of the system. This paper also presents the influence of inclination of heat pipes on effectiveness and NTU parameters which are important in heat exchanger design.

A C2-continuous high-resolution upwind convection scheme

A bounded upwinding scheme for numerical solution of hyperbolic conservation laws and Navier-Stokes equations is presented. The scheme is based on convection boundedness criterion and total variation diminishing stability criteria and developed by employing continuously differentiable functions. The accuracy of the scheme is verified by assessing the error and observed convergence rate on 1-D benchmark test cases. A comparative study between the new scheme and conventional total variation diminishing/convection boundedness criterion-based upwind schemes to solve standard nonlinear hyperbolic conservation laws is also accomplished. The scheme is then examined in the simulation of Newtonian and non-Newtonian fluid flows of increasing complexity; a satisfactory agreement has been observed in terms of the overall behavior. Finally, the scheme is used to study the hydrodynamics of a gas-solid flow in a bubbling fluidized bed. © 2013 John Wiley & Sons, Ltd.

Simulação numérica de escoamentos de fluidos com efeitos de rotação

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)