Este trabalho, realizado com o auxílio da termodinâmica computacional, dá continuação a um estudo sistemático, que trata dos fundamentos dos tratamentos termoquímicos metalúrgicos relacionados à nitretação e à nitrocarburação. O estudo concentrou-se nas ligas da família ABNT 10XX, do aço ABNT 4140 e do aço ferramenta AISI M2. A composição das ligas utilizadas na simulação, contudo, não correspondem exatamente à comercial. Uma simplificação foi necessária para que a simulação do processo fosse facilitada. O uso da termodinâmica computacional como ferramenta para determinação das fases presentes (ou que podem vir a se formar) em um sistema a certa temperatura e pressão simulando o processo de nitretação mostrou-se uma valiosa ferramenta de análise.Apesar de, na prática, estes prognósticos ainda dependerem de uma cinética micro ou macroscópica favorável para ser verdadeiros, a comparação dos resultados obtidos com a literatura mostrou uma boa concordância nos resultados, validando, assim, o emprego desta técnica de análise. Assim, além deste estudo ser um avanço na compreensão das fases que participam das microestruturas (nitretos e carbonitretos) que se formam nas ligas estudadas, também pode ser utilizado como um balizador na definição de parâmetros que permitam a seleção destas fases finais desejadas, economizando tempo e recursos na realização de ensaios experimentais. / This work, accomplished with the help of computational thermodynamics, provides the development of a systematic study concerning the basis of the thermochemical treatments related to nitriding and nitrocarburising. This study focused on the alloys of the ABNT 10XX family, on the ABNT 4140 steel and on the AISI M2 tool steel. However, the compositions used in the simulation do not correspond exactly to the commercial composition of the alloy, for they had their composition simplified in order to facilitate the simulation process. The use of computational thermodynamics as a tool to determine the phase present (or that can be formed) in a system under a specific temperature and pressure, simulating the nitriding process, proved to be a valuable analysis tool. Nevertheless, in practice, these prognostics still depend on a micro or macro kinetics in order to be true; the comparison of the results obtained with those found in the literature showed good agreement of results, validating, thus, the use of this analysis technique. So, besides this study being an advance when it comes to understanding phases (nitrides and carbonitrides) that are formed in the studied alloys, it can also be used as a mark in defining parameters that allow the selection of final desired structures, saving time and resources.

O presente trabalho pesquisa o fenômeno da transferência simultânea de calor e umidade em solos insaturados que envolvem dutos enterrados. Estes dutos e o solo envolvente compõem o chamado sistema trocador-armazenador de calor no solo, que é muitas vezes utilizado como fonte complementar de aquecimento em estufas solares agrícolas. O funcionamento do sistema trocador-armazenador de calor é baseado na energia armazenada no solo durante os horários de máxima insolação, sendo que durante a noite parte desta energia é recuperada. Durante o dia o ar quente do meio interno da estufa solar é bombeado para dentro do feixe de tubos enterrados, que devolvem o ar mais frio na outra extremidade. Por outro lado, durante a noite o ar mais frio do meio interno da estufa é bombeado para dentro dos dutos, os quais efetuam troca térmica com o solo envolvente, que neste horário possui a temperatura mais elevada do sistema. O objetivo geral deste trabalho é resolver o problema transiente periódico e tridimensional da transferência simultânea de calor e umidade em solos não-saturados, que compõem o sistema trocador-armazenador de calor, utilizando a simulação numérica. Como objetivos secundários têm-se: o melhoramento do sistema de troca térmica, a quantificação da parcela de calor transportado pela difusão da umidade no solo e a análise dos campos de temperatura e de conteúdo de umidade, sendo que a análise dos campos de umidade permite verificar se existe formação de frentes de secagem significativas na vizinhança dos dutos durante os sucessivos períodos de aquecimento e resfriamento a que o sistema é submetido. Na resolução do problema em questão é empregado o modelo clássico de Philip e De Vries para a transferência simultânea de calor e massa em meios porosos insaturados Neste modelo, as equações de conservação de energia e massa obtidas trazem explicitamente as influências combinadas dos gradientes de temperatura e de conteúdo de umidade nos processos de transporte de calor e umidade. O sistema de equações diferenciais governantes do problema em questão é resolvido numericamente utilizando o Método dos Volumes Finitos e na discretização destas equações é usada uma integração temporal totalmente implícita. Todas as propriedades difusivas e termofísicas empregadas são consideradas variáveis com a temperatura e o conteúdo de umidade. Os dutos de seção transversal circular do sistema trocadorarmazenador de calor no solo, são modelados como dutos de seção transversal quadrada de área equivalente para que coordenadas cartesianas possam ser utilizadas nos modelos analisados. Neste trabalho são simulados quatro modelos computacionais associados ao sistema trocador-armazenador de calor no solo. Estes modelos são compostos por: um duto isolado, um duto com convecção, dois dutos isolados e dois dutos com convecção. A variação da temperatura do ar na entrada do(s) escoamento(s), assim como a temperatura do meio ambiente, para os modelos com convecção, é dada por uma senóide com uma amplitude de 14 ºC. No modelo de um duto isolado, são realizadas simulações utilizando várias combinações dos parâmetros do modelo em questão e os resultados, assim obtidos, são comparados com aqueles encontrados na literatura Visando melhorar o sistema de troca térmica dos modelos computacionais investigados, são selecionados valores e intervalos de valores recomendados para os parâmetros do modelo de um duto isolado. Para este modelo, com um diâmetro de 0,1 m, são escolhidos valores (ou intervalos de valores) recomendados: de 4 m/s para a velocidade do escoamento interno dentro do duto, de 0,25 para o conteúdo de umidade do solo, de 5 até 20 metros para o comprimento do duto e de 0,20 até 0,30 m para a distância entre centros do dutos. As simulações dos quatro modelos computacionais realizadas utilizando as várias combinações dos valores recomendados para os parâmetros destes modelos, mostrou que não há diferença significativa entre os valores de calor volumétrico armazenado no solo empregando a resolução acoplada das equações de energia e de massa e a resolução da equação da temperatura. Mesmo para os modelos de um e de dois dutos com convecção a diferença percentual encontrada foi insignificante. Finalmente, são apresentados e analisados os campos de temperatura e de conteúdo de umidade para os quatro modelos computacionais avaliados. Os perfis de temperatura e de conteúdo de umidade em diferentes horários mostraram que, durante o dia, o solo absorve calor dos escoamentos internos de ar e, uma vez que, junto à superfície dos dutos tem-se regiões de maior temperatura, há, conseqüentemente, uma migração da umidade nestas regiões. Durante a noite, ocorre o contrário, o solo fornece calor aos escoamentos dentro dos dutos, e, desta forma, as regiões próximas aos dutos apresentam níveis de conteúdo de umidade superiores ao inicial. Ainda, os perfis de conteúdo de umidade para todas as situações analisadas mostraram que, não há formação de frentes de secagem significativas nas proximidades dos dutos que compõem os quatro modelos computacionais avaliados.

Este trabalho tem o objetivo de realizar uma análise de sensibilidade dos parâmetros operacionais e de projeto de uma usina termelétrica. A análise de sensibilidade serve para indicar qual a importância dessas grandezas no resultado do sistema. A planta modelada é a da Usina Termelétrica Presidente Médici – Fase B, localizada em Candiota, RS, onde cada módulo é capaz de produzir até 160 MW de potência. Para esta planta é apresentado o equacionamento de cada componente do sistema. Destacam-se duas modelagens para a turbina a vapor, usadas para a simulação de seu comportamento em carga parcial, o da Elipse de Stodola e o de Schegliáiev. Ambos os modelos são comparados. Também são apresentados os modos de controle de carga por válvula de estrangulamento, por válvula de bocais e por pressão deslizante. A análise de sensibilidade é feita por meio de três métodos distintos, cujos resultados são comparados. O Método Diferencial utiliza as derivadas parciais para estimar a sensibilidade. O Método de Monte Carlo realiza uma série de avaliações do sistema, com os dados de entrada gerados randomicamente em cada avaliação e análise estatística dos resultados para avaliar a sensibilidade das respostas. É desenvolvido neste trabalho o Método da Transformada de Fourier, que gera os dados de forma senoidal e utiliza a transformada de Fourier para medir as sensibilidades. A medida da sensibilidade é feita através de dois índices. O índice de importância é a derivada parcial de um resultado em relação a um parâmetro, adimensionalizado pelas médias das variáveis, e o índice de sensibilidade é a composição fracional da variância de um resultado. Como conclusão, observa-se que os modelos para operação em carga parcial da turbina de Schegliáiev e de Stodola produzem praticamente os mesmos resultados. Os três métodos de análise de sensibilidade conduzem aos mesmos índices de sensibilidade no caso estudado. A ordenação da importância dos parâmetros estudados é idêntica quando se usa tanto o índice de importância quanto o de sensibilidade, mesmo seus valores sejam diferentes para cada parâmetros.

O objetivo deste trabalho é estudar os efeitos não lineares que resultam na dinâmica de um dímero rígido e amortecido, provido de rotação deslizando sobre um substrato periódico unidimensional. Simulações numéricas são realizadas para dinâmicas com amortecimento nos graus de liberdade translacional e rotacional, para lançamentos do dímero com velocidades translacionais finitas (estado transiente); com o dímero sujeito a uma força externa aplicada na coordenada do centro de massa e com temperaturas finitas (estado estacionário). As equações de movimento descritas em coordenadas de centro de massa e rotacional são caracterizadas por um acoplamento roto-translacional, este que é ativado em regimes distintos durante o deslizamento do dímero, resultando em transferências de energia entre suas coordenadas. A motivação deste estudo é tentar compreender os efeitos rotacionais emergentes na dinâmica do menor objeto possível que seja capaz de rotar, de forma a possibilitar simples contribuições para o entendimento da origem do atrito em escalas nanométricas. No estado de transiente foi possível obter soluções analíticas para os diferentes regimes da velocidade do centro de massa do dímero, bem como uma compreensão extensiva de sua dinâmica, realizada através de um estudo sobre o comportamento geral de sua inclinação angular e com comparações instrutivas com trabalhos anteriores realizados com um dímero vibrante. No estado estacionário foi possível determinar o comportamento no limiar de deslizamento (atrito estático) através do seu estudo como função da inclinação angular inicial. Três diferentes dinâmicas surgem neste estado, uma relacionada com um adátomo, outra com caráter pendular, e a última em que o dímero é capaz de deslizar realizando rotações completas. Nas três dinâmicas surgem efeitos não lineares, sendo que nas duas últimas acontece devido ao acoplamento rototranslacional, em que a natureza das dinâmicas são qualitativamente compreendidas através de uma análise do potencial efetivo como função da inclinação angular. Argumentos e considerações concisas são desenvolvidas com o âmbito de se formular simples expressões aproximadas que predizem com eficácia alguns dos limites que foram previstos. Também foi estudado o efeito da inclusão de temperaturas finitas no sistema, resultando numa atenuação no efeito não linear. / The aim of this work is to study the effects of rotation of the non linear friction of a damped dimer sliding on a 1D periodic substrate. Numerical simulations are performed with: a damping in the translational and rotational coordinate, throwing the dimer with a finite initial translational velocity (transient state); with the dimer subjected to an external force applied in the center of mass coordinate and with finite temperatures (steady state). The equations of motion in terms of center of mass and rotational coordinate show a rototranslational coupling, whose is activated for distinct regimes while the dimer is sliding, resulting in an energy transfer between the coordinates. The motivation of this work is to understand the rotational effects that emerges from the dynamics of the smallest object that can rotate, in order to achieve simple contributions to the understanding of the friction origin in nanometric scale. The first system was the transient state, and we were able to obtain analytical solutions of the center of mass velocity while the dimer is sliding in the three distinct regimes that arise. As well as an extensive comprehension of this dynamics, performed by studying the whole behavior of the angular inclination and instructive comparisons with the vibrating dimer. In the steady state system, an analysis of the behavior of the dimer at the threshold sliding (static friction) as function of the angular inclination allowed a complete understanding of the threshold effects. Three distinct dynamics arise from this state, one related to the adatom, another with a pendulous nature, and the last with the dimer performing complete rotations while is sliding. Non linear effects arise from the thee dynamics, for the two former cases the effects arise from roto-translational coupling, which are qualitatively understood by analysing the effective potential as function of the angular inclinations. Reasonable and proper arguments are developed aiming to construct simple and approximated expressions that predict with good efficiency some of the limit parameters values, that describes the transitions between different dynamics. Also, numerical simulations performed with finite temperature shows an attenuation in the non linear effects.

Neste trabalho estudamos um sistema de equações diferenciais parabólicas que modelam um processo de difusão-reação em duas dimensões da mistura molecular e reação química irreverssível de um só passo entre duas espécies químicas A e B para formar um produto P. Apresentamos resultados analíticos e computacionais relacionados à existência e unicidade da solução, assim como estimativas do erro local e global utilizando elementos finitos. Para os resultados analíticos usamos a teoria de semigrupos e o principio do m´aximo, e a simulação numérica é feita usando diferenças finitas centrais e o esquema simplificado de Ruge-Kutta. As estimativas do erro local para o problema semi-discretizado são estabelecidas usando normas de Sobolev, e para estimar o erro global usamos shadowing finito a posteriori. Os resultados computacionais obtidos mostram que o comportamento da solução está dentro do esperado e concorda com resultados da referências. Assim mesmo as estimativas do erro local e global são obtidas para pequenos intervalos de tempo e assumindo suficiente regularidade sobre a velocidade do fluído no qual realiza-se o processo. Destacamos que a estimativa do erro global usando shadowing finito é obtida sob hipóteses a posteriori sobre o operador do problema e o forte controle da velocidade numa vizinhança suficientemente pequena.

O presente trabalho é dedicado ao estudo de métodos de simulação para ciclos de Rankine. O trabalho é iniciado com a modelagem de um ciclo de Rankine simples e segue evoluindo para configurações mais complexas tal como o ciclo de Rankine com reaquecimento e regeneração. São adotadas as considerações mais convencionais da prática de projeto de centrais termelétricas cujos sistema térmicos baseiam-se no ciclo de Rankine, incluindo-se queda de pressão em tubulações do circuito além de outras perdas. Em seguida, são estabelecidas as expressões matemáticas que possibilitam a determinação das propriedades termodinâmicas da água em seus mais diversos estados ao longo do ciclo. Por último, são desenvolvidos métodos de simulação, chamados neste trabalho de Substituição Sucessiva e Bloco Único, que caracterizam-se pela resolução simultânea do conjunto de equações algébricas dos ciclos elaborados. As simulações são efetuadas através de programas escritos na linguagem Fortran. Os métodos de simulação são aplicados para a obtenção dos resultados considerados mais importantes na análise de sistemas térmicos de potência, tais como rendimento térmico do ciclo, título na saída da turbina, vazões mássicas pelo sistema, potência nas bombas e calor trocado no gerador de vapor e no condensador Na maioria das simulações, estes resultados apresentam-se como funções da: (1) potência elétrica requerida, eficiência isentrópica e pressões na turbina; (2) eficiência térmica, pressão e temperatura no gerador de vapor; (3) pressão e grau de subresfriamento do líquido saturado no condensador e (4) eficiência isentrópica das bombas. São obtidos os mesmos resultados para os métodos de simulação utilizados. O método da Substituição Sucessiva apresentou menor tempo computacional, principalmente para configurações de ciclo mais complexas. Uma aplicação alternativa do método de Bloco Único demonstrou ser inconveniente para ciclos de configurações mais complexas devido ao elevado tempo computacional, quando todas as equações de cálculo das propriedades termodinâmicas são incluídas no sistema de equações a ser resolvido. Melhores rendimentos térmicos e título na saída da turbina foram obtidos para configurações de ciclo de Rankine com reaquecimento e regeneração.

Os sistemas de aquecimento através da energia solar são uma fonte interessante de energia acessível e distribuída. No Brasil, a utilização de coletores solares planos já é bem conhecida e difundida. Os coletores solares de tubo evacuado até pouco tempo tinham um custo considerado elevado para o nosso mercado, mas atualmente esta tecnologia está cada vez mais difundida no mercado nacional com custos menores. Para se entender se o custo adicional deste coletor quando comparado aos tradicionais coletores planos é vantajoso, são necessários estudos comparativos para determinar as reais vantagens e quando as mesmas podem proporcionar uma maior economia financeira à instalação. Esta dissertação compara as tecnologias atuais de coletores planos com os de tubo evacuado disponíveis no mercado nacional, mostrando comparações em diferentes regiões do país e analisando se as vantagens encontradas são interessantes do ponto de vista do custo da instalação. Para realizar estas analises foram desenvolvidos testes para caracterização dos tubos destes coletores, de forma a determinar os parâmetros de desempenho dos mesmos. Também foi preparada uma bancada de testes para os coletores solares e foram realizados ensaios de comportamento térmico. Conhecendo os resultados foi possível adequar e utilizar um software de simulação desenvolvido no Laboratório de Energia Solar da Universidade Federal do Rio Grande do Sul. Este software visa a simulação de sistemas de aquecimento de água com diferentes configurações em qualquer região do país. Com a implementação dos coletores de tubo evacuado, e com simulações do seu desempenho em diferentes regiões, é possível comparar o seu comportamento em relação ao coletor plano. Nesta dissertação as curvas de desempenho foram comparadas em diferentes situações, onde foi possível verificar em que momentos o coletor de tubos pode ser melhor que o plano. Também foram realizadas comparações entre a área ocupada e o custo destes coletores para instalações em 3 diferentes regiões do país. Com estas comparações ficou evidente que para climas mais frios o coletor de tubos evacuados pode apresentar melhor desempenho que os coletores planos, mas o seu custo ainda é muito maior que os coletores planos nacionais. / The solar water heating systems are a interesting source of accessible and distributed energy. In Brazil, the flat plate solar collector use is well known. The vaccum tube solar collectors until little time ago had an elevated cost for our market, but, today, the technology is available in our national market with lower costs. And to understand if this costs, that are higher than the flat plane collectors are worth it, some studies are necessary to understand the real benefits of this technology. This work aims to compare both technologies available in the national market, showing comparisons in different regions of the country and studying if the advantages are interesting in the economical point of view. To achieve that some tests were developed in order to obtain the performance parameters of this collectors. Also an external bench was developed to test this collectors according to the Brazilian standard. With the results and knowing the behavior of this collectors it was possible to program the simulation software that was developed in the Laboratório de Energia Solar of the Universidade Federal do Rio Grande do Sul. This software aims at the simulation of water heating systems in different regions of the country. With the implementation of the vacuum collectors in the software, it was possible to simulate different climate conditions over the country and compare the behavior with the flat plate collectors. In this work the performance curves were compared in different climate situations, so that was possible to understand when the vaccum collector can be better. Also some comparisons where made between the occupied area and the costs in three different regions of the country. With the results of this comparisons it is evident that for cold climates the vaccum collector can obtain a better performance than the flat plate collectors, but its cost is still higher.

O principal objetivo dessa dissertação consiste nas simulações numéricas de um escoamento de um fluido viscoplástico usando o método de elementos finitos. A modelagem mecânica empregada está baseada nas equações conservativas de massa e momentum acopladas ao modelo constitutivo recentemente proposto por Souza Mendes e Dutra – aqui chamado de SMD. A modelagem mecânica foi aproximada via método de Galerkin Mínimos-Quadrados (GLS), que é construído para superar a falha numérica do método clássico de Galerkin para formulações mistas sem danificar a sua consistência. Primeiro, a fim de checar o código computacional, um escoamento em uma cavidade forçada. Em seguida, dois escoamentos de fluidos viscoplásticos foram considerados: um fluido SMD escoando através de uma expansão seguida de uma contração abrupta e outro em torno de um cilindro circular. Os resultados numéricos gerados forneceram informações físicas dos fenômenos de viscoplasticidade estudados e se mostraram de acordo com a literatura. / The main objective of this paper consists of numerical simulations of viscoplastic fluid flows using the finite element method. The mechanical modeling employed is based on the mass and momentum conservative equations coupled with the constitutive model recently proposed by Souza Mendes and Dutra - herein only named by SMD. The mechanical modeling was approximated via the Galerkin Least-Squares method (GLS) which is built to overcome the numerical shortcoming of the classical Galerkin method for mixed formulations without damaging its consistency. First, in order to check the computational code, a Newtonian leaky cavity flow. Next, two viscoplastic fluid flows has been considered: a SMD fluid flow through a sudden expansion followed by a contraction geometry and around a circular cylinder. The generated numerical results gave some physical information on the viscoplastic phenomena studied and are in accordance with the literature.

A estimativa de parâmetros térmicos relacionados à Soldagem por Fricção Linear (também conhecida como Soldagem por Fricção e Mistura Mecânica) de uma placa de alumínio AA 2195- T8 é estudada nesta dissertação. A determinação dos parâmetros é tratada como um problema de otimização, no qual a função objetivo corresponde a uma função erro entre temperaturas medidas numericamente e temperaturas calculadas, para certos valores da taxa de calor gerada pela fricção, do coeficiente de transferência calor por convecção natural e do coeficiente de transferência de calor entre a placa e a base de suporte. A distribuição transiente de temperaturas sobre a placa é determinada pela solução da equação da condução de calor tridimensional transiente, a qual é resolvida pelo método dos volumes finitos. A minimização da função objetivo é realizada com o auxílio do algoritmo Otimização Extrema Generalizada (GEO), um método evolucionário que pode lidar com diversos tipos de problemas de otimização. A avaliação da sensibilidade da leitura de temperatura é feita com relação ao posicionamento dos sensores, com o objetivo de descobrir as melhores posições para aquisição de dados, e também a determinação do parâmetro t (parâmetro utilizado no método GEO que deve ser ajustado para cada tipo de problema) que melhor se adapta para cada um dos conjuntos de medição de temperatura. Além disso, realiza-se um estudo no qual é analisado um caso em que a fonte de calor possui forma de distribuição desconhecida, sem a utilização de equações para a sua descrição, desta maneira mais parâmetros devem ser estimados pela análise inversa, assim o processo de otimização é mais complexo (cinco valores de taxa de calor e dois coeficientes de transferência de calor, ao invés de uma taxa de calor e dois coeficientes). Ainda, para simular medições de dados reais, os valores de temperatura obtidos a partir da solução numérica para valores específicos da taxa de calor e dos coeficientes de transferência de calor são perturbados com ruídos de acordo com desvios-padrão típicos dessa forma de medição. Este trabalho demonstra que a aproximação inversa proposta pode ser um modo muito efetivo para avaliar e predizer os parâmetros que governam o processo de transferência de calor na Soldagem por Fricção Linear, um importante passo para o seu controle em tempo real. / The estimation of thermal parameters related to Friction Stir Welding of a AA 2195-T8 plate is studied in this dissertation. The determination of the parameters is carried out by means of an optimization problem, in which the objective function corresponds to an error function between the numerically measured temperature and the temperature computed for each estimated values of the heat rate input, the heat transfer coefficient on the bottom surface, and the natural convection heat transfer coefficient. The time-dependent temperature distribution on the plate is determined by the solution of the three-dimensional transient state conduction equation, which is solved by the control-volum method. The minimization of the objective function is accomplished with the aid of the Generalized Extremal Optimization (GEO) method, an evolutionary method that can deal with several types of optimization problems. The evaluation of the temperature reading sensitivity is carried out, in relation to the readers locations on the plate, in order to determine the best positions that can be used to acquire datas, and the determination of the t parameter (this parameter is used in the GEO method and it must be adjusted for each type of problem) that best fits with each one of the temperature readers assembly, allowing the determination of the best temperature readers assembly. In addition, it is executed a study that considers a case where the heat rate input has unknown profile distribution, without using equations for its description, in this manner more parameters have to be estimated by the inverse analysis, so the optimization process is more complicated (five heat rate inputs and two heat transfer coefficients, instead of one heat rate input and two heat transfer coefficients). Furthermore, to simulate real-data measurements, the temperature inputs, obtained from a numerical solution for specific values of the heat rate input and the heat transfer coefficients, were perturbed with noises according to the standard deviation of the measurement procedure. This work demonstrates that the proposed inverse approach can be a very effective way to evaluate and predict the parameters that govern the heat transfer process in Friction Stir Welding, an important step to real-time control of this process.

A otimização de processos tem como objetivo a maximização, ou minimização, de uma função objetivo que descreve o problema a ser resolvido. Para tanto, dependendo do problema e visando a obtenção da melhor resposta no menor tempo, podem ser utilizadas técnicas estocásticas ou determinísticas. Atualmente, o emprego de otimizadores não mais se restringe às aplicações da indústria de processos e manufaturas, mas também em outras áreas, como, por exemplo, a de energias renováveis. Esse fato se deve à conscientização acerca das consequências da utilização combustíveis fósseis e da busca pelo maior custo-benefício em projetos da área. Este trabalho tem como objetivo estudas o posicionamento de aerogeradores em um parque eólico utilizando otimização de processos. Para tal, foi desenvolvido um modelo baseado em Katic et al. (1986) para a descrição do déficit de velocidade do vento, aplicando o conceito de blending functions para combinar as perdas pelo efeito de esteira em todas as direções de vento incidente. O resultado é um modelo de rápida execução e qualitativamente correto. Comparou-se, então, soluções de otimizadores de busca direta global do tipo estocástico (PSO) e determinístico (DIRECT), escolhidos em função da não-regularidade da função objetivo. Considerou-se, também, o valor de 15% de perda de energia em relação à geração ideal como sendo o máximo aceitável para as respostas. O DIRECT alcançou soluções viáveis para até 8 turbinas, enquanto o PSO solucionou o problema para até 13 aerogeradores. Os resultados do PSO com 80 partículas e 100 iterações, para posicionamento de 13 aerogeradores, foi comparado com os resultados mostrados por Mosetti et al. (1994), mostrando uma necessidade de área de terreno cerca de cinco vezes menor para cada turbina. Avaliou-se, então, a influência os parâmetros algoritmo PSO, tanto os internos quantos os externos (número de partículas e iterações), em cenários de média e alta complexibilidade para a mesma quantidade de avaliações da função objetivo, respectivamente. As soluções encontradas mostraram que tais parâmetros tem grande efeito nas respostas dos otimizadores. / The aim of process optimization is maximize, or minimize, an objective function that describes a problem to be solved. To do that, depending on the problem and aiming the best solution in a shortest time, stochastic or deterministic techniques can be used. Currently, the use of optimizers no longer is restricted to process and manufacturing industrial problems, but in many others fields, like, for example, the renewable energy. This is due to the awareness about the consequences of using fossils fuels and looking for the best cost-effective in the projects of this area. The objective of this work is study the set of wind turbines in a wind park using process optimization. To do that, it was developed a model based on Katic et al. (1986) to describe the wind speed decay, using the knowledge of blending functions to add the wake losses in all wind directions. The result is a fast performance and qualitatively correct model. Then, global direct search optimizers solutions, a stochastic (PSO) and a deterministic (DIRECT), were compared. Those were chosen due the non-regular objective function. It was considerer 15% as the maximum acceptable energy losses, compared to the ideal, of the positioning answeres. The DIRECT reached viable solutions up to 8 turbines while the PSO solved the problem up to 13 turbines. The results of the PSO with 80 particles and 100 iteractions to positioning 13 turbines, was compared with the results shown by Mosetti et al. (1994) resulting in about 5 times less land required to each turbine. Then, it was evaluated the effect os the PSO parameters, inner and outer (number of particles and iteractions), in medium and high complexity scenarios to the same amount of objective functions evaluations. The solutions reached shown that those parameters have large influence in the optimizers solutions.