Numerical simulations of thermally stratified flows in an open domain