An iterative coupling of a finite-difference boundary layer code "STAN5"[1] and a boundary element code "BEASY”[2] have been used to study two dimensional convective-conductive problems. The coupling method essentially leaves undisturbed the codes of these programs and uses their custom features in their applicable domains. The boundary layer program is used in the convective domain, while "BEASY” is used in the solid domain. Compatibility conditions are established between these domain solutions and an iterative scheme is used until these conditions are mutually satisfied. The indirect communication of compatibility (boundary) criteria between these programs was achieved by two editor/data-transferor like programs. Combined entire package controlled by systems dependent communication program. Numerical scheme used by boundary layer program incorporates a stepwise method. The method presents certain advantages, which makes it suitable for coupling with any other domain method such as the boundary element method. A number of applications have been studied to verify the results for each program as well as the coupling scheme. These include flows over Cylinder inclining surfaces, airfoil and wedge shaped models with internal regions subjected to convective or isothermal environments. The results of the conjugate coupling showed the following: 1) Total job-time reduction due to high convergence rate, and time required for the model preparation. 2) A data storage reduction. 3) A substantial collection of system characteristics. The result of iterative coupling scheme is a useful engineering tool that can be used to study conjugate system consisting of boundary layer and conductive domains of arbitrary shapes. [1] Program and method originally developed by Patankar/Spalding later modifications done by Crawford/Kays and P.R. Iyanger [2] Boundary Element Analysis SYstem developed by Computational Mechanics