Descripción del proyecto
EL PROYECTO CONSITE EN EL DESARROLLO DE UN CODIGO DE TRANSPORTE CON FVM Y OD PARA PROBLEMAS DE AUTOVALOR CON FUENTE DE FISION Y DE FUENTE FIJA, BASANDONOS EN UNA LIBRERIA DE DISCRETIZACION DE ECUACIONES EN DERIVADAS PARCIALES POR EL METODO DE VOLUMENES FINITOS (FVM), SE PROGRAMARA LA PARTE GEOMETRICA DE LAS ECUACIONES DE TRANSPORTE NEUTRONICO EN ESTADO ESTACIONARIO Y EN 3 DIMENSIONES, PARA EL CALCULO DE AUTOVALOR CON FUENTE DE FISION LAS LIBRERIAS DE FVM DEBEN DE MODIFICARSE, PUESTO QUE ESTAN PREPARADAS PARA LA RESOLUCION DE PROBLEMAS DE FUENTE FIJA O FLUJO, RESOLVIENDO UN SISTEMA DE ECUACIONES LINEALES APOYANDOSE EN LIBRERIAS COMO MUMPS O PETSC, PARA LOS PROBLEMAS DE CRITICIDAD, EN LOS QUE SE DEBE DE RESOLVER UN PROBLEMA DE AUTOVALORES/AUTOVECTORES, SE DEBEN DE MODIFICAR ESTAS LIBRERIAS PARA RESTRUCTURAR EL PLANTEAMIENTO DEL PROBLEMA DE TAL FORMA QUE QUEDE COMO UN PROBLEMA DE AUTOVALORES Y ADEMAS INCLUIR EN SUS METODOS DE CALCULO ALGUNA LIBRERIA NUMERICA QUE RESUELVA ESTE TIPO DE PROBLEMAS, COMO SLEPC, LA PARTE ANGULAR DEL PROBLEMA SE RESOLVERA POR EL METODO DE ORDENADAS DISCRETAS (O METODO SN) UTILIZANDO UNA CUADRATURA DE ORDEN 16, EL PROBLEMA SE PLANTEARA DISCRETIZANDO TAMBIEN LA VARIABLE ENERGETICA CON LA APROXIMACION MULTIGRUPO CON UN NUMERO GENERICO DE GRUPOS Y TENIENDO EN CUENTA NO SOLO LOS TERMINOS DE DOWNSCATERING SINO TAMBIEN LOS DE UPSCATERING, LOS FICHEROS DE ENTRADA DE LA PARTE GEOMETRICA SE REALIZARAN CON UN MALLADOR DE LIBRE DISTRIBUCION Y CODIGO LIBRE COMO GMSH, QUE PERMITE TANTO GENERAR MALLAS ESTRUCTURADAS COMO NO ESTRUCTURADAS, LAS CONDICIONES DE CONTORNO Y OPCIONES NUMERICAS SE ESPECIFICARAN EN UN FICHERO DE ENTRADA EN FORMATO ASCII,LOS CODIGOS DE TRANSPORTE GENERICOS EN GEOMETRIAS 3D, COMO EL QUE SE PRETENDE DESARROLLAR, SON GRANDES CONSUMIDORES DE RECURSOS COMPUTACIONALES, SE PROPONE ADEMAS LA PARALELIZACION TANTO DE LA PARTE ESPACIAL COMO DE LA ANGULAR, TANTO A NIVEL DE SOLVER (PETSC/SLEPC), COMO A NIVEL DE SUBDOMINIO CON LA LIBRERIA MPI, ESTE PUNTO ES CRUCIAL PARA EL SALTO A SU USO PRODUCCION DEL CODIGO, LA PARALELIZACION DE LA PARTE ESPACIAL YA SE HA PROBADO EN EL DESARROLLO DE UN CODIGO DE DIFUSION NEUTRONICA 3D POR EL METODO DE VOLUMENES FINITOS (FVM), LA PARTE ANGULAR TAMBIEN SE PUEDE DESCOMPONER EN VARIAS TAREAS (PARES) PARA ENVIAR A DIFERENTES PROCESADORES,TODO PROBLEMA DE AUTOVALORES CONLLEVA ASOCIADO UN PROBLEMA ADJUNTO (CONJUGADO TRASPUESTO), EL PROBLEMA DIRECTO, POR CONSIDERACIONES FISICAS, TIENE SOLUCION REAL, POR LO QUE SU PROBLEMA ADJUNTO CONSISTE SIMPLEMENTE EN LA TRASPOSICION DEL PROBLEMA DIRECTO, EL CUAL TAMBIEN SE PROGRAMARA, TIENE IMPORANTANCIA PARA EL CALCULO DEL MODULO TRANSITORIO,LA RESOLUCION DE LA PARTE ANGULAR DE LA ECUACION DEL TRANSPORTE POR ORDENADAS DISCRETAS (SN) NO PERMITE LA OBTENCION DE LOS MODOS (ARMONICOS NATURALES), ASI SE IMPLEMENTARA TAMBIEN UN METODO BASADO EN ARMONICOS ESFERICOS (PN) O HIBRIDO PN-SN, CUYA DISCRETIZACION SI PERMITE LA OBTENCION DE MODOS,ADEMAS, TAMBIEN SE RESOLVERA LA ECUACION DEL TRANSPORTE EN MULTIGRUPO Y CON PRECURSORES DE NEUTRONES DIFERIDOS MEDIANTE UN METODO CUASI-ESTATICO VARIACIONAL (PERTURBACIONES) PESADO CON EL FLUJO ADJUNTO EL CUAL PERMITE LA DESCOMPOSICION DEL FLUJO NEUTRONICO EN UNA PARTE ESPACIAL Y OTRA TEMPORAL (AMPLITUD), LA PARTE ESPACIAL CAMBIA MUY POCO (DEBIDO A LOS CAMBIOS DE LAS SECCIONES EFICACES) EN CADA PASO DE TIEMPO, ASI Y TODO SE IMPLEMENTARA UN METODO DE ACTUALIZACION (CADA N PASOS DE TIEMPO) PARA LA SOLUCION ESPACIAL, TRANSPORTE RADIACIÓN\ANÁLISIS MODAL\VOLÚMENES FINITOS\ORDENADAS DISCRETAS\SEGURIDAD NUCLEAR\CÁLCULO DE DOSIS