É realizado um estudo sobre o problema do segundo espectro de frequências no modelo de Timoshenko para uma viga bi-apoiada. A equação característica, associada às condições de contorno bi-apoiada, permite determinar dois autovalores que dependem de um inteiro n não negativo, sendo que o de menor módulo está associado ao valor predito pela teoria de Euler-Bernoulli, e o de maior módulo é usualmente referido como sendo do segundo espectro. Este último espectro foi salientado por Traill-Nash e Collar (1953) e desconsiderado por Timoshenko, que utilizou funções trigonométricas que satisfazem as condições de contorno para obter os modos de vibração. Neste trabalho, o modelo de Timoshenko é formulado matricialmente e permite que o estudo dos autovalores e autofunções seja realizado com o uso de uma base da equação modal matricial de segunda ordem completa, gerada por uma solução matricial fundamental. Devido a uma representação analítica desta solução fundamental, o estudo é realizado com o uso da solução de um problema de valor inicial escalar, cujo comportamento torna-se completamente oscilatório acima de um valor crítico. São obtidos resultados que permitem classificar os autovalores como simples e duplos, em ambos os espectros, e determinar seus correspondentes modos. Observa-se que autovalores que correspondem a dois inteiros diferentes e localizados em espectros distintos, porém relativamente próximos, estão associados a modos de vibração descritos por funções trigonométricas que apresentam número de nós bastante diferente. / A study about the problem of the second spectrum of frequencies that arises in the Timoshenko model for a bi-supported beam is accomplished . The characteristic equation associated with the bi-supported boundary conditions allows to determine two eigenvalues that depend of a non-negative integer n, being that the eigenvalue with smaller module is associated to the value predicted by theory of Euler-Bernoulli, and the one with greater module is usually referred as belonging to the second spectrum. This latter spectrum was emphasized by Traill-Nash and Collar (1953) but disregarded by Timoshenko who used trigonometric functions that satisfy the boundary conditions in order to obtain the vibration modes. In this work the model of Timoshenko is formulated in matrix terms and allows that the study of the eigenvalues and eigenfunctions be performed using a basis of a complete second order matrix modal equation, generated by a fundamental matrix solution. Due to an analytical representation of this fundamental solution, the study is done with the solution of a scalar initial value problem, whose behavior becomes completely oscillatory above a critical value. Results are obtained in such a way that allow to classify the eigenvalues as simple and double in both spectra and to determine their corresponding modes. It is observed that eigenvalues that correspond to different integers and localized on distinct spectra, but relatively close together, are associated with vibration modes described by trigonometric functions which have quite different number of nodes.

O presente trabalho apresenta a verificação da Hipótese Topológica (HT) sobre o Modelo Esférico ferromagnético. Esta hipótese sugere uma nova abordagem a transições de fase termodinâmicas (TF), apoiada na idéia recentemente proposta de que o mecanismo na origem das transições de fase seja uma quebra de difeomorficidade adequada do espaço de configurações dos sistemas no ponto da transição. Embora se saiba que nem toda transição topológica (TT) esteja ligada a uma TF, foi provado recentemente para uma classe de modelos que TTs são necessárias na presença de transições de fase de primeira e segunda ordem, pelo menos. Da termodinâmica do modelo esférico campo médio, sabe-se que este apresenta uma TF na ausência de campo magnético externo, e que esta é destruída para campo finito. Da evolução da topologia de subvariedades do espaço de configurações, definidas pela (função) energia potencial, resulta que, no caso sem campo externo, ocorre uma TT nestas variedades em uma energia potencial igual à energia potencial (média) crítica da transição de fase. Este resultado corrobora a HT no sentido da necessidade de uma TT em correspondência a uma TF. No entanto, a mudança observada na topologia não é particularmente “grande”, se considerada dentro do critério sugerido por trabalhos anteriores. Do caso com campo externo, resulta que uma TT idêntica à encontrada no caso sem campo persiste. Da termodinâmica do modelo, observa-se, no entanto, que a região de energias potenciais em que a TT ocorre não é acessível ao sistema físico, o que é consistente com a ausência de TF no modelo neste caso. Entretanto, o formalismo topológico sobre o espaço de configurações não se mostrou suficiente na caracterização deste comportamento do sistema, e esta conclusão só é alcançada dado o conhecimento prévio da termodinâmica do modelo. Apresentam-se também resultados parciais da aplicação da HT ao modelo esférico a dimensão finita. / In the present work the verification of the Topological Hypothesis (TH) on the ferromagnetic Spherical Model is presented. This hypothesis suggests a new approach to phase transitions (PTs), based on the recently proposed idea that the mechanism at the origin of these phenomena might be a suitable break of diffeomorphicity of the systems configuration spaces at the respective transition points. Although it is known that any topological transition (TT) will in general not be linked to a PT, it was recently proved for a class of models the necessity of occurence of a TT in the presence of (at least) first and second order PTs. From the thermodynamics of the mean field Spherical Model, we see that it presents a PT in vanishing external magnetic field, and that it is destroyed by the external field. From the topological evolution of suitable submanifolds in configuration space, defined by the potential energy function, results that, in zero external field, a TT occurs in these submanifolds in a value of potential energy which is the same as the critical (mean) potential energy. This result is in accordance with the TH in the extent of the necessity of a TT in correspondence with a PT. However, the observed TT is not particularly "strong", when considered in the criteria suggested by previous related works. From the finite field case, results that a TT identical to the one in the zero field case persists. From the thermodynamics of the model in the presence of a field, however, it is seen that the potential energy range in which this TT happens is not accessible to the physical system, which is still consistent with the absence of a PT in this case. Nevertheless, the formalism of Differential Topology has been proven insufficient in the characterization of this behaviour, and the conclusion above was only reached given the previous knowledge of the model’s thermodynamics. Some partial results of the application of the TH to the Spherical Model in finite dimension are also presented here.

Este trabalho faz uma pesquisa das vibrações de uma viga elástica não linear de Timoshenko sobre a influência de força axial e com uso do método espectral de Galerkin. O modelo não-linear de Timoshenko é obtido através do principio estendido de Hamilton. A função de energia é derivada de maneira geral, incluindo o caso linear, e com identificação das condições de contorno de natureza conservativa. A determinação das autofunções do sistema linear é realizada através de um problema de autovalor descrito por uma equação diferencial linear de segunda ordem com coeficientes matriciais que dependem não - linearmente no autovalor. A ortogonalidade das autofunções é obtida para as condições de contorno clássicas. As correspondentes autofunções são obtidas na base de Euler e na base gerada pela função matricial de Green de valor inicial. O método de Galerkin foi formulado matricialmente e a existência e unicidade foi obtida para uma viga bi-apoiada com o uso da função da energia e dados regulares. / This work investigates the vibrations of a nonlinear elastic Timoshenko beam, subject to an axial force, using the spectral method of Galerkin. The nonlinear model of Timoshenko is obtained through extended Hamilton s principle. The energy function is derived in a general form, including the linear case, and with identi cation of the boundary conditions of conservative nature. The determination of the eigenfunctions of the linear system is done through an eigenvalue problem described by a second-order di erential equation with matrix coe cients that depend non-linearly upon the eigenvalue. The orthogonality of the eigenfunctions is obtained for classical boundary conditions. The corresponding eigenfunctions are determined in the Euler s base and with the base generated by the initial value matrix Green function. The method of Galerkin was formulated in matrix terms and the existence and uniqueness of solutions was obtained for the bi-supported beam with use of the function of the energy and smooth data.

É amplamente conhecida a crescente utilização de modelos hidrológicos com o fim de representar simplificadamente os complexos processos do ciclo hidrológico em bacias hidrográficas. O ajuste destes modelos passa por uma etapa de calibração onde os valores de seus parâmetros são definidos. Novas ferramentas de otimização vem sendo desenvolvidas com o passar do tempo e surge a necessidade de avaliar seu desempenho na calibração automática de modelos hidrológicos que utilizamos no dia-a-dia como, por exemplo, o IPH II. Nesse contexto, nessa pesquisa foi desenvolvido um algoritmo baseado no comportamento de colônias de formigas, denominado ACO-IPH. Os algoritmos baseados no comportamento de colônia de formigas, desde a sua conceição foram projetados para fins de otimização e testados, até a atualidade, em diversas aplicações, onde apresentaram desempenhos aceitáveis, embora as aplicações na calibração de modelos hidrológicos sejam escassas. O algoritmo ACO-IPH foi inicialmente testado em funções matemáticas com ótimos conhecidos apresentando um bom desempenho em termos de eficiência e eficácia. Posteriormente foi acoplado ao modelo hidrológico IPH II como uma alternativa na calibração automática mono-objetivo e verificado seu desempenho e as simplificações e limitações de sua utilização. Uma vez acoplado, o algoritmo foi testado na calibração automática do modelo IPH II utilizando uma série sintética de vazões e, posteriormente, utilizando séries reais de vazões em quatro bacias hidrográficas. O desempenho do algoritmo foi avaliado em termos de eficácia e eficiência, para três diferentes funções-objetivo e utilizando como critério de comparação o algoritmo SCE-UA. Os diferentes testes executados na calibração automática do modelo hidrológico IPH II mostraram que o ACO-IPH consegue encontrar boas soluções, porém não consegue apresentar um desempenho, em termos de eficiência e eficácia, comparável a algoritmos do estado da arte atual, como o algoritmo SCE-UA. A perda de eficiência do ACO-IPH em comparação ao SCE-UA é a mais importante, podendo demorar, em média, até 10 vezes mais até a convergência. / It is widely known the increasing use of hydrologic models in order to represent the complex processes of the hydrological cycle in watersheds. To adjust these models a calibration step is necessary, where the values of its parameters are defined. New optimization tools have been recently developed that could be used for the automatic calibration of hydrologic models as the IPH II. In this study an algorithm was developed based on the behavior of ant colonies, named ACO-IPH. Since its conception, the mathematical algorithms based on ants foraging behavior were designed for optimization purposes and tested, until today, in several applications, where it showed acceptable performance, although the applications for the calibration of hydrologic models are scarce. The ACO-IPH algorithm was initially tested on mathematical functions where the global optimum was known showing a good performance in terms of efficiency and effectiveness.Then it was coupled to the hydrologic model IPH II as an alternative for the single-objective automatic calibration. Several indices were estimated to assess the ACO-IPH performance after analyzing the simplifications and limitations of its use. Once coupled to the IPH II model, the algorithm was tested in the calibration processes using synthetic streamflow series and later using real series of observed streamflows in four watersheds. The algorithm's performance was assessed in terms of effectiveness and efficiency for three different objective functions and using as a basis of comparison the SCE -UA algorithm. The different tests performed in automatic calibration of hydrological model IPH II showed that ACO-IPH can find good solutions, but cannot reach the efficiency and effectiveness of other evolutionary algorithms like the SCE-UA. The loss of effectiveness of ACO-IPH compared to the SCE-UA is the most relevant. The ACO-IPH took, on average, 10 times longer toconvergence.

Esta dissertação, tem como propósito, estimar a elasticidade intertemporal para a oferta de trabalho masculina no Brasil. Também, verifica-se a presença de fatos estilizados relacionados ao ciclo de vida do indivíduo. Este trabalho fundamenta-se na teoria que tem como referência Heckman (1974 e 1976), e a idéia principal é supor que o indivíduo procura aumentar sua oferta de trabalho em períodos onde o salário real é alto. A partir de um modelo com dados de painel, formulado por MaCurdy (1981), procedemos a estimação. Para a verificação dos fatos estilizados foram feitos gráficos. No caso da estimação, houve um tratamento explícito para as questões da forma funcional, heterogeneidade não observada e da endogeneidade da variável salário. Os gráficos confirmam, em parte, a presença dos fatos estilizados, as variáveis salário por hora e rendimentos do trabalho têm uma relação côncava com a variável idade do indivíduo, no mercado de trabalho brasileiro. Contudo, o modelo econométrico estimado não confirma a efetividade da hipótese de elasticidade intertemporal, uma vez que, as estimativas geradas foram baixas e negativas, embora significativas .

O objetivo deste trabalho é a introdução e desenvolvimento de uma metodologia analítico-simbólica para a obtenção de respostas dinâmicas e forçadas (soluções homogêneas e não homogêneas) de sistemas distribuídos, em domínios ilimitados e limitados, através do uso da base dinâmica gerada a partir da resposta impulso. Em domínios limitados, a resposta impulso foi formulada pelo método espectral. Foram considerados sistemas com condições de contorno homogêneas e não homogêneas. Para sistemas de natureza estável, a resposta forçada é decomposta na soma de uma resposta particular e de uma resposta livre induzida pelos valores iniciais da resposta particular. As respostas particulares, para entradas oscilatórias no tempo, foram calculadas com o uso da fun»c~ao de Green espacial. A teoria é desenvolvida de maneira geral permitindo que diferentes sis- temas evolutivos de ordem arbitrária possam ser tratados sistematicamente de uma forma compacta e simples. Realizou-se simulações simbólicas para a obtenção de respostas dinâmicas e respostas for»cadas com equações do tipo parabólico e hiperbólico em 1D,2D e 3D. O cálculo das respostas forçadas foi realizado com a determinação das respostas livres transientes em termos dos valores iniciais das respostas permanentes. Foi simulada a decomposição da resposta forçada da superfície livre de um modelo acoplado oceano-atmosfera bidimensional, através da resolução de uma equação de Klein-Gordon 2D com termo não-homogêneo de natureza dinâmica, devido a tensão de cisalhamento na superfície do oceano pela ação do vento. / The objective of this work is the introduction and the development of an analytical-symbolic methodology for obtaining dynamic and forced responses of distributed systems, in unlimited and limited domains, through the use of the dynamic basis generated by the impulse response. In limited domains, the impulse response was formulated with the spec- tral method. Systems were considered with homogeneous and nonhomogeneous boundary conditions. For systems of stable nature, the forced response is decomposed into the sum of a particular response and a free response that is induced by the initial values of the particular response. The particular responses, for oscillatory time inputs, were calculated with the use of the spatial Green's function. The theory is developed in a general way that allows that di®erent arbitrary order evolution systems can be systematically treated in a compact and simple way. Symbolic simulations ware performed for obtaining dynamics responses and forced responses with equations of parabolic and hyperbolic type in 1D,2D and 3D. The calculation of the forced responses were accomplished with the determina- tion of the transient free responses in terms of the initial values of the permanent responses. It was simulated the decomposition of the forced response of the free surface of a two-dimensional coupled model ocean-atmosphere, through the inte- gration of a 2D Klein-Gordon equation with a nonhomogeneous terms of dynamic nature, due to the shear stress in the surface of the ocean by the action of the wind.

O agronegócio da soja, caracterizado como investimento em atividade produtiva, envolve problemas de incertezas com o comprometimento definitivo de recursos, com a impossibilidade de reverter decisões tomadas e, também, com o esforço de prever a renda a ser obtida com a colheita. A questão relevante para o sojicultor é como avaliar o rendimento esperado de grãos em um contexto de incertezas, base para o cômputo da taxa de retorno do investimento na produção de soja. Este trabalho teve como objetivo estabelecer uma metodologia que associa um modelo estocástico de precipitação a um modelo de predição do rendimento da cultura de grãos de modo a permitir a atribuição de probabilidades aos rendimentos esperados de soja. Foram utilizados dados meteorológicos diários referentes à precipitação pluvial, temperatura máxima, mínima e média do ar, umidade relativa do ar, radiação solar global e velocidade do vento a 2 metros, bem como o rendimento de soja de Cruz Alta(RS) do período 1/1/1975-31/12/1997 A precipitação pluvial diária foi ajustada à cadeia de Markov de primeira ordem e dois estados e na função gama com variação sazonal para os meses de outubro a abril. Os resultados indicaram que a metodologia associando a geração de séries de precipitação pluvial, evapotranspiração da soja e evapotranspiração real, aplicadas como variáveis de entrada no modelo de predição do rendimento de grãos, possibilitou a definição da densidade de probabilidade do rendimento esperado de soja e, por conseguinte, a quantificação do risco na forma da distribuição de probabilidade acumulada do rendimento. O risco de ocorrência de um valor de rendimento de soja, menor que um valor mínimo de aceitação, quantiifcado pela sua distribuição de probabilidade acumulada, fundamentada em um fenômeno da natureza - a precipitação pluvial - que não depende da vontade humana. A obtenção do rendimento desejado envolve um processo que depende de organização, isto é, depende de fatores subjetivos, de decisões e escolhas. Enfim, um processo que comporta alternativas e depende de iniciativas humanas.

Submitted by Israel Vieira Neto (israel.vieiraneto@ufpe.br) on 2015-03-06T16:24:43Z No. of bitstreams: 2 TESE Carlos Alberto Gomes de Amorim Filho.pdf: 1106273 bytes, checksum: d3ecb6d6b657e535fd2dbe040fd08d5b (MD5) license_rdf: 1232 bytes, checksum: 66e71c371cc565284e70f40736c94386 (MD5) / Made available in DSpace on 2015-03-06T16:24:43Z (GMT). No. of bitstreams: 2 TESE Carlos Alberto Gomes de Amorim Filho.pdf: 1106273 bytes, checksum: d3ecb6d6b657e535fd2dbe040fd08d5b (MD5) license_rdf: 1232 bytes, checksum: 66e71c371cc565284e70f40736c94386 (MD5) Previous issue date: 2013-09-17 / O crescimento econômico desordenado e as persistentes agressões ambientais, juntamente com a ausência de ações eficazes de gerenciamento e preservação dos recursos hídricos, vêm limitando a disponibilidade dos mesmos. Isto restringe o atendimento aos usos e exacerba os conflitos entre os usuários de água. É nesse contexto que se encontra a bacia hidrográfica do rio Tapacurá, cujo reservatório é responsável pelo abastecimento público de cerca de 1,5 milhões de habitantes da Região Metropolitana do Recife e, apesar de sua grande importância, o mesmo tem sofrido com problemas de poluição e de disponibilidade hídrica. Neste estudo foi desenvolvido um modelo hidro-econômico para a bacia do rio Tapacurá capaz de determinar uma alocação econômica ótima de água para os seus usuários. O modelo incorpora oferta de água que simula condições climáticas atuais e futuras, e as demandas presentes. Busca-se dessa forma contribuir para o melhor planejamento e gerenciamento dos recursos hídricos desta região, de modo a minimizar os potenciais conflitos pela água que ali possam eventualmente ocorrer e possibilitar seu desenvolvimento sustentável. Os resultados mostraram que, utilizando dados de um ano hidrológico médio para um cenário com as condições climáticas atuais, as restrições de quantidade de água não comprometeram o atendimento das demandas requeridas pelos usuários. No entanto, ao serem consideradas também as restrições de qualidade da água, houve uma redução na quantidade total da água alocada aos usuários, o que resulta em uma diminuição do benefício econômico para esses usuários. Apesar do modelo ter sido desenvolvido inicialmente para ser aplicado num horizonte de tempo de um ano com passos mensais, o mesmo foi estendido e aplicado para um período de 3 anos, possibilitando uma gestão da água interanual. Além das simulações para um cenário atual, o modelo foi executado para cenários futuros, considerando mudanças climáticas. Nesse caso, o modelo foi aplicado em duas situações: utilizando dados de um ano representativo do período de 2010 a 2039, em que foi considerado uma diminuição das vazões de 4,98% em relação a um ano médio do cenário atual, e do período de 2040 a 2069, em que a redução considerada foi de 14,28%. Para o primeiro período (2010-2039), às restrições de quantidade não comprometeram o atendimento das demandas requeridas pelos usuários. Entretanto, ao ser executado o modelo com as restrições de qualidade, tem-se uma redução nos benefícios dos usuários. Já para o período de 2040 a 2069, o atendimento às restrições de disponibilidade hídrica resultam em uma pequena redução dos benefícios líquidos dos usuários, contudo, ao serem adicionadas as restrições de qualidade essa redução passa a ser maior do que a do período de 2010 a 2039.

Submitted by Eduarda Figueiredo (eduarda.ffigueiredo@ufpe.br) on 2015-03-13T12:14:36Z No. of bitstreams: 2 Tese Andre Inacio do Carmo.pdf: 1646519 bytes, checksum: 919f7ad3394f397b2cb42070d27de624 (MD5) license_rdf: 1232 bytes, checksum: 66e71c371cc565284e70f40736c94386 (MD5) / Made available in DSpace on 2015-03-13T12:14:36Z (GMT). No. of bitstreams: 2 Tese Andre Inacio do Carmo.pdf: 1646519 bytes, checksum: 919f7ad3394f397b2cb42070d27de624 (MD5) license_rdf: 1232 bytes, checksum: 66e71c371cc565284e70f40736c94386 (MD5) Previous issue date: 2012-06-20 / FACEPE / A contaminação de solos na Região Metropolitana da Cidade do Recife, PE decorrente de vazamento em postos de gasolina tem acontecido, principalmente quando os postos de gasolina são mais antigos e quando não há monitoramento adequado das condições dos tanques de armazenamento. Neste sentido foi estudada a situação de três solos desta região com o objetivo de compreender os mecanismos envolvidos na sorção e no transporte do Naftaleno, um dos componentes dos Hidrocarbonetos Aromáticos Policíclicos (HPAs) dos combustíveis de petróleo. Foram realizados ensaios cinéticos e isotérmicos, em batelada, bem como, ensaios de transporte em colunas de solo deformado, em laboratório. Nos três solos, as cinéticas de adsorção foram não lineares e as isotermas de sorção lineares. Os valores do KD em condições naturais (na presença da matéria orgânica) foram 3,7; 26,8 e 11,2 cm3 g-1, respectivamente, para o solo Franco Arenoso 1, Franco Arenoso 2 e Franco Argilo Siltoso, e 0,43; 2,1 e 0,91 cm3 g-1, respectivamente na ausência da matéria orgânica. O solo Franco Arenoso 1 apresentou menor capacidade de adsorção ao Naftaleno, e o solo Franco Arenoso 2, apresentou a maior capacidade. Verificou-se ainda que para os três solos, a interação da matéria orgânica com o Naftaleno foi bastante significativa, sendo que no solo Franco Arenoso 2 foi observada uma forte influência dos óxidos de ferro (goethita) na adsorção do Naftaleno. Os ensaios de transporte com o Naftaleno foram realizados na vazão de 0,45 cm3 min-1 e na concentração de 20 mg L-1. As variáveis hidrodispersivas foram obtidas pelo ajuste do modelo convecção dispersão (CDE) às curvas de eluição do Naftaleno por intermédio do programa CXTFIT 2.0. Em condições naturais, na presença da matéria orgânica, o solo Franco Arenoso 1 foi o mais dispersivo dos três solos e o menos reativo ao Naftaleno, ou seja, foi o solo que apresentou a maior susceptibilidade à contaminação a esse composto. Na ausência de matéria orgânica, o solo Franco Arenoso 2 foi o mais dispersivo dos três solos e o mais reativo ao Naftaleno, tendo apresentado assim, um bom potencial de adsorção, dentre os solos estudados. O modelo CDE representou adequadamente os dados experimentais das curvas de eluição do Naftaleno. Os resultados obtidos evidenciaram que o Naftaleno, no solo Franco Arenoso 1, mesmo com a presença da matéria orgânica, ofereceu um risco real de contaminação das águas subterrâneas da Região Metropolitana do Recife, PE. Esses resultados permitiram estimar o risco do transporte reativo do Naftaleno em três solos urbanos da Região Metropolitana do Recife, PE.

Submitted by Eduarda Figueiredo (eduarda.ffigueiredo@ufpe.br) on 2015-03-13T14:38:36Z No. of bitstreams: 2 TESE Priscila Alves Lima.pdf: 1194515 bytes, checksum: b995d6c97f419dbf8969c321268318bb (MD5) license_rdf: 1232 bytes, checksum: 66e71c371cc565284e70f40736c94386 (MD5) / Made available in DSpace on 2015-03-13T14:38:36Z (GMT). No. of bitstreams: 2 TESE Priscila Alves Lima.pdf: 1194515 bytes, checksum: b995d6c97f419dbf8969c321268318bb (MD5) license_rdf: 1232 bytes, checksum: 66e71c371cc565284e70f40736c94386 (MD5) Previous issue date: 2012-10-10 / A aplicação de esterco bovino em áreas de lavoura representa a adição de nutrientes às plantas e também uma alternativa de reciclagem. Contudo, as doses excessivas, ocasionam acúmulo e potencialmente perdas de fósforo por lixiviação, reduzindo a eficiência à nutrição de plantas e representam potencial contaminante do ambiente. O objetivo deste trabalho foi avaliar o transporte de P em solo arenoso adubado com esterco bovino. O solo estudado foi Neossolo Regolítico coletado em três camadas, 0-20; 20-40 e 40-60 cm no munícipio de Esperança-PB, com histórico de adubação orgânica aproximadamente de 40 anos. Foi aplicado o método de deslocamento miscível em colunas de solo com 20 cm de comprimento; 2,6 cm de diâmetro utilizando bomba peristáltica com vazão de 0,286 cm3 min-1. Por intermédio das curvas de eluição do KBr e do P identificou o mecanismo do processo de transporte do P e determinou os parâmetros hidrodispersivos por meio dos modelos convecção-dispersão (CDE), CDE com duas regiões de água Móvel e Imóvel (CDE-MIM) por intermédio do programa CXTFIT 2.0. Os parâmetros estimados verificaram que o transporte do P, a sorção foi mais lenta que a dessorção e o modelo CDE a dois Sítios de Sorção Tipo-1 representou adequadamente os dados experimentais das curvas de eluição do P. Os resultados obtidos mostram que a adubação com esterco bovino, em solo arenoso, oferece risco real de contaminação das águas subterrâneas da região.