Um método SN híbrido direto para cálculos de sistemas combustível-moderador em geometria unidimensional

Desenvolvemos nesta dissertação um método híbrido direto para o cálculo do fator de desvantagem e descrição da distribuição do fluxo de nêutrons em sistemas combustível-moderador. Na modelagem matemática, utilizamos a equação de transporte de Boltzmann independente do tempo, considerando espalhamento linearmente anisotrópico no modelo monoenergético e espalhamento isotrópico no modelo multigrupo, na formulação de ordenadas discretas (SN), em geometria unidimensional. Desenvolvemos nesta dissertação um método híbrido direto para o cálculo do fator de desvantagem e descrição da distribuição do fluxo de nêutrons em sistemas combustível-moderador. Na modelagem matemática, utilizamos a equação de transporte de Boltzmann independente do tempo, considerando espalhamento linearmente anisotrópico no modelo monoenergético e espalhamento isotrópico no modelo multigrupo, na formulação de ordenadas discretas (SN), em geometria unidimensional. Descrevemos uma análise espectral das equações de ordenadas discretas (SN)a um grupo e a dois grupos de energia, onde seguimos uma analogia com o método de Case. Utilizamos, neste método, quadraturas angulares diferentes no combustível (NC) e no moderador (NM), onde em geral assumimos que NC > NM . Condições de continuidade especiais que acoplam os fluxos angulares que emergem do combustível (moderador) e incidem no moderador (combustível), foram utilizadas com base na equivalência entre as equações SN e PN-1, o que caracteriza a propriedade híbrida do modelo proposto. Sendo um método híbrido direto, utilizamos as NC + NM equações lineares e algébricas constituídas pelas (NC + NM)/2 condições de contorno reflexivas e (NC + NM)/2 condições de continuidade para determinarmos as NC + NM constantes. Com essas constantes podemos calcular os valores dos fluxos angulares e dos fluxos escalares em qualquer ponto do domínio. Apresentamos resultados numéricos para ilustrar a eficiência e a precisão do método proposto.

Efeitos do controle da jaqueira, Artocarpus heterophyllus L., sobre a comunidade de pequenos mamíferos e a rede de dispersão de sementes na Ilha Grande, RJ

Atualmente observa-se uma expressiva perda de biodiversidade global resultante de atividades antrópicas, sendo a introdução de espécies exóticas uma das mais impactantes. A jaqueira Artocarpus heterophyllus é uma espécie exótica introduzida no Brasil durante o período colonial, sendo considerada invasora em diversas localidades. Na Mata Atlântica invade áreas de mata aberta e de borda, habitualmente associadas a ambientes antrópicos. Na Ilha Grande é encontrada em grande abundância em decorrência do histórico de ocupação humana. Para compreender como a mastofauna responde a presença da jaqueira, o Laboratório de Ecologia de Mamíferos da Universidade do Estado do Rio de Janeiro (UERJ) vem desenvolvendo um estudo ao longo de seis anos nos arredores da Vila Dois Rios, localizada na face oceânica da Ilha Grande. A partir dos resultados prévios iniciou-se uma segunda etapa do estudo no mesmo local que buscou avaliar diferentes métodos de controle das jaqueiras. O presente estudo é uma continuação direta desses dois trabalhos anteriores e teve como objetivo acompanhar as respostas da comunidade de pequenos mamíferos no período imediatamente posterior ao controle. Durante 18 meses foram amostradas bimestralmente 18 grades, 10 aonde foi efetuado o controle das jaqueiras e 8 aonde não foi constatada a presença desta árvore. Em cada grade foram colocadas 11 armadilhas de captura viva sendo banana a isca utilizada. Os mamíferos capturados foram medidos e suas fezes coletadas. A quantidade de jacas em cada área também foi anotada bimensalmente. As fezes foram analisadas em laboratório e as sementes encontradas identificadas. Os resultados obtidos indicam que a influência de A. heterophyllus sobre a estrutura da comunidade de pequenos mamíferos foi menor após o tratamento de controle. A única espécie que parece ainda responder a abundância de jaqueiras é o roedor Trinomys dimidiatus, que apresentou densidades mais elevadas nas áreas em tratamento, porém mais próximas a resultados obtidos para espécies congêneres em áreas pouco antropizadas. Utilizando uma abordagem de redes complexas observamos que, embora T. dimidiatus seja a espécie mais abundante em termos de número de indivíduos, o gambá Didelphis aurita parece ser a espécie de mamífero mais importante para dispersão de sementes nativas, aparecendo como espécie com maior número de conexões com espécies de sementes nas redes contruídas para as áreas sem jaqueiras e com jaqueiras antes e após o tratamento. Finalmente, a partir dos dados obtidos criamos um modelo matemático para a população de T. dimidiatus dos arredores da Vila Dois Rios, baseado em um crescimento logístico. Os resultados do modelo proposto se mostraram correlacionados com os dados de abundância reais, de modo que ele parece ser um simulador adequado da população local.

Simulações de problemas inversos com aplicações em engenharia nuclear usando técnicas de transporte de partículas neutras monoenergéticas na formulação unidimensional de ordenadas discretas

Neste trabalho, três técnicas para resolver numericamente problemas inversos de transporte de partículas neutras a uma velocidade para aplicações em engenharia nuclear são desenvolvidas. É fato conhecido que problemas diretos estacionários e monoenergéticos de transporte são caracterizados por estimar o fluxo de partículas como uma função-distribuição das variáveis independentes de espaço e de direção de movimento, quando os parâmetros materiais (seções de choque macroscópicas), a geometria, e o fluxo incidente nos contornos do domínio (condições de contorno), bem como a distribuição de fonte interior são conhecidos. Por outro lado, problemas inversos, neste trabalho, buscam estimativas para o fluxo incidente no contorno, ou a fonte interior, ou frações vazio em barras homogêneas. O modelo matemático usado tanto para os problemas diretos como para os problemas inversos é a equação de transporte independente do tempo, a uma velocidade, em geometria unidimensional e com o espalhamento linearmente anisotrópico na formulação de ordenadas discretas (SN). Nos problemas inversos de valor de contorno, dado o fluxo emergente em um extremo da barra, medido por um detector de nêutrons, por exemplo, buscamos uma estimativa precisa para o fluxo incidente no extremo oposto. Por outro lado, nos problemas inversos SN de fonte interior, buscamos uma estimativa precisa para a fonte armazenada no interior do domínio para fins de blindagem, sendo dado o fluxo emergente no contorno da barra. Além disso, nos problemas inversos SN de fração de vazio, dado o fluxo emergente em uma fronteira da barra devido ao fluxo incidente prescrito no extremo oposto, procuramos por uma estimativa precisa da fração de vazio no interior da barra, no contexto de ensaios não-destrutivos para aplicações na indústria. O código computacional desenvolvido neste trabalho apresenta o método espectronodal de malha grossa spectral Greens function (SGF) para os problemas diretos SN em geometria unidimensional para gerar soluções numéricas precisas para os três problemas inversos SN descritos acima. Para os problemas inversos SN de valor de contorno e de fonte interior, usamos a propriedade da proporcionalidade da fuga de partículas; ademais, para os problemas inversos SN de fração de vazio, oferecemos a técnica a qual nos referimos como o método físico da bissecção. Apresentamos resultados numéricos para ilustrar a precisão das três técnicas, conforme descrito nesta tese.

Problemas inversos em processos difusivos com retenção

Um Estudo para a solução numérica do modelo de difusão com retenção, proposta por Bevilacqua et al. (2011), é apresentado, bem como uma formulação implícita para o problema inverso para a estimativa dos parâmetros envolvidos na formulação matemática do modelo. Através de um estudo minucioso da análise de sensibilidade e do cálculo do coeficiente de correlação de Pearson, são identificadas as chances de se obter sucesso na solução do problema inverso através do método determinístico de Levenberg-Marquardt e dos métodos estocásticos Algoritmo de Colisão de Partículas (Particle Collision Algorithm - PCA) e Evolução Diferencial (Differential Evolution - DE). São apresentados os resultados obtidos através destes três métodos de otimização para três casos de conjunto de parâmetros. Foi observada uma forte correlação entre dois destes três parâmetros, o que dificultou a estimativa simultânea dos mesmos. Porém, foi obtido sucesso nas estimativas individuais de cada parâmetro. Foram obtidos bons resultados para os fatores que multiplicam os termos diferenciais da equação que modela o fenômeno de difusão com retenção.

Custo efetividade do uso do Peguinterferon alfa 2a combinado com Ribavirina no tratamento de respondedores virológicos lentos coinfectados com VHC/HIV

No mundo, as hepatites decorrentes de infecções virais têm sido uma das grandes preocupações em saúde pública devido a seu caráter crônico, curso assintomático e pela sua capacidade de determinar a perda da função hepática. Com o uso em larga escala de medicamentos antirretrovirais, a doença hepática relacionada à infecção pelo vírus da hepatite C (VHC) contribuiu para uma mudança radical na história natural da infecção pelo vírus da imunodeficiência humana (HIV). Não se sabe ao certo o peso da coinfecção VHC/HIV no Brasil, mas evidências apontam que independentemente da região geográfica, esses indivíduos apresentam maiores dificuldades em eliminar o VHC após o tratamento farmacológico, quando comparados a monoinfectados. No âmbito do SUS, o tratamento antiviral padrão para portadores do genótipo 1 do VHC e do HIV é a administração de peguinterferon associado à Ribavirina. Quanto ao período de tratamento e aos indivíduos que devem ser incluídos, os dois protocolos terapêuticos mais recentes possuem divergências. A diretriz mais atual preconiza o tratamento de indivíduos respondedores precoces somados a respondedores virológicos lentos, enquanto a diretriz imediatamente anterior exclui na 12 semana indivíduos que não respondem completamente. Com base nessa divergência, esse estudo objetivou avaliar o custo-efetividade do tratamento contra o VHC em indivíduos portadores do genótipo 1, coinfectados com o HIV, virgens de tratamento antiviral, não cirróticos e imunologicamente estabilizados, submetidos às regras de tratamento antiviral estabelecidos pelas duas mais recentes diretrizes terapêuticas direcionadas ao atendimento pelo SUS. Para tal, foi elaborado um modelo matemático de decisão, baseado em cadeias de Markov, que simulou a progressão da doença hepática mediante o tratamento e não tratamento. Foi acompanhada uma coorte hipotética de mil indivíduos homens, maiores de 40 anos. Adotou-se a perspectiva do Sistema Único de Saúde, horizonte temporal de 30 anos e taxa de desconto de 5% para os custos e consequências clínicas. A extensão do tratamento para respondedores lentos proporcionou incremento de 0,28 anos de vida ajustados por qualidade (QALY), de 7% de sobrevida e aumento de 60% no número de indivíduos que eliminaram o VHC. Além dos esperados benefícios em eficácia, a inclusão de respondedores virológicos lentos mostrou-se uma estratégia custo-efetiva ao alcançar um incremental de custo efetividade de R$ 44.171/QALY, valor abaixo do limiar de aceitabilidade proposto pela Organização Mundial da Saúde OMS - (R$ 63.756,00/QALY). A análise de sensibilidade demonstrou que as possíveis incertezas contidas no modelo são incapazes de alterar o resultado final, evidenciando, assim, a robustez da análise. A inclusão de indivíduos coinfectados VHC/HIV respondedores virológicos lentos no protocolo de tratamento apresenta-se, do ponto de vista fármaco-econômico, como uma estratégia com relação de custoefetividade favorável para o Sistema Único de Saúde. Sua adoção é perfeitamente compatível com a perspectiva do sistema, ao retornar melhores resultados em saúdeassociados a custos abaixo de um teto orçamentário aceitável, e com o da sociedade, ao evitar em maior grau, complicações e internações quando comparado à não inclusão.

Teoria de controle ótimo com aplicações a sistemas biológicos

Neste trabalho apresentamos as etapas para a utilização do método da Programação Dinâmica, ou Princípio de Otimização de Bellman, para aplicações de controle ótimo. Investigamos a noção de funções de controle de Lyapunov (FCL) e sua relação com a estabilidade de sistemas autônomos com controle. Uma função de controle de Lyapunov deverá satisfazer a equação de Hamilton-Jacobi-Bellman (H-J-B). Usando esse fato, se uma função de controle de Lyapunov é conhecida, será então possível determinar a lei de realimentação ótima; isto é, a lei de controle que torna o sistema globalmente assintóticamente controlável a um estado de equilíbrio. Como aplicação, apresentamos uma modelagem matemática adequada a um problema de controle ótimo de certos sistemas biológicos. Este trabalho conta também com um breve histórico sobre o desenvolvimento da Teoria de Controle de forma a ilustrar a importância, o progresso e a aplicação das técnicas de controle em diferentes áreas ao longo do tempo.

Análise da resposta dinâmica de passarelas de pedestres considerando-se uma modelagem probabilística do caminhar humano.

Passarelas de pedestres com arquitetura moderna, esbeltas e leves são uma constante nos dias atuais, apresentando grandes vãos e novos materiais. Este arrojo arquitetônico tem gerado inúmeros problemas de vibrações excessivas, especialmente sobre passarelas mistas (aço-concreto). As normas e recomendações de projeto consideram, ainda, que as forças induzidas pelo caminhar humano são determinísticas. Todavia, o caminhar humano e as respectivas forças dinâmicas geradas apresentam comportamento randômico. Deste modo, o presente trabalho de pesquisa objetiva contribuir com os projetistas estruturais, a partir do emprego de uma abordagem probabilística para avaliação do estado limite de utilização deste tipo de estrutura, associado a vibrações excessivas que podem vir a causar desconforto humano. Para tal, utiliza-se como modelo estrutural uma passarela de pedestres mista (aço-concreto) construída no campus do Instituto de Traumatologia e Ortopedia (INTO), na cidade do Rio de Janeiro. Com base na utilização dos métodos probabilísticos, torna-se possível determinar a probabilidade dos valores das acelerações de pico da estrutura ultrapassarem ou não os critérios de conforto humano estabelecidos em normas e recomendações de projeto. Os resultados apontam para o fato de que os valores das acelerações de pico calculadas com base exclusivamente nos métodos determinísticos podem ser superestimados em algumas situações de projeto.

Modelagem matemática da transferência de calor numa placa plana sob o efeito de uma fonte pontual externa de radiação térmica.

Este trabalho apresenta uma modelagem matemática para o processo de aquecimento de um corpo exposto a uma fonte pontual de radiação térmica. O resultado original que permite a solução exata de uma equação diferencial parcial não linear a partir de uma seqüência de problemas lineares também é apresentado. Gráficos gerados com resultados obtidos pelo método de diferenças finitas ilustram a solução do problema proposto.

Análise das dosagens do fator de crescimento endotelial vascular no plasma e fluidos peritoneais de pacientes com câncer epitelial de ovário

A carcinogênese epitelial ovariana tem sido foco de estudos científicos em todo o mundo desenvolvido. A angiogênese tumoral ovariana é um processo multifatorial que resulta em vários produtos pró-angiogênicos. Entre eles, o fator de crescimento vascular endotelial (VEGF) é predominante. Os objetivos deste estudo foram relacionar as dosagens do VEGF dos fluidos peritoneais, do plasma periférico e do infundíbulo pélvico aos níveis de citorredução em pacientes operadas de adenocarcinoma epitelial de ovário (CEO); formular um modelo probabilístico de citorredução e utilizar estas dosagens para estimar a probabilidade do desfecho de citorredução. Além disto, foi criada uma nova variável chamada carga de VEGF. Pelo procedimento step-wise a citorredução foi melhor descrita pela carga de VEGF, mas faltou Normalidade aos resíduos, não sendo possível a adequação de um modelo matemático. Porém, a curva ROC, forneceu uma área sob a curva de 0,84, com sensibilidade de 71,4 % e especificidade variando de 69,5 a 73,9%. O ponto de corte ótimo foi 15,52 log de picograma de carga de VEGF. A odds-ratio (OR) calculada para citorredução ótima descrita pela carga de VEGF foi de 11 (IC= 2,59 ; 46,78). No grupo com estágio avançado (III e IV), a OR foi de 6 (IC= 1,15; 31,22). Apesar do pequeno número de casos, esta nova variável pode vir a ser um auxiliar na determinação de situações onde cirurgia citorredutora deixa de ser a pedra fundamental do tratamento primário do CEO e a indução quimioterápica passe a ter o principal papel na citorredução química antes da cirugia nestes casos.

Impacto do aterro controlado Morro do Céu na atmosfera

O principal objetivo deste estudo é estimar o alcance dos poluentes emitidos pelo Aterro Controlado do Morro do Céu (ACMC), localizado em Niterói e as possíveis áreas afetadas. Utilizou-se um modelo matemático para o estudo do transporte e difusão dos poluentes que fornece a concentração nas vizinhanças, a partir dos dados de emissões, meteorológicos e topográficos. Foram analisados e tratados os dados meteorológicos e a topografia da região de alcance da pluma. As amostras foram coletadas por uma bomba de ar a uma vazão de 500 mL min-1 por um período de 10 min, em dois dias e horários diferentes. A câmara de fluxo empregada foi de 30 L e de PVC, no formato cilíndrico, cuja boca foi emborcada no solo por 5 cm e a amostragem realizada por uma válvula posicionada no fundo. Foram empregados como meios de amostragem cartuchos de carvão ativo de leito duplo e as análises químicas foram realizadas por cromatografia de fase gasosa com detecção por ionização na chama. Os resultados encontrados indicam que o ACMC pode conter quantidade significativa de vários compostos químicos, principalmente compostos orgânicos voláteis, impactando a vegetação, prejudicando a saúde e a qualidade de vida da população do seu entorno. Foram obtidos valores de compostos orgânicos voláteis acima de 400 g m-3 a menos que 500 m do ACMC. Considerando que o ACMC não será desativado em um tempo próximo, creches, escolas e postos de saúde devem ser deslocados para uma área onde o impacto dos poluentes atmosféricos seja menor, contribuindo para a melhoria da qualidade de vida da população

Um modelo físico-matemático para escoamentos em meios porosos com transição insaturado-saturado.

Neste trabalho é apresentada uma nova modelagem matemática para a descrição do escoamento de um líquido incompressível através de um meio poroso rígido homogêneo e isotrópico, a partir do ponto de vista da Teoria Contínua de Misturas. O fenômeno é tratado como o movimento de uma mistura composta por três constituintes contínuos: o primeiro representando a matriz porosa, o segundo representando o líquido e o terceiro representando um gás de baixíssima densidade. O modelo proposto possibilita uma descrição matemática realista do fenômeno de transição insaturado/saturado a partir de uma combinação entre um sistema de equações diferenciais parciais e uma desigualdade. A desigualdade representa uma limitação geométrica oriunda da incompressibilidade do líquido e da rigidez do meio poroso. Alguns casos particulares são simulados e os resultados comparados com resultados clássicos, mostrando as consequências de não levar em conta as restrições inerentes ao problema.

Análise dos Modos de Plasmon em Fibras Ópticas com Contraste Arbitrário de Índices de Refração.

Neste trabalho, é feita a análise dos modos de plasmon que se propagam em um filme metálico que cobre uma fibra óptica generalizada. Os modos de plasmon estudados são: Fuga pela Cobertura (lcv), Ligado Simétrico (Sb), Fuga pelo Núcleo (lcr) e Ligado Assimétrico (ab). Os filmes metálicos, para efeito de comparação, utilizados neste trabalho, são: a prata, o ouro e o paládio. Desenvolveu-se um modelo matemático do fenômeno eletromagnético e um software, que gerou um banco de dados que facilitasse a análise de estruturas, com diversas combinações de parâmetros. Com o banco de dados, foram obtidos diversos gráficos, que permitiram: analisar os modos de plasmon, verificar a atenuação das ondas e o comportamento do campo eletromagnético em cada região da estrutura. As confrontações entre as estruturas com filmes de: prata, ouro e paládio, permitiram concluir que aquelas elaboradas com os filmes de prata e de ouro são as que apresentam menores perdas, portanto, as recomendadas na confecção de sensores. Como a prata é mais acessível que o ouro, aconselha-se a sua utilização. A análise e os resultados deste trabalho são originais na literatura especializada.

Influência da Velocidade, do Espaçamento e do Número de Veículos sobre a Resposta Dinâmica de Pontes Rodoviárias de Concreto Armado.

Nesta dissertação são avaliados os efeitos dinâmicos provenientes da travessia de comboios de veículos sobre o tabuleiro irregular de obras de arte rodoviárias de concreto armado. O modelo matemático empregado para simular o comportamento do sistema veículo-ponte considera a participação da massa e da rigidez das viaturas na definição das freqüências do sistema e, conseqüentemente, a força de interação entre os veículos e a ponte é afetada pela flexibilidade desta. A ponte é modelada a partir do emprego de elementos finitos de barra unidimensionais e discretizado com massas concentradas e flexibilidade distribuída. O modelo de veículo empregado baseia-se no veículo TB-12 preconizado pela norma brasileira NBR 7188. Este veículo é simulado por sistemas de massas, molas e amortecedores sendo descrito por graus de liberdade à translação e rotação no plano. As irregularidades da pista são definidas por um modelo não-determinístico com base na densidade espectral do pavimento. O carregamento sobre a ponte é constituído por sucessões de veículos deslocando-se com velocidade constante sobre a obra. Devido à própria natureza das irregularidades da pista e do comboio de veículos, atenção especial é concentrada na fase permanente da resposta do sistema. São estudadas as respostas de dois modelos estruturais existentes, com base em tabuleiros isostáticos, em concreto armado, com e sem balanços, em seção do tipo T e duplo T, respectivamente, em termos de deslocamentos e esforços nas seções onde ocorrem os efeitos máximos. As conclusões do trabalho versam sobre a influência da velocidade, espaçamento e do número de veículos, referentes a situações distintas de carregamento, no que tange a resposta dinâmica das pontes rodoviárias de concreto armado. A magnitude dos efeitos dinâmicos associados à interação dos veículos com o pavimento irregular também é investigada.

Estudo do equilíbrio líquido-líquido e da estabilidade de sistemas envolvendo biodiesel, álcoois e água

O estudo da estabilidade de emulsões de biocombustíveis destaca-se, dentre outros aspectos, pela necessidade de identificação das características destas. Esta necessidade reforça-se pela dificuldade de detecção de biodiesel e suas misturas quando estão fora de especificação, além de se aplicar aos processos de separação e purificação do biodiesel após a transesterificação. O entendimento deste problema poderá ser alcançado a partir de estudos de estabilidade associados aos modelos termodinâmicos de equilíbrio líquido-líquido (modelo NRTL) em sistemas ternários, que serão a base da modelagem proposta neste trabalho. O levantamento de dados envolve técnicas de observação direta em titulação, e ensaios de estabilidade e tensiometria, com algumas alterações das técnicas clássicas. Desta forma, o objetivo principal do trabalho é propor uma maneira de identificar e caracterizar as interfaces das emulsões formadas e estudar a estabilidade das fases envolvidas. Além disso, estuda-se a introdução de alcoóis como agentes tensoativos para viabilizar a utilização de misturas de combustíveis, mesmo na região de mais de uma fase. Para o desenvolvimento do trabalho, emulsões envolvendo os constituintes básicos serão combinadas entre si em diferentes proporções. A etapa inicial constitui-se de um mapeamento das misturas e proporções que formam as emulsões, baseadas nas curvas de equilíbrio. Em seguida, faz-se o estudo de caracterização e estabilidade destes sistemas. Com estas informações é possível mapear o sentido de incremento de tensão na interface e instabilidade da separação de fases das emulsões em relação aos diagramas ternários correspondentes. Por fim, propõe-se uma estratégia para estimar os parâmetros termodinâmicos do NRTL para as espécies presentes nos sistemas estudados e prever o comportamento dos mesmos e de sistemas similares. Este modelo matemático proposto tem como objetivo complementar a técnica experimental adotada para construção das curvas de equilíbrio, sendo bem sucedido

Otimização do problema de localização de instalações aplicado ao comércio e distribuição de combustíveis

Um dos problemas mais relevantes em organizações de grande porte é a escolha de locais para instalação de plantas industriais, centros de distribuição ou mesmo pontos comerciais. Esse problema logístico é uma decisão estratégica que pode causar um impacto significativo no custo total do produto comercializado. Existem na literatura diversos trabalhos que abordam esse problema. Assim, o objetivo desse trabalho é analisar o problema da localização de instalações proposto por diferentes autores e definir um modelo que seja o mais adequado possível ao mercado de distribuição de combustíveis no Brasil. Para isso, foi realizada uma análise do fluxo de refino e distribuição praticado neste segmento e da formação do respectivo custo de transporte. Foram consideradas restrições como capacidade de estoque, gama de produtos ofertados e níveis da hierarquia de distribuição. A partir dessa análise, foi definido um modelo matemático aplicado à redução dos custos de frete considerando-se a carga tributária. O modelo matemático foi implementado, em linguagem C, e permite simular o problema. Foram aplicadas técnicas de computação paralela visando reduzir o tempo de execução do algoritmo. Os resultados obtidos com o modelo Single Uncapacited Facility Location Problem (SUFLP) simulado nas duas versões do programa, sequencial e paralela, demonstram ganhos de até 5% em economia de custos e redução do tempo de execução em mais de 50%.