Termo
Termo
Termo
Y TERMOHIDRAULICA
3
ENERGIA LIBERADA EN LA FISIÓN
MeV
Total: ~ 205
4
ENERGIA PRODUCIDA EN LAS CAPTURAS
5
Barra de combustible para un PWR
de 900 MWe
6
ENERGIA LIBERADA EN LA FISION
7
TERMOHIDRAULICA
• El análisis térmico tanto del núcleo del reactor como del resto de componentes de la central.
• El núcleo del reactor se debe diseñar para producir la energía térmica deseada sin
sobrepasar las condiciones de diseño de temperaturas de los materiales.
• El fallo de dichos materiales por sobrepasar esos límites da lugar a fallo de componentes y a
la posible liberación de material radiactivo al refrigerante.
Rango de
operación
estacionaria
Retraso de la respuesta
de los operadores
Tiempo
Accidentes
• Un desequilibrio puede estar motivado por accidentes, los cuales casi siempre
dan lugar a desequilibrios en los balances térmicos.
• También se desean altas temperaturas del refrigerante a la salida del núcleo, ya que
esto da mejores rendimientos termodinámicos.
• La vaina sirve como contención para evitar el paso de los productos de fisión al refrigerante.
• La vaina casi nunca alcanza la fusión, a no ser que el flujo crítico sobrepase su valor de
diseño.
• Esta sometida a fuertes tensiones debidas al refrigerante, a los gases de fisión y al
hinchamiento del combustible.
• Estas tensiones son dependientes de la historia térmica, es decir, variaciones de la
temperatura por diversos niveles de potencia, y de la irradiación del combustible.
• En los LWR se limita el flujo de calor que puede transferirse de la vaina al refrigerante. Por
encima de ciertos valores, se forma una película de vapor en la vaina que impide dicha
transferencia, incrementándose la temperatura de la vaina y la posterior rotura. (Se vera)
• Reacción agua-circonio con producción de hidrógeno a partir de 1260 ºC. Límite de
temperatura en operación de 980 ºC, que coincide con 400 W/cm
Limitaciones en el refrigerante
• Aunque se desea que la temperatura sea lo mas alta posible a la salida, sin embargo el
aumento de la temperatura debe mantenerse en un rango para evitar ondas térmicas, como
podría ocurrir en una parada.
19
a. Experimentos termohidráulicos
Años 50:
Los años cincuenta presenciaron el comienzo de
la penetración de la energía nucleoeléctrica en el
mercado de la electricidad.
Fue un período de gran entusiasmo y de intensa
investigación y desarrollo, que ofreció esperanzas de
proporcionar al mundo una fuente sustitutiva de energía
barata y prácticamente inagotable.
En 1960 había 17 reactores de potencia en
explotación, con una capacidad total de generación de
electricidad de 1200 MW(e), en cuatro países: Francia,
la URSS, el Reino Unido y los Estados Unidos.
20
a. Experimentos termohidráulicos
Años 60:
Expertos en reactores nucleares cuestionaron el funcionamiento del ECCS
(Emergency Core Cooling System).
Un fallo del ECCS desencadena: una pérdida masiva de refrigerante, un fallo en la
contención y la fusión del núcleo.
Dadas la controversia y la falta de datos experimentales del ECCS en caso
de grandes roturas, se diseñó la instalación LOFT (Loss of Fluid Tests).
Como parte del diseño, se realizaron semiscale tests (tests de efectos
separados), que demostraron la falta de predictibilidad de los códigos existentes.
En consecuencia, se adoptaron los ECCS interim acceptance criteria,
recomendaciones muy conservadoras para simular los fenómenos termohidráulicos
en el diseño del ECCS de las nuevas CCNN.
21
a. Experimentos termohidráulicos
Años 70:
Tras un periodo de incertidumbre provocada por la desconfianza hacia los
ECCS interim acceptance criteria, en 1978 se puso en marcha la instalación LOFT
(Idaho, EEUU), con combustible real y una escala similar a la de las plantas reales.
Se realizaron 8 experimentos, 6 de ellos termohidráulicos (oecd-nea.org/tools/ie/):
- LP-FW-1: Loss-of-feedwater, primary feed and bleed recovery procedure.
- LP-SB-1: Hot leg SB LOCA, early pump trip.
- LP-SB-2: Hot leg SB LOCA, delayed pump trip.
- LP-SB-3: Cold leg SB LOCA, core uncovery, secondary feed and bleed recovery procedure,
accumulator injection at low-pressure differential.
- LP-02-6: 200 % large-break LOCA, US licensing case.
- LP-LB-1: 200 % large-break LOCA, UK licensing case.
22
a. Experimentos termohidráulicos
Años 70:
Además de LOFT, otras
instalaciones experimentales
sirvieron para estudiar
fenómenos individuales y
fenómenos integrales y
desarrollar códigos
termohidráulicos.
En las instalaciones de
BETHSY (CEA, Grenoble, Francia)
se simularon la mayoría de los
accidentes más relevantes de los
reactores PWR.
23
b. Desarrollo de códigos TH
Años 70 y 80:
En este periodo se desarrollaron diferentes códigos termohidráulicos (TH),
entre los cuales el más popular en España fue RELAP5. Ésta fue su evolución:
1. RELAPSE-1: PWR con presionador, volumen caliente y volumen frío en circuito 1.
2. RELAP2: versión anterior + BWR con separación de burbujas en dos fases.
3. RELAP3: 20 volúmenes de control y 50 uniones definidas por el usuario +
modelo de conducción de varilla de combustible + correlaciones de la
transferencia de calor.
4. RELAP4: mejoras en las ecs. de conservación y el tratamiento numérico + ECCS.
5. RELAP5: 6 ecs. de conservación + formato libre de entrada de datos + obtención
de gráficos + flexibilidad de simulación + NPA (Nuclear Plant Analyser).
Veamos… https://www.youtube.com/watch?v=Z4BBMW-4KVg
24
b. Desarrollo de códigos TH
Etapas:
25
c. Desarrollo de modelos de CCNN
El modelo de central nuclear:
Es el modelo completo de planta en el que se integran tres partes:
1. La termohidráulica:
Resuelve ecuaciones de balance mediante elementos básicos: volúmenes,
uniones, válvulas, bombas, conexiones laterales, estructuras de calor y otros
componentes (separadores, acumuladores, presionador, etc.).
2. El reactor:
Calcula la generación térmica en el combustible, incluyendo: potencia de
fisión, calor de decaimiento, oxidación de las vainas y su distribución espacial.
3. La instrumentación y el control:
Define los mecanismos de procesamiento electrónico que modulan los
requisitos para activar cualquier acción automática.
Contiene 1.500 ecs. o 15.000 – 20.000 registros en RELAP5 efecto usuario.
26
c. Desarrollo de modelos de CCNN
Research on thermal hydraulic behavior of SB LOCAs (small break
loss of coolant accident) in AP1000 (WEST-PWR):
Años 90 - 2000:
Tras la diversificación de los códigos surgen
movimientos de unificación de códigos TH
europeos y americanos, dando lugar al código
consolidado USRNC; y de éste con otros códigos,
como los códigos neutrónicos y de accidentes
severos.
Además, se amplían los modelos para los
diseños de reactores avanzados y se desarrollan
interfaces para facilitar y acelerar su uso.
Flujo de calor (transferencia): Cantidad de calor por unidad de tiempo, lo que significa que se mide en
J/s o W.
𝐶
q=
∆𝑡
Densidad de flujo de calor: flujo de calor por unidad de superficie, se mide en W/m2, q’’
Modos de transferencia de calor:
- A través de la materia (conducción / difusión)
- Por flujo de fluido (convección)
- Radiación
Difusión (Conducción):
Se produce a nivel atómico, la energía se trasfiere por la vibración de átomos y electrones.
La transferencia de calor a través de una pared homogénea se
describe macroscópicamente como (Ley de Fourier):
𝑘 𝐴 ∆𝑇
q=
∆𝑥
Conductividad
térmica
(W/m2K/m=W/mK)
El perfil de temperatura muestra el gradiente de esta durante la conducción de calor a
través de una pared; es más brusca en los malos conductores
Convección
𝑞 = ℎ𝐴Δ𝑇
coeficiente de convección (W/m2K),
Convección natural: El fluido se desplaza por diferencias en densidades del fluido causadas
por diferencias en temperatura o cambios de fase
P.e. una capa caliente de fluido gana flotabilidad y asciende, es sustituida por el fluido circundante más frío
que, a continuación, se calienta de nuevo, se dispersa y se enfría. El proceso continúa mientras haya una
fuente de calor y un disipador.
Ejemplo de convección natural por
diferencias en temperatura:
Emision de energía interna en forma de r.e.m. (no precisa materia para desplazarse).
q=ε𝜎𝐴𝑇 4
Ejemplo:
Si la vasija se queda sin agua (LOCA) la radiación podría ser la única forma de perdida del calor de
decaimiento, esto podría conducir a altas temperaturas en el núcleo y a su fusión.
Transferencia de calor con cambio de fase
∆𝑇 𝑞 = ℎ𝐴Δ𝑇
𝑞=𝑘𝐴
∆𝑥
• Pastilla, conducción
𝑞 ′′′ 𝑟02 𝑞′
𝑇𝐶𝐿 − 𝑇𝐹 = =
4𝑘𝐹 4𝜋𝑘𝐹
𝑞′' = ℎ𝑒𝑓 𝑇𝐹 − 𝑇𝐶
𝑞 ′′ 𝑟0𝑞′′′
r0 𝑇𝐹 − 𝑇𝐶 = =
ℎ 𝑒𝑓 2ℎ
𝑒
• Vaina
𝑑𝑇
𝑞 ′′ = −𝑘𝑉 𝛻𝑇= −𝑘 𝑉 𝐿𝑒𝑦 𝑑𝑒 𝐹𝑜𝑢𝑟𝑖𝑒𝑟
𝑑𝑟
𝑞′ 𝑟1 𝑞′ 𝑟2 𝑟02 𝑞 ′′′ 𝑟2
𝑇𝐶−𝑇𝑆= − 𝑙𝑛 = 𝑙𝑛 = 𝑙𝑛
2𝜋𝑘𝑉 𝑟2 2𝜋𝑘𝑉 𝑟1 2𝑘𝑉 𝑟1
• Refrigerante
q’
𝑞 ′′ = ℎ𝐶 𝑇𝑆 − 𝑇𝐹𝑙 ⇒ = ℎ𝐶 𝑇𝑆 − 𝑇𝐹𝑙
2𝜋𝑟2
𝑞′ 𝑞 ′′′ 𝑟02
𝑇𝑆 − 𝑇𝐹𝑙 = =
2𝜋𝑟2ℎ𝐶 2𝑟2ℎ𝐶
L Martes Miercoles
24/04/2023 6 -7 CH 4:30-6 CH
01/05/2023 4:30-5:30 MH
08/05/2023 6 -8:30 MH + CH
22/05/2023 4:30 -6 SR
Convección
Conducción Resumen de transferencia
𝑞 ′′′ 𝑟02 𝑞′
𝑇𝐶𝐿 − 𝑇𝐹 = =
4𝑘𝐹 4𝜋𝑘𝐹
𝑞 ′′ 𝑟 0 𝑞 ′′′ 𝑞′
𝑇𝐹 − 𝑇𝐶 = = =
ℎ𝑒 2ℎ 2𝜋𝑟 ℎ
𝑒𝑓 0 𝑒𝑓
𝑟02 𝑞 ′′′
𝑟2 𝑞′ 𝑟2
𝑇𝐶 −𝑇𝑆 = 𝑙𝑛 = 𝑙𝑛
2𝑘𝑉 𝑟1 2𝜋𝑘𝑉 𝑟1
𝑞 ′′′ 𝑟02 𝑞′
𝑇𝑆 − 𝑇𝐹𝑙 = =
r0 2𝑟2ℎ𝐶 2𝜋𝑟2ℎ𝐶
r1
r2
𝑝𝑜𝑡𝑒𝑛𝑐𝑖𝑎 𝑡𝑒𝑟𝑚𝑖𝑐𝑎
𝑞 ′′′ =
𝑣𝑜𝑙𝑢𝑚𝑒𝑛 𝑑𝑒 𝑐𝑜𝑚𝑏𝑢𝑠𝑡𝑖𝑏𝑙𝑒
Perfil de temperatura
𝑞′′ = ℎ(𝑇𝑤 − 𝑇𝑏 )
Tª del fluido
Coef. de convección Tª de la vaina
Transmisión de calor en pared
Temperatura de pared, 𝑇𝑤
Es la temperatura de la pared
Regímenes exterior de la vaina
ℎ𝐷ℎ 𝑐𝑜𝑛𝑣𝑒𝑐𝑐𝑖ó𝑛 𝑘
𝑁𝑢 = = → ℎ = 𝑁𝑢 ·
𝑘 𝑐𝑜𝑛𝑑𝑢𝑐𝑐𝑖ó𝑛 𝐷ℎ
1. Régimen de flujo
2. Velocidad del flujo
3. Propiedades físicas
4. Diámetro hidráulico (Dh)
𝑁𝑢𝐿𝑊𝑅 = 35 − 220, 𝑁𝑢 𝑆𝐹𝑅 = 5 − 7
Termohidráulica Nuclear (THN) | Máster Universitario en Ciencia y Tecnología Nuclear | Universidad Politécnica de Madrid 56
Transmisión de calor en pared
Temperatura de pared, 𝑇𝑤
Es la temperatura de la pared
Regímenes exterior de la vaina
Termohidráulica Nuclear (THN) | Máster Universitario en Ciencia y Tecnología Nuclear | Universidad Politécnica de Madrid 57
Transmisión de calor en pared
Regímenes
Termohidráulica Nuclear (THN) | Máster Universitario en Ciencia y Tecnología Nuclear | Universidad Politécnica de Madrid 58
Transmisión de calor en pared
Nucleación
Para que comience la ebullición es preciso que ocurra el fenómeno de
nucleación que dará pie a la ebullición nucleada: crecimiento de las
burbujas que inician el mecanismo de TDC por ebullición. 𝑃𝐿
Nucleación homogénea 𝜎
𝑃𝑉
Se produce sobre superficies lisas sin porosidades cuando un núcleo
de vapor de radio 𝑟 ∗ adquiere la suficiente energía como para que la
presión 𝑃𝑉 del vapor logre vencer las fuerzas de tensión superficial 𝜎 y la 𝑟∗
presión del líquido 𝑃𝐿
2𝜎
𝜋𝑟 ∗2 𝑃𝑉 = 𝜋𝑟 ∗2 𝑃𝐿 + 2𝜋𝑟 ∗ 𝜎 → 𝑃𝑉 − 𝑃𝐿 = ∗
𝑟
𝑟 ∗ es del orden de nm. La temperatura en la superficie debe alcanzar un
valor muy superior respecto a 𝑇𝑠𝑎𝑡 , hasta 220 ºC (frente a 100 ºC).
Termohidráulica Nuclear (THN) | Máster Universitario en Ciencia y Tecnología Nuclear | Universidad Politécnica de Madrid 59
Transmisión de calor en pared
Nucleación
Para que comience la ebullición es preciso que ocurra el fenómeno
de nucleación que dará pie a la ebullición nucleada: crecimiento de
las burbujas que inician el mecanismo de TDC por ebullición.
Nucleación heterogénea
𝑃𝐿
La nucleación heterogénea ocurre debido a la presencia de
microcavidades que almacenan gas, constituyendo los centros
de nucleación. La presión de ese gas facilita la ebullición. 𝜎
2𝜎
𝜋𝑟 ∗2 (𝑃𝑉 + 𝑃𝑔𝑎𝑠 ) = 𝜋𝑟 ∗2 𝑃𝐿 + 2𝜋𝑟 ∗ 𝜎 → 𝑃𝐵 − 𝑃𝐿 = ∗
𝑟 𝑃𝑉 + 𝑃𝑔
𝑃
𝑇𝑤 − 𝑇𝑠𝑎𝑡 = 25 exp − 𝑞 ′′0.25
6.2
𝑃 es la presión (MPa)
Termohidráulica Nuclear (THN) | Máster Universitario en Ciencia y Tecnología Nuclear | Universidad Politécnica de Madrid 62
Transmisión de calor en pared
Regímenes
63
Transmisión de calor en pared
Regímenes
Termohidráulica Nuclear (THN) | Máster Universitario en Ciencia y Tecnología Nuclear | Universidad Politécnica de Madrid 64
Transmisión de calor en pared
Regímenes 𝑇𝑠𝑎𝑡
𝑇𝑠𝑎𝑡
Ebullición nucleada saturada
𝑇𝑏 𝑇𝑠𝑎𝑡
• Cuando aumenta el tamaño de las burbujas producidas, estas
generan mayor agitación en la superficie caliente
𝑄𝑠𝑒𝑛
𝑇𝑠𝑎𝑡
• Esto hace que aumente la eficiencia del mecanismo de TDC en la
superficie
𝑇𝑠𝑎𝑡
• El bulk solo se calienta por la convección del fluido, ya que el calor
latente que transportan las burbujas escapa del subcanal con el vapor 𝑄𝑙𝑎𝑡
𝑇𝑤
𝑄𝑠𝑒𝑛
Termohidráulica Nuclear (THN) | Máster Universitario en Ciencia y Tecnología Nuclear | Universidad Politécnica de Madrid 65
Transmisión de calor en pared
Regímenes
Termohidráulica Nuclear (THN) | Máster Universitario en Ciencia y Tecnología Nuclear | Universidad Politécnica de Madrid 66
Transmisión de calor en pared
Regímenes
Termohidráulica Nuclear (THN) | Máster Universitario en Ciencia y Tecnología Nuclear | Universidad Politécnica de Madrid 67
Transmisión de calor en pared
Regímenes
Correlación de Chen (convección forzada y capa líquida) 𝑞 ′′ = ℎ2Φ 𝑇𝑤 − 𝑇𝑠𝑎𝑡 → ℎ2Φ = ℎ𝑁𝐵 + ℎ𝐶
𝒉𝑵𝑩 es el coeficiente de transmisión de calor por ebullición nucleada
Termohidráulica Nuclear (THN) | Máster Universitario en Ciencia y Tecnología Nuclear | Universidad Politécnica de Madrid 68
Transmisión de calor en pared
Regímenes
Correlación de Chen (convección forzada y capa líquida) 𝑞 ′′ = ℎ2Φ 𝑇𝑤 − 𝑇𝑠𝑎𝑡 → ℎ2Φ = ℎ𝑁𝐵 + ℎ𝐶
𝒉𝑪 es el coeficiente de transmisión de calor por convección forzada
Termohidráulica Nuclear (THN) | Máster Universitario en Ciencia y Tecnología Nuclear | Universidad Politécnica de Madrid 69
Transmisión de calor en pared
Regímenes
Termohidráulica Nuclear (THN) | Máster Universitario en Ciencia y Tecnología Nuclear | Universidad Politécnica de Madrid 70
Transmisión de calor en pared
Regímenes
Termohidráulica Nuclear (THN) | Máster Universitario en Ciencia y Tecnología Nuclear | Universidad Politécnica de Madrid 71
Transmisión de calor en pared
Regímenes
Ebullición en película
Ebullición en película
Film
q’’ boiling
Convección natural
- Natural convection
- Occurs at a small temperature
difference
Hasta (TW - TS) = 5ºC
Transferencia de calor en pared – curva de ebullición
Film
Ebullición nucleada q’’ boiling
q’’ Film
boiling
Ebullición nucleada
q’’ Film
boiling
Ebullición parcial
q’’ Film
boiling
Ebullición en película y radiación
Punto c:
q’’ Film
“Crisis de ebullición o Departure boiling
from Nucleate Boiling (DNB)”.
La TW aumenta considerablemente.
DNB y DNBR
En operación, DNBR esta alrededor de 2.0 y dependerá de los parámetros del sistema primario de
refrigeración.
En un PWR tipico, DNBR debe de ser siempre >1.36. Lo que asegura con un 95% de probabilidad que
no ocurrirá una crisis de ebullición
Calculo de CHF/DNB
• Estudios detallados del fenómeno han dado lugar a modelos mecanísticos que pueden proporcionar un mejor
entendimiento del proceso de flujo crítico de calor.
• Sin embargo, para diseño se siguen usando correlaciones empíricas que tienen diversos inconvenientes:
1. Solo son válidas para el rango de parámetros para el que han sido desarrolladas.
La extrapolación da lugar a importantes errores.
Resumen
Los regímenes de flujo y los modos de transmisión de calor están muy interrelacionados, pero no se
relacionan uno a uno.
Existen modelos o correlaciones que permiten determinar el valor del coeficiente de transmisión de
calor de la vaina al refrigerante y que son dependientes del régimen de transmisión de calor.
El concepto de flujo crítico de calor (CHF) es de gran importancia en seguridad nuclear ya que de
alcanzarse se podría dar degradación de la vaina y otros fenómenos derivados de la elevación brusca de
la temperatura de esta.
Los valores de flujo crítico de calor se determinan también con modelos empíricos que dependen
principalmente del tipo de reactor y de los parámetros que tienen influencia en los modos de transmisión de
calor que se da en este.
Relación con los límites de seguridad del núcleo del reactor
• Durante el funcionamiento normal y durante los incidentes, debe mantenerse la integridad del
núcleo del reactor (que las pastillas de combustible no fusionen y que el revestimiento no
sufra daños).
• La integridad del núcleo en caso de funcionamiento normal y de incidentes moderadamente
frecuentes está garantizada por los límites de diseño que relacionan la temperatura de las
pastillas de combustible y la crisis de ebullición.
• Los límites serían relativamente fáciles de definir si la distribución de energía en el núcleo
fuera uniforme. Desgraciadamente, no es así.
• La distribución de la potencia es desigual y depende de los parámetros nucleares,
termodinámicos y de fabricación,
• Un enfoque práctico consiste en describir la distribución de potencia multiplicando los valores
medios por factores de pico.
Relación con los límites de seguridad del núcleo del reactor
Densidad de potencia lineal media (q’): potencia térmica liberada por unidad de longitud de la
barra de combustible (kW/m) Cociente entre la potencia térmica nominal del reactor y la
longitud total de todas las varillas de combustible.
Esta densidad no es igual para todas las varillas de un elemento combustible. Ni entre elementos
individuales y ni con la altura de una varilla de combustible en el núcleo.
Los F se calculan a partir de mediciones de flujo de neutrones térmicos que se realizan, por regla
general, cada 1000 MWd/t de quemado del núcleo, lo que equivale aproximadamente a una vez al
mes en un PWR típico.
Los límites operativos de los F no se traspasan siempre que el reactor funcione de forma que:
* la diferencia de flujo axial (la diferencia entre la potencia de la mitad superior e inferior del núcleo)
esté dentro de la banda operacional;
* se retiren diferentes grupos de barras de control con un solapamiento programado,
* las barras de control están por encima del límite de inserción,
* las barras de control individuales de un grupo se retiran a la misma altura; en ningún caso deben
estar separadas más de un cierto número de pasos.
Relación con los límites de seguridad del núcleo del reactor