`_).
The differences to the original problem are that the grid is not adaptive and no
stabilisation method is used.
The problem can be described using the following equations (the Boussinesq equations):
.. code-block:: none
-div(2 * eta * eps(u)) + nabla(p) = -rho * beta * g * T in Omega
-div(u) = 0 in Omega
dT/dt + div(u * T) - div(k * nabla(T)) = gamma in Omega
The mesh is a simple square (0,1)x(0,1):
.. image:: _static/square(0,1)x(0,1)-50x50.png
:width: 300 px
The temperature and the velocity vectors plot:
.. image:: _static/tutorial_dealii_7-results.png
:height: 400 px
Animation:
.. image:: _static/tutorial_dealii_7-animation.gif
:height: 400 px
]]>