source: palm/trunk/TUTORIAL/SOURCE/exercise_topography.tex @ 1538

Last change on this file since 1538 was 1515, checked in by boeske, 9 years ago

several updates in the tutorial

  • Property svn:keywords set to Id
File size: 11.2 KB
Line 
1%$Id: exercise_topography.tex 1515 2015-01-02 11:35:51Z hoffmann $
2\input{header_tmp.tex}
3%\input{../header_lectures.tex}
4
5\usepackage[utf8]{inputenc}
6\usepackage{ngerman}
7\usepackage{pgf}
8\usepackage{subfigure}
9\usepackage{units}
10\usepackage{multimedia}
11\usepackage{hyperref}
12\newcommand{\event}[1]{\newcommand{\eventname}{#1}}
13\usepackage{xmpmulti}
14\usepackage{tikz}
15\usetikzlibrary{shapes,arrows,positioning,decorations.pathreplacing}
16\def\Tiny{\fontsize{4pt}{4pt}\selectfont}
17
18%---------- neue Pakete
19\usepackage{amsmath}
20\usepackage{amssymb}
21\usepackage{multicol}
22\usepackage{pdfcomment}
23\usepackage{xcolor}
24
25\institute{Institute of Meteorology and Climatology, Leibniz UniversitÀt Hannover}
26\selectlanguage{english}
27\date{last update: \today}
28\event{PALM Seminar}
29\setbeamertemplate{navigation symbols}{}
30\setbeamersize{text margin left=.5cm,text margin right=.2cm}
31\setbeamertemplate{footline}
32  {%
33    \begin{beamercolorbox}[rightskip=-0.1cm]&
34     {\includegraphics[height=0.65cm]{imuk_logo.pdf}\hfill \includegraphics[height=0.65cm]{luh_logo.pdf}}
35    \end{beamercolorbox}
36    \begin{beamercolorbox}[ht=2.5ex,dp=1.125ex,%
37      leftskip=.3cm,rightskip=0.3cm plus1fil]{title in head/foot}%
38      {\leavevmode{\usebeamerfont{author in head/foot}\insertshortauthor} \hfill \eventname \hfill \insertframenumber \; / \inserttotalframenumber}%
39    \end{beamercolorbox}%
40%    \begin{beamercolorbox}[colsep=1.5pt]{lower separation line foot}%
41%    \end{beamercolorbox}
42  }%\logo{\includegraphics[width=0.3\textwidth]{luhimuk_logo.eps}}
43
44\title[Exercise - Topography]{Exercise - Topography}
45\author{PALM group}
46
47% Notes:
48% jede subsection bekommt einen punkt im menu (vertikal ausgerichtet.
49% jeder frame in einer subsection bekommt einen punkt (horizontal ausgerichtet)
50\begin{document}
51% Folie 1
52\begin{frame}
53\titlepage
54\end{frame}
55
56\section{Exercise}
57\subsection{Exercise}
58
59% Folie 2
60\begin{frame}
61   \frametitle{Exercise}
62   Please carry out \textbf{two runs} with following conditions.
63   \begin{itemize}
64      \item<2->{Single cube}
65      \begin{itemize}
66         \item[1.)]{First run ''generic'' using {\tt topography = 'single\_building'}}
67         \item[2.)]{Second run ''raster'' using {\tt topography = 'read\_from\_file'}
68                   with ASCII file ...\_topo}
69      \end{itemize}
70      \item<3->{Neutral boundary layer in a channel}
71      \item<4->{Constant bulk velocity}
72      \item<5->{No Coriolis force}
73      \item<6->{Simulation features:}
74      \begin{itemize}
75         \item{domain size: (80 m)$^3$ (x/y/z)}
76         \item{grid size: 2 m equidistant}
77         \item{cube: size (40 m)$^3$, location centered in the domain center}
78         \item{simulated time: 7200 s}
79         \item{initial velocity: u = 1, v = 0 m/s}
80      \end{itemize}
81   \end{itemize}
82   \onslide<7->\textbf{Please use the same building (size, location) for both runs!}
83\end{frame}
84
85% Folie 3
86\begin{frame}
87   \frametitle{Questions to be Answered}
88   \small
89   \begin{itemize}
90      \item<2->{Can you identify flow convergence / divergence patterns near the cube?}
91      \begin{itemize}
92         \item{What kind of output do you need to answer this?}
93      \end{itemize}
94      \item<3->{How does the horizontally and temporally averaged momentum flux profile look
95            like?}
96      \begin{itemize}
97         \item{How long should the averaging time interval be?}
98      \end{itemize}
99      \item<4->{Is it really a large-eddy simulation?}
100      \begin{itemize}
101         \item{Are the subgrid-scale fluxes much smaller than the resolved-scale fluxes?}
102         \item{How do the total kinetic energy and the maximum velocity components change
103               with time?}
104      \end{itemize}
105   \end{itemize}
106   \onslide<5->\textbf{Final question:}
107   \begin{itemize}
108      \item{Do the results of both runs agree?}
109   \end{itemize}
110\end{frame}
111
112% Folie 4
113\begin{frame}
114   \frametitle{Hints (I)}
115   \scriptsize
116   \begin{itemize}
117      \item<2->{\textbf{Domain size}}
118      \begin{itemize}
119         \scriptsize
120         \item{Is controlled by grid size (\textbf{dx}, \textbf{dy}, \textbf{dz}) and number of grid points
121               (\textbf{nx}, \textbf{ny}, \textbf{nz}). Since the first grid point along one of the directions has
122               index 0, the total number of grid points used are \textbf{nx}+1, \textbf{ny}+1, \textbf{nz}+1.
123               The total domain size in case of cyclic horizontal boundary conditions is
124               (\textbf{nx}+1)$\cdot$\textbf{dx}, (\textbf{ny}+1)$\cdot$\textbf{dy}.}
125         \end{itemize}
126      \item<3->{\textbf{Initial profiles}}
127      \begin{itemize}
128         \scriptsize
129         \item{Constant with height. See parameter \textbf{initializing\_actions} for available
130               initialization methods. See \textbf{ug\_surface}, \textbf{vg\_surface} for initial values of
131               velocity.}
132      \end{itemize}
133      \item<4->{\textbf{Boundary conditions}}
134      \begin{itemize}
135         \scriptsize
136         \item{For channel boundary condition, see \textbf{bc\_uv\_t}.}
137      \end{itemize}
138      \item<5->{\textbf{Forcing}}
139      \begin{itemize}
140         \scriptsize
141         \item{For constant bulk velocity, see \textbf{conserve\_volume\_flow}.}
142         \item{For Coriolis force, see \textbf{omega}.}
143      \end{itemize}
144      \item<6->{\textbf{Topography}}
145      \begin{itemize}
146         \scriptsize
147         \item{For generic topography, see \textbf{building\_height}, \textbf{building\_length\_x} and
148               \textbf{building\_length\_y}.}
149         \item{For raster topography, please use a text editor to manually create an
150               ASCII ''raster\_topo'' file that contains the same building.}
151      \end{itemize}     
152    \end{itemize}
153 \end{frame}
154
155% Folie 5
156\begin{frame}
157   \frametitle{Hints (II)}
158   \footnotesize
159   \begin{itemize}
160      \item<2->{\textbf{Simulation time}}
161      \begin{itemize} 
162         \footnotesize       
163         \item{See parameter \textbf{end\_time}.}
164      \end{itemize}
165      \item<3->{\textbf{Variables}}
166      \begin{itemize}
167         \footnotesize
168         \item{Output variables are chosen with parameters \textbf{data\_output} (3d-data or
169               2d-cross-sections) and \textbf{data\_output\_pr} (profiles).}
170         \item{Time series are activated using \textbf{dt\_dots}.}         
171      \end{itemize}
172      \item<4->{\textbf{Output intervals}}
173      \begin{itemize}
174         \footnotesize
175         \item{Output intervals are set with parameter \textbf{dt\_data\_output}. This parameter
176               affects all output (cross-sections, profiles, etc.). Individual temporal
177               intervals for the different output quantities can be assigned using
178               parameters \textbf{dt\_do3d}, \textbf{dt\_do2d\_xy}, \textbf{dt\_do2d\_xz}, \textbf{dt\_do2d\_yz}, \textbf{dt\_dopr},
179               etc. }
180      \end{itemize}
181      \item<5->{\textbf{Time averaging}}
182      \begin{itemize}
183         \footnotesize
184         \item{Time averaging is controlled with parameters \textbf{averaging\_interval},
185               \textbf{averaging\_interval\_pr}, \textbf{dt\_averaging\_input}, \textbf{dt\_averaging\_input\_pr}.}
186      \end{itemize}
187   \end{itemize}
188\end{frame}
189
190% Folie 6
191\begin{frame}
192   \frametitle{Further Hints}
193   \scriptsize
194   Please see under \\
195   \textbf{http://palm.muk.uni-hannover.de/wiki/doc/app/netcdf}  \\ 
196   \par\medskip 
197   where the complete PALM netCDF-data-output and the respective steering parameters are described.
198   \par\medskip
199   For topography, see \\
200   \textbf{http://palm.muk.uni-hannover.de/wiki/doc/app/inipar\#topo}\\
201   \par\medskip
202   and especially for raster topography, see also
203   \textbf{http://palm.muk.uni-hannover.de/wiki/doc/app/iofiles\#TOPOGRAPHY\_DATA} \\
204   \par\medskip
205   as well as the presentation ''Using topography (I)''.
206\end{frame}
207
208% Folie 7
209\begin{frame}
210   \frametitle{Proceeding}
211   Please proceed as follows:
212   \begin{itemize}
213      \item<2->[1.]{Please run with the ''generic'' topography case first.}
214      \item<3->[2.]{Check your results to answer all questions – except the final question.}
215      \item<4->[3.]{After this run has finished, use ncview, ncdump etc. to check the precise
216               location of the building (look at 2D array zusi that is contained in 2D xy
217               cross-sections and 3D volume data).}
218      \item<5->[4.]{Use this information to manually create the ''raster\_topo'' file.}
219      \item<6->[5.]{Run the ''raster'' topography case.}
220      \item<7->[6.]{Compare both simulation results to answer the final question.}
221   \end{itemize} 
222\end{frame}
223
224% Folie 8
225\begin{frame}
226   \frametitle{How to Start?}
227   \footnotesize
228   \begin{itemize}
229      \item<2->{Create two \textbf{INPUT} directories for both new runs: \\
230            {\tt cd $\sim$/palm/current\_version} \\           
231            {\tt mkdir -p JOBS/generic/INPUT} \\
232            {\tt mkdir -p JOBS/raster/INPUT}} 
233      \item<3->{Create the parameter files and {\tt raster\_topo} file and set the required
234            parameters in \\
235            {\tt JOBS/generic/INPUT/generic\_p3d} \\
236            {\tt JOBS/raster/INPUT/raster\_p3d}}
237      \item<4->{Start the runs one by one with mrun-commands \\
238            {\tt mrun -d generic -K parallel ...} \\
239            {\tt mrun -d raster -K parallel ...}}
240      \item<5->{and analyze the output files in \\
241            {\tt JOBS/generic/OUTPUT} \\
242            {\tt JOBS/raster/OUTPUT}}
243   \end{itemize}
244\end{frame}
245
246\section{Results}
247\subsection{Results}
248
249% Folie 9
250\begin{frame}
251   \frametitle{Flow convergence / divergence (I)}
252   \includegraphics[width=0.45\textwidth]{exercise_topography_figures/cross_sections/u_xy.eps} \hspace{0.8cm}
253   \includegraphics[width=0.45\textwidth]{exercise_topography_figures/cross_sections/v_xy.eps} 
254\end{frame}
255
256% Folie 10
257\begin{frame}
258   \frametitle{Flow convergence / divergence (II)}
259   \includegraphics[width=0.45\textwidth]{exercise_topography_figures/cross_sections/w_xy.eps} \hspace{0.8cm}
260   \includegraphics[width=0.45\textwidth]{exercise_topography_figures/cross_sections/w_xz.eps} 
261\end{frame}
262
263% Folie 11
264\begin{frame}
265   \frametitle{Streamlines}
266   \includegraphics[width=0.45\textwidth]{exercise_topography_figures/streamlines/streamlines_xy.eps} \hspace{0.8cm}
267   \includegraphics[width=0.45\textwidth]{exercise_topography_figures/streamlines/streamlines_xz.eps} \hspace{0.8cm}
268\end{frame}
269
270
271% Folie 12
272\begin{frame}
273   \frametitle{Vertical profiles of $\overline{w'u'}$, $\overline{w'v'}$}
274   \includegraphics[width=0.45\textwidth]{exercise_topography_figures/profiles/wu_time_pr.eps} \hspace{0.8cm}
275   \includegraphics[width=0.45\textwidth]{exercise_topography_figures/profiles/wv_time_pr.eps}
276\end{frame}
277
278% Folie 13
279\begin{frame}
280   \frametitle{LES? - Fluxes}
281   \begin{center}
282      \includegraphics[width=0.6\textwidth]{exercise_topography_figures/profiles/wu_comp_pr.eps}
283   \end{center}   
284\end{frame}
285
286% Folie 14
287\begin{frame}
288   \frametitle{LES? - Time Series (I)}
289   \begin{center}
290      \includegraphics[width=0.95\textwidth]{exercise_topography_figures/timeseries/E_ts.eps} \\
291      \includegraphics[width=0.95\textwidth]{exercise_topography_figures/timeseries/umax_ts.eps} 
292   \end{center}
293\end{frame}
294
295% Folie 15
296\begin{frame}
297   \frametitle{LES? - Time Series (II)}
298   \begin{center}
299      \includegraphics[width=\textwidth]{exercise_topography_figures/timeseries/vmax_ts.eps} \\
300      \includegraphics[width=\textwidth]{exercise_topography_figures/timeseries/wmax_ts.eps}
301   \end{center}
302\end{frame}
303
304\end{document}
Note: See TracBrowser for help on using the repository browser.