A high-order discontinuous Galerkin method for simulating incompressible fluid-thermal-structural interaction problems

•An accurate, high-order viscous discontinuous Galerkin solver for the computation of fluid-thermal-structure interaction of Boussinesq flows is proposed.•Formal h and p convergence studies employing the Method of Manufactured Solutions verify that the expected order of accuracy (up to sixth) is achieved.•Four benchmark cases further validate the new scheme.