Spin-S bilayer Heisenberg models: Mean-field arguments and numerical calculations