Optimización de algoritmos de resolución de la ecuación de transporte con orientación a su paralelización.
La física de reactores está basada en la resolución de la ecuación de transporte de neutrones, para la cual existen diversos esquemas numéricos y códigos de cálculo que los implementan. La programación en paralelo permite aumentar la velocidad de estos códigos mediante la utilización de múltiples...
Guardado en:
| Autor principal: | |
|---|---|
| Formato: | Tesis NonPeerReviewed |
| Lenguaje: | Español |
| Publicado: |
2015
|
| Materias: | |
| Acceso en línea: | http://ricabib.cab.cnea.gov.ar/504/1/1Garc%C3%ADa_Manuel.pdf |
| Aporte de: |
| Sumario: | La física de reactores está basada en la resolución de la ecuación de transporte de
neutrones, para la cual existen diversos esquemas numéricos y códigos de cálculo que
los implementan. La programación en paralelo permite aumentar la velocidad de estos
códigos mediante la utilización de múltiples unidades de procesamiento, mejorando
la capacidad de modelado. En este trabajo se estudió la paralelización de distintos
algoritmos asociados a la ecuación de transporte.
Para la etapa de celda se consideraron los métodos de probabilidades de colisión y
de respuesta heterogenea (HRM). El cálculo de
flujos de respuesta en el primer método
fue paralelizado por grupos de energía, y en el segundo por grupos y elementos. Para el
método multigrupo correspondiente a HRM se desarrolló un algoritmo de paralelización
de las iteraciones térmicas y de las operaciones matriciales asociadas al esquema ite-
rativo. También se obtuvieron algoritmos paralelos para los cálculos de quemado. Los
métodos paralelizados en esta etapa fueron implementados en el código CONDOR,
desarrollado en INVAP.
El método de difusión en diferencias finitas fue analizado para los cálculos de núcleo.
La paralelización se realizó sobre el método de Gauss-Seidel, con un esquema rojo-negro
por filas. Este algoritmo fue implementado en el código CITVAP, basado en CITATION
2.0 de Oak Ridge.
La programación de los algoritmos desarrollados se realizó en OpenMP, una herra-
mienta de procesamiento en paralelo en el modelo de memoria compartida. También
se consideró un modelo híbrido MPI-OpenMP para el cálculo de
flujos de respuesta en
HRM. |
|---|