1
Fork 0

Add olaflow porosity parameters calibration

This commit is contained in:
Edgar P. Burkhart 2022-05-06 11:39:37 +02:00
parent b5a72a03c5
commit 7ee68c1d5c
Signed by: edpibu
GPG key ID: 9833D3C5A25BD227
4 changed files with 40 additions and 0 deletions

View file

@ -27,3 +27,42 @@ Results are plotted using python \autoref{fig:resola}.
\includegraphics{resola.pdf}
\caption{Results from olaFlow model.}\label{fig:resola}
\end{figure}
\subsection{Porosity parameters}
Several parameters should be calibrated in order to accurately model the porous media: $a$, $b$ and $c$ are friction
parameters in Forcheimer's equation; D50 is the median diameter of the elements constituting the porous media; $p$ is
the porosity of the media.
7 cases were run with the values in \autoref{tab:porotest}.
\begin{table}
\centering
\begin{tabular}{cccccc}
\toprule
\textbf{Case} & $a$ & $b$ & $c$ & D50 (\si{\m}) & $\phi$ \\
\midrule
\textbf{0} & \num{50} & \num{1.2} & \num{0.34} & \num{4} & \num{0.4} \\
\textbf{1} & \intersemibold\num{0} & \num{1.2} & \num{0.34} & \num{4} & \num{0.4} \\
\textbf{2} & \intersemibold\num{5000} & \num{1.2} & \num{0.34} & \num{4} & \num{0.4} \\
\textbf{3} & \num{50} & \intersemibold\num{0} & \num{0.34} & \num{4} & \num{0.4} \\
\textbf{4} & \num{50} & \intersemibold\num{3.0} & \num{0.34} & \num{4} & \num{0.4} \\
\textbf{5} & \num{50} & \num{1.2} & \num{0.34} & \intersemibold\num{2} & \num{0.4} \\
\textbf{6} & \num{50} & \num{1.2} & \num{0.34} & \num{4} & \intersemibold\num{0.25} \\
\bottomrule
\end{tabular}
\caption{Test cases for porosity parameters.}\label{tab:porotest}
\end{table}
Some results are displayed in \autoref{fig:diff} and \autoref{fig:diff2}. No major differences are noticable
between cases (excepted for case 3, where $b=0$).
\begin{figure}
\centering
\includegraphics{diff.pdf}
\caption{Tests for porosity parameters; water - air border at \SI{-50}{\m}.}\label{fig:diff}
\end{figure}
\begin{figure}
\centering
\includegraphics{diff2.pdf}
\caption{Tests for porosity parameters; water - air border at \SI{-20}{\m}.}\label{fig:diff2}
\end{figure}

View file

@ -170,6 +170,7 @@
%inter-unit-product = {.},
per-mode = power-positive-first,
uncertainty-mode = separate,
reset-text-family = false,
}
\SendSettingsToPgf
\DeclareSIUnit\met{met}

BIN
report/fig/diff.pdf Normal file

Binary file not shown.

BIN
report/fig/diff2.pdf Normal file

Binary file not shown.