The topic sound interesting and your first step promising. For the coupling you could look in the OpenLB example section “thermal”. There, a coupled model is used, ADE + NSE with a coupling approach based on Boussinesq-approximation. It should go similar to it.. You can keep me up to date on it by email, it would be good to share your results with then OpenLB community..
Best
Mathias