Métodos numéricos integrales para ecuaciones Fokker-Planck

El trabajo presenta un método numérico de resolucion para ecuaciones de fokker-planck no lineales. El modelo mejora el tratamiento de problemas cinéticos con mas de una cantidad conservada, subsanando las deficiencias que, en este sentido, presentan los esquemas habituales en diferencias. El método...

ver descrição completa

Detalhes bibliográficos
Autor: Donoso Vargas, José Manuel
Formato: tesis doctoral
Fecha de publicación:2002
País:España
Recursos:Universidad Complutense de Madrid (UCM)
Repositorio:Docta Complutense
Idioma:español
OAI Identifier:oai:docta.ucm.es:20.500.14352/62848
Acesso em linha:https://hdl.handle.net/20.500.14352/62848
Access Level:acceso abierto
Palavra-chave:Física nuclear
2207 Física Atómica y Nuclear
Descrição
Resumo:El trabajo presenta un método numérico de resolucion para ecuaciones de fokker-planck no lineales. El modelo mejora el tratamiento de problemas cinéticos con mas de una cantidad conservada, subsanando las deficiencias que, en este sentido, presentan los esquemas habituales en diferencias. El método se basa en la ecuación integral de evolución para la distribución a través de la probabilidad de transición o propagador. En el segundo capitulo se propone un sencillo algoritmo para la obtención del propagador a tiempos cortos basado la construcción de una ecuación auxiliar local de fokker-plack con solución conocida. La distribución se implementa en el tiempo mediante una matriz de evolución con elementos positivos y acotados, que conserva la norma unidad y la positividad de . El ajuste del propagador se consuma durante el ejercicio de la computación, con el uso recursivo de correcciones a segundo orden de potencias en el paso temporal. El método se ilustra el capitulo tercero, en el que se resuelven ecuaciones simples. En los dos capítulos finales el método se aplica a la ecuación de fokker-planck de la física del plasma, para simetrías esférica y cilíndrica en el espacio de velocidades. En ambos casos se dan las expresiones de los propagadores a tiempos cortos. Estos propagadores configuran un eficiente operador colisional discreto, de propiedades conservativas consistentes con el problema continuo, con independencia absoluta en el número de iteraciones.