• Refine Query
  • Source
  • Publication year
  • to
  • Language
  • 47
  • 2
  • Tagged with
  • 49
  • 49
  • 32
  • 18
  • 10
  • 9
  • 9
  • 7
  • 7
  • 6
  • 5
  • 5
  • 5
  • 4
  • 4
  • About
  • The Global ETD Search service is a free service for researchers to find electronic theses and dissertations. This service is provided by the Networked Digital Library of Theses and Dissertations.
    Our metadata is collected from universities around the world. If you manage a university/consortium/country archive and want to be added, details can be found on the NDLTD website.
41

"Teoria e modelamento computacional de aquecimento de plasma por ondas de alfvén no tokamak TCABR" / Theory and computer modelling of Alfvén wave heating in TCABR Tokamak

Sanabria, Edgar Rodolfo Rondán 10 August 2006 (has links)
Neste trabalho apresentamos o estudo da possivilidade de melhores regimes para o uso dos experimentos de aquecimento e geracao de corrente e fluxo de plasma no tokamak TCABR. Apresentamos um estudo dos efeitos de rotacao de plasma em baixa frequencia (low-frequency (LF)), penetração de campo eletromagnético, absorção e forças ponderomotoras no “Tokamak Chauffage Alfvén Brésilien" (TCABR) com ênfase na faixa de frequências de 0, 5–10, 0kHz. Os campos de LF são dirigidos pelo limitador magnético ergódico (ergodic magnetic limiter (EML)) no TCABR. Foi feito um estudo analítico das ondas de Alfvén e ressonância usando modelos simples. Um estudo num´erico tembém foi realizado utilizando três códigos, quais sejam, o código cinético toroidal, o código cilíndrico e o código ALTOK. / In this work we present the study of the determination the best regimes and parameters¶for the heating experiments and current generation and plasma flow in the tokamak TCABR. Study of effects of plasma rotation in low frequency (LF), field penetration, absorption and ponderomotive forces in “Tokamak Chauffage Alfvén Brésilien" (TCABR)is investigated with emphasis in the frequency range of 0, 5–10, 0kHz. The fields of LF are driven by the ergodic magnetic limiter (EML) in TCABR. A qualitative analytical study of the Alfvén waves and their resonances is performed using simple models. A numeric study was carried out using through three codes, called the kinetic totoidal code, the cylindrical code and the ALTOK code.
42

Cálculo da fração de vazio em escoamentos bifásicos (gás/líquido) a partir da identificação de bolhas em imagens digitais / Two-phase flow void fraction estimation based on bubble segmentation and dimensioning using neural nets and modified randomized hough transform

Serra, Pedro Luiz Santos 21 June 2017 (has links)
A Agência Internacional de Energia Atômica (IAEA - \"International Atomic Energy Agency\") vem incentivando o desenvolvimento de sistemas passivos de refrigeração em plantas nucleares visando a simplificação e o incremento da confiabilidade em funções essenciais de segurança nos projetos de uma próxima geração de reatores nucleares refrigerados a água. O principal fundamento desses sistemas é o emprego da circulação natural como sistema de segurança aplicável em operações de desligamento do reator para manutenção ou na ocorrência de acidentes. A circulação natural é um fenômeno que surge em virtude do gradiente de temperatura em pontos diferentes do circuito de refrigeração. Em condições extremas de estabilidade têm-se o estabelecimento do escoamento bifásico gás/líquido podendo configurar-se segundo diferentes regimes. A fração de vazio é reconhecida como um dos parâmetros chave na predição da ocorrência de instabilidades do escoamento bifásico. Apresenta-se neste trabalho uma inovadora metodologia para estimativa da fração de vazio a partir de imagens digitais capturadas diretamente de circuitos experimentais que geram o escoamento bifásico. O método é baseado na aquisição de imagens, com controle da profundidade de campo, de uma seção do Circuito de Circulação Natural (CCN) presente no IPEN/CNEN-SP. A imagem é segmentada com base na inferência fuzzy de diferentes parâmetros de segmentação e ajustada ao foco utilizado na sua aquisição. Ela é varrida de um modo inédito e iterativo, utilizando máscaras de diferentes tamanhos integrando um conjunto de redes neurais com a Transformada Randomizada de Hough. Cada diferente tamanho de máscara é escolhido de acordo com os tamanhos das bolhas que são os objetos de interesse. O volume da bolha é estimado baseado em sua projeção plana capturada nas imagens digitais. O cálculo da fração de vazio considera o volume da seção geométrica do escoamento no tubo de vidro cilíndrico e a profundidade de campo utilizada e nos parâmetros geométricos inferidos para cada bolha detectada. Os resultados mostraram que a integração entre o conjunto de redes neurais e a Transformada Randomizada de Hough aumentaram a robustez das estimativas do sistema. / The International Atomic Energy Agency (IAEA) has been encouraging the use of passive cooling systems in new designs of nuclear power plants. Next nuclear reactor generations are intended to possess simpler and robust safety functions. Natural circulation based systems hold an undoubtedly prominent position among these. Natural circulation phenomenon occurrence depends only on the existence of refrigerant liquid temperature gradient in different sections of the plant refrigerator circuit. The study of limit conditions for these systems has led to instability behavior analysis where many different two-phase flow patterns are present. Void fraction is a key parameter in thermal transfer analysis of theses flow instability conditions. This works presents a new method to estimate void fraction from digital images captured at an experimental two-phase flow circuit. The method is primarily based on depth-of-field controlled image acquisition of a section of a closed loop of natural circulating water through cylindrical glass tubes. This loop is called Natural Circulation Facility (NCF) and is located at Nuclear and Engineering Research Institute in Brazil (IPEN/CNEN-SP). Image is segmented based on fuzzy inference of different segmentation parameters and adjusted to image acquisition focus. The image is then scanned in an inedited way using different-sized masks integrating a set of different artificial neural networks with a modified Randomized Hough Transform. Each different mask size is chosen in accordance to bubble sizes which are objects of interest. The bubble volume is estimated based on two-dimensional projection sizing based on digitally acquired images. Void fraction calculation takes into account the volume of the geometrical section of flow inside cylindrical glass tube considering used depth-of-field. It is also based on the summed bubble geometrical parameters inferred for each detected bubble. The results have shown that integration between artificial-neural-net sets and Randomized Hough Transforms increase system estimations robustness.
43

Desenvolvimento do programa ATHELCORE para análise termohidráulica do núcleo completo de reatores do tipo LMFBR

Francisco Antonio Braz Filho 01 June 1991 (has links)
O conhecimento da distribuição de temperatura e velocidade no núcleo de um reator é necessária para a verificação dos limites do projeto estrutural, a análise de segurança, o desempenho do combustível, etc. O presente trabalho consiste no desenvolvimento do programa ATHELCORE, o qual representa uma expansão do modelo físico utilizado no programa ATHEL, para a análise termohidráulica do núcleo completo de reatores do tipo LMFBR, considerando a transferência de calor entre feixes. O tempo de execução do programa e a memória requeridos são reduzidos quando comparados a programas similares, tais como o COBRA-WC e o COMMIX, fornecendo uma descrição detalhada do campo de velocidade e temperatura. Os resultados obtidos do programa, em relação aos dados experimentais, são considerados satisfatórios.
44

Desenvolvimento de um modelo para dimensionamento da capacidade produtiva de fábricas de combustível nuclear para reatores de pesquisa / Development of a model for dimensioning the production capacity of nuclear fuel factories for research reactors

Negro, Miguel Luiz Miotto 06 October 2017 (has links)
A demanda por combustível nuclear para reatores de pesquisa está aumentando em nível mundial, enquanto várias de suas fábricas têm pequeno volume de produção. Este trabalho estabeleceu um modelo conceitual com duas estratégias para o aumento da capacidade produtiva dessas fábricas. Foram abordadas as fábricas que produzem elementos combustíveis tipo placa carregados com LEU U3Si2-Al, tipicamente usados em reatores nucleares de pesquisa. A primeira estratégia baseia-se na literatura da área de administração da produção e é uma prática frequente nas fábricas em geral. A segunda estratégia aproveita a possibilidade de desmembrar setores produtivos, comum em instalações de produção de combustível nuclear. Ambas as estratégias geraram diferentes cenários de produção, os quais devem ser seguros em relação à criticalidade. Foram coletados dados de uma fábrica real de combustível nuclear para reatores de pesquisa. As duas estratégias foram aplicadas a esses dados com a finalidade de testar o modelo proposto, o que configurou um estudo de caso. A aplicação das estratégias aos dados coletados deu-se por meio de simulação de eventos discretos em computador. Foram criados diversos modelos de simulação para abranger todos os cenários gerados, de forma que o teste indicou um aumento da capacidade produtiva de até 207% sem necessidade de aquisição de novos equipamentos. Os resultados comprovam que o modelo atingiu plenamente o objetivo proposto. Como principal conclusão pode-se apontar a eficácia do modelo proposto, fato que foi validado pelos dados da fábrica. / Although many nuclear fuel factories have small production volumes, the demand for nuclear fuel for research reactors is increasing worldwide. This work established a conceptual model with two strategies to increase the production capacity of these factories. We addressed factories that produce plate-type fuel elements loaded with LEU U3Si2-Al, which are typically used in nuclear research reactors. The first strategy is based on production management literature and is a regular practice in general manufacturing plants. The second strategy takes advantage of the fact that productive sectors can be separated in nuclear fuel production facilities. Both strategies have generated different production scenarios that are assumed to be safe in relation to nuclear criticality. We collected data from a real plant that produces nuclear fuel for research reactors and applied the model to that data, aiming to test the proposed model by setting up a case study. Through the use of computer software, we applied the two strategies to this data by means of discrete events simulation and created several simulation models in order to cover all generated scenarios. Our tests indicated an increase of up to 207% in productive capacity without the need of acquiring new equipment, thus showing that the model has fully achieved its proposed objective. One of the main conclusions that we point out is the models effectiveness, which was validated by the factory data.
45

Desenvolvimento de um modelo para dimensionamento da capacidade produtiva de fábricas de combustível nuclear para reatores de pesquisa / Development of a model for dimensioning the production capacity of nuclear fuel factories for research reactors

Miguel Luiz Miotto Negro 06 October 2017 (has links)
A demanda por combustível nuclear para reatores de pesquisa está aumentando em nível mundial, enquanto várias de suas fábricas têm pequeno volume de produção. Este trabalho estabeleceu um modelo conceitual com duas estratégias para o aumento da capacidade produtiva dessas fábricas. Foram abordadas as fábricas que produzem elementos combustíveis tipo placa carregados com LEU U3Si2-Al, tipicamente usados em reatores nucleares de pesquisa. A primeira estratégia baseia-se na literatura da área de administração da produção e é uma prática frequente nas fábricas em geral. A segunda estratégia aproveita a possibilidade de desmembrar setores produtivos, comum em instalações de produção de combustível nuclear. Ambas as estratégias geraram diferentes cenários de produção, os quais devem ser seguros em relação à criticalidade. Foram coletados dados de uma fábrica real de combustível nuclear para reatores de pesquisa. As duas estratégias foram aplicadas a esses dados com a finalidade de testar o modelo proposto, o que configurou um estudo de caso. A aplicação das estratégias aos dados coletados deu-se por meio de simulação de eventos discretos em computador. Foram criados diversos modelos de simulação para abranger todos os cenários gerados, de forma que o teste indicou um aumento da capacidade produtiva de até 207% sem necessidade de aquisição de novos equipamentos. Os resultados comprovam que o modelo atingiu plenamente o objetivo proposto. Como principal conclusão pode-se apontar a eficácia do modelo proposto, fato que foi validado pelos dados da fábrica. / Although many nuclear fuel factories have small production volumes, the demand for nuclear fuel for research reactors is increasing worldwide. This work established a conceptual model with two strategies to increase the production capacity of these factories. We addressed factories that produce plate-type fuel elements loaded with LEU U3Si2-Al, which are typically used in nuclear research reactors. The first strategy is based on production management literature and is a regular practice in general manufacturing plants. The second strategy takes advantage of the fact that productive sectors can be separated in nuclear fuel production facilities. Both strategies have generated different production scenarios that are assumed to be safe in relation to nuclear criticality. We collected data from a real plant that produces nuclear fuel for research reactors and applied the model to that data, aiming to test the proposed model by setting up a case study. Through the use of computer software, we applied the two strategies to this data by means of discrete events simulation and created several simulation models in order to cover all generated scenarios. Our tests indicated an increase of up to 207% in productive capacity without the need of acquiring new equipment, thus showing that the model has fully achieved its proposed objective. One of the main conclusions that we point out is the models effectiveness, which was validated by the factory data.
46

"Teoria e modelamento computacional de aquecimento de plasma por ondas de alfvén no tokamak TCABR" / Theory and computer modelling of Alfvén wave heating in TCABR Tokamak

Edgar Rodolfo Rondán Sanabria 10 August 2006 (has links)
Neste trabalho apresentamos o estudo da possivilidade de melhores regimes para o uso dos experimentos de aquecimento e geracao de corrente e fluxo de plasma no tokamak TCABR. Apresentamos um estudo dos efeitos de rotacao de plasma em baixa frequencia (low-frequency (LF)), penetração de campo eletromagnético, absorção e forças ponderomotoras no “Tokamak Chauffage Alfvén Brésilien” (TCABR) com ênfase na faixa de frequências de 0, 5–10, 0kHz. Os campos de LF são dirigidos pelo limitador magnético ergódico (ergodic magnetic limiter (EML)) no TCABR. Foi feito um estudo analítico das ondas de Alfvén e ressonância usando modelos simples. Um estudo num´erico tembém foi realizado utilizando três códigos, quais sejam, o código cinético toroidal, o código cilíndrico e o código ALTOK. / In this work we present the study of the determination the best regimes and parameters¶for the heating experiments and current generation and plasma flow in the tokamak TCABR. Study of effects of plasma rotation in low frequency (LF), field penetration, absorption and ponderomotive forces in “Tokamak Chauffage Alfvén Brésilien” (TCABR)is investigated with emphasis in the frequency range of 0, 5–10, 0kHz. The fields of LF are driven by the ergodic magnetic limiter (EML) in TCABR. A qualitative analytical study of the Alfvén waves and their resonances is performed using simple models. A numeric study was carried out using through three codes, called the kinetic totoidal code, the cylindrical code and the ALTOK code.
47

Modelo para análise de circuitos de resfriamento em circulação natural

Silas Joaquim dos Santos 01 June 1988 (has links)
Este trabalho apresenta uma análise de termosifões e seus números adimensionais. O modelo matemático considera pressão constante, escoamento monofásico e incompressível. O modelo simula tanto termosifões abertos quanto fechados, modelando as fontes quentes como núcleos de PWR ou aquecedores elétricos e as fontes frias como trocadores de calor ou reservatórios. O programa computacional STRATS foi desenvolvido baseado neste modelo. Os resultados mostram uma boa concordância tanto com os resultados analíticos quanto com os dados experimentais apresentados na literatura.
48

Implications of advanced computational methods for reactivity initiated accidents in nuclear reactors. / Implicações do uso de métodos computacionais avançados na análise de acidentes iniciados por reatividade em reatores nucleares.

Busquim e Silva, Rodney Aparecido 26 May 2015 (has links)
Advanced computational tools are applied to simulate a nuclear power plant (NPP) control rod assembly ejection (CRE) accident. The impact of these reactivity-initiated accidents (RIAs) on core reactivity behavior, 3D power distribution and stochastic reactivity estimation are evaluated. The three tools used are: the thermal-hydraulic (TH) RELAP5 (R5) code, the neutronic (NK) PARCS (P3D) code, and the coupled version P3D/R5, with specially developed linkage using the environment code MATLAB. This study considers three different-size cores: NPP1 (2772 MWt); NPP2 (530 MWt); and NPP3 (1061 MWt). The three cores have the same general design and control rod assembly (CRA) positions, and the ejected CRA has similar worth and at the same rod ejection pace. The CRE is assessed under both hot zero power (HZP) and hot full power (HFP) conditions. The analyses indicate that RIA modeling and simulation should be carried out through a systematic coding and configuration approaches, otherwise the results will not capture the true transient behavior of the core under analysis. The simulation of one code depends on the appropriate configuration of parameters generated by the other code and on the correct determination of the TH/NK mapping weight factors for the various mesh regions in each of the models. From the design point of view, the standalone codes predict milder magnitude of power and reactivity increase compared to the coupled P3D/R5 simulation. The magnitudes of reduced peak power and reactivity become larger as the core size shrinks. The HFP simulation shows that the three NPPs have the same transient peak value, but the post-transient steady power is lower for a smaller core. The HZP analysis indicates that the transient peak is lower for the smaller core, but the post-transient power occurs at the same level. The three-dimensional (3D) power distributions are different among the HFP and HZP cases, but do not depend on the size of the core. The results indicate: i) HFP: core power increases in the area surrounding the ejected rod/bank assembly, and this increase becomes lower as the NPPs shrinks however, the power is well-distributed after the transient; and ii) HZP: the area surrounding the CRA stays hotter, but the 3D peak assembly factor becomes lower, during and after the transients, as the NPPs shrinks. These features confirm that the smaller cores yield a safer response to a given inserted reactivity compared to larger cores. A stochastic extended Kalman filter (EKF) algorithm is implemented to estimate the reactivity based on the reactor power profile, after the addition of random noise. The inverse point kinetics (IPK) deterministic method is also implemented and the results of the application of EKF and IPK are compared to the P3D/R5 simulation. The following sophisticated strategies made the EKF algorithm robust and accurate: the system is modeled by a set of continuous time nonlinear stochastic differential equations; the code uses a time step directly based on the power measured and applies that to the model for online discretization and linearization; filter tuning goes automatically up from the first time step; and the state noise covariance matrix is updated online at each time step. It was found that the IPK reactivity has higher noise content compared to the EKF reactivity for all cases. Thus, the EKF presents superior and more accurate results. Furthermore, under a small reactivity insertion, the IPK reactivity varies widely from positive to negative values: this variation is not observed within the EKF. A sensitivity analysis for three distinct standard deviation (SD) noise measurements suggests that EKF is superior to IPK method, independent of the noise load magnitude. As the noise content increases, the error between the IPK and P3D/R5 reactivity also increases. A sensitivity analysis for five distinct carry-over effects of different random noise loads indicates that the random addition of different noise loads to the reactor power does not change the overall performance of both algorithms. / Este trabalho aplica métodos computacionais avançados para simular a ejeção de barras de controle (CRE) em uma planta térmica nuclear (NPP). São avaliados o impacto da ocorrência de acidentes iniciados por reatividade (RIAs) na reatividade total, na distribuição da potência em três dimensões (3D) e na determinação da reatividade. As ferramentas utilizadas são: o código termo-hidráulico (TH) RELAP5 (R5), o código neutrônico (NK) PARCS (P3D), a versão acoplada P3D/R5, e o ambiente computacional MATLAB. Este estudo considera três reatores nucleares de diferentes tamanhos: NPP1 (2772 MWT); NPP2 (530 MWt); e NPP3 (1061 MWt). Os três núcleos possuem projeto similar e idêntica posição dos grupos das barras de controle (CRA), além do mesmo valor de reatividade diferencial das CRA ejetadas e idêntica velocidade de ejeção. A ocorrência da CRE é avaliada sob condições de hot zero power (HZP) e de hot full power (HFP). As análises indicam que a modelagem e a simulação de RIAs devem ser realizadas sistematicamente, caso contrário os resultados não irão refletir o comportamento em regime transitório do núcleo. A simulação de um modelo em um código depende da apropriada configuração de parâmetros gerados pelo outro código e da determinação adequada do mapeamento TH/NK para as várias malhas dos modelos. Do ponto de vista de projeto, a utilização de códigos independentes resulta em cálculos de potência e reatividade conservadores em comparação com os resultados utilizando-se P3D/R5. Os picos de potência e de reatividade são menores à medida que o núcleo encolhe. A simulação em condições de HFP resulta em valores de pico de potência similares durante transitório para as três NPPs, mas a potência de pós-transitórios é menor para o menor núcleo. A análise em condições de HZP também indica que o valor máximo durante o transitório é menor para o menor núcleo, mas o pós-transitórios ocorre aos mesmos níveis de potência das demais NPPS. A distribuição de potência em 3D também apresenta resultados distintos para condições de HFP e HZP, mas tais resultados são independentes do tamanho do núcleo: i) HFP: há um aumento da potência do núcleo em torno da CRE, mas tal comportamento diminui para núcleos menores - no entanto, a potência é bem distribuída após o transitório; e ii) HZP: há aumento de potência na área do CRE, mas o pico de potência em 3D é menor durante e depois dos transitórios para núcleos menores. Tais características indicam que os núcleos menores respondem de forma mais segura quando da inserção de reatividade em comparação a reatores de maiores dimensões. O método estocástico de filtragem de Kalman estendido (EKF) foi codificado para estimar a reatividade com base no perfil de potência da NPP, após a adição de ruído aleatório. O método determinístico da cinética pontual inversa (IPK) também foi implementado e os resultados da aplicação dos algoritmos do EKF e IPK foram comparados com os resultados da simulação do P3D/R5. As seguintes estratégias, implementadas neste trabalho, possibilitaram a aplicação robusta e precisa do EKF: o sistema foi modelado por um conjunto de equações diferenciais não-lineares estocásticas de tempo contínuo; o algoritmo obtém o passo de tempo diretamente da potência medida e aplica-o ao modelo para a discretização e linearização online; o ajuste do filtro ocorre automaticamente a partir do primeiro passo de tempo; e a matriz de covariância do ruído no estado é atualizada online. Verificou-se que a reatividade calculada pelo método IPK possui maior nível de ruído quando comparada ao EKF para todos os casos estudados. Portanto, o EKF apresenta resultados superiores e mais precisos. Além disso, sob uma pequena inserção de reatividade, a reatividade calculada pelo método IPK varia consideravelmente de valores positivos para negativos: esta variação não é observada com o EKF. Uma análise de sensibilidade para três desvios padrão (SD) sugere que o algoritmo EKF é superior ao método IPK, independente da magnitude do ruído. Com o aumento da magnitude do ruído, o erro entre as reatividades calculadas pelo IPK e pelo P3D/R5 aumenta. A análise de sensibilidade para cinco ruídos aleatórios indica que a adição de ruído na potência do reator não altera o desempenho global de ambos os algoritmos.
49

Implications of advanced computational methods for reactivity initiated accidents in nuclear reactors. / Implicações do uso de métodos computacionais avançados na análise de acidentes iniciados por reatividade em reatores nucleares.

Rodney Aparecido Busquim e Silva 26 May 2015 (has links)
Advanced computational tools are applied to simulate a nuclear power plant (NPP) control rod assembly ejection (CRE) accident. The impact of these reactivity-initiated accidents (RIAs) on core reactivity behavior, 3D power distribution and stochastic reactivity estimation are evaluated. The three tools used are: the thermal-hydraulic (TH) RELAP5 (R5) code, the neutronic (NK) PARCS (P3D) code, and the coupled version P3D/R5, with specially developed linkage using the environment code MATLAB. This study considers three different-size cores: NPP1 (2772 MWt); NPP2 (530 MWt); and NPP3 (1061 MWt). The three cores have the same general design and control rod assembly (CRA) positions, and the ejected CRA has similar worth and at the same rod ejection pace. The CRE is assessed under both hot zero power (HZP) and hot full power (HFP) conditions. The analyses indicate that RIA modeling and simulation should be carried out through a systematic coding and configuration approaches, otherwise the results will not capture the true transient behavior of the core under analysis. The simulation of one code depends on the appropriate configuration of parameters generated by the other code and on the correct determination of the TH/NK mapping weight factors for the various mesh regions in each of the models. From the design point of view, the standalone codes predict milder magnitude of power and reactivity increase compared to the coupled P3D/R5 simulation. The magnitudes of reduced peak power and reactivity become larger as the core size shrinks. The HFP simulation shows that the three NPPs have the same transient peak value, but the post-transient steady power is lower for a smaller core. The HZP analysis indicates that the transient peak is lower for the smaller core, but the post-transient power occurs at the same level. The three-dimensional (3D) power distributions are different among the HFP and HZP cases, but do not depend on the size of the core. The results indicate: i) HFP: core power increases in the area surrounding the ejected rod/bank assembly, and this increase becomes lower as the NPPs shrinks however, the power is well-distributed after the transient; and ii) HZP: the area surrounding the CRA stays hotter, but the 3D peak assembly factor becomes lower, during and after the transients, as the NPPs shrinks. These features confirm that the smaller cores yield a safer response to a given inserted reactivity compared to larger cores. A stochastic extended Kalman filter (EKF) algorithm is implemented to estimate the reactivity based on the reactor power profile, after the addition of random noise. The inverse point kinetics (IPK) deterministic method is also implemented and the results of the application of EKF and IPK are compared to the P3D/R5 simulation. The following sophisticated strategies made the EKF algorithm robust and accurate: the system is modeled by a set of continuous time nonlinear stochastic differential equations; the code uses a time step directly based on the power measured and applies that to the model for online discretization and linearization; filter tuning goes automatically up from the first time step; and the state noise covariance matrix is updated online at each time step. It was found that the IPK reactivity has higher noise content compared to the EKF reactivity for all cases. Thus, the EKF presents superior and more accurate results. Furthermore, under a small reactivity insertion, the IPK reactivity varies widely from positive to negative values: this variation is not observed within the EKF. A sensitivity analysis for three distinct standard deviation (SD) noise measurements suggests that EKF is superior to IPK method, independent of the noise load magnitude. As the noise content increases, the error between the IPK and P3D/R5 reactivity also increases. A sensitivity analysis for five distinct carry-over effects of different random noise loads indicates that the random addition of different noise loads to the reactor power does not change the overall performance of both algorithms. / Este trabalho aplica métodos computacionais avançados para simular a ejeção de barras de controle (CRE) em uma planta térmica nuclear (NPP). São avaliados o impacto da ocorrência de acidentes iniciados por reatividade (RIAs) na reatividade total, na distribuição da potência em três dimensões (3D) e na determinação da reatividade. As ferramentas utilizadas são: o código termo-hidráulico (TH) RELAP5 (R5), o código neutrônico (NK) PARCS (P3D), a versão acoplada P3D/R5, e o ambiente computacional MATLAB. Este estudo considera três reatores nucleares de diferentes tamanhos: NPP1 (2772 MWT); NPP2 (530 MWt); e NPP3 (1061 MWt). Os três núcleos possuem projeto similar e idêntica posição dos grupos das barras de controle (CRA), além do mesmo valor de reatividade diferencial das CRA ejetadas e idêntica velocidade de ejeção. A ocorrência da CRE é avaliada sob condições de hot zero power (HZP) e de hot full power (HFP). As análises indicam que a modelagem e a simulação de RIAs devem ser realizadas sistematicamente, caso contrário os resultados não irão refletir o comportamento em regime transitório do núcleo. A simulação de um modelo em um código depende da apropriada configuração de parâmetros gerados pelo outro código e da determinação adequada do mapeamento TH/NK para as várias malhas dos modelos. Do ponto de vista de projeto, a utilização de códigos independentes resulta em cálculos de potência e reatividade conservadores em comparação com os resultados utilizando-se P3D/R5. Os picos de potência e de reatividade são menores à medida que o núcleo encolhe. A simulação em condições de HFP resulta em valores de pico de potência similares durante transitório para as três NPPs, mas a potência de pós-transitórios é menor para o menor núcleo. A análise em condições de HZP também indica que o valor máximo durante o transitório é menor para o menor núcleo, mas o pós-transitórios ocorre aos mesmos níveis de potência das demais NPPS. A distribuição de potência em 3D também apresenta resultados distintos para condições de HFP e HZP, mas tais resultados são independentes do tamanho do núcleo: i) HFP: há um aumento da potência do núcleo em torno da CRE, mas tal comportamento diminui para núcleos menores - no entanto, a potência é bem distribuída após o transitório; e ii) HZP: há aumento de potência na área do CRE, mas o pico de potência em 3D é menor durante e depois dos transitórios para núcleos menores. Tais características indicam que os núcleos menores respondem de forma mais segura quando da inserção de reatividade em comparação a reatores de maiores dimensões. O método estocástico de filtragem de Kalman estendido (EKF) foi codificado para estimar a reatividade com base no perfil de potência da NPP, após a adição de ruído aleatório. O método determinístico da cinética pontual inversa (IPK) também foi implementado e os resultados da aplicação dos algoritmos do EKF e IPK foram comparados com os resultados da simulação do P3D/R5. As seguintes estratégias, implementadas neste trabalho, possibilitaram a aplicação robusta e precisa do EKF: o sistema foi modelado por um conjunto de equações diferenciais não-lineares estocásticas de tempo contínuo; o algoritmo obtém o passo de tempo diretamente da potência medida e aplica-o ao modelo para a discretização e linearização online; o ajuste do filtro ocorre automaticamente a partir do primeiro passo de tempo; e a matriz de covariância do ruído no estado é atualizada online. Verificou-se que a reatividade calculada pelo método IPK possui maior nível de ruído quando comparada ao EKF para todos os casos estudados. Portanto, o EKF apresenta resultados superiores e mais precisos. Além disso, sob uma pequena inserção de reatividade, a reatividade calculada pelo método IPK varia consideravelmente de valores positivos para negativos: esta variação não é observada com o EKF. Uma análise de sensibilidade para três desvios padrão (SD) sugere que o algoritmo EKF é superior ao método IPK, independente da magnitude do ruído. Com o aumento da magnitude do ruído, o erro entre as reatividades calculadas pelo IPK e pelo P3D/R5 aumenta. A análise de sensibilidade para cinco ruídos aleatórios indica que a adição de ruído na potência do reator não altera o desempenho global de ambos os algoritmos.

Page generated in 0.0817 seconds