Trabajo de grado - Pregrado
Implementación de una interfaz gráfica para la simulación de plasmas electrostáticos y su aplicación en el estudio de ondas solitarias observadas en al región aural y otras regiones de la magnetosfera
Registro en:
Universidad Tecnológica de Pereira
Repositorio UTP
Autor
Osorio Tascón, Juan Camilo
Gómez Tabares, Oscar Julián
Institución
Resumen
Este trabajo de grado cuenta con un modelo teórico desarrollado, modelo computacional, interfaz gráfica y simulaciones de plasma electrostático.
dentro de sus anexos se encuentran los códigos desarrollados para este trabajo los cuales pueden ser consultado en el sitio :
Github://gitfront.io/r/user-1533735/bc0c57d4d1f428ba39910dfde17342a16c0dd00c/PEWAV/ Para este trabajo se construyó una interfaz gráfca a partir del lenguaje de programación Python para simulación de plasmas electrostáticos tomando como base el código OneGPLA.py desarrollado anteriormente en el grupo de investigación GPLA. Esto se hizo con la fínalidad de brindar una herramienta académica para el estudio de plasmas bajo condiciones que no son fácilmente reproducibles ya sea por la falta de los instrumentos o infraestructura necesaria para ello, como por su poca probabilidad de obtenerla aquí en la tierra, por ejemplo, el caso del plasma astrofísico. Es así que al implementar la interfaz gráfca: Plasma Electrostatic Wave (PEWAV), el cual utiliza el algoritmo enfocado en la dinámica de múltiples partículas llamado Particle-in-Cell (PIC), se decidió aplicarlo en las inestabilidades electrónicas del plasma para corroborar que bajo ciertas condiciones dichas inestabilidades son el mecanismo principal de generación de las ondas solitarias electrostáticas (OSE) que han sido observadas en algunas regiones de la magnetosfera como la zona de media y alta altitud. Pregrado Ingeniero(a) Físico(a) Índice general
Índice general 3
Índice de tablas 5
Índice de guras 6
Resumen y Agradecimientos 8
1. Introducción 9
1.1. Objetivos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.1.1. Objetivos Generales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.1.2. Objetivos especí cos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2. Modelo Teórico 11
2.1. Magnetosfera . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.2. Teoría Cinética del Plasma . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.3. Longitud de Debye . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.4. Frecuencia del Plasma . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.5. Velocidad Térmica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.6. Amortiguamiento de Landau . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.7. Inestabilidades en el Plasma . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.7.1. Inestabilidad Two-Stream . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.7.2. Inestabilidad Bump-of-Tail . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.7.3. Inestabilidad Beam-Plasma . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.8. Teoría de Fluido del Plasma . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.8.1. Ecuación de Continuidad . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.8.2. Ecuación de Momento . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.8.3. Ecuación de estado . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.9. Ondas Solitarias y la Ecuación Korteweg-de Vries . . . . . . . . . . . . . . . . . 19
3. Modelo Computacional 25
3.1. Método Particle-in-Cell . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
3.1.1. Método de ponderación . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
3.1.2. Integración de las ecuaciones de campo . . . . . . . . . . . . . . . . . . . 26
3.1.3. Integración de las ecuaciones de movimiento . . . . . . . . . . . . . . . . 27
3.1.4. Condiciones de frontera . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.1.5. Normalización . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.1.6. Condiciones iniciales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
3
4 ÍNDICE GENERAL
4. Interfaz Grá ca PEWAV 30
4.1. Librerías Usadas de Python . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
4.1.1. Numpy . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
4.1.2. Pandas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
4.1.3. Tkinter . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
4.1.4. Matplotlib . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
4.1.5. H5py . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
4.2. Adaptación de los códigos OneGPLA . . . . . . . . . . . . . . . . . . . . . . . . 33
4.3. Composición . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
4.3.1. ROneGPLAfunctions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
4.3.2. ROneGPLA_AIO . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
4.3.3. PEWAV . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
4.4. Manual de Usuario . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
4.4.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
4.4.2. Ejecución de la Interfaz: Introducción al Entorno . . . . . . . . . . . . . 37
4.4.3. Elementos de la Interfaz : Condiciones Iniciales . . . . . . . . . . . . . . 38
4.4.4. Elementos de la Interfaz: Grá cas . . . . . . . . . . . . . . . . . . . . . . 40
5. Simulaciones de plasmas electrostáticos en PEWAV 42
5.1. Parámetros . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
6. Resultados 45
6.1. Dos especies de electrones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
6.1.1. Inestabilidad two-stream . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
6.1.2. Inestabilidad bump-on-tail . . . . . . . . . . . . . . . . . . . . . . . . . . 46
6.1.3. Inestabilidad warm two-stream . . . . . . . . . . . . . . . . . . . . . . . 46
6.1.4. Grá cas 2S . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
6.2. Tres especies de electrones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
6.2.1. Inestabilidad two-stream . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
6.2.2. Inestabilidad weak-beam . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
6.2.3. Inestabilidad bump-on-tail . . . . . . . . . . . . . . . . . . . . . . . . . . 56
6.2.4. Grá cas 3S . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
7. Conclusiones, Discusiones y Trabajo Futuro 69
7.1. Conclusiones y Discusiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
7.2. Trabajo Futuro . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
8. Anexos 72
Bibliografía 73