A stabilized weighted interior penalty method for thermal convection model in heterogeneous porous media | Synapse