Preguntas frecuentes sobre ZIVIS (Blog de Zivis)

Ya podemos enseñaros los primeros resultados de Zivis. Los hemos obtenido de
las trayectorias que habéis simulado en vuestros ordenadores correspondientes
al plasma A0.

Se trata de un plasma ECRH (resonancia electrón-ciclotrón). Brevemente, este
plasma se calienta mediante radiación microondas cuya energía absorben los
electrones y se la comunican a los iones (chocando con ellos). Están
caracterizados por una temperatura electónica alta, temperatura iónica baja y
densidad también baja.

Como ya se ha comentado en este blog, A0 corresponde a "apagar" el campo
eléctrico del TJ-II. De este modo, comparando A1 y A0, uno estima cuál es el
efecto del campo eléctrico en el confinamiento. Lo de apagar el campo
eléctrico es algo que nos podemos permitir hacer en la simulación, claro, en
la realidad, siempre está.



Lo que veis son, para A0 y vistos desde arriba, los puntos de choque con la
cámara de vacío de los iones que habéis simulado. De esa imagen pueden
sacarse dos tipos de información relevantes:

- Qué partes de la cámara están sometidas a más impactos y podrían resultar
dañadas.

- Dónde se van a generar más impurezas en el plasma: los iones que golpean las
paredes pueden producir que éstas emitan partículas al plasma. Esto es
perjudicial porque estas impurezas enfrían el plasma.


Uno puede ver varios patrones en los puntos de choque. Por ejemplo, esa serie
de bandas radiales tiene una explicación sencilla: todas las trayectorias
que calculáis, antes o después, acaban como la que veis en la imagen:



Después de dar muchas vueltas (típicamente entre 100 y 1000!) se han ido
frenando y se quedan atrapadas en un llamado "pozo magnético", entre dos
bobinas, donde el campo magnético es más intenso. Ahí empiezan a oscilar a la
vez que "derivan" hacia fuera, hasta que se encuentran una pared de la cámara
de vacío, como veis más abajo.



La idea, por cierto, es que esto ocurra cuanto más tarde mejor.

El siguiente paso será estudiar la diferencia entre A0 y A1.


JLVG (del equipo de física de Zivis)

Calculados con lgv, a eso de las cinco y media del sabado 14.

118345 ficheros A0, que hacen 1183450 trayectorias

120099 ficheros A1, que hacen 1200990 trayectorias

9812 ficheros B0, que harian 98120

18 B1, que hacen 180 y otros

61231 ficheros B1 cortas, que hacen 61231 trayectorias.

Calculados con lgv2 (quizas haya que descartarlos)

1839 ficheros B0, que hacen 18390 trayectorias

16318 ficheros B1cortas, que hacen 16318 solo, claro

calculados con lgv3, el que hemos lanzado esta mañana.

279 B0, que hacen 2790

2626 B1, que hacen, claro, 2626

A ver como se nos da el dia.

Etiquetas: fisica, tecnica

Segun me dicen los de la parte fisica, la diferencia entre las tareas x0 y x1 es que en las 1 se enciende un campo electrico que ayuda al confinamiento. Luego cuando lo entienda con mas detalle os lo cuento. Pero el resultado es que B1 genera confinamientos de orden de decimas de segundo. En los batch que hemos estado mandando se calculaban 10 trayectorias, y se ve que un segundo de simulacion es mucha CPU.

Ahora estamos mandando una serie nueva de B1 con solo una trayectoria y solamente con recogida de estadistica, sin escribirla en el disco. De acuerdo a lo que tarden veremos si saltamos a a otra configuracion o si nos atrevemos con B1 a 10 trayectorias en las maquinas mas potentes.

ADDENDA:

El responsable cientifico me ha contado algo mas. Una diferencia entre A y B es que la temperatura ionica de B es mas baja, asi que uno podria esperar que tanto B1 como B0 dieran mas tiempo de confinamiento que A1 y A0. Por otro lado no es que en las "1" se encienda algo en el ordenador, sino en el software: se activa un flag que hace que las trayectorias tomen en consideracion el campo electrostatico que produce el propio plasma, y que se toma de datos medidos experimentalmente en el TJIII (el dispositivo que estamos simulando). Al tomar en consideracion este campo, los resultados son más realistas. Y si el potencial genera un buen pozo, pues el tiempo de confinamiento puede ser facilmente un orden de magnitud superior, y por tanto bastante más en los casos extremos.

Lo que haremos hoy martes y mañana miercoles es generar trayectorias B1 una por una, de forma que no saturen las maquinas, y cuando tengamos unas cuantas podremos tomar decisiones acerca de como estudiar esta configuracion. Si a final de la tarde vemos que las maquinas pequeñas no pueden con ellas, seguiremos mandandoles A0 y A1 a estas para mejorar la estadistica de la configuracion A.

Las simulaciones nuevas, de una sola trayectoria, se llaman "lgv_B1_cortas_". A ver si con esas hay menos errores. El campo "progreso" saltará de 0% a 100% de repente, porque lo actualizabamos al final de cada trayectoria.

Por cierto que ahora me estan entrando ganas de añadir dibujo al programa para que se vea por donde viaja la particula hasta que choca.