The energy equilibrium equation and discrete ordinate methods are combined to establish the one-dimensional steady heat transfer mathematical model of multi-layer thermal insulations (MTIs) in metallic thermal protection systems. The inverse problem of heat transfer is solved by the genetic algorithm and data from the steady heat transfer experiment of fibrous thermal insulations. The density radiation attenuation coefficient, the albedo of fibrous thermal insulations and the surface emissivity of reflective screens are optimized. Finally, the one-dimensional steady heat transfer model of MTIs with optimized thermal physical parameters is verified by experimental data of the effective MTI conductivity.