Univariate polynomial equation providing on-lattice higher-order models of thermal lattice Boltzmann theory.

: A univariate polynomial equation is presented. It provides on-lattice higher-order models of the thermal lattice Boltzmann equation. The models can be accurate up to any required level and can be applied to regular lattices, which allow efficient and accurate approximate solutions of the Boltzmann equation. We derive models approaching the complete Galilean invariant and providing accuracy of the fourth-order moment and beyond. We simulate one-dimensional thermal shock tube problems to illustrate the accuracy of our models. Moreover, we show the remarkably enhanced stability obtained by our models and our discretized equilibrium distributions.