Estudo de um modelo convectivo-difusivo-reativo em combustão no método de elementos finitos

Neste trabalho, apresenta-se um estudo numérico de um modelo convectivo-difusivo-reativo em combustão baseado no Método de Elementos Finitos. Primeiramente, apresenta-se o desenvolvimento das equações de balanço (quantidade de movimento, massa, espécie e energia) que modelam um processo de mistura molecular e reação química, irreversível, de passo único e exotérmica entre duas espécies químicas F (Combustível) e O (Oxidante). Tais espécies reagem e formam um produto P, conforme vFF +vOO ! vPP + calor, onde vF , vO e vP são os coeficientes estequiométricos molares. No modelo, considera-se que a reação é de primeira ordem com respeito a cada um dos reagentes e que a taxa de reação específica segue a cinética de Arrhenius. Em seguida, o modelo é estudado numericamente considerando-se um domínio retangular e condições de contorno do tipo Neumann. Tanto a Técnica das Diferenças Finitas como a Técnica de Elementos Finitos são utilizadas na discretização espacial das equações do modelo. Para a integração no tempo, utiliza-se a método de Runge-Kutta simplificado de três estágios. Os diferentes códigos computacionais obtidos, tanto pela Técnica de Diferenças Finitas como de Elementos Finitos, são comparados frente ao problema de interesse. Observa-se que ambas as técnicas apresentam resultados equivalentes. Além disso, os códigos desenvolvidos são robustos (capazes de lidar com vários conjuntos de parâmetros), de baixo custo e precisos. Por fim, apresenta-se uma revisão do trabalho de Zavaleta [48], no qual obtem-se uma estimativa local do erro na aproximação do problema estudado pela Técnica de Elementos Finitos.

Um modelo matemático baseado em wavelets para análise do método térmico de recuperação de óleo pesado aplicando irradiação eletromagnética

This work proposes a model to investigate the use of a cylindrical antenna used in the thermal method of recovering through electromagnetic radiation of high-viscosity oil. The antenna has a simple geometry, adapted dipole type, and it can be modelled by using Maxwell s equation. The wavelet transforms are used as basis functions and applied in conjunction with the method of moments to obtain the current distribution in the antenna. The electric field, power and temperature distribution are carefully calculated for the analysis of the antenna as electromagnetic heating. The energy performance is analyzed based on thermo-fluid dynamic simulations at field scale, and through the adaptation in the Steam Thermal and Advanced Processes Reservoir Simulator (STARS) by Computer Modelling Group (CMG). The model proposed and the numerical results obtained are stable and presented good agreement with the results reported in the specialized literature

Modelagem do desempenho ambiental dos projetos de exploração e produção aplicando equações estruturais

A legislação ambiental e os principais agentes que se relacionam com a empresa se constituem em fatores exógenos que não podem ser negligenciados ao formular-se e avaliar-se a política ambiental corporativa. As influências exógenas e seus efeitos sobre a gestão ambiental e o gerenciamento de projetos de exploração e produção (E&P) e, por essa via, sobre o desempenho ambiental, foram objetos de estudo desta tese. Embora o desempenho ambiental seja um assunto relevante, a pesquisa sobre esse tema ainda é escassa. Tal carência desponta ainda mais acentuada quando se aborda o desempenho ambiental de projetos na indústria de petróleo e gás. O principal objetivo deste estudo foi avaliar a relação entre a legislação ambiental vigente, as ações de órgãos reguladores, fornecedores, empresas terceirizadas e comunidades locais e o desempenho ambiental dos projetos de E&P na indústria de petróleo e gás e, também, analisar os efeitos do sistema de gestão ambiental e o gerenciamento dos projetos sobre tal desempenho. Na fase abdutiva, foi conduzido um estudo de caso com abordagem qualitativa em uma grande empresa brasileira do setor de petróleo e gás, na fase dedutiva, foi realizada uma pesquisa survey explanatória de corte transversal com abordagem quantitativa, incluindo 113 projetos de E&P de cinco unidades executoras da empresa. Foi formulado um modelo conceitual, com cinco construtos e sete hipóteses de pesquisa, representativo dos efeitos de fatores externos sobre o desempenho ambiental dos projetos de E&P. Os dados foram tratados aplicando a Análise Fatorial Exploratória e a Modelagem de Equações Estruturais com aplicação dos softwares IBM® SPSS® Statistics 20.0 e IBM® SPSS® Amos 18.0. O modelo de equações estruturais foi reespecificado e estimado utilizando o método de Máxima Verossimilhança e o procedimento bootstrap com 2000 reamostragens, até alcançar adequados valores dos índices de ajustamento. O modelo mostrou boa aderência às evidências empíricas, representando uma teoria explicativa dos fatores que influenciam o desempenho ambiental dos projetos de E&P na empresa estudada. As estatísticas descritivas apontaram adequado desempenho dos projetos de E&P com relação aos efluentes descartados, volume de água reutilizada, redução de resíduos e práticas de reciclagem. Identificou-se que projetos de maior porte alcançam melhor desempenho ambiental em relação aos de menor tamanho. Não foram achadas diferenças significativas entre os desempenhos de projetos executados por unidades operacionais distintas. Os resultados da modelagem indicaram que nem a legislação ambiental, nem os agentes externos exercem influência significativa sobre a sistemática da gestão dos projetos de E&P. Os agentes externos atuam sobre a gestão ambiental da empresa exercitando capacidades colaborativas, obstrutivas e propositivas. A legislação ambiental é percebida como entrave ao desenvolvimento dos projetos ao longo de seu ciclo de vida, principalmente, pelas deficiências dos órgãos ambientais. Identificou-se que o sistema de gestão ambiental influencia diretamente o Programa de Desenvolvimento e Execução de Projetos de E&P, que, por sua vez, provoca efeitos diretos e indiretos sobre o desempenho ambiental. Finalmente, comprovou-se que o Sistema de Gestão Ambiental da empresa é determinante para o desempenho ambiental dos projetos de E&P, tanto pelos seus efeitos diretos, como pelos indiretos, estes últimos mediados pela sistemática de gestão dos projetos de E&P

Modelagem matemática para o transporte de partículas sujeitas a múltiplos mecanismos de retenção

Discrepancies between classical model predictions and experimental data for deep bed filtration have been reported by various authors. In order to understand these discrepancies, an analytic continuum model for deep bed filtration is proposed. In this model, a filter coefficient is attributed to each distinct retention mechanism (straining, diffusion, gravity interception, etc.). It was shown that these coefficients generally cannot be merged into an effective filter coefficient, as considered in the classical model. Furthermore, the derived analytic solutions for the proposed model were applied for fitting experimental data, and a very good agreement between experimental data and proposed model predictions were obtained. Comparison of the obtained results with empirical correlations allowed identifying the dominant retention mechanisms. In addition, it was shown that the larger the ratio of particle to pore sizes, the more intensive the straining mechanism and the larger the discrepancies between experimental data and classical model predictions. The classical model and proposed model were compared via statistical analysis. The obtained p values allow concluding that the proposed model should be preferred especially when straining plays an important role. In addition, deep bed filtration with finite retention capacity was studied. This work also involves the study of filtration of particles through porous media with a finite capacity of filtration. It was observed, in this case, that is necessary to consider changes in the boundary conditions through time evolution. It was obtained a solution for such a model using different functions of filtration coefficients. Besides that, it was shown how to build a solution for any filtration coefficient. It was seen that, even considering the same filtration coefficient, the classic model and the one here propposed, show different predictions for the concentration of particles retained in the porous media and for the suspended particles at the exit of the media

Avaliação das equações de predição das exigências de proteína para matrizes pesadas na fase de crescimento

Objetivou-se avaliar as equações para predição das exigências protéicas de matrizes pesadas em crescimento, desenvolvidas em estudos anteriores na Faculdade de Ciências Agrárias e Veterinárias da UNESP, Jaboticabal, tendo como padrão de comparação as recomendações estabelecidas pelo manual de criação da linhagem. Foram utilizadas 288 matrizes pesadas da linhagem Hubbard Hy-Yield, com cinco semanas de idade, em delineamento inteiramente casualizado, com dois tratamentos e seis repetições, sendo 24 aves por unidade experimental até 14 semanas e, de 15 a 20 semanas de idade, 20 aves por unidade experimental. Um tratamento foi representado pelas recomendações protéicas para a linhagem (testemunha) e o outro pelas seguintes equações de predição das exigências de proteína: de 5 a 8 semanas PB (g/ave/dia) = 2,02.P0,75+ 0,350.G, 9 a 14 semanas - PB = 2,02.P0,75 + 0,406.G e de 15 a 20 semanas PB = 2,02.P0,75 + 0,463.G, em que P é o peso corporal (kg) e G, o ganho de peso diário (g). Os resultados indicaram que as equações de predição para proteína, apesar de não comprometerem o peso ou a uniformidade corporal em relação ao tratamento testemunha, determinaram menor porcentagem de peito e maior de gordura abdominal e maior teor de gordura na carcaça, sendo atribuído à menor ingestão de proteína proporcionada pelas equações de predição e, conseqüentemente, deficiência em relação a alguns aminoácidos que não foram suplementados.

Método de avaliação do gerenciamento de resíduos em canteiro de obras da construção civil : proposta baseada em empresas construtoras da cidade de Natal-RN

According to great concern between the developed industrial activities and resultant impacts over the environment, an association of several factors have occurred, procedures to the efficient management of the rotation between economical development and the environment have been improved. A research in field have been realized inside building sites of companies in order to provide knowledge about the implemented and accomplish actions according to the resolution from CONAMA nº307. Trough the interview among the representations of the companies and photographic survey in loco, such as, what makes the companies implement this management, reutilization and recycling, transport and disposition. The present study had as objective: analyze the insertion of the used tools to residuals management, proposing improvements, in a way that it can be easily identified during the procedures execution in the building sites of the building companies of the city of Natal/RN. To reach the goal, in the first place a revision of the pertinent literature was performed; there for, it can be seen the relation between residues management and environment sustainability, once it happens in a continued way it may prevent the waste and reduces the risk that the activities way bring to the employees, community and environment; once found the great difficult faced with regard to labors, material, equipment, project, planning, costumer s interference, furnisher. And still, it could be verified wich materials generate greater indexes of residues in the works and the main occurrences of waste and loss. However a greater transparency is needed coming from the high administration in the commitment with the continued actions, to make it so, there must be a cultural change inside the company. There for there will be a greater productivity and quality of the under taking such as costumer s satisfaction

Um algoritmo para estimação de estado em alimentadores de distribuição de energia elétrica com base no método da soma de potências

Most algorithms for state estimation based on the classical model are just adequate for use in transmission networks. Few algorithms were developed specifically for distribution systems, probably because of the little amount of data available in real time. Most overhead feeders possess just current and voltage measurements at the middle voltage bus-bar at the substation. In this way, classical algorithms are of difficult implementation, even considering off-line acquired data as pseudo-measurements. However, the necessity of automating the operation of distribution networks, mainly in regard to the selectivity of protection systems, as well to implement possibilities of load transfer maneuvers, is changing the network planning policy. In this way, some equipments incorporating telemetry and command modules have been installed in order to improve operational features, and so increasing the amount of measurement data available in real-time in the System Operation Center (SOC). This encourages the development of a state estimator model, involving real-time information and pseudo-measurements of loads, that are built from typical power factors and utilization factors (demand factors) of distribution transformers. This work reports about the development of a new state estimation method, specific for radial distribution systems. The main algorithm of the method is based on the power summation load flow. The estimation is carried out piecewise, section by section of the feeder, going from the substation to the terminal nodes. For each section, a measurement model is built, resulting in a nonlinear overdetermined equations set, whose solution is achieved by the Gaussian normal equation. The estimated variables of a section are used as pseudo-measurements for the next section. In general, a measurement set for a generic section consists of pseudo-measurements of power flows and nodal voltages obtained from the previous section or measurements in real-time, if they exist -, besides pseudomeasurements of injected powers for the power summations, whose functions are the load flow equations, assuming that the network can be represented by its single-phase equivalent. The great advantage of the algorithm is its simplicity and low computational effort. Moreover, the algorithm is very efficient, in regard to the accuracy of the estimated values. Besides the power summation state estimator, this work shows how other algorithms could be adapted to provide state estimation of middle voltage substations and networks, namely Schweppes method and an algorithm based on current proportionality, that is usually adopted for network planning tasks. Both estimators were implemented not only as alternatives for the proposed method, but also looking for getting results that give support for its validation. Once in most cases no power measurement is performed at beginning of the feeder and this is required for implementing the power summation estimations method, a new algorithm for estimating the network variables at the middle voltage bus-bar was also developed

Dispositivos planares integrados utilizando método dinâmico com metamateriais e PBG

This work presents a theoretical, numerical and computation analysis of parameters of a rectangular microstrip antenna with metamaterial substrate, fin line as a coupler and also integrated devices like integrated filter antenna. It is applied theory to full-wave of Transverse Transmission Line - TTL method, to characterize the magnitude of the substrate and obtain the general equations of the electromagnetic fields. About the metamaterial, they are characterized by permittivity and permeability tensor, reaching to the general equations for the electromagnetic fields of the antenna. It is presented a study about main representation of PBG(Photonic Band Gap) material and its applied for a specific configuration. A few parameters are simulated some structures in order to reduce the physical dimensions and increase the bandwidth. The results are presented through graphs. The theoretical and computational analysis of this work have shown accurate and relatively concise. Conclusions are drawn and suggestions for future work

Análise de antenas de microfita com patches circulares sobre substratos anisotrópicos usando o método dos potenciais de Hertz

This work consists on the theoretical and numerical analysis of some properties of circular microstrip patch antennas on isotropic and uniaxial anisotropic substrates. For this purpose, a full wave analysis is performed, using Hertz Vector Potentials method in the Hankel Transform domain. In the numerical analysis, the moment method is also used in order to determine some characteristics of the antenna, such as: resonant frequency and radiation pattern. The definition of Hertz potentials in the Hankel domain is used in association with Maxwell´s equations and the boundary conditions of the structures to obtain the Green´s functions, relating the components of the current density on the patch and the tangential electric field components. Then, the Galerkin method is used to generate a matrix equation whose nontrivial solution is the complex resonant frequency of the structure. In the analysis, a microstrip antenna with only one isotropic dielectric layer is initially considered. For this structure, the effect of using superconductor patches is also analyzed. An analysis of a circular microstrip antenna on an uniaxial anisotropic dielectric layer is performed, using the Hertz vector potentials oriented along the optical axis of the material, that is perpendicular to the microstrip ground plane. Afterwards, the circular microstrip antenna using two uniaxial anisotropic dielectric layers is investigated, considering the particular case in which the inferior layer is filled by air. In this study, numerical results for resonant frequency and radiation pattern for circular microstrip antennas on isotropic and uniaxial anisotropic substrates are presented and compared with measured and calculated results found in the literature

Aplicação do método LTT às estruturas retangulares e triangulares em multicamadas e empilhadas em substratos PBG para comunicações móveis

Recently the planar antennas have been studied due to their characteristics as well as the advantages that they offers when compared with another types of antennas. In the mobile communications area, the need for this kind of antennas have became each time bigger due to the intense increase of the mobile communications that needs of antennas which operate in multifrequency and wide bandwidth. The microstrip antennas presents narrow bandwidth due the loss in the dielectric generated by radiation. Another limitation is the radiation pattern degradation due the generation of surface waves in the substrate. In this work some used techniques to minimize the disadvantages (previously mentioned) of the use of microstrip antennas are presented, those are: substrates with PBG material - Photonic Bandgap, multilayer antennas and with stacked patches. The developed analysis in this work used the TTL - Transverse Transmission Line method in the domain of Fourier transform, that uses a component of propagation in the y direction (transverse to the direction real of propagation z), treating the general equations of electric and magnetic field as functions of y and y . This work has as objective the application of the TTL method to microstrip structures with single and multilayers of rectangular and triangular patches, to obtaining the resonance frequency and radiation pattern of each structure. This method is applied for the treatment of the fields in stacked structures. The Homogenization theory will be applied to obtaining the effective permittivity for s and p polarizations of the substrate composed of PBG material. Numerical results for the triangular and rectangular antennas with single layer, multilayers resonators with triangular and rectangular patches are presented (in photonic and isotropic substrates). Conclusions and suggestions for continuity of this work are presented

Método dinâmico aplicado para antenas cilíndricas

Nowadays there has been a major breakthrough in the aerospace area, with regard to rocket launches to research, experiments, telemetry system, remote sensing, radar system (tracking and monitoring), satellite communications system and insertion of satellites in orbit. This work aims at the application of a circular cylindrical microstrip antenna, ring type, and other cylindrical rectangular in structure of a rocket or missile to obtain telemetry data, operating in the range of 2 to 4 GHz, in S-band. Throughout this was developed just the theoretical analysis of the Transverse transmission line method which is a method of rigorous analysis in spectral domain, for use in rockets and missiles. This analyzes the spread in the direction "ρ" , transverse to dielectric interfaces "z" and "φ", for cylindrical coordinates, thus taking the general equations of electromagnetic fields in function of e [1]. It is worth mentioning that in order to obtain results, simulations and analysis of the structure under study was used HFSS program (High Frequency Structural Simulator) that uses the finite element method. With the theory developed computational resources were used to obtain the numerical calculations, using Fortran Power Station, Scilab and Wolfram Mathematica ®. The prototype was built using, as a substrate, the ULTRALAM ® 3850, of Rogers Corporation, and an aluminum plate as a cylindrical structure used to support. The agreement between the measured and simulated results validate the established processes. Conclusions and suggestions are presented for continuing this work

Modelagem numérica de escoamento em padrão anular com transferência de calor e massa

Annular flow is the prevailing pattern in transport and energy conversion systems and therefore, one of the most important patterns in multiphase flow in ducts. The correct prediction of the pressure gradient and heat transfer coefficient is essential for optimizing the system s capacity. The objective of this work is to develop and implement a numerical algorithm capable of predicting hydrodynamic and thermal characteristics for upflow, vertical, annular flow. The numerical algorithm is then complemented with the physical modeling of phenomena that occurs in this flow pattern. These are, turbulence, entrainment and deposition and phase change. For the development of the numerical model, axial diffusion of heat and momentum is neglected. In this way the time-averaged equations are solved in their parabolic form obtaining the velocity and temperature profiles for each axial step at a time, together with the global parameters, namely, pressure gradient, mean film thickness and heat transfer coefficient, as well as their variation in the axial direction. The model is validated for the following conditions: fully-developed laminar flow with no entrainment; fully developed laminar flow with heat transfer, fully-developed turbulent flow with entrained drops, developing turbulent annular flow with entrained drops, and turbulent flow with heat transfer and phase change

Modelagem e simulação de um separador bifásico com inversão de fases induzida

The present work has the main goal to study the modeling and simulation of a biphasic separator with induced phase inversion, the MDIF, with the utilization of the finite differences method for the resolution of the partial differencial equations which describe the transport of contaminant s mass fraction inside the equipment s settling chamber. With this aim, was developed the deterministic differential model AMADDA, wich was admensionalizated and then semidiscretizated with the method of lines. The integration of the resultant system of ordinary differential equations was realized by means of a modified algorithm of the Adam-Bashfort- Moulton method, and the sthocastic optimization routine of Basin-Hopping was used in the model s parameter estimation procedure . With the aim to establish a comparative referential for the results obtained with the model AMADDA, were used experimental data presented in previous works of the MDIF s research group. The experimental data and those obtained with the model was assessed regarding its normality by means of the Shapiro-Wilk s test, and validated against the experimental results with the Student s t test and the Kruskal-Wallis s test, depending on the result. The results showed satisfactory performance of the model AMADDA in the evaluation of the MDIF s separation efficiency, being possible to determinate that at 1% significance level the calculated results are equivalent to those determinated experimentally in the reference works

Avaliação e modelagem da absorção de H2S do gás natural em coluna de leito estagnado

Natural gas, although basically composed by light hydrocarbons, also presents contaminant gases in its composition, such as CO2 (carbon dioxide) and H2S (hydrogen sulfide). The H2S, which commonly occurs in oil and gas exploration and production activities, causes damages in oil and natural gas pipelines. Consequently, the removal of hydrogen sulfide gas will result in an important reduction in operating costs. Also, it is essential to consider the better quality of the oil to be processed in the refinery, thus resulting in benefits in economic, environmental and social areas. All this facts demonstrate the need for the development and improvement in hydrogen sulfide scavengers. Currently, the oil industry uses several processes for hydrogen sulfide removal from natural gas. However, these processes produce amine derivatives which can cause damage in distillation towers, can cause clogging of pipelines by formation of insoluble precipitates, and also produce residues with great environmental impact. Therefore, it is of great importance the obtaining of a stable system, in inorganic or organic reaction media, able to remove hydrogen sulfide without formation of by-products that can affect the quality and cost of natural gas processing, transport, and distribution steps. Seeking the study, evaluation and modeling of mass transfer and kinetics of hydrogen removal, in this study it was used an absorption column packed with Raschig rings, where the natural gas, with H2S as contaminant, passed through an aqueous solution of inorganic compounds as stagnant liquid, being this contaminant gas absorbed by the liquid phase. This absorption column was coupled with a H2S detection system, with interface with a computer. The data and the model equations were solved by the least squares method, modified by Levemberg-Marquardt. In this study, in addition to the water, it were used the following solutions: sodium hydroxide, potassium permanganate, ferric chloride, copper sulfate, zinc chloride, potassium chromate, and manganese sulfate, all at low concentrations (»10 ppm). These solutions were used looking for the evaluation of the interference between absorption physical and chemical parameters, or even to get a better mass transfer coefficient, as in mixing reactors and absorption columns operating in counterflow. In this context, the evaluation of H2S removal arises as a valuable procedure for the treatment of natural gas and destination of process by-products. The study of the obtained absorption curves makes possible to determine the mass transfer predominant stage in the involved processes, the mass transfer volumetric coefficients, and the equilibrium concentrations. It was also performed a kinetic study. The obtained results showed that the H2S removal kinetics is greater for NaOH. Considering that the study was performed at low concentrations of chemical reagents, it was possible to check the effect of secondary reactions in the other chemicals, especially in the case of KMnO4, which shows that your by-product, MnO2, acts in H2S absorption process. In addition, CuSO4 and FeCl3 also demonstrated to have good efficiency in H2S removal

Equações preditivas para as pressões respiratórias máximas de crianças brasileiras

I ntroduction: The assessment of respiratory muscle strength is important in the diagnosis and monitoring of the respiratory muscles weakness of respiratory and neuromuscular diseases. However, there are still no studies that provide predictive equations and reference values for maximal respiratory pressures for children in our population. Aim: The purpose of this study was to propose predictive equations for maximal respiratory pressures in healthy school children. Method: This is an observational cross-sectional study. 144 healthy children were assessed. They were students from public and private schools in the city of Natal /RN (63 boys and 81 girls), subdivided in age groups of 7-8 and 9-11 years. The students presented the BMI, for age and sex, between 5 and 85 percentile. Maximal respiratory pressures were measured with the digital manometer MVD300 (Globalmed ®). The maximal inspiratory pressure (MIP) and maximal expiratory pressures (MEP) were measured from residual volume and total lung capacity, respectively. The data were analyzed using the SPSS Statistics 15.0 software (Statistical Package for Social Science) by assigning the significance level of 5%. Descriptive analysis was expressed as mean and standard deviation. T'Student test was used for unpaired comparison of averages of the variables. The comparison of measurements obtained with the predicted values in previous studies was performed using the paired t'Student test. The Pearson correlation test was used to verify the correlation of MRP's with the independent variables (age, sex, weight and height). For the equations analysis the stepwise linear regression was used. Results: By analyzing the data, we observed that in the age range studied MIP was significantly higher in boys. The MEP did not differ between boys and girls aged 7 to 8 years, the reverse occurred in the age between 9 and 11 years. The boys had a significant increase in respiratory muscle strength with advancing age. Regardless sex and age, MEP was always higher than the MIP. The reference values found in this study are similar to a sample of Spanish and Canadian children. The two models proposed in previous studies with children from other countries were not able to consistently predict the values observed in this studied population. The variables sex, age and weight correlated with MIP, whereas the MEP was also correlated with height. However, in the regression models proposed in this study, only gender and age were kept exerting influence on the variability of maximal inspiratory and expiratory pressures. Conclusion: This study provides reference values, lower limits of normality and proposes two models that allow predicting, through the independent variables, sex and age, the value of maximal static respiratory pressures in healthy children aged between 7 and 11 years old