Articulo de referencia

criterio de Lawson

Criterio de Lawson para importantes experimentos de fusión por confinamiento magnético El criterio de Lawson es un indicador de rendimiento utilizado en la investigación de la f...

Criterio de Lawson para importantes experimentos de fusión por confinamiento magnético

El criterio de Lawson es un indicador de rendimiento utilizado en la investigación de la fusión nuclear . Compara la tasa de energía generada por las reacciones de fusión dentro del combustible con la tasa de pérdidas de energía al medio ambiente. Cuando la tasa de producción es mayor que la tasa de pérdida, el sistema produce energía neta. Si el combustible captura suficiente energía, el sistema se vuelve autosostenible y se dice que se ha encendido .

El concepto fue desarrollado por primera vez por John D. Lawson en un documento clasificado de 1955 [ 1 ] que fue desclasificado y publicado en 1957. [ 2 ] Tal como fue formulado originalmente, el criterio de Lawson da un valor mínimo requerido para el producto de la densidad del plasma (electrón) n e y el " tiempo de confinamiento de energía ".τmi{\displaystyle \tau _{E}}eso conduce a una producción neta de energía.

Análisis posteriores sugirieron que un indicador de mérito más útil es el producto triple de la densidad, el tiempo de confinamiento y la temperatura del plasma T. Este producto triple también tiene un valor mínimo requerido, y el nombre "criterio de Lawson" podría referirse a dicho valor.

El 8 de agosto de 2021, investigadores de la Instalación Nacional de Ignición del Laboratorio Nacional Lawrence Livermore en California confirmaron haber producido la primera ignición exitosa de una reacción de fusión nuclear que superó los criterios de Lawson en el experimento. [ 3 ] [ 4 ]

Balance energético

El concepto central del criterio de Lawson consiste en examinar el balance energético de cualquier central de fusión que utilice plasma caliente. Esto se muestra a continuación:

Potencia neta = Eficiencia × (Fusión − Pérdida por radiación − Pérdida por conducción)

  1. La potencia neta es el exceso de energía que supera la necesaria internamente para que el proceso se lleve a cabo en cualquier central de fusión nuclear.
  2. La eficiencia se refiere a la cantidad de energía necesaria para hacer funcionar el dispositivo y a la eficacia con la que este recolecta energía de las reacciones.
  3. La fusión es la tasa de energía generada por las reacciones de fusión.
  4. La pérdida por radiación es la energía que se pierde en forma de luz (incluidos los rayos X ) que sale del plasma.
  5. La pérdida por conducción es la energía que se pierde cuando las partículas abandonan el plasma, llevándose consigo esa energía.

Lawson calculó la tasa de fusión suponiendo que el reactor de fusión contiene una nube de plasma caliente con una curva gaussiana de energías de partículas individuales, una distribución de Maxwell-Boltzmann caracterizada por la temperatura del plasma. Basándose en esta suposición, estimó el primer término, la energía de fusión producida, utilizando la ecuación de fusión volumétrica. [ 5 ]

Fusión = Densidad numérica del combustible A × Densidad numérica del combustible B × Sección transversal (Temperatura) × Energía por reacción

  1. La fusión es la tasa de energía de fusión producida por el plasma.
  2. La densidad numérica es la densidad de partículas por unidad de volumen de los combustibles respectivos (o de un solo combustible, en algunos casos).
  3. La sección transversal es una medida de la probabilidad de un evento de fusión, que se basa en la temperatura del plasma.
  4. La energía por reacción es la energía liberada en cada reacción de fusión.

Esta ecuación se suele promediar sobre una población de iones con distribución normal . El resultado es la cantidad de energía que genera el plasma en cualquier instante.

Lawson luego estimó [ 5 ] las pérdidas por radiación utilizando la siguiente ecuación:

PAGB=1.4×1034norte2T1/2Wdometro3{\displaystyle P_{B}=1.4\times 10^{-34}\cdot n^{2}\cdot T^{1/2}\;{\frac {\mathrm {W} }{\mathrm {cm} ^{3}}}}

donde n es la densidad numérica de la nube y T es la temperatura . Para su análisis, Lawson ignora las pérdidas por conducción. En realidad, esto es casi imposible; prácticamente todos los sistemas pierden energía a través del transporte de partículas, donde la masa abandona el plasma y se lleva consigo su energía.

Al igualar las pérdidas por radiación y las tasas de fusión volumétricas, Lawson estimó la temperatura mínima para la fusión en la reacción deuterio - tritio (DT).

12D+13T24Hmi(3.5METROmiV)+01norte(14.1METROmiV){\displaystyle _{1}^{2}\mathrm {D} +\,_{1}^{3}\mathrm {T} \rightarrow \,_{2}^{4}\mathrm {He} \left(3.5\,\mathrm {MeV} \right)+\,_{0}^{1}\mathrm {n} \left(14.1\,\mathrm {MeV} \right)}

para ser de 30 millones de grados (2,6 keV), y para la reacción deuterio - deuterio (DD)

12D+12D13T(1.0METROmiV)+11pag(3.0METROmiV){\displaystyle _{1}^{2}\mathrm {D} +\,_{1}^{2}\mathrm {D} \rightarrow \,_{1}^{3}\mathrm {T} \left(1.0\,\mathrm {MeV} \right)+\,_{1}^{1}\mathrm {p} \left(3.0\,\mathrm {MeV} \right)}

ser 150 millones de grados (12,9 keV). [ 2 ] [ 6 ]

Extensiones en n τ E

El tiempo de confinamientoτmi{\displaystyle \tau _{E}}mide la tasa a la que un sistema pierde energía a su entorno. Cuanto más rápida sea la tasa de pérdida de energía,PAGloss{\displaystyle P_{\mathrm {loss} }}Cuanto menor sea el tiempo de confinamiento de energía, menor será la densidad de energía.W{\displaystyle W}(contenido energético por unidad de volumen) dividido por la densidad de pérdida de potenciaPAGloss{\displaystyle P_{\mathrm {loss} }}(tasa de pérdida de energía por unidad de volumen):

τmi=WPAGloss{\displaystyle \tau _{E}={\frac {W}{P_{\mathrm {loss} }}}}

Para que un reactor de fusión funcione en estado estacionario, el plasma de fusión debe mantenerse a una temperatura constante. Por lo tanto, es necesario añadir energía térmica al mismo ritmo que el plasma la pierde para mantener las condiciones de fusión. Esta energía puede provenir de las propias reacciones de fusión, según el tipo de reacción, o bien, de un calentamiento adicional mediante diversos métodos.

A modo de ejemplo, aquí se derivará el criterio de Lawson para la reacción DT, pero el mismo principio puede aplicarse a otros combustibles de fusión. También se asumirá que todas las especies tienen la misma temperatura, que no hay iones presentes aparte de los iones del combustible (sin impurezas ni cenizas de helio), y que D y T están presentes en la mezcla óptima 50-50. La densidad iónica es entonces igual a la densidad electrónica y la densidad de energía de electrones e iones juntos viene dada, según la ley de los gases ideales , por

W=3norteT{\displaystyle W=3nT}

dóndeT{\displaystyle T}es la temperatura en electronvoltios (eV) ynorte{\displaystyle n}es la densidad de partículas.

La tasa de volumenF{\displaystyle f}(reacciones por volumen por tiempo) de las reacciones de fusión es

F=nortednortetσv=14norte2σv{\displaystyle f=n_{\mathrm {d} }n_{\mathrm {t} }\langle \sigma v\rangle ={\frac {1}{4}}n^{2}\langle \sigma v\rangle }

dóndeσ{\displaystyle \sigma }es la sección transversal de fusión ,v{\displaystyle v}es la velocidad relativa y{\displaystyle \langle \rangle }denota un promedio sobre la distribución de velocidad maxwelliana a la temperaturaT{\displaystyle T}.

La tasa volumétrica de calentamiento por fusión esF{\displaystyle f}vecesmidoh{\displaystyle E_{\mathrm {ch} }}, la energía de los productos de fusión cargados (los neutrones no pueden ayudar a calentar el plasma). En el caso de la reacción DT,midoh=3.5METROmiV{\displaystyle E_{\mathrm {ch} }=3.5\,\mathrm {MeV} }.

El criterio de Lawson, o valor mínimo de (densidad electrónica * tiempo de confinamiento de energía) requerido para el auto calentamiento, para tres reacciones de fusión. Para DT, nτ E se minimiza cerca de la temperatura de 25 keV (300 millones de kelvin).

El criterio de Lawson requiere que el calentamiento por fusión supere las pérdidas:

FmidohPAGloss{\displaystyle fE_{\rm {ch}}\geq P_{\rm {loss}}}

Sustituyendo cantidades conocidas se obtiene:

14norte2σvmidoh3norteTτmi{\displaystyle {\frac {1}{4}}n^{2}\langle \sigma v\rangle E_{\rm {ch}}\geq {\frac {3nT}{\tau _{E}}}}

Reordenando la ecuación se obtiene:

La cantidadT/σv{\displaystyle T/\langle \sigma v\rangle }es una función de la temperatura con un mínimo absoluto. Reemplazar la función con su valor mínimo proporciona un límite inferior absoluto para el producto.norteτmi{\displaystyle n\tau _{E}}Este es el criterio de Lawson.

Para la reacción deuterio - tritio , el valor físico es al menos

norteτmi1.51020smetro3{\displaystyle n\tau _{E}\geq 1.5\cdot 10^{20}{\frac {\mathrm {s} }{\mathrm {m} ^{3}}}}

El mínimo del producto se produce cercaT=26kmiV{\displaystyle T=26\,\mathrm {keV} }.

Extensión hacia el "producto triple"

Una figura de mérito aún más útil es el "producto triple" de densidad, temperatura y tiempo de confinamiento, nT τ E . Para la mayoría de los conceptos de confinamiento, ya sean inerciales , de espejo o toroidales, la densidad y la temperatura pueden variar en un rango bastante amplio, pero la presión máxima alcanzable p es constante. En tal caso, la densidad de potencia de fusión es proporcional a p 2v >/ T  2 . Por lo tanto, la potencia de fusión máxima disponible de una máquina dada se alcanza a la temperatura T donde <σ v >/ T  2 es máxima. Siguiendo la derivación anterior, se obtiene fácilmente la siguiente desigualdad:

norteTτmi12midohT2σv{\displaystyle nT\tau _{\rm {E}}\geq {\frac {12}{E_{\rm {ch}}}}\,{\frac {T^{2}}{\langle \sigma v\rangle }}}
La condición de triple producto de fusión para tres reacciones de fusión

La cantidadT2σv{\displaystyle {\frac {T^{2}}{\langle \sigma v\rangle }}}También es una función de la temperatura con un mínimo absoluto a una temperatura ligeramente inferior a laTσv{\displaystyle {\frac {T}{\langle \sigma v\rangle }}}.

Para la reacción DT, el mínimo se produce en T  =  14  keV. El promedio <σ v > en esta región de temperatura se puede aproximar como [ 7 ]

σv=1.11024T2metro3s,TinortekmiV,{\displaystyle \left\langle \sigma v\right\rangle =1.1\cdot 10^{-24}T^{2}\;{\frac {{\rm {m}}^{3}}{\rm {s}}}\,{\rm {,}}\quad {\rm {T\,in\,keV}}{\rm {,}}}

por lo que el valor mínimo del valor del producto triple en T  =  14  keV es aproximadamente

norteTτmi12142kmiV21.11024metro3s1423500kmiV31021keV s/metro3(3.51028K s/metro3){\displaystyle {\begin{matrix}nT\tau _{E}&\geq &{\frac {12\cdot 14^{2}\cdot {\rm {keV}}^{2}}{1.1\cdot 10^{-24}{\frac {{\rm {m}}^{3}}{\rm {s}}}14^{2}\cdot 3500\cdot {\rm {keV}}}}\approx 3\cdot 10^{21}{\mbox{keV s}}/{\mbox{m}}^{3}\\\end{matrix}}(3.5\cdot 10^{28}{\mbox{K s}}/{\mbox{m}}^{3})}

Esta cifra aún no se ha alcanzado en ningún reactor, aunque las últimas generaciones de máquinas se han acercado. JT-60 reportó 1,53 × 10²¹ keV · sm⁻³ . [ 8 ] Por ejemplo, el TFTR ha alcanzado las densidades y vidas útiles energéticas necesarias para lograr Lawson a las temperaturas que puede generar, pero no puede generar esas temperaturas al mismo tiempo. ITER pretende hacer ambas cosas.

En cuanto a los tokamaks , existe una motivación especial para utilizar el triple producto. Empíricamente, se ha encontrado que el tiempo de confinamiento de energía τ E es casi proporcional a n 1/3 / P  2/3 . En un plasma encendido cerca de la temperatura óptima, la potencia de calentamiento P es igual a la potencia de fusión y, por lo tanto, es proporcional a n 2 T 2 . El triple producto escala como 

norteTτminorteT(norte1/3/PAG2/3)norteT(norte1/3/(norte2T2)2/3)T1/3{\displaystyle {\begin{matrix}nT\tau _{E}&\propto &nT\left(n^{1/3}/P^{2/3}\right)\\&\propto &nT\left(n^{1/3}/\left(n^{2}T^{2}\right)^{2/3}\right)\\&\propto &T^{-1/3}\\\end{matrix}}}

El producto triple depende solo débilmente de la temperatura como T  -1/3 . Esto hace que el producto triple sea una medida adecuada de la eficiencia del esquema de confinamiento.

Confinamiento inercial

El criterio de Lawson se aplica tanto a la fusión por confinamiento inercial (FCI) como a la fusión por confinamiento magnético (FCM), pero en el caso inercial resulta más útil expresarlo de una forma diferente. Una buena aproximación para el tiempo de confinamiento inercialτmi{\displaystyle \tau _{E}}es el tiempo que tarda un ion en recorrer una distancia R a su velocidad térmica.

vth=kBTmetroi{\displaystyle v_{th}={\sqrt {\frac {k_{\rm {B}}T}{m_{i}}}}}

donde m i denota la masa iónica media. El tiempo de confinamiento inercialτmi{\displaystyle \tau _{E}}por lo tanto, se puede aproximar como

τmiRvth=RkBTmetroi=RmetroikBT .{\displaystyle {\begin{matrix}\tau _{E}&\approx &{\frac {R}{v_{th}}}\\\\&=&{\frac {R}{\sqrt {\frac {k_{\rm {B}}T}{m_{i}}}}}\\\\&=&R\cdot {\sqrt {\frac {m_{i}}{k_{\rm {B}}T}}}{\mbox{ .}}\\\end{matrix}}}

Sustituyendo la expresión anterior en la relación ( 1 ), obtenemos

norteτminorteRmetroikBT12midohkBTσvnorteR12midoh(kBT)3/2σvmetroi1/2norteR(kBT)3/2σv .{\displaystyle {\begin{matrix}n\tau _{E}&\approx &n\cdot R\cdot {\sqrt {\frac {m_{i}}{k_{B}T}}}\geq {\frac {12}{E_{\rm {ch}}}}\,{\frac {k_{\rm {B}}T}{\langle \sigma v\rangle }}\\\\n\cdot R&\gtrapprox &{\frac {12}{E_{\rm {ch}}}}\,{\frac {\left(k_{\rm {B}}T\right)^{3/2}}{\langle \sigma v\rangle \cdot m_{i}^{1/2}}}\\\\n\cdot R&\gtrapprox &{\frac {\left(k_{\rm {B}}T\right)^{3/2}}{\langle \sigma v\rangle }}{\mbox{ .}}\\\end{matrix}}}

Este producto debe ser mayor que un valor relacionado con el mínimo de T  3/2 /<σv>. El mismo requisito se expresa tradicionalmente en términos de densidad de masa ρ  =  < nm i >:

ρR1gramo/dometro2{\displaystyle \rho \cdot R\geq 1\mathrm {g} /\mathrm {cm} ^{2}}

La satisfacción de este criterio a la densidad del DT sólido (0,2 g/cm³ ) requeriría un pulso láser de energía implausiblemente grande. Suponiendo que la energía requerida escala con la masa del plasma de fusión ( E láser ~ ρR³ ~ ρ⁻² ), comprimir el combustible a 10³ o 10⁴ veces la densidad sólida reduciría la energía requerida en un factor de 10⁶ o 10⁸ , situándola en un rango realista. Con una compresión de 10³ , la densidad comprimida será de 200 g/cm³ , y el radio comprimido puede ser tan pequeño como 0,05 mm. El radio del combustible antes de la compresión sería de 0,5 mm. El pellet inicial será quizás el doble de grande, ya que la mayor parte de la masa se ablacionará durante la compresión.  

La potencia de fusión multiplicada por la densidad es un buen parámetro para determinar la temperatura óptima para el confinamiento magnético, pero para el confinamiento inercial, la fracción de quemado del combustible probablemente sea más útil. El quemado debería ser proporcional a la tasa de reacción específica ( < σv > ) multiplicada por el tiempo de confinamiento (que escala como T⁻¹ ) dividido por la densidad de partículas n : 

fracción de combustión norte2σvT1/2/norte(norteT)σv/T3/2{\displaystyle {\begin{matrix}{\mbox{burn-up fraction }}&\propto &n^{2}\langle \sigma v\rangle T^{-1/2}/n\\&\propto &\left(nT\right)\langle \sigma v\rangle /T^{3/2}\\\end{matrix}}}

Por lo tanto, la temperatura óptima para la fusión por confinamiento inercial maximiza <σv>/ T 3/2 , que es ligeramente superior a la temperatura óptima para el confinamiento magnético.

Sistemas no térmicos

El análisis de Lawson se basa en la tasa de fusión y la pérdida de energía en un plasma termalizado. Existe una clase de máquinas de fusión que no utilizan plasmas termalizados, sino que aceleran directamente iones individuales hasta alcanzar las energías requeridas. Los ejemplos más conocidos son el migma , el fusor y el polywell .

Aplicado al fusor, el análisis de Lawson se utiliza como argumento para afirmar que las pérdidas por conducción y radiación son los principales impedimentos para alcanzar la potencia neta. Los fusores utilizan una caída de voltaje para acelerar y colisionar iones, lo que da como resultado la fusión. [ 9 ] La caída de voltaje se genera mediante jaulas de alambre, y estas jaulas conducen las partículas.

Los Polywells son mejoras de este diseño, diseñados para reducir las pérdidas por conducción eliminando las jaulas de alambre que las causan. [ 10 ] Sin embargo, se argumenta que la radiación sigue siendo un impedimento importante. [ 11 ]

Véase también

Notas

^a Es sencillo relajar estas suposiciones. La pregunta más difícil es cómo definirnorte{\displaystyle n}cuando los iones y los electrones difieren en densidad y temperatura. Considerando que se trata de un cálculo de producción y pérdida de energía por iones, y que cualquier concepto de confinamiento de plasma debe contener las fuerzas de presión del plasma, parece apropiado definir la densidad efectiva (electrónica).norte{\displaystyle n}a través de la presión (total)pag{\displaystyle p}comonorte=pag/2Ti{\displaystyle n=p/2T_{\mathrm {i} }}. El factor de2{\displaystyle 2}se incluye porquenorte{\displaystyle n}Por lo general se refiere a la densidad de los electrones solamente, peropag{\displaystyle p}aquí se refiere a la presión total. Dadas dos especies con densidades iónicasnorte1,2{\displaystyle n_{1,2}}números atómicosZ1,2{\displaystyle Z_{1,2}}temperatura iónicaTi{\displaystyle T_{\mathrm {i} }}y temperatura electrónicaTmi{\displaystyle T_{\mathrm {e} }}Es fácil demostrar que la potencia de fusión se maximiza mediante una mezcla de combustible dada pornorte1/norte2=(1+Z2Tmi/Ti)/(1+Z1Tmi/Ti){\displaystyle n_{1}/n_{2}=(1+Z_{2}T_{\mathrm {e} }/T_{\mathrm {i} })/(1+Z_{1}T_{\mathrm {e} }/T_{\mathrm {i} })}Los valores paranorteτ{\displaystyle n\tau },norteTτ{\displaystyle nT\tau }y la densidad de potencia debe multiplicarse por el factor(1+Z1Tmi/Ti)(1+Z2Tmi/Ti)/4{\displaystyle (1+Z_{1}T_{\mathrm {e} }/T_{\mathrm {i} })\cdot (1+Z_{2}T_{\mathrm {e} }/T_{\mathrm {i} })/4}. Por ejemplo, con protones y boro (Z=5{\displaystyle Z=5}) como combustible, otro factor de3{\displaystyle 3}debe incluirse en las fórmulas. Por otro lado, para los electrones fríos, todas las fórmulas deben dividirse por4{\displaystyle 4}(sin factor adicional paraZ>1{\displaystyle Z>1}).

Referencias

  1. Lawson, JD (diciembre de 1955). Algunos criterios para un reactor termonuclear útil (PDF) (Informe técnico). Centro de Investigación de Energía Atómica, Harwell, Berkshire, Reino Unido.
  2. 1 2 Lawson, JD (diciembre de 1955). "Algunos criterios para un reactor termonuclear productor de energía". Actas de la Sociedad Física, Sección B. 70 ( 1): 6– 10. Bibcode : 1957PPSB...70....6L . doi : 10.1088/0370-1301/70/1/303 .
  3. "Los científicos lograron la fusión nuclear autosostenible... pero ahora no pueden replicarla" . Sciencealert . 16 de agosto de 2022.
  4. Abu-Shawareb, H.; Acree, R.; Adams, P.; Adams, J.; Addis, B.; Aden, R.; Adrian, P.; Afeyan, BB; Aggleton, M.; Aghaian, L.; Aguirre, A.; Aikens, D.; Akre, J.; Albert, F.; Albrecht, M. (2022-08-08). "Se supera el criterio de Lawson para la ignición en un experimento de fusión inercial" . Physical Review Letters . 129 (7) 075001. Bibcode : 2022PhRvL.129g5001A . doi : 10.1103/PhysRevLett.129.075001 . hdl : 10044/1/99300 . ISSN 0031-9007 . PMID 36018710 .  
  5. 1 2 Spitzer, Lyman; Seeger, Raymond J. (noviembre de 1963). "Física de gases totalmente ionizados". American Journal of Physics . 31 (11): 890– 891. Bibcode : 1963AmJPh..31..890S . doi : 10.1119/1.1969155 . ISSN 0002-9505 . 
  6. "Convertidor de energía" . www.phys.ksu.edu . Universidad Estatal de Kansas . Consultado el 17 de febrero de 2023 .
  7. Wesson, J. (2004). "Tokamaks" . Oxford Engineering Science Series (48) (3.ª ed.). Oxford: Clarendon Press. 
  8. "Se ha marcado el triple producto de fusión más alto del mundo en plasmas de modo H de alto βp" . Archivado del original el 6 de enero de 2013.
  9. Hirsch, Robert L. (octubre de 1967). "Confinamiento electrostático inercial de gases de fusión ionizados". Journal of Applied Physics . 38 (11): 4522– 4534. Bibcode : 1967JAP....38.4522H . doi : 10.1063/1.1709162 . ISSN 0021-8979 . 
  10. Bussard, Robert W (2 de octubre de 2006). «El advenimiento de la fusión nuclear limpia: potencia y propulsión espacial de alto rendimiento». 57.º Congreso Internacional de Astronáutica . Reston, Virginia: Instituto Americano de Aeronáutica y Astronáutica. doi : 10.2514/6.iac-06-d2.8.05 . ISBN 978-1-62410-042-0.
  11. Rider, Todd H. (1997-04-01). "Limitaciones fundamentales en sistemas de fusión de plasma que no están en equilibrio termodinámico" . Física de plasmas . 4 (4): 1039– 1046. Bibcode : 1997PhPl....4.1039R . doi : 10.1063/1.872556 . hdl : 1721.1/11412 . ISSN 1070-664X . 
  • Derivación matemática , archivada en 2019 desde el original.