@@ -147,7 +147,7 @@ Target nuclei (incident projectile, resulting fragments) resulting nuclei
---
-## Transmutation reactions
+# Transmutation reactions
Reactions that result in a change of the isotope
@@ -183,7 +183,7 @@ Pb208(n,2n)Pb207
---
-## Transmutation of lead to gold
+# Transmutation of lead to gold
@@ -212,7 +212,7 @@ Pb208(n,2n)Pb207
---
-## Q values
+# Q values
Amount of energy absorbed (-ve) or release (+ve) during the nuclear reaction
@@ -229,7 +229,7 @@ Online Q value calculator at [NNDC](https://www.nndc.bnl.gov/qcalc/)
---
-## Fusion fuels
+# Fusion fuels
@@ -325,7 +325,7 @@ Q values of fusion fuel reactions
---
-## Reaction rate equation
+# Reaction rate equation
- The reaction rate ($RR$) can be found by knowing the number of neutrons per unit volume ($n$), the velocity of neutrons ($v$), the material density ($p$), Avogadro's number ($N_{a}$), the microscopic cross section at the neutron energy ($\sigma_{e}$) and the atomic weight of the material ($M$).
@@ -337,7 +337,7 @@ $$ RR = \frac{nv\rho N_{a}\sigma_{e} }{M} = \phi N_{d} \sigma_{e} = \phi \Sigma_
---
-## Macroscopic cross section
+# Macroscopic cross section
@@ -419,7 +419,7 @@ Reactions have characteristics
---
-## Energy distribution
+# Energy distribution
@@ -460,7 +460,7 @@ Source [IAEA nuclear data services](https://nds.iaea.org/dataexplorer/?target_el
---
-## Nuclear data libraries
+# Nuclear data libraries
There are several groups that produce and distribute nuclear data
@@ -538,7 +538,7 @@ $\xi = 1+ \frac{(A-1)^2}{2A} ln \frac{(A-1)}{(A+1)}$
---
-## Why lithium
+# Why lithium
@@ -581,7 +581,7 @@ $\xi = 1+ \frac{(A-1)^2}{2A} ln \frac{(A-1)}{(A+1)}$
---
-## Other materials
+# Other materials
## Tungsten
@@ -601,7 +601,7 @@ $\xi = 1+ \frac{(A-1)^2}{2A} ln \frac{(A-1)}{(A+1)}$
---
-## Neutron spectra through materials
+# Neutron spectra through materials
By knowing the materials present can you identify which blanket results in which spectrum
diff --git a/3_prompt_response.md b/3_prompt_response.md
index 380b486..62793a1 100644
--- a/3_prompt_response.md
+++ b/3_prompt_response.md
@@ -176,7 +176,7 @@ TODO helium production plot in Iron -->
-# Tritium Breeding Ratio
+# [Tritium Breeding Ratio](https://prezi.com/embed/rnzt6pjj-xfu/?bgcolor=ffffff&lock_to_path=0&autoplay=1&autohide_ctrls=1&landing_data=bHVZZmNaNDBIWnNjdEVENDRhZDFNZGNIUE43MHdLNWpsdFJLb2ZHanI0eWk1QlBaUER3dVArS1hRQTAxNXdDZWNRPT0&landing_sign=ABm-Z3JCWCuKHnLF1Q-0yjuTsqyWAQdv3CEpUjcYcXk)
diff --git a/all_slides.md b/all_slides.md
new file mode 100644
index 0000000..c9f0e1b
--- /dev/null
+++ b/all_slides.md
@@ -0,0 +1,1109 @@
+---
+marp: true
+# theme: gaia
+# theme: uncover
+# _class: lead
+paginate: true
+backgroundColor: #fff
+title: Neutronics Analysis of Fusion Systems
+description: Presentation slides for the fusion energy neutronics workshop
+author: Jonathan Shimwell
+keywords: fusion,neutronics,neutron,photon,radiation,simulation,openmc,dagmc
+style: |
+ .columns {
+ display: grid;
+ grid-template-columns: repeat(2, minmax(0, 1fr));
+ gap: 1rem;
+ },
+ .columns3 {
+ display: grid;
+ grid-template-columns: repeat(3, minmax(0, 1fr));
+ gap: 1rem;
+ },
+ h1 {
+ text-align: center
+ }
+---
+
+
+
+# Neutronics Analysis of Fusion Systems
+
+![Neutron](images/cover.png)
+
+
+Slides and scripts to reproduce plots with simulations are available on GitHub
+https://github.com/fusion-energy/fusion-neutronics-presentation-slides
+
+---
+
+
+# Why is neutronics useful
+
+
+![bg vertical height:15cm left:10%](images/why_neutronics.png)
+- **Radioactivity** - Neutrons activate material, making it radioactive leading to handling and waste storage requirements.
+- **Hazardous** - Neutrons are Hazardous to health and shielded will be needed to protect the workforce.
+- ***Produce fuel*** - Neutrons will be needed to convert lithium into tritium to fuel the reactor.
+- ***Electricity*** - 80% of the energy release by each DT reaction is transferred to the neutron.
+- ***Structural integrity*** - Neutrons cause damage to materials such as embrittlement, swelling, change conductivity …
+- ***Diagnose*** - Neutrons are an important method of measuring a variety of plasma parameters (e.g. Q value).
+
+---
+
+# Topics Covered
+
+- Nuclear data
+- Prompt responses
+- Delayed responses
+- Simulation approaches
+
+---
+
+
+
+
+
+# Nuclear data
+
+
+
+- 1 stable isotope of gold Au$_{79}^{197}$
+- 3 natural isotopes of lead
+ - Pb$_{82}^{204}$ ⚛ -3 protons, -4 neutrons
+ - Pb$_{82}^{206}$ ⚛ -3 protons, -6 neutrons
+ - Pb$_{82}^{207}$ ⚛ -3 protons, -7 neutrons
+ - Pb$_{82}^{208}$ ⚛ -3 protons, -8 neutrons
+- 2 reactions for converting gold to lead
+ - Pb204 (n,3npa) Au197
+ - Pb204 (n,nta) Au197
+- No cross section data found in ENDF
+
+
+
+
+---
+
+# Q values
+
+Amount of energy absorbed (-ve) or release (+ve) during the nuclear reaction
+
+| Reaction | Energy release [MeV] |Threshold reaction |
+| -------- | ------- |------- |
+| Be9(n,2n) | -1.6 | Yes |
+| Pb208(n,2n) | -7.3 | Yes |
+| Li6(n,t) | 4.8 | No |
+| Li7(n,nt) | -2.4 | Yes |
+
+Mass and Binding energy converted to kinetic energy
+
+Online Q value calculator at [NNDC](https://www.nndc.bnl.gov/qcalc/)
+
+---
+
+# Fusion fuels
+
+
+
+
+![](images/fusion-cross-sections.png)
+
+
+
+
+
+Q values of fusion fuel reactions
+
+| Reaction | Energy release (MeV) |
+| --------- | ------- |
+| D + T -> He$^{4}$ + n | 17.6 |
+| D + D -> He$^{3}$+n | 3.3 |
+| D + D -> T + p | 4.0 |
+| D + He$^{3}$->He$^{4}$+p | 18.3 * |
+
+* No neutron emitted
+
+---
+
+# Aneutronic Fusion fuels
+
+
+
+
+- Neutrons are not emitted in the primary fuel reaction
+- Neutrons can be emitted by reactions with the products
+- Energy capture via direct conversion or divertor?
+
+
+
+- A DT plasma has several fusion reactions.
+- DT is the most likely reaction.
+- DD and TT reactions also occur with lower probabilities.
+- All reactions and emit different energy neutrons.
+
+
+
+
+![](images/dd_tt_dt.png)
+
+
+
+
+---
+
+# Microscopic Cross Section
+
+
+
+
+- Measured in Barns (1 barn = $10^{-28}m^{2}$)
+- Energy dependant
+- Cross section evaluations exist for:
+ - different incident particles
+ - different nuclides
+ - different interactions.
+- Important neutron reactions plotted
+ - Tritium breeding
+ - Neutron multiplication
+
+
+
+---
+
+# Reaction rate equation
+
+
+- The reaction rate ($RR$) can be found by knowing the number of neutrons per unit volume ($n$), the velocity of neutrons ($v$), the material density ($p$), Avogadro's number ($N_{a}$), the microscopic cross section at the neutron energy ($\sigma_{e}$) and the atomic weight of the material ($M$).
+- This reduces down to the neutron flux ($\phi$), nuclide number density ($N_{d}$) and microscopic cross section $\sigma_{e}$.
+- This can be reduced one more stage by making use of the Macroscopic cross section ($\Sigma_{e}$).
+
+
+$$ RR = \frac{nv\rho N_{a}\sigma_{e} }{M} = \phi N_{d} \sigma_{e} = \phi \Sigma_{e} $$
+
+---
+
+# Macroscopic cross section
+
+
+
+
+
+- Lithium metatitanate has a material density of 3.4 g/cm3
+- When plotting materials the Macroscopic cross section accounts for number density of the different isotopes
+- Units for Macroscopic cross section are cm$^{-1}$
+
+
+
+
+![](images/macroscopic_cross_sections.png)
+
+
+
+
+---
+
+# Multigroup cross sections
+
+
+
+
+- Discretize a continuous distribution
+- Histogram of average cross section in each energy bin
+- Continuous cross section has rules for interpolation that can be accounted for.
+- Groups are not equally spaced.
+- Structures are optimized for different energy ranges (fission, fast fission, fusion etc)
+
+
+
+- The scattering angle varies depending on the energy of the incident neutron
+- Low energy neutrons have isotropic scattering (even probability in all directions)
+- High energy neutrons are more likely to have a low deflection angle and are forwards bias.
+
+
+
+
+![](images/angle_energy_cross_section.png)
+
+
+
+
+---
+
+# Energy distribution
+
+
+
+
+- There is also data on neutrons released in reactions such as (n,2n).
+- The (n,2n) reaction is a threshold reaction and requires energy.
+- No run away chain reaction possible.
+
+
+
+
+![](images/angle_energy_be9.png)
+
+
+
+
+---
+
+# Experimental data
+
+
+
+
+- Availability of experimental data varies for different reactions and different isotopes.
+
+- Typically the experimental data is then interpreted to create evaluation libraries, such as ENDF, JEFF, JENDL, CENDL.
+
+
+
+
+
+[![](images/exfor_be_n_2n.png)](https://nds.iaea.org/dataexplorer/)
+
+Source [IAEA nuclear data services](https://nds.iaea.org/dataexplorer/?target_elem=Au&target_mass=197&reaction=n%2Cg)
+
+
+
+
+---
+
+# Nuclear data libraries
+
+There are several groups that produce and distribute nuclear data
+
+- TENDL 2023 🇪🇺 2850 neutron
+- JENDL 5 🇯🇵 795 neutron
+- ENDF/B-VIII.0 🇺🇸 557 neutron
+- JEFF 3.3 🇪🇺 562 neutrons
+- BROND 3.1 🇷🇺 372 neutrons
+- FENDL 3.2b 🌐 191 neutron
+- CENDL 3.2 🇨🇳 272 neutron
+
+---
+
+# Path length
+
+
+
+
+- Path length = 1 / $\Sigma_{T}$
+- A 14MeV neutron will lose energy via scattering interactions
+- As the neutron energy decreases the path length also decreases
+- Path length at thermal energy is more constant
+
+![](images/neutron-scatter.png)
+
+
+
+---
+
+# Energy loss
+
+The average logarithmic energy decrement (or loss) per collision ($\xi$) is related to the atomic mass ($A$) of the nucleus
+
+
+
+- Lithium has a particularly high cross section for tritium production
+- Li6 has a very high cross section at low neutron energies
+- Li7 has a reasonable cross section at high neutron energies
+- Other reaction channels are relativity low
+- Often alloyed with Si or other elements to improve material properties (e.g. flammability)
+
+
+
+
+![](images/all_tritium_multi.png)
+
+* Elements up to Iron plotted
+
+
+
+---
+
+# Why beryllium
+
+
+
+
+- Beryllium has the lowest threshold energy for any isotope with a n,2n reaction.
+- This means even low energy 3MeV neutrons can undergo (n,2n) reactions.
+- Often alloyed with Ti or other elements to improve material properties (e.g. swelling due to retention)
+- Lead is also a popular choice for a neutron multiplier
+
+
+
+
+![](images/all_neutron_multi.png)
+* Elements up to Iron plotted
+
+
+
+
+---
+
+# Other materials
+
+
+## Tungsten
+
+- High atomic number = good gamma attenuation
+
+- High neutron capture resonances = good neutron attenuation
+
+## Water
+
+- High hydrogen content = excellent neutron moderator
+
+## Helium 4
+
+- Low interaction cross sections and low density = transparent to neutrons and gammas
+
+
+---
+
+# Neutron spectra through materials
+
+
+By knowing the materials present can you identify which blanket results in which spectrum
+
+- FLiBe, Molten salt, typically 90% enriched Li6
+- HCPB, helium cooled pebble bed, typically 60% enriched Li6
+- HCLL, helium cooled lithium lead, typically 90% enriched Li6
+- WCCB, Water cooled ceramic breeder, typically 60% enriched Li6
+- WCLL, water cooled lithium lead, typically 90% enriched Li6
+- Liquid Lithium, typically natural enrichment
+
+
+---
+
+
+![bg](images/who-wants-to-be-million-spectra.png)
+
+---
+
+![](images/blanket_spectra_legend.png)
+
+---
+
+# Prompt responses
+
+
+
+
+ - Neutron wall loading
+ - Heating
+
+
+
+
+
+
+
+ - Tritium breeding
+
+
+ - Dose
+
+
+
+---
+
+# Neutron wall loading
+
+
+
+
+- Energy carried by uncollided source neutrons incident on a unit area of first wall per unit time
+- Units typically used $MW m^{-2}$
+- Useful for estimating neutronics results and scaling or comparing results
+- For simple source distributions and geometry, can calculate analytically
+- Complex source distributions or geometries require more sophisticated methods (e.g Monte Carlo)
+
+
+
+
+
+
+![](scripts/wall-loading.png)
+
+
+
+
+
+
+---
+
+# Neutron wall example
+
+- Significant poloidal variation of neutron wall loading occur in toroidal magnetic confinement fusion reactors
+
+![width:650px](images/wall-loading.png)
+source http://dx.doi.org/10.13182/FST13-751
+
+---
+
+# Nuclear Heating
+
+- Energy deposition calculated from the flux using “Kinetic Energy Released in MAterials” (KERMA) factors
+- Energy lost by a neutron from a collision is assumed to be deposited locally
+- Gamma photons produced by neutrons are transported to determine where their energy is deposited (need coupled neutron-photon transport)
+- The power density distribution is used in thermal-hydraulics calculations and subsequent structural analysis (e.g. thermal stress)
+- Total heating is used for sizing cooling systems
+- Nuclear energy multiplication (Mn) is ratio of energy deposited by neutrons and gamma photons in the reactor to neutron energy incident on FW
+
+
+---
+
+# Nuclear Heating depends on material and location
+
+
+
+
+
+- At same location with same neutron flux, nuclear heating depends on material
+- High-Z materials usually yield higher nuclear heating than low-Z materials
+- Gamma heating represents ~85% of nuclear heating in high-Z materials and only ~40% in low-Z materials
+- Nuclear heating drops rapidly as we move away from FW
+
+
+
+- Different types of dose, absorbed, equivalent and effective.
+- Effective dose is typically used for dose maps.
+- Dose coefficients units of $Sv.cm^2$
+- Neutron flux ($particles.cm^{-2}s^{-1}$)
+- Resulting dose in Sv per second
+
+
+
+
+
+![](images/dose.png)
+
+
+
+
+---
+
+
+---
+
+
+# Delayed response
+
+
+
+
+
+ - Activation
+ - Activity build up and decay
+ - Emission spectra
+ - Shut down dose
+
+
+
+
+
+
+ - Analysis needed to lift or cool components
+
+ - Impact of burn up on TBR
+
+
+
+
+
+
+
+- New isotopes created during irradiation
+
+
+
+
+
+- Radioactive isotopes decay and will eventually reach a point where decay rate is equal to activation rate.
+
+
+
+
+
+- Decay is more noticeable once the plasma is shutdown.
+
+
+
+- The activity is related to the irradiation time and the nuclide half life.
+
+
+
+- High energy neutron activation
+![](images/isotope_chart_high_activation.png)
+
+
+
+
+- Low energy neutron activation
+![](images/isotope_chart_low_activation.png)
+
+
+
+
+---
+
+
+# Activation products from fission
+
+- Fission of large atoms (e.g. U235)
+- Results in two fission products far from stability
+
+![bg 60%](images/isotope_chart_fission_activation.png)
+
+
+---
+
+# Emission during decay
+
+
+
+
+
+- Characteristic gamma energies and intensities emitted
+- Reduces with half life of unstable isotope
+- Problematic sources in fusion Co60
+- Neutrons also emitted by isotopes such as N17 found which is formed by Oxygen irradiation in water
+
+
+
+
+![](images/gamma_spec2.png)
+
+
+
+
+---
+
+# Shut down dose rate
+
+
+
+
+- Post irradiation gamma and even neutron emission from radioactive isotopes continues.
+- Gamma and neutrons emitted cause dose field that makes human maintenance difficult.
+- This causes components to generate self heating
+- Reduced strength of components due to temperature, lift carefully
+- Activated coolant pumped outside of the bio-shield
+
+
+
+Image source [Eurofusion](https://scipub.euro-fusion.org/wp-content/uploads/eurofusion/WPBBPR17_17590_submitted.pdf)
+
+
+
+---
+
+![bg](images/million-solid.png)
+
+
+
+---
+
+# Overview of neutronics simulation software
+
+- Inventory codes
+- Monte Carlo Radiation transport
+- Geometry conversion software
+
+---
+
+# Inventory codes
+
+Solving the Bateman equation
+
+| Name of software | Group / community / country |
+| ----- | -----|
+| ACAB | UNED, Spain |
+| ALARA | Wisconsin, US |
+| Aburn | North China Electrical Power |
+| OpenMC | MIT, ANL, community |
+| Origen | LANL, US |
+| Serpent | VTT, Finland |
+| Fispact | CCFE / UKAEA |
+| Fornax | Silver Fir Software, US |
+
+
+---
+
+# Radiation transport
+
+Sampling the Boltzman transport equation
+
+- Stochastic / Monte Carlo is most widely used method in fusion
+- Track individual particle histories through phase space
+- Random sampling of particle behavior at each event
+- Accumulate contributions to the mean behavior from each history
+- Variance reduction used to speed up simulation
+
+
+---
+
+# Monte Carlo Simulations
+
+| Name of software | Group / community / country |
+| ----- | -----|
+| FLUKA | CERN |
+| GEANT | CERN |
+| MCNP | LANL |
+| OpenMC | MIT, ANL and open source community |
+| Serpent | VTT, Finland |
+| TopMC | China |
+| TRIPOLI | France |
+| SCONE | Cambridge UK |
+| MC DC | US |
+
+---
+
+# Geometry for Monte Carlo
+
+
+
+
+
+---
+
+
+# Geometry conversion
+
+[Link to flowchat](https://www.mermaidchart.com/raw/bfea01f7-56e6-4780-9687-0a6c99e58b74?theme=light&version=v0.1&format=svg)
+
+![mermaid](images/cad-toneutronics-routes.png)
+
+
+---
+
+# Software distribution
+
+Open source codes such as OpenMC and DAGMC are distributed via GitHub, conda.
+
+Some codes used in neutronics are controlled codes under export control
+
+Distribution in the US by RSICC and in the EU by the NEA databank.
+
+![RSICC](images/rsicc.png)
+[RSICC](http://rsicc.ornl.gov/Default.aspx)
+
+
+
+---
+
+# Questions
+
+# 📧 mail@jshimwell.com
+# ![width:60](images/github.png) @shimwell
+
+
+---
\ No newline at end of file
diff --git a/all_slides.pdf b/all_slides.pdf
new file mode 100644
index 0000000..3a954ed
Binary files /dev/null and b/all_slides.pdf differ
diff --git a/index.html b/index.html
index bd31ab9..5175e96 100644
--- a/index.html
+++ b/index.html
@@ -285,7 +285,7 @@
Reactions
text-align: center
}
;background-color:#fff;background-image:none;" data-marpit-pagination-total="55">
-
Neutron induced reactions
+
Neutron induced reactions
@@ -330,7 +330,7 @@
Neutron induced reactions
text-align: center
}
;background-color:#fff;background-image:none;" data-marpit-pagination-total="55">
-
Transmutation reactions
+
Transmutation reactions
Reactions that result in a change of the isotope
@@ -380,7 +380,7 @@
Isotope transmutation
text-align: center
}
;background-color:#fff;background-image:none;" data-marpit-pagination-total="55">
-
Transmutation of lead to gold
+
Transmutation of lead to gold
@@ -435,7 +435,7 @@
Transmutation of lead to gold
text-align: center
}
;background-color:#fff;background-image:none;" data-marpit-pagination-total="55">
-
Q values
+
Q values
Amount of energy absorbed (-ve) or release (+ve) during the nuclear reaction
@@ -498,7 +498,7 @@
Q values
text-align: center
}
;background-color:#fff;background-image:none;" data-marpit-pagination-total="55">
-
Fusion fuels
+
Fusion fuels
@@ -562,7 +562,7 @@
Fusion fuels
text-align: center
}
;background-color:#fff;background-image:none;" data-marpit-pagination-total="55">
-
Aneutronic Fusion fuels
+
Aneutronic Fusion fuels
@@ -640,7 +640,7 @@
Aneutronic Fusion fuels
text-align: center
}
;background-color:#fff;background-image:none;" data-marpit-pagination-total="55">
-
Energy of neutrons from DT fuel
+
Energy of neutrons from DT fuel
@@ -682,7 +682,7 @@
Energy of neutrons from DT fuel
text-align: center
}
;background-color:#fff;background-image:none;" data-marpit-pagination-total="55">
-
Microscopic Cross Section
+
Microscopic Cross Section
@@ -735,7 +735,7 @@
Microscopic Cross Section
text-align: center
}
;background-color:#fff;background-image:none;" data-marpit-pagination-total="55">
-
Reaction rate equation
+
Reaction rate equation
The reaction rate () can be found by knowing the number of neutrons per unit volume (), the velocity of neutrons (), the material density (), Avogadro's number (), the microscopic cross section at the neutron energy () and the atomic weight of the material ().
This reduces down to the neutron flux (), nuclide number density () and microscopic cross section .
@@ -769,7 +769,7 @@
Reaction rate equation
text-align: center
}
;background-color:#fff;background-image:none;" data-marpit-pagination-total="55">
-
Macroscopic cross section
+
Macroscopic cross section
@@ -810,7 +810,7 @@
Macroscopic cross section
text-align: center
}
;background-color:#fff;background-image:none;" data-marpit-pagination-total="55">
-
Multigroup cross sections
+
Multigroup cross sections
@@ -853,7 +853,7 @@
Multigroup cross sections
text-align: center
}
;background-color:#fff;background-image:none;" data-marpit-pagination-total="55">
-
Cross section regions
+
Cross section regions
Reactions have characteristics
@@ -938,7 +938,7 @@
Angular distribution
text-align: center
}
;background-color:#fff;background-image:none;" data-marpit-pagination-total="55">
-
Energy distribution
+
Energy distribution
@@ -1024,7 +1024,7 @@
Experimental data
text-align: center
}
;background-color:#fff;background-image:none;" data-marpit-pagination-total="55">
-
Nuclear data libraries
+
Nuclear data libraries
There are several groups that produce and distribute nuclear data
TENDL 2023 2850 neutron
@@ -1166,7 +1166,7 @@
Energy loss
text-align: center
}
;background-color:#fff;background-image:none;" data-marpit-pagination-total="55">
-
Why lithium
+
Why lithium
@@ -1212,7 +1212,7 @@
Why lithium
text-align: center
}
;background-color:#fff;background-image:none;" data-marpit-pagination-total="55">
-
Why beryllium
+
Why beryllium
@@ -1257,7 +1257,7 @@
Why beryllium
text-align: center
}
;background-color:#fff;background-image:none;" data-marpit-pagination-total="55">
-
Other materials
+
Other materials
Tungsten
@@ -1303,7 +1303,7 @@
Helium 4
text-align: center
}
;background-color:#fff;background-image:none;" data-marpit-pagination-total="55">
-
Neutron spectra through materials
+
Neutron spectra through materials
By knowing the materials present can you identify which blanket results in which spectrum
FLiBe, Molten salt, typically 90% enriched Li6
@@ -1639,7 +1639,7 @@
Nuclear Heating depend
text-align: center
}
;background-color:#fff;background-image:none;" data-marpit-pagination-total="55">
-