Abstract: Numerical simulations of boundary layers play a significant role in the study and interpretation of physical experiments for theoretical explanations of boundary layer disturbances. The influence of thermal boundary layer on the control of heat transport across flat plates is particularly examined. The Crank-Nicolson differential method, which is widely favoured for finite-difference modelling of boundary layer equations, is reviewed. The stability of this method is compared with other numerical approaches in order to establish the appropriate scheme for sustainable applications, involving the design of any conjugate system with heat transfer between the solid and fluid interface. Specific applications to the analysis of cabin comfort in automobiles are anticipated.