A maioria dos líquidos encontrados na natureza são não newtonianos e o estudo do seu comportamento reológico tem uma importância significante em diferentes áreas da engenharia. Entre eles, existe uma classe de fluidos que exibem pequena deformação aparente quando sujeitos a um nível de tensões inferior a uma tensão limite de escoamento, referido como comportamento viscoplástico. Nesta classe de materiais, alguns apresentam também comportamento elástico quando submetidos a baixas taxas de cisalhamento. A presente Tese tem como objetivo o estudo numérico de escoamentos bidimensionais em regime permanente de fluidos elasto-viscoplásticos através de uma expansão-contração planar. O modelo mecânico é definido pelas equações de conservação de massa e de balanço de momentum acopladas ao modelo elasto-viscoplástico proposto nesta Tese. Esta modelagem é aproximada por um método de elementos finitos multi-campos estabilizado baseado na metodologia de Galerkin mínimos-quadrados que possui como variáveis primais os campos de tensão extra polimérica, velocidade e pressão. As condições de compatibilidade entre os sub-espaços de elementos finitos para tensão extra-velocidade e velocidade-pressão são violadas, permitindo assim a utilização de interpolações de igual ordem. O método estabilizado foi implementado no código de elementos finitos para fluidos não newtonianos em desenvolvimento no Laboratório de Mecânica dos Fluidos Aplicada e Computacional (LAMAC) da Universidade Federal do Rio Grande do Sul. Nesta Tese é adotada uma metodologia alternativa para a definição das zonas rígidas do escoamento como sendo a posição onde a taxa de cisalhamento é igual a um valor dado pela relação de parâmetros reológicos do fluido, especificamente a tensão limite de escoamento e a viscosidade newtoniana para baixas taxas de cisalhamento. Nas simulações numéricas realizadas, o tempo de relaxação adimensional, o número de salto, o coeficiente power-law, a vazão adimensional e a massa específica adimensional são variados de forma a avaliar de que modo influenciam na dinâmica de escoamentos elastoviscoplásticos. Os resultados obtidos estão qualitativamente de acordo com a literatura, atestando a estabilidade da formulação empregada. / Non-Newtonian fluids are the majority of liquids found in nature and the study of their rheological behavior has a significant importance on different areas of engineering. Among them, there is a class of materials that exhibits little apparent deformation when subjected to a stress level behind an yield stress, referenced as viscoplastic material. In this class of materials, some fluids also exhibit elastic behavior at low shear rates. The present work aimed to a numerical study of two-dimensional steady state laminar flows of elasto-viscoplastic fluids through a planar expansion-contraction cavity. The mechanical model was defined by the mass conservation and momentum balance equations coupled to the elasto-viscoplastic model porposed in this work. This modeling has been approximated by a stabilized multi-field finite element method based on the Galerkin least-squares methodology, having as primal variables the elastic extra-stress component, velocity and pressure fields. In this way, the compatibility conditions between the extra-stress-velocity and pressure-velocity (Babuška- Brezzi condition) finite element subspaces are violated, allowing to use equal-order finite element interpolations. The stabilized method has been implemented in the finite element code for non-Newtonian fluids under development at the Laboratory of Applied and Computational Fluid Mechanics (LAMAC) of the Federal University of Rio Grande do Sul. An alternative methodology is adopted to define the yield surface as the position where the strain rate is equal to a value given by the relation of the rheological parameters of the fluid, namely the yield stress and the viscosity at low shear rates. In the performed numerical simulations, the non-dimensional relaxation time, the jump number, the power-law coefficient, the non-dimensional flow rate, and the non-dimensional density are varied in order to evaluate their influence on the elasto-viscoplastic fluid dynamics. All results found are in qualitatively accordance with the affine literature, and attesting the good stability features of the formulation.

Em problemas de viscoelasticidade computacional, a discretização espacial para a solução global das equações de equilíbrio é acoplada à discretização temporal para a solução de um problema de valor inicial local do fluxo viscoelástico. É demonstrado que este acoplamento espacial-temporal (ou global-local) éconsistente se o tensor de deformação total, agindo como elemento acoplador, tem uma aproximação de ordem p ao longo do tempo igual à ordem de convergência do método de integração de Runge-Kutta (RK). Para a interpolação da deformação foram utilizados polinômios baseados em soluções obtidas nos tempos tn+1, tn, . . ., tn+2−p, p ≥ 2, fornecendo dados consistentes de deformação nos estágios do RK. Em uma situação onde tal regra para a interpolação da deformação não é satisfeita, a integração no tempo apresentará, consequentemente, redução de ordem, baixa precisão e, por conseguinte, eficiência inferior. Em termos gerais, o propósito é generalizar esta condição de consistência proposta pela literatura, formalizando-a matematicamente e o demonstrando através da utilização de métodos de Runge-Kutta diagonalmente implícitos (DIRK) até ordem p = 4, aplicados a modelos viscoelásticos não-lineares sujeitos a deformações finitas. Através de exemplos numéricos, os algoritmos de integração temporal adaptados apresentaram ordem de convergência nominal e, portanto, comprovam a validade da formalização do conceito de interpolação consistente da deformação. Comparado com o método de integração de Euler implícito, é demonstrado que os métodos DIRK aqui aplicados apresentam um ganho considerável em eficiência, comprovado através dos fatores de aceleração atingidos. / In computational viscoelasticity, spatial discretization for the solution of the weak form of the balance of linear momentum is coupled to the temporal discretization for solving a local initial value problem (IVP) of the viscoelastic flow. It is shown that this spatial- temporal (or global-local) coupling is consistent if the total strain tensor, acting as the coupling agent, exhibits the same approximation of order p in time as the convergence order of the Runge-Kutta (RK) integration algorithm. To this end we construct interpolation polynomials based on data at tn+1, tn, . . ., tn+2−p, p ≥ 2, which provide consistent strain data at the RK stages. If this novel rule for strain interpolation is not satisfied, time integration shows order reduction, poor accuracy and therefore less efficiency. Generally, the objective is to propose a generalization of this consistency idea proposed in the literature, formalizing it mathematically and testing it using diagonally implicit Runge-Kutta methods (DIRK) up to order p = 4 applied to a nonlinear viscoelasticity model subjected to finite strain. In a set of numerical examples, the adapted time integrators obtain full convergence order and thus approve the novel concept of consistency. Substantially high speed-up factors confirm the improvement in the efficiency compared with Backward Euler algorithm.

O comportamento estrutural de veículos de transporte coletivo, quando submetido a eventos de impacto, é de grande relevância na engenharia automobilística, pois principalmente o ônibus, vem se tornando um dos meios de transporte de elevada importância, sendo que seu uso vem aumentando a cada ano. No entanto, as estruturas de ônibus que circulam nas estradas brasileiras não são capazes de resistir, sem causar danos aos passageiros, a eventos de impacto em acidentes, o que pode ser verificado em pesquisas divulgadas pela Agência Nacional de Transportes Terrestres (ANTT) em 2007, que mostra a evolução de acidentes e vítimas envolvendo ônibus, que cada ano aumentam consideravelmente. Neste contexto, pretende-se desenvolver um absorvedor de impacto que aumente a capacidade de absorção da energia de impacto e evite a invasão da região de sobrevivência dos passageiros, e desta forma diminuindo o número de vítimas em acidentes. Serão abordados neste trabalho eventos de impacto semi-frontal, justificado pela gravidade dos acidentes desse tipo nas estradas, sendo este o tipo de acidente que deixa um maior número de vítimas, pois geralmente parte da estrutura lateral do ônibus é removida, expondo os passageiros no período de impacto. Esse efeito é conhecido como “Efeito Abridor de Latas”. O modelo numérico deste trabalho será formado por elementos unifilares, sendo a estrutura do ônibus formada por vigas flexíveis e/ ou rígidas, unidas através de juntas não lineares rotacionais e translacionais. A rigidez de cada junta não linear é obtida a partir de métodos analíticos que descrevem o comportamento do tubo de parede fina quando submetido a impacto. / The structural behavior of transportation vehicles when subjected to events of impact is of great relevance in automotive engineering, mainly because the bus has become a highly important alternative means of transportation, and its use is increasing every year. However, the structures of the buses circulating on Brazilian roads are not able to withstand strong impact, without causing passengers’ injury in events where accidents occur, which can be shown in surveys released by the National Transport Agency (ANTT) in 2007, presenting that the great incidence of accidents with victims involving buses is increasing considerably every year. In this context, the development of an impact absorber that had the capability to increase the impact energy absorption and prevent the invasion of part of bus on the region where passengers are located (reducing the number of fatalities in accidents) are the objectives of this work. In this work, a semi-frontal impact will be presented and studied due to the severity of such accident, as the side of the bus structure is removed, exposing the passengers during the impact. This is the effect known as “Can opener effect”. The numerical model of this study will consist of a single-wire element, with a bus structure formed by flexible and/or rigid beams joined by nonlinear rotational and translational joints. The rigidity of each nonlinear joint is obtained by analytical methods describing the behavior of a thin-walled tube subjected to an impact force.

Fluidos Não-newtonianos estão presentes em grande quantidade na natureza e em muitas aplicações industriais. Tais fluidos apresentam uma variedade de comportamentos específicos. Dentre eles podemos citar: elasticidade, plasticidade, viscosidade variando com a taxa de deformação, tixotropia. Esta tese de doutorado apresenta um recente modelo para Fluidos Estruturados Tixotrópicos (de Souza Mendes, 2009), no qual baseasse nas tradicionais equações de balanço de momentum e massa, em uma equação constitutiva de Maxwell-b modificada e uma equação de evolução para o nível de estrutura do material. Na equação constitutiva de Maxwell-b modificada, tempo de relaxação e viscosidade estrutural são funções da estrutura do material. A equação evolutiva para o parâmetro de estruturação prevê efeitos de tixotropia, dependentes do tempo de equilíbrio do material. As equações que constituem o modelo mecânico para Fluidos Estruturados Tixotrópicos são aproximadas numericamente através do Método de Elementos Finitos de Galerkin Mínimos-quadrados (GLS). Tal método adicionou termos dependentes de malha ao método clássico de Galerkin, de forma a aumentar a estabilidade da formulação original, proporcionando compatibilizar os sub-espaços de velocidade e pressão. São investigados escoamentos creeping flow de Fluidos Estruturados Tixotrópicos em dois problemas clássicos de engenharia: em torno de um cilindro confinado entre duas placas e em uma expansão planar abrupta. Parâmetros relevantes ao modelo são variados em faixas pertinentes ao problema físico simulado e os resultados demonstraram uma boa capacidade de previsão. / Non-Newtonian fluids are present in large quantities in nature and in many industrial applications. Such fluids exhibit a variety of specific behaviors. Among them we can mention: elasticity, plasticity, viscosity varying with the strain rate, thixotropy, etc. This thesis presents a new model for thixotropic Structured Fluids (de Souza Mendes, 2009), in which were a basis in traditional balance equations of mass and momentum, the constitutive equation of Maxwell-type and an evolution equation for the level structure of the material. In the modified constitutive equation of Maxwell-b, relaxation time and structural viscosity are functions of material structure. The structure parameter evolutive equation provides effects of thixotropy, dependent on the equilibration time of the material. The equations which constitute the mechanical model for Structured Thixotropic fluids are approximated numerically using the Galerkin Least-squares Finite Element Method (GLS). Such a method added mesh dependent terms to the classical method of Galerkin in order to increase the stability of the original formulation, providing compatibility between sub-spaces of velocity and pressure. Are investigated Structured Thixotropic Fluids creeping flow in two classical problems in engineering: around a cylinder confined between two plates and in a planar abrupt expansion. Relevant parameters to the model are varied in pertinent ranges to physical problem simulated and the results demonstrated a good predictive capability.

O presente trabalho utiliza o método Constructal Design para desenvolver o estudo numérico da configuração de materiais de alta condutividade térmica em forma de “Y” que minimiza a resistência ao fluxo de calor, quando áreas ocupadas pelos materiais de alta e baixa condutividades são mantidas constantes. Para a solução numérica da equação diferencial da difusão do calor e suas respectivas condições de contorno, foi utilizado o software MATLAB ®, mais especificamente a ferramenta PDETOOL, Partial Differential Equations Tool. O objetivo deste trabalho é a minimização da resistência térmica do sistema gerador de calor com baixa condutividade térmica com a utilização de vias em formato de Y com material de alta condutividade térmica e volume constante, sendo variáveis os comprimentos e espessuras do material dos ramos simples e bifurcados. Todas as possibilidades geométricas foram avaliadas e a geometria ótima foi aquela que conduziu a menor resistência térmica. Duas condições são apresentadas, a primeira tem os ramos e a base da geometria “Y” com igual condutividade térmica. Os resultados para esta configuração mostram que existem valores específicos para os graus de liberdade que minimizam a resistência térmica. Nesse caso, os ramos se degeneraram e a configuração ótima tem a forma de um “V”. A segunda configuração apresenta combinações de condutividade térmica diferentes, para os ramos e a bases. Para estes casos obteve-se um valor otimizado próximo de 1 para a razão entre os comprimentos dos ramos simples e bifurcados, indicando que a configuração otimizada tem realmente a forma de um “Y” o que demonstra a dependência entre a geometria e as condições impostas pelo meio. Embora o design inicial do Y possa assumir diversas configurações, quando comparado o primeiro design com o design final, no caso do Y com iguais condutividades térmicas se conseguiu uma melhora superior a 28% e no caso do Y com condutividades diferentes mais de 30 %. Finalmente, este trabalho mostra que a geometria otimizada é aquela que melhor distribui as imperfeições, isto é, os pontos quentes (pontos de temperatura máxima), o que está de acordo com o princípio da ótima distribuição das imperfeições. / The present work used the method Constructal Design to develop numerical analyses of pathways of high thermal conductivity in "Y" shape which minimizes the thermal resistance when areas occupied by the materials of high and low conductivities are kept constant. For the numerical solution of the differential equations of heat diffusion and their boundary conditions, we used the MATLAB ® software, specifically the PDETOOL tool. The aim was to minimize the thermal resistance of the heat generator system with low thermal conductivity with the use of Y-shaped pathways with high thermal conductivity and constant volume, with variable lengths and thicknesses of material from stem and forked branches. All geometric possibilities were evaluated and the optimal geometry was that which resulted in lower thermal resistance. Two conditions were studied. In the first one the stem and branches of the "Y" have equal thermal conductivity. The results for this configuration show that there are specific values for the degrees of freedom to minimize the thermal resistance. In this case, the branches have degenerated and the optimum configuration has the shape of a "V". The second configuration offers different combinations of thermal conductivity, for branches and bases. For these cases we obtained a optimized value close to 1 for the ratio between the lengths of stem and bifurcated branches, indicating that the optimized configuration actually has the shape of a "Y" which shows the dependency of geometry and condition imposed by the environment. Although the initial design of Y can take various configurations, when compared the first design to the final design, in the case of Y with equal thermal conductivity it this improvement was achieved an improvement greater than 28% and in the case of Y with different conductivities over 30%. Finally, this study showed that the optimized geometry is the one that better distributes imperfections, this is, hot spots (points of maximum temperature), which is in accordance with the principle of the optimal distribution of imperfections.

As imagens são sinais que possuem muita informação. Os objetos representados em ima- gens podem sofrer deformações fazendo com que suas características mudem, o que difi- culta o reconhecimento do objeto. A chave para o sucesso de um sistema de reconheci- mento de padrões em imagens é escolher adequadamente a sua abordagem e um modelo para as feições presentes nas imagens. Uma das dificuldades é extrair e selecionar as feições que são mais discriminantes entre as diferentes classes usadas para modelar um objeto. Os modelos ativos de formas (Active Shape Models ASM) adaptam-se às defor- mações de um objeto. O objeto a ser modelado, pode ser representado com um modelo de pontos distribuídos (Point Distribution Models PDM). O PDM consiste de pontos de interesse ou marcadores, que permitem a extração de feições em localizações específicas do objeto. Após tratar a incerteza da localização e oclusão dos marcadores é possível ex- trair as feições mais representativas, obtendo-se um desempenho alto em termos da taxa de reconhecimento. Nesta tese são introduzidas novas formas para extrair e selecionar feições com modelos ativos de formas, que melhoram a taxa de classificação onde há objetos deformáveis. Esta tese é inovadora no sentido de aperfeiçoar o uso de ASMs na classificação de faces humanas, e na sua aplicação no monitoramento visual de outros tipos de objetos deformáveis. / Images are signals that have a lot of information. The objects depicted in images may suf- fer deformations causing changes in their characteristics, which hinders the recognition of the object. The key to the success of a system of pattern recognition in images is to choose properly your approach and a model for the features in the images. One difficulty is to extract and select the features that are most discriminating between different classes used to model an object. The Active Shape Models (ASM) adapt to deformations of an object. The object to be modeled can be represented with a Points Distribution Model (PDM). The PDM consists of points of interest or landmarks that allow the extraction of features in specific locations of the object. After treating the uncertainty of the location and occlu- sion of the landmarks it is possible to extract the most representative features, obtaining a high performance in terms of recognition rate. This thesis introduces new ways to extract and select features with ASMs, which improve the classification rate where deformable objects are present. This thesis is innovative in the sense that improves the use of ASMs in the classification of human faces, and that can be applied in visual monitoring of other types of deformable objects.

Na presente tese é resolvida a equação de difusão de nêutrons estacionária, bem como problemas de cinética, em geometria unidimensional cartesiana multi-região considerando o modelo de multigrupos de energia. Um dos objetivos e inovação neste trabalho é a obtenção de uma solução aproximada com estimativa de erro, controle de precisão e na forma de uma expressão analítica. Com esse tipo de solução não há a necessidade de recorrer a esquemas de interpolação, geralmente necessários em caso de discretizações do domínio. O fluxo de nêutrons é expandido em uma série de Taylor cujos coeficientes são encontrados utilizando a equação diferencial e as condições de contorno e interface. O domínio é dividido em várias células, cujo tamanho e o grau do polinômio são ajustáveis de acordo com a precisão requerida. Para resolver o problema de autovalor é utilizado o método da potência. A metodologia é aplicada em um benchmark que consiste na solução da equação de difusão como condição inicial e na solução de problemas de cinética para diferentes transientes. Os resultados são comparados com sucesso com resultados da literatura. A convergência da série é garantida pela aplicação de um raciocínio baseado no critério de Lipschitz para funções contínuas. Cabe ressaltar que a solução obtida, em conjunto com a análise da convergência, mostra a solidez e a precisão dessa metodologia. / In the present dissertation the one-dimensional neutron diffusion equation for stationary and kinetic problems in a multi-layer slab has been solved considering the multi-group energy model. One of the objectives and innovation in this work is to obtain an approximate solution with error estimation, accuracy control and in the form of an analytical expression. With this solution there is no need for interpolation schemes, which are usually needed in case of discretization of the domain. The neutron flux is expanded in a Taylor series whose coefficients are found using the differential equation and the boundary and interface conditions. The domain is divided into several layers, whose size and the polynomial order can be adjusted according to the required accuracy. To solve the eigenvalue problem the conventional power method has been used. The methodology is applied in a benchmark problem consisting of the solution of the diffusion equation as an initial condition and solving kinetic problems for different transients. The results are compared successfully with the ones in the literature. The convergence of the series is guaranteed by applying a criterion based on the Lipschitz criterion for continuous functions. Note that the solution obtained, together with the convergence analysis, shows the robustness and accuracy of this methodology.

Neste trabalho, são inicialmente discutidas as dificuldades referentes à obtenção de resultados numéricos para a operação de uma turbina Savonius independentes do grau de discretização, tamanho de domínio de cálculo e de máximo tempo físico simulado. Também são relatadas as divergências entre as metodologias numéricas e experimentais adotadas por diversos autores, que dificultam análises e comparações dos resultados obtidos por meio dessas metodologias com os resultados próprios obtidos. Devido a esses fatos, no presente trabalho, uma série de procedimentos experimentais e numéricos são realizados para efetuar análises do escoamento sobre uma turbina eólica Savonius. Nos experimentos em canal aerodinâmico, perfis de velocidade e parâmetros da turbulência são obtidos pela técnica de anemometria de fio quente. Medições com o uso de tubos de Pitot e manômetros eletrônicos são efetuadas para avaliar a variação da pressão e os perfis de velocidade média em posições selecionadas. Além de dados para análise, informações úteis para uso como condições de contorno nas simulações numéricas também são obtidas. Os fenômenos são reproduzidos através de simulações numéricas pelo Método de Volumes Finitos, que solucionam as equações da continuidade, de Navier-Stokes com médias de Reynolds e do modelo de turbulência k-ω SST. Análises experimentais e numéricas considerando o escoamento sobre um cilindro, que mantém semelhanças com o escoamento sobre o rotor, também são realizadas. Simulações numéricas do escoamento sobre o cilindro são efetuadas, fornecendo resultados representativos do escoamento real, quando geometrias tridimensionais são aplicadas na modelagem numérica. Nas simulações do escoamento sobre o rotor Savonius em condição estática, resultados representativos do escoamento real são obtidos com o uso de uma modelagem que leva em consideração a rugosidade das pás do rotor, estacionado na posição angular de 90°. Para posições angulares menores, não se obteve uma boa concordância entre os resultados experimentais e numéricos. A realização deste trabalho fornece informações úteis para a análise do fenômeno e tem potencial para contribuir com futuros trabalhos desse tema. / This research work initially presents a discussion about the difficulties related to obtaining numerical results for the operation of a turbine Savonius independent of the degree of discretization, calculation domain size and maximum physical time of the simulation. The differences between the numerical and experimental methodologies adopted by various authors difficult the analysis and comparisons of the results obtained through these methods with the results obtained by the methodology. Due to these facts, in this research work, a series of experimental and numerical procedures are performed to conduct analyzes of flow over a Savonius wind turbine. In the experiments on aerodynamic channel, velocity and turbulence profiles parameters are obtained by the technique of hot wire anemometry. Measurements using Pitot tubes and electronic manometers are made to evaluate the variation of pressure and mean velocity profiles at selected positions. In addition to data analysis, useful information for use as boundary conditions in the numerical simulations are also obtained. The phenomena are reproduced through numerical simulations by Finite Volume Method, that solve the equations of continuity, Reynolds-averaged Navier–Stokes equations and the equation of the turbulence model k-ω SST. Experimental and numerical analyzes considering the flow over a cylinder, which holds similarities with the flow over the rotor, are also performed. Numerical simulations of the flow over the cylinder are made, providing results representative of the actual flow when three-dimensional geometries are applied in numerical modeling. In flow simulations over the Savonius rotor in static condition, at 90°, representative results of the actual flow are obtained using a model that takes into account the roughness of the rotor blades. For smaller angular positions a good agreement between experimental and numerical results was not obtained. This work provides useful information for the analysis of the phenomenon and has the potential to contribute to future work on this theme.

Sistemas de Distribuição de Energia Elétrica (SDEE) estão, pelas próprias condições de funcionamento, sujeitos a perdas de energia e quedas nos níveis de tensão. Uma das formas de reduzir as perdas e ao mesmo tempo melhorar o perfil de tensão do SDEE consiste na utilização de bancos de capacitores e reguladores de tensão ao longo dos alimentadores. Para isto a alocação ótima destes equipamentos torna-se de fundamental importância, logo este trabalho tem por objetivo a utilização de um método quadrático de programação por metas multiobjetivo para a ótima alocação desses equipamentos visando a redução de perdas (meta 1) com o menor custo de investimento (meta 2) sujeito a restrições do fluxo de potência, limites técnicos de operação, mudança no controle dos capacitores chaveados e dos taps dos reguladores de tensão. O constante aumento na penetração da geração distribuída traz novos desafios para o controle e a operação dos SDEE. A metodologia proposta considera os efeitos destas fontes de geração e os seus impactos na alocações ótima dos bancos de capacitores e reguladores de tensão. Com o intuito de resolver o problema de otimização em programas comerciais como Gurobi e Cplex uma metodologia linear de cálculo de fluxo de potência é utilizada. Os algoritmos foram desenvolvidos no software Matlab, e posteriormente foi realizada a conversão para linguagem GAMS (Generic Algebraic Modeling System) para utilização dos solvers na plataforma NEOS (Network Enabled Optimization System). / Power distribution systems (DS) are susceptible to energy loss and voltage drops, mainly because of the way that they are operated and built. One way to minimize system loss and improve the voltage profile of DS is with the use of capacitor banks and voltage regulators (VRs) along the distribution feeders. In this scenario, the optimal placement of these devices becomes crucial. In this paper, a quadratic multi objective goal programming method is proposed to obtain the optimal placement of these devices, where the goals are energy loss (goal 1) and the total investment cost (goal 2). The model is subject to system operation constraints, such as power flow demands, limits of voltages and currents, number of stages of switched capacitors and the tap control of VRs. The increase of distributed generation penetration brings new challenges for control and operation of DS. The proposed methodology considers the effects of these power sources and their impact on optimal allocations of capacitors and VRs. With the aim of solving the optimization problem in commercial programs, a linear power flow is used. The algorithms were developed with Matlab software and the conversion was made to a General Algebraic Modeling System (GAMS) language and solved in the NEOS (Network Enabled Optimization System) platform.

A velocidade de queima laminar adiabática é um importante parâmetro da combustão que dita o comportamento de chamas pré-misturadas. Dos métodos disponíveis para a medição desse parâmetro, o método do fluxo de calor destaca-se pela simplicidade e precisão. No presente trabalho, esse método é utilizado para medir a velocidade de queima de biogás (modelado como CH4 com diferentes níveis de diluição com CO2) e de gás de síntese (modelado como uma mistura de CH4, H2, CO, CO2 e N2) em ar a 298 K e 1 atm. Tais gases são de crescente interesse para a sociedade em função de aspectos ambientais, porém, suas velocidades de queima não foram amplamente estudadas ainda. Os resultados obtidos são comparados com as previsões de cinco mecanismos cinéticos (GRI-Mech 3.0, Davis et al., Konnov, San Diego e USC Mech II) a fim de avaliar a sua capacidade preditiva. Os resultados experimentais e numéricos das velocidades de queima de biogás e ar apresentam uma boa concordância e as incertezas encontradas foram condizentes com as relatadas na literatura. Os resultados experimentais desse gás foram parametrizados em uma correlação empírica de fácil utilização em modelos numéricos. As medições da velocidade de queima de gás de síntese e ar, por outro lado, apresentaram valores inferiores às previsões numéricas de todos os mecanismos estudados. Os dados experimentais da literatura, para a mesma mistura, diferem tanto em valores quanto em comportamento dos resultados do presente trabalho. Tal comportamento está provavelmente relacionado a alguma contaminação no CO utilizado, já que quando esse gás está presente observa-se uma chemi-luminescência não relatada na literatura. / The adiabatic laminar burning velocity is an important combustion parameter that dictates premixed flames characteristics. Among the measuring methods available in literature, the heat flux method stands out for its simplicity and accuracy. In the present work, this method is used to measure the adiabatic laminar burning velocity of biogas (modeled as CH4 with different dilution levels with CO2) and syngas (modeled as a CH4, H2, CO, CO2 and N2 mixture) in air at 298 K and 1 atm. Such gases are of growing society interest due to environmental aspects, however, their adiabatic laminar burning velocity have not been widely studied yet. The experimental results are compared to predictions of five kinetic mechanisms (GRI-Mech 3.0, Davis et al., Konnov, San Diego e USC Mech II) to evaluate their predictive capacity. Experimental and numerical results of biogas/air mixtures adiabatic laminar burning velocity show good agreement and the found uncertainties are in agreement with literature. Experimental results of this gas were fitted in an empiric correlation of simple numerical application. Experimental results of the laminar burning velocity of syngas/air, on the other hand, show lower values than the numerical predictions of all studied kinetic mechanisms. Literature available data for the same mixture differ both in values and behavior of the present work results. Such behavior is probably related to some contamination on the CO used since a chemi-luminescence not reported in literature can be noted when this gas is present.