Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Chapter 4 modifications #50

Merged
merged 17 commits into from
Jun 27, 2022
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
14,453 changes: 7,157 additions & 7,296 deletions bibfile.bib

Large diffs are not rendered by default.

19 changes: 10 additions & 9 deletions chapters/chapter1/intro.tex
Original file line number Diff line number Diff line change
Expand Up @@ -238,17 +238,18 @@ \subsection{Plasma-facing materials}
Unfortunately, graphite being porous, hydrogen (and therefore tritium) retention was high \sidecite{sugiyama_measurement_2004}.
Plus carbon can react with \gls{plasma} particles forming methane.
Methane is then deposited on locations hard to access in the reactor \gls{trapping} tritium even more.
For these two safety reasons, \gls{cfc} was replaced with tungsten or beryllium (or both).
For these two safety reasons, \gls{cfc} was replaced with \Gls{W} or \gls{Be} (or both).

\Gls{W} has a very high melting point (\SI{3422}{\celsius}) and retains less tritium \sidecite{pajuste_tritium_2021}.
However, tungsten being a high-Z element, eroded tungsten will make the \gls{plasma} radiate and cool it down.
For this reason, the \acrshort{iter} divertor will be made of tungsten but the first wall (which has a large surface area) will be made of beryllium.
However, \Gls{W} being a high-Z element, eroded \Gls{W} will make the \gls{plasma} radiate and cool it down.
For this reason, the \acrshort{iter} divertor will be made of \Gls{W} but the first wall (which has a large surface area) will be made of \gls{Be}.

\subsection{Divertor}\labsec{divertor section}

In a fusion reactor, heat and particles (fusion ashes) need to be extracted.
In most \glspl{tokamak}, the escaping \gls{plasma} is diverted towards a dedicated component that is heat-resistant.
Such a configuration is called a \emph{\gls{divertor}} configuration.
The \gls{divertor} allows the creation of a \gls{x-point} and a \gls{separatrix} decoupling the \gls{plasma} edge from the walls.
The X-divertor is a common configuration (used in \acrshort{west}, \acrshort{jet}, \acrshort{iter}) but more advanced configurations exist such as the Super-X \gls{divertor} (\acrshort{mast-u}) \sidecite{havlickova_effect_2015}, X-Point Target (\acrshort{sparc}) \sidecite{rodriguez-fernandez_overview_2022, kuang_divertor_2020} or the Snowflake configurations \sidecite{ryutov_snowflake_2015}.

\begin{figure} [h]
Expand All @@ -275,7 +276,7 @@ \subsection{Divertor}\labsec{divertor section}
A prototype of the \gls{ivt} is shown on \reffig{inner target photo}.
These elements are themselves made of rows, called \glspl{pfuLabel}, of small unit bricks of a few dozens of millimetres called \glspl{monoblock}.
In \acrshort{iter}, the \gls{ivt} has 16 \glspl{pfuLabel} and the \gls{ovt} has 22 \glspl{pfuLabel}.
\Glspl{monoblock} are typically made of a tungsten substrate with a cooling pipe running through.
\Glspl{monoblock} are typically made of a \Gls{W} substrate with a cooling pipe running through.
This cooling channel is necessary to keep the component's temperature below its operating limit and exhaust heat.

Several \gls{monoblock} designs are currently studied for DEMO with varying dimensions, different materials for the cooling pipe or the interlayer, etc. \sidecite{vizvary_european_2020, huang_tungsten_2016, hirai_use_2016, domptail_design_2020}.
Expand All @@ -296,7 +297,7 @@ \subsection{Divertor}\labsec{divertor section}
\labfig{inner target photo}
\end{figure}

Studies on ITER-like \glspl{monoblock} have demonstrated the resistance of the \gls{monoblock} design to high heat loads while investigating the effect of tungsten recrystallisation \sidecite{durif_impact_2019, durif_modelisation_2019, visca_manufacturing_2018} (see \reffig{monoblock_temperature_exp_model}).
Studies on ITER-like \glspl{monoblock} have demonstrated the resistance of the \gls{monoblock} design to high heat loads while investigating the effect of \Gls{W} recrystallisation \sidecite{durif_impact_2019, durif_modelisation_2019, visca_manufacturing_2018} (see \reffig{monoblock_temperature_exp_model}).

\begin{figure} [h]
\centering
Expand Down Expand Up @@ -367,7 +368,7 @@ \subsection{Breeding}

The \gls{tbr} is defined by the number of tritium atoms produced by generated neutrons.
In order to ensure tritium self-sufficiency, the \gls{tbr} of the blanket must be greater than or equal to one \sidecite{abdou_blanketfirst_2015}.
A \gls{tbr} greater than one can only be obtained by neutron multiplication with lead or beryllium.
A \gls{tbr} greater than one can only be obtained by neutron multiplication with lead or \gls{Be}.

Moreover, the \gls{tbr} must account for:
\begin{itemize}
Expand Down Expand Up @@ -919,7 +920,7 @@ \section{Problem definition}
Numerous simulation tools have been developed throughout the years (see \reftab{code comparison}).
Most of these codes are not able to run multimaterial and/or multidimensional simulations.
These features are however essential to fully simulate monoblocks (see \refsec{divertor section}).
Many of them rely on the \gls{fdm} whereas \acrshort{hit} \sidecite{candido_integrated_2020}, Abaqus \sidecite{benannoune_multidimensional_2020} and \gls{achlys} \sidecite{stephen-dixon_aurora-multiphysicsachlys_2021} rely on the \gls{fem}.
Many of them rely on the \gls{fdm} whereas \acrshort{hit} \sidecite{candido_integrated_2020}, Abaqus \sidecite{benannoune_multidimensional_2020} and \gls{achlys} \sidecite{dixon_aurora-multiphysicsachlys_2021} rely on the \gls{fem}.
Moreover, some do not have an integrated heat transfer solver - essential for an accurate estimation of temperature fields and therefore thermally activated processes.
Some of these codes rely on proprietary software like Abaqus or COMSOL for \acrshort{hit} - limiting their accessibility and scalability (parallel computing).
The code \gls{achlys} meets all these requirements and is the only one available open-source but was only released in mid 2021.
Expand All @@ -941,13 +942,13 @@ \section{Problem definition}
\\
\acrshort{mhims} \cite{hodille_study_2016} & $\checkmark$ & & & & & & Fortran \\
\\
\acrshort{tessim} \cite{schmid_transport_2014} & $\checkmark$ & & & & & & \\
\acrshort{tessim} \cite{schmid_transport_2014} & $\checkmark$ & & & $\checkmark$ & & & \\
\\
\acrshort{hit} \cite{candido_integrated_2020} & $\checkmark$ & $\checkmark$ & $\checkmark$ & $\checkmark$ & $\checkmark$ & & COMSOL\\
\\
Abaqus \cite{benannoune_multidimensional_2020} & $\checkmark$ & $\checkmark$ & $\checkmark$ & $\checkmark$ & $\checkmark$ & & Fortran\\
\\
\gls{achlys} \cite{stephen-dixon_aurora-multiphysicsachlys_2021} & $\checkmark$ & $\checkmark$ & $\checkmark$ & $\checkmark$ & $\checkmark$ & $\checkmark$ & C++\\
\gls{achlys} \cite{dixon_aurora-multiphysicsachlys_2021} & $\checkmark$ & $\checkmark$ & $\checkmark$ & $\checkmark$ & $\checkmark$ & $\checkmark$ & C++\\
\\
\end{tabular}
\caption{Comparison of some hydrogen transport modelling tools.}
Expand Down
4 changes: 2 additions & 2 deletions chapters/chapter2/model_description.tex
Original file line number Diff line number Diff line change
Expand Up @@ -109,12 +109,12 @@ \subsubsection{Analytical simplification for an implanted source of H} \labsec{t

Plasma implantation of hydrogen particles can be modelled with a volumetric source.
Typically, the depth of the implantation profile is a few nanometres depending on the incident particles energy and incident angle.
These profiles can be simulated by Monte Carlo codes like SRIM \sidecite{ziegler_srim_2010} and have a gaussian-like shape (see \reffig{srim_implantation_profile_example}).
These profiles can be simulated by Monte Carlo codes like \gls{srim} \sidecite{ziegler_srim_2010} and have a gaussian-like shape (see \reffig{srim_implantation_profile_example}).

\begin{figure}
\centering
\includegraphics[width=\linewidth]{Figures/Chapter1/srim_implantation_range.pdf}
\caption{\SI{100}{eV} deuterium implantation profile in tungsten computed from SRIM. Reproduced from \cite{shimada_improved_2019}.}
\caption{\SI{100}{eV} deuterium implantation profile in tungsten computed from \gls{srim}. Reproduced from \cite{shimada_improved_2019}.}
\labfig{srim_implantation_profile_example}
\end{figure}

Expand Down
Loading