From 7dfcc480ba1e19bd3232349fc733caef94034292 Mon Sep 17 00:00:00 2001 From: stainer_t Date: Mon, 8 Sep 2025 13:48:49 +0200 Subject: Initial commit from Polytechnique Montreal --- doc/IGE335/Section3.05_shi.tex | 137 +++++++++++++++++++++++++++++++++++++++++ 1 file changed, 137 insertions(+) create mode 100644 doc/IGE335/Section3.05_shi.tex (limited to 'doc/IGE335/Section3.05_shi.tex') diff --git a/doc/IGE335/Section3.05_shi.tex b/doc/IGE335/Section3.05_shi.tex new file mode 100644 index 0000000..92f7f22 --- /dev/null +++ b/doc/IGE335/Section3.05_shi.tex @@ -0,0 +1,137 @@ +\subsection{The {\tt SHI:} module}\label{sect:SHIData} + +The {\tt SHI:} module perform self-shielding calculations in DRAGON, using +the generalized Stamm'ler method.\cite{SHIBA} This approach is based on an heterogeneous-homogeneous equivalence principle. In this case, an {\sl equivalent dilution parameter} $\sigma_{{\rm e},g}$ is computed for each resonant isotope, in each resonant region and +each resonant energy group $g$. This dilution parameter is used to interpolate pretabulated effective cross sections for the infinite homogeneous medium, previously obtained with the {\sl flux calculator} of the {\tt GROUPR} module in code NJOY.\cite{njoy2010} +Each resonant isotope, identified as such by the \dusa{inrs} +parameter defined in \Sect{LIBData}, is to be recalculated. The general format of +the data for this module is: + +\begin{DataStructure}{Structure \dstr{SHI:}} +\dusa{MICLIB} \moc{:=} \moc{SHI:} $\{$ \dusa{MICLIB} $|$ \dusa{OLDLIB} $\}$ +\dusa{TRKNAM} $[$ \dusa{TRKFIL} $]$ \moc{::} \dstr{descshi} +\end{DataStructure} + +\noindent +where + +\begin{ListeDeDescription}{mmmmmmmm} + +\item[\dusa{MICLIB}] {\tt character*12} name of the \dds{microlib} that will +contain the microscopic and macroscopic cross sections updated by the +self-shielding module. If +\dusa{MICLIB} appears on both LHS and RHS, it is updated; otherwise, the +internal library \dusa{OLDLIB} is copied into +\dusa{MICLIB} and \dusa{MICLIB} is updated. + +\item[\dusa{OLDLIB}] {\tt character*12} name of a read-only \dds{microlib} +that is copied into \dusa{MICLIB}. + +\item[\dusa{TRKNAM}] {\tt character*12} name of the required \dds{tracking} +data structure. + +\item[\dusa{TRKFIL}] {\tt character*12} name of the sequential binary tracking +file used to store the tracks lengths. This file is given if and only if it was +required in the previous tracking module call (see \Sect{TRKData}). + +\item[\dstr{descshi}] structure describing the self-shielding options. + +\end{ListeDeDescription} + +Each time the \moc{SHI:} module is called, a sub-directory is updated in the +\dds{microlib} data structure to hold the last values defined in the +\dstr{descshi} structure. The next time this module is called, +these values will be used as floating defaults. + +\subsubsection{Data input for module {\tt SHI:}}\label{sect:descshi} + +\begin{DataStructure}{Structure \dstr{descshi}} +$[$ \moc{EDIT} \dusa{iprint} $]$ \\ +$[$ \moc{GRMIN} \dusa{lgrmin} $]~~[$ \moc{GRMAX} \dusa{lgrmax} $]$ \\ +$[$ \moc{MXIT} \dusa{imxit} $]~~[$ \moc{EPS} \dusa{valeps} $]$ \\ +$[~\{$ \moc{LJ} $|$ \moc{NOLJ} $\}~]$ $[~\{$ \moc{GC} $|$ \moc{NOGC} $\}$ +$[~\{$ \moc{TRAN} $|$ \moc{NOTR} $\}~]$ +$[$ \moc{LEVEL} \dusa{ilev} $]$ \\ +$[~\{$ \moc{PIJ} $|$ \moc{ARM} $\}~]$ \\ +{\tt ;} +\end{DataStructure} + +\noindent where + +\begin{ListeDeDescription}{mmmmmmmm} + +\item[\moc{EDIT}] keyword used to modify the print level \dusa{iprint}. + +\item[\dusa{iprint}] index used to control the printing of this module. The +amount of output produced by this tracking module will vary substantially +depending on the print level specified. + +\item[\moc{GRMIN}] keyword to specify the minimum group number considered +during the self-shielding process. + +\item[\dusa{lgrmin}] first group number considered during the +self-shielding process. By default, \dusa{lgrmin} is set to the first group +number containing self-shielding data in the library. + +\item[\moc{GRMAX}] keyword to specify the maximum group number considered +during the self-shielding process. + +\item[\dusa{lgrmax}] last group number considered during the self-shielding +process. By default, \dusa{lgrmax} is set is set to the last group +number containing self-shielding data in the library. + +\item[\moc{MXIT}] keyword to specify the maximum number of iterations during +the self-shielding process. + +\item[\dusa{imxit}] the maximum number of iterations. The default is +\dusa{imxit}=20. + +\item[\moc{EPS}] keyword to specify the convergence criterion for the +self-shielding iteration. + +\item[\dusa{valeps}] the convergence criterion for the self-shielding iteration. +By default, \dusa{valeps}=$1.0\times 10^{-4}$. + +\item[\moc{LJ}] keyword to activate the Livolant and Jeanpierre +normalization scheme which modifies the self-shielded averaged neutron fluxes in +heterogeneous geometries. By default the Livolant and Jeanpierre +normalization scheme is not activated. + +\item[\moc{NOLJ}] keyword to deactivate the Livolant and Jeanpierre +normalization scheme which modifies the self-shielded averaged neutron fluxes in +heterogeneous geometries. This is the default option. + +\item[\moc{GC}] keyword to activate the Goldstein-Cohen approximation in +cases where Goldstein-Cohen parameters are stored on the internal library. These +parameters may not be available with some libraries (e.g., {\tt APLIB1}, {\tt +APLIB2} or MATXS-type libraries). The Goldstein-Cohen parameters can always be +imposed using the \moc{IRSET} keyword of the \moc{LIB:} module (see +\Sect{LIBData}). This is the default option. + +\item[\moc{NOGC}] keyword to deactivate the Goldstein-Cohen approximation in +cases where Goldstein-Cohen parameters are stored on the internal library. + +\item[\moc{TRAN}] keyword to activate the transport correction option for +self-shielding calculations (see \moc{CTRA} in \Sectand{MACData}{LIBData}). This is the default option. + +\item[\moc{NOTR}] keyword to deactivate the transport correction option for +self-shielding calculations (see \moc{CTRA} in \Sectand{MACData}{LIBData}). + +\item[\moc{LEVEL}] keyword to specify the self-shielding model. + +\item[\dusa{ilev}] $=0$: original Stamm'ler model (without distributed effects); + $=1$: use the Nordheim (PIC) distributed self-shielding model\cite{toronto04}; + $=2$: use both Nordheim (PIC) distributed self-shielding model and Riemann integration + method\cite{hasan}. By default, \dusa{ilev}$\,=0$. + +\item[\moc{PIJ}] keyword to specify the use of complete collision +probabilities in the self-shielding calculations of {\tt SHI:}. +This is the default option for \moc{EXCELT:} and \moc{SYBILT:} trackings. +This option is not available for \moc{MCCGT:} trackings. + +\item[\moc{ARM}] keyword to specify the use of iterative flux techniques +in the self-shielding calculations of {\tt SHI:}. +This is the default option for \moc{MCCGT:} trackings. + +\end{ListeDeDescription} +\eject -- cgit v1.2.3