O conhecimento da distribuição da umidade no solo para diferentes tipos de solos e diferentes vazões pode ser aplicado na otimização do dimensionamento de um sistema de irrigação, no manejo da água na zona radicular, bem como, auxiliar na aplicação eficiente de fertilizantes na rizosfera. Portanto, objetivou-se com a presente pesquisa desenvolver um modelo numérico bidimensional capaz de simular a distribuição da umidade no perfil de solo para diferentes sistemas de irrigação. Sendo que para tal, o modelo utiliza o método das diferenças finitas, mediante uma discretização da solução da equação de Richards e o método iterativo de Picard, o qual foi utilizado para garantir a conservação da massa. O modelo numérico ainda torna possível considerar a extração de água pelas plantas e a evaporação da água na superfície do solo. E para avaliar o seu desempenho foi feita uma análise de sensibilidade. O modelo matemático foi desenvolvido no Departamento de Engenharia de Biossistemas, pertencente à Escola Superior de Agricultura \"Luiz de Queiroz\" - ESALQ/USP. A linguagem de programação utilizada foi Fortran 90. A estrutura computacional foi feita de forma a permitir ao usuário a entrada de informações tais como: a) dados do perfil do solo, no que se refere às suas propriedades físico-hídricas, b) informações sobre o tipo de irrigação, tempo de simulação, tempo de aplicação de água via irrigação. O programa permite ao usuário a opção entre simular a extração de água pela as plantas ou não. Para a validação do modelo foram utilizados os dados obtidos por Rivera (2004), em um experimento conduzindo no antigo Departamento de Engenharia Rural da Escola Superior de Agricultura \"Luiz de Queiroz\". O desempenho do modelo foi avaliado com base nos parâmetros estatísticos: índice de concordância de Willmott; coeficiente de eficiência, raiz quadrada do erro médio, erro médio e o erro máximo absoluto. A análise de sensibilidade do modelo foi avaliada em função dos parâmetros: densidade de fluxo de água, condutividade hidráulica do solo saturado, alfa e n (parâmetros de ajuste da curva de retenção). Pelos resultados obtidos, o modelo apresentou um desempenho satisfatório na simulação dos perfis de umidade em relação aos dados medidos em condições de experimento em solo não saturado. Em relação à análise de sensibilidade, o modelo apresentou maior sensibilidade aos parâmetros de fluxo de entrada e ao parâmetro n da equação de retenção, indicando que tais parâmetros de entrada, necessitam ser determinados com maior precisão. Houve baixa sensibilidade ao parâmetro \"alfa\" da curva de retenção. / The soil moisture distribution for different soil types under different flow regimes can be used to optimize the design of an irrigation system, water manage water in the root zone, or increase the precision of the application of fertilizers in the rhizosphere. Therefore, the objective of this research was to develop a two-dimensional numerical model capable of simulating the soil moisture distribution in soil profile for different irrigation systems. The model is a finite difference solution of the Richards equation solution, in which the Picard iteration method is used to ensure the conservation of mass. The numerical model also takes water extraction by plants and water evaporation on the soil surface into account. The mathematical model was developed in the Department of Biosystems Engineering, College of Agriculture \"Luiz de Queiroz\", ESALQ / USP. It was coded in Fortran 90, and structured to allow interactive input of a) physical and hydraulic properties soil profile, and b) type of irrigation, irrigation application time, and simulation. The simulation of water extraction by the plants is optional. And to evaluate its performance a sensitivity analysis was done. Data obtained in an experiment carried out in the older Rural Engineering Department, College of Agriculture \"Luiz de Queiroz\" by Rivera (2004), were used to validate the model. Model performance was evaluated using the Willmott index of agreement; efficiency coefficient, root mean square error, average error and the maximum absolute error. Sensitivity analyses were performed on the density of water flow, saturated hydraulic conductivity, and on alpha and n (setting parameters of the soil water retention curve). From the results obtained, the model satisfactory simulated the measured soil moisture distributions in the experimental soil moisture profiles. The model was most sensitive to variations in the parameters of the input stream and to n, a parameter in the soil water retention equation (van Genuchten model, 1980). This results is indicative that these input parameters need to be determined with greater accuracy. The model output was not very sensitive to alpha, a shape parameter in the soil water retention curve.

O conhecimento da distribuição da umidade no solo para diferentes tipos de solos e diferentes vazões pode ser aplicado na otimização do dimensionamento de um sistema de irrigação, no manejo da água na zona radicular, bem como, auxiliar na aplicação eficiente de fertilizantes na rizosfera. Portanto, objetivou-se com a presente pesquisa desenvolver um modelo numérico bidimensional capaz de simular a distribuição da umidade no perfil de solo para diferentes sistemas de irrigação. Sendo que para tal, o modelo utiliza o método das diferenças finitas, mediante uma discretização da solução da equação de Richards e o método iterativo de Picard, o qual foi utilizado para garantir a conservação da massa. O modelo numérico ainda torna possível considerar a extração de água pelas plantas e a evaporação da água na superfície do solo. E para avaliar o seu desempenho foi feita uma análise de sensibilidade. O modelo matemático foi desenvolvido no Departamento de Engenharia de Biossistemas, pertencente à Escola Superior de Agricultura \"Luiz de Queiroz\" - ESALQ/USP. A linguagem de programação utilizada foi Fortran 90. A estrutura computacional foi feita de forma a permitir ao usuário a entrada de informações tais como: a) dados do perfil do solo, no que se refere às suas propriedades físico-hídricas, b) informações sobre o tipo de irrigação, tempo de simulação, tempo de aplicação de água via irrigação. O programa permite ao usuário a opção entre simular a extração de água pela as plantas ou não. Para a validação do modelo foram utilizados os dados obtidos por Rivera (2004), em um experimento conduzindo no antigo Departamento de Engenharia Rural da Escola Superior de Agricultura \"Luiz de Queiroz\". O desempenho do modelo foi avaliado com base nos parâmetros estatísticos: índice de concordância de Willmott; coeficiente de eficiência, raiz quadrada do erro médio, erro médio e o erro máximo absoluto. A análise de sensibilidade do modelo foi avaliada em função dos parâmetros: densidade de fluxo de água, condutividade hidráulica do solo saturado, alfa e n (parâmetros de ajuste da curva de retenção). Pelos resultados obtidos, o modelo apresentou um desempenho satisfatório na simulação dos perfis de umidade em relação aos dados medidos em condições de experimento em solo não saturado. Em relação à análise de sensibilidade, o modelo apresentou maior sensibilidade aos parâmetros de fluxo de entrada e ao parâmetro n da equação de retenção, indicando que tais parâmetros de entrada, necessitam ser determinados com maior precisão. Houve baixa sensibilidade ao parâmetro \"alfa\" da curva de retenção. / The soil moisture distribution for different soil types under different flow regimes can be used to optimize the design of an irrigation system, water manage water in the root zone, or increase the precision of the application of fertilizers in the rhizosphere. Therefore, the objective of this research was to develop a two-dimensional numerical model capable of simulating the soil moisture distribution in soil profile for different irrigation systems. The model is a finite difference solution of the Richards equation solution, in which the Picard iteration method is used to ensure the conservation of mass. The numerical model also takes water extraction by plants and water evaporation on the soil surface into account. The mathematical model was developed in the Department of Biosystems Engineering, College of Agriculture \"Luiz de Queiroz\", ESALQ / USP. It was coded in Fortran 90, and structured to allow interactive input of a) physical and hydraulic properties soil profile, and b) type of irrigation, irrigation application time, and simulation. The simulation of water extraction by the plants is optional. And to evaluate its performance a sensitivity analysis was done. Data obtained in an experiment carried out in the older Rural Engineering Department, College of Agriculture \"Luiz de Queiroz\" by Rivera (2004), were used to validate the model. Model performance was evaluated using the Willmott index of agreement; efficiency coefficient, root mean square error, average error and the maximum absolute error. Sensitivity analyses were performed on the density of water flow, saturated hydraulic conductivity, and on alpha and n (setting parameters of the soil water retention curve). From the results obtained, the model satisfactory simulated the measured soil moisture distributions in the experimental soil moisture profiles. The model was most sensitive to variations in the parameters of the input stream and to n, a parameter in the soil water retention equation (van Genuchten model, 1980). This results is indicative that these input parameters need to be determined with greater accuracy. The model output was not very sensitive to alpha, a shape parameter in the soil water retention curve.

Made available in DSpace on 2016-08-31T13:24:20Z (GMT). No. of bitstreams: 1 LuizAAS_DISSERT.pdf: 2790394 bytes, checksum: 1bac58cbcfc89d893d66e58cfa50081f (MD5) Previous issue date: 2009-03-06 / In this work, a mathematical model was developed for the simulation of soil water movement under surface drip irrigation, being used the method of the finite volumes for the resolution of the partial differential equation of water flow in porous media. The soil water flow was treated three-dimensionally for the determination of your movement in the infiltration and redistribution phases. The model allows the determination dimensions and shape of the wetted bulb, of the partial and final dimension of the saturated disc radius and the determination of the saturated soil volume, when to deal about of isolated bulb. In the case of overlapping bulbs, the model allows the determination of the width and depth of the formed wetted strip, making possible the choice of the drippers spacing that better it is adapted to the conditions of the irrigation project and management. A computational program was developed with base in the model, and the obtained results were validated starting from comparisons with field data and results of simulations of other models. These comparisons demonstrated that the model presents reliable results and it can be used as tool for the designing of drip irrigation installations. The sensibility analysis accomplished starting from the variation of some parameters of the soil and of the emitter (initial moisture, saturated soil hydraulic conductivity and trickle discharge), it demonstrated that, with relationship to the final moisture inside of the bulb, the model is so much relatively sensitive to positive variations as negatives of those parameters. / Neste trabalho, foi desenvolvido um modelo matemático para a simulação do movimento de água no solo sob irrigação por gotejamento superficial, utilizando-se o método dos volumes finitos para a resolução da equação diferencial parcial de escoamento de água em meios porosos. O fluxo de água no solo foi tratado tridimensionalmente para a determinação do seu movimento nas fases de infiltração e redistribuição. O modelo permite a determinação da forma e das dimensões do bulbo molhado, da dimensão parcial e final do raio do disco saturado e a determinação do volume de solo saturado, quando tratar-se de bulbo isolado. No caso de bulbos sobrepostos, o modelo permite a determinação da largura e profundidade da faixa molhada formada, possibilitando a escolha do espaçamento entre gotejadores que melhor se adeque às condições do projeto e manejo da irrigação. Um programa computacional foi desenvolvido com base no modelo, e os resultados obtidos foram validados a partir de comparações com dados de campo e resultados de simulações de outros modelos. Estas comparações demonstraram que o modelo apresenta resultados confiáveis e pode ser utilizado como ferramenta para o dimensionamento de instalações de irrigação por gotejamento. A análise de sensibilidade realizada a partir da variação de alguns parâmetros do solo e do emissor (umidade inicial, condutividade hidráulica do solo saturado e vazão do gotejador), demonstrou que, com relação à umidade final dentro do bulbo, o modelo é relativamente sensível tanto à variações positivas quanto negativas desses parâmetros.