Heat Bubble

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