Abstract A fully coupled high order discontinuous Galerkin (DG) solver for viscoelastic Oldroyd B fluid flow problems is presented. Contrary to known methods combining DG the discretization of convective terms material model with standard finite element (FEM) and using elastic viscous stress splitting (EVSS) its derivatives, a local (LDG) formulation first described hyperbolic convection‐diffus...