Neutronenfluss

Der Neutronenfluss (englisch Neutron flux), abweichend davon auch Neutronenflussdichte^[1][2][3][4] genannt,^[5] ist eine physikalische Größe der Reaktorphysik. Es ist eine skalare Größe. Anschaulich gibt er die Summe aller Wege an, die von den in einem Raumbereich vorhandenen freien Neutronen in einem Zeitintervall zurückgelegt werden, geteilt durch das Volumen des Raumbereichs und die Dauer des Zeitintervalls.^[6] Sein übliches Formelzeichen ist ${\displaystyle \Phi }$ (großes Phi), seine Dimension ${\displaystyle \mathrm {L} ^{-2}\mathrm {T} ^{-1}}$. Die übliche Maßeinheit ist cm⁻² s⁻¹, der Deutlichkeit wegen meist geschrieben als n cm⁻² s⁻¹ („Neutronen pro Quadratzentimeter und Sekunde“).

Der Neutronenfluss ist weder ein Fluss noch eine Flussdichte im Sinne der sonst üblichen physikalischen Nomenklatur.^[7][8]

Er kann durch die Anzahldichte ${\displaystyle n}$ der Neutronen und den Mittelwert ${\displaystyle v}$ der Geschwindigkeitsbeträge der Neutronen ausgedrückt werden:^{[6][9][10][11][7][12][13]}

${\displaystyle \Phi =n\cdot v}$.

Der Neutronenfluss lässt sich mittels Neutronendetektoren messen.

Freie Neutronen treten in Sternen, in Supernovae und – verursacht durch kosmische Strahlung oder Gewitter – in unserer natürlichen irdischen Umgebung auf. Von besonderer Bedeutung sind der Neutronenfluss und einige mit ihm verwandte Größen in Kernreaktoren, Abschirmungen usw.

Definition

Der Neutronenfluss in einem Raumbereich ist definiert als:

${\displaystyle \Phi (t,{\vec {r}}):={\frac {1}{V}}\sum _{i}v_{i}}$

Hier sind:

• ${\displaystyle i}$: die Zählvariable für die Neutronen (zu einer eventuellen Auswahl der gezählten Neutronen siehe "Neutronenwinkeldichte und weitere differentielle Größen")
• ${\displaystyle t}$: der Zeitpunkt, zu dem die Anzahl der Neutronen und deren Geschwindigkeit erfasst wird,
• ${\displaystyle {\vec {r}}}$: der Ortsvektor, der die Lage des Raumbereichs angibt,
• ${\displaystyle V}$: das Volumen des betrachteten Raumbereichs,
• ${\displaystyle v_{i}}$: der Betrag der Geschwindigkeit ${\displaystyle {\vec {v}}_{i}}$ des ${\displaystyle i}$-ten Neutrons.

Unabhängige Variable

Der Neutronenfluss hängt außer vom Ort auch von der kinetischen Energie ${\displaystyle E}$ der Neutronen ab und kann sich mit der Zeit ${\displaystyle t}$ ändern:

${\displaystyle \Phi =\Phi ({\vec {r}},E,t)}$.

Der Vektor ${\displaystyle {\vec {r}}}$ ist der Ortsvektor, der die Lage eines – im Grenzfall infinitesimal kleinen – Raumbereichs in einem dreidimensionalen Koordinatensystem angibt. Folglich hängt der Neutronenfluss im Allgemeinen von fünf unabhängigen Variablen ab.

Anschauliche Erklärung

Der Neutronenfluss im Raumpunkt ${\displaystyle {\vec {r}}}$ lässt sich auch veranschaulichen durch die Vorstellung einer kleinen Kreisscheibe mit dem Mittelpunkt ${\displaystyle {\vec {r}}}$, die von den Neutronen durchdrungen wird. Man stellt sich vor, die Kreisscheibe werde bei festem Mittelpunkt für jedes einzelne Neutron so im Raum gedreht, dass das Neutron in Richtung der Normalen einfällt. Der Neutronenfluss ist dann die Zahl der Neutronen, die in einem Zeitintervall die Scheibe durchdringen, geteilt durch das Zeitintervall und durch die Fläche der Scheibe. Das Drehen der Scheibe in alle Richtungen ergibt als Einhüllende eine Kugel mit dem Durchmesser der Scheibe.^[14][15]

Neutronenwinkeldichte und weitere differentielle Größen

Die Auswahl der zum Zeitpunkt ${\displaystyle t}$ zu zählenden Neutronen im Volumen ${\displaystyle V}$ um den Ort ${\displaystyle {\vec {r}}}$ kann danach erfolgen, ob die Neutronen eine Energie ${\displaystyle E_{i}}$ innerhalb eines Energiebereichs ${\displaystyle E-\delta E\leq E_{i}\leq E+\delta E}$ haben und zugleich ihre Flugrichtung ${\displaystyle {\hat {e}}_{i}}$ in der Umgebung einer vorgegebenen Flugrichtung mit dem Einheitsvektor ${\displaystyle {\vec {\Omega }}}$ liegt (${\displaystyle {\hat {e}}_{i}\cdot {\vec {\Omega }}\geq 1-\delta e}$). Daraus ergibt sich die differentielle Größe Neutronenwinkeldichte (englisch angular neutron density) oder kurz Winkeldichte ${\displaystyle n}$:

${\displaystyle n({\vec {r}},E,t,{\vec {\Omega }}):=\lim _{V\rightarrow 0}\lim _{\delta E\rightarrow 0}\lim _{\delta e\rightarrow 0}{\frac {1}{V}}{\frac {1}{\delta E}}{\frac {1}{4\pi (\delta e)^{2}}}\Phi (t,{\vec {r}},V,E,\delta E,{\vec {\Omega }},\delta e)}$.

Eine genaue Analyse der Transportprozesse von Neutronen in einem Kernreaktor (entsprechend auch von anderen Teilchen) erfordert die Definition noch weiterer, mit dem Neutronenfluss verknüpfter Größen. Die Namen dieser Größen sind nicht in allen reaktorphysikalischen Lehrbüchern gleich. Wir folgen hier der Namensgebung und der Definition einiger dieser Größen nach dem Lehrbuch Nuclear reactor theory von Bell und Glasstone.^[16]

Neutronenflussspektrum

Die Größe Neutronenflussspektrum, Neutronenspektrum,^[17] Neutronenenergiespektrum, energiediffenzieller Neutronenfluss oder energieabhängiger Neutronenfluss ist die partielle Ableitung des Neutronenflusses nach der Energie:

${\displaystyle \varphi ({\vec {r}},E,t)={\frac {\partial \Phi ({\vec {r}},E,t)}{\partial E}}}$.

Die Maßeinheit ist dementsprechend z. B. n cm⁻²s⁻¹eV⁻¹, Neutronen pro Quadratzentimer, Sekunde und Elektronenvolt. Der Fluss der Neutronen mit Energien zwischen ${\displaystyle E}$ und ${\displaystyle E+\partial E}$ ist ${\displaystyle \varphi (E)\,\partial E}$.

Flugrichtung der Neutronen

Um Neutronen nach ihrer Flugrichtung zu unterscheiden, sind zwei weitere unabhängige Variable erforderlich, die in dem Vektor Neutronenflugrichtung (Einheitsvektor der Neutronenflugrichtung) zusammengefasst werden.

Folgendes sei vorangestellt: Es gibt Größen in der Physik, die vom Ort und vom Impuls abhängen, zum Beispiel die Wellenfunktion der Quantenmechanik. Für solche Größen werden neben den Ortskoordinaten ${\displaystyle {\vec {r}}=(x,y,z)}$ auch die Impulskomponenten ${\displaystyle {\vec {p}}=(p_{x},p_{y},p_{z})}$ als unabhängige Variable in die Symbolik aufgenommen. Sei ${\displaystyle f}$ so eine Größe, die außerdem von der Zeit ${\displaystyle t}$ abhängt, dann symbolisiert man die Abhängigkeit der Größe ${\displaystyle f}$ von sieben unabhängigen Variablen mit ${\displaystyle f({\vec {r}},{\vec {p}},t)}$.

In der Reaktorphysik könnte man ebenfalls mit diesen sieben unabhängigen Variablen rechnen, den drei Ortskoordinaten, den drei Impulskomponenten und der Zeit. Aber anstelle der drei Impulskomponenten wählt man drei andere unabhängige Variable, die kinetische Energie ${\displaystyle E}$ der Neutronen und zwei Raumwinkelvariable ${\displaystyle {\vec {\Omega }}=(\theta ,\varphi )}$, die die Flugrichtung der Neutronen erfassen. Es sind der Polarwinkel ${\displaystyle \theta }$ und der Azimutwinkel ${\displaystyle \varphi }$ der (auf den Wert 1 normierten) Flugrichtung in Kugelkoordinaten, wie in der Abbildung gezeigt.

Der Vektor ${\displaystyle {\vec {\Omega }}}$, der Einheitsvektor in Richtung der Neutronenbewegung, wird definiert durch

${\displaystyle {\vec {\Omega }}={\frac {{\vec {v}}(E)}{|{\vec {v}}(E)|}}={\frac {{\vec {v}}(E)}{v(E)}}}$.

Dabei symbolisieren ${\displaystyle {\vec {v}}(E)}$ den Geschwindigkeitsvektor, ${\displaystyle |{\vec {v}}(E)|}$ und ${\displaystyle v(E)}$ den Betrag des Geschwindigkeitsvektors, jeweils in Abhängigkeit von der kinetischen Energie ${\displaystyle E}$.

Man beachte den Unterschied zwischen diesem zweikomponentigen Vektor ${\displaystyle {\vec {\Omega }}}$ und dem gewöhnlichen Raumwinkel, der ein Skalar ist und mit ${\displaystyle \Omega }$ symbolisiert wird, was zu Fehlinterpretationen führen kann.

Mit diesen unabhängigen Variablen symbolisiert man eine reaktorphysikalische Größe ${\displaystyle f}$ dann entsprechend mit ${\displaystyle f({\vec {r}},E,{\vec {\Omega }},t)}$. Diese Wahl der unabhängigen Variablen hat sich als zweckmäßig erwiesen und ist aus physikalischer Sicht einer Wahl der Impulskomponenten gleichwertig. Aus der kinetischen Energie und dem Raumwinkelvektor können die Impulskomponenten berechnet werden und umgekehrt.

Neutronenwinkelfluss

Der Neutronenwinkelfluss (Angular neutron flux) oder kurz Winkelfluss ist die abhängige Variable der Neutronentransportgleichung und damit eine der wichtigsten physikalischen Größen der Reaktortheorie überhaupt. Sie wird definiert als

${\displaystyle \psi ({\vec {r}},E,{\vec {\Omega }},t)=v\cdot n({\vec {r}},E,{\vec {\Omega }},t)}$.

Neutronenfluss oder Neutronenskalarfluss

Aus dem Neutronenwinkelfluss ergibt sich der oben definierte Neutronenfluss einfach als Integral des Winkelflusses über alle Raumrichtungen:

${\displaystyle \Phi ({\vec {r}},E,t)=\int _{4\pi }d{\vec {\Omega }}\psi ({\vec {r}},E,{\vec {\Omega }},t)}$.

Weil der Neutronenfluss durch Integration über die vektorielle Größe Neutronenwinkelfluss gebildet wird, bezeichnet man den Neutronenfluss im Kontext der Neutronentransporttheorie auch als Neutronenskalarfluss oder skalaren Neutronenfluss.

Neutronenfluenz

Durch Integration des Neutronenflusses oder des Neutronenflussspektrums über die Zeit, z. B. die Dauer einer Bestrahlung, ergibt sich entsprechend die totale bzw. energieabhängige Neutronenfluenz. Sie ist wichtig zur Berechnung z. B. der durch Neutronen verursachten Strahlenschäden oder der Ausbeute einer Neutronenaktivierung.

Bedeutung für Kernreaktoren

Die Größe Neutronenfluss dient im Zusammenhang mit Kernreaktoren hauptsächlich dazu, die Frage zu beantworten, ob und warum eine Spaltstoffanordnung (Kernreaktor) kritisch wird, und um Kernreaktionsraten im Reaktor zu berechnen. Ohne Kenntnis räumlich und energetisch hinreichend „aufgelöster“ Kernreaktionsraten können beide Fragen nicht beantwortet werden. Die entscheidende Größe Kernreaktionsratendichte ist das Produkt aus Neutronenfluss und der materialabhängigen Größe makroskopischer Wirkungsquerschnitt. Der makroskopische Wirkungsquerschnitt hängt von Teilchendichten und nuklearen Eigenschaften aller Nuklide ab, mit denen sich die Neutronen den Raumbereich „teilen“.

In einem Kernreaktor ist der Neutronenfluss im Allgemeinen abhängig von Ort, Neutronenenergie und Zeit. In einem gleichbleibend kritisch gehaltenen Reaktor kann man die Zeitabhängigkeit über eine kleine Zeitspanne vernachlässigen. Man spricht dann vom stationären Zustand des Reaktors. Der totale oder der thermische Neutronenfluss ist relativ einfach zu messen, z. B. mit in den Reaktorkern eingebauten Spaltkammern.

Aus dem Produkt von Neutronenfluss und einem speziellen Wirkungsquerschnitt, einer material- und ortsabhängigen Größe, ergibt sich die Leistungsdichte in einem Volumenbereich des Reaktors. Nimmt man (etwas vereinfachend) an, dass die Form der räumlichen Verteilung des Neutronenflusses über den Reaktorkern immer gleich bleibt, dann genügt die Messung des Neutronenflusses an einer Stelle, um auf die Gesamtleistung des Reaktors zu schließen. Eine Temperaturmessung wäre dagegen als Leistungsmaß nicht geeignet, denn die Temperatur an einer Stelle ist Resultat der Leistung an der Stelle selbst, der Leistung in benachbarten Bereichen in der jüngeren und mittleren Vergangenheit sowie der Kühlleistung.

Das Fluss-Messsignal wird daher allgemein zur Steuerung und Überwachung des Reaktors verwendet. Auch bei abgeschaltetem (unterkritischem) Reaktor wird die Flussmessung ständig in Betrieb gehalten. Eine zu diesem Zweck eingebaute radioaktive Neutronenquelle sorgt stets für einen geringen Neutronenfluss; dadurch wird die Funktion der Messinstrumentierung dauernd überwacht.

Mittelwert und Diskretisierung

Ein gemessener Neutronenflusswert ist stets ein Mittelwert über einen gewissen Raumbereich, ein Energieintervall und eine Zeitdauer. Auch für Berechnungen des Neutronenflusses ist eine Diskretisierung der unabhängigen Variablen notwendig, da die entsprechenden Gleichungen nur numerisch lösbar sind. Die Intervallgrößen für die Diskretisierung werden in der Regel vor dieser Berechnung festgelegt.

Räumliche Diskretisierung

Größe und Form des interessierenden Raumbereichs können sehr unterschiedlich sein. Bei einer ausgedehnten Anordnung, etwa einem Reaktorkern, setzt die verfügbare Rechenkapazität der Auflösung in kleine Raumbereiche praktische Grenzen. Wie fein die Ortsauflösung gewählt werden muss, hängt zusammen mit der mittleren freien Weglänge der Neutronen; diese wird durch die Neutronenenergie und das jeweilige Medium bestimmt. Typische Abmessungen der einzelnen „Zelle“ bei Reaktorberechnungen liegen im Zentimeterbereich.^[20]

Energiediskretisierung

Die Energieauflösung wird in Reaktorberechnungen je nach Fragestellung gewählt, vom Gesamtbereich der möglichen Neutronenenergie (10⁻⁴ bis 2·10⁷) eV – also gar keiner „Auflösung“ – bis hin zu einigen hundert Energieintervallen („Neutronengruppen“). Bei thermischen Reaktoren reicht für manche Fragestellungen die Unterteilung in zwei Gruppen aus (thermischer Fluss und schneller Fluss; siehe erste Abbildung). Für schnelle Reaktoren sind Berechnungen z. B. oft mit 26 Gruppen durchgeführt worden.^[21]

Durch Integration des Neutronenflussspektrums über das Energieintervall der jeweiligen Neutronengruppe werden die sogenannten Gruppenflüsse berechnet. Dies sind zum Beispiel im Fall von zwei Energiegruppen die Integrale

${\displaystyle \Phi _{\text{thermisch}}=\int _{E_{u}}^{E_{1}}\varphi ({\vec {r}},E,t)\ \mathrm {d} E}$ (thermischer Fluss) und

${\displaystyle \Phi _{\text{schnell}}=\int _{E_{1}}^{E_{o}}\varphi ({\vec {r}},E,t)\ \mathrm {d} E}$ (schneller Fluss).

Im Beispiel der ersten Abbildung wurden die Energiegrenzen ${\displaystyle E_{u}=10^{-4}\ \mathrm {eV} }$, ${\displaystyle E_{1}=0{,}625\ \mathrm {eV} }$ und ${\displaystyle E_{o}=20\ \mathrm {MeV} }$ verwendet.^[22]

Das Neutronenflussspektrum wird u. a. als Wichtungsfunktion benötigt, um die sog. Gruppenkonstanten für Neutronendiffusionsberechnungen zu erhalten.

Die Abbildung zeigt das Neutronenflussspektrum eines frischen Brennelements eines Druckwasserreaktors (DWR). Es handelt sich um das gleiche Modell, zu dem auch die am Artikelanfang abgebildeten Neutronenflussverteilungen gehören, und wurde ebenfalls mit dem Programmsystem HELIOS 1.8^[23][24] berechnet. Dargestellt ist das Neutronenflussspektrum einer ausgewählten Kühlmittelregion (Region 1, die Region an der linken oberen Ecke) und der Mittelwert des Neutronenflussspektrums, gemittelt über das gesamte Brennelement. Man erkennt:

• Die Energie beim Maximum des Neutronenflussspektrums liegt nahe an der thermischen Energie von 0,0253 eV (20 °C).
• Das Neutronenflussspektrum bricht bei etwa 20 MeV ab, da durch Kernspaltung keine Neutronen mit einer höheren Energie entstehen.

Beispiel: Neutronenfluss in einem Druckwasserreaktor

Die erste Abbildung des Artikels zeigt die horizontale Ortsabhängigkeit des thermischen und des schnellen Neutronenflusses über den Querschnitt eines Brennelements eines Druckwasserreaktors in einem stationären Zustand, gewonnen aus einer Neutronentransportrechnung^[24] mit dem Zell- und Abbrandprogramm HELIOS 1.8^[23] und Wirkungsquerschnitten aus der Kerndatenbibliothek ENDF/B-VI data files (Rose and Dunford, 1990).^[25]

Das Gitter des Brennelements enthält zwei Typen von „Zellen“: Brennstabzellen und Führungsrohrzellen. Die Brennstabzelle enthält drei Materialien: Brennstoff (rot), Zirkalloy-Hülle (grün) und Wasser (blau), die Führungsrohrzelle nur das Führungsrohr und Wasser. Die nebenstehende Abbildung zeigt links die beiden Zelltypen und rechts daneben ihre Diskretisierung in 6 bzw. 8 Regionen.

Das Brennelement ist eine 18×18-Anordnung mit 300 Brennstäben und 24 Regelstabführungsrohren, in die Neutronenabsorber eintauchen können. In dieser Rechnung sind die Regelstäbe gezogen. Da das Brennelement symmetrisch ist, reichte es aus, ein Viertel des Brennelements zu berechnen, das die Grafik zeigt. Das Symmetriezentrum (die Mitte des Brennelement-Querschnitts) ist auf dem Diskretisierungsgitter links oben als kleines schwarzes Quadrat markiert.

An den Orten der gezogenen Regelstäbe ist der thermische Fluss stark erhöht, im mittleren Teil der ersten Abbildung des Artikels deutlich an den sechs rot bis orange eingefärbten Maxima zu erkennen. Der schnelle Fluss dagegen ist an diesen Orten vergleichsweise niedrig. Die Zahlen innerhalb der Abbildung geben die Flusswerte in n cm⁻²s⁻¹ bei einer spezifischen Wärmeleistung von 37,4 W/g Schwermetall an.^[26] In diesem Modellfall liegt der schnelle Fluss im Intervall (2,2·10¹⁴ – 2,3·10¹⁴) Neutronen cm⁻²s⁻¹, der thermische Fluss im Intervall (2,9·10¹³ – 4,4·10¹³) Neutronen cm⁻²s⁻¹.

Kernreaktoren

Die Neutronenflusswerte (genauer: Neutronenwinkelflusswerte) in einem Kernreaktor gehorchen der Neutronentransportgleichung,^[6][10] einer Bilanzgleichung für Neutronen. Diese Integro-Differentialgleichung numerisch zu lösen gehört zu den anspruchsvollsten Aufgaben der Physik und numerischen Mathematik überhaupt.

Programmsysteme zur numerischen Lösung der Neutronentransportgleichung, zum Beispiel für ein Kernreaktor-Brennelement in einer räumlichen und energetischen Auflösung, wie sie in der Abbildung am Anfang des Artikels dargestellt ist, werden von nur sehr wenigen spezialisierten Firmen auf der Welt entwickelt. Es sind entweder staatlich dominierte Firmen, zum Beispiel in Frankreich unter dem Dach der EDF oder dem CEA (Programmsystem APOLLO), oder privatrechtliche Firmen wie Studsvik Scandpower (Programmsysteme CASMO und HELIOS^[23]). Die privatrechtlichen Firmen sind meist durch Outsourcing von Arbeitsgruppen entstanden, die mit der Programmentwicklung in einem staatlichen Institut oder an einer Universität begonnen haben.^[27] Die Entwicklung solcher Programmsysteme erfordert Dutzende von Personenjahren und kann an einer Universität gegenwärtig nicht geleistet werden. Hinzu kommt für jeden Entwickler eine Einarbeitungszeit in die physikalische Theorie und die mathematisch-numerischen Lösungsverfahren von mehreren Jahren, ehe er mit dem Programmieren überhaupt beginnen kann. Auch ein Anwender eines solchen Programms sollte für die Einarbeitung etwa drei bis fünf Jahre veranschlagen.

Für die Entwicklung neuer Kernreaktortypen, zum Beispiel eines Salzschmelzenreaktors für den Leistungsbetrieb oder für den Routinebetrieb eines existierenden Kernkraftwerks, ist eine ganze Kette von Programmen erforderlich, u. a. auch weit weniger „aufwendige“ Neutronentransportprogramme, die Spezialaufgaben numerisch lösen. Das erste Glied dieser Kette bleibt aber stets ein hochleistungsfähiges Neutronentransportprogramm. Der gängige Typname Zell- und Abbrandprogramm für ein solches Programm umschreibt dessen Leistungsumfang nur andeutungsweise.^[28] Diese Zell- und Abbrandprogramme wurden (Stand Jahr 2007) ausschließlich in Fortran unterschiedlicher Versionen geschrieben. Ein Programmsystem dieses Typs kostete im Jahr 2000 mindestens 100000 Dollar.

Das Verhalten der Neutronen im Reaktor kann jedoch angenähert auch als Diffusionsvorgang beschrieben werden.^[6][10] Ein solches Neutronendiffusionsprogramm berechnet ausschließlich den Neutronenfluss, nicht den Neutronenwinkelfluss. Ein Programm dieses Typs kann, grob überschlagen, in zwei bis fünf Personenjahren (je nach räumlicher Dimension: 1D, 2D, 3D) entwickelt werden. Die übliche Programmiersprache für ein Neutronendiffusionsprogramm ist ebenfalls Fortran. Programme dieses Typs werden für öffentliche Institute und Universitäten kostenlos von der Datenbank der NEA Computer Program Library, speziell unter der Kategorie C. STATIC DESIGN STUDIES bereitgestellt (Einarbeitungszeit einige Monate).

Für spezielle Aufgaben werden auch Monte-Carlo-Rechenprogramme wie MCNP („Monte-Carlo N-Particle Transport Code“) eingesetzt.^[29][30]

Fusionsreaktoren

Zur Berechnung der Neutronenflussverteilung und damit zusammenhängender Größen in Fusionsreaktoren und entsprechenden Versuchsanlagen wird allgemein die Monte-Carlo-Methode eingesetzt.^[31][32][33]

Natürlicher Neutronenfluss

Neutroneninduzierte Kernreaktionen in AGB-Sternen (englisch Asymptotic giant branch) sind für die meisten natürlichen Elemente verantwortlich, die durch Nukleosynthese entstanden sind. Das sind Elemente massereicher als Eisen. Der Neutronenfluss ist relativ niedrig und liegt in der Größenordnung von 10⁵ bis 10¹¹ n cm⁻²s⁻¹. Er führt zur Nukleosynthese durch den s-Prozess (englisch Slow-Neutron-Capture-Prozess).

Im Gegensatz dazu ist der Neutronenfluss nach der Explosion eines massereichen Sterns (Supernova) sehr hoch und erreicht die Größenordnung von 10³² n cm⁻²s⁻¹.^[34] Das führt zur Nukleosynthese durch den r-Prozess (englisch Rapid-Neutron-Capture-Prozess).

In der Erdatmosphäre erzeugt die kosmische Strahlung, die vorwiegend aus hochenergetischen Protonen (1 GeV und höher) besteht, durch Spallation der Atomkerne von Stickstoff und Sauerstoff freie Neutronen. Zugleich erzeugt sie sekundäre Protonen und geladene und neutrale Pionen, die ihrerseits in weiteren Reaktionen Neutronen freisetzen können. Der so entstehende Neutronenfluss hängt stark vom primären Protonenstrom und dem Ort der Reaktion in der Atmosphäre ab.

Bei einem Blitz werden Atome des atmosphärischen Stickstoffs und Sauerstoffs ionisiert, aber auch das in einer Wasserdampfwolke immer vorhandene schwere Wasser. Im elektrischen Feld der Blitzentladung können die Ionen (wie die freigesetzten Elektronen auch) beschleunigt werden und Kernreaktionen auslösen, die ihrerseits Neutronen freisetzen. Insbesondere kann es auch zu Fusionsreaktionen kommen.^[35] Für diese Vorgänge stehen Theorie und Experiment noch am Anfang. Der atmosphärische Neutronenfluss während eines Gewitters erreicht bis zu etwa 90 n cm⁻²s⁻¹.^[36][37][38]

Weitere Zahlenbeispiele für den Neutronenfluss

Vom „Komitee Forschung mit Neutronen“^[39] werden Forschungsneutronenquellen (FNQ) unterschieden in

• Spallationsquellen. Es gibt 2 Spallations-Neutronenquellen in Europa.
• Reaktoren mit hohem Neutronenfluss (${\displaystyle \Phi _{\mathrm {th} }}$ > 10¹⁵ n cm⁻²s⁻¹). Es gibt einen solchen Hochflussreaktor in Europa, und zwar im ILL.
• Reaktoren mit mittlerem Neutronenfluss (10¹⁴ n cm⁻²s⁻¹ < ${\displaystyle \Phi _{\mathrm {th} }}$ < 10¹⁵ n cm⁻²s⁻¹). Es gibt 3 solche FNQ in Europa.
• Reaktoren mit niedrigem Neutronenfluss (${\displaystyle \Phi _{\mathrm {th} }}$ ≤ 10¹⁴ n cm⁻²s⁻¹). Es gibt 7 solche FNQ in Europa.

${\displaystyle \Phi _{\mathrm {th} }}$ meint den thermischen Neutronenfluss.

Die Spallations-Neutronenquelle SINQ erreicht einen totalen Neutronenfluss von 10¹⁴ n cm⁻²s⁻¹, der Forschungsreaktor FRM-II 8×10¹⁴ n cm⁻²s⁻¹.^[40]

In einem Fusionsreaktor wird die Wand des Plasmagefäßes ebenfalls einem Neutronenfluss von etwa 10¹⁴ n cm⁻²s⁻¹ ausgesetzt sein,^[41] der hier ganz überwiegend aus Neutronen der hohen Energie von etwa 14 MeV besteht.

Weblinks

• NEA Computer Program Library
• C. STATIC DESIGN STUDIES

Einzelnachweise und Anmerkungen

1. ↑ Dieter Emendörfer, Karl-Heinz Höcker: Theorie der Kernreaktoren. Bibliographisches Institut, Mannheim / Wien / Zürich 1970 (380 S.). 
2. ↑ Dieter Smidt: Reaktortechnik. 2. Auflage. Braun, Karlsruhe 1976, ISBN 3-7650-2019-2 (XVI, 325 S.).  Smidt gibt auf S. 20 für Neutronenfluss fast wörtlich dieselbe Definition wie Emendörfer/Höcker, Theorie der Kernreaktoren Bd. 1 auf S. 63 für Neutronenflussdichte, nämlich: „der Weg, der von allen Neutronen einer Volumeinheit in einer Zeiteinheit insgesamt zurückgelegt wird.“
3. ↑ Albert Ziegler: Lehrbuch der Reaktortechnik. Springer, Berlin / Heidelberg 1983, ISBN 3-540-12198-6 (XI, 242 S.). 
4. ↑ Albert Ziegler, Hans-Josef Allelein (Hrsg.): Reaktortechnik: Physikalisch-technische Grundlagen. 2., neu bearbeitete Auflage. Springer Vieweg, Berlin 2013, ISBN 978-3-642-33846-5 (634 S., eingeschränkte Vorschau in der Google-Buchsuche [abgerufen am 21. Januar 2018]). 
5. ↑ Die Einordnung der Größe Neutronenfluss als Dichte hat sich vor allem in der deutschsprachigen Fachliteratur eingebürgert. Die Größe Neutronenfluss ist per definitionem eine intensive Größe, so dass sich ein Appendix -dichte erübrigt.
6. ↑ ^{a b c d} Samuel Glasstone, Milton C. Edlund: The elements of nuclear reactor theory. MacMillan, London 1952 (VII, 416 S.).  Diese Monografie nimmt eine herausragende Stellung ein, weil sie wie keine andere die damals junge Generation der Reaktorphysiker in West und Ost und die späteren Lehrbuchschreiber geprägt hat. Sie ist im 6. Druck vom Februar 1957 vollständig online einsehbar.babel.hathitrust.org. Volltextsuche ist möglich.
7. ↑ ^{a b} James J. Duderstadt, Louis J. Hamilton: Nuclear reactor analysis. Wiley, New York 1976, ISBN 978-0-471-22363-4 (xvii, 650 S.).  Die Autoren schreiben auf S. 106: „… the tradition in nuclear engineering of referring to this quantity as the neutron "flux" is very misleading.“ ( … die Tradition in der Kerntechnik, diese Größe als Neutronen"fluss" zu bezeichnen, ist sehr irreführend.)
8. ↑ Was die hier zitierten Monografien betrifft, sind Duderstadt und Hamilton (S. 106) die einzigen Autoren, die den Fehlgriff des Namens Neutronenfluss explizit ansprechen. Implizit tun das auch andere englischsprachige Autoren, z. B. Weinberg und Wigner s. o., indem sie den Namen Track length (Bahnlänge) ins Spiel bringen. Auch dieser Name entfernt sich weit von der üblichen Namensgebung für physikalische Größen. Er suggeriert, diese Größe habe die Dimension einer Länge. Daran ändert auch die Tatsache nichts, dass von einigen Autoren korrekt angehängt wird, dass diese Bahnlänge noch durch das Volumen des Raumbereichs und das Zeitintervall zu teilen sei, in dem die Bahnlänge „gemessen“ worden ist. Die Autoren deutschsprachiger Monografien (zumindest die hier zitierten) definieren und interpretieren die Namen Neutronenfluss und Bahnlänge unkritisch. Manche Autoren „verwässern“ ihre Definitionen sogar noch, zum Beispiel durch Maßeinheiten: „Diese (die Neutronenflussdichte) kann aufgefasst werden als die gesamte, von allen Neutronen je cm³ und s zurückgelegte Bahnlänge.“(Ziegler 1983, S. 58)
9. ↑ Alvin M. Weinberg, Eugene Paul Wigner: The physical theory of neutron chain reactors. Univ. of Chicago Press, Chicago 1958, ISBN 0-226-88517-8 (XII, 800 S.).  Die Autoren schreiben auf S. 23: „The quantity nv (i.e., the flux of incident particles) has the dimension of cm⁻² sec⁻¹. It is sometimes also called "track length,“ since it is the total distance traveled during unit time by all particles contained in unit volume." (Die Größe nv (d. h. der Fluss der einfallenden Teilchen) hat die Dimension cm⁻² sec⁻¹. Sie wird manchmal auch „Bahnlänge“ genannt, da es sich um die Gesamtdistanz handelt, die während der Zeiteinheit von allen in der Volumeneinheit enthaltenen Teilchen zurückgelegt wird.)
10. ↑ ^{a b c} Karl Wirtz, Karl Heinz Beckurts: Elementare Neutronenphysik. Springer, Berlin 1958 (VIII, 243 S., eingeschränkte Vorschau in der Google-Buchsuche [abgerufen am 3. Januar 2018]). 
11. ↑ Aleksej D. Galanin: Theorie der thermischen Kernreaktoren. Teubner, Leipzig 1959 (XII, 382 S.).  Die Monografie ist im gleichen Jahr original in russischer Sprache erschienen und ein Jahr später bei Pergamon Press in englischer Sprache unter dem Titel „Thermal reactor theory“.
12. ↑ Rudi J. J. Stamm'ler, Máximo J. Abbate: Methods of steady-state reactor physics in nuclear design. Acad. Press, London 1983, ISBN 0-12-663320-7 (XVI, 506 S.). 
13. ↑ Paul Reuss: Neutron physics. EDP Sciences, Les Ulis 2008, ISBN 978-2-7598-0041-4 (xxvi, 669, eingeschränkte Vorschau in der Google-Buchsuche).  In dieser Monografie werden auf S. 98 die Größen Neutronendichte, Neutronenfluss und Kernreaktionsrate sehr klar definiert.
14. ↑ K. H. Beckurts, K. Wirtz: Neutron Physics. Springer 1964, ISBN 978-3-642-87616-5, Seite 82–83
15. ↑ A. Ziegler, H.-J. Allelein (Hrsg.): Reaktortechnik. Physikalisch-technische Grundlagen. 2. Auflage. Springer 2013, ISBN 978-3-642-33845-8, Seite 58
16. ↑ George I. Bell, Samuel Glasstone: Nuclear reactor theory. Van Nostrand Reinhold, New York 1970, S. 2 ff. (XVIII, 619 S.). 
17. ↑ "Neutronenspektrum" ist doppeldeutig. Es könnte sich auch um das Neutronenanzahldichtespektrum handeln.
18. ↑ Der Schnitt kann parallel zu allen Außenflächen des Quaders geführt werden. Das Maximum des Neutronenflusses liegt im Mittelpunkt des Quaders. Wird der Schnitt nicht durch den Mittelpunkt geführt, verringert sich der Neutronenfluss gegenüber dem durch den Mittelpunkt geführten nur um eine Konstante, die mit dem Abstand des Schnitts vom Mittelpunkt kleiner wird.
19. ↑ Da die Neutronendiffusionsgleichung eine homogene lineare Differentialgleichung ist, ist mit einer Neutronenfluss-Lösung ${\displaystyle \Phi }$ auch der mit einer beliebigen Konstanten multiplizierte Neutronenfluss eine Lösung der Neutronendiffusionsgleichung.
20. ↑ Im Fall des in der Abbildung dargestellten Rechenmodells wurde der „Zentimeterbereich“ teilweise unterschritten.
21. ↑ H. Giese: KfK Analysis of the Superphenix-1 control rod experiments part 2: Rod worth calculations. Kernforschungszentrum Karlsruhe Report KfK-4896 (1992)
22. ↑ Bei Reaktorberechnungen mit nur zwei Gruppen ist die Grenzenergie 0,625 eV ein „Quasistandard“.
23. ↑ ^{a b c} Rudi J. J. Stamm'ler et al.: HELIOS Methods: Version 1.8. Studsvik Scandpower 2003 (192 S.). 
24. ↑ ^{a b} RK: PWR-Calculations with the Code-System HELIOS 1.8, Studsvik 2005 International User’s Group Meeting, Charlotte, NC, USA, June 1-3, 2005.
25. ↑ www-nds.iaea.org
26. ↑ Der Neutronenfluss ist in guter Näherung proportional der spezifischen Wärmeleistung. Die spezifische Wärmeleistung des Reaktors ist der Quotient aus der Gesamt-Wärmeleistung und der anfangs eingesetzten Schwermetallmasse des Reaktorkerns. Der Reaktor des Kernkraftwerks Emsland zum Beispiel wird bei einer thermischen Leistung von 3850 MW betrieben. Die anfangs eingesetzte Schwermetallmasse beträgt 103 t. Daraus ergibt sich eine mittlere spezifische Leistung von 37,4 W/g.
27. ↑ Studsvikstudsvik.episerverhosting.com im Fall des Programmsystems CASMO, Uni Oslo im Fall des Programmsystems HELIOS. Beide Firmen fusionierten Anfang der 2000er Jahre zur Firma Studsvik ScandPower
28. ↑ Bei der Entwicklung eines neuen Reaktortyps benötigt der Anteil dessen, was man unter Neutronenphysik zusammenfassen kann, nur ca. 20 % des Gesamtaufwands.
29. ↑ R. A. Forster, L. J. Cox, R. F. Barrett et al.: MCNP Version 5. Nuclear Instruments in Physics Research Section B Band 213 (2004) Seite 82–86
30. ↑ Z. Xu, J. Rhodes, K. Smith: CASMO-5 versus MCNP-5 benchmark of radial power profile in a fuel pin. Int. Conf. on Mathematics, Computational Methods and Reactor Physics, Saratoga Springs, 2009 (PDF; 371 kB)
31. ↑ Y. Li, L. Lu, A. Ding, H. Hu, Q. Zheng, S. Zheng, Y. Wu: Benchmarking of MCAM 4.0 with the ITER 3D model. Fusion Engineering and Design Band 82 (2007) Seite 2861–2866
32. ↑ S. P. Simakov, U. Fischer, K. Kondo and P.Pereslavtsev: Status of the McDeLicious Approach for the D-Li Neutron Source Term Modeling in IFMIF Neutronics Calculations. Fusion Science and Technology Band 62 (2012) Seite 233–239
33. ↑ P. Pereslavtsev, L. Lu, U. Fischer, O. Bitz: Neutronic analyses of the HCPB DEMO reactor using a consistent integral approach. Fusion Engineering and Design Band 89 (2014) Seite 1979–1983
34. ↑ E. Margaret Burbidge, G.R. Burbidge, William Fowler, Fred Hoyle: Synthesis of the Elements in Stars. In: Rev. Mod. Phys. Band 29, Nr. 4, 1957, S. 548–650, doi:10.1103/RevModPhys.29.547. 
35. ↑ Neutrons Born In Lightning. PhysOrg, 2005
36. ↑ Christoph Köhn, Ute Ebert: Calculation of beams of positrons, neutrons and protons associated with terrestrial gamma-ray flashes. In: Journal of Geophysical Research: Atmospheres, 2015, 120, S. 1620–1635. doi:10.1002/2014JD022229
37. ↑ A. V. Gurevich, A. M. Almenova: Observations of high-energy radiation during thunderstorms at Tien-Shan. In: Physical Review D. Americal Physical Society, 2016, 94 (2), S. 023003. doi:10.1103/PhysRevD.94.023003
38. ↑ Christoph Köhn, Gabriel Diniz, Mushin N. Harakeh: Production mechanisms of leptons, photons, and hadrons and their possible feedback close to lightning leaders. In: Journal of Geophysical Research: Atmospheres, 2017, 122, S. 1365–1383, doi:10.1002/2016JD025445
39. ↑ sni-portal.de
40. ↑ Forschungneutronenquelle Heinz Maier-Leibnitz
41. ↑ Weston M. Stacey: Fusion. An Introduction to the Physics and Technology of Magnetic Confinement Fusion. Wiley-VCH, 2010, ISBN 978-3-527-40967-9, eingeschränkte Vorschau in der Google-Buchsuche