Parallel solver for unsteady Rayleigh Benard problem in limited domains