Abstract A class of Bernstein-Bézier basis based high-order finite element methods is developed for the Galerkin-characteristics solution convection-diffusion problems. The formulation derived using a semi-Lagrangian discretization total derivative in considered spatial performed method on unstructured meshes. Lagrangian interpretation this approach greatly reduces time truncation errors Euleri...