1980020753
1980020753
1980020753
C
APPENDIX
38 IF(H)39,39,40
39 HE -HE
GO TO 42
40 UE UK USTR*OL
SUHPXOL
CALL UMP(UEpOZEPHE_FUME)
IF(FUME-EUM)42p42,41
41 DL = OL$FUMIIFUN-FUNE)
GO TO 13
4Z DCQ DCO + OL*VN
X-X+DL
IF(X -OLV)47,50pSO
47 H HE
O2 DIE
GO TO(43,13,43p45,13p43)pISGV "
43 IF(HS-.l*ESP)44p13,13
44 OL Z**OL
ISTAB =ISTAB-1
GO TO 13
45 IF(HOIFF)4%,ZSp13
49 STP II(ISTAS-q)
50 HR HE
DZR OZE
RETURN
ENO
FUNCTION RUNO(A_8)
RUNO (AINT(A*Be$IGN(.SpA)))IB
RETURN
ENO
FUNCTION COSG(A)
COSG COS(A**O174532925199)
RETURN
END
FUNCTION SING(A)
SING - SIN(A**O17_532925199)
RETURN
END
FUNCTION CSLG(ApB}
CSLG ALOGIABS(SING(A-B)))
RETURN
ENO
SUBROUTINE ORAW(WCpWSjWLpDKeO_jFL_AGpRA)
CALL DRAW(WC_WSeWLpORAK(J)eORAR(J)eFLA(J)PABGRp_KP)
WHEN RKP = Oe DK - KAPPA; OTHERWISE, DK K.
COSP COSG(AG*FL)
FK DK
IF(RA)Z,1P2
1 FK FK*(1.-COSP)I(I=COSP)
2 SINP $1NG(AG*FL)
113
APPENDIX
WL = -DMtALOG(1.eFK)
IF(FLoEO.O.) FK 1.
BETA -(COSP-1.)IFK + COSP
BQMI - BETA$*2 - I.
WUBEO= $ORT(ABS(BOMI))
U= (I.+BETA}$SINPI(I.COSP)
IF(BOMZIS, Sp4
3 WF ALOG(ABS[(WUBEOeU)I(WUBEO-U)))
GD TD 6
4 WF Z.* ATAN(UlWUBEO)
GO TO 6
WF Go
WC -(WUBEO*WF - SINP - BETAeFL*AGeI.74_329E-Z)eDM
W$ (COSP-I.)e(1.-(lolFK 1.)*ALOG(1.eFK))*DM
RETURN
ENO
SUBROUTINE TRAPRO
DIMEHSION FLS(Z)pFLA(E)pDRAK(Z),DRAM(Z)_AC(4p3)PD(3)p
I_$I(2)pWCI(Z)_FINT(3)pA(4)jHK(Z)PR(3)tFKERN(30)
CDMMON/EA/ILESpIDRUplSTApNNESE
CO_MDNIPRAL/DLTsDLTUpALNpALV(14)pNALPITP_NAMP(12)pCML(14)pCRL(Z4)
I _CPV[Z)pALTXt4,E)pDARFJITITZpITIT2
COMMON pl(1Z1)pP(1ZI);XP(1ZZ)wYP[lZ1)_PUFF(14)*AGAM(|4)jX(122)_
1y(_z1),DS(1ZZ)_VF(lZ1)pARG(1Z1),ANI(ZB)tALFR(Zg)_IZZIKFU_NQ_NUPRO'
2JAB_JSTsCMpETA_ABFApPI_BOGENpDARGpPURES(ZS)pGAP(450)_ALFA(29)
DATA ABSZIO.I
CALL WANDEL(NUPROpNAHPp12,_)
ALFR(JAB+_)=O.
ABZT'_NI(JAB)
IF[ABS[ABZT-ABSZi.LT**I) GO TO 14
IB_INT(.ZS_ABZT+.I)
_Oo2*IB
NKRm2*M_
ABSZ-FLOAT(NKR)
ABGRs360olABSZ
HABGR=SeABGR
DO B M=I, IB
ARI=FLOAT(MQ+I-Z*MI*NABGR
B FKERN(M].ABGR*COSGIARI)I(SING(ARI)*PI)
14 MAGAM=INT(AGAM(3))
NQ_NXR_I
IF(MAGAM.EQ.O) GO TO 22
NZT-NZPZ(3_O)
WRITE(IDRU_SEINZT
Z2 DO 23 I'1,2q
Z3 ALFA(1)=ALFR(I)
1-1
Jl
Z4 FLS(J)- PURES(I)eABFA
114
,
APPENDIX
CALL
CALL
WCI(J)- WCI(J)+WC
WSI(J)- WSI(J)+WS
WLI WLIWL
FLA(J)- PURE$(I+I) * ABFA
IF(FLA(J))ZSpZS_Zb
Z5 DRAK(J)I 0
ORAM(J) 1.
GOTO 34
26 WI - COSG(ABGR*FLA(J))
IF(PURES(I+Z)-I.)ZTs3OpZ9
27 ORAK(J) .I*PURE$(I+3)
Z80RAM(J)- .I*PURE$(I+4)
GOTO 34
29
DRAW(WCpWSpWL,.b,-1.PFL$(J)gABGRal)
DRAW(WCI(J)_W$I(J)pWLIp-,bp-I.pFLS(J),ABGR,1)
DRAK(J)'((.I*PURES(I4})**(-IO./PURE$(I3))-I.)*(1.+WI)/(1._WI)
DRAK(J)- RUND(DRAK(J)_IO00.)
DRAM(J) .I*PURES(I3)
GD TO 34
30 AA .05*(I.-WI )*PURE$(I3)
WILN - ALDG(.I*PURE$(I4))
FMIT ._
MIT 0
31 FM -WILN/ALDG(AA/FMIT +1.1
MIT-MITI
IFIABS(FM-FMIT)-I.E-6) 33,3Z_3Z
3Z FMIT FM
GD TO 31
33 DRAMiJ) - RUND(FMplO00.)
ORA$ - -05_PURE$(I3)*(WII,)/FM
DRAK(J) RUNO(ORASplO00.)
34 Io I+5
J" J+l
IF(J-3)24_38p38
38 MER 0
W$1(Z)- -WSI(Z}
O(1) - WLI *(WSI(Z)W$I(I})
D(Z) --WLI *(WCI(Z)+WCI(I))
D(3) - WCI(1)*WSI(Z)-WCI(Z)*W$I{1)
ITMOOeINT(PURE$(11))
RUFIO0.
IF(ITMOD*GE.4.AND.ITMOO.LE.6)RUF-IO00.
ITMRITMOD
SHKS .I*PURES(IZ|
HKST'.l*ABS(PURE${13))
35 DO 36 J'lp3
36 AC(I_J)- O.
ALIV O.
SINAI O.
115
APPENDIX"
COSAI " I.
FNI " O.
J'l
37 C$AIP " COSG(Z._ALFA(J)}
SNAIP - %ING(Z**ALFAIJ))
IF(J-JST-1)40,39p40
39 AC(2,1)" SINAI
AC{Z,Z)" -I.-COSAI
ACi2_3)" -1.
AC(3pI)" -SNAIP
AC(3,2)- 1.CSkIP
AC(3p3)- 1.
AC(4_I)- COSAI-CSAIR
AC(4p2)u $INAI-SNAIP
AC(4_3)- O
ALI5 - ALIV
ALISP - ALFA{J)
GOTO 41
40 FIT - CSLG(HABGR*FNI-qO.JALIV)
FIIP- CSLGfHABGR_FNI-90.pALFA(J)|
PB-FNI*HABGR*BOGEN
ACil, 1).-FIIP*SNAIPFII*$INAI+(COSAI-CSAIP)*PB AC(Z#I)
AC(1, Z}.-FII*(loCOSAI)FIIP*(1._CSAIP)+(EINAI-SNAIP)*PBAC{I_2|
AC{1,3)- FIIP- FII AC{I_3)
41 IF(J-JAB-1)42p43t43
42 ALIV mALFA(J!
$1NAI_SNAIP
COSAI_C$AIP
FNI ANIIJ)
J'Jl
GO TO 37
43 DO 47 J-lp2
IF(FLA(J)! 47p47p49
49 CALL DRAW(WCpWSjWLPDRAK(J)pDRAM(J)_FLA(J)pABGRpO)
AC(I_I)- WC+ AC(1,1)
IFiJ-Z)45,_4p45
_4 W$ J -_$
WL "- WL
45 AC(lJZ) - W$ AC(1_2)
AC(1,3) --WL AC(1,3)
47 CONTINUE
DO 5Z J-1_4
A(J} - O.
DO 52 11_3
52 A(J) - A(J)+Dil)*AC(JpI)
SOLUTION OF TRAN$CENOENTAL EQUATION
53 I'O
FV 9E9
PHISH - ._ *(ALI$ALISP)
116
APPENDIX
6O
6Z
65
CSLI - CSLG(PHISN, ALIS)
CSLIP- CSLG(PHISH, ALISP)
FPuA (1 ]+A( Z)*CSL IA(3 )*CSLI P+A( 4 )*BOGEN* (90. +PHISH)
IF(I.GE.20) GO TO 6&
IF(ABS(FP)-ABS(FV).LT.-.SE-9) GO
I=20
PHISH " PHISH - POIF
GO TO 60
POIF -FP / (A(2)IIPHISH-ALIS)
I = I+l
FV=FP
PHISH PHISH + POIF
TO 62
A(3)/(PHISH-ALISP))
IF(PHISH.LT.ALIS.AND.PHISH.GT.ALISP) GO TO 60
WRITE(IDRUP6)HERJITMQD
66 FORHAT (66HOTRANSCENOENTAL EQUATION HAS OIVERGED. CHECK TRA1ANO
eTRAZ CAROS.,1ZH ITERATIONJIgs8H HODE _I1)
STOP
66 ANI(JST) (PHISH90.)IHABGR
70 DO 71 I=1;3
71 FINT[I)=AC|I_I)+AC|gsI)*CSLI+AC(3_I)*CSLIPAC(_pI)*BOGEN*(PHISH+
FgO.)
69 HK(1)=(FINT(1)*WLI-FINT[3)*WCI{Z))ID(2)
HK(Z)'(FZNT(1)*WLI+FINT(3)*WCI(1))/D(2)
HK$ HK(1)+HK(2)
IF(ITMOD.EQ.O.OR.ABS(HKS-SHKS).LT.HKST] GO TO 76
IF(HAGAH.LT.Z.ANO.HAGAN-HER.NE.1) GO TO 100
GO TO 76
74 ITHO0 =0
IF(HAGAHEQ.O) GO TO 300
76 NZTuNZpZ(ZIJA8+6)
WRITE(IDRU_77]NZT_NUPROpNEReITNR
77 FORHAT (AZs62HTRANSCENOENTAL EQUATION RESULTS AIRFOIL _I6p
*IZH ITERATION_IZ_EH MODE ,I1)
NZT=NZPZ(1;O)
WRITEIIDRU_78]NZT
78 FORMAT (AlpbqH NU ALPHA* OffiEGA t OMEGA K HU K H
* LAHBOA LAHBDA*)
JHm 1
DO 85 JN'lpJAB
79 NZT=NZPZ(E;O)
IF(JN.NE,IAND.JNNEJAB) GO TO 83
Xl .5.(1. COSG(FLA(JHI*ABGR))
WHK (IORAK(JH)*(1.-Xl)FXl)**(-DRAfl(JH))
MSTR ORAHIJHI*DRAK(JH)/X1
MRITE(IDRU_82)NZTpANI(JN)pALFA(JN)pWSTR_WHKeDRAK(JH)_DRAH(JH)_HK(j
1H)_FLA[JH);FLS(JH)
BE FORMAT (AI_Fb.Z_FBZ_FB.3eFT3s3FB3_F8.Z_FT.Z)
JH=Z
GO TO 85
117
APPENDIX
83 WRITE(IDRU_SZ)NZTp&NI(JN)_&LFA(JN)
85 CONTINUE
IF(ITHOD.EQ,O) GO TO 300
100 IF(MER)IO3_IOZeZ03
102 DAL =o1
GO TO 10,
103 IF(HKS-HKSV.EQ.O.)GO TO 74
DAL - ($HK$-HKS)eD&LI(HKS-HKSV)
DALD-D&L
IF(ITP,EQ,O]DAL=RUND(DAL_RUF)
IF(_AGAH.EQ.O0)GO TO 100_
NZT'NZPZ(ZPO)
WRITE(IDRU_IOO3)NZTpMERpHKSpDALOpDAL
1003 FORMAT (AZelOH ITERATION_IZ_3XsSHK S =_Fqobp3XpTHDELTA =pF1Z.Bp
s3X_gHROUNOED -_F6oZ)
1004 ZF(DAL.EQ00,)GO TO 7_
ZF(MERoGEo3oANDoAB$(DALVI.LE.ABS(DAL))GO TO 76
ZOk DALV-DAL
IF(IT_ODoGEo_)GO TO 113
DO 111 J-1pJAB
IF(ITMOD.NE.Z.&MD.JoLE.JST) ALFA(J)=ALFA(J)D&L
IF(ZTMOD.NE.ZoAND.J.GT.JST) ALFA(J):ALFA(J)-D&L
111CONTZNUE
GO TO 11Z
113 IF(ITMOD.GE07]GO TO 114
IF(ITMOD.NEoS)DRAK(I)-DRAK(I)+DkL
IF(ZT_ODoNEo6)DRkKIZ)=DRkK(Z)*DAL
GO TQ lie
11_ ZF(ZTMOD.NEoB)ALFA|JST)-&LF&IJST)DAL
IF(;TMOD.NE.7)ALFA(JST+I)-&LFA(JST+I)-DAL
112 HKSVoHKS
_ER =MER*I
GOTO 35
300 kK- .5*(COSG(PHISH-&LFAiJSTeZ))I$ING(PHISH-ALFA(JST*I))
1 -COSG(PHISH-ALFA(JST))ISZNGiPH|$H-ALFA(JST)))
AKP .AK_180.19.86960_
PHIM = O.
NU=I
1- 1
ANU sO.
JH=O
VI- O.
302 JH'JH+I
FF1 COSG(kBGR*FLA(JH))
FFZ - DRAK(JH)I(Z.FF1)
FG1 o COSG(kBGRtFLS(JH))
FG3 - .bI(FGI-1.)
304 V|= VZ - CSLG(PHZN-90.p ALFA(Z))
GO TO 310
118
306
310
31Z
31_
320
324
326
1
i
328
i
APPENDIX
ARGN = ANU
IF(ANU.GT..5_ ABSZ)ARGN= ABSZ - ANU
CSP = CQSG(ARGNeABGR)
F-O.
IF(ARGN.LT.FLA(JH)IFeORAM(JH)eALOGI(CSP-FFZ)eFF2+Z.)
GoO.
IF(ARGN.LT.FL$(JH))G--HK(JH)eALOG(1.-((CSP-FG1)eFG3)*_Z)
P(NU)= F+G+CSLG(ANU*HABGR-90. pALFA(I)) VI
PI(NU)-P(NU)-AKeABS($ING((ANUeMABGR - 90.) - PHISH))
NU = NU + 1
ANU" ANU+ 1.
IFiANU-ANI(1))306,306p312
IF(ANU- ABSZ)31_p320P320
PHIM ANIiI]_HABGR "
VI=VI+CSLGiPHIM-9OotALFA(1))
I I+l
IF(I-1-JST)30_P302,304
PS=O.
B2-O.
00 326 I=lpNKR
PS=P$+P[I)
BI Z*lI-1)
B2 B2 + SING(BI_ABGR)*P(I!
Vl Zo_EXP(PSIABSZ)
SXI " ,00000000
SY-O.
00328 N=IpNQ
O=O.
00326 M'I, IB
MN N + I + MQ - Z_M
MM Z_N - MN
IF(RN.GT.NKR) MN - MN - NKR
IF(MM.LT.I) MM MM + NKR
O Q+ FKERN(M)e(PI(MN)-PI(MM))
ANU- N-1
ZP ANU_HABGR - 90.
ZL COSG(ZP - PHISH)
ZL = ASS([1,-ZL)/(1.+ZL))
IF(ZLoNE,O=IZLuALOGiZL)
&RG(N) = 0 - AKP_SINGIZP-PHISH)_ZL + ZP
VF(N) VI*EXP(-P(N))
WV COSG(ZP)IVF(N)
XP(N)= WVOSING(ARG(N))
YP(N)=-WV_COSG(ARG(H))
5Xl= SXl XP(N)
SY SY + YP(N)
SX " SXI
XPK " SX/(kSSZ -- 1.)
YPK SYI(ABSZ -l.)
I19
APPENDIX
DO 329 N-ZpNKR
XP(N)sXP(N)-XPK
329 YP{N).YP(N)-YPK
CALL CINT(XPeXpNQ_IZZ)
CALL CINT(YP,Y_NQPIZZ)
RQV O.
DO 330 Nm2pNKR
ROeX(N)tX(N)+Y(N)*Y(N)
IF(ROGTRQV)LsN
330 RQV = RQ
D0 327 I 1,3
IEPPL L-21
327 _(T)=SQPT(X(IEPPL)*X(IEPPL)+Y(IEPPL)*Y(ZEPPL))
333 TAU (R(3)-R(1))/(4.*(R(2)+R(2)-Ril)-R(3)))
XNAS X(L)+TAUtiX(L+I)-X(L-1),Zo_TAU_(X(L+I)+X(L-1)-X(L)-X(L)))
YNAS - y(L|+TAU*(Y[L+I)-Y(L-1)+Z.tTAU*(Y(L+I)+Y(L-1)-Y(L)-Y(L)))
SO XNAStXNAS + YNASeYNA$
AT"XNA$1SQ
Bs YNA$/$O
STREF - lolSQRT($Q)
ETA ABSZ*STREFIPI
CM .SeETASSTREFeB2
DARG 19,098_q S[3o_YNA$1XNA$ -- (YNASIXNA$)e_3)
IF(ABS(SX)ABS(SYi.LTo.OOOI_ABSZ) GO TO 335
$XSTREFeSX*200.
SY-STREFeSY_ZO0,
NZTNZPZ(ZpO)
WRITE (IDRUp33) NZTeSXpSY
33 FORffiAT (AZ_IHWARNING - $X spFb,3/3Xe_HSY -_F6,3)
335 D0 331 N=ZeNO
XR=X(N)
X(N)" 1.-B*Y(N)-ATeXR
Y(N)= BIXR -AT_Y(N)
ARG(N) ARG(N) - DARG
WQ (XP(N)+XP(N-Z)-XPK-XPK)ee2 (YP(N)YP(N-1)-YPK-YPK)e_2
331 DS(N-1) - STREFeSQRTiWO)e(1,+.bbb6667_((XP(N)eYP(N-1)
1-XP(N-1)tYP(N))IMQ)_2)
NHKW-NOI12
DLT - Y(NHKW)I(BOGEN_(I,-X(NHK_)))
DLTU--Y(NJ'K_)I(8_SEY*[I.-X(_H_)))
X(1) = 1.
33Z ARG(1) ARG(1) - DARG
3_6 ITPI
ALN=DARG
TF(PURES(13).GE.O.)GO TO 11
PURE$(1ZI-IO._HKS
PURE$(13)=,O0001
11 RETURN
120
........... ---_.
APPENDIX
END
SUBROUTINE CINT(P_Z, NQtIZZ)
DIMENSION P(NQ)_Z(NO)
FCINT=O.
IF(IZZ,NE.O) FCINT=.O8333333333
ZIZ)=O
PVVoP(NO-1)
PV=P(1)
PL'P(2)
DO 10 N-2, NQ
NZ-N+Z-(NQ-Z)*(NINQ)
PZ=P(NZ)
Z(N)=Z(N-1)+PLePV(PL+PV-PZ-PVV)OFCINT
PVVaPV
PVePL
10 PL=PZ
RETURN
END
SUBROUTINE QIP(X,Y;A)
DIMENSION X(3)pY(3)pA(3)
Cl - (Y(2)-Y(1))I(X(Z)-X(1))
A(3)-(Y(3)-Y(1)-CIO(X(3)-X(1))I/((X(3)-X(1))*(X(E)-X(1)))
A(1)=Y(1)-CltX(1)A(3)_X(1)*X(2)
A(2)-C1-A(3)*(X(1)+X(2))
RETURN
END
SUBROUTINE OIA(XpY_NPpO)
DIMENSION U(3I_V(3)_W(3I_Z(3I_A(3)_B(3IpX(1Z1)_Y(1Z1)
D,O,
XV-X(1)
00 4 N'2pNP
NR NP-I
XS=X(N)
1 IF(X(NR)-XS)3_ZwZ
2 NR u NR - 1
GO TO 1
3 IF(NR-N.LE.3)GO TO 5
DNN=Y(N)-Y(NR)-(Y(NRel)-Y(NR))O(X(N)-X(NR))/(X(NR+I)-XiNR))
IF(DNN.LE.D)GO TO 4
NMuN
NRMmNR
O=ONN
4 XV'XS
5 TA'I.
YF-Oo
IIT'O
? IIT'IIT*I
IF(IIT.GT,8)GO TO 11
DO 6 l'ls3
121
APPENDIX
ND-NM+2-I
NU"NRM-2I
U(I)sTA*(X(NO)YF*Y(NO))
V(I)sTA*(Y(NO)-YF*X(NO))
W{I)uTA*(X(NU)_YF*Y(NU))
6 Z(I),TA*(Y(NU)-YF*X(NU})
CALL QIP(UpVPA)
CALL OIPIWPZ_B)
IF(ABS(A(3)-B(3)).LT..O001) GO TO
XST - (B(Z)-A(2))_.51(A(3)=B(3))
YS - A(2)Z,eA(3)*XST
IF(XST.GE.U(1).AND.XST.LE.U(3))GO
IF(XST,GT.U(3))NR=NM-1
IF(XSTLToU(1))Nn=NMZ
GO TO 7
ZI IF(XST.GE.W(1).ANDeXSToLEoW(3))GD
IF(XST.GT.W(3))NRMmNRM+I
IF(XST.LToW(1))NRM-NRM-1
GC TG T
22 IF(ABS(YS)-.OOOIIg_Bp8
8 YF-YF+YSITA
TA-I./SQRT(I.YF*YF)
GO TD 7
11
TO Zl
TO 2Z
D AI1)-B(1){A(Z)-B(Z)(A(3)-B(3))*XST)eXST
11 RETURN
END
SUBROUTINE DIAGR[ISTIFTpISTIM_ENV)
CDMM_NIBRALIDLT_OLTUpALN_ALV(14)sNAL_ITP_NAMP(ZZ)pCML(16)pCRL(l_)
1 ,CPV(2)pALTX(4p2I_DARFpITITI_ITIT2
COMMON Pl(1ZI)pP(IZ1)eXP(121)_YP(IZ1)_PUFF(l_)_AGAM(14)pX(lZl)p
1Y ( 1Z1 ) _ DS ( 1ZZ ) VF ( 121 ) p ARG ( 121 ), ANI ( 28 ), ALFA (29) p I ZZp KFU_ NO_ NUPRO
EXTERNAL TRUEp YG%p YG5
LOGICAL PLpTRUEpYG%pYG5
DATA RZEPALMIN,ALMAXI177.8p-5._15.I
PL-ISTIFT.GE.5
IF(ISTIMoGE,Z)GO TO 10
XMX'O,
IF(MENV.EO,1IGO TO 6
IF(PUFF(1).NE.O.)RZE-IOO**PUFF(1)
QUER-Z.3_RZE
HOCH-1.SeRZE
CALL DEFINE(PLpQUERpHOCH_-.Z_Ilp'Z_I,3)
GO TO 10
4 ALMI-O,
ALMAsO,
IF(PUFF(1).NE.O.)ALMINsPUFF(t)
IF(PUFF(Z).NE.O.)ALMAXsPUFF(Z)
HOCH-IO.eiALMAX--ALMIN)
CALL DEFINEIPL_Z6O._HOCH_-,3_4.3_ALMIN_AL_AX)
122
.T
l
APPENDIX
10 DO 348 _=ZsNAL
VMX'Oo
00 369 Nel_NQ
V'.5*A8S(VPR(NpH))
VQ'V*V
IF(VMX.LT. VQ)VMX=VQ
IF(XMXLTV)XHX'V
349 P(N)=V
IF|HENV.NE.I)CALL POLZUG(ISTIFTpXeP_NCp. TRUE.,YG4)
348 PI(N)-4**V_X-1.
IF(MENV.EQ,I)GO TO ZO
DO 350 I'IpNQ
350 P(I)-Y(I).Z
CALL $MOOTH(ISTIFTjXpPpXP,NQnTRUE)
GO T_ 30
ZO CALL POLZUG(ISTZFTpPlpALV_NAL_,TRUEJTRUE)
ALHI=AMINI(ALMIsALV(1))
ALMA-AMAXI(AL_A_ALV(NAL))
30 IF(ISTIM.EQI.0R.ISTIMGE.31RETURN
IF(MENV.EQ.I)GO TO 40
XMX-AMINl(XMX_l.25)
CALL ACHS(Z_O.P-.OI_.OZpO.jXMX,.OSj.ZStYGS)
CALL ACHS(lpO,_-.OI_-.OZ_O._I,P.lp.SpTRUE)
GO TO 50
40 CALL ACHS(ZJO,p-.Zp-.5_O.p4.e.Z,.SgTRUE)
CALL ACHSIZ_O.e-.O49-.l_ALNIwALMApI.,5._TRUE)
50 CALL GCLOSE
RETURN
ENO
LOGICAL FUNCTION YG6(XjY)
A-X
YG4-Y.GE.oZ.AND.Y.LE.I.Z5
RETURN
END
LOGICAL FUNCTION YGS(XpY)
YGS=X.LT.-.OOOI.OR.AES(Y-.Z).GTOZ5
RETURN
END
SUBROUTINE PUDECK
CONMQN Pl(121)_P(121)pXP(ZZ1)_YP(121)pPUFF(16|_AGAN[16)pX(121)_
1Y(121)_O$(1ZZ)_VF(121)pARG(1Z1)pANX(Z8)eALFA(Z9)_IZZpKFUpNQpNUPRO,
ZJABeJSTpCMpETApABFAePZpBOGEN
COMnONIPRALIDLTeOLTUtALNpALV(14)_NAL_ITP_NA_P(12)tC_L(14)eCRL(14)
1 ,CPV(Z)_ALTX(4_Z)_OARF_ITITI,ITITZ
CO_HONIEAIILESeIDRU_ISTA_NNESE
DIHENS|ON INX(4)_INY(4)_INA[4)
EOUIVALENCE(P(1)_INX(1))_(P(5]_INY|I))e(F(9)_INA(1))
$TZ-IO0000.
IF(X(1)GT.80.)STZ-IO00.
123
APPENDIX
124
DO 6 N=lpNQp6
NZE-HINO(k,NQ-N+E)
DO 4 H=E,NZE
L=N+M-1
INX(H)=INT(STZ*X(L)+SIGN(.SpX(L))) r
INY(M)-INT(STZ*Y(L)+SIGN(,SPYIL)))
ARGS-ARG(L)*BOGEN
IF(ITP.EQ.Z)ARGS-XP(L}
4 INA(M)-INT(IOOOO,eARGSSIGN(.SpARG$))
b WRITE(ISTA_B)iNARP(M)pMn3_T)pNe(INX(M)gINY(M)JINA'(M)PR=I*NZE)
8 FORMAT(SAleI3PlZIb]
RETURN
END
SUBROUTINE STRAAK(TE_RUApYB,MXZpISTIFT)
DIMENSION XA(Z)_YA(Z)
COMMON Pl(121)tP(1Z1)tXP(121)tYP(1Z1)pPUFF(14)tAGAM(16)pX(121)_
1Y(1Z1)pOS[1ZZ)pVF(121),ARG(XZ1)pANI(Ze)_ALFA(Zg)pIZZeKFUpNQ
COMMONIPLTMIMPL_NGCtXZEH,YZEH_MSPLI
EXTERNAL TRUE
LOGICAL TRUEPPL
DATA YV, YXAIIOOOO.tlO000.1
T " ABSKTE)
PL'ISTIFT.GE.5
YBL=AB$[YB)
QUER=T+10.
QRM=OUER-S,
HOCHgYBL
YMX O,
YMN " O.
DO 2 I m I_NO
IF[Y(I).GT.YMX) YRX - Y(I)
IF(Y(I),LT,YRN) YRN " Y(I)
Z CONTINUE
IF(HGC)3pb93
3 IF[TE)Splb,8
5 CALL GCLOSE
ZF(T-Z.Ilbpbp6
b CALL DEFINE(PL_QUER_HOCHp-5.pQRMpOopHOCH)
YXA 5o-T*YMN
MSTR=O
GO TO 10
8 YXA YXA RUA
IF(YB,GT.O.) YXA RUND(YV-T*YRNI,5*RUA_I.IRUA)
10 YV YXA TeYRX
IF(YVoGE,YBL.AND.MSTR.NE,O}GO TO
_STR=I
DO 12 I=Z,NQ
PI(I) - T*X(I)
12 P(I) TtY(I)YXA
APPENDIX
i
J
]
i
i
14
16
=;F7,4)
NZT=NZPZ(lJO)
MRITE(IORUp4)NZTp(Mp_-N,NZL)
4 FORMAT (AIP?(IZ_ZTH X(CM) Y(CM) ))
NZTsNZPZ(lpO)
IF(NZT,EQoNNESE)GD TO 1
Z=I+l
00 6 K " N,NZL
L s K-N+1
XOR(L) = X(I)*.I*ABS(T(K))
6 YOR(L) Y[I)*.I*AB$(T(K))
WRITE(IORU, IO) NZT_(XOR(M)_YOR(M)pH-1,L)
10 FORMAT (A1,7(Fg. ZpF8.Z_ZX))
IF(I.LT.NQ) GOTO 5
8 CONTZNUE
RETURN
END
CALL SMOOTH(ZSTIFTtPl_P,XP_NQpTRUE)
IF(HXZI16p16p14
XA(Z) T
XA(2) O.
YA(1} YXA
YA(Z) YXA
CALL PQLZUG([STIFTpXApYApZp.TRUEopTRUE]
RETURN
ENO
SUBROUTINE 3TRDR (T_NT)
OINENSION XDR(8)p YOR(8),T(4Z)
COMMONIEAI ILE$# IORUe ISTAe NNESE
COMMON Pl( 1Z1)pP (1Z1 ]_ XP(1Z1) tYP(1ZZ)P PUFF (14) pAGAM(16) _X (121),
1Y(1Z1)pOS{ 12Z) eVF(1Z1)eARG(1Z1)_ ANZ( ZS)p ALFA[ Z(;) _ IZZ_KFU_ NQp NUPRO_
ZJA8_ JSTpCM_ ETAeABFA_ PI_,BOGEN_ DARG
EQUIVALENCE (XDR (1) J YP(lob) )p (YOR (1) sYP(116) )
DO 8 N Z,NT,8
I = 0
NZL N+7
IF (NZL.GT. NT)NZL NT
NZTNZPZ(3pO)
1 WRITE(IDRU_,Z)NZT_,NUPRO_,DARGJCM
2 FORMAT (AlsBHAIRFOZL ,I_ie.3XeSHALPHAO ,,eFS.2tSH DEGREESe3Xp
*SNCMO
125
P,F..FE_N_S ,-
i. Abbott, Ira H.; Yon Doenhoff, Albert E.; and Stivers, Louis S., Jr. :
Summary of Airfoil Data. NACA Rep. 824, 1945. (Supersedes NACA
WR L-560. )
2. Abbott, Ira H.; and Von Doenhoff, Albert E.: Theory of Wing Sections.
Dover Publ., Inc., c.1959.
3_ Schlichting, Hermann (J. Kestin, transl.): Boundary Layer Theory. McGraw-
Hill Book Co., Inc., 1955.
4. Smetana, Frederick O.; Summey, Delbert C.; smith, Neill S.; and Carden,
Ronald K.: Light Aircraft Lift, Drag, and Moment Prediction - A Review
and Analysis. NASA CR-2523, 1975.
5. Eppler, R.: Direct Calculation of Airfoils From Pressure Distribution.
NASA TT F-15,417, 1974. (Translated from Ingenieur-Archiv, vol. 25,
no. i, 1957, pp- 32-57.)
6. Lighthill, M. J.: A New Method of Two-Dimensional Aerodynamic Design.
R. & M. No. 2112, British A.R.C., 1945.
7. Hess, John L.: The Use of Higher-Order Surface Singularity Distributions
To Obtain Improved Potential Flow Solutions for Two-Dimensional Lifting
Airfoils. Comput. Methods AppI. Mech. & Eng., vol. 5, no. i, Jan. 1975,
pp. 11-35.
8. Eppler, R.: Practical Calculation of Laminar and Turbulent Bled-Off
Boundary Layers. NASA TM-75328, 1978. (Translated from Ingenieur-
Archiv, vol. 32, 1963, pp. 221-245.)
9. EFpler, Richard; and Somers, Dan M.: Low Speed Airfoil Design and Analysis.
Advanced Technology Airfoil Research - Volume I, NASA CP-2045, Part l,
1979, pp. 73-99.
I0. Kowalik, J.; and Osborne, M. R.: Methods for Unconstrained Optimization
Problems. American Elsevier Pub. Co., Inc., 1968.
ii. Eppler, R. (Francesca Neffgen, transl.): Laminar Airfoils for Reynolds
Numbers Greater Than 4 x l06. B-819-35, Apr. 1969. (Available from
NTIS as N69-28178.) (Translated from Ingenieur-Archiv, vol. 38, 1969,
pp. 232-240.)
12. Eppler, Richard: Turbulent Airfoils for General Aviation. J. Aircr.,
vol. 15, no. 2, Feb. 1978, pp. 93-99.
13. Squire, H. B.; and Young, A. D.: The Calculation of the Profile Drag of
Aerofoils. R. & M. NO. 1838, British A.R.C., 1938.
14. Woods, L. C.: The Theory of Subsonic Plane Flow. Cambridge Univ. Press,
c.1961.
126
z-plane
v=l
X
li
{-plane
Figure I.- Mathematical planes.
127
w
w
_ pressure recovery
t
i I ] clo_ure
i i i contribution
I I I ,
X X C X
W 5
Figure 2.- Typical pressure-recovery distribution.
128
............. Z. , ......
i,le
i
Cw - I.
2_ s
*_oi
\\\
Figure 3.- Arc limits on _ unit circle.
129
V
U
cc
ct,deg
l 2.7 ;us and ls
l.O-
0.5
i _ t i j
0.5 1.0
x/c
Figure 4.- Airfoil 1095. (_ relative to zero-lift line.)
130
1.0
V
0.5
.ORIGTNAI_ PAGE IS
OF POOR QUALITY
a, deg
3 69 ;us and Is
/
0.5 1.0
x/c
Figure 5.- Airfoil 1096. (_ relative to zero-lift line.)
i
; 131
_, deg
12
0
IO
8
V
U
1
12
i
I i i i
Figure 6.- Airfoil 1098. (_
i | i i l i
0.5 1.0
x/c
relative to zero-lift line.)
132
2.0_
V
_._IG_-_'L QUA__'_-_
p_ Eoo_
I$
Figure 7.- Airfoil 1097.
x/c
(_ relative to zero-lift line.)
133
Pn
Panel between Pn and Pn+l
Shape: Cubic in q({)
Vorticity: Parabolic
Flow condition: Inner tangential velocity
equals zero at Pn and Pn+l
Figure 8.- Panel method.
i
134
(_,n)
..... _ W
WrI
Figure 9.- Induced-velocity vector.
135
Figure i0.- Induced tangential velocity on panel.
136
P2
P
nq-1
Trailing-edge angle equals zero
Flow conditions: Equal velocity on both sides of
trailing edge; normal velocity equals
zero at trailing edge
Figure Ii.- Sharp trailing edge.
P1-P
nq
P2
P1
P
nq-1
Pnq
Linear vorticity and source distributions at the base
No flow singularity at PI or Pnq
Flow conditions: Inner tangential and inner normal
velocities equal zero at the middle
of the base
Figure 12.- Blunt trailing edge.
137
138
R62
104
lO3
lO2
lO l
_nt
separation
Transition
criteri on
(r=O)
_, deg
13
Laminar
separation
1.50 H32 1.55 1.60
Figure 13.- Boundary-layer development; airfoil 1098, upper surface,
R = 1 x 106 . (U relative to zero-lift line.)
i
Y
I
I
I
I
I
v I
X
Ssep
u/f
oo
Figure 14.- Lift-coefficient correction due to separation.
Subroutine GRP
Subroutine GRS
Subroutine H12B
Subroutine UMP
Subroutine GRUP
Subroutine CDCF
Figure 15.- Block diagram of boundary-layer subroutines.
139
!
\
\
\
\
\
\
\
% %
J
J
I
!
o
, r , J _ i , i
|
140