diff --git a/tex/header/preamble.tex b/tex/header/preamble.tex index 31a420ee1c888c61e4236c0e3b37150a47001a4b..90f7d127cef8cd2472fb11727c0cd1a32482613b 100644 --- a/tex/header/preamble.tex +++ b/tex/header/preamble.tex @@ -14,7 +14,11 @@ \usepackage{silence} % Package that lets ignore warnings and errors for using packages in subdirectories \WarningFilter{latex}{You have requested package `.maindir/tex/header/} -\usepackage[tikz]{.maindir/tex/header/stdmath} % , scriptLetters +\usepackage{.maindir/tex/header/stdmath} % , scriptLetters +\usepackage{listings} +\lstMakeShortInline[columns=fixed]| +\usepackage[section]{placeins} % makes sure that the figures stay within their section +\usepackage{siunitx} \setmainfont{Latin Modern Roman} \setsansfont{Latin Modern Sans} diff --git a/tex/main.tex b/tex/main.tex index 1b0a2ef9cd464573face7249793c0e18eca7003e..3c0a44371ed9a26de8531275bb14c4f4583d0a1b 100644 --- a/tex/main.tex +++ b/tex/main.tex @@ -5,6 +5,7 @@ \author{Felix Hilsky} \maketitle % use \textinput as described in /header and _TEMPLATE +\textinput{setup} \begin{appendices} % \textinput{some-appendix-section} \end{appendices} diff --git a/tex/setup.tex b/tex/setup.tex index f3e7145e97c003f546236c6d0ca90c2baa4fca1d..fd7236f671697830860c6b4c1f64545a9b07f5bf 100644 --- a/tex/setup.tex +++ b/tex/setup.tex @@ -3,6 +3,7 @@ \input{.maindir/tex/header/preamble-section} % inputs the preamble only if necessary \docStart + \section{Mesh} The mesh has two zones: one around the cylinder and a wider one for the remainder of the domain. The zone around the cylinder is a square with arked edges. @@ -29,19 +30,44 @@ In order to keep the cell sizes similar at the corners of the square, the number \end{figure} \section{Boundary values and parameters} +% The turbulence parameters $α_t$, $ν_t$, $k$, $ε$ and $ω$ have a wall function as the boundary value at the cylinder: |alphatJayatillekeWallFunction|, |nutkWallFunction| with value $0$, |kqRWallFunction| with value \num{4e-5} and |epsilonWallFunction| with value \num{2e-8} and |omegaWallFunction| with value \num{5e-4} respectively. +The turbulence parameters all have the condition |inletOutlet| for the inlet and the outlet. +This is supposed to be an uniform condition at inflow (at inlet) and a zero gradient condition on the outlet. +The pressure is calculated and the actual pressure boundary conditions are given for the adjusted pressure $p_{rgh}$. +All other values have cyclic conditions on the top and bottom +(|type cyclic;|) +and empty conditions for front and back to ensure a 2D-problem. -\section{Solvers} - -\section{Schemes} +The other boundary conditions are listed in table \ref{fig:boundary}. +\begin{table}[ht] + \centering + \caption{The boundary conditions for all value and boundaries except top, bottom, front and back.} + \label{tab:boundaryconditions} + \begin{tabular}{ccccc} + Value & std. value & inlet & outlet & cylinder \\ + $α_t$ & 0 & inletOutlet & inletOutlet & |alphatJayatillekeWallFunction| \\ + $ν_t$ & 0 & inletOutlet & inletOutlet & |nutkWallFunction| \\ + $k$ & \num{4e-5} & inletOutlet & inletOutlet & |kqRWallFunction| \\ + $ε$ & \num{2e-8} & inletOutlet & inletOutlet & |epsilonWallFunction| \\ + $ω$ & \num{5e-4} & inletOutlet & inletOutlet & |omegaWallFunction| \\ + $T$ & \num{300} & Dirichlet $T_{∞}$ & $∇T = 0$ & Dirichlet $T_w$ \\ + $U$ & & Dirichlet $U_{in}$ & |fluxCorrectedVelocity| & |noSlip| \\ + $p_{rgh}$ & & |fixedFluxPressure| & Dirichlet $0$ & $∇p_{rgh} = 0$ + \end{tabular} +\end{table} -% do not use \input nor \begin{document nor \end{document} +\section{Solvers} +% check at end +For $p_{rgh}$ preconditioned conjugate gradient (|PCG|) solver with |DIC|-smoothed |GAMG|-preconditioner. +The preset |GAMG| solver with |DICGaussSeidel| preconditioner was used before and the change did not affect the runtime significantly. -% write the text here -% for including a tikz picture, place it in the directory tikz and -% use \tikzinput{tikzpicture} -% for including another tex (text) file, place it in the directory -% tex and use \textinput{textfile} -% for including some other file use \maininput{file/path/relative/to/main/repodir} +For the other values a smooth solver with symmetric Gauss-Seidel smoother was used. +The preset PBiCGStab solver was used before and the runtime was not significantly affected by the change. +\section{Schemes} +For the time derivative a |backward| scheme was used. +Gauss linear was used as a gradient Scheme while |Gauss limitedLinear 1| was used as all divergence Schemes (in the vector version for $U$.) +% I hope the other ones don't matter +The wall distance (|wallDist|) is calculated via a |meshWave|. \docEnd