A computer program for the finite element calculation of axisymmetric subcritical compressible flow

A computer program, developed at NLR, for the calculation of the axisymmetric subcritical compressible flow around and inside bodies of revolution is described. The solution is formulated in terms of a velocity potential which is approximated quadratically using TRIM6 and isoparametric elements for the discretization of the flow field. The program is written in ALGOL for a CDC Cyber 72 computer. A survey is given of the formulae. The flow of the computation is illustrated by means of a flow chart. The procedures are described and a complete list of input and output quantities is included. The computer program was applied to different configurations. A typical example involved 902 unknowns to be determined. The semi-bandwidth was 69 and the number of iteration steps required was 6.