  In this work, we develop adaptive schemes using goal-oriented error control
for a highly nonlinear flow temperature model with temperature dependent
density. The dual-weighted residual method for computing error indicators to
steer mesh refinement and solver control is employed. The error indicators are
used to employ adaptive algorithms, which are substantiated with several
numerical tests. Therein, error reductions and effectivity indices are
consulted to establish the robustness and efficiency of our framework.
