The density current test case
Geometry:
Computational domain extends in horizontal direction from -25.6 to 25.6 km and in vertical direction from 0 to 12.8 km.
Profile:
Temperature
\(\theta = 300.0 \,K\)
\( \Delta T= \begin{cases}
0.0 \text{if}\quad L \gt 1.0, \\
-15.0 \cos(L \pi/2) \text{if}\quad L \le 1.0
\end{cases}\)
where
\(L= ((\frac{x-x_{c}}{x_{r}})^{2} + (\frac{z-z_{c}}{z_{r}})^{2})^{1/2}\)
and
\(x_{c}=0.0 \,km ; x_{r}=2.0 \,km ; z_{c}=2.0 \,km\) an \(z_{r}=2.0\,km\)
results
top plane: | density |
upper middle plane: | horizontal veocity |
lower middle plane: | vertical velocity |
bottom plane: | temperature perturbation |