Estudo do galgamento em estruturas marítimas verticais utilizando o modelo numérico AMAZON

O galgamento oceânico é um fenómeno associado à interação de ondas com estruturas marítimas. Nesta dissertação estudou-se o galgamento em estruturas marítimas verticais através da aplicação do modelo numérico AMAZON, que é um modelo bidimensional, baseado na resolução das equações não lineares de águas pouco profundas. Para validar a sua aplicação a estruturas marítimas verticais impermeáveis, os resultados de caudal médio galgado obtidos com o modelo AMAZON foram comparados com os resultados de ensaios em modelo físico bidimensional realizados anteriormente e disponíveis na base de dados do projeto europeu CLASH. Os resultados foram ainda comparados com os das fórmulas empíricas e semi-empíricas de Van der Meer & Bruce de 2014, de Goda de 2008 e de Hedges & Reis de 2015. Antes da sua validação, o modelo foi calibrado utilizando um dos ensaios do CLASH. Para tal, realizou-se uma análise de sensibilidade a alguns parâmetros do modelo, concretamente da dimensão mínima da malha, da escala de aplicação do AMAZOM (à escala reduzida ou à do protótipo), da forma de modelação da parede vertical, já que, dado este modelo ser integrado na vertical, não é possível considerar paredes verticais sem utilizar um artifício de cálculo, dos valores de alguns dos parâmetros de entrada do modelo (profundidade mínima molhada e a profundidade mínima de fricção) e da dimensão do domínio de cálculo. Os resultados da calibração do AMAZON mostram que, por causa da precisão do modelo para os menores caudais médios de galgamento, se deve aplicar o modelo à escala do protótipo. Para estruturas verticais, a aproximação que o modelo faz, por defeito, para a secção vertical da estrutura através de um talude fictício a 10:1 apresenta melhores resultados que a utilização de um talude a 20:1. Contudo, ambos requerem a adoção de uma dimensão da malha bastante reduzida nesta zona e na sua vizinhança (0.03 m). No caso em estudo, os dois parâmetros de entrada do modelo não têm impacto significativo nos resultados, ao contrário do domínio de cálculo que, tal como recomendado na literatura de referência, se deve iniciar (fronteira de entrada) a aproximadamente um comprimento de onda do pé da estrutura. Os resultados obtidos mostram que, em geral, o AMAZON subestima o valor medido do caudal médio galgado com erros relativos variáveis entre -91% e -49%; apenas para um caso o modelo sobrestima o caudal médio, com um erro relativo de 43%. Quando comparado com as formulações empíricas e semi-empíricas, a fórmula de Van der Meer & Bruce apresenta melhores aproximações aos caudais médios medidos que constam na base de dados do CLASH, seguida da fórmula de Goda e dos resultados do modelo AMAZON, embora estes dois últimos apresentem valores bastante próximos. Assim, uma vez que o modelo está em constante desenvolvimento, espera-se que no futuro seja uma ferramenta mais poderosa. Contudo, enquanto não se fizerem melhorias no modelo e/ou não se estender a outros ensaios a análise feita no âmbito desta dissertação, os resultados sugerem que não há vantagem na sua utilização face às formulações empíricas e semi-empíricas disponíveis, tendo também em conta o significativo tempo computacional exigido para executar o AMAZON com a dimensão da malha escolhida.

Estimativa da área foliar da mangueira (Mangifera indica L.) cvs. Tommy Atkins e Haden, utilizando dimensões lineares

Os trabalhos existentes sobre a estimativa da área foliar da mangueira são raros e de resultados pouco convergentes. O objetivo deste trabalho foi determinar as equações que melhor descrevem a relação entre as dimensões lineares (comprimento e largura máxima) e a área da folha da mangueira, cultivares Tommy Atkins e Haden, possibilitando a estimativa rápida e não-destrutiva da área foliar (AF) utilizando apenas o comprimento (C), a largura máxima (L), ou ambos. A área foliar da mangueira pode ser estimada multiplicando-se o parto do comprimento pela largura, pelo fator 0,74 (para a cultivar Tommy Atkins) ou 0,78 (para a cultivar Haden), bem como pelas equações: AF= 4,96349 C - 33,42976 (R²=0,86); AF= 17,02964 L - 18,88065 (R²=0,85) e AF= 0,73499 (CxL) + 0,59459 (R²=0,92) para a cultivar Tommy Atkins ou AF= 5,35282C - 33,17061(R²=0,88); AF= 19,09951L - 24,61777(R²=0,89) e AF= 0,76015(CxL) + 0,43257(R²=0,99) para a cultivar Haden.

Equações alométricas para estimativa de biomassa e carbono em árvores de reflorestamentos de restauração

Este estudo teve como objetivos desenvolver e ajustar modelos alométricos para estimativa da biomassa seca e do carbono total em árvores plantadas em áreas de restauração florestal. Os dados dos ajustes são provenientes de 107 árvores de 44 espécies plantadas na região do Médio Vale do Paranapanema, SP, localizado nos Biomas Mata Atlântica e Cerrado. A biomassa seca e a massa de carbono foram obtidas por meio de amostragem destrutiva da parte aérea e subterrânea das árvores. Para o teste e ajuste dos modelos foi realizada a estratificação do conjunto inicial dos dados em ritmos de crescimento das espécies amostradas. Os ajustes foram feitos usando-se oito modelos lineares de cada variável dependente e dois obtidos pelo processo Stepwise-Forward. Os melhores modelos para estimativa da biomassa seca e carbono orgânico apresentaram coeficientes de determinação ajustado acima de 0,95 e erros-padrão percentuais abaixo de 32%. Os modelos baseados nos ritmos de crescimento das espécies mostraram melhores resultados estatísticos, atingindo R² = 0,985 e Syx% = 16,15 para biomassa seca das espécies de crescimento lento. Modelos provenientes do procedimento Stepwise geraram as melhores equações para as estimativas de biomassa seca e carbono total, e a estratificação dos dados dos diferentes ritmos de crescimento das espécies amostradas mostrou-se que foi adequada para melhorar o desempenho dos modelos.

Infiltração de água no solo em áreas cultivadas com cana-de-açúcar sob diferentes sistemas de colheita e modelos de ajustes de equações de infiltração

As técnicas de colheita da cultura de cana nem sempre seguem preceitos conservacionistas, fato que pode influenciar atributos físicos do solo com consequente alteração da velocidade de infiltração de água no solo. O trabalho realizado teve como objetivo avaliar o efeito de diferentes sistemas de colheita (manual e mecanizada) da cana-de-açúcar (crua e queimada), sobre as características de infiltração e verificar a adequação de modelos matemáticos para a estimativa da taxa de infiltração de água no solo. Foi utilizado um simulador de chuva portátil calibrado para aplicar uma precipitação de 60 mm h-1. Os testes de infiltração de água no solo foram realizados em três sistemas de colheita: colheita manual de cana queimada, colheita mecanizada de cana crua e colheita mecanizada de cana queimada. A qualidade do ajuste dos modelos foi avaliada por meio de regressões não lineares entre os valores estimados e os valores médios observados em cada tratamento estudado. Nos tratamentos onde foi empregada a colheita mecânica, observou-se uma redução da taxa de infiltração final em relação à área sob colheita manual de cana queimada. A equação de Horton foi a mais adequada para a estimativa da taxa de infiltração de água no solo estudado.

Estimativa da área foliar de Panicum maximum usando dimensões lineares do limbo foliar

Com o objetivo de obter uma equação que, por meio de parâmetros lineares dimensionais das folhas, permitisse a estimativa da área foliar de Panicum maximum, estudaram-se relações entre a área foliar real (Sf) e os parâmetros dimensionais do limbo foliar, como o comprimento ao longo da nervura principal (C) e a largura máxima (L), perpendicular à nervura principal. As equações lineares simples, exponenciais e geométricas obtidas podem ser usadas para estimação da área foliar. Do ponto de vista prático, sugere-se optar pela equação envolvendo apenas o produto C x L. Desse modo, a estimativa da área foliar de P. maximum pode ser feita pela fórmula Sf = 0,6058 x (C x L), que equivale a tomar 60,58% do produto entre o comprimento ao longo da nervura principal e a largura máxima, com coeficiente de determinação de 0,8586.

Estimativa da área foliar de Cissampelos glaberrima usando dimensões lineares do limbo foliar

Com o objetivo de obter uma equação matemática que, através de parâmetros lineares dimensionais das folhas, permitisse a estimativa da área foliar de Cissampelos glaberrima, estudaram-se relações entre a área foliar real (Sf) e os parâmetros dimensionais do limbo foliar, como o comprimento ao longo da nervura principal (C) e a largura máxima (L) perpendicular à nervura principal. As equações lineares simples, exponenciais e geométricas obtidas podem ser usadas para estimação da área foliar da falsa parreira-brava. Do ponto de vista prático, sugere-se optar pela equação linear simples envolvendo o produto C x L, usando-se a equação de regressão Sf = 0,7878 x (C x L), que equivale a tomar 78,78% do produto entre o comprimento ao longo da nervura principal e a largura máxima, com coeficiente de correlação de 0,9307.

Estimativa da área foliar de Typha latifolia usando dimensões lineares do limbo foliar

Com o objetivo de obter uma equação que, através de parâmetros lineares dimensionais das folhas, permita a estimativa da área foliar de Typha latifolia, estudaram-se relações entre a área foliar real (Sf) e parâmetros dimensionais do limbo foliar, como o comprimento ao longo da nervura principal (C) e a largura máxima (L), perpendicular à nervura principal. As equações lineares simples, exponenciais e geométricas obtidas podem ser usadas para estimação da área foliar da taboa. Do ponto de vista prático, sugere-se optar pela equação linear simples que envolve o produto C x L, usando-se a equação de regressão Sf = 0,9651 x (C x L), que equivale a tomar 96,51% do produto entre o comprimento ao longo da nervura principal e a largura máxima, com um coeficiente de determinação de 0,9411.

Estimativa da área foliar de Tridax procumbens usando dimensões lineares do limbo foliar

Com o objetivo de obter uma equação que, através de parâmetros lineares dimensionais das folhas, permita a estimativa da área foliar de Tridax procumbens, estudaram-se relações entre a área foliar real (Sf) e os parâmetros dimensionais do limbo foliar, como o comprimento ao longo da nervura principal (C) e a largura máxima (L), perpendicular à nervura principal. As equações lineares simples, exponenciais e geométricas obtidas podem ser usadas para estimação da área foliar da erva-de-touro. Do ponto de vista prático, sugere-se optar pela equação linear simples envolvendo o produto C x L, usando-se a equação de regressão Sf = 0,6008 x (C x L), que equivale a tomar 60,08% do produto entre o comprimento ao longo da nervura principal e a largura máxima, com um coeficiente de determinação de 0,8731.

Estimativa da área foliar de Brachiaria plantaginea usando dimensões lineares do limbo foliar

Com o objetivo de obter uma equação que, por meio de parâmetros lineares dimensionais das folhas, permita a estimativa da área foliar de Brachiaria plantaginea, estudaram-se relações entre a área foliar real (Sf) e os parâmetros dimensionais do limbo foliar, como o comprimento ao longo da nervura principal (C) e a largura máxima (L), perpendicular à nervura principal. As equações lineares simples, exponenciais e geométricas obtidas podem ser usadas para estimação da área foliar do capim-marmelada. Do ponto de vista prático, deve-se optar pela equação linear simples, envolvendo o produto C x L, usando-se a equação de regressão Sf = 0,7338 x (C x L), o que equivale a tomar 73,38% do produto entre o comprimento ao longo da nervura principal e a largura máxima, com um coeficiente de determinação de 0,8754.

Estimativa da área foliar de cinco espécies do gênero Amaranthus usando dimensões lineares do limbo foliar

Este trabalho foi realizado com o objetivo de obter equações que estimem a área foliar de espécies de plantas daninhas do gênero Amaranthus utilizando as dimensões lineares do limbo foliar, bem como comparar as equações obtidas para plantas que cresceram em casa de vegetação e em campo. Para isso, plantas de A. deflexus, A. hybridus, A. retroflexus, A. spinosus e A. viridis foram cultivadas em casa de vegetação ou coletadas em infestações naturais de áreas agrícolas e não-agrícolas do município de Piracicaba-SP. Para cada condição foram analisados 100 limbos foliares, com relação a comprimento ao longo da nervura principal (C), largura máxima (L) e área foliar real (Ar). Os conjuntos de dados foram ajustados à equação linear que passa pela origem Ar = a.(C.L) e, para todas as variáveis, comparados por meio de intervalos de confiança, a 5% de significância. Concluiu-se que: há espécies de plantas daninhas do gênero Amaranthus com tamanhos diferenciais de folhas, sendo A. deflexus aquela que, em média, possui as menores dimensões; a equação linear passando pela origem (Ar = a.(C.L)) foi adequada para ajustar a relação entre as medidas lineares do limbo e a área foliar real; e plantas que cresceram em casa de vegetação e em campo podem apresentar diferenças quanto ao tamanho de folhas ou quanto ao parâmetro a da equação de ajuste da reta.

Paralelização de métodos de resolução de sistemas lineares esparsos com o DECK em um Cluster de PCs

O objetivo desta dissertação é a paralelização e a avaliação do desempenho de alguns métodos de resolução de sistemas lineares esparsos. O DECK foi utilizado para implementação dos métodos em um cluster de PCs. A presente pesquisa é motivada pela vasta utilização de Sistemas de Equações Lineares em várias áreas científicas, especialmente, na modelagem de fenômenos físicos através de Equações Diferenciais Parciais (EDPs). Nessa área, têm sido desenvolvidas pesquisas pelo GMC-PAD – Grupo de Matemática da Computação e Processamento de Alto Desempenho da UFRGS, para as quais esse trabalho vem contribuindo. Outro fator de motivação para a realização dessa pesquisa é a disponibilidade de um cluster de PCs no Instituto de Informática e do ambiente de programação paralela DECK – Distributed Execution and Communication Kernel. O DECK possibilita a programação em ambientes paralelos com memória distribuída e/ou compartilhada. Ele está sendo desenvolvido pelo grupo de pesquisas GPPD – Grupo de Processamento Paralelo e Distribuído e com a paralelização dos métodos, nesse ambiente, objetiva-se também validar seu funcionamento e avaliar seu potencial e seu desempenho. Os sistemas lineares originados pela discretização de EDPs têm, em geral, como características a esparsidade e a numerosa quantidade de incógnitas. Devido ao porte dos sistemas, para a resolução é necessária grande quantidade de memória e velocidade de processamento, característicos de computações de alto desempenho. Dois métodos de resolução foram estudados e paralelizados, um da classe dos métodos diretos, o Algoritmo de Thomas e outro da classe dos iterativos, o Gradiente Conjugado. A forma de paralelizar um método é completamente diferente do outro. Isso porque o método iterativo é formado por operações básicas de álgebra linear, e o método direto é formado por operações elementares entre linhas e colunas da matriz dos coeficientes do sistema linear. Isso permitiu a investigação e experimentação de formas distintas de paralelismo. Do método do Gradiente Conjugado, foram feitas a versão sem précondicionamento e versões pré-condicionadas com o pré-condicionador Diagonal e com o pré-condicionador Polinomial. Do Algoritmo de Thomas, devido a sua formulação, somente a versão básica foi feita. Após a paralelização dos métodos de resolução, avaliou-se o desempenho dos algoritmos paralelos no cluster, através da realização de medidas do tempo de execução e foram calculados o speedup e a eficiência. As medidas empíricas foram realizadas com variações na ordem dos sistemas resolvidos e no número de nodos utilizados do cluster. Essa avaliação também envolveu a comparação entre as complexidades dos algoritmos seqüenciais e a complexidade dos algoritmos paralelos dos métodos. Esta pesquisa demonstra o desempenho de métodos de resolução de sistemas lineares esparsos em um ambiente de alto desempenho, bem como as potencialidades do DECK. Aplicações que envolvam a resolução desses sistemas podem ser realizadas no cluster, a partir do que já foi desenvolvido, bem como, a investigação de précondicionadores, comparação do desempenho com outros métodos de resolução e paralelização dos métodos com outras ferramentas possibilitando uma melhor avaliação do DECK.

Reconhecimento automático de locutor utilizando medidas de invariantes dinâmicas não-lineares

As técnicas utilizadas em sistemas de reconhecimento automático de locutor (RAL) objetivam identificar uma pessoa através de sua voz, utilizando recursos computacionais. Isso é feito a partir de um modelamento para o processo de produção da voz. A modelagem detalhada desse processo deve levar em consideração a variação temporal da forma do trato vocal, as ressonâncias associadas à sua fisiologia, perdas devidas ao atrito viscoso nas paredes internas do trato vocal, suavidade dessas paredes internas, radiação do som nos lábios, acoplamento nasal, flexibilidade associada à vibração das cordas vocais, etc. Alguns desses fatores são modelados por um sistema que combina uma fonte de excitação periódica e outra de ruído branco, aplicadas a um filtro digital variante no tempo. Entretanto, outros fatores são desconsiderados nesse modelamento, pela simples dificuldade ou até impossibilidade de descrevê-los em termos de combinações de sinais, filtros digitais, ou equações diferenciais. Por outro lado, a Teoria dos Sistemas Dinâmicos Não-Lineares ou Teoria do Caos oferece técnicas para a análise de sinais onde não se sabe, ou não é conhecido, o modelo detalhado do mecanismo de produção desses sinais. A análise através dessa teoria procura avaliar a dinâmica do sinal e, assumindo-se que tais amostras provêm de um sistema dinâmico não-linear, medidas qualitativas podem ser obtidas desse sistema. Essas medidas não fornecem informações precisas quanto ao modelamento do processo de produção do sinal avaliado, isto é, o modelo analítico é ainda inacessível. Entretanto, pode-se aferir a respeito de suaO problema analisado ao longo deste trabalho trata da busca de novos métodos para extrair informações úteis a respeito do locutor que produziu um determinado sinal de voz. Com isso, espera-se conceber sistemas que realizem a tarefa de reconhecer um pessoa automaticamente através de sua voz de forma mais exata, segura e robusta, contribuindo para o surgimento de sistemas de RAL com aplicação prática. Para isso, este trabalho propõe a utilização de novas ferramentas, baseadas na Teoria dos Sistemas Dinâmicos Não-Lineares, para melhorar a caracterização de uma pessoa através de sua voz. Assim, o mecanismo de produção do sinal de voz é analisado sob outro ponto de vista, como sendo o produto de um sistema dinâmico que evolui em um espaço de fases apropriado. Primeiramente, a possibilidade de utilização dessas técnicas em sinais de voz é verificada. A seguir, demonstra-se como as técnicas para estimação de invariantes dinâmicas não-lineares podem ser adaptadas para que possam ser utilizadas em sistemas de RAL. Por fim, adaptações e automatizações algorítmicas para extração de invariantes dinâmicas são sugeridas para o tratamento de sinais de voz. A comprovação da eficácia dessa metodologia se deu pela realização de testes comparativos de exatidão que, de forma estatisticamente significativa, mostraram o benefício advindo das modificações sugeridas. A melhora obtida com o acréscimo de invariantes dinâmicas da forma proposta no sistema de RAL utilizado nos testes resultou na diminuição da taxa de erro igual (EER) em 17,65%, acarretando um intrínseco aumento de processamento. Para sinais de voz contaminados com ruído, o benefício atingido com o sistema proposto foi verificado para relações sinal ruído (SNRs) maiores que aproximadamente 5 dB. O avanço científico potencial advindo dos resultados alcançados com este trabalho não se limita às invariantes dinâmicas utilizadas, e nem mesmo à caracterização de locutores. A comprovação da possibilidade de utilização de técnicas da Teoria do Caos em sinais de voz permitirá expandir os conceitos utilizados em qualquer sistema que processe digitalmente sinais de voz. O avanço das técnicas de Sistemas Dinâmicos Não-Lineares, como a concepção de invariantes dinâmicas mais representativas e robustas, implicará também no avanço dos sistemas que utilizarem esse novo conceito para tratamento de sinais vocais.

Ondas e instabilidades em plasmas com inomogeneidades suaves na densidade e no campo magnético de equilíbrio

No presente trabalho, estudamos a absorção e amplificação de ondas eletromagnéticas que se propagam em plasmas com densidade e temperatura fracamente inomogêneas, imersos em um campo magnético também inomogêneo, tendo como base a teoria cinética, dentro do contexto da aproximação local. Esse estudo se dá efetivamente a partir da obtenção do tensor dielétrico do plasma, que deve ser empregado na relação de dispersão. Iniciamos com uma revisão dos conceitos básicos sobre plasmas homogêneos e inomogêneos. Os fundamentos da teoria cinética também foram abordados. Apresentamos uma revisão de trabalhos anteriores que enfocam o mesmo tema, embora descrevendo separadamente os dois tipos de inomogeneidades. A partir desses trabalhos, obtivemos um tensor dielétrico geral, que descreve de forma simultânea as inomogeneidades do campo magnético de equilíbrio e da função distribuição de equilíbrio.Tal tensor foi obtido a partir de um sólido desenvolvimento teórico, que garante a correta descrição da troca de energia entre as ondas e as partículas do plasma. Abordamos os aspectos gerais das instabilidades de deriva, direcionando o estudo à faixa de frequência das ondas híbridas inferiores, e às instabilidades LHDI e MTSI (IWI). Utilizamos perfis lineares de inomogeneidades de campo magnêtico ambiente e densidade para modelar a região da magnetosfera conhecida como neutral sheet. Particularizamos o tensor dielétrico para o estudo específico das instabilidades LHDI e MTSI (IWI), para o tipo de perfil citado acima. Apresentamos uma nova rela»c~ao de dispers~ao para plasmas inomogêneos, que incorpora explicitamente as derivadas espaciais do tensor dielétrico do plasma. Usamos o tensor que unifica os tratamentos das inomogeneidades do campo e densidade nessa relação de dispersão, e obtivemos uma descrição unificada das instabilidades LHDI e MTSI (IWI).

Dinâmica estocástica de íons sujeitos a um conjunto quase-monocromático de ondas do tipo híbrida inferior

Neste trabalho, estudamos a interação de íons com um conjunto quase-monocromático de ondas eletrostáticas de frequência na faixa das frequências híbridas inferiores, propagando-se perpendicularmente a um campo magnético uniforme. Consideramos que as fases das ondas são aleatoriamente distribuídas (ondas incoerentes), tratando o caso de ondas de fases coerentes (ondas coerentes) como um caso particular. Derivamos o Hamiltoniano adequado a esse sistema, e deduzimos as equações de movimento, cujas soluções são analisadas numericamente, mostrando a ocorrência de difusão estocástica no espaçoo de fase ângulo-ação, para amplitudes de onda suficientemente grandes. Também fazemos estimativas sobre a amplitude mínima (threshold) para o aparecimento de ilhas de primeira ordem no espaço de fase. Estimamos, também, o limiar para as ilhas de segunda ordem e de ordens maiores, bem como o limiar de estocasticidade. A análise mostra que para o caso de várias ondas o comportamento estocástico ocorre antes do limiar de estocasticidade comparado com o caso de uma onda. No caso de ondas coerentes, observa-se que o limiar de estocasticidade diminui com o aumento do número de ondas que comp˜oem o conjunto de ondas, proporcionalmente ao inverso da raiz quadrada deste número, portanto, tendendo a ser nulo no limite em que o número de ondas no pacote tende a infinito. No caso de ondas incoerentes, observa-se também uma diminuição do limiar de estocasticidade com o aumento do número de ondas, mas nesse caso, saturando com valor até um terço do valor do limiar de estocasticidade para o caso de uma onda. Observa-se também que o limite superior da região de estocasticidade no espaço de fase aumenta com o aumento do número de ondas. No caso de ondas coerentes, esse aumento é proporcional à raiz cúbica do número de ondas que compõem o conjunto de ondas. No caso de ondas incoerentes o limite superior da região de estocasticidade têm um aumento de até o dobro em relação ao caso de uma onda. A análise também mostra que o mecanismo da estocasticidade para o caso de várias ondas é diferente do mecanismo atuante no caso de uma onda. No caso de uma onda, a estocasticidade ocorre por superposição de ilhas de ordens maiores do que um, com o aumento da intensidade da onda. No caso de várias ondas, a presençaa de ondas de frequências próximas à frequência de ressonância causa pequenas perturbações na trajetória principal das partículas, causada pela onda central, espalhando-a pelo espaço de fase de forma mais eficiente que o mecanismo de estocasticidade para o caso de uma onda.

Solução de uma classe de equações diferenciais parciais em espaços de Niko'skii através do método de Galerkin

Neste trabalho apresenta-se uma solu c~ao para um problema abstrato de Cauchy. Basicamente, d a-se uma formula c~ao abstrata para certos tipos de equa c~oes diferenciais parciais n~ao lineares de evolu c~ao em espa cos de Nikol'skii, tais espa cos possuem boas propriedades de regularidade e resultados de imers~ao compacta, num certo sentido s~ao intermedi arios entre os espa cos de Holder e os espa cos de Sobolev. Aplicando o m etodo de Galerkin, prova-se resultados de exist^encia global de solu c~oes fracas, como tamb em a exist^encia de solu c~oes fracas com a propriedade de reprodu c~ao. E impondo mais hip oteses sobre os operadores envolvidos demonstra-se unicidade de solu c~oes fracas.