The present study aims to propose a coupled numerical model for analyzing the influence of the configuration of the vapor-cooled shield tubes on insulation performance. The numerical model is developed to calculate the effective thermal conductivity of multi-layer insulation (MLI) and rigorously evaluated. Specifically, the tank evaporation simulation software BoilFAST is manually coupled with the CFD commercial software of ANSYS FLUENT (2020 R2) to accurately determine the mass flow rate and temperature of boil-off gas in a tank, which is then used as inputs for the CFD simulations. This coupling methodology can be used for reliable calculations of boil-off losses for liquid hydrogen and heat transfer inside the MLI structures. Results demonstrate that the estimated effective thermal conductivity agrees with the earlier published data and that increasing the diameter and number of tubes enhances heat exchange between the VCS tubes and MLI structures. Also, the heat flux into the tank is relatively reduced by 32.04 % by increasing the diameter of the tube and 49.14 % by adding the number of tubes. Moreover, the maximum temperature difference of the VCS is estimated to be less than 1 K, indicating nearly uniform distributions of the VCS.