In this paper we present computations of the nonweak solutions of class C11 of the strong form of transient Navier–Stokes equations for compressible flow in Lagrangian frame of reference using space–time least squares finite element formulation (STLSFEF) with primitive variables ρ, u, T. For high speed compressible flows the solutions reported here possess the same orders of continuity as the governing differential equations (GDEs). It is demonstrated that with this approach accurate numerical solutions of Navier–Stokes equations without any assumptions or approximations are possible. In the approach presented here SUPG, SUPG/DC, or SUPG/DC/LS operators are neither used nor needed. The role of diffusion, i.e., viscosity (physical or artificial) and thermal conductivity on shock structure is demonstrated. Merits of using ρ, u, T as primitive variables over ρ, u, p are discussed. Compression of air in a rigid cylinder by a rigid, massless and frictionless piston is used as a model problem. True time evolutions of class C11 are reported beginning with the first time step until steady shock conditions are achieved. Comparisons with analytical solutions are presented when possible.