目前主流的热工水力系统分析程序采用一阶差分格式来离散基本守恒方程。一阶差分截断误差大,导致系统分析程序在模拟自然循环问题时计算精度不高。为提高计算精度,本文采用两流体双压力两相流模型高精度数值解法对Welander自然循环问题进行计算分析。数值结果表明,高精度数值解法能有效减小数值扩散,提高预测精度。