Submitted by Claudinei Pracidelli (cpracide@ipen.br) on 2015-03-17T10:41:16Z No. of bitstreams: 0 / Made available in DSpace on 2015-03-17T10:41:16Z (GMT). No. of bitstreams: 0 / Tese (Doutorado em Tecnologia Nuclear) / IPEN/T / Instituto de Pesquisas Energeticas e Nucleares - IPEN-CNEN/SP

Made available in DSpace on 2014-10-09T12:35:04Z (GMT). No. of bitstreams: 0 / Made available in DSpace on 2014-10-09T14:00:44Z (GMT). No. of bitstreams: 0 / Este trabalho apresenta a realização de experimentos no reator nuclear IPEN/MB-01 submetido à presença de um refletor com água pesada instalado na sua face oeste. Após a instalação do refletor no reator foram conduzidos três tipos de experimentos: A calibração das barras de controle, a verificação da influência do aumento da temperatura do moderador na reatividade e a medição das taxas de reações por meio da irradiação de fios e folhas de ativação. Devido às propriedades nucleares de interação de nêutrons com água pesada, notadamente sua elevada capacidade de espalhamento e sua baixa capacidade de absorção, é possível notar alterações no funcionamento do reator observadas pelas mudanças dos padrões de retiradas e inserções de barras de controle no núcleo. Essas alterações são apresentadas no decorrer do trabalho. / Dissertação (Mestrado) / IPEN/D / Instituto de Pesquisas Energeticas e Nucleares - IPEN-CNEN/SP

Submitted by Claudinei Pracidelli (cpracide@ipen.br) on 2015-03-17T10:41:16Z No. of bitstreams: 0 / Made available in DSpace on 2015-03-17T10:41:16Z (GMT). No. of bitstreams: 0 / Este trabalho apresenta um novo método para o cálculo dos coecientes de sensibilidade, através da união do metodo diferencial e da teoria da perturbação generalizada, que são os dois métodos tradicionalmente utilizados em física de reatores para a obtenção de tais grandezas. Esses dois métodos apresentam algumas deciências tornando os cálculos dos coeficientes de sensibilidade lentos ou computacionalmente exaustivos, mas unindo-os e possível eliminar as deciências apresentadas por ambos e obter uma nova equação para o coe- ciente de sensibilidade. O método proposto neste trabalho foi aplicado em um reator do tipo PWR , onde foi feita análise de sensibilidade da produção e da razão de conversão do 239Pu, para um ciclo de 120 dias de queima. O código utilizado para a análise de queima e análise de sensibilidade, o CINEW, foi desenvolvido durante este trabalho e os resultados obtidos foram comparados com os códigos amplamente utilizados em física de reatores, como o CINDER e o SERPENT. As conclusões obtidas foram que o novo método matemático para a obtenção dos coeficientes de sensibilidade e o CINEW, além de fornecer agilidade numérica também presentam eciência e segurança. Pois o novo método matemático para a obtenção dos coeficientes quando comparados com os métodos tradicionais utilizados para a análise de sensibilidade, mostram resultados satisfatórios, mesmo quando o método utiliza aproximações matemáticas que diferem do método proposto, e com a vantagem de não apresentar as deciências apresentadas pelos métodos diferencial e da teoria da perturbação generalizada. As análises de queima obtidas pelo CINEW foram comparadas com o CINDER, que mostraram uma diferença aceitável, apesar do CINDER apresentar alguns problemas computacionais que advém da época em que foi feito. A originalidade deste trabalho e a aplicação do método proposto em problemas que envolvem dependência temporal e a elaboração do primerio código nacional que faz análise de queima e análise de sensibilidade. / Tese (Doutorado em Tecnologia Nuclear) / IPEN/T / Instituto de Pesquisas Energeticas e Nucleares - IPEN-CNEN/SP

Standardmäßig erfolgt die Datenaufbereitung der Neutronenwirkungsquerschnitte für Reaktorkernrechnungen mit 2D-Zellcodes. Ziel dieser Arbeit war es, einen 3D-Zellcode zu entwickeln, mit diesem Code 3D-Effekte zu untersuchen und die Notwendigkeit einer 3D-Datenaufbereitung der Neutronenwirkungsquerschnitte zu bewerten. Zur Berechnung des Neutronentransports wurde die Methode der Erststoßwahrscheinlichkeiten, die mit der Ray-Tracing-Methode berechnet werden, gewählt. Die mathematischen Algorithmen wurden in den 2D/3D-Zellcode TransRay umgesetzt. Für den Geometrieteil des Programms wurde das Geometriemodul eines Monte-Carlo-Codes genutzt. Das Ray-Tracing in 3D wurde auf Grund der hohen Rechenzeiten parallelisiert. Das Programm TransRay wurde an 2D-Testaufgaben verifiziert. Für einen Druckwasser-Referenzreaktor wurden folgende 3D-Probleme untersucht: Ein teilweise eingetauchter Regelstab und Void (Vakuum oder Dampf) um einen Brennstab als Modell einer Dampfblase. Alle Probleme wurden zum Vergleich auch mit den Programmen HELIOS (2D) und MCNP (3D) nachgerechnet. Die Abhängigkeit des Multiplikationsfaktors und der gemittelten Zweigruppenquerschnitte von der Eintauchtiefe des Regelstabes bzw. von der Höhe der Dampfblase wurden untersucht. Die 3D berechneten Zweigruppenquerschnitte wurden mit drei üblichen Näherungen verglichen: Lineare Interpolation, Interpolation mit Flusswichtung und Homogenisierung. Am 3D-Problem des Regelstabes zeigte sich, dass die Interpolation mit Flusswichtung eine gute Näherung ist. Demnach ist hier eine 3D-Datenaufbereitung nicht notwendig. Beim Testfall des einzelnen Brennstabs, der von Void umgeben ist, erwiesen sich die drei Näherungen für die Zweigruppenquerschnitte als unzureichend. Demnach ist eine 3D-Datenaufbereitung notwendig. Die einzelne Brennstabzelle mit Void kann als der Grenzfall eines Reaktors angesehen werden, in dem sich eine Phasengrenzfläche herausgebildet hat.

Standardmäßig erfolgt die Datenaufbereitung der Neutronenwirkungsquerschnitte für Reaktorkernrechnungen mit 2D-Zellcodes. Ziel dieser Arbeit war es, einen 3D-Zellcode zu entwickeln, mit diesem Code 3D-Effekte zu untersuchen und die Notwendigkeit einer 3D-Datenaufbereitung der Neutronenwirkungsquerschnitte zu bewerten. Zur Berechnung des Neutronentransports wurde die Methode der Erststoßwahrscheinlichkeiten, die mit der Ray-Tracing-Methode berechnet werden, gewählt. Die mathematischen Algorithmen wurden in den 2D/3D-Zellcode TransRay umgesetzt. Für den Geometrieteil des Programms wurde das Geometriemodul eines Monte-Carlo-Codes genutzt. Das Ray-Tracing wurde auf Grund der hohen Rechenzeiten parallelisiert. Das Programm TransRay wurde an 2D-Testaufgaben verifiziert. Für einen Druckwasser-Referenzreaktor wurden folgende 3D-Probleme untersucht: Ein teilweise eingetauchter Regelstab und Void (bzw. Moderator mit geringerer Dichte) um einen Brennstab als Modell einer Dampfblase. Alle Probleme wurden zum Vergleich auch mit den Programmen HELIOS (2D) und MCNP (3D) nachgerechnet. Die Abhängigkeit des Multiplikationsfaktors und der gemittelten Zweigruppenquerschnitte von der Eintauchtiefe des Regelstabes bzw. von der Höhe der Dampfblase wurden untersucht. Die 3D berechneten Zweigruppenquerschnitte wurden mit drei üblichen Näherungen verglichen: linearer Interpolation, Interpolation mit Flusswichtung und Homogenisierung. Am 3D-Problem des Regelstabes zeigte sich, dass die Interpolation mit Flusswichtung eine gute Näherung ist. Demnach ist hier eine 3D-Datenaufbereitung nicht notwendig. Beim Testfall des einzelnen Brennstabs, der von Void (bzw. Moderator geringerer Dichte) umgeben ist, erwiesen sich die drei Näherungen für die Zweigruppenquerschnitte als unzureichend. Demnach ist eine 3D-Datenaufbereitung notwendig. Die einzelne Brennstabzelle mit Void kann als der Grenzfall eines Reaktors angesehen werden, in dem sich eine Phasengrenzfläche herausgebildet hat.

McMaster University DOCTOR OF PHILOSOPHY (2022) Hamilton, Ontario (Engineering Physics) TITLE: Development and Evaluation of Polaris CANDU Geometry Modelling and of TRACE_Mac/PARCS_Mac Coupling with RRS for CANDU Analysis AUTHOR: Simon Younan, M.A.Sc. (McMaster University), B.Eng. (McMaster University) SUPERVISOR: Dr. David Novog NUMBER OF PAGES: xiv, 163 / In the field of nuclear safety analysis, as computers have become more powerful, there has been a trend away from low-fidelity models using conservative assumptions, to high-fidelity best-estimate models combined with uncertainty analysis. A number of these tools have been developed in the United States, due to the popularity of light water reactors. These include the SCALE analysis suite developed by ORNL, as well as the PARCS and TRACE tools backed by the USNRC. This work explores adapting the capabilities of these tools to the analysis of CANDU reactors. The Polaris sequence, introduced in SCALE 6.2, was extended in this work to support CANDU geometries and compared to existing SCALE sequences such as TRITON. Emphasis was placed on the Embedded Self-Shielding Method (ESSM), introduced with Polaris. Both Polaris and ESSM were evaluated and found to perform adequately for CANDU geometries. The accuracy of ESSM was found to improve when the precomputed selfshielding factors were updated using a CANDU representation. The PARCS diffusion code and the TRACE system thermalhydraulics code were coupled, using the built-in coupling capability between the two codes. In addition, the Exterior Communications Interface (ECI), used for coupling with TRACE, was utilized. A Python interface to the ECI library was developed in this work and used to couple an RRS model written in Python to the coupled PARCS/TRACE model. A number of code modifications were made to accommodate the required coupling and correct code deficiencies, with the modified versions named PARCS_Mac and TRACE_Mac. The coupled codes were able to simulate multiple transients based on prior studies as well as operational events. The code updates performed in this work may be used for many future studies, particularly for uncertainty propagation through a full set of calculations, from the lattice model to a full coupled system model. / Thesis / Doctor of Philosophy (PhD) / Modern nuclear safety analysis tools offer more accurate predictions for the safety and operation of nuclear reactors, including CANDU reactors. These codes take advantage of modern computer hardware, and also a shift in philosophy from conservative analysis to best estimate plus uncertainty analysis. The goal of this thesis was to adapt a number of modern tools to support CANDU analysis and uncertainty propagation, with a particular emphasis on coupling of multiple interacting models. These tools were then demonstrated, and results analyzed. The simulations performed in this work were successful in producing results comparable to prior studies along with experimental and operational data. This included the simulation of four weeks of reactor operation including “shim mode” operation. Sensitivity and uncertainty analyses were performed over the course of the work to quantify the precision and significance of the results as well as to identify areas of interest for future research.

The goal of this thesis is to study the transient behaviour of the SLOWPOKE-2 reactor using Monte-Carlo simulations with the G4-STORK code. G4-STORK is a 3-dimensional Monte-Carlo code derived from the GEANT4 physics simulation toolkit. Methods were developed for the proper treatment of delayed neutrons and a lumped capacitance model was used to track the time-dependent fuel properties (temperature, density) based on the fission power. By validating the methods in G4-STORK with experimental measurements we hope to extend our understanding of reactor transients as well as further develop our methods to model the transients of the next generation reactor designs. A SLOWPOKE-2 reactor such as the one at RMC was chosen for simulation due to its compact size, and well-known transient response of control rod removal and measured temperature feedback. Static simulations in G4-STORK find a neutron flux of order 10^12 cm−2 s−1 which agrees with experiment and a control rod worth of (4.9 ± 2.0) mk compared to the experimentally measured worth of 5.45 mk. Transient simulations from rod pluck-out find similar trends to the experimental findings as our results suggest a negative temperature feedback due to the doppler broadening of the U-238 absorption spectrum which contributes to the overall safety mechanism seen in the SLOWPOKE reactor. It is determined that the methods in G4-STORK provide a reasonable ability to simulate reactor transients and it is recommended that a full-core thermal-hydraulics model be coupled to G4-STORK to achieve a higher level of accuracy. / Thesis / Master of Applied Science (MASc)

<p>When the material or geometry of a reactor varies with time, the neutron flux will respond in the form of a reactor transient. These transients can occur during normal operations when control rods are moved or the reactor is refuelled (CANDU). During a reactor accident, the transient response is especially important because the reactor properties vary quickly with large amplitudes. Therefore, better understanding these conditions allows for improved identification, prevention and mitigation of reactor transients. However, current nuclear simulation codes are generally limited in their ability to model transient behaviour.</p> <p>The NStable code was created to model time-dependent neutron populations in multiplying mediums using the Geant4 Monte Carlo toolkit. The neutron population is allowed to evolve in time, but is periodically renormalized so that the total number of neutrons is constrained within a manageable range. This ensures that the simulation is viable even in highly sub- or supercritical environments. Since Geant4 was not intrinsically designed to track a neutron population over "long" time periods (up to 10 s), the population renormalization mechanisms needed to be created and integrated with Geant4. Additionally, nuclear reactor analysis functionality was added to calculate important quantities such as k_eff.</p> <p>The NStable code was validated using three established nuclear simulation codes: MCNP 5, DRAGON 3.06J, and TART 2005. The validation cases compared spatial distributions and criticality estimates for either homogeneous spheres (uranium-235 or a uranium-heavy water mixture) or the standard CANDU 6 lattice cell. For all three systems, the criticality estimates in NStable agreed with the appropriate validation code within 10 mk (TART for the spheres and DRAGON for the CANDU 6 lattice). Finally, the NStable code was also used to simulate a temperature transient in a UHW sphere where the temperature linear increased by 700 K over 50 ms. In response to the increasing temperature, k_eff decreased by 100 mk over the same period. In the future, transient modelling in NStable should be investigated further to reproduce actual experimental results, and to couple NStable with a thermohydraulics code to simulate a full transient response.</p> / Master of Applied Science (MASc)

Os experimentos realizados e apresentados nesta dissertação resultaram na caracterização neutrônica de núcleos na configuração cilíndrica com 30 varetas combustíveis de diâmetro com um espaço, criado pela retirada de 16 varetas centrais, preenchido com água leve (H2O) ou água pesada (D2O) no reator nuclear de pesquisa IPEN/MB-01. Nestes núcleos, efetuou-se experimentos de correlação de canais nucleares, calibração de barras de controle e irradiação de detectores de ativação de diversos materiais em forma de folhas, cujas faixas energéticas de atuação abrangem grande parte do espectro de nêutrons do núcleo do reator, para a obtenção de parâmetros nucleares, como excessos de reatividade, reatividades totais, atividades saturadas por núcleo alvo, razões espectrais, razões de cádmio e fluxo de nêutrons multigrupo. Com a irradiação de fios de ativação de ouro na parte radial foram obtidas as formas espaciais dos fluxos de nêutrons térmicos e epitérmicos. Os resultados mostraram as características espectrais dessa nova configuração com o espaço das 16 varetas combustíveis preenchido com os dois materiais moderadores. No espaço com a água leve houve um aumento significativo de 294% do fluxo de nêutrons térmicos em comparação com a configuração padrão retangular de 28×26 varetas combustíveis. Com a água pesada aumentou-se a reatividade do sistema com &rho; = (783 ± 54) pcm a mais de excesso de reatividade que na configuração com água leve. Os resultados calculados foram simulados nos códigos computacionais MCNP5, SANDBP e CITATION, onde se obtiveram resultados acurados e precisos para as atividades saturadas por núcleo alvo, as distribuições energéticas e espaciais dos fluxos de nêutrons da parte ativa e de parte do refletor e as comparações diretas das seções de choque entre as razões espectrais experimentais e calculadas. / The experiments performed and presented in this thesis results in the neutronic characterization of the core with cylindrical configuration with 30 fuel rods diameter and a space, created by the removal of 16 central rods, filled with light water (H2O) or heavy water (D2O) in the IPEN/MB-01 nuclear research reactor. In these cores were performed experiments of nuclear channels correlation, control rod worth and irradiation of activation detectors of various materials in foils shapes, whose energy performances cover much of the reactor core neutron spectrum, to obtain nuclear parameters, such as, reactivity excesses, total reactivities, saturated activities per target nucleus, spectral ratios, cadmium ratios and multigroup neutron flux. Activation gold wires detectors were irradiated in radial part to obtain the spatial forms of thermal and epithermal neutron fluxes. The results show the spectral characteristics of this new configuration with the space of 16 fuel rods filled with the two moderator materials. In the space with light water there was a significant increase of 294% of the thermal neutron flux compared to standard rectangular configuration of 28×26 fuel rods. With heavy water the system reactivity was increased, more &rho; = (783 ± 54) pcm in excess reactivity than in the light water configuration. The calculated results were simulated in computational codes MCNP5, SANDBP and CITATION, where accurate and precise results were obtained for saturated activities per target nucleus, the energy and spatial distributions of the neutron fluxes for the active part and part of the reflector and the direct comparisons of cross sections between the experimental and calculated spectral ratios.

Neste trabalho apresentamos um estudo das reações nucleares de alta energia que são fundamentais na definição do termo fonte dos reatores nucleares subcríticos acionados por fonte externa. Estas reações nucleares, também conhecidas como \"spallation\", consistem na interação de hádrons de alta energia com os núcleons do núcleo atômico. A fenomenologia destas reações consiste em duas etapas, sendo que à primeira, o próton interage através de espalhamentos múltiplos, em um processo denominado cascata intra-nuclear seguido da etapa na qual o núcleo excitado oriundo da cascata intranuclear ou evapora partículas de forma a atingir estados energéticos moderados ou fissiona, em um processo conhecido como competição entre evaporação e fissão. Neste trabalho os principais modelos nucleares, os modelos de Bertini e Cugnon, são revistos, pois estes modelos são fundamentais para propósito de projeto devido à falta de dados nucleares avaliados para estas reações. A implementação e validação dos métodos de cálculo para o projeto destas fontes são realizadas. A implementação da metodologia é realizada utilizando o programa MCNPX ( \"Monte Carlo N-Particle eXtended\"), dedicado para cálculos de transporte destas partículas e a validação é realizada mediante uma cooperação internacional junto a um projeto coordenado de pesquisa da Agencia Internacional de Energia Atômica e trabalhos disponíveis. O objetivo é qualificar os cálculos relacionados às reações nucleares e os canais de desexcitação envolvidos. O CRISP, um código nacional para a descrição da fenomenologia das reações envolvidas, também foi estudado e os modelos implementados no código foram revistos e melhorados de forma a dar continuidade ao seu processo de qualificação. Devido às limitações dos principais modelos na descrição de produção de nuclídeos leves, a reação de multi-fragmentação foi estudada. As discrepâncias nos cálculos de produção destes nuclídeos são atribuídas à falta do canal de multi-fragmentação estatística do núcleo. A implementação deste canal foi realizada para a aplicação em reações de altas energias junto ao código CRISP de forma a reproduzir a produção de nuclídeos leves, bem como sua validação mediante a comparação com dados experimentais disponíveis para este fenômeno, obtendo com isso uma melhor reprodução de todo o espectro de produção de nuclídeos do processo. / This work presents a study of high energy nuclear reactions which are fundamental to dene the source term in accelerator driven systems. These nuclear reactions, also known as spallation, consist in the interaction of high energetic hadrons with nucleons in the atomic nucleus. The phenomenology of these reactions consist in two step. In the rst, the proton interacts through multiple scattering in a process called intra-nuclear cascade. It is followed by a step in which the excited nucleus, coming from the intranuclear cascade, could either, evaporates particles to achieve a moderate energy state or ssion. This process is known as competition between evaporation and ssion. In this work the main nuclear models, Bertini and Cugnon are reviewed, since these models are fundamental for design purposes of the source term in ADS, due to lack of evaluated nuclear data for these reactions. The implementation and validation of the calculation methods for the design v of the source is carried out to implement the methodology of source design using the program MCNPX (Monte Carlo N-Particle eXtended), devoted to calculation of transport of these particles and the validation performed by an international cooperation together with a Coordinated Research Project (CRP) of the International Atomic Energy Agency and available jobs, in order to qualify the calculations on nuclear reactions and the de-excitation channels involved, providing a state of the art of design and methodology for calculating external sources of spallation for source driven systems. The CRISP, is a brazilian code for the phenomenological description of the reactions involved and the models implemented in the code were reviewed and improved to continue the qualication process. Due to failure of the main models in describing the production of light nuclides, the multifragmentation reaction model was studied. Because the discrepancies in the calculations of production of these nuclides are attributes to the lack of reaction channel and the implementation of this channel was carried out for applications in high energy reactions with the CRISP code to reproduce the production of light nuclides, as well, as its validation by comparison with experimental data available for this phenomenon. Thus, obtaining a better reproduction of the whole spectrum of production of nuclides in the process.