The Welander problem was investigated using high-order numerical methods. High-order numerical methods reduced numerical diffusion significantly for stability analysis. High-order numerical methods were proved to be promising for natural circulation analysis.