Solution of a time fractional inverse advection-dispersion problem by discrete mollification

Abstract

ABSTRACT: We consider an inverse problem for a time fractional advection-dispersion equation in a 1-D semi-infinite setting. The fractional derivative is interpreted in the sense of Caputo and advection and dispersion coefficients are constant. The inverse problem consists on the recovery of the boundary distribution of solute concentration and dispersion flux from measured (noisy) data known at an interior location. This inverse problem is ill-posed and thus the numerical solution must include some regularization technique. Our approach is a finite difference space marching scheme enhanced by adaptive discrete mollification. Error estimates and illustrative numerical examples are provided.

RESUMEN: Consideramos un problema inverso para una ecuación de advección-dispersión con derivada temporal fraccionaria, en una configuración unidimensional. La derivada fraccionaria se interpreta en el sentido de Caputo y las coeficientes de advección y de dispersión son constantes. El problema inverso involucra la reconstrucción simultánea de la concentración de soluto y del flujo de dispersión en una de las fronteras del dominio físico, a partir de lecturas de datos perturbados en un punto interior del dominio. Mostramos que el problema inverso es mal condicionado y por tanto una solución numérica del problema requiere de alguna técnica de regularización. Proponemos un esquema de diferencias finitas de marcha en el espacio, que utiliza molificación discreta como técnica de regularización. Se incluyen estimativos de error y ejemplos numéricos ilustrativos.