Solución en diferencias finitas de la ecuación de Boussinesq con porosidad drenable variable y condición de radiación fractal en la frontera

Abstract

El drenaje subterráneo es utilizado para eliminar excedentes de agua en la zona radical y suelos salinos para lixiviar las sales. La dinámica del agua es estudiada con la ecuación de Boussinesq, sus soluciones analíticas son obtenidas asumiendo que la transmisibilidad del acuífero y la porosidad de drenable son constantes y que la superficie libre se abate de manera instantánea sobre los drenes. La solución en el caso general requiere de soluciones numéricas. Se ha mostrado que la condición de frontera en los drenes es una condición de radiación fractal y la porosidad drenable es variable y relacionada con la curva de retención de humedad, y ha sido resuelta con el método del elemento finito, que en un esquema unidimensional puede hacerse equivalente al método de diferencias finitas. Aquí se propone una solución en diferencias finitas de la ecuación diferencial considerando la porosidad drenable variable y la condición de radiación fractal.