Sobre a Equação da Difusão de Nêutrons dependente do tempo:
a) Qual é o significado físico de cada um dos termos envolvidos na equação acima?
b) Sob quais condições físicas ou geométricas este modelo de difusão deixa de ser válido?
c) Assuma as hipóteses de estado estacionário e de um reator crítico sem fonte externa independente (
). Modele a difusão substituindo o operador Laplaciano pelo buckling geométrico (
) e mostre algebricamente que o fator de multiplicação pode ser escrito como:
d) Com base na expressão deduzida no item (c), analise qualitativamente o que ocorre com o valor de
nas situações abaixo, justificando com base nos termos da equação:
- O tamanho físico do reator é reduzido (sem alterar a composição material).
- Barras de controle (material absorvedor forte) são inseridas no núcleo.
- O enriquecimento do combustível é aumentado.
b) Sob quais condições físicas ou geométricas este modelo de difusão deixa de ser válido?
c) Assuma as hipóteses de estado estacionário e de um reator crítico sem fonte externa independente (
Passo 1
O enunciado dá a equação da difusão de nêutrons dependente do tempo:
e pede:
a) Explicar fisicamente cada termo da equação.
b) Dizer em que condições físicas/geométricas o modelo de difusão deixa de ser válido.
c) Assumindo estado estacionário e que a fonte é a fissão , substituir
pelo buckling geométrico
e mostrar algebricamente que o fator de multiplicação pode ser escrito como
a) Explicar fisicamente cada termo da equação.
b) Dizer em que condições físicas/geométricas o modelo de difusão deixa de ser válido.
c) Assumindo estado estacionário e que a fonte é a fissão
As hipóteses dadas a usar são:
(estado estacionário) e
(substituição pelo buckling geométrico).
d) A partir da expressão de (c), analisar qualitativamente o que ocorre com quando: (i) o tamanho do reator é reduzido (mesma composição), (ii) inserem-se barras de controle absorvedoras, (iii) aumenta-se o enriquecimento do combustível.
d) A partir da expressão de (c), analisar qualitativamente o que ocorre com
Passo 2
Vamos entender o que é cada variável/função que aparece na equação (isso ajuda a interpretar fisicamente os termos).
-
: fluxo de nêutrons (fluxo escalar). Fisicamente é a taxa de passagem de nêutrons por unidade de área por unidade de tempo; unidades típicas: nêutrons·cm ·s . Intuitivamente, está ligada à densidade de nêutrons por em média. -
: velocidade média dos nêutrons (cm·s ). -
: termo fonte volumétrica de nêutrons (nêutrons·cm ·s ); pode representar fonte externa ou a produção por fissão. -
: seção de choque macroscópica de absorção (unidade: cm ). Multiplicada por dá taxa volumétrica de remoção por absorção (nêutrons·cm ·s ). -
: coeficiente de difusão (unidade: cm). Relaciona gradientes de fluxo à corrente difusiva. -
(Laplaciano): operador que mede a curvatura espacial de ; unidades de 1/cm quando aplicado a .
Passo 3
Interpretação física termo a termo (didático, com unidades e intuição):
- Termo temporal
: representa a taxa de variação temporal do fluxo de nêutrons, convertido para a mesma escala das taxas volumétricas dos outros termos. Como tem unidade n·cm ·s e tem unidade cm·s , tem unidade n·cm ·s (mesma unidade de etc.). Fisicamente indica se o número de nêutrons está aumentando ou diminuindo no volume por unidade de tempo. Em regime estacionário esse termo é zero. - Termo fonte
: taxa de geração local de nêutrons por unidade de volume por unidade de tempo. Pode ser uma fonte externa (por exemplo, um feixe) ou a produção por fissão (quando , onde é o número médio de nêutrons liberados por fissão e é a seção macroscópica de fissão). - Termo de absorção
: representa a perda local de nêutrons devido à absorção (inclui fissão como processo de remoção se contado separadamente; atenção ao que está incluído nas seções). Unidade: n·cm ·s . É um termo de “sink”. - Termo de difusão
: vem de combinar a Lei de Fick (corrente de nêutrons ) com a conservação (vazamento). De fato . Esse termo representa o ganho ou perda local devido ao transporte por espalhamento entre vizinhanças (ou seja, vazamento para fora da região ou entrada de nêutrons vindos de regiões de alto fluxo). Intuitivamente, se é maior aqui do que na vizinhança, há perda por difusão; se é menor, há ganho. Dimensionalmente também dá n·cm ·s .
Passo 4
Breve resumo das hipóteses físicas implícitas na equação de difusão (útil para a parte b):
- Aproximação de difusão assume distribuição angular de nêutrons quase isotrópica (condição que surge quando espalhamento é muito mais provável que absorção e os gradientes espaciais são suaves).
- A energia dos nêutrons foi “integrada” ou tratada por grupos de energia; a equação escrita é tipicamente mono-energética ou para um grupo de energia.
- A equação é uma aproximação do transporte de Boltzmann, válida quando o comprimento livre médio é pequeno comparado às escalas de variação do fluxo.
Passo 5
Quando o modelo de difusão deixa de ser válido — lista com justificativa prática (vai direto aos pontos que costuma cair em prova):
- Meios com pouco espalhamento (muito absorvedores ou quase vácuo): se o coeficiente de espalhamento
não é muito maior que , a distribuição angular não é isotrópica e o modelo falha. Intuição: sem espalhamento suficiente, os nêutrons mantêm direcionalidade (transporte balístico), não se acomodam à hipótese de difusão. - Regiões próximas a superfícies/contornos onde ocorrem fortes anisotropias (dentro de algumas mean free paths da superfície): perto de uma interface com vácuo ou com material muito diferente, o fluxo angular fica anisotrópico e a aproximação difusiva quebra.
- Heterogeneidades finas comparáveis ao comprimento livre médio (combustível fintado, canôs, perfis pequenos): quando a escala geométrica de variação é da ordem do comprimento livre médio, a média local perde validade.
- Presença de fontes pontuais muito intensas ou feixes direcionados: a suposição de isotropia local cai.
- Reatores rápidos / nêutrons de alta energia onde espalhamento anisotrópico e as reações dependem fortemente da energia: a equação mono-energética de difusão não descreve bem o problema; é necessário tratamento por transporte multigrupo ou equação de transporte.
- Regiões vazias (voids) e meios óticos finos: difusão pressupõe muitas colisões locais; em voids não há colisões suficientes para estabelecer um comportamento difusivo.
Em resumo numérico: a difusão costuma valer bem quandoe as variações espaciais de ocorrem em comprimentos muito maiores que o comprimento livre médio. Fora desse regime, use-se a equação de transporte (Boltzmann).
Passo 6
Agora a parte (c) — derivação algébrica pedida, passo a passo e bem mastigada.
Hipóteses: estado estacionário (
) e fonte igual à produção por fissão
. Começamos da equação dada:
Em estado estacionário
, então o lado esquerdo é zero e obtemos
Substituímos
:
Agora substituímos a hipótese do buckling geométrico
(nota: isso corresponde à separação de variáveis na forma da moda fundamental; o sinal negativo indica curvatura do modo fundamental):
Ou seja
Rearranjando para colocar a produção
em evidência:
Definimos o fator de multiplicação local
como a razão entre produção e remoção (absorção + fuga). Assim
Esta é exatamente a expressão que o enunciado pediu para mostrar:
Observação prática e útil para prova: a
aparece em todos os termos e, se
, pode-se cancelar
e escrever a forma mais compacta
Mas o enunciado pediu a forma com
explícita, então mostramos a forma requisitada e indicamos a simplificação.
Passo 7
Explicação do que é
(buckling geométrico) e como ele depende do tamanho/forma:
-
é um parâmetro geométrico que captura a curvatura espacial do modo fundamental de fluxo. Valores maiores de correspondem a modos com curvatura maior (variações espaciais mais rápidas) — tipicamente ocorrem em geometrias menores. - Intuição: quanto menor o reator (mantendo mesma composição), maior será a fuga relativa de nêutrons por superfície; isso se traduz em maior
. Assim atua como termo de fuga efetiva (perda por vazamento). - Exemplos qualitativos (não precisa decorar para entender): reduzir dimensões aumenta
; aumentar o tamanho diminui .
Passo 8
Agora a análise qualitativa pedida em (d), usando a expressão
(forma simplificada que resulta do cancelamento de
). Vou analisar cada situação e justificar com base nos termos da equação.
(i) Reduzir o tamanho físico do reator (mesma composição):
- Ao reduzir o tamanho, a fuga aumenta →
aumenta. - Como
está no denominador, um aumento de torna o denominador maior → diminui. - Interpretação: o reator tende a ficar menos multiplicador (mais subcrítico) porque mais nêutrons escapam antes de causar fissão.
(ii) Inserção de barras de controle (material fortemente absorvedor):
- Barras de controle aumentam a absorção efetiva no núcleo, logo aumentam
(ou adicionam um termo local de absorção). - Um aumento de
aumenta o denominador → diminui. - Interpretação: absorvedores removem nêutrons que poderiam causar fissões, reduzindo a multiplicação.
(iii) Aumento do enriquecimento do combustível:
- Enriquecimento normalmente aumenta
(e possivelmente se não mudar muito), ou seja aumenta a produção por fissão. - Um aumento de
aumenta o numerador → aumenta. - Interpretação: mais fissões por unidade de fluxo significa maior multiplicação; reactor pode tornar-se mais crítico ou supercrítico se o aumento for suficiente.
Passo 9
Comentário final de prova e dica prática:
- Quando pedirem analisar efeitos qualitativos, pense sempre em termos de numerador (produção por fissão) e denominador (absorção + fuga). Aumentar produção →
. Aumentar absorção ou fuga → . - Se precisar calcular números, lembre-se de que
pode ser cancelada na expressão de se não for zero, o que facilita a conta.
Resposta
a) Termos:
— taxa temporal de variação do fluxo (nêutrons·cm
·s
);
— fonte volumétrica de nêutrons (externa ou por fissão);
— perda por absorção (sink) (
é a seção macroscópica de absorção);
— termo de difusão (vazamento/redistribuição espacial devido ao espalhamento), com
coeficiente de difusão e
a curvatura espacial do fluxo.
b) O modelo de difusão falha quando a hipótese de isotropia local não vale: meios com pouco espalhamento ( ), perto de superfícies/vácuo (dentro de algumas mean free paths), em heterogeneidades finas comparáveis ao comprimento livre médio, em voids, para fontes muito direcionais, ou em regimes de nêutrons rápidos/energia fortemente dependente — nesses casos é necessário usar a equação de transporte.
c) Em estado estacionário e com e
obtemos algebricamente:
b) O modelo de difusão falha quando a hipótese de isotropia local não vale: meios com pouco espalhamento (
c) Em estado estacionário e com
portanto
e, se
, simplifica para
.
d) (i) Diminuir o tamanho do reator → →
→ denominador
→
.
(ii) Inserir barras de controle (absorvedores) → → denominador
→
.
(iii) Aumentar o enriquecimento → → numerador
→
.
d) (i) Diminuir o tamanho do reator →
(ii) Inserir barras de controle (absorvedores) →
(iii) Aumentar o enriquecimento →