Extension of PPVN.f to Axisymmetric Problems

A pre-processor program should usually be used to generate the following data lines 

As already mentioned, the present code corresponds to the solution of steady-state non-isothennal Navier-Stokes equations in two-dimensional Cartesian domains by the continuous penalty method. As an example, we consider modifications required to extend the program to the solution of creeping (Stokes) non-isothermal flow in axisymmetric domains  

Step 1 To solve a Stokes flow problem by this program the inertia term in the elemental stiffness matrix should be eliminated. Multiplication of the density variable by zero enforces this conversion (this variable is identified in the program listing). 

Step 2 General structure of stiffness matrices derived for the model equations of Stokes flow in (x, 3O and (r, z) formulations (see Chapter 4) are compared. 

216 COMPUTER SIMULATIONS - FINITE ELEMENT PROGRAM Stiffness matrix corresponding to flow equations in (x, 3O formulation 

---