No description available.

English abstract Hydrothermal systems related to highly evolved granitic magmas host diverse mineralization styles and provide an important source of economic metals. This master thesis concentrates on description and interpretation of geological structure, petrographic and textural variability, alteration zoning and calculation of time-integrated fluid fluxes recorded in highly evolved granites and tin-mineralized greisens of the Horní Blatná massif in the Western Krušné hory pluton. The massif is a composite intrusion, which consists of a large number of intrusive units emplaced during two stages. The first stage is represented by sparsely porphyritic fine-grained low-lithian annite granites that can be correlated with marginal granites (G2) of the Fichtelgebirge (Smrčiny) batholith or with intermediate granites (Walfischkopf type) of the Western Krušné hory (Erzgebirge) pluton. Intrusive batches of the second stage progressively evolve from medium- to coarse-grained serial high-lithian annite and zinnwaldite granites with topaz and rare tourmaline towards aphyric fine-grained zinnwaldite (or trilithionite) granites. This suite corresponds to the EIB2 and EIB3 facies of the younger intrusive complex in the Western Krušné hory (Erzgebirge) pluton and it can be compared to the G3 Waldstein and G4 units in...

O método ISPH (do inglês, Incompressible Smoothed Particle Hydrodynamics) é um método de aproximação livre de malha que, através de um conjunto finito de partículas e uma formulação completamente Lagrangeana, permite a solução de diversos tipos de escoamentos. Entretanto, sua aplicação para escoamentos bifásicos ainda é um desafio, principalmente no que refere-se à manutenabilidade da interface entre fluidos. Diante disso, nesta tese é apresentado o desenvolvimento de um código numérico baseado no método ISPH, sendo propostas duas técnicas de tratamento de interface. Para tanto é realizado um estudo a cerca do método, considerando diferentes metodologias, e analisando pontos específicos, tais como a solução do campo de pressões. São apresentados resultados que mostram a eficácia do método, tanto em escoamentos monofásicos, quanto em escoamentos multifásicos, onde, neste caso, são destacadas as melhorias obtidas através das técnicas de tratamento de interface propostas. Por fim, é realizado um estudo do comportamento de misturas bifásicas, com referência ao fenômeno da inversão de fase / Incompressible Smoothed Particle Hydrodynamics (ISPH) method is a meshless approximation that has been used to simulate several types of fluid flows, through a finite particle set and fully lagrangian formulation. The application of ISPH method in two-fluid flow simulations however, has presented many challenges, specially related to the presence of the interface between different fluids. Thus, we present in this study the development of a numerical code based on ISPH, introducing novel interface treatment techniques. A thorough study about this method is provided, considering different methodologies and analysing specific points such as the position of the interface and the obtained pressure field. Results have been presented to show the methods developed in this thesis efficiently simulate two-fluid flows, illustrating the improvements achieved by the proposed interface treatment techniques. Finally, a study of biphasic mixture behavior is carried out with reference to phase inversion phenomena

This thesis deals with a theoretical and an experimental exploration of squeeze film levitation (SFL) of light objects. The investigations aimed to find the important design parameters controlling this levitation mechanism and also to suggest an alternative way to implement SFL. The study, through computer modelling and experimental validation, focused on Poisson's contraction effect for generating SFL. A finite element model (ANSYS) was verified by experimental testing of five different plate designs. Each plate was subjected to a uniaxial plain stress by an arrangement of two hard piezoelectric actuators (PZT) bonded to the bottom of the plate and driven with DC or AC voltages. It was observed that pulsation of a dimple or crest shaped elastic deformation along the longitudinal axis in the central area of the plate was created because of Poisson's contraction. This Poisson's effect generated the squeeze-film between the plate and the levitated object. The separation distance between a floating lightweight object and the plate was analysed using computational fluid dynamics (ANSYS CFX) through creation of a modelling model for the air-film entrapped between the two interacting surfaces - a typical three-dimensional fluid-solid interaction system (FSI). Additionally, the levitation distance has been experimentally measured by a Laser Sensor. A satisfactory agreement has been found between model predictions and experimental results. Two levitation systems, one based on a horn transducer (Langevin type) and the other one in the form of a plain rectangular plate made of Aluminium and firmly fastened at both ends with a surface-mounted piezoelectric actuator, were compared in this thesis. Both devices were based on SFL mechanism. Evidently, the performances of both designs were greatly influenced by the design structure and in particular by the driving plate characteristics such as plate size and geometry as well as the driving boundary conditions. To this end, physical experiments were carried out and it was found that the device utilising horn-type transducer yields better levitation performance. Ultimately, the research explained the confusion between three approaches to non-contact levitation through literature review and also pointed out some essential parameters like piezoelectric actuators location, material of the driving structure, coupled-field between the actuators and the driving structure and the fluid-solid interface that was existed between the excited plate and the levitated object.

O método de Floquet foi utilizado para estudar a transição da esteira de dois cilindros posicionados lado a lado. Quando os cilindros estão afastados, cada cilindro possui sua própria esteira e para espaçamentos grandes (T/D > 4) a principal influência do cilindro ao lado é a sincronização das esteiras em anti-fase. A esteira em anti-fase corresponde a um escoamento simétrico em relação à linha central, igualmente espaçada em relação aos dois cilindros. Para espaçamentos com T/D < 3, a esteira em fase visualizada por Williamson (1985) é observada. A análise de Floquet da esteira em anti-fase simétrica encontrou um modo que quebra a simetria e induz a transição para esteira em fase. Reduzindo o espaçamento diversos modos complexos apresentam multiplicadores de Floquet que crescem e se tornam dominantes para T/D = 2 coincidindo com o surgimento de esteira irregular. A aproximação dos cilindros também afeta o surgimento de tridimensionalidades no escoamento. Os modos A e B de Williamson (1996b) observados para o cilindro isolado ainda existem, mas observa-se uma supressão do modo B. Por outro lado, o modo A se torna instável para Re menores. A esteira em fase também apresenta os modos A e B, mas neste caso, o modo B se torna instável para Re menores que o observado por Williamson. / The Floquet method was used to study the transition in the wake of two side-by-side circular cylinders. When the cylinders are far apart, each cylinder presents its own wake and for large spacing (T/D > 4) the main influence from the neighboring cylinder is an anti-phase synchronization of the wakes. The antiphase wake corresponds to a symmetric flow in relation to the mid-line between cylinders. When T/D < 3 the in-phase wake described by Williamson (1985) is also observed. The Floquet analysis of the symmetric anti-phase wake found an unstable mode that breaks the symmetry and triggers the transition to in-phase vortex-shedding. As the spacing is further reduced, several complex modes have Floquet multipliers that grow with the decreasing spacing and when T/D = 2 they become dominant suggesting an important role on the irregular wake regime. The presence of the neighboring cylinder also affects the transition to three-dimensional wake. Modes A and B described by Williamson (1996b) for the isolated circular cylinder are still observed but mode B is suppressed as the spacing decreases. On the other hand, mode A becomes unstable for lower Re numbers. The in-phase wake also presents modes A and B but in this case, mode B becomes unstable for lower Re numbers when compared to the single circular cylinder.

This thesis is aimed to calculate the effective permeability tensor and to simulate the fluid flow in naturally fractured reservoirs. This requires an understanding of the mechanisms of fluid flow in naturally fractured reservoirs and the detailed properties of individual fractures and matrix porous media. This study has been carried out to address the issues and difficulties faced by previous methods; to establish possible answers to minimise the difficulties; and hence, to improve the efficiency of reservoir simulation through the use of properties of individual fractures. The methodology used in this study combines several mathematical and numerical techniques like the boundary element method, periodic boundary conditions, and the control volume mixed finite element method. This study has contributed to knowledge in the calculation of the effective permeability and simulation of fluid flow in naturally fractured reservoirs through the development of two algorithms. The first algorithm calculates the effective permeability tensor by use of properties of arbitrary oriented fractures (location, size and orientation). It includes all multi-scaled fractures and considers the appropriate method of analysis for each type of fracture (short, medium and long). In this study a characterisation module which provides the detail information for individual fractures is incorporated. The effective permeability algorithm accounts for fluid flows in the matrix, between the matrix and the fracture and disconnected fractures on effective permeability. It also accounts for the properties of individual fractures in calculation of the effective permeability tensor. The second algorithm simulates flow of single-phase fluid in naturally fractured reservoirs by use of the effective permeability tensor. This algorithm takes full advantage of the control volume discretisation technique and the mixed finite element method in calculation of pressure and fluid flow velocity in each grid block. It accounts for the continuity of flux between the neighbouring blocks and has the advantage of calculation of fluid velocity and pressure, directly from a system of first order equations (Darcy???s law and conservation of mass???s law). The application of the effective permeability tensor in the second algorithm allows us the simulation of fluid flow in naturally fractured reservoirs with large number of multi-scale fractures. The fluid pressure and velocity distributions obtained from this study are important and can considered for further studies in hydraulic fracturing and production optimization of NFRs.

We present predictions for the statistical error due to finite sampling in the presence of thermal fluctuations in molecular simulation algorithms. Specifically, we present predictions for the error dependence on hydrodynamic parameters and the number of samples taken. Expressions for the common hydrodynamic variables of interest such as flow velocity, temperature, density, pressure, shear stress and heat flux are derived using equilibrium statistical mechanics. Both volume-averaged and surface-averaged quantities are considered. Comparisons between theory and computations using direct simulation Monte Carlo for dilute gases, and molecular dynamics for dense fluids, show that the use of equilibrium theory provides accurate results. / Singapore-MIT Alliance (SMA)

The objective of this dissertation is to investigate the non-Darcy flow effects on field base data by considering gas viscosity, gas deviation factor and gas density as variables. To achieve it, different correlations from the literature and field data have been combined to Sawyer-Brown Method, thus a contribution has been achieved. Production history of selected gas field has been implemented to a numerical simulator. To find out non-Darcy effects quantitatively, Darcy flow conditions have also been run in the simulator for each scenario in addition to non-Darcy flow correlation runs. Extracted data from simulation runs have been analyzed on the basis of Sawyer-Brown Method by introducing several correlations to consider gas viscosity, gas deviation factor and gas density as variables. Engineering and scientific research on non-Darcy flow is still being conducted in order for better understanding the nonlinear flow behavior of fluids through porous media. The deviations from Darcy&rsquo / s Law are attributed to the occurrence of all or alternating combinations of factors that can be categorized as the anisotropy of porosity and permeability, multi-phase flow of fluids in varying phases, magnitude of pressure drop and the subsequent phase change in fluids, and the change in flow regime at elevated rates of flow in porous media. Throughout this dissertation, the factors causing deviations from Darcy flow behavior have been investigated.

Problems involving coupled thermo-poro-chemo-mechanical processes are of great importance in geothermal and petroleum reservoir systems. In particular, economic power production from enhanced geothermal systems, effective water-flooding of petroleum reservoirs, and stimulation of gas shale reservoirs are significantly influenced by coupled processes. During such procedures, stress state in the reservoir is changed due to variation in pore fluid pressure and temperature. This can cause deformation and failure of weak planes of the formation with creation of new fractures, which impacts reservoir response. Incorporation of geomechanical factor into engineering analyses using fully coupled geomechanics-reservoir flow modeling exhibits computational challenges and numerical difficulties. In this study, we develop and apply efficient numerical models to solve 3D injection/extraction geomechanics problems formulated within the framework of thermo-poro-mechanical theory with reactive flow. The models rely on combining Displacement Discontinuity (DD) Boundary Element Method (BEM) and Finite Element Method (FEM) to solve the governing equations of thermo-poro-mechanical processes involving fracture/reservoir matrix. The integration of BEM and FEM is accomplished through direct and iterative procedures. In each case, the numerical algorithms are tested against a series of analytical solutions. 3D study of fluid injection and extraction into the geothermal reservoir illustrates that thermo-poro-mechanical processes change fracture aperture (fracture conductivity) significantly and influence the fluid flow. Simulations that consider joint stiffness heterogeneity show development of non-uniform flow paths within the crack. Undersaturated fluid injection causes large silica mass dissolution and increases fracture aperture while supersaturated fluid causes mineral precipitation and closes fracture aperture. Results show that for common reservoir and injection conditions, the impact of fully developed thermoelastic effect on fracture aperture tend to be greater compare to that of poroelastic effect. Poroelastic study of hydraulic fracturing demonstrates that large pore pressure increase especially during multiple hydraulic fracture creation causes effective tensile stress at the fracture surface and shear failure around the main fracture. Finally, a hybrid BEFEM model is developed to analyze stress redistribution in the overburden and within the reservoir during fluid injection and production. Numerical results show that fluid injection leads to reservoir dilation and induces vertical deformation, particularly near the injection well. However, fluid withdrawal causes reservoir to compact. The Mandel-Cryer effect is also successfully captured in numerical simulations, i.e., pore pressure increase/decrease is non-monotonic with a short time values that are above/below the background pore pressure.

No description available.