En este trabajo se ha realizado una simulación numérica detallada del fenómeno de cambio de fase sólido-líquido, por ser esta fenomenología de gran interés en diferentes áreas industriales. La simulación realizada implica problemas de no linealidad, fuertes acoplamientos y frontera móvil. Como resultado de esto, solo para las configuraciones más simples se pueden utilizar herramientas analíticas, mientras que para resolver la mayoría de problemas de interés se requiere el uso de métodos numéricos. Estos métodos consisten en discretizar las ecuaciones que definen la fenomenología que nos ocupa en pequeñas celdas o volúmenes de control. En este estudio se ha optado por utilizar el Método de Volúmenes Finitos (FVM) para la discretización de las ecuaciones gobernantes utilizando mallas cartesianas. Se utiliza una malla desplazada; esto quiere decir que las componentes del vector velocidad se calculan en las caras de los volúmenes de control, lo que permite un acoplamiento adecuado entre la ecuación de conservación de la masa y momentum. Tanto el código como las soluciones numéricas han sido convenientemente verificados. La verificación del código consiste en comprobar que éste está libre de errores de programación y que el comportamiento de los esquemas numéricos implementados está acorde con su comportamiento teórico. Para la verificación de la solución numérica se han utilizado métodos de extrapolación de Richardson o realizando un estudio de refinamiento de malla y observando la evolución de algunas magnitudes características del problema como pueden ser la fracción de líquido en el dominio o el número de Nusselt en la pared caliente.Una vez que el código y las soluciones numéricas han sido convenientemente verificados, la validación final de la simulación del proceso es la que se obtiene de comparar los resultados predecidos con los datos experimentales.El problema de la fusión del galio en una cavidad rectangular calentada por un lado ha sido ampliamente utilizado por investigadores con el fin de evaluar los métodos numéricos para la resolución del cambio de fase. Si bien comentar que este material tiene las ventajas de que sus propiedades termofísicas están bien establecidas, que tiene una temperatura de cambio de fase cercana a la temperatura ambiente y que es un material con gran interés industrial, también presenta algunas desventajas como es el hecho de que tiene un comportamiento anisotrópico en cuanto a la conductividad térmica de la fase sólida. No obstante, el gran número de trabajos experimentales que se encuentran en la literatura relacionados con este material nos ha conducido a escoger este material para la realización de un estudio detallado del cambio de fase sólido-líquido en este trabajo.Si bien en la literatura aparece este problema con diferentes configuraciones, en este trabajo nos hemos centrado en el estudio del caso en que la relación de aspecto (alto/ancho) es de 0,5.Existen diferentes métodos numéricos pararesolver problemas de cambio de fase sólido líquido: métodos que siguen la frontera móvil, métodos que fijan la frontera móvil, etc. En este trabajo utilizaremos el Método Entálpico pues nos permite utilizar una malla fija en todo el dominio, la condición de Stefan queda impuesta de manera implícita, permite la coexistencia de más de un frente de cambio de fase y permite que la interfase tenga un cierto grosor. El caso que nos ocupa tiene la singularidad de encontrarse en un rango de Prandtl muy bajo. Esto provoca que se alcance el régimen turbulento para números de Rayleigh relativamente bajos. Esto nos ha hecho pensar sobre la conveniencia de realizar un estudio más detallado para determinar para que valor de Rayleigh se produce la transición de régimen permanente a flujo oscilatorio y de éste a un régimen caótico. Los problemas que han ido surgiendo a lo largo de la realización de este trabajo nos han conducido a la utilización del método multibloc, también conocido como método de descomposición de subdominios. El método se emplea considerando flujos incompresibles y mallas desplazadas. En este trabajo se explicarán las modificaciones que han sido necesarias para la utilización de este método en la fenomenología del cambio de fase sólido líquido. Se han utilizado dos aproximaciones: una conservativa y otra basada en la presión. Para la obtención de la solución de referencia se ha empleado los métodos de subdominios basados en la presión, pues con el método conservativo se han observado discrepancias entre la solución obtenida con un único subdominio y la obtenida con varios subdominios. Finalmente, se ha realizado un estudio parámetrico del caso, para el cual se han considerado diferentes relaciones de aspecto, diferentes condiciones de contorno y variaciones de ±10% en las propiedades termofísicas con respecto al caso de referencia, con el objetivo de ver como afectan estas modificaciones sobre la fenomenología que nos ocupa. / In this thesis a detailed numerical simulation of liquid-solid phase change phenomena has been made, because this phenomenology is of great interest in different industrial areas. The simulation done implies problems of nonlinearity, strong couplings and movable interphase. Like a result, only for the simplest configurations analytical tools can be used, whereas to solve the most of interest problems numerical methods are needed. These methods consist in the discretisation of the equations that define the phenomenology in small cells or control volumes. In this study the Finite Volume Method (FVM) has been used for the governing equations discretisation using Cartesian meshes. A displaced mesh is used; this means that the components of the speed vector are calculated in the faces of the control volumes, which allows a correct coupling between the continuity and momentum equations. So the code as the numerical solutions have been properly verified. The code verification consist of verifying that this is free of programming errors and that the behaviour of the numerical schemes is agreed with the theoretical one. For the verification of the numerical solution the Richardson Extrapolation Method or a mesh refinement study have been used. Once the code and the numerical solutions have been properly verified, the final validation of the process simulation is obtained comparing the numerical results with experimental ones.The Gallium melting problem in a square cavity heated by a side has been widely used by investigators with the objective of evaluatingthe numerical methods used for solving the phase change phenomena. Although this material have the advantages of its thermophysical properties are well established, the phase change temperature is near the room temperature and is a material with a great industrial interest, also presents some disadvantages like having an anisotropic thermal conductivity of the solid phase. However, the great number of experimental works that can been found in the literature, lead us to choose this material for doing a detailed study of the solid liquid phase change in this work. Although in the literature this problem appears with different configurations, in this thesis we have centered in the study of the case with an aspect ratio (height/width) of 0.5. Different numerical methods exist to solve solid-liquid phase change problems: methods that follow the moving interphase, methods that fix the moving interphase, etc. In this work we have used the Enthalpy Method because it allows us to use a fixed mesh in all the domain, the Stefan condition is imposed implicitly, it allows the coexistence of more than one front of phase change and allows the interphase has a thickness.This case has the singularity of being in a low Prandtl number range. This causes that a turbulent state has been reached for relatively low Rayleigh numbers. This has made think us on the convenience of making a detailed study to determine the transition of permanent regime to oscillating flow and from oscillating flow to chaotic regime. The problems than have been found in the accomplishment of this work have lead us to the use of multiblock method. This method is used considering incompressible flows and moved meshes. We will explain the modifications that have been necessary for using this method in the phase change phenomena. Two approaches have been used: a conservative one and a pressure based one. For obtaining the reference solution has been use the pressure based method, because the conservative method presents discrepancies between the obtained solution with an only subdomain and the obtained with several subdomains. Finally, a parametric study has been done, considering different aspect ratios, boundary conditions and variations of ±10% in the thermophysical properties with respect to the reference case, with the objective to see as these changes affect on the phenomenology that occupies to us.

In recent years, in parallel with advances in computer technology, a multitude of computer tools have been developed through which it is possible to obtain a detailed description of the phenomena occurring in the core of nuclear reactors. The final ob-jective of these new tools is to perform safety analysis using best estimate techniques. The best estimate techniques, as opposed to the conservative ones, allow the operation of the reactor with narrower safety margins, and thus greater core economy. In this context, in this work is developed an multiphysics computer platform that inte-grates simulation codes that cover most of the physics that take place in nuclear reac-tors. For the integration of the different feedback phenomena between thermal-hydraulics, neutronics and heat transfer, a series of couplings have been developed between the codes that compose the platform. All the developments carried out are intended to realistically represent the design and behavior of the nuclear facility, in-cluding the control system, fuel elements and fuel rods. The computer platform includes some of the state-of-the-art codes for reactor analysis. The thermal-hydraulics is covered with a developed coupled code, consisting of the semi-implicit coupling between the TRACE system code and the subchannel code COBRA-TF (CTF), whose parallel version has been created in this work. In transients where three-dimensional neutron calculations are necessary, the explicit coupling be-tween the three-dimensional PARCS core simulator and the subchannel code CTF has been developed. For the analysis of the integrity of the fuel rods, the FRAPCON and FRAPTRAN codes are used, coupling the latter explicitly with CTF. All the developed tools have been included in the same computer platform that en-compasses and coordinates the simulations under the user's guidelines. The platform has enough flexibility to perform safety studies in a multitude of operational or acci-dental scenarios, and it is hoped that in the future it may be used for supporting li-cense calculations. The developed tools have been verified through a series of practical applications in different transient and accidental scenarios in light water reactors. The results obtained have been compared with actual plant measurements and with the results of other simulation codes showing adequate predictive capacity. The work carried out in this doctoral thesis is part of the research line financed by the Ministerio de Economía y Competitividad in the NUC-MULTPHYS project (ENE2012-34585) and the interdisciplinary collaboration projects of the Universitat Politècnica de Valencia COBRA_PAR (PAID-05-11-2810) and Open-NUC (PAID-05-12). / En los últimos años, paralelamente a los avances en tecnología informática, se están desarrollando multitud de herramientas informáticas mediante las que es posible obte-ner una descripción detallada de los fenómenos que tienen lugar en el núcleo de los reactores nucleares. El objeto de estas nuevas herramientas es el de realizar análisis de seguridad en reactores nucleares utilizando técnicas de mejor estimación. Las técnicas de mejor estimación, en contraposición con las conservadoras, permiten la operación del reactor con márgenes de seguridad más estrechos, y por tanto mayor economía del núcleo. En este contexto, en la presente tesis doctoral se desarrolla una plataforma informática que integra códigos informáticos que cubren la mayor parte de las físicas que tienen lugar en los reactores nucleares. Para la integración de los diferentes fenómenos de realimentación entre termohidráulica, neutrónica, mecánica y transmisión de calor se han desarrollado una serie de acoplamientos entre los códigos que componen la plata-forma. Todos los desarrollos realizados tienen por objetivo representar de forma rea-lista el diseño y comportamiento de la instalación nuclear, incluyendo el sistema de control, los elementos y las varillas de combustible. En la plataforma informática se incluyen algunos de los códigos de última generación (estado de arte) para el análisis del comportamiento de reactor. En el plano termohi-dráulico se utiliza el código acoplado desarrollado, formado por el acople semi-implícito entre el código de sistema TRACE y el de subcanal COBRA-TF (CTF), cuya versión paralela ha sido creada en este trabajo. En transitorios en los que resultan ne-cesarios los cálculos de neutrónica tridimensional, se ha desarrollado el acople explíci-to entre el simulador tridimensional de núcleos PARCS y el código de subcanal CTF. Para el análisis de la integridad de las varillas de combustible se emplean los códigos FRAPCON y FRAPTRAN, acoplando este último de forma temporalmente explícita con CTF. Todos los desarrollos realizados se han incluido en una misma plataforma informática que los engloba y coordina las simulaciones bajo las directrices del usuario. La plata-forma posee suficiente flexibilidad para realizar estudios de seguridad en multitud de escenarios operacionales o accidentales, y se desea que en un futuro pueda ser utilizada en cálculos de apoyo a licencia. Las herramientas desarrolladas han sido verificadas mediante una serie de aplicaciones prácticas en distintos transitorios y escenarios acci-dentales en reactores de agua ligera. Los resultados obtenidos se han comparado con medidas reales de planta y con los resultados de otros códigos de simulación mostran-do una adecuada capacidad predictiva. El trabajo realizado en la presente tesis doctoral se enmarca dentro de la línea de in-vestigación financiada por el Ministerio de Economía y Competitividad en el proyec-to NUC-MULTPHYS (ENE2012-34585) y los proyectos de colaboración interdisci-plinar de la Universitat Politècnica de Valencia COBRA_PAR (PAID-05-11-2810) y Open-NUC (PAID-05-12) / En els últims anys, paral·lelament als avanços en tecnologia informàtica, s'estan desenvolupant multitud de ferramentes informàtiques mitjançant les quals és possible obtindre una descripció detallada dels fenòmens que tenen lloc en el nucli dels reactors nuclears. L'objecte final d'aquestes noves ferramentes és el de realitzar anàlisis de segu-retat a reactors nuclears utilitzant tècniques de millor estimació. Les tècniques de mi-llor estimació, en contraposició amb les conservadores, permeten l'operació del reactor amb marges de seguretat més estrets, i per tant una major economia del nucli. En aquest context, en el present treball de tesi es desenvolupa una plataforma in-formàtica que integra codis informàtics que cobreixen la major part de les físiques que tenen lloc als reactors nuclears. Per a la integració dels diferents fenòmens de reali-mentació entre termohidràulica, neutrònica i transmissió de calor s'han desenvolupat una sèrie d'acoblaments entre els codis que componen la plataforma. Tots els desenvo-lupaments realitzats tenen per objectiu representar de forma realista el disseny i com-portament de la instal·lació nuclear, incloent el sistema de control, els elements i les varetes de combustible. A la plataforma informàtica s'inclouen alguns dels codis d'última generació (estat de l'art) per a l'anàlisi del comportament de reactor. En el pla termohidràulic s'utilitza el codi acoblat desenvolupat, format per l'acoblament semi-implícit entre el codi de sis-tema TRACE i el de subcanal COBRA-TF (CTF), en una versió paral·lela creada en aquest treball. En transitoris en els que resulten necessaris els càlculs de neutrònica tridimensional, s'ha desenvolupat l'acoblament explícit entre el simulador tridimensio-nal de nuclis PARCS i el codi de subcanal CTF. Per a l'anàlisi de la integritat de les varetes de combustible s'empren els codis FRAPCON i FRAPTRAN, acoblant aquest últim de forma temporalment explícita amb CTF. Tots els desenvolupaments realitzats s'han inclòs en una mateixa plataforma informàti-ca que els engloba i coordina les simulacions sota les directrius de l'usuari. La plata-forma posseeix suficient flexibilitat per realitzar estudis de seguretat en multitud d'es-cenaris operacionals o accidentals, i es desitja que en un futur pugui ser utilitzada en càlculs de suport a llicència. Les ferramentes desenvolupades han sigut verificades mitjançant una sèrie d'aplicacions pràctiques en diferents transitoris i escenaris acci-dentals en reactors d'aigua lleugera. Els resultats obtinguts s'han comparat amb mesu-res reals de planta i amb els resultats obtinguts amb altres codis de simulació, mostrant una adequada capacitat predictiva. El treball realitzat en la present tesi doctoral s'emmarca dins de la línia d'investigació finançada pel Ministeri d'Economia i Competitivitat en el projecte NUC-MULTPHYS (ENE2012-34585) i els projectes de col·laboració interdisciplinar de la Universitat Politècnica de València COBRA_PAR (PAID-2810.11.05) i Open-NUC (PAID-05-12). / Abarca Giménez, A. (2017). Desarrollo y verificación de una plataforma multifísica de altas prestaciones para análisis de seguridad en ingeniería nuclear [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/88399

No description available.