Abstract

We develop a virtual element method to solve a convective Brinkman-Forchheimer problem coupled with a heat equation. This coupled model may allow for thermal diffusion and viscosity as a function of temperature. Under standard discretization assumptions and appropriate assumptions on the data, we prove the well posedness of the proposed numerical scheme. We also derive optimal error estimates under suitable regularity assumptions for the solution and appropriate assumptions on the data. We conclude with a series of numerical tests performed on different families of meshes that complement the theoretical findings. © 2025 Elsevier Ltd