De Vahl Davis natural convection test

The classical de Vahl Davis benchmark test is defined for the natural convection of the air ($\Pr =0.71$) in the square closed cavity. The only physical free parameter of the test is the thermal Rayleigh number. In the original paper [de Vahl Davis, 1983] de Vahl Davis tested the problem up to the Rayleigh number , however in the latter publications, the results of more intense simulations were presented with the Rayleigh number up to . Lage and Bejan [Lage and Bejan, 1991] showed that the laminar domain of the closed cavity natural convection problem is roughly below $\text{Gr1}{{\text{0}}^{9}}$. It was reported [Janssen and Henkes, 1993; Nobile, 1996] that the natural convection becomes unsteady for $\text{Ra}=2\cdot {{10}^{8}}$. Here we present a MLSM solution of the de Vahl Davis case.