La difusión es el movimiento neto de partículas desde regiones de mayor concentración hacia regiones de menor concentración. Cuando el proceso estocástico subyacente no es browniano, la difusión se dice anómala. En particular, la superdifusión anómala se modela mediante operadores de diferenciación espacial de orden no-entero, que son no-locales y de carácter integro-diferencial. En este proyecto, estudiamos propiedades analíticas y desarrollamos y analizamos herramientas computacionales para el tratamiento de problemas con operadores de este tipo sobre dominios acotados. Motivados por diversas aplicaciones, tratamos con problemas no-lineales para operadores no-locales. En estas aplicaciones, el relajamiento causado por la incorporación de modelos no-locales permitiría capturar fenómenos que sus contrapartes locales no logran capturar completamente. La no-localidad conlleva desafíos tanto analítica como computacionalmente. El desarrollo de métodos numéricos se ve dificultado por la presencia de núcleos hipersingulares y la necesidad de integrar sobre dominios no acotados. Mostramos que las soluciones de los problemas correspondientes sean poco regulares, fundamentalmente debido a un pobre comportamiento cerca de la frontera del dominio. Esta regularidad de soluciones es un elemento fundamental a tener en cuenta para el análisis de los métodos numéricos que desarrollamos y juega un papel preponderante en nuestra obtención de estimaciones de error.