Towards a numerical benchmark for 3D mixed-convection low Mach number flows in a rectangular channel heated from below