Many multiphysics problem can be described by the coupling of several models through physical surfaces. Relying on existing model-specific solvers is very desirable, however they must be coupled in a way that ensures an accurate and stable coupled simulation. In this contribution, we present a multistep coupling scheme which relies on the history of the exchanged quantities to enable a high-order accurate coupling with time adaptation. Explicit and implicit variants are discussed in details. Numerical experiments conducted with an opensource demonstrator on a conjugate heat transfer problem show that high-order convergence is attained, and that stability is favourable compared to other classical approaches.