Design options have been evaluated for the Modular Helium Reactor (MHR) for higher temperature operation. An alternative configuration for the MHR coolant inlet flow path is developed to reduce the peak vessel temperature (PVT). The coolant inlet path is shifted from the annular path between reactor core barrel and vessel wall through the permanent side reflector (PSR). The number and dimensions of coolant holes are varied to optimize the pressure drop, the inlet velocity, and the percentage of graphite removed from the PSR to create this inlet path. With the removal of ~10% of the graphite from PSR the PVT is reduced from 541 0C to 421 0C. A new design for the graphite block core has been evaluated and optimized to reduce the inlet coolant temperature with the aim of further reduction of PVT. The dimensions and number of fuel rods and coolant holes, and the triangular pitch have been changed and optimized. Different packing fractions for the new core design have been used to conserve the number of fuel particles. Thermal properties for the fuel elements are calculated and incorporated into these analyses. The inlet temperature, mass flow and bypass flow are optimized to limit the peak fuel temperature (PFT) within an acceptable range. Using both of these modifications together, the PVT is reduced to ~350 0C while keeping the outlet temperature at 950 0C and maintaining the PFT within acceptable limits. The vessel and fuel temperatures during low pressure conduction cooldown and high pressure conduction cooldown transients are found to be well below the design limits. The reliability and availability studies for coupled nuclear hydrogen production processes based on the sulfur iodine thermochemical process and high temperature electrolysis process have been accomplished. The fault tree models for both these processes are developed. Using information obtained on system configuration, component failure probability, component repair time and system operating modes and conditions, the system reliability and availability are assessed. Required redundancies are made to improve system reliability and to optimize the plant design for economic performance. The failure rates and outage factors of both processes are found to be well below the maximum acceptable range.

In order to protect gas turbine blades from hot gases exiting the combustor, several intricate external and internal cooling concepts are employed. High pressure stage gas turbine blades feature serpentine passages where rib turbulators are installed to enhance heat transfer between the relatively colder air bled off from the compressor and the hot internal walls. Most of the prior studies have been restricted to Reynolds number of 90000 and several studies have been carried out to determine geometrically optimized parameters for achieving high levels of heat transfer in this range of Reynolds number. However, for land-based power generation gas turbines, the Reynolds numbers are significantly high and vary between 105 and 106. Present study is targeted towards these high Reynolds numbers where traditional rib turbulator shapes and prescribed optimum geometrical parameters have been investigated experimentally. A steady-state liquid crystal thermography technique is employed for measurement of detailed heat transfer coefficient. Five different rib configurations, viz., 45 deg., V-shaped, inverse V-shaped, W-shaped and M-shaped have been investigated for Reynolds numbers ranging from 150,000 to 400,000. The ribs were installed on two opposite walls of a straight duct with aspect ratio of unity. For very high Reynolds numbers, the heat transfer enhancement levels for different rib shapes varied between 1.3 and 1.7 and the thermal hydraulic performance was found to be less than unity. / Master of Science

The modelling of multi-phase water flow is an important modern-day design tool used by engineers to develop practical systems which are beneficial to society . Multi-phase water flow can be found in many important industrial applications such as large scale conventional and nuclear power systems, heat transfer machinery, chemical process plants, and other important examples. Because of many inherent complexities in physical two-phase flow processes, no generalised system of equations has been formulated that can accurately describe the two-phase flow of water at all flow conditions and system geometries. This has led to the development of many different models for the simulation of two-phase flow at specific conditions. These models vary greatly in complexity. The simplest model that can be used to simulate two-phase flow is termed the homogeneous equilibrium (HEM) two-phase flow model. This model has been found useful in investigations of choking and flashing flows, and as an initial investigative model used before the formulation of more complex models for specific applications. This flow model is fully de ned by three conservation equations, one each for mass, momentum and energy. To close the model, an equation of state (EOS) is required to deliver system pressure values. When solving the HEM, a general practice is to employ an equation of state that is derived from a fundamental expression of the second law of thermodynamics. This methodology has been proven to deliver accurate results for two-phase system simulations. This study focused on an alternative formulation of the equation of state which was previously developed for the time dependent modelling of HEM two-phase flow systems, termed the rate form of the equation of state (RFES). The goal of the study was not to develop a new formulation of the EOS, but rather to implement the RFES in a transient simulation model and to verify that this implementation delivers appropriate results when compared to the conventional implementation methodology. This was done by formulating a transient pipe and reservoir network model with the HEM, and closing the model using both the RFES and a benchmark EOS known to deliver accurate system property values. The results of the transient model simulations were then compared to determine whether the RFES delivered the expected results. It was found that the RFES delivered sufficiently accurate results for a variety of system transients, pressure conditions and numerical integration factors. / MIng (Nuclear Engineering), North-West University, Potchefstroom Campus, 2014

The modelling of multi-phase water flow is an important modern-day design tool used by engineers to develop practical systems which are beneficial to society . Multi-phase water flow can be found in many important industrial applications such as large scale conventional and nuclear power systems, heat transfer machinery, chemical process plants, and other important examples. Because of many inherent complexities in physical two-phase flow processes, no generalised system of equations has been formulated that can accurately describe the two-phase flow of water at all flow conditions and system geometries. This has led to the development of many different models for the simulation of two-phase flow at specific conditions. These models vary greatly in complexity. The simplest model that can be used to simulate two-phase flow is termed the homogeneous equilibrium (HEM) two-phase flow model. This model has been found useful in investigations of choking and flashing flows, and as an initial investigative model used before the formulation of more complex models for specific applications. This flow model is fully de ned by three conservation equations, one each for mass, momentum and energy. To close the model, an equation of state (EOS) is required to deliver system pressure values. When solving the HEM, a general practice is to employ an equation of state that is derived from a fundamental expression of the second law of thermodynamics. This methodology has been proven to deliver accurate results for two-phase system simulations. This study focused on an alternative formulation of the equation of state which was previously developed for the time dependent modelling of HEM two-phase flow systems, termed the rate form of the equation of state (RFES). The goal of the study was not to develop a new formulation of the EOS, but rather to implement the RFES in a transient simulation model and to verify that this implementation delivers appropriate results when compared to the conventional implementation methodology. This was done by formulating a transient pipe and reservoir network model with the HEM, and closing the model using both the RFES and a benchmark EOS known to deliver accurate system property values. The results of the transient model simulations were then compared to determine whether the RFES delivered the expected results. It was found that the RFES delivered sufficiently accurate results for a variety of system transients, pressure conditions and numerical integration factors. / MIng (Nuclear Engineering), North-West University, Potchefstroom Campus, 2014

One of the greatest challenges facing the nuclear power industry is the final disposition of nuclear waste. To meet the needs of the nuclear power industry, a new fuel assembly design, called DUPLEX, has been developed which provides higher fuel burnups, burns transuranic waste while reducing minor actinides, reduces the long term radiotoxicity of spent nuclear fuel, and was developed for use in current light water reactors. The DUPLEX design considered in this thesis is based on a seed and blanket unit (SBU) configuration, where the seed region contains standard UO2 fuel, and the blanket region contains an inert matrix (Pu,Np,Am)O2-MgO-ZrO2 fuel. The research efforts of this thesis are first to consider the higher burnup effects on DUPLEX assembly thermal-hydraulic performance and thermal safety margin over the assembly’s expected operational lifetime. In order to accomplish this, an existing burnup-dependent thermal-hydraulic methodology for conventional homogeneous fuel assemblies has been updated to meet the modeling needs specific to SBU-type assemblies. The developed framework dramatically expands the capabilities of the latest thermal-hydraulic evaluation framework such that the most promising and unique DUPLEX fuel design can be evaluated. As part of this updated methodology, the posed DUPLEX design is evaluated with respect to the minimum departure from nucleate boiling ratio, peak fuel temperatures for both regions, and the peak cladding temperatures, under ANS Condition I, II, and III transient events with the thermal-hydraulic code VIPRE-01. Due to difficulty in the fabrication and handling of minor actinide dioxides, documented thermal conductivity values for the considered IMF design are unavailable. In order to develop a representative thermal conductivity model for use in VIPRE-01, an extensive literature survey on the thermal conductivity of (Pu,Np,Am)O2-MgO-ZrO2 component materials and a comprehensive review of combinatory models was performed. Using the updated methodology, VIPRE-01 is used to perform steady-state and transient thermal hydraulic analyses for the DUPLEX fuel assembly. During loss-of-flow accident scenarios, the DUPLEX design is shown to meet imposed safety criteria. However, using the most conservative thermal conductivity modeling approach for (Pu,Np,Am)O2-MgO-ZrO2, the blanket region fuel temperatures remain only slightly below the design limit.

Two new fuel assembly designs for light water reactors using advanced mixed-oxide fuels have been proposed to reduce the radiotoxicity of used nuclear fuel discharged from nuclear power plants. The research efforts of this thesis are the first to consider the effects of burnup on advanced mixed-oxide fuel assembly performance and thermal safety margin over an assembly?s expected operational burnup lifetime. In order to accomplish this, a new burnup-dependent thermal-hydraulic analysis methodology has been developed. The new methodology models many of the effects of burnup on an assembly design by including burnup-dependent variations in fuel pin relative power from neutronic calculations, assembly power reductions due to fissile content depletion and core reshuffling, and fuel material thermal-physical properties. Additionally, a text-based coupling method is developed to facilitate the exchange of information between the neutronic code DRAGON and thermal-hydraulic code VIPRE-01. The new methodology effectively covers the entire assembly burnup lifetime and evaluates the thermal-hydraulic performance against ANS Condition I, II, and III events with respect to the minimum departure from nucleate boiling ratio, peak cladding temperatures, and fuel centerline temperatures. A comprehensive literature survey on the thermal conductivity of posed fuel materials with burnup-dependence has been carried out to model the advanced materials in the thermal-hydraulic code VIPRE-01. Where documented conductivity values are not available, a simplified method for estimating the thermal conductivity has been developed. The new thermal conductivity models are based on established FRAPCON-3 fuel property models used in the nuclear industry, with small adjustments having been made to account for actinide additions. Steady-state and transient thermal-hydraulic analyses are performed with VIPRE- 01 for a reference UO2 assembly design, and two advanced mixed-oxide fuel assembly designs using the new burnup-dependent thermal-hydraulic analysis methodology. All three designs maintain a sufficiently large thermal margin with respect to the minimum departure from nucleate boiling ratio, and maximum cladding and fuel temperatures during partial and complete loss-of-flow accident scenarios. The presence of a thin (Am,Zr)O2 outer layer on the fuel pellet in the two advanced mixed-oxide fuel assembly designs increases maximum fuel temperatures during transient conditions, but does not otherwise greatly compromise the thermal margin of the new designs.

Complex flow structures consisting of branching, multi-scale, hierarchically arranged flow paths can be a beneficial in certain applications by providing lower hydraulic and thermal resistances than conventional flow arrangements. In this study, an experimental approach was used to investigate the hydrodynamic and thermal effects of the complexity, or degree of branching, in microscale complex flow structures. The primary focus of this work was to develop new concepts to advance the current capabilities of complex flow structures through management of complexity. The effects of complexity were determined from experiments performed on a set of microfluidic test sections which were identical except for the complexity of the underlying microchannel configuration. Comparison of the relative hydrodynamic and thermal performance indicates that complexity has a strong effect on both the pressure drop and heat transfer. When the pumping power is taken into account, the results suggest that higher complexity arrangements improve the overall thermal-hydraulic performance. This conclusion was confirmed by the trends observed in the coefficient of performance, a measure of the device thermal efficiency. To address the limitations of conventional fixed-complexity designs, the concept of a variable-complexity flow structure is developed. With a variable-complexity design, the configuration of a branched flow structure can be dynamically controlled to improve performance as operational conditions vary. This concept was successfully demonstrated by developing and testing an active variable-complexity microfluidic device in which pneumatically controlled microvalves were used to create different flow channel configurations. The variable-complexity concept was further refined by developing a microfluidic device with a passive variable-complexity design in which the flow channel configuration changed autonomously based on local temperatures. By using microvalves containing a temperature sensitive polymer, the flow configuration of the device was made thermally adaptive. Experiments were performed to characterize the behavior of the polymer microvalves and the overall device performance. The results showed that the device was capable of tracking changes in external heat sources by adapting and reconfiguring its internal flow structure. The experiments also showed how this variable-complexity design can reduce the pumping power expenditure by automatically directing flow only to areas where it is required. / text

Regenerative heat exchangers are established as a means of heat recovery in many industrial applications. The fixed-bed regenerators are mostly used to transfer heat from hot flue gas to cold air. In this work, several mathematical models of regenerators and several calculation methods were compared, while the preferred method is Willmott’s open method from 1964. Analysis of the influence of geometrical and operational parameters was carried out for the linear regenerator model.

In accident analysis regarding nuclear power plants, it is very common to use thermal hydraulic system codes, such as TRACE, developed by U.S. NRC. In the case of licensing a power plant, this is one of the necessities. TRACE is a relatively new thermal hydraulic system code and a lot of knowledge is needed to implement it in a correct way, especially in accident analysis, where it is a requirement that the rules and statements in Appendix-K, dealing with criteria for ECCS-models, are modelled. In this thesis an improved model of a Swedish Boiling Water Reactor within TRACE is realized and tested. Afterwards, once a working and representative model has been obtained, a sensitivity study in conducted in order to investigate the sensitivity of TRACE for a couple of thermal hydraulic parameters. The sensitivity study is focussing on the eect of the peak cladding temperature, as well as the coolability of the nuclear fuel in terms of quenching and quench-front velocities. It is found to be hard to say unilaterally what the eect of changing a certain number of parameters on the reactor behaviour is. As it turns out to be, although strongly related, the peak cladding temperatures and the quench phenomena can behave dierently

This PhD study, developed at Universitat Politècnica de València (UPV), aims to cover the first phase of the benchmark released by the expert group on Uncertainty Analysis in Modeling (UAM-LWR). The main contribution to the benchmark, made by the thesis' author, is the development of a MATLAB program requested by the benchmark organizers. This is used to generate neutronic libraries to distribute among the benchmark participants. The UAM benchmark pretends to determine the uncertainty introduced by coupled multi-physics and multi-scale LWR analysis codes. The benchmark is subdivided into three phases: 1. Neutronic phase: obtain collapsed and homogenized problem-dependent cross sections and criticality analyses. 2. Core phase: standalone thermohydraulic and neutronic codes. 3. System phase: coupled thermohydraulic and neutronic code. In this thesis the objectives of the first phase are covered. Specifically, a methodology is developed to propagate the uncertainty of cross sections and other neutronic parameters through a lattice physics code and core simulator. An Uncertainty and Sensitivity (U&S) analysis is performed over the cross sections contained in the ENDF/B-VII nuclear library. Their uncertainty is propagated through the lattice physics code SCALE6.2.1, including the collapse and homogenization phase, up to the generation of problem-dependent neutronic libraries. Afterward, the uncertainty contained in these libraries can be further propagated through a core simulator, in this study PARCSv3.2. The module SAMPLER -available in the latest release of SCALE- and DAKOTA 6.3 statistical tool are used for the U&S analysis. As a part of this process, a methodology to obtain neutronic libraries in NEMTAB format -to be used in a core simulator- is also developed. A code-to-code comparison with CASMO-4 is used as a verification. The whole methodology is tested using a Boiling Water Reactor (BWR) reactor type. Nevertheless, there is not any concern or limitation regarding its use in any other type of nuclear reactor. The Gesellschaft für Anlagen und Reaktorsicherheit (GRS) stochastic methodology for uncertainty quantification is used. This methodology makes use of the high-fidelity model and nonparametric sampling to propagate the uncertainty. As a result, the number of samples (determined using the revised Wilks' formula) does not depend on the number of input parameters but only on the desired confidence and uncertainty of output parameters. Moreover, the output Probability Distribution Functions (PDFs) are not subject to normality. The main disadvantage is that each input parameter must have a pre-defined PDF. If possible, input PDFs are defined using information found in the related literature. Otherwise, the uncertainty definition is based on expert judgment. A second scenario is used to propagate the uncertainty of different thermohydraulic parameters through the coupled code TRACE5.0p3/PARCSv3.0. In this case, a PWR reactor type is used and a transient control rod drop occurrence is simulated. As a new feature, the core is modeled chan-by-chan following a fully 3D discretization. No other study is found using a detailed 3D core. This U&S analysis also makes use of the GRS methodology and DAKOTA 6.3. / Este trabajo de doctorado, desarrollado en la Universitat Politècnica de València (UPV), tiene como objetivo cubrir la primera fase del benchmark presentado por el grupo de expertos Uncertainty Analysis in Modeling (UAM-LWR). La principal contribución al benchmark, por parte del autor de esta tesis, es el desarrollo de un programa de MATLAB solicitado por los organizadores del benchmark, el cual se usa para generar librerías neutrónicas a distribuir entre los participantes del benchmark. El benchmark del UAM pretende determinar la incertidumbre introducida por los códigos multifísicos y multiescala acoplados de análisis de reactores de agua ligera. El citado benchmark se divide en tres fases: 1. Fase neutrónica: obtener los parámetros neutrónicos y secciones eficaces del problema específico colapsados y homogenizados, además del análisis de criticidad. 2. Fase de núcleo: análisis termo-hidráulico y neutrónico por separado. 3. Fase de sistema: análisis termo-hidráulico y neutrónico acoplados. En esta tesis se completan los principales objetivos de la primera fase. Concretamente, se desarrolla una metodología para propagar la incertidumbre de secciones eficaces y otros parámetros neutrónicos a través de un código lattice y un simulador de núcleo. Se lleva a cabo un análisis de incertidumbre y sensibilidad para las secciones eficaces contenidas en la librería neutrónica ENDF/B-VII. Su incertidumbre se propaga a través del código lattice SCALE6.2.1, incluyendo las fases de colapsación y homogenización, hasta llegar a la generación de una librería neutrónica específica del problema. Luego, la incertidumbre contenida en dicha librería puede continuar propagándose a través de un simulador de núcleo, para este estudio PARCSv3.2. Para el análisis de incertidumbre y sensibilidad se ha usado el módulo SAMPLER -disponible en la última versión de SCALE- y la herramienta estadística DAKOTA 6.3. Como parte de este proceso, también se ha desarrollado una metodología para obtener librerías neutrónicas en formato NEMTAB para ser usadas en simuladores de núcleo. Se ha realizado una comparación con el código CASMO-4 para obtener una verificación de la metodología completa. Esta se ha probado usando un reactor de agua en ebullición del tipo BWR. Sin embargo, no hay ninguna preocupación o limitación respecto a su uso con otro tipo de reactor nuclear. Para la cuantificación de la incertidumbre se usa la metodología estocástica Gesellschaft für Anlagen und Reaktorsicherheit (GRS). Esta metodología hace uso del modelo de alta fidelidad y un muestreo no paramétrico para propagar la incertidumbre. Como resultado, el número de muestras (determinado con la fórmula revisada de Wilks) no depende del número de parámetros de entrada, sólo depende del nivel de confianza e incertidumbre deseados de los parámetros de salida. Además, las funciones de distribución de probabilidad no están limitadas a normalidad. El principal inconveniente es que se ha de disponer de las distribuciones de probabilidad de cada parámetro de entrada. Si es posible, las distribuciones de probabilidad de entrada se definen usando información encontrada en la literatura relacionada. En caso contrario, la incertidumbre se define en base a la opinión de un experto. Se usa un segundo escenario para propagar la incertidumbre de diferentes parámetros termo-hidráulicos a través del código acoplado TRACE5.0p3/PARCSv3.0. En este caso, se utiliza un reactor tipo PWR para simular un transitorio de una caída de barra. Como nueva característica, el núcleo se modela elemento a elemento siguiendo una discretización totalmente en 3D. No se ha encontrado ningún otro estudio que use un núcleo tan detallado en 3D. También se usa la metodología GRS y el DAKOTA 6.3 para este análisis de incertidumbre y sensibilidad. / Aquest treball de doctorat, desenvolupat a la Universitat Politècnica de València (UPV), té com a objectiu cobrir la primera fase del benchmark presentat pel grup d'experts Uncertainty Analysis in Modeling (UAM-LWR). La principal contribució al benchmark, per part de l'autor d'aquesta tesi, es el desenvolupament d'un programa de MATLAB sol¿licitat pels organitzadors del benchmark, el qual s'utilitza per a generar llibreries neutròniques a distribuir entre els participants del benchmark. El benchmark del UAM pretén determinar la incertesa introduïda pels codis multifísics i multiescala acoblats d'anàlisi de reactors d'aigua lleugera. El citat benchmark es divideix en tres fases: 1. Fase neutrònica: obtenir els paràmetres neutrònics i seccions eficaces del problema específic, col¿lapsats i homogeneïtzats, a més de la anàlisi de criticitat. 2. Fase de nucli: anàlisi termo-hidràulica i neutrònica per separat. 3. Fase de sistema: anàlisi termo-hidràulica i neutrònica acoblats. En aquesta tesi es completen els principals objectius de la primera fase. Concretament, es desenvolupa una metodologia per propagar la incertesa de les seccions eficaces i altres paràmetres neutrònics a través d'un codi lattice i un simulador de nucli. Es porta a terme una anàlisi d'incertesa i sensibilitat per a les seccions eficaces contingudes en la llibreria neutrònica ENDF/B-VII. La seua incertesa es propaga a través del codi lattice SCALE6.2.1, incloent les fases per col¿lapsar i homogeneïtzar, fins aplegar a la generació d'una llibreria neutrònica específica del problema. Després, la incertesa continguda en la esmentada llibreria pot continuar propagant-se a través d'un simulador de nucli, per a aquest estudi PARCSv3.2. Per a l'anàlisi d'incertesa i sensibilitat s'ha utilitzat el mòdul SAMPLER -disponible a l'última versió de SCALE- i la ferramenta estadística DAKOTA 6.3. Com a part d'aquest procés, també es desenvolupa una metodologia per a obtenir llibreries neutròniques en format NEMTAB per ser utilitzades en simuladors de nucli. S'ha realitzat una comparació amb el codi CASMO-4 per obtenir una verificació de la metodologia completa. Aquesta s'ha provat utilitzant un reactor d'aigua en ebullició del tipus BWR. Tanmateix, no hi ha cap preocupació o limitació respecte del seu ús amb un altre tipus de reactor nuclear. Per a la quantificació de la incertesa s'utilitza la metodologia estocàstica Gesellschaft für Anlagen und Reaktorsicherheit (GRS). Aquesta metodologia fa ús del model d'alta fidelitat i un mostreig no paramètric per propagar la incertesa. Com a resultat, el nombre de mostres (determinat amb la fórmula revisada de Wilks) no depèn del nombre de paràmetres d'entrada, sols depèn del nivell de confiança i incertesa desitjats dels paràmetres d'eixida. A més, las funcions de distribució de probabilitat no estan limitades a la normalitat. El principal inconvenient és que s'ha de disposar de les distribucions de probabilitat de cada paràmetre d'entrada. Si és possible, les distribucions de probabilitat d'entrada es defineixen utilitzant informació trobada a la literatura relacionada. En cas contrari, la incertesa es defineix en base a l'opinió d'un expert. S'utilitza un segon escenari per propagar la incertesa de diferents paràmetres termo-hidràulics a través del codi acoblat TRACE5.0p3/PARCSv3.0. En aquest cas, s'utilitza un reactor tipus PWR per simular un transitori d'una caiguda de barra. Com a nova característica, cal assenyalar que el nucli es modela element a element seguint una discretizació totalment 3D. No s'ha trobat cap altre estudi que utilitze un nucli tan detallat en 3D. També s'utilitza la metodologia GRS i el DAKOTA 6.3 per a aquesta anàlisi d'incertesa i sensibilitat.¿ / Mesado Melia, C. (2017). Uncertainty Quantification and Sensitivity Analysis for Cross Sections and Thermohydraulic Parameters in Lattice and Core Physics Codes. Methodology for Cross Section Library Generation and Application to PWR and BWR [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/86167 / TESIS