UNA NUEVA ALTERNATIVA NUMÉRICA PARA LA SOLUCIÓN DE LA ECUACIÓN UNIDIMENSIONAL DE RICHARDS: ESTUDIO DE DRENAJE E INFILTRACIÓN DE FLUIDOS EN LA ZONA NO SATURADA

Evelyn Álvarez Sierra, Lilliam Álvarez Díaz, Maria Teresa Alonso González

Resumen


En el trabajo se comparan dos métodos numéricos para resolver el modelo unidimensional de infiltración y drenaje de agua en la zona no saturada en medios porosos, el cual es modelado respecto al contenido de humedad utilizando la ecuación no-lineal de Richards. El primer método está basado en el método clásico de Diferencias Finitas y el segundo en el método de Líneas combinado con el código DASSL para la solución de las ecuaciones diferenciales-algebraicas resultantes. Se muestra que el último método proporciona una vía numérica eficiente para la solución de problemas de EDPs que tienen un comportamiento singular de shock u ondas viajeras, como es el caso de la ecuación de Richard, los cuales se pueden resolver numéricamente con éxito sólo utilizando esquemas muy estables. Los métodos numéricos discutidos en el trabajo se aplican a dos tipos de suelos reales: Yolo Light Clay y Brindabella Silty Clay Loam, usando las propiedades hidráulicas referidas en Broadbridge y White [16]. Para validar el modelo, se comparan los perfiles de humedad con los resultados reportados por Warrick, Lomen e Islas [18].Igualmente, se demuestra con datos reales la ventaja de resolver numéricamente el comportamiento del flujo unidimensional en la zona no saturada de un medio poroso. Finalmente, el modelo propuesto y los resultados numéricos obtenidos posibilitan brindar un pronóstico sobre la utilización de recursos hídricos para el caso de riego en agricultura y también para el transporte de contaminantes. 


Palabras clave


Ecuación de Richards, Drenaje, Infiltración de Fluidos, Zonas Saturadas

Texto completo:

RESUMEN ARTICULO COMPLETO

Referencias


L. A. Richards, “Capillary conduction of liquids in porous media,” in Physics 1, no 5, 1931, pp. 318-333.

R. G. Hills et al. “Modeling one-dimensional infiltration into very dray soils,” Water Resources Research, 1989, pp. 181-187.

P. J. Ross, “Efficient numerical methods for infiltration using Richards equation,” Water Resources Research, vol. 26, 1990, pp.279-290.

M. R. Kirkland and R.G. Hills, “Algorithms for solving Richards equation for variably saturated soils”, Water Resources Research, 1992.

B. N. Wilson and D. C. Slack. A comparison of three Infiltration Models, 1981.

A. G. Williams et al. Approaches to Modeling Richards Equation, 1996.

V. Leitâo et al. “The dual boundary elements formulation for elastopalastic fractury mechanics,” International Journal for Numerical Methods in Engineering.

K. E. Brenam et al. Numerical solution of initial value problems of differential-algebraic equations, Elsevier Publs., 1989.

L. R. Petzold. A description of DASSL, a differentialalgebraic solver, SANDIA Report, SAND-828637, 1982.

E. Griepentrog and R. März. Differential-Algebraic Equations and their Numerical Treatment, Teubner, Leipzig, 1986.

Hindmarsh, LSODE and LSODI. Two initial value ordinary differential equation solvers, ACM SIGNUM Newsletter, 15, 1980, pp10-11.

L. Shampine and M. W. Reichelt. The Matlab ODE, Suite, 1999.

C. Cunha and L. Alvarez. Numerical solution of the Rapoport-Leas equation using the method of lines and the DASSL code, Relatorio de Pesquisa No.50/1997, IMECC, Brazil; Journal on Computational Mathematics and Math. Phys, vol. 11, 1998, pp. 1851-1858.

J. M. Sanz – Serna. Converge analysis of one-step schemes in the method of lines, Applied Math and Comput. 1993, pp 183-195.

P. Broadbrige and I. White. “Constant rate rainfall infiltration: A versatile nonlinear model, 1. Analytical solution,” Water Resources Res., 24, NO.1145-154, 1988(a).

I. White and P. Broadbrige. “Constant rate rainfall infiltration: A versatile nonlinear model, 2: Applications of Solutions,” Water Resources Research, 24, 1988, pp. 155-162.

A. A. Samarski. Introducción a los Métodos Numéricos, Ed. Mir, 1982.

A. W. Warrick et al. “An analytical solution to Richards Equation for draining soil profile,” Water Resources Research, vol. 26, No. 2, 1990, pp. 253-258.


Enlaces refback

  • No hay ningún enlace refback.


ESTADÍSTICAS DEL ARTICULO
Resumen : 133
ARCHIVO PDF RESUMEN : 33
ARCHIVO PDF ARTICULO COMPLETO : 53



Copyright (c) 2018 Revista Investigación & Desarrollo

Licencia de Creative Commons
Este obra está bajo una licencia de Creative Commons Reconocimiento-NoComercial-CompartirIgual 4.0 Internacional.