As fases fenológicas e produtividade de soja para a região de Londrina, norte do estado do Paraná, foram estimadas utilizando um modelo de desenvolvimento e crescimento de soja denominado SOYGRO V5.42. O modelo simula o desenvolvimento fenológico (duração dos estádios vegetativos e reprodutivos) da cultura da soja e balanço hídrico, em função dos elementos meteorológicos diários, do tipo de solo e da cultivar. Os testes aplicados para avaliar a sensibilidade do modelo, mostram que o modelo é sensível à densidade de plantio, e às variáveis de solo principalmente à variação de umidade. O modelo, porém não mostrou sensibilidade em relação ao albedo e coeficiente de evaporação do solo no estágio 1. O modelo requer dados meteorológicos diários, informações sobre os parâmetros de solo, principalmente sobre as condições hídricas e dados genéticos da cultivar. Permite ainda entrar com manejo de irrigação e aplicação de fertilizantes. Os resultados apresentados pelo modelo, em geral, mostraram uma tendência em superestimar as produtividades, e aumentar os períodos semeadura-emergência, semeadura-florescimento e semeadura-maturidade fisiológica / not available

Quantificação de incerteza e risco tem aplicação direta no projeto de limites de cava e na programação de produção. Pequenas variações nas condições de contorno de um projeto podem ter impacto significativo sobre o retorno final do mesmo na ordem de milhões de dólares. Preço de mercado do bem mineral, o custo de capital, taxas de atratividade, etc, são fatores usualmente testados para analisar a viabilidade de um empreendimento, porém, raramente é considerada a incerteza relacionada ao atributo geológico em questão. O propósito de um planejamento de lavra tem sido prover subsidio para o engenheiro de minas decidir sobre a capacidade de minerar determinadas unidades de lavra do depósito a partir de teores estimados. Salienta-se porém que existe, a partir dos dados amostrais, incertezas a respeito do modelo geológico que devem ser consideradas ao se projetar uma cava e desenvolver a lavra. A simulação geoestatistica tem a capacidade de produzir múltiplos modelos equiprováveis, os quais podem ser avaliados independentementecomo possíveis cenários do depósito mineral. Simulação condicional, ao contrário de técnicas de interpolação e estimativa, provê respostas para questões associadas a risco devido à variações nos teores do modelo geológico. Ao gerar múltiplos cenários tem-se acesso à probabilidade e conseqüentemente às variações de retorno financeiro e rotas de extração de minério de um projeto. o presente trabalho tem como objetivo a investigação de novas metodologias que contribuam para a construção de cenários de planejamento de lavra e avaliação do impacto provocado nestes pela incerteza fornecida a partir de modelos simulados. As respostas buscadas dentro da abordagem sugerida por esse trabalho, compreendem a definição de uma metodologia para controle e planejamento de lavra a médio e longo prazo por incorporação da flutuabilidade local associada ao minério, avaliando a sensibilidade do retorno financeiro e o traçado de lavra em relação ao método de geração do modelo geológico de teores. Para solucionar o problema em questão, sugere-se a geração de modelos estocásticos e a alimentação de múltiplos modelos selecionados por critérios considerados relevantes para a delimitação do espaço de incerteza. A aplicação de funções de transferência tais como mudança de suporte e seqüenciamento de produção sobre os modelos simulados é a chave para a obtenção de respostas a respeito do impacto sobre a lucratividade de um projeto. Ao alimentar à essas funções modelos equiprováveis que preservam as características estatísticas e a conectividade espacial dos dados amostrais mas que invariavelmente possuem diferentes distribuições de teores tem-se a dimensão em termos econômicos do risco associado à incerteza geológica. Foi confrontado o retorno financeiro produzido por modelos simulados e verificou-se para o caso especifico que os métodos de simulação geoestatistica empregados não produziram diferenças significativas quando comparados casos selecionados segundo os mesmos critérios. Porém, entre cenários extremos gerados pelo mesmo algoritmo de simulação foram verificadas desigualdades relevantes. Verificou-se também a transferência de minério entre diferentes classes de qualidade ao proceder-se com variações nas dimensões dos blocos de lavra.

A simulação paralela de eventos é uma área da computação que congrega grande volume de pesquisas, pela importância em facilitar o estudo de novas soluções nas mais diferentes áreas da ciência e tecnologia, sem a necessidade da construção de onerosos protótipos. Diversos protocolos de simulação paralela podem ser encontrados, divididos em dois grandes grupos de acordo com o algoritmo empregado para a execução em ordem dos eventos: os conservadores e os otimistas; contudo, ambos os grupos utilizam trocas de mensagens para a sincronização e comunicação. Neste trabalho, foi desenvolvido um novo protocolo de simulação paralela, fazendo uso de memória compartilhada, o qual foi implementado e testado sobre um ambiente de estações de trabalho, realizando, assim, simulação paralela com uso de memória compartilhada distribuída. O protocolo foi desenvolvido tendo como base de funcionamento os protocolos conservadores; utilizou diversas características dos mesmos, mas introduziu várias mudanças em seu funcionamento. Sua execução assemelha-se às dos protocolos de execução síncrona, utilizando conceitos como o lookahead e janelas de tempo para execução de eventos. A principal mudança que o novo protocolo sofreu foi proporcionada pelo acesso remoto à memória de um LP por outro, produzindo diversas outras nas funções relativas à sincronização dos processos, como o avanço local da simulação e o agendamento de novos eventos oriundos de outro LP. Um ganho adicional obtido foi a fácil resolução do deadlock, um dos grandes problemas dos protocolos conservadores de simulação paralela. A construção de uma interface de comunicação eficiente com uso de memória compartilhada é o principal enfoque do protocolo, sendo, ao final da execução de uma simulação, disponibilizado o tempo de simulação e o tempo de processamento ocioso (quantia utilizada em comunicação e sincronização). Além de uma implementação facilitada, propiciada pelo uso de memória compartilhada ao invés de trocas de mensagens, o protocolo oferece a possibilidade de melhor ocupar o tempo ocioso dos processadores, originado por esperas cada vez que um LP chega a uma barreira de sincronização. Em nenhum momento as modificações efetuadas infringiram o princípio operacional dos protocolos conservadores, que é não possibilitar a ocorrência de erros de causalidade local. O novo protocolo de simulação foi implementado e testado sobre um ambiente multicomputador de memória distribuída, e seus resultados foram comparados com dois outros simuladores, os quais adotaram as mesmas estratégias, com idênticas ferramentas e testados em um mesmo ambiente de execução. Um simulador implementado não utilizou paralelismo, tendo seus resultados sido utilizados como base para medir o speedup e a eficiência do novo protocolo. O outro simulador implementado utilizou um protocolo conservador tradicional, descrito na literatura, realizando as funções de comunicação e sincronização através de trocas de mensagens; serviu para uma comparação direta do desempenho do novo protocolo proposto, cujos resultados foram comparados e analisados.

Apresentar um modelo para simular um sistema de armazenamento de calor no solo em estufas para plasticultura é o objetivo do presente trabalho. O sistema consiste num feixe de tubos enterrados no solo. A convecção forçada de ar no seu interior realiza a troca térmica necessária para manter as estufas sob faixas desejadas de temperatura. O objetivo do modelo é investigar os efeitos no calor armazenado e a influência das variáveis, tais como diâmetro, comprimento, espaçamento entre os tubos e a velocidade de ar no canal provocam no sistema. O solo é tratado como um meio difusivo e avalia-se a contribuição do termo de condensação e evaporação da água contida no ar em escoamento nos tubos. A equação da energia é resolvida para o solo e para o ar. Os tubos de seção transversal circular são modelados como tubos de seção transversal quadrada com o objetivo de que as simulações possam ser processadas em coordenadas cartesianas. O programa resolve situações tridimensionais, transientes e emprega o Método dos Volumes Finitos para integrar as equações diferenciais governantes. O modelo original é baseado no modelo de Gauthier et al., 1997, tendo sido os resultados do mesmo foram usados para a validação do presente estudo. Um circuito de água quente é também projetado e apresentado para o aquecimento das estufas. A água circula através de mangueiras sobre o solo e é aquecida por um sistema de queimadores a gás liqüefeito de petróleo ou óleo combustível, transferindo assim calor para o interior da mesma. O projeto de aquecimento foi realizado através de um programa de parceria entre a Ufrgs, Sebrae, Fapergs e a Agropecuária Clarice.

O presente trabalho apresenta uma comparação entre resultados computacionais, obtidos através do "software" EnergyPlus e experimentais do comportamento térmico de um ambiente condicionado e não condicionado. Para tanto, monitoraram-se os dados climáticos de radiação, velocidade do vento e temperatura, no período de 11 a 20 de janeiro de 2002 e produziu-se um arquivo climático. Simultaneamente, fez-se a aquisição das temperaturas de uma sala-teste, localizada no terceiro pavimento de um prédio na cidade de Porto Alegre, bem como das salas adjacentes. As temperaturas do ar de insuflamento e de retorno dos condicionadores de ar, localizados na sala-teste, foram medidas durante o dia, em seis dias do período de monitoramento. Mediu-se também a velocidade do ar de retorno e determinou-se a potência sensível de refrigeração. Os ganhos de calor interno da sala também foram medidos e utilizaramse as formulações apresentadas pela ASHRAE, 2001, para determiná-los em relação às pessoas e à infiltração. Tais dados foram declarados ao programa como variáveis de entrada. As simulações do comportamento térmico da sala-teste foram implementadas informando-se ao programa a temperatura das salas ou os coeficientes de uma equação representativa das mesmas. Por considerar que a primeira representava melhor as condições térmicas das salas contíguas, utilizou-se esta modalidade para análise As simulações foram conduzidas, alterando-se opções fornecidas pelo programa: modelo de céu isotrópico e anisotrópico, coeficiente de convecção simples e detalhado. Os resultados da carga térmica e temperatura ambiente da sala-teste obtidos nas simulações foram comparados com os dados experimentais do período de monitoramento. A melhor concordância foi obtida para o modelo anisotrópico, coeficiente de convecção detalhado. Constatou-se uma grande discrepância entre as cargas térmicas, para o modelo de convecção simples. Assim, conclui-se que o "software" EnergyPlus representa bem o comportamento térmico de uma edificação "termicamente pesada" para coeficiente de convecção detalhado, necessitando pesquisa para as demais edificações.

Combinando métodos computacionais da eletrostática e dinâmica molecular, este trabalho teve como objetivo descrever processos de solvatação em sistemas quimicamente importantes. Foram obtidas propriedades termodinâmicas necessárias para o entendimento do processo de solvatação. Estão descritos e testados modelos para descrever a interação soluto-solvente, possibilitando, assim, aprimorar a descrição físico-química dos processos de solvatação. Utilizaram-se programas desenvolvidos em nosso grupo e programas comerciais que permitem os cálculos de dinâmica molecular e química quântica. Uma nova abordagem para o cálculo de energia livre de solvatação foi desenvolvida proporcionando a obtenção acurada e eficiente dessa propriedade, dentro do enfoque da dinâmica molecular. Nessa nova abordagem, novas metodologias para a geração de cavidades moleculares foram propostas e avaliadas. As energias livres de solvatação obtidas estão em boa concordância com os valores experimentais.

O uso da mecânica de fluidos computacional no estudo de processos envolvendo o escoamento de fluidos poliméricos está cada vez mais presente nas indústrias de transformação de polímeros. Um código computacional voltado a esta função, para que possa ser aplicado com sucesso, deve levar a predições mais próximas possível da realidade (modelagem), de uma forma relativamente rápida e eficiente (simulação). Em relação à etapa de modelagem, o ponto chave é a seleção de uma equação constitutiva que represente bem as características reológicas do fluido, dentre as diversas opções existentes. Para a etapa de simulação, ou seja, a resolução numérica das equações do modelo, existem diversas metodologias encontradas na literatura, cada qual com suas vantagens e desvantagens. Neste tópico se enquadra o trabalho em questão, que propõe uma nova metodologia para a resolução das equações governantes do escoamento de fluidos viscoelásticos. Esta se baseia no método dos volumes finitos, usando o arranjo co-localizado para as variáveis do problema, e na utilização de aproximações de alta ordem para os fluxos médios lineares e não-lineares e para outros termos não lineares que surgem da discretização das equações constitutivas. Nesta metodologia, trabalha-se com os valores médios das variáveis nos volumes durante todo o processo de resolução, sendo que os valores pontuais são obtidos ao final do procedimento via deconvolução. A solução do sistema de equações não lineares, resultante da discretização das equações, é feita de forma simultânea, usando o método de Newton São mostrados então, resultados da aplicação da metodologia proposta em problemas envolvendo escoamentos de fluidos newtonianos e fluidos viscoelásticos. Para descrever o comportamento reológico destes últimos, são usadas duas equações constitutivas, que são o modelo de Oldroyd-B e o modelo de Phan-Thien-Tanner Simplificado. Por estes resultados pode-se ver que a metodologia é muito promissora, apresentando algumas vantagens frente às metodologias convencionais em volumes finitos. A implementação atual da metodologia desenvolvida está restrita a malhas uniformes e, consequentemente, soluções para problemas com geometrias complexas, que necessitam de refinamento localizado da malha, foram obtidas somente para baixos números de Weissenberg, devido a limitação do custo computacional. Esta restrição pode ser contornada, tornando o seu uso competitivo.

Cada vez mais as siderúrgicas focadas na produção de aços especiais fazem uso de dispositivos eletromagnéticos que agem de forma benéfica no escoamento do aço no interior do molde e também nas trocas térmicas entre o aço líquido e o molde de cobre. Porém, o uso desses dispositivos pode acarretar perturbações indesejáveis na região de interface entre o metal e a camada protetora de escória no topo do molde. Essas perturbações na interface podem proporcionar arraste e emulsificação de escória. O aprisionamento de escória no molde pode dar origem a problemas como defeitos superficiais e internos nos tarugos solidificados. Baseado em práticas correntes nas plantas siderúrgicas este trabalho pretende verificar as condições que propiciam o arraste de escória na região superior do molde, nos casos em que se aplica agitação eletromagnética. Este estudo foi realizado em duas etapas, sendo que a primeira etapa destina-se o estudo e desenvolvimento de um modelo, fazendo uso de uma ferramenta numérica, que represente o funcionamento de agitadores eletromagnéticos usados na siderurgia. Dando seqüência ao estudo, foi então desenvolvido um modelo fluidodinâmico do molde do lingotamento contínuo, dando enfoque ao comportamento da interface metal – escória. O efeito do agitador sob o escoamento foi acoplado ao modelo fluidodinâmico, através de resultados do modelo eletromagnético calculados no modelo anterior. O objetivo desta etapa do trabalho é estudar as características de deformação do menisco e tentar avaliar a ocorrência de arraste de escória quando da variação de alguns parâmetros, como por exemplo, intensidade de agitação e espessura de camada líquida de escória. Foram analisados escoamentos em moldes circulares com duas diferentes de seções transversais. Nestes escoamentos foram testadas três intensidades de agitação eletromagnéticas para cada um deles. Outro parâmetro testado foi a espessura da camada de escória no menisco. Este estudo demonstrou que a agitação eletromagnética excessiva pode acarretar arraste de escória, como indicado na literatura. Foi possível também verificar o tipo de deformação que ocorre na interface e as características que podem proporcionar circunstâncias críticas para o arraste. Além da intensidade de agitação, os resultados também sugerem que, quanto menor profundidade de submersão da válvula, maior a probabilidade de ocorrência de arraste. / Each time more the steelmaking plants, focused on the special steels productions, have used electromagnetic devices which act well not only on the steel flow inside the mold but also on the thermal changes between the liquid steel and the copper mold. However, the use of these devices may cause undesirable instability in the interface steel – slag in the top of the mold. These instabilities in the interface may cause entrapment and emulsification of the slag. The entrapment of the slag in the mold may originate problems like superficial and internal defects in the solidified billet. Based on current practices in the steelmaking plants, this paper aims to verify the conditions that promote the slag entrapment in the upper region of the mold, in the cases in which the electromagnetic stirring is applied. This study was developed in two stages. The first stage aimed to study and to construct a numerical model which represents the electromagnetic stirrer work. As a sequence of the study, a fluid-dynamic model to the flow in the continuous casting mold was developed, with focus on the behavior of the steel – slag interface. The electromagnetic stirrer effects, under the steel flow in the mold, were coupled with the fluid-dynamic model, through the electromagnetic model results, previously calculated. The aim of this work is to verify the possibilities of slag entrapment occurrence when there is variation of the parameters like the electric current intensity (electromagnetic stirring), slag thickness, slag viscosity and submersion nozzle depth. In this work the fluid flow was analyzed in a circular mold with two different sections (180mm and 230mm). Each section mold was tested with a type of mold powder. In this flow three different stirrer intensities were tested. This study shows that the excessive electromagnetic stirring may cause slag entrapment, as indicated in the literature. It was also possible to verify the kind of deformation that happens on the steel-slag interface, and the characteristics that bring critical circumstances to the slag entrapment. Besides the stirring intensity, the results also suggest that the smaller the submersion nozzle depth the greater the probability of the slag entrapment occurrence.

Nesse trabalho é apresentada uma solução analítica para a equação de transporte bidimensional em um domínio retangular considerando o espalhamento Compton, utilizando o método LTSN e a aproximação PN na variável angular. Esse procedimento leva a uma formulação para a taxa de dose total devido µa radiação gama. Nessa derivação, a energia depositada pelo elétron livre é considerada. Para alcançarmos esse objetivo dois problemas são resolvidos: o primeiro para o transporte de fótons, assumindo como seção de choque diferencial o núcleo de espalhamento de Klein-Nishina bem como modelo multigrupo em energia; e o segundo para o transporte de elétrons livres, considerando a seção de choque de espalhamento diferencial de Rutherford. Simulações numéricas e validações com os resultados obtidos pelo método de Monte Carlo são apresentadas para ambos os problemas. / In this work we present the analytical solution for the two-dimensional transport equation in a rectangle considering Compton scattering, using the LTSN method and the PN approach for the angular variable. This procedure leads to a formulation for the total dosis due the gamma radiation. In this derivation, the energy deposited by the free electron is taken into account. To reach this objective we solve two model problems: the ¯rst one for the transport of photons, assuming the Klein-Nishina scattering kernel as the scattering di®erential cross section as well as the multigroup model in the energy variable; and the second one the transport of free electrons considering the screened Rutherford scattering kernel. Numerical simulations and validations with Monte Carlo results are reported for both problems.

Este trabalho analisa as etapas de fabricação, por Forjamento em Matriz Aberta, de um eixo vazado em aço ABNT 4140 (~DIN número 1.7225, grau 42CrMo4), com auxílio da simulação numérica computacional utilizando o software PEP/LARSTRAN. Com base nos resultados de simulação obtidos foram feitas alterações no modelo e adaptações nas simulações numéricas, visando a obtenção de um processo confiável. A estratégia final desenvolvida por simulação se tratou de um forjamento de todo o perímetro do diâmetro externo da pré-forma e subsequente avanço longitudinal do ferramental para repetir o mesmo ciclo de forjamento. Esta estratégia de forjamento foi utilizada para realizar um forjamento de um eixo piloto, no qual ocorreu o forjamento do eixo a 1200 °C, a uma velocidade de 80 mm/s, com recalque de 20 mm por batida. O resultado obtido foi um eixo com alongamento Δl = 77% sem deformação significativa do diâmetro interno. Os resultados obtidos estão relativamente de acordo com a simulação numérica computacional. / This paper analyzes hot forging steps of a hollow shaft by the open die forging process. The study was aided by numerical simulation using the software PEP/LARSTRAN. The material used was a steel ABNT 4140 (~DIN number 1.7225, grade 42CrMo4). The proposed process was simulated and optimized sequentially. That means several changes in the numerical model and simulation in order to obtain a feasible process. The final strategy developed by simulation was forge initially the entire perimeter of the preform with its outer diameter and subsequently forward longitudinally with the same tools to repeat the complete forging cycle. The tool set was a flat upper die and a V-profile lower die. A pilot shaft was forged to test the final forging strategy developed by simulation. The physical experiment was carried out with the workpiece at 1200 °C. The tool speed was 80 mm/s and the stroke was 20 mm per blow. It was obtained 77% final elongation of the part without significant deformation of internal diameter. The experimental results are in relatively good agreement with the numerical simulations.