A robust high-order compact method for large eddy simulation