Desenvolvimento de um método para correção de curvas de desempenho em bombas centrífugas submersas operando com fluidos viscosos

Electrical Submersible Pump (ESP) is used as an artificial lift technique. However, pumping viscous oil is generally associated with low Reynolds number flows. This condition leads to a performance degradation respect to the performance expected from the regular operation with water that most of the centrifugal pumps are originally designed for. These issues are considered in this investigation through a numerical study of the flow in two different multistage, semi-axial type ESPs. This investigation is carried out numerically using a Computational Fluid Dynamics (CFD) package, where the transient RANS equations are solved numerically. The turbulence is modeled using the SST model. Head curves for several operating conditions are compared with manufacturer’s curves and experimental data for a three-stage ESP, showing good agreement for a wide range of fluid viscosities and rotational speeds. Dimensionless numbers (n, n, n e Re) are used to investigate performance degradation of the ESPs. In addition, flow phenomena through the impellers of the ESPs are investigated using flow field from numerical results. Results show that performance degradation is directly related to rotational Reynolds number, Re. In addition, it was verified that performance degradation occurs for constant normalized specific speedn, which shows that performance degradation occurs similarly for different centrifugal pumps. Moreover, experimental data and numerical results agreed with a correlation from literature between head and flow correction factors proposed by Stepanoff (1967). A definition of modified Reynolds number was proposed and relates the head correction factor to viscosity. A correlation between head correction factor and the modified Reynolds number was proposed, which agreed well with numerical and experimental data. Then, a method to predict performance degradation based on the previous correlations was proposed. This method was compared with others from literature. In general, results and conclusions from this work can also be useful to bring more information about the flow of highly viscous fluids in pumps, especially in semi-axial, multistage ESPs.

Medidas de concentração de radônio proveniente de argamassas de cimento portland, gesso e fosfogesso

Portland cement being very common construction material has in its composition the natural gypsum. To decrease the costs of manufacturing, the cement industry is substituting the gypsum in its composition by small quantities of phosphogypsum, which is the residue generated by the production of fertilizers and consists essentially of calcium dihydrate and some impurities, such as fluoride, metals in general, and radionuclides. Currently, tons of phosphogypsum are stored in the open air near the fertilizer industries, causing contamination of the environment. The 226 Ra present in these materials, when undergoes radioactive decay, produces the 222Rn gas. This radioactive gas, when inhaled together with its decay products deposited in the lungs, produces the exposure to radiation and can be a potential cause of lung cancer. Thus, the objective of this study was to measure the concentration levels of 222Rn from cylindrical samples of Portland cement, gypsum and phosphogypsum mortar from the state of Paraná, as well as characterizer the material and estimate the radon concentration in an environment of hypothetical dwelling with walls covered by such materials. Experimental setup of 222Rn activity measurements was based on AlphaGUARD detector (Saphymo GmbH). The qualitative and quantitative analysis was performed by gamma spectrometry and EDXRF with Au and Ag targets tubes (AMPTEK), and Mo target (ARTAX) and mechanical testing with x- ray equipment (Gilardoni) and the mechanical press (EMIC). Obtained average values of radon activity from studied materials in the air of containers were of 854 ± 23 Bq/m3, 60,0 ± 7,2 Bq/m3 e 52,9 ± 5,4 Bq/m3 for Portland cement, gypsum and phosphogypsum mortar, respectively. These results extrapolated into the volume of hypothetical dwelling of 36 m3 with the walls covered by such materials were of 3366 ± 91 Bq/m3, 237 ± 28 Bq/m3 e 208 ± 21 Bq/m3for Portland cement, gypsum and phosphogypsum mortar, respectively. Considering the limit of 300 Bq/m3 established by the ICRP, it could be concluded that the use of Portland cement plaster in dwellings is not secure and requires some specific mitigation procedure. Using the results of gamma spectrometry there were calculated the values of radium equivalent activity concentrations (Raeq) for Portland cement, gypsum and phosphogypsum mortar, which were obtained equal to 78,2 ± 0,9 Bq/kg; 58,2 ± 0,9 Bq/kg e 68,2 ± 0,9 Bq/kg, respectively. All values of radium equivalent activity concentrations for studied samples are below the maximum level of 370 Bq/kg. The qualitative and quantitative analysis of EDXRF spectra obtained with studied mortar samples allowed to evaluate quantitate and the elements that constitute the material such as Ca, S, Fe, and others.

Protocolo de roteamento para redes oportunistas utilizando otimização por enxame de partículas

In this research work, a new routing protocol for Opportunistic Networks is presented. The proposed protocol is called PSONET (PSO for Opportunistic Networks) since the proposal uses a hybrid system composed of a Particle Swarm Optimization algorithm (PSO). The main motivation for using the PSO is to take advantage of its search based on individuals and their learning adaptation. The PSONET uses the Particle Swarm Optimization technique to drive the network traffic through of a good subset of forwarders messages. The PSONET analyzes network communication conditions, detecting whether each node has sparse or dense connections and thus make better decisions about routing messages. The PSONET protocol is compared with the Epidemic and PROPHET protocols in three different scenarios of mobility: a mobility model based in activities, which simulates the everyday life of people in their work activities, leisure and rest; a mobility model based on a community of people, which simulates a group of people in their communities, which eventually will contact other people who may or may not be part of your community, to exchange information; and a random mobility pattern, which simulates a scenario divided into communities where people choose a destination at random, and based on the restriction map, move to this destination using the shortest path. The simulation results, obtained through The ONE simulator, show that in scenarios where the mobility model based on a community of people and also where the mobility model is random, the PSONET protocol achieves a higher messages delivery rate and a lower replication messages compared with the Epidemic and PROPHET protocols.

Ambiente para aquisição e processamento de sinal eletromiográfico de superfície

In this work, a platform to the conditioning, digitizing, visualization and recording of the EMG signals was developed. After the acquisition, the analysis can be done by signal processing techniques. The platform consists of two modules witch acquire electromyography (EMG) signals by surface electrodes, limit the interest frequency band, filter the power grid interference and digitalize the signals by the analogue-to- digital converter of the modules microcontroller. Thereby, the data are sent to the computer by the USB interface by the HID specification, displayed in real-time in graphical form and stored in files. As processing resources was implemented the operations of signal absolute value, the determination of effective value (RMS), Fourier analysis, digital filter (IIR) and the adaptive filter. Platform initial tests were performed with signal of lower and upper limbs with the aim to compare the EMG signal laterality. The open platform is intended to educational activities and academic research, allowing the addition of other processing methods that the researcher want to evaluate or other required analysis.

Uma janela tropical: análise do desempenho luminoso de ambientes com aberturas sombreadas para o clima de Natal/RN

This master thesis introduces assessment procedures of daylighting performance in office rooms with shaded opening, recommendations for Natal-RN (Latitude 05,47' S, Longitude 35,11' W). The studies assume the need of window exterior shading in hot and humid climate buildings. The daylighting performance analyses are based on simulated results for three levels of illuminance (300,500 e 1000 lux) between 08h00 e 16h00, in rooms with 2,80 m height, 6 m large and 4 m, 6 m e 8 m depths, with a centered single opening, window wall ratio (20%, 40% e 60%), four orientations (North, East, South and West), and two types of sky (clear and partially cloudy). The sky characteristics were statistically determined based on hourly data from INPE-CRN solar and daylighting weather station. The lighting performance is resulted from dynamic computer simulation of 72 models using Troplux 3.12. The simulation results were assessed using a new parameter to quantify the use of interior daylighting, the useful percentage of daylight (PULN), which corresponds to the time fraction with satisfactory light, in accordance with the illuminance design. The passive zone depths are defined based on the PULN. Despite the failures of illuminance data from the weather station, the analyses ratified the high potential of daylighting for shaded rooms. The most influential variables on the lighting performance are the opening size and the illuminance of design, while the orientation is a little influential

Implementação paralelizada de métodos de resolução de sistemas algébricos na simulação de reservatórios de gás

Desde a década de 1960, devido à pertinência para a indústria petrolífera, a simulação numérica de reservatórios de petróleo tornou-se uma ferramenta usual e uma intensa área de pesquisa. O principal objetivo da modelagem computacional e do uso de métodos numéricos, para a simulação de reservatórios de petróleo, é o de possibilitar um melhor gerenciamento do campo produtor, de maneira que haja uma maximização na recuperação de hidrocarbonetos. Este trabalho tem como objetivo principal paralelizar, empregando a interface de programação de aplicativo OpenMP (Open Multi-Processing), o método numérico utilizado na resolução do sistema algébrico resultante da discretização da equação que descreve o escoamento monofásico em um reservatório de gás, em termos da variável pressão. O conjunto de equações governantes é formado pela equação da continuidade, por uma expressão para o balanço da quantidade de movimento e por uma equação de estado. A Equação da Difusividade Hidráulica (EDH), para a variável pressão, é obtida a partir deste conjunto de equações fundamentais, sendo então discretizada pela utilização do Método de Diferenças Finitas, com a escolha por uma formulação implícita. Diferentes testes numéricos são realizados a fim de estudar a eficiência computacional das versões paralelizadas dos métodos iterativos de Jacobi, Gauss-Seidel, Sobre-relaxação Sucessiva, Gradientes Conjugados (CG), Gradiente Biconjugado (BiCG) e Gradiente Biconjugado Estabilizado (BiCGStab), visando a uma futura aplicação dos mesmos na simulação de reservatórios de gás. Ressalta-se que a presença de heterogeneidades na rocha reservatório e/ou às não-linearidades presentes na EDH para o escoamento de gás aumentam a necessidade de métodos eficientes do ponto de vista de custo computacional, como é o caso de estratégias usando OpenMP.

Aplicações da teoria da percolação à modelagem e simulação de reservatórios de petróleo

In this thesis we study some problems related to petroleum reservoirs using methods and concepts of Statistical Physics. The thesis could be divided percolation problem in random multifractal support motivated by its potential application in modelling oil reservoirs. We develped an heterogeneous and anisotropic grid that followin two parts. The first one introduce a study of the percolations a random multifractal distribution of its sites. After, we determine the percolation threshold for this grid, the fractal dimension of the percolating cluster and the critical exponents ß and v. In the second part, we propose an alternative systematic of modelling and simulating oil reservoirs. We introduce a statistical model based in a stochastic formulation do Darcy Law. In this model, the distribution of permeabilities is localy equivalent to the basic model of bond percolation

Simulação de transporte eletrônico em dispositivos unimoleculares baseados em indicadores de pH

Nas últimas décadas, diversos pesquisadores têm tentado empregar moléculas em dispositivos eletrônicos de nanoescala. Por este motivo, diferentes parâmetros eletro/ópticos, que regem o transporte eletrônico em moléculas orgânicas, precisam ser analisados. Neste trabalho foi desenvolvido um estudo de transporte de carga para o composto Vermelho de Propila, popularmente utilizado como indicador de pH. A motivação para estudá-lo resulta de sua estrutura constituída por subunidades doadora-aceitadora, acopladas via grupo azo (N=N), uma característica bem conhecida em retificadores moleculares. A metodologia utilizada para tratar o sistema em equilíbrio é baseada em métodos de Mecânica Molecular e Hartree-Fock. Sendo que, para simular o sistema em não-equilíbrio, foi empregado o formalismo de Landauer-Büttiker. Através desses métodos, as curvas características do sistema molecular foram traçadas e comparadas. O resultado da comparação permitiu explicar os fenômenos que regem o transporte eletrônico na nanoestrutura. Além disso, foram analisados os efeitos de contatos metálicos, ligados a molécula na presença de campo elétrico externo.

Estabilidade de métodos numéricos LCR para a simulação de escoamento

Pós-graduação em Matematica Aplicada e Computacional - FCT

Nuevos métodos para la obtención de soluciones iniciales en el problema de transporte

Una vez formulado el modelo matemático del problema de transporte, el siguiente paso consiste en resolver el modelo, es decir, en obtener los mejores valores númericos para las variables de decisión. La forma en que se obtengan estos valores depende del tipo específico del modelo matemático utilizado, Por tanto, una vez definido el modelo, se podrá elegir un método apropiado para resolverlo.

Desenvolvimento de um método espectronodal livre de erros de truncamento espacial para problemas adjuntos de transporte de partículas neutras monoenergéticas na formulação de ordenadas discretas em geometria unidimensional

Um método numérico nodal livre de erros de truncamento espacial é desenvolvido para problemas adjuntos de transporte de partículas neutras monoenergéticas em geometria unidimensional com fonte fixa na formulação de ordenadas discretas (SN). As incógnitas no método são os fluxos angulares adjuntos médios nos nodos e os fluxos angulares adjuntos nas fronteiras dos nodos, e os valores numéricos gerados para essas quantidades são os obtidos a partir da solução analítica das equações SN adjuntas. O método é fundamentado no uso da convencional equação adjunta SN discretizada de balanço espacial, que é válida para cada nodo de discretização espacial e para cada direção discreta da quadratura angular, e de uma equação auxiliar adjunta não convencional, que contém uma função de Green para os fluxos angulares adjuntos médios nos nodos em termos dos fluxos angulares adjuntos emergentes das fronteiras dos nodos e da fonte adjunta interior. Resultados numéricos são fornecidos para ilustrarem a precisão do método proposto.

Métodos espectronodais para cálculos de transporte de partículas neutras com fonte fixa na formulação de ordenadas discretas e multigrupo de energia

Um método espectronodal é desenvolvido para problemas de transporte de partículas neutras de fonte fixa, multigrupo de energia em geometria cartesiana na formulação de ordenadas discretas (SN). Para geometria unidimensional o método espectronodal multigrupo denomina-se método spectral Greens function (SGF) com o esquema de inversão nodal (NBI) que converge solução numérica para problemas SN multigrupo em geometria unidimensional, que são completamente livre de erros de truncamento espacial para ordem L de anisotropia de espalhamento desde que L < N. Para geometria X; Y o método espectronodal multigrupo baseia-se em integrações transversais das equações SN no interior dos nodos de discretização espacial, separadamente nas direções coordenadas x e y. Já que os termos de fuga transversal são aproximados por constantes, o método nodal resultante denomina-se SGF-constant nodal (SGF-CN), que é aplicado a problemas SN multigrupo de fonte fixa em geometria X; Y com espalhamento isotrópico. Resultados numéricos são apresentados para ilustrar a eficiência dos códigos SGF e SGF-CN e a precisão das soluções numéricas convergidas em cálculos de malha grossa.

Modelagem e simulação do escoamento imiscível em meios porosos fractais descritos pela equação de Kozeny-Carman Generalizada

O presente trabalho trata do escoamento bifásico em meios porosos heterogêneos de natureza fractal, onde os fluidos são considerados imiscíveis. Os meios porosos são modelados pela equação de Kozeny-Carman Generalizada (KCG), a qual relaciona a porosidade com a permeabilidade do meio através de uma nova lei de potência. Esta equação proposta por nós é capaz de generalizar diferentes modelos existentes na literatura e, portanto, é de uso mais geral. O simulador numérico desenvolvido aqui emprega métodos de diferenças finitas. A evolução temporal é baseada em um esquema de separação de operadores que segue a estratégia clássica chamada de IMPES. Assim, o campo de pressão é calculado implicitamente, enquanto que a equação da saturação da fase molhante é resolvida explicitamente em cada nível de tempo. O método de otimização denominado de DFSANE é utilizado para resolver a equação da pressão. Enfatizamos que o DFSANE nunca foi usado antes no contexto de simulação de reservatórios. Portanto, o seu uso aqui é sem precedentes. Para minimizar difusões numéricas, a equação da saturação é discretizada por um esquema do tipo "upwind", comumente empregado em simuladores numéricos para a recuperação de petróleo, o qual é resolvido explicitamente pelo método Runge-Kutta de quarta ordem. Os resultados das simulações são bastante satisfatórios. De fato, tais resultados mostram que o modelo KCG é capaz de gerar meios porosos heterogêneos, cujas características permitem a captura de fenômenos físicos que, geralmente, são de difícil acesso para muitos simuladores em diferenças finitas clássicas, como o chamado fenômeno de dedilhamento, que ocorre quando a razão de mobilidade (entre as fases fluidas) assume valores adversos. Em todas as simulações apresentadas aqui, consideramos que o problema imiscível é bidimensional, sendo, portanto, o meio poroso caracterizado por campos de permeabilidade e de porosidade definidos em regiões Euclideanas. No entanto, a teoria abordada neste trabalho não impõe restrições para sua aplicação aos problemas tridimensionais.

Uma Metodologia de Estudo de Simulação Tridimensional de Escoamento Turbulento Estratificado no Reservatório de Plantas Hidrelétricas.

Uma simulação numérica que leva em conta os efeitos de estratificação e mistura escalar (como a temperatura, salinidade ou substância solúvel em água) é necessária para estudar e prever os impactos ambientais que um reservatório de usina hidrelétrica pode produzir. Este trabalho sugere uma metodologia para o estudo de escoamentos ambientais, principalmente aqueles em que o conhecimento da interação entre a estratificação e mistura pode dar noções importantes dos fenômenos que ocorrem. Por esta razão, ferramentas de simulação numérica 3D de escoamento ambiental são desenvolvidas. Um gerador de malha de tetraedros do reservatório e o modelo de turbulência algébrico baseado no número de Richardson são as principais ferramentas desenvolvidas. A principal dificuldade na geração de uma malha de tetraedros de um reservatório é a distribuição não uniforme dos pontos relacionada com a relação desproporcional entre as escalas horizontais e verticais do reservatório. Neste tipo de distribuição de pontos, o algoritmo convencional de geração de malha de tetraedros pode tornar-se instável. Por esta razão, um gerador de malha não estruturada de tetraedros é desenvolvido e a metodologia utilizada para obter elementos conformes é descrita. A geração de malha superficial de triângulos utilizando a triangulação Delaunay e a construção do tetraedros a partir da malha triangular são os principais passos para o gerador de malha. A simulação hidrodinâmica com o modelo de turbulência fornece uma ferramenta útil e computacionalmente viável para fins de engenharia. Além disso, o modelo de turbulência baseado no número de Richardson leva em conta os efeitos da interação entre turbulência e estratificação. O modelo algébrico é o mais simples entre os diversos modelos de turbulência. Mas, fornece resultados realistas com o ajuste de uma pequena quantidade de parâmetros. São incorporados os modelos de viscosidade/difusividade turbulenta para escoamento estratificado. Na aproximação das equações médias de Reynolds e transporte de escalar é utilizando o Método dos Elementos Finitos. Os termos convectivos são aproximados utilizando o método semi-Lagrangeano, e a aproximação espacial é baseada no método de Galerkin. Os resultados computacionais são comparados com os resultados disponíveis na literatura. E, finalmente, a simulação de escoamento em um braço de reservatório é apresentada.

Simulação numérica da dispersão advecção de pesticidas no solo sob efeito da temperatura.

É apresentado um modelo de dispersão-advecção de evolução unidimensional que simula a lixiviação de pesticidas em lisímetros ou colunas de solo sob efeito da temperatura média diária do perfil do solo. O modelo matemático e todo o seu conjunto de hipóteses será denominado DAPESTE. Nas simulações numéricas do modelo DAPESTE serão utilizados o método dos elementos finitos para a semi-discretização da variável espacial e o método de Eüler regressivo para a discretização da variável temporal. Serão utilizados métodos de elementos finitos apropriados para problemas de dispersão-advecção nos quais o transporte advectivo predomina sobre o transporte dispersivo. O modelo DAPESTE supõe que as difusividades do pesticida nas fases gasosa e aquosa do solo dependem da temperatura média diária do solo a qual varia periodicamente com a profundidade e com o tempo. O coeficiente de dispersão hidrodinâmico do modelo DAPESTE dependerá da temperatura do solo. O coeficiente de partição água-ar, variando com a temperatura, será determinado pela equação de Clausius-Clapeyron. A equação de van?t Ho® será usada para determinar o coeficiente de sorção do pesticida no solo como função da temperatura. Com a equação de Arrhenius será estimado o efeito da temperatura na taxa de degradação do pesticida. Estas relações de dependência entre os parâmetros do modelo e a temperatura do solo auxiliam sobremaneira na compreensão do destino de pesticidas no solo sob diferentes cenários de temperaturas médias diárias, especialmente a meia vida do pesticida e a concentração lixiviada no perfil do solo.