Using the problem symmetry, only one half of the domain was
considered in the model. Besides symmetry boundary conditions,
the right edge of the domain was let free and a vertical displacement
was applied on the top boundary corresponding to the symmetry
plane: