El objetivo de este trabajo es describir una implementación eficiente de modelos de gases reales cúbicos y multiparamétricos en un solucionador de Galerkin discontinuo existente para extender sus capacidades a la simulación de flujos turbulentos de gases reales. Los modelos termodinámicos adoptados son van der Waals, Peng-Robinson y Span-Wagner, que difieren entre sí en términos de precisión y costo computacional. Los flujos numéricos convectivos a través de las interfaces de los elementos se calculan con un solucionador de Riemann linealizado consistente termodinámicamente, mientras que para las condiciones de contorno se emplea una expresión linealizada de los invariantes generalizados de Riemann. Las propiedades de transporte se tratan como cantidades dependientes de la temperatura y la densidad a través de correlaciones multiparamétricas. Se adopta una integración temporal implícita; la matriz Jacobiana y las derivadas termodinámicas se obtienen con la herramienta de diferenciación automática Tapenade. La precisión del solucionador se evalúa calculando tanto casos de prueba de gases reales en estado estacionario como no estacionario disponibles en la literatura, y se investiga el efecto del tamaño de la malla y el grado polinómico de aproximación en la precisión de la solución. Se observa un buen acuerdo con datos experimentales y de referencia numérica, y fenómenos específicos no clásicos son bien reproducidos por el solucionador.