This paper presents finite-element solution of the viscous compressible flow governed by Navier-Stokes equations in two dimensions. The numerical algorithm is two step Taylor-Galerkin method using unstructured triangular grids. To increase accuracy and stability, combined moving node method and grid refinement method have been used.
To validate the present algorithm, computational results are compared with the existing experimental data and other numerical solutions. They showed good agreement one another. Four test problems are solved for demonstration of this numerical approach. They include subsonic flow over flat plate, Carter flat plate problem, laminar shock boundary layer interaction, and finally laminar flow around a NACA0012 airfoil with zero angle of attack and free stream Mach number 0.85. The results have indicated that adaptive triangular grid is convenient for present preliminary viscous flow calculation but needs more research for complicated geometries.