Feature-oriented grid topology design for aerospace configurations

The last few years have seen a substantial increase in the geometric complexity for 3D flow simulation. In this paper we describe the challenges in generating computation grids for 3D aerospace configuations and demonstrate the progress made to eventually achieve a push button technology for CAD to visualized flow. Special emphasis is given to the interfacing from the grid generator to the flow solver by semi-automatic generation of boundary conditions during the grid generation process. In this regard, once a grid has been generated, push button technology of most commercial flow solvers has been achieved. This will be demonstrated by the ad hoc simulation for the Hopper configuration.

Simulação e análise do escoamento no habitáculo de um automóvel

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

Otimização de vazão de poços injetores em projeto de injeção de água utilizando simulação por linhas de fluxo

Until the early 90s, the simulation of fluid flow in oil reservoir basically used the numerical technique of finite differences. Since then, there was a big development in simulation technology based on streamlines, so that nowadays it is being used in several cases and it can represent the physical mechanisms that influence the fluid flow, such as compressibility, capillarity and gravitational segregation. Streamline-based flow simulation is a tool that can help enough in waterflood project management, because it provides important information not available through traditional simulation of finite differences and shows, in a direct way, the influence between injector well and producer well. This work presents the application of a methodology published in literature for optimizing water injection projects in modeling of a Brazilian Potiguar Basin reservoir that has a large number of wells. This methodology considers changes of injection well rates over time, based on information available through streamline simulation. This methodology reduces injection rates in wells of lower efficiency and increases injection rates in more efficient wells. In the proposed model, the methodology was effective. The optimized alternatives presented higher oil recovery associated with a lower water injection volume. This shows better efficiency and, consequently, reduction in costs. Considering the wide use of the water injection in oil fields, the positive outcome of the modeling is important, because it shows a case study of increasing of oil recovery achieved simply through better distribution of water injection rates

Análise comparativa entre simuladores de linha de fluxo e diferenças finitas para um reservatório de petróleo submetido a injeção de água

Water injection is the most widely used method for supplementary recovery in many oil fields due to various reasons, like the fact that water is an effective displacing agent of low viscosity oils, the water injection projects are relatively simple to establish and the water availability at a relatively low cost. For design of water injection projects is necessary to do reservoir studies in order to define the various parameters needed to increase the effectiveness of the method. For this kind of study can be used several mathematical models classified into two general categories: analytical or numerical. The present work aims to do a comparative analysis between the results presented by flow lines simulator and conventional finite differences simulator; both types of simulators are based on numerical methods designed to model light oil reservoirs subjected to water injection. Therefore, it was defined two reservoir models: the first one was a heterogeneous model whose petrophysical properties vary along the reservoir and the other one was created using average petrophysical properties obtained from the first model. Comparisons were done considering that the results of these two models were always in the same operational conditions. Then some rock and fluid parameters have been changed in both models and again the results were compared. From the factorial design, that was done to study the sensitivity analysis of reservoir parameters, a few cases were chosen to study the role of water injection rate and the vertical position of wells perforations in production forecast. It was observed that the results from the two simulators are quite similar in most of the cases; differences were found only in those cases where there was an increase in gas solubility ratio of the model. Thus, it was concluded that in flow simulation of reservoirs analogous of those now studied, mainly when the gas solubility ratio is low, the conventional finite differences simulator may be replaced by flow lines simulator the production forecast is compatible but the computational processing time is lower.

Simulador de escoamento multifásico em poços de petróleo (SEMPP)

The multiphase flow occurrence in the oil and gas industry is common throughout fluid path, production, transportation and refining. The multiphase flow is defined as flow simultaneously composed of two or more phases with different properties and immiscible. An important computational tool for the design, planning and optimization production systems is multiphase flow simulation in pipelines and porous media, usually made by multiphase flow commercial simulators. The main purpose of the multiphase flow simulators is predicting pressure and temperature at any point at the production system. This work proposes the development of a multiphase flow simulator able to predict the dynamic pressure and temperature gradient in vertical, directional and horizontal wells. The prediction of pressure and temperature profiles was made by numerical integration using marching algorithm with empirical correlations and mechanistic model to predict pressure gradient. The development of this tool involved set of routines implemented through software programming Embarcadero C++ Builder® 2010 version, which allowed the creation of executable file compatible with Microsoft Windows® operating systems. The simulator validation was conduct by computational experiments and comparison the results with the PIPESIM®. In general, the developed simulator achieved excellent results compared with those obtained by PIPESIM and can be used as a tool to assist production systems development

Simulação computacional da interação fluido-estrutura em bombas de cavidades progressivas

The pumping through progressing cavities system has been more and more employed in the petroleum industry. This occurs because of its capacity of elevation of highly viscous oils or fluids with great concentration of sand or other solid particles. A Progressing Cavity Pump (PCP) consists, basically, of a rotor - a metallic device similar to an eccentric screw, and a stator - a steel tube internally covered by a double helix, which may be rigid or deformable/elastomeric. In general, it is submitted to a combination of well pressure with the pressure generated by the pumping process itself. In elastomeric PCPs, this combined effort compresses the stator and generates, or enlarges, the clearance existing between the rotor and the stator, thus reducing the closing effect between their cavities. Such opening of the sealing region produces what is known as fluid slip or slippage, reducing the efficiency of the PCP pumping system. Therefore, this research aims to develop a transient three-dimensional computational model that, based on single-lobe PCP kinematics, is able to simulate the fluid-structure interaction that occurs in the interior of metallic and elastomeric PCPs. The main goal is to evaluate the dynamic characteristics of PCP s efficiency based on detailed and instantaneous information of velocity, pressure and deformation fields in their interior. To reach these goals (development and use of the model), it was also necessary the development of a methodology for generation of dynamic, mobile and deformable, computational meshes representing fluid and structural regions of a PCP. This additional intermediary step has been characterized as the biggest challenge for the elaboration and running of the computational model due to the complex kinematic and critical geometry of this type of pump (different helix angles between rotor and stator as well as large length scale aspect ratios). The processes of dynamic generation of meshes and of simultaneous evaluation of the deformations suffered by the elastomer are fulfilled through subroutines written in Fortan 90 language that dynamically interact with the CFX/ANSYS fluid dynamic software. Since a structural elastic linear model is employed to evaluate elastomer deformations, it is not necessary to use any CAE package for structural analysis. However, an initial proposal for dynamic simulation using hyperelastic models through ANSYS software is also presented in this research. Validation of the results produced with the present methodology (mesh generation, flow simulation in metallic PCPs and simulation of fluid-structure interaction in elastomeric PCPs) is obtained through comparison with experimental results reported by the literature. It is expected that the development and application of such a computational model may provide better details of the dynamics of the flow within metallic and elastomeric PCPs, so that better control systems may be implemented in the artificial elevation area by PCP

Incorporação do vínculo de suavidade no ajuste de histórico de reservatórios de petróleo

The history match procedure in an oil reservoir is of paramount importance in order to obtain a characterization of the reservoir parameters (statics and dynamics) that implicates in a predict production more perfected. Throughout this process one can find reservoir model parameters which are able to reproduce the behaviour of a real reservoir.Thus, this reservoir model may be used to predict production and can aid the oil file management. During the history match procedure the reservoir model parameters are modified and for every new set of reservoir model parameters found, a fluid flow simulation is performed so that it is possible to evaluate weather or not this new set of parameters reproduces the observations in the actual reservoir. The reservoir is said to be matched when the discrepancies between the model predictions and the observations of the real reservoir are below a certain tolerance. The determination of the model parameters via history matching requires the minimisation of an objective function (difference between the observed and simulated productions according to a chosen norm) in a parameter space populated by many local minima. In other words, more than one set of reservoir model parameters fits the observation. With respect to the non-uniqueness of the solution, the inverse problem associated to history match is ill-posed. In order to reduce this ambiguity, it is necessary to incorporate a priori information and constraints in the model reservoir parameters to be determined. In this dissertation, the regularization of the inverse problem associated to the history match was performed via the introduction of a smoothness constraint in the following parameter: permeability and porosity. This constraint has geological bias of asserting that these two properties smoothly vary in space. In this sense, it is necessary to find the right relative weight of this constrain in the objective function that stabilizes the inversion and yet, introduces minimum bias. A sequential search method called COMPLEX was used to find the reservoir model parameters that best reproduce the observations of a semi-synthetic model. This method does not require the usage of derivatives when searching for the minimum of the objective function. Here, it is shown that the judicious introduction of the smoothness constraint in the objective function formulation reduces the associated ambiguity and introduces minimum bias in the estimates of permeability and porosity of the semi-synthetic reservoir model

Sistema de controle de escoamento de gás natural em dutos de distribuição com detecção de vazamentos utilizando redes neurais

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Análise da forma da partícula na aplicação do Método dos Elementos Discretos (DEM)

O presente trabalho analisa o efeito da forma da partícula por meio do método dos Elementos Discretos (DEM). São analisadas quatro partículas com formas e esfericidades diferentes – uma esférica e três não-esféricas. O ângulo de repouso é o principal instrumento de avaliação dos escoamentos. A análise inicia com a calibração da partícula esférica usando dados da literatura, em seguida, estes parâmetros calibrados são empregados para a simulação dos escoamentos das partículas não-esféricas. São feitas comparações dos esforços computacionais e, estas informações, são usadas para verificar as vantagens que a partícula esférica oferece sobre as outras três formas. Nesse cenário, procedimentos foram desenvolvidos para ajudar no processo de calibração dos ângulos de repousos, tendo como base, o conhecimento da sensibilidade de alguns parâmetros DEM. Os resultados mostram a influência das partículas não-esféricas e, principalmente, que é possível obter com a partícula esférica, escoamentos similares aos das partículas não-esféricas, com ganho computacional.

Modelagem de motor de combustão internae simulação do processo de queima de combustível

The internal combustion engine is a heat engine widely used in the automotive industry. In order to better understand its behavior many models in the literature have been proposed in the last years. The 0-D thermodynamic model is a fairly simple tool but it is very useful to understand the phenomenon of combustion inside the chamber of internal combustion engines. In the first phase of this work, an extensive literature review was made in order to get information about this kind of analysis and, after this, apply them in a model able to calculate the instantaneous temperature and pressure in one zone of the combustion chamber of a diesel engine. Therefore some considerations were made with the aim of increasing the accuracy of the model in predicting the correct behavior of the engine, adding the combined effects of heat transfer, leakage and injection. In the second phase, the goal was to study the internal flow of a three-dimensional model of an internal combustion engine. In order to achieve this goal the software Solidworks was used to create the geometries of an engine and the suite of softwares Ansys was used to create the moving mesh (ICEM CFD and CFX-Pre) and to solve the CFD problem (Ansys CFX code). The model was able to perform the air flow simulation during the four-stroke cycle of an engine: admission, compression, expansion and exhaust. The results obtained from both models were suitable and they open a new range of possibilities for future researches on the field

Simulação numérica de fluxo de águas subterrâneas na bacia do Rio Corrente (BA)

The numeric simulation is an important tool applied in understanding the dynamics of groundwater flow. In a hydrogeological model the processes responsible for groundwater flow are described by numerical formulations that allow the simplification, representation and understanding of the dynamics of the Aquifer System. In this work, a steady state groundwater flow simulation of Urucuia Aquifer System (UAS) part of the Corrente river basin was conducted, using the finite element method through software FEFLOW, to understand the dynamics of groundwater flow and quantify the hydrologic balance. The aquifer system Urucuia lodges in the São Francisco hydrogeological province and corresponds to a set of interconnected aquifers that occur in rocks from Urucuia group in the Urucuia sub-basin described by Campos e Dardenne (1997). The system is a porous media one, in a shape of a thick table mountain, consisting essentially of sandstones. The Corrente river basin is located in UAS in Western State of Bahia and it's one of the main units to maintaining permanent flow (Q95) and average natural flow of the São Francisco river. The simulation performed in this work obtained the following results for the modelled region: horizontal hydraulic conductivity of 3 x 10-4 m/s and vertical one 6 x 10-5 m/s; maximum recharge of 345 mm and minimum of 85 mm/a. It was concluded that: (1) regional groundwater flow has eastbound; with an exception of the extreme northeast portion, where the flow has opposite direction; (2) there are smaller water side dividers with an approximate direction EW, that guide the flow of water to the drainage that cut the aquifer; and (3) the UAS at Corrente river basin can be understood as a free regional aquifer system, isotropic and homogeneous. Regionally, the small lithological variations present in the Urucuia group can be neglected and do not exhibit significant influences on the dynamics of ground water flow

Análise da influência de variantes do método de solução de escoamentos utilizando os métodos de elementos finitos com compressibilidade artificial e pseudo-características

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

A numerical study of temperature effects on hot wire measurements inside turbulent channel flows

A way to investigate turbulence is through experiments where hot wire measurements are performed. Analysis of the in turbulence of a temperature gradient on hot wire measurements is the aim of this thesis work. Actually - to author's knowledge - this investigation is the first attempt to document, understand and ultimately correct the effect of temperature gradients on turbulence statistics. However a numerical approach is used since instantaneous temperature and streamwise velocity fields are required to evaluate this effect. A channel flow simulation at Re_tau = 180 is analyzed to make a first evaluation of the amount of error introduced by temperature gradient inside the domain. Hot wire data field is obtained processing the numerical flow field through the application of a proper version of the King's law, which connect voltage, velocity and temperature. A drift in mean streamwise velocity profile and rms is observed when temperature correction is performed by means of centerline temperature. A correct mean velocity pro�le is achieved correcting temperature through its mean value at each wall normal position, but a not negligible error is still present into rms. The key point to correct properly the sensed velocity from the hot wire is the knowledge of the instantaneous temperature field. For this purpose three correction methods are proposed. At the end a numerical simulation at Re_tau =590 is also evaluated in order to confirm the results discussed earlier.

Studio e ottimizzazione della ipersostentazione del velivolo V-173

Studio della ipersostentazione di un velivolo non convenzionale degli anni ‘40, il Vought V-173 Flying Pancake. Per realizzare l’intero lavoro si sono adoperati dei software di fluidodinamica computazionale come Flow Simulation di SolidWorks e il programma JavaFoil, inoltre, si è scelto l’ambiente Matlab di supporto ai calcoli e per rappresentare i risultati ottenuti a seguito delle simulazioni, realizzando script e funzioni per l’approssimazione di questi e la loro successiva raffigurazione grafica. In particolar modo, a partire dal modello tridimensionale in SolidWorks del V-173 si sono ricreate le curve CL-α e CD-α a partire dai punti ottenuti dalle simulazioni per diverse configurazioni del velivolo. In una prima fase si è valutata l’aerodinamica del velivolo ‘pulito’ senza la presenza delle eliche e di ipersostentatori, successivamente si sono seguite due strade diverse per valutare il comportamento del velivolo: nel primo caso si sono eseguiti studi dell’aerodinamica del Pancake in presenza degli ipersostentatori già presenti sul velivolo (i plain flap), proponendo soluzioni alternative per migliorare l’ipersostentazione del V-173 tramite spillamenti (slots) e diverse configurazioni di vortex generator che energizzassero lo strato limite e ottimizzassero le prestazioni del velivolo con particolare attenzione alla fase di atterraggio. In secondo luogo, tenendo in considerazione che il Pancake è un aeroplano bielica, si è voluta studiare l’influenza delle due eliche sulla sua aerodinamica: dopo aver riprodotto nel modo più verosimile entrambe le eliche utilizzando SolidWorks si è fatto uno studio di massima ricavando risultati che potessero indicare la compatibilità tra elica e velivolo a seguito dei risultati sperimentali ottenuti con Flow Simulation.