A fully-implicit high-order system thermal-hydraulics model for advanced non-LWR. Applied fully-coupled JFNK solution scheme to avoid the operator-splitting errors. The effects of different spatial and temporal discretization schemes are investigated. High-order convergence rates of discretization schemes are verified.