Page 215 - Computational Fluid Dynamics for Engineers

6.6 Applications of the Panel Method 203

remains unchanged, subroutine POISSON (modified subroutine POISSON) has
some minor changes as discussed below.
The statements in the beginning of this subroutine starting from

ELEMENTS OF DIAGONAL VECTORS IN THE BLOCK MATRICES A, B AND C

require changes, as indicated below, up to the statement
ELEMENTS OF W-VECTOR, SEE EQ. (4.5.2)

The remaining statements stay unchanged. The listing is given in Appendix B,
Chapter 6.

Appendix 6B Panel Program for an Airfoil

6B.1 M A I N Program

In the revised MAIN program given below, we read the value of r (= GAMMA),
the total number of grid points (NODTOT) in the interval 0 < 6 < 2TT. The
listing is given in Appendix B, Chapter 6.

6B.2 Subroutine COEF

Since the Kutta condition is not used, the statements regarding this condition
are removed as shown in Appendix B, Chapter 6.
The index KUTTA, which represents the total number of equations, is re-
placed by NEQNS = NODTOT, and Eq. (6.5.2) is incorporated in the statement
under "ELEMENTS OF VECTOR B for I = 1 to N."

6B.3 Subroutine VPDIS

With values of GAMMA given as the input and passed from the MAIN program,
the revised subroutine is given in Appendix B, Chapter 6.

Appendix 6C Panel Program for Multielement Airfoils

6C.1 M A I N Program

The number of elements (NBODY) the number of panels for the i^-th element
(NODT(K)), airfoil coordinates (X(I,K), Y(I,K)) for the K-th element, the an-
gle of attack (ALPHA) and the reference length (CHORD) are input to the
revised MAIN Program. The listing of MAIN and other subroutines are given
in Appendix B, Chapter 6.