Cold 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 6.4 km.
Integration time: t=1800s
Profile: Temperature

\(\theta = 3000\,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}=4.0\,km ; z_{c}=3.0 \,km\) and \(z_{r}=2.0 \,km\)

results

top plane: density
upper middle plane: horizontal veocity
lower middle plane: vertical velocity
bottom plane: temperature perturbation