A multi-region algorithm for N-TH coupling calculation and its application to nuclear thermal propulsion reactor