Fórmulas de derivación numérica para el operador de difusión angular de Fokker-Planck

Las ecuaciones cinéticas del transporte son bien conocidas por los ingenieros nucleares. En una de ellas, la llamada ecuación de Fokker-Planck, aparece un operador diferencial de difusión, similar al de la ecuación del calor, para el cual existe, desde 1985, una discretización ciertamente curiosa, puesto que emplea nodos y pesos gaussianos, más ligados a la cuadratura que a la diferenciación. En esta charla estableceremos un marco que permite interpretar esa fórmula de derivación numérica como la perturbación de una fórmula estándar, y daremos resultados que demuestran orden 2 de convergencia.

Trabajo realizado en colaboración de los profesores Javier Segura (Universidad de Cantabria, España) y Barry Ganapol (University of Arizona, USA).