Initial commit from Polytechnique Montreal

1 files changed, 304 insertions, 0 deletions

diff --git a/doc/IGE344/SectDTHM.tex b/doc/IGE344/SectDTHM.tex
new file mode 100644
index 0000000..c546b00
--- /dev/null
+++ b/doc/IGE344/SectDTHM.tex
@@ -0,0 +1,304 @@
+\subsection{Contents of \dir{thm} data structure}\label{sect:thmdir}
+
+This data structure contains the thermal-hydraulics information required in a multi-physics calculation
+
+\subsubsection{The main \dir{thm} directory}\label{sect:thmdirmain}
+
+The following records and sub-directories will be found in the first level of a \dir{thm} directory:
+
+\begin{DescriptionEnregistrement}{Main records and sub-directories in \dir{thm}}{8.0cm}
+\CharEnr
+  {SIGNATURE\blank{3}}{$*12$}
+  {parameter $\mathsf{SIGNA}$ containing the signature of the data structure}
+\IntEnr
+  {STATE-VECTOR}{$40$}
+  {array $\mathcal{S}^{th}_{i}$ containing various integer parameters that are required to describe this data structure}
+\RealEnr
+  {REAL-PARAM\blank{2}}{$40$}{}
+  {array $\mathcal{R}^{th}_{i}$ containing various floating-point parameters that are required to describe this data structure}
+\RealEnr
+  {MESHZ\blank{7}}{$\mathcal{S}^{th}_{1}$}{m}
+  {initial axial meshes as recovered from the fuelmap}
+\RealEnr
+  {REF-RAD\blank{5}}{$(\mathcal{S}^{th}_{7}-1)\times\mathcal{S}^{th}_{1}$}{m}
+  {initial radial meshes as recovered from the first call to {\tt THM:} module}
+\RealEnr
+  {NB-FUEL\blank{5}}{$\mathcal{S}^{th}_{2}$}{}
+  {number of active fuel rods in a single assembly or number of active fuel pins in the cluster}
+\RealEnr
+  {NB-TUBE\blank{5}}{$\mathcal{S}^{th}_{2}$}{}
+  {number of active guide tubes in a single assembly}
+\RealEnr
+  {FRACT-PU\blank{4}}{$\mathcal{S}^{th}_{2}$}{}
+  {plutonium mass enrichment}
+\OptRealEnr
+  {KCONDF\blank{6}}{$\mathcal{S}^{th}_{16}+3$}{$\mathcal{S}^{th}_{12}\ne 0$}{}
+  {coefficients of the user-defined correlation for the fuel thermal conductivity}
+\OptCharEnr
+  {UCONDF\blank{6}}{$12$}{$\mathcal{S}^{th}_{12}\ne 0$}
+  {string variable set to {\tt KELVIN} or to {\tt CELSIUS}}
+\OptRealEnr
+  {KCONDC\blank{6}}{$\mathcal{S}^{th}_{17}+3$}{$\mathcal{S}^{th}_{13}\ne 0$}{}
+  {coefficients of the user-defined correlation for the clad thermal conductivity}
+\OptCharEnr
+  {UCONDC\blank{6}}{$12$}{$\mathcal{S}^{th}_{13}\ne 0$}
+  {string variable set to {\tt KELVIN} or to {\tt CELSIUS}}
+\RealEnr
+  {ERROR-T-FUEL}{1}{K}
+  {absolute error in fuel temperature}
+\RealEnr
+  {ERROR-D-COOL}{1}{g/cc}
+  {absolute error in coolant density}
+\RealEnr
+  {ERROR-T-COOL}{1}{K}
+  {absolute error in coolant temperature}
+\RealEnr
+  {ERROR-P-COOL}{1}{Pa}
+  {absolute error in coolant pressure}
+\RealEnr
+  {MIN-T-FUEL\blank{2}}{1}{K}
+  {minimum fuel temperature}
+\RealEnr
+  {MIN-D-COOL\blank{2}}{1}{g/cc}
+  {minimum coolant density}
+\RealEnr
+  {MIN-T-COOL\blank{2}}{1}{K}
+  {minimum coolant temperature}
+\RealEnr
+  {MIN-P-COOL\blank{2}}{1}{Pa}
+  {minimum coolant pressure}
+\RealEnr
+  {MAX-T-FUEL\blank{2}}{1}{K}
+  {maximum fuel temperature}
+\RealEnr
+  {MAX-D-COOL\blank{2}}{1}{g/cc}
+  {maximum coolant density}
+\RealEnr
+  {MAX-T-COOL\blank{2}}{1}{K}
+  {maximum coolant temperature}
+\RealEnr
+  {MAX-P-COOL\blank{2}}{1}{Pa}
+  {maximum coolant pressure}
+ \DirEnr
+  {HISTORY-DATA}
+  {sub-directory containing the historical values taken by the thermal-hydraulics parameters (mass flux, density, pressure, enthalpy, temperature)
+  in the coolant and in the fuel rod for the whole geometry}
+\OptRealEnr
+  {RAD-PROF\_R\blank{2}}{$\mathcal{S}^{th}_{18}$}{$\mathcal{S}^{th}_{18}\ne 0$}{m}
+  {abscissas of the user-defined radial form factor table}
+\OptRealEnr
+  {RAD-PROF\_F\blank{2}}{$\mathcal{S}^{th}_{18}$}{$\mathcal{S}^{th}_{18}\ne 0$}{ }
+  {form-factor values of the user-defined radial form factor table}
+\OptRealEnr
+  {TIME-SR1\blank{2}}{$\mathcal{S}^{th}_{19}$}{$\mathcal{S}^{th}_{19}\ne 0$}{s}
+  {tabulation abscissa in time}
+\OptRealEnr
+  {POWER-SR1\blank{2}}{$\mathcal{S}^{th}_{19}$}{$\mathcal{S}^{th}_{19}\ne 0$}{ }
+  {tabulation power factor corresponding to each tabulation abscissa in time}
+\end{DescriptionEnregistrement}
+
+The signature for this data structure is $\mathsf{SIGNA}$=\verb*|L_THM|. The array $\mathcal{S}^{h}_{i}$
+contains the following information: 
+
+\begin{itemize}
+\item $\mathcal{S}^{th}_{1}$ contains the number of axial meshes $N_{\rm z}$.
+\item $\mathcal{S}^{th}_{2}$ contains the number of channels in the radial plane $N_{\rm ch}$.
+\item $\mathcal{S}^{th}_{3}$ contains the maximum number of iterations for computing the
+conduction integral.
+\item $\mathcal{S}^{th}_{4}$ contains the maximum number of iterations for computing the
+center pellet temperature.
+\item $\mathcal{S}^{th}_{5}$ contains the maximum number of iterations for computing the
+coolant parameters (velocity, pressure, enthapy, density) in case of a transient calculation.
+\item $\mathcal{S}^{th}_{6}$ contains the number of discretisation points in fuel.
+\item $\mathcal{S}^{th}_{7}$ contains the number of total discretisation points in the whole fuel rod (fuel+cladding) $N_{\rm dtot}$.
+\item $\mathcal{S}^{th}_{8}$ contains the type of calculation performed by the \moc{THM:} module:
+
+\begin{displaymath} \mathcal{S}^{th}_{8} = \left\{
+\begin{array}{rl}
+ 0 & \textrm{steady-state calculation} \\
+ 1 & \textrm{transient calculation.} \\
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{th}_{9}$ contains the current time index.
+\item $\mathcal{S}^{th}_{10}$ flag to set the gap correlation:
+
+\begin{displaymath} \mathcal{S}^{th}_{10} = \left\{
+\begin{array}{rl}
+ 0 & \textrm{built-in correlation is used} \\
+ 1 & \textrm{set the heat exchange coefficient of the gap as a user-defined constant.} \\
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{th}_{11}$ flag to set the heat transfer coefficient between the clad and fluid:
+
+\begin{displaymath} \mathcal{S}^{th}_{11} = \left\{
+\begin{array}{rl}
+ 0 & \textrm{built-in correlation is used} \\
+ 1 & \textrm{set the heat exchange coefficient between the clad and fluid as a user-defined constant.} \\
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{th}_{12}$ flag to set the fuel thermal conductivity:
+
+\begin{displaymath} \mathcal{S}^{th}_{12} = \left\{
+\begin{array}{rl}
+ 0 & \textrm{built-in correlation is used} \\
+ 1 & \textrm{set the fuel thermal conductivity as a function of a simple user-defined correlation.} \\
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{th}_{13}$ flag to set the clad thermal conductivity:
+
+\begin{displaymath} \mathcal{S}^{th}_{13} = \left\{
+\begin{array}{rl}
+ 0 & \textrm{built-in correlation is used} \\
+ 1 & \textrm{set the clad thermal conductivity as a function of a simple user-defined correlation.} \\
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{th}_{14}$ type of approximation used during the fuel conductivity evaluation:
+
+\begin{displaymath} \mathcal{S}^{th}_{14} = \left\{
+\begin{array}{rl}
+ 0 & \textrm{use a rectangle quadrature approximation} \\
+ 1 & \textrm{use an average approximation.} \\
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{th}_{15}$ type of subcooling model:
+
+\begin{displaymath} \mathcal{S}^{th}_{15} = \left\{
+\begin{array}{rl}
+ 0 & \textrm{use the Bowring correlation} \\
+ 1 & \textrm{use the Saha-Zuber correlation.} \\
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{th}_{16}$ contains the number of terms in the user-defined correlation for the fuel thermal conductivity (if $\mathcal{S}^{th}_{12}=1$).
+\item $\mathcal{S}^{th}_{17}$ contains the number of terms in the user-defined correlation for the clad thermal conductivity (if $\mathcal{S}^{th}_{13}=1$).
+\item $\mathcal{S}^{th}_{18}$ type of radial form factor for the power:
+
+\begin{displaymath} \mathcal{S}^{th}_{18} = \left\{
+\begin{array}{rl}
+ 0 & \textrm{flat radial form factor} \\
+ N_{\rm rad} & \textrm{number of point in the radial form factor table.} \\
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{th}_{19}$ number of points in the user-defined time-power table.
+
+\item $\mathcal{S}^{th}_{20}$ type of fluid:
+\begin{displaymath} \mathcal{S}^{th}_{20} = \left\{
+\begin{array}{rl}
+ 0 & \textrm{light water (H$_2$O)} \\
+ 1 & \textrm{heavy water (D$_2$O).} \\
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{th}_{21}$ flag indicating if the gap is considered:
+\begin{displaymath} \mathcal{S}^{th}_{21} = \left\{
+\begin{array}{rl}
+ 0 & \textrm{gap is considered} \\
+ 1 & \textrm{is not.} \\
+\end{array} \right.
+\end{displaymath}
+
+\item $\mathcal{S}^{th}_{22}$ flag indicating the pressure drop option:
+\begin{displaymath} \mathcal{S}^{th}_{22} = \left\{
+\begin{array}{rl}
+ 0 & \textrm{no pressure drop} \\
+ 1 & \textrm{pressure drop is computed.} \\
+\end{array} \right.
+\end{displaymath}
+
+\end{itemize}
+
+The array $\mathcal{R}^{th}_{i}$ contains the following information: 
+
+\begin{itemize}
+\item $\mathcal{R}^{th}_{1}$ contains the current time step in s. 
+\item $\mathcal{R}^{th}_{2}$ contains the fraction of reactor power released in fuel.
+\item $\mathcal{R}^{th}_{3}$ contains the critical heat flux in W/m$^2$.
+\item $\mathcal{R}^{th}_{4}$ contains the inlet coolant velocity in m/s.
+\item $\mathcal{R}^{th}_{5}$ contains the outlet coolant pressure in Pa.
+\item $\mathcal{R}^{th}_{6}$ contains the inlet coolant temperature in K.
+\item $\mathcal{R}^{th}_{7}$ contains the fuel porosity.
+\item $\mathcal{R}^{th}_{8}$ contains the fuel pellet radius
+\item $\mathcal{R}^{th}_{9}$ contains the internal clad rod radius in m.
+\item $\mathcal{R}^{th}_{10}$ contains the external clad rod radius in m.
+\item $\mathcal{R}^{th}_{11}$ contains the guide tube radius in m.
+\item $\mathcal{R}^{th}_{12}$ contains the hexagonal side in m. Used only for cluster geometries.
+\item $\mathcal{R}^{th}_{13}$ contains the temperature maximum absolute error (in K) allowed in the solution of the conduction equations.
+\item $\mathcal{R}^{th}_{14}$ contains the maximum relative error allowed in the matrix resolution of the conservation equations of the coolant.
+\item $\mathcal{R}^{th}_{15}$ contains the relaxation parameter for the multiphysics parameters (temperature of fuel and coolant  and density of coolant).
+\item $\mathcal{R}^{th}_{16}$ contains the time in s.
+\item $\mathcal{R}^{th}_{17}$ contains the heat transfer coefficient of the gap (if $\mathcal{S}^{th}_{10}=1$).
+\item $\mathcal{R}^{th}_{18}$ contains the heat transfer coefficient between the clad and fluid (if $\mathcal{S}^{th}_{11}=1$).
+\item $\mathcal{R}^{th}_{19}$ contains the surface temperature weighting factor of effective fuel temperature for the Rowlands approximation.
+\item $\mathcal{R}^{th}_{20}$ reactor power, as defined after the {\tt POWER-LAW} keyword.
+\item $\mathcal{R}^{th}_{21}$ maximum of variable variations in local parameters (used for time step adjustment strategy).
+\item $\mathcal{R}^{th}_{22}$ contains the rugosity of the fuel rod in m, used in M\"uller-Steinhagen correlation for coolant friction.
+\item $\mathcal{R}^{th}_{23}$ contains the angle in radians of the fuel channel with respect of the vertical axis.
+\end{itemize}
+
+\subsubsection{The \moc{HISTORY-DATA} sub-directory}\label{sect:thmdirhistorydata}
+
+In the \moc{HISTORY-DATA} directory, the following sub-directories will be found:
+\begin{DescriptionEnregistrement}{Sub-directories in \moc{HISTORY-DATA} directory}{7.0cm} \label{tabl:tabhistorydatadir}
+ \DirEnr
+ {TIMESTEP0000}
+ {sub-directory containing all the values of the thermal-hydraulics parameters computed by the \moc{THM:} module {\sl in steady-state conditions}.}
+ \DirEnr
+ {TIMESTEP{\sl numt}}
+ {sub-directories containing all the values of the thermal-hydraulics parameters computed by the \moc{THM:} module in transient conditions at a given time index {\sl numt}. {\sl numt} can take values between 1 and 9999 in I4.4 format.}
+\end{DescriptionEnregistrement}
+
+\noindent In the \moc{TIMESTEP0000} and in each of the \moc{TIMESTEP}{\sl numt} sub-directories, the following records will be found:
+\begin{DescriptionEnregistrement}{Records in \moc{TIMESTEP} directories}{7.0cm} \label{tabl:tabtimestepdir}
+ \RealEnr
+ {TIME\blank{8}}{$1$}{$s$}
+ {time}
+ \DirlEnr
+ {\moc{CHANNEL}\blank{5}}{$N_{\rm ch}$}
+ {list of $N_{\rm ch}$ sub-directories containing all the values of the thermal-hydraulics parameters computed by the \moc{THM:} module and sorted channel by channel.}
+\end{DescriptionEnregistrement}
+
+\noindent In each of the \moc{CHANNEL} sub-directories, the following records will be found:
+\begin{DescriptionEnregistrement}{Records in each \moc{CHANNEL} directory}{7.0cm} \label{tabl:tabchanneldir}
+ \RealEnr
+ {VINLET\blank{6}}{$1$}{$m.s^{-1}$}
+ {inlet velocity}
+ \RealEnr
+ {TINLET\blank{6}}{$1$}{$K$}
+ {inlet temperature}
+ \RealEnr
+ {PINLET\blank{6}}{$1$}{$Pa$}
+ {inlet pressure}
+ \RealEnr
+ {VELOCITIES\blank{2}}{$N_{\rm z}$}{$m.s^{-1}$}
+ {velocity in each of the $N_{\rm z}$ bundles of the channel numbered {\sl numc}}
+ \RealEnr
+ {PRESSURE\blank{4}}{$N_{\rm z}$}{$Pa$}
+ {pressure in each bundle of the channel}
+ \RealEnr
+ {ENTHALPY\blank{4}}{$N_{\rm z}$}{$J.kg^{-1}$}
+ {enthalpy in each bundle of the channel}
+ \RealEnr
+ {DENSITY\blank{5}}{$N_{\rm z}$}{$kg.m^{-3}$}
+ {density in each bundle of the channel}
+ \RealEnr
+ {LIQUID-DENS\blank{1}}{$N_{\rm z}$}{$kg.m^{-3}$}
+ {density of liquid phase in each bundle of the channel}
+ \RealEnr
+ {TEMPERATURES}{$N_{\rm z}, N_{\rm dtot}$}{$K$}
+ {distribution of the temperature in the fuel-pin for each bundle of the channel}
+ \RealEnr
+ {CENTER-TEMPS}{$N_{\rm z}$}{$K$}
+ {center fuel pellet temperature in each bundle of the channel}
+ \RealEnr
+ {RADII}{$N_{\rm z}, N_{\rm dtot}-1$}{$m$}
+ {fuel and clad radii}
+\end{DescriptionEnregistrement}
+
+\clearpage