SYMBOL DESCRIPTION

THETA Crank Angle (in degree), zero at TDC
P Pressure (atm)
Vv Instantaneous engine (or charge) volume (cm)
T Mass averaged temperature CK)
U Mass averaged specific internal energy (cal/gm)
PHEE Over-all equivalence ratio
DPDTH Step-averaged pressure derivative, a (atm/C.A.)
WP Step-averaged work rate, P * (cal/C.A.)
DFJDTH Fuel injection rate (gm/C.A.)
SMDDTH Step-averaged over-all dilution rate (gm/C.A.)
SMVDTH Step-averaged over-all vitiation (air) rate (gm/C.A.)
DOQRDTH Step-averaged over-all rate of radiative heat transfer to the
cylinder wall (cal/C.A.)
SVBDTH Step-averaged over-all premixed burning rate (gm/C.A.)
SDBDTH Step-averaged over-all droplet burning rate (gm/C.A.)
SLDDTH Step-averaged over-all mass rate of-ignition of liquid fuel (gm/C.A.)
SUMQL Cumulated (time integrated) heat loss(convective and radiative heat
transfer) to the cylinder wall. (cal.)
SUMFB Cumulated mass of fuel burned to equilibrium products
DFBDTH Over-all mass burning rate (gm/C.A.)
WKNET Cumulated net work (cal), i.e., fe * de
oN
RR Mass averaged gas constant (em? - atm/gm - °K)
UCHAR Characteristics charge velocity (cm/sec)
3
DVDTH Engine Volume rate (cm /C.A.)
DELE12 Over-all numerical error involved in integration of the system energy

conservation equations, (normalized by total chemical energy released
i.e. XMFJ * HCOM).

19


APPENDIX: FORTRAN LISTING OF DIESEL ENGINE PROGRAM

The subroutines follow the main program are listed in alphabetical
order.

43


9€

THETAIDEG CAPs = 500 = PEATNDs 7.A5GEC0L VICHSD= 1.5166002 TKKI= 1,6196003 U(CAL/GH)® S.B77E001 PHEE= 5.96:

DRATH CET Ms Mb eWSGE-OL MPECAL/CAD= S6A22E9C9 DFINTHIGH/CAD> 0, SHDOTHIGN/CAD2 3.248£-02 S4VOTHIGH/CAD= 3,228E-C3
DACUTHIGOL/CADS 3.830L600 OORDTH(GAL/CAD= 6.7245-01 SVODTHIGH/CAD= 0. SOBDTHIGH/CAD= 2.089E-93 SLOOTHIGY/CAI= 6.
SUACLICALD= 7.935E061 SUMFU(GH)= 7.657E-02 SFADTH(GH/CAI= 2.0N9E-03 WKNET(CALD=-3.091E402 RRICH3-ATH/GM-K)= 2.5966 +00
UCHROICH/SECIS GS31E03 OVATHIGM3/CADs 2.229€400 MELEL2= 4.657E-03 ITERPCE 1

XMNOF(GHD= 46974E-03 YNOE= 1,959E-03 XNOE(PPHD= 1,648£403 SOOTE(GH)= 3,039E-04

THETO(DEG .CAP= 7.60 | PUATHD= 7.575C001 VACH3I= 1.569802 TIKI= 166356403 UICAL/GMD= 5.109E401 PHEE= 5.9632~31
OPUTHEATH/CAD==.,023F6G0 MPICAL/CEI= S.392E60C OFIOTHIGH/CAD= 0, SMNOTHIGH/GA)'S 3.2106-G2 SHVOTH(GH/CA)s 3.2186-03
OGEOTH(CALJCAD= 3.813£000 DOROTHICAL/TAY= B.653E-01 SVUUTHIGH/CADS 0. SOBDTHISH/CAD= 1.821£-93 SLOOTHIGM/CAD= 2.

SUNCLICALD = O6922E4G1 “SUMFD(GH)= 7.99 TF 02 SFODTHIGH/CAD= 1,421E-03 WKNETICALI=-2.99LE602 RRICNI-ATH/GH-K)= 2.8976 000
UCHACICRSSECIS &.533E003 OVOTHICMS/CAD= 3,101E+00 NFLF 122 &.3720-03 ITERPCE 1
KRNOE(GHD= $o239E-03 YNOES 2.060£-03 XNOE(PPMD= 1.944E+03 SOOTE(GHD= 3.154E-64

THETACQEG.CAI= = 9.60 PLAT) = 723935401 VICH3D= 16406462 TUK

1.643F 403 ULCAL/GND= &.233E40L PHEE> 5-6563Z-01

DPOTHCATM/CAI=©1.595E400 WPI CAL/CAD= 647995400 OFJOTHIGH/CAD= 9, SHODTHIGH/CA) = 3.172£-02 SHVOTHIGH/CAIs 3.172:
OSCOTHICAL/CAY= 3,754E609 DORDTHICAL/CA WS3E-0% SVOODTHIGH/CAI= 9, SOGOTHIGH/GAD= 1,586E-03 SLOOTH(GH/CAI=
SUMILECELD= 968475471 SUMFOIGMD= 8, 3125-52 SFOUNTHIGM/CAI= 1.5862-03 WKNET(CAL)=-2.8625¢02 RRICHS-ATM/GM-K}= 2,837E000
USAAF (CH/SEGIE GA5S4VEFO3 OVOTHICHS/CAD= 3.977E4G0 MELEL2> 4.d49E-03 ITERPC= 1

AMUOE(GADS 5.452F-03 YNOEZ 2.161E-03 XNOE(PPM)= 2,d20F403 SOOTE(GHI= 3.225F-04

THETRIDEG.CAP= 12,00  PUATHD = 7.032Fe OL VECHS)= 26729402 TEKI= 1,646E403 UICAL/GHD= 3.256001 PHEE= 5. 9635-08
CEATUIATHZCRDs91,.A76F 009 WP(CAL/GAD= 769876400 OFIDTHCON/ZAI= GC. SHOOTHIGH/CA}® 3.086£-02 SHVOTH{GM/CAD= 3.006E-63

ACOTHICAL/CAI= 3.661€000 DOKOFHICAL/SAD= B.152E-01 SVAUTHIGM/CAD= 0. SOBOTHIGH/CAD= 2.326£-93 SLOOTHIGH/C!
SUNCLICALD= 1,)75F462 SUMFAIGMD= 6.590F=G2 SFYOTHIGH/CEI= 1.528E-03 WKNETICALD=-2.708E402 RRICH3-ATH/GH-K= 2.8966 400
UCHBRICH/SECI= &.537E403 OVOTHICM3/CAD= 4.845E40G DELLI2= 3.491E-03 ITERPC= 1

XMNCE(GH)= S.615E-03 YNOEs 2.233E-03 XNOE(PPMI= 2.078E+03 SOOTEIGMD= 3.2616-04

a.

THETAUDEG.CAD= 42400  PEATHD= 6.82BZe01 VICMSI= 1.786002 TAKD= 1.644E403 UMCALJGHD= 2,7466001 PHEES 5.95
or THEATM/CQ1S=-2,958E000 HAICAL/CAD= A.UO2E40O NFJOTHIGH/CAI= 0. SHOOTH(GH/CA) = 3.035E-02 SHVDTHIGHs

an
Vs 3.035E-a3

THICAL/CAD J.4C4E000 DQROTHICAL/CAI= 7,940E-O1 SVENTHIGH/CAD= 2. SOGOTHIGH/SAD= 1,210€-03 SLOOTHIGH/CAI= ce
Sancutceuys As119E4G2 SUMFAIGMD= B.7112-02 SFAOTHIGH/CAI= 1.216b-03 WKNETICALI=-2.623£002 RRICMI-ATH/GH-K)= 2.3966 400
UCHAPICH/SECI= &.SS8E6C3 OVOTHICHS/CAD= 5S.276E400 OFLEI2Z= 3.315€-03 ITERPC= 1

XMNOE (GIS 56644E-03 YNOE= 2,229E-03 XNOE(PPMI= 2,103E¢63 SOOTEIGMI= 3.272E-64

THETAUDEG.CAD= 13.90 — PUATHD= G.613E4U1 VICHSI= L.A33E002 TEKI> 1.66160G3 UICAL/GHD= 2,214E401 PHEES 5.453£-01

OPOTHIATH/CAI=-2,117E469 W>(CAL/CAd= 8.926E400 OF UNTHEGH/CAD= a. SHNOTHIGH/GA) = 2.975E-02 SHVDTHIGH/CAD= 2.975€-03
MECCTHICAL/CAI= 3.541F600 NOPNIHICAL/CAD= 7,725E-01 SVONTH(ON/CRIS Oe SOUDTHISH/CAI= 1.105£~03 SLONTHIGH/EA)= 9.
SUVSLICALI= 2,162€462 SUMFO(GHI= 8,622f-22 SFAOTHIGI/CAD= 1.21650-03 WKNETICALD=-2.534£+02 RRICMS-ATH/GH-K)= 2.8958 606
USHAPICH/SECDS 540003 OVOTHICMI/CAD= 5.704EFGG OFLELZ= 26936£-03 ITERPC= 1

XMNOE(GMD= %.742E-03 YNOES 2.2526-03 XNME(PPMD= 2.120.463 SOOTEIGHI= 3.277E-C4

THETAUDEG.CADS = 15.09 — PUATH)= G16KEHU1 VICHID= 1,.9566002 TIKD= 1.633003 ULSAL/GMD= 1.139E401 PHEES 5.963E-01
CPOTHCATA/CAD=22,235E409 MPICAL/CAD= F.644E400 OFINTHIGH/CAI= 1. SHODTHIGH/CAD= 2.84GE-02 SNVDTH(GH/CAD= 2.846E-03
GACITHICAL CONS 3.402E000 MQKOTHICAL/SAI= 7,257E-G1 SVGOTHIGH/CAI= 0. SOBOTH(SH/CAI= 9.021E~04 SLOOTHIGH/CAI= Co
SUMGLICALD= 1,246E092 SUMFACGHD= 9,012°-02 SFADTHIGH/CAI= 9.C21E-04 WKNET(CALD=-2.344E402 RRICHS-ATH/GN-KI= 2.694E 000
UDMAPICH/SECD= G544F 403 OVNTHICHI/CAD= 6.553400 OFLFIZ= 2.708€-03 ITERPCE 1

VENCEIGHD® SoA2GE“G3 YNOT= 2. 282E=93 XNOE(PPMI= 2.154£463 SOOTE(GHI= 3.280E-04


QPF= C*(DLPF+*O2PF +O 3PF)
PRAPF= -F1*F 1 FF AROS

PUPP= F1"*O*(030/P1)

PuPT Cv

PUPF= F1I*{-24+09F)+ PFIPFFUX

PRP FASCF#CRR2 (5) FPHE#RR2 (8) /P1

PRP ~FLFC*(FR2CG) FPHE#R22(7)) ¥1000./T2

PRXPF= F822(2)+ CF(RP2(3)+ RRZ2(G)FTI+t RR2(5) FALCGP)
PRPF= F1i*PRXFF+ PFIPFFERX

RETURN

END

94


UF J=(TFJ-TBEPF) #CPFLEUFLBP
HFJ= UFJ+ PINISFOENFATOXS
© CRITICAL FUEL PROPERTIES
UFCRIT= (TCRITF-T8PF) *CPFL+ UFLBF
C ENGINE DISPLACEMENT AND CLEARANCE VOLUME
VDIS= 0.5*3,14159*DE*DE*REW
vos vOIS/(CR-1.9)
FARORD= XMFJ/XMAT
QFRORD= 1.5/FARORD
EGR1iI=1.6-EGR
XMEGR=EGR*XMAIZEGRI
XMAEG XMEGR/ (1. 0#F ARORD)
XMFEGR= F ARORD*#XMAEGR
c= 1.0

THETAT= THEVC
CALL VOLM(THET4I,OSOTH,DVOTH,V1)
" WP= PL*DOVOTHFATOXS
IF (XMRD.GT.0.0) GO TO 4
C EVALUATE RESTOUAL SAS VIA 4 SIMPLIFIED OTTO CYCLE (OPTION)
v2= V1/CR
Gi2= 1.35
G34= 1.272
CVi= 6.2514 '
XMROL= £29
00 5 ITER=1, 29 .
XMAO= XMRO+XMAT +XMEGR
Ti= PLFVLFAMWAT/XMAO/R
T2= TA*(CR)¥F(G 12-1. 0)
C ASSUMMING 25 PER CENT HEAT LOSS
T3= T2+ HCOM/CV1/ (OF RORDN41.0)*9.75
(XMAO +XMF IS) FRFTZSOMWAT SV 2
Pus PIF(CR)**(-534)
V4P= V1*C(PQ/PE)F*(1./634)3
RSD= V2sV42
XMRO= RSO*(XMAO+XMF UD
IF (ABS(XMRD-XMRD01).LT.0.01*XNRD) GO TO 4&
XMRD1= XMPD
5 CONTINUE
WRITE (6,117)
stop
& CONTINUS
XMAO= XMRD+ XMAT+ XMEGR
XMARDQ= XMRO/ (1. 6+F AR ORD)
XMFRO= FARORD*XMARD
G EQUIVALENCE RATIO FOR #A1IR# (I.E. INTAKE PURE AIR # EGR + RESIOUES }
FARA= (XMFEGR+XMFRDI/ (XMAT#XMARD+XMAEGR)
AFRA= 1.9/FARG
PHEA= FARASFARIS
S INITIALIZATION
TIM=1.0E-15
BTAA=1.0
Ga=1.0
XMA=XMAO
XMPA= XMA/(1.C+A4F RA)
XMPAA= XMA/C1,0 +FARAD
T1i= PLFVIFAMWOT/XMAOZR

48


Table 3.

Fuel Properties

WwW

4K
FDEN
CEN
API
XMUF
SURFT
TBPF
TCRITF
CPFL
CVFV
CPFV
DHFV
HFFL
HCOM

No. of C atoms in fuel molecule

No. of H atoms in fuel molecule

Density of fuel (gm/cm>)

Centan No.

API Gravity at 60°F

Viscosity (gm/cm-sec) at TF]

Surface tension (gm/sec*) at TF]

Boiling point at 1 atm (°K)

Critical temperature (°x)

Liquid specific heat (cal/gm)

Vapor constant volume specific heat (cal/gm)
Vapor constant pressure specific heat (cal/gm
Latent heat at 298°K (cal/gm)

Heat of formation at 298°K (cal/gm)

Heat of combustion (cal/gm)

ia; ie? a, as a a a a) ie? es a: Mas es Bas Ba

NA

NA
Cption
NA

NA

NA

NA

NA

NA

NA

NA

NA

NA
Option*
Cption*

* If this value is not specified it will be calculated by program.

12



FOECK HOMNOR
SUBROUTINE HCMNOR (Ps ToXNCoXN29X02,XNCNEG,CNCRyXNOMTIL)
C THIS FOUTINE CALCULATES HOMOGENECUS NC PRCDLCTICN RATE
RT= 1.987*T
C RATE CONSTANTS

XKF 1s 1.35E14* EXP (-75400./RT)

XK31= JA 2ZELSFEXP(-4NO./ RT)
L.SSELIFTFEXP(-7080./2T)
3.2E9FT FEXP(-39100.7°T)

= S.O*EXP(-58300./2T)

XK1= CNORF0.229782* (XKFL*XKECT)

XK2 = UXKBLFXKB2/ (XKF1FXKF2) )

XK3= ¥KOL/SXKF2

CORAL= 1.0- XK2* (XNONEQ*XNCNEOS UXN2* XC2))
CORR 1.04 XK3*(XNONEG/YO2Z)

GORR=s CCRRLI/JCORS2

XNOOTL= XKLFXN2*SQRTUXO2*P/T) ¥CORP
RETURN

END



68

7a

OMDAVX= G.5*(OMOOTH(ID+OMO0TO(ID)
OMVAVX= 0.5*(DMVDTH(I) #OMVOTO(T))
OLBAVX= 0.5*(DLSO0TH(I) +0L80TO(T))

XXX1= HAO*(OAPOTH(I) 40MOAVX ~OMVAYX 2+ HFOSOFOLBAVX

1 + HFGO*DVGDTH(I)- QBMAVX
XXX2= 0 XXX1L*OTHETASXMBM(T)
RHP= RBM(T)

THP= TOMCT)
00 68 ITER=1yNTER
USMIT= XXX2- RHP*THP

CALL TPROPITHP» Py VXHP,CVHPsSPHP yUHP, HHP yAWMHP yGAMHP,

4 RHP,PHECT))

XXX= UBMIT=- UHP

IF (ABS(XXX/UBMIT) «LT.0.062) GO TO 69
THP=THP+#XXX/CVHP

CONTINUE

WRITE (6997) I, ITERPCsTHP,UBMII», UHP

sToP

CONTINUF

T3M(L)= THP

CONTINUE

XMNOCT)= XMNOOCI)+ (ONOOTO(T)+ONODTH{(I)) *DTHET2
IF (XMNOCT?.LT* 960) XMNOCTI= £29

XMSOOT(I) = XMSOTOC(I) + (OMSOTO(T) #OMSOTH(T)) *OTHET2

IF (TBM(T).LT.TRANGL «OR. TSM(I) .GT.TRANG2) WRITE(6,109)

LTBM(L) ,PHE(I),OTBOTH(I), OTBOTOC(I), TS8MO(T) ,OPOTH
2 ,OLECTH(T) ,OV30TH(I) »OMDDTH(T) yOMVOTH(T)

IF (PHE(I) LT.PHERGL «OR. PHE(I) .GT.PHERG2) WRITE(6,169)

1 TOM(T),PHECT)

YNO(I) = XMNO(I) /XMBM(T)

XNONEQ(I) = CAMW3M(T)/30.)*YNO(T

YSC(T)= XMSOOTUII/XMSM(T)

BTA8M (T= XMBM(I)* RBMCI)* TBMCID
GBM(T) = BTABMCIY**(,6666667) *(TBM(I)-TW)
SUMG=SUMG +GBM(I)

SUMU= SUMU+ XM3M(T) *UBM(T)

SUMRT= SUMRT+8TABMIT)

SUMR= SUMR4#XMBM (IT) *RBM(T)

CONTINUE
3 CONTINUE

IF (XMO.LE.0.d) GO TO 35

BTAASBTAA/SUMRT

GA= GA/SUMG

CONTINUE

IF (XMFG.LE.G.0) GO TO 36

BIAFG= @TAFG/SUMRT

GEG GFG/SUMG

CONTINUE

IF (II.LT.t) GO TO 74

O00 75 I=1,IT

IF (XMFAM(T) LE.O.0) GO TO 75

BIABM(T)= BTABM(TI/SUMRT

GSM(I)= GBMCT) /SUMG
5 CONTINUE

CONTINUE

SUMM= XMA4+ XMFG+ SUMBM

SUMSKM= XMAt XMFL# XMFLB+ XMFG+ SUMBM

CHK= ABS CSUMCKM-SUMF J-XM40)

IF (CHK.GT. 169€-5*XM40) WRITE (3,118) SUMCKM,CHK

1 »XMQ,XMFL4XMFLB»XMFGySUMEM

58

I,THETAT,

I,THETAT,


XFIOTO(I)= XFIOTH(T
XAIDTO(I) = XALOTH(T)
NO 83 L=1,NSZ2G
XKKOCI,LY= XKKCTyL)
83 CONTINUE
53 CONTINUE
54 CONTINUE
SMDDTO= SMOOTH
SMVDTO= SMVDTH
SNBOTO= SOBOTH
XA&9TO= XAAOTH
XFADTO= XFADTH
DFJOTO= OF DTH
OFVOTO= DFVOTH
SLOOTO= SLODTH
SVB0TO= SVBOTH
APBOTO= APBOTH
ogboto= 0acoTH
DORDTO= DOROTH
THETAT=THETAI#0 THETA
CALL VOLM(THETAI,DSOTH,OVOTH,V)
Coeees CORRECTOPS. +. SORECTORS
53 CONTINUE
© MASS CONSERVATION
SUMFJ= SUMFJO# (DFIDTHFOFJOTO) FOTHET2
XMFV = XMFVO + (OFVOTH+OFVOTO? *DTHET2
XMFLO= XMFLOO# (SLOOTH+#SLOOTO) *OTHET2
XMFPR= XMFPBOF (SV3DTH+SV3NTO) FOTHET2
XMAPB= XMAPBO+ -CAPBOTH+APSOTO) FOTHET2

XMO = XMDO + (SMOOTH#SMODTO) *OTHET2
XMYo = XMVO + [SMVDTH#SMVDTO) *OTHET2
XMFIQB= XMFOBOt (SOSOTH+SOBOTO) FOTHET2
XMFA = XMFAO + (XFAOTH#XFAOTO) FO THET2
XMPAA= XMPAAO# (XAADTHFEXANDTO) *DTHET2

QFRA= XMPAA/XMFA
FaRa= 1.0/0FRA
PHEA= FARA/FARIS
SUM3M=0.0
XMFLB= 0.90
If (II.LE.0) GO TO 8
00 9 T=1,IT
XMF3AMCT)= XMFBMOCID+ (XFIOTH(ID+XFIDOTOC(I)) *O0THET2
XMPAM(I)= XMPAMO(I)+ (XAIDTH(I)+XAIDOTO(I)) *OTHET2
XMBN(T Y= XMPAMCID4+XMFBM(I )
AFR(T)= XMPAM(TI/SXMFRM (TY
PHE{T)= 1.0/7 (AF RCID*FAROS)
SUM3M= SUMBM+ XMBM(TI)~
XMFLI(I)= 0.0
OO 84 L=1,NSZ2G :
B2= OLO(I»LIFOLOCIyL I= (XKKCTyL) #XKKOCI,LI) #0T2
IF (02.LT.1.06-16). O2= 0.6
DI(I,LI= SART(N2)
OMPLICI,LI= XNOI(TyLY*(0.5235987757 *FDEN*OT (I yL)**3)
XMFLICIDV= XMFLICID+ OMFLI(I,L)
@L CONTINUE
XMFLB= XMFLB+ XMFLICI)
9 CONTINUE
3 CONTINUF
XMQ = XMAQ +XMV-XMD-XMAPB
XMFL= SUMF J-XMFV-XMFLD

56


1G2 FORMAT (1X p#L FUEL #,1P9E12.3,0P2F7.4)

104 FORMAT (1X ,*G FUEL *,1P9E12.3, 0P2F7.4

105 FORMAT(1X,* AIR *,1P9E12.3,0P2F7.4)

166 FORMAT (1X ,14,2X,1P9E12.5,0P2F7 4)

197 FORMAT (1H4)

20% FORMAT( SX,1P11511.3)

109 FORMAT (1X,15,1P 1210.2)

112 FORMAT (//50X,*M OLE FRACTIONS *//1Xy *SYSTEM* 44Xy FH2¥*y BX 9 FH20% 97X49 F50
1 48Xs FCC] Fy 7X, FN2F BX ys FNOF 9 BX sy FOHF BX, FHF, IX FOF GX *NF 9X yg FO2Fy
2 9Xe*CHE*)

111 FORMAT (1X,*FORMATION OF A NEW SYSTEM*,1I5,1P10210.3)

112 FORMAT (/5Xy*L¥_ 3X5 FD (CM) *y3X9*NO. DROPS* 2X5 *MASS (GM) *)

114 FORMAT (/1X_9* RATE INFORMATION*®/1X,*SYSTEM OTSDTH*s3X_*OMDDTH* 4X,
1 *OMVOTH*s4Xe*OVBOTH*,4Xo*OLBDTH® 44Xs *ONODTH*44X5 *ONODTI* 5 4X_
2 *ONODT2*,4&X, *OMSTTH*44Xs*OMSOT1%94X, FOMSOT2*)

115 FORMATI(/26X,*OROPLET STATISTICS *//* SYSTEM*s3(13537X))

116 FORMAT (//1X,* wro-- INDICATED ENGINE PERFORMACE AND ENGINE EMISS
FION sace- FS 7/720X5FT MoE oP. (PSI) =F
1 19E12.3/20X%,%1.HoPo (HPP =F, E12. 3/20X, FL .S oF Ce (LB/HP-HRI =F E1203
2 /29X 4 *SOOT(MG/M3STP)=*%5512.3/20Xs*NITRIC OXIDE(PPM)=*5E12.3)

117 FORMAT (1X,* RESIDUAL GAS CALCULATION DOES NOT CONVERGE */)

118 FORMAT(/* ERROReeeee MASS CONSERVATION DOES NOT SATISFIEDsesseo*y
1 1P7E10.3/)

119 FORMAT (2X yI5,1P4(2X,2E10.3))

126 FORMAT (/1X_* ERROR IN ITERATION #,1P10E10.3)

999 CALL EXIT

END

62


*FOFCK FYAP
SUBCOUTINE EVAP(P ,T,TAy TL yCO,CMEGAgT IM, RHOS, XMFGy XMAFA,KCUNT,CEVAP
1,AIGNGs XMSOTUD

C THIS SCUTINS CALCULATES THE VAPORIZATION CF FUEL CRCPLETS IN A SFRAY

C ANC ALSG ESTIMATES THE IGNITICN CF EAGH SUBSYSTEM
COMMON/INTEGR/S INO1,INC2,1ISOCT,IWETAT yNFRTyICTAILsNTEReNAXIL,y
1 MSHEL »NSZGeNTTST,IGUGyII,III
CCMMON/EONVES XJgATO OTS XK ySIG SR RUCTOTHyCTEOT »GTROT SPIT
COMMON/SCA/ THEMAX,THEIJ,THETG,OINJs THENBs THEVGe THEVCyDTING,

1 QTHMIN,OTRMEC,OTEMOX sCTHETASOT
COMMCN /SETRES XMOTSXNF do XMEGR yXMRDgXMAG py XMARO y XMAEGRyXMFROyXMFEGRy
1 AMWA,AFRA,FARASPHEA

_ COMMONZFUELP Y's XX eYFLSZELeYSTOAsZSTOAAPIyFCEN sFAROS ps AMWF, SURFT
oXMUF CEN, GXST »HCOM

* conuon JSHEDT/ HFCO2,HFHZCs0HKVH2C,HFHECG ,CFFLSCPFVsCHFV,HHV,T3PF

1 sHF Jy TFEJoHFFL SUPPL

COMMON /SCOTP/ SDENsDSCOT, ASCTP,ASOTA,XMSOCT (100

CCMMONZENGEO/ NEsRFWyRLeVOySCsS USTs CReRPMs TH yONCZZ¢NNOZZ,0BCWLe

1 VFJO

COMMCN/DPEOPRCS HUTTM, ITPHY sXTLYO2,CPREE, PHON, TORITF yXL 2984
XMBM (100) ,O0VQNTH(109) ,SLOCTH(1C0),XMFLI(109)4
OLC100,7),OMFLI (160.2) »OMLCT (109,43) ,XNDT (1000305
XKK0100,3)

COMMCN/JETPY DIF yXSC20) pVS (2094) VS (C2094) 4 2HOF (2044) pAFT(20 94) y
PTIM(20, MV sETAT (2925) sDYS (205404 GOSZ(3),SZBAR (3D,
C259 Gr Bde XES(25 44) pVSS (2544) sCYSS (25,4 a TAFV (25 94) 9
BEV C2554) XMFVS {25 54) sXMSOT (2544) yg TOC25 9 4) POV 25 gb) yg
TOIG(25494)

ane |

Pal

1.6
FR**0,333333
GENBQ= CEN**0.69

METS = Neer
XMUAS L.25E-SF (TAF*L.S/S(TAGIIOL)
XL= XL298FCCTCRITFATLUIS(TCRITE -298.9) ** 0.39
TF (TIM.GT.OTINJ) GO TC 24
C ISSUING 4 NEW SURSYSTEM
TIl= III+1
Oc 22 J=1, NSHEL
TAFVCLIT,J)= TA
AFV(TIT, J =1.0510
XMSCT(LII,J)=0.9
IOIG(IIT,J= .
TOCITIyv= T
POC(IIT y= P
90 23 L=1, NSZG
OULIITyJ,L)= SZBAPIL)
CONTINUE
CONTINUE
CONTINUE
XMFG1= 0.6
XMAFAL= 06.9
MSOTI= 0.0
ec £ Te1,101
TI4T= TIM-(T-1) *0T
OG 2 J=1,NSHEL
IFCINTG(T,J).6£.3) GO TO 2
C LOCATE THF NFHW COCROINATES FCR FACH ELEMENT
Ne 11 TO=2,NAXTL
IF CTIMI GE SETIM{ID-1,5) «AND. TIMIGLT. PTIM(ICs JI) GO TO 12

Nw

ae)

co
+

66


SUMS= SUMS+ AAM (ID *XMBMCTI*CROMC ID +TBM (I) *PRPTBM(TD)
74 CONTINUE
79. CONTINUE
S LIQUID FUEL SYSTEM (VAPORING)
IF (XMFL.«LE.0.9 OR. IWETBST.EQ.0) GO TO 27
C THe TEMPERATURE OF VAPORIZING LIQUID FUEL IS ASSUMED TO BE THE STABLE
C WET-3ULR TEMPERATURE
CALL WETBT(P,TAsAMWA,AMWF A,B,C, TFLyDHFYP)
UFL= UFFL+ CPFL*ITFL-298.)
HFL= UFL+ P/FDEN*ATOXS
DTLOTH= (TFL-TFLOI/OTHETA
€ LIQUID FUEL SYSTEM (BURNING)

TFL3= TFL
UFL3= UFL
IF CICRITA.NE.C) GO TO 27
HFOS= HEL
TFOS= TFL

27 CONTINUE
FUEL SAS SYSTEM
QAFG= 0.0
BPFS= 0.6
ORFOTH=0.5
IF (XMFG.LE.0.0) GO TO 28
VFG=VFBTAFG
CALL RADHZ(TFS, XMFGyVFG,)ORFOTH)
OGROTH= OGROTH+ DRFOTH .
OCFITH= GFG*OGCOTH
HFG= UFG+RFGFTFS
C*****# ASSUMING THE LIQUIG FUEL IS VAPORIZED TO THE LOCALL FUEL GAS TEMP.
HFV= HFG
OQRFOTH+0CF OTH
CHEV -HFG) *OFVOTH- QFG
MFG* (CVFV#RFG)
(QF GFTFG) /PS(CVFVtRFEG)
RBRFG= XX1/XX2
SUML= SUM1+# RFG*TFG*OFGOTH
SUM3= SUMZ+ 3BFG*XMFG*RFG
SUM5= SUMS+ AAFS*XMFG*RFG
28 CONTINUE
4IR SYSTEM
SAA= G.6
BBA= 0.0
ORAOTH=2.0
IF (XMA.LE.G.U) GO TO 26
CALL TPROPIT4yP sVXAy CVA, CPA,UA,HAy AMWAdGAMA,RAyPHEA)
CALL SPECY(P 4TAsPHEA ,AMWA,XH2A yXHZOA yXCON, XCO2ZNgXN2N gXNOAYXOHAY
1 XH&,XOD,XNA,XO2A,XCHGAD
VA=VFETAA
RHOA= XMA/VA
CALL RANHZ(TA, XMF A,VA,ORADTH)
DOROTH= OGROTH+ DRADTH
OCADTH=G4*OQCOTH
QEVOTH= CHFV-UFL)*OF VOTH (UFJ-UFLI*OFJOTH+ XMFL*CPFLFOTLOTH
QS= DRANTHHOCADTHFOEVNTH- (HF J-UFJ)*DFIDTH
(ROFTA/P-9UPP-TAFPRPP)FXMA
CPUP THR AFT AFPRPT) FXMA
XMA* (PUPF ETO FORDE)

oO

ce

x Ged
SNOODTH= 0.0
SMSOTH= 0.0

54


SGIVEN
oc
REY
RL

cr

sc

TM
RPM
scr
ONne7z
NNOZZ
OBOuL
THEYC
THEVO
xmay
EGP
XMRO
ONEGA-
NSTROK
THEI
orns
XMF J
PINS
TFD

"

Xx

FOEN

API
xHUF
SURFT

TOPF

a

ry

a

(2)

O.1397E +925
O.762E+Cty
O.30G8E +02,
O.17E+02y
G.2032E+00,
O.5E+035
O.1S5E 464, -
DetE+Ol,
De2545-01y
be
G.82E+OL,
06155403,
O.145403,
O.24E+01,
Oy
O.ly
0.66403,
-2882303761516078135
-0.2€+02,
Q.17E+02y
Os1Et se,
O.231€+03,
0.35E+035
O.7E+01,
O.12E+02,
600 +
O.4G7E+02,
0. 345E+02,
O.12E-OLy
O27E+02,

O.3776EHO3y

OUTPUT.

23


*
o
m
>
az

GK WAIN
PROGRAM MAIN CINOUT,OUTPUT, TAPES=INPUT, TAPES=OUTPUT, TAPES)
SINGLE-CYLINDERED OIESEL ENGINE ENGINE ANALYSIS
CYCLE THERMODYNAMICS MOOIFIED (12/2/74 5CJK)
CYCLE THERMODYNAMISS MODIFIED (12/5/74,CJK)
EQUILIBRIUM COMPOSITION CALCULATIONS MODIFIED (1/30/75,CJK
VITIATION ADDED AND IGNITION DELAY MODIFIED (2/12/75,CJK)
Q-S THIN FLAME IN FINITE SPACE FOR HET. NO RATE INCLUDED (3/27/75,CJX)
NUMERICAL SCHEME REFINED (4/30/775,CJK)
ENTIRE PROGRAM CLEANEN AND STORED IN A PERMANENT FILE (5/12/755CJK)
VERSION B eseee IMPROVE COMPUTATIONAL EFFICIENCY 3Y 26 PERCENT (6/1/75,CJUK
OIMENSION [00(8)
DIMENSION OT389TH(100) ,ONOOTH(100)+DMSOTH(19C) »DAPDTH(100),5
1 OPHOTH(100) »XFIDTH(105) »XATOTH(100), QBM(100),OMDDOTH(109),
2 OMVDTH{1L59) ,OCMOTH{100) ,RTABM(100),G693M(100),44M(100),B8M(100),5
3 YNO(160)5,YS(100),XMNO(100)
DIMENSION TBMO(100),DTBDTO(1CG6) »XMNOO(166) ,ONODTO(106),XMSOTO(106)
2 +OMSDT9(100),XMFBMO(100),XMPAMO(i00)
3 sXFIDTO(100),XAIDTOCL09) sDIO(10055) »XKKO(10055)
4 9Q8M0 (100), DMO0TO(100),OMVOTOC100),0L89TO(100)
COMMON/INTEGRS INO1,INO2,ISOOT, TNETBT NPRT sIDIATL NTER NAXILs

AMAVQANRANDA

1 NSHEL,NSZG,NITST,I3UG,II,III
SOMMON/SPECYIZ = XH201009 yXH20 (100) »XCO(1GE),XSO2(10 Ode XN2(1E0),
1 ¥NOC400) ,XOH (109) yXH(105),XOC109) »XNC€190) »XO2(10 0)
2 9 XCH4 (100) ,XNONEQ(100)
COMMON/PRATETY ONOOT1(100) ,ONOOTZ(100) ,OMSOT1(100) yOMSOT2(100),
4 BRMOTH(100),0LB0TH(100) »XMOF (100)
COMMON/TPROPIY P,T3M (100) ,USM(100) sHBM(100),V3M(160) ,CPBM(100),
1 AMWBM (100) ,R3M(100) ,PHEC LOC) sAPRCLO0) yXMFEM(100),
2 XMP4M(190)
COMMON /QIRP/S XMATyXMF Jy XMEGRy XMAD y XMAO yg XMARDy XMAEGRy XMFROQXMFESRy
1 AMWASAFRA,FORA,PHEA
COMMON /FUELP/Wy XX y ¥FLy ZFLy VS TOs ZETOA, API, FIEN, FAROS y AMWE y SURFT
1 yXMUF »CENsOXSTsHCOM

COMMON/ORGPC/ HUTIM, DTPHY, XTLY 02 sCPBSAR,EHOD, TORITF yXL298,y
4 XMBM(100) ,OVBOTH(10G),OLDOTH(100),XMFLIC100);
2 OTC190653),OMFLT(100,3) sOMLOT(19), 3) sXNOIT(10693)5
3 XKK (100,53)

COMMON /SHEAT/ HFCO2,HFH20,0HVH20,HFH20S »CPFLyCPFVs,DHFV,HHV,TBPF
4 sHF Jy TESyHFFLQUFFL

COMMON /SOOTPY SDEN,DSOOTs ASOTP ,ASOTAyXMSOOT(199)
COMMON/RANGS/ TPRANG1,TRANG2,PRANG1 yPRANG2s PHERGL» PHERG2

CCMMON /PORVI/ PUPP, PUPT,PRPF, PUPF yPRPET,PRPP

COMMON/PORVIZ/ PUPPIM(100) ,PUPT3M(10G) »PRPFBM(100) ,PUPFEM(1060),5
1 PRPTOMN(100) »PRPPSM (109)

COMMON/GONVF YS XJeATOsOTRyXKySIGyRe RUZ OTOTHs D THOT yOTROTsPI
COMMON/SCA/ THEMAX,THETJS,THEIG ,OINJyTHENB, THEVO,THEVCyOTINGg

1 OTHMIN,OTHMEDsOTHMAX,»OTHETA,OT
COMMON/ENGEO/ DEgRFWyRLy VC 4SC oS, UST yCRyRPM,TWsONOZZ»NNOZZ+OGOWLY
1 VF JO

COMMON/SPTAB/ NPTByNTTB,NPHTB, PT B(7) yTTB(13) sPHETBC110 4
XH2T3(07,13511) »XH20T3(07513511) -XTOT3C751351195
XS0273 (7413511) ,XN2TB(7,13,11) »XNOTB(7413,11)5
XDHT3¢7 913,41) »XHTB( 7915511) »XOTBI7y13511)5
XNT8(7 513,14) ,XO2TEC 7913511) »XCH4TAC7 512511)
NAMELIST /GIVEN/ DEyRFWeRLeCRy SST Wy RPM »SSR+DNOZZ»+NNOZZ,DBOWL,
THE VE, THEVOsXMAL 9EGR »XMRD»OMEGA,NSTROK,
THE TJ sD INS XMFISFINS TR
WoXXsFOENyCEN,API yXMUF,SURFT,
T3PF, TORITF sCPFLy CPF VyCVFVyDHFVyHFFL,HCOOM

anit «

PWN

44


°

SAVEVCITEP)Y= YLCNTTST)
YLCNTITST) = YNO8
HxX= APS(CRK)
ITF (CCAPSCHK ALT 61 .0£-6) eANDSCYLINITST) «GE 005)) WOR.
4 (AISCHK.LT.0.005*YNC3I) SC TO 2
IMPROVING INITIAL GUESS
IF (ILCW.EG.D .AND» CHK.LT.9.0) GC TC 14
TF (TUPS.E9.0 eANT. CHK.GT.9.G) GO TO 15
IF (CEKLT.C.9 AND. CHK.GE.CHKLOW) GC TC 14
TE (CRK.GT.5e0 AND. CHK.LE.CHKUPR) GC TC 15
co Te 16
44 CONTINUS
TLOW= 1
y1oucw= Yitt)
YU8LOH= SAVE2( ITER)
CHKLOWS CHK
GQ TO 17
15 CONTTNUF
TUP?= 4
yigues= vit)
CHKUPR= CHK
17 CONTINUE
TUPLC= ILAW*TUP?
IF (IUPLO.LE.G) GO TO 16
YiCLb= (YLOLCNFCHKUOR= YLCUPRFCHKLOW) / (CRKUPR-CHKLOW)
TF (Y18LOW.LELD.0Y YLCLI= OFF Y1OLOK+ YIOUPE)
59 10 6
15 CONTIBUE
IF (ITEP.ER.1) GO TO 2
ITERL= ITEO-14
TE CAPSUSAVE2 (ITEP) HSAVE2(ITEFLD) eLE. 120&-20) GO TC 2
TE (APS(SAVEL(TTERI-SAVELCITEFI)) oLT. 1.0&-20) GO TC 4
FCOYL= (SAVE 2CITESY-SAVE2CITEFLII/S(SAVELC ITER -SOVELCITERID)
YICLI= YiClie CHK NSO
69 TC 6
4G CONTINUE
WPITE (6,101) ITER, PyAPW,RROO 4yYNC8SXMOCT SRL yRF RB, RMAX
2G 8 T=1,NITST
WPITE (6,199) TyP0F (CT) sTOF CI) sYFCF (I), YO2CF (I), YN2CF (IT) yWOOTA(T)
1 HONT(ID, YNCOF (IT), ¥Z¢19,SAVEL(I),SAVFZC1)
8 CONTINUE
sToe
% CONTINUE
CVER=ALL KC PROSUSTION PER NECPLET
CNCcOTI=0.0
OO 5 Tet,NItst1
O.S*(ROF(T+1) +R0F (TD)
FOF (I +1)-R20F (1)
NTI= ONCCTI+ WOOTAV(T) FAT FRI FORT

PMOITTS 12.566 37062¥ONCOTT

166 FCRMAT(1X%,1591P12E1962)

D1 FORMATIS/1X, %.. FRPOR IN HETNCR ITESATION 26* y I3,1P9E1062)
PETHON
ey

75


I. INTRODUCTION

This document describes the computer program and serves as a user's
manual for a numerical model on THE PREDICTION OF POLLUTANT FORMATION
IN DIESEL ENGINES. The model is capable of predicting both engine
performance and formation of nitric oxide and soot in a direct injected
diesel engine combustion process. The program is sufficiently versatile to
investigate the influence of several engine parameters, such as geometry,
injection timing, compression ratio, load, swirl, etc., and thus can be
considered as a tool readily accessible to the design engineers. The
physical and chemical processes, the detailed mathematical analysis and
the assumptions upon which the model is based are discussed in the full
report on this investigation. This manual contains brief descriptions of the
program subroutines, program input and output. One sample case is also
given to illustrate the use of the program which is designed to facilitate
program check-out.

The computer program is coded in standard FORTRAN IV, and a version
is currently operational in CDC 6600 computer.


FMECK RAQDHZ
SURROUTINE RADHZCTHPs X¥MFSyVI,DOROTH)
G FATE OF RADIATIVE H.T. FROM INDIVIDUAL HOT ZCNE TC THE WALL
C BASED GN LOTOL ZCNE TEMPERATURE ©
COMMON/SCONVFS XJysATOSDTReXKySIGsRsRUSCTOTHSCTHOT,OTROT,PI
COMMCN/ENGEOS DEyRFWsRLy VC SC ySyUST oC Ry RPM, TK »CNOZZyNNOZZ,0B0WL y
1 VF JO
ALFAH=0.9
ALFAS=0.25
XLI=S+SC
AH= 1.576796227*DE*DE
ASLVI= 3.14159*DE*XLI
ALOH= CAH*FALFAHFASLYTI*ALFAS) /(AHFASLVT)
RHOFI=0.0
IF (VI .GT. 1.0E&-10) RPHOFI= XMFA/VI
RCP= 5.752*RHOFT
VZE= VI/NNOZ2
XXX= €.0FVZE/PT
PLEXXX**(0.3323337)
AZE= 2.1G159*OL*PL
EPSA= ALPH*(1.0-EXP(-RCPFPFL))
XXX=(THP¥FGATH FFG) FSTC
QPF= NNOZZFEPSAFAZEF XXX
CCRITH=QRFFOTOTH
RETURN
ENO

78


FRECK HOMNOZ
SUB20UTINE HCMNOZ(CNOR)

© TRIS PCUTINE CALCULATES HOMNEEKECUS NC PRCDUCTICN FATE
COMMON/INTEGRS INO1,INC2,ISOCT,IWETBTsNERTsICIATL »NTERsNAXILy

1 NSHEL »NSZGyNITST»IBUG,IL,III
COMMCN/SSPFLYTZ = XH2(100),XH2C(100)5xXCC(100),xCC20100) »XN2(100),
1 XNOC100) -XCHC100) yXH(100)5xXC (100) ,XxN(1009,x02(100)
2 +XCH4 (100) »¥NCNEQ(100
COMMON/PRATEI/S ONODT1(¢1009,0NCOT2(106),0MSOT1(109),9MSDT2(100),
1 ORMOTH(109) ,OLEOTH(100) ,XMOF (100)
COMMCNSTFOOPTS P4TRM (1099 USM C100) HEMC10G9,VIN(100) sCPAM(1O0CI»
tL AMWEM (100) ,29 (100) ,PHE (100) ,AFR2(100),¥MFEM(100),

XMPAM (109)

COMMON/OPOPC/ HUTIM, DTPHY,XTLYO2,CPBAK,°HON,TCRITF,XL298,
1 XMBM (100), PVEOTH (100) ,OLDDOTH(195) +XMFLI (100),
2 91(10053) ,OMFLI (100, 3) OMLOT(10%43) »XNOTC1005 3)
3 XKK (10043)

oc 1 T=1,It
PMONTIITI= 0.9
TF (TNOL1EQ.0 .OP. TOM(TY.LT.1900. «CR. XO2(T).LE.1-€-15) GO TO 1

= TRM(T)

£Ts 1.927"
© PATF CONSTANTS

XKF 1 B6E1UFEXP(=754906/RT
Je 2EL3 EXP (-4C0./FT)

1.336108T SEXP (-7080./2T)
' =a*T *EXP(-39100./9T)
MKEO3= 5.0*ENP(-59300./RT)
XKL= CNO9*O. 220782 *(XKF1*XKECS)
XK?2 = (X¥KBLFXKRAS(XKFLFXKFE2))
XH %SXKPL/XKE2 .
CORPL= Led XK2*CXNONEG(T) FXNCNEC(T) I / CXN2ET) #XO20T))
CORP2= 1.04 XK3B*(XNONEC(TIZXC2(T))
COR®= CORPL/CORR?
VMCOTL= YKL*XNI CL) *SORT(XCZ(1) #P/T) FCCRR
YNCATI= (320./4MWAM (TD) ®XNCOTL
OMCOTLCTIS CXMRMCTD=XMOF (TD) *YNODTL

RETUON
END

77


FOECK STZE9
OUTINE SIZEDCDELPyRHOAyFCENy SURFT»XMUF yCNCZZySZBAF4GDSZ yNSZG)

C THIS ROUTINE CALCULATES THE INITIAL OROP-STIZE OISTRIAUTICN OF A FUEL JET
DIMENSION GOSZ(€ 5)4SZ8ARt 5) .
CM= 37.5
OEL=0.85

OPXX= DELP*1.01325E6
QMAX= CM*(ONCZZ**.5) #(SUQFT FE, 25) * (OF XX ¥* (00375) ) ¥ (FOENFE (-6125))
1 *(XMUL FQ, 25) *(OHOAF* (-0629)
DMIN= O.01*NMAX
DELD= (CMAX=CMINI/SNSZG
DO 1 T=1,NSZ2¢
Gi= OMIN+(I-1)*DELD
N2= 01+0£L0
XX1= DELFALOGIN1/ (CDMAX-D1))
ITF (O2+LTeDMOX) XX2= NELFALCG(N2/ (DMAX-D2))
GOSZCIV= 065%4(1 .0-ERFCXX1))
IF(O2.LT.OMAX) GOSZ(IIV= O.5* (ERFUXX2)-ERFCXX1))
SZSAR(TI= 0.5% (01402)
IF( 02.GE.OMAX) SZBAR(T)= 0.5% (OMAX+D1)
‘ 4 CONTINUE

SUM=0.0
00 2 IT=1,NS2G
SUM=SUM+#+GDSZ (1)

2 CONTINUE

C NORMALIZATION

OO 2 [=41,NSZG
GOSZ(TI= GOSZ(ID/SUM

3 CONTINUE
RETURN
ENO

80


ve

THETEINEG CAPS ~11400 PEATMD® &630E G1 VICHID= 1.72AF 102 os st sue eos VICAL/OHD= 726708001 PHEES 3.221
DEITHCATH/CODe 262966490 WOCCALSGADS-5 U18E 400 NFIDTHCCH/CAD 03 SHNOTHIGH/GA)® 6.666F-03 SNVOTHIGH/CA)= 8.606E-C4
QICOTHCCML/CAIs Z.cCEMEOIO OCPOTHICAL/CAD= 2,195E-91 SVENTHUGH/Ga Ie aye secs SODDTHIGH/CAI= 1.420E-03 SLODTHIG4/CAD= 7.162E-G3
SUPCLECOLDS 2672FE 461 SUMFIIGH)s 2.998F=02 SFIDTHIGH/CA= 6.G43E-03 WKNETICAL)=-2.821£¢02 RRICH3-ATM/GN-K)= 2.857006
UCHERICH/SECIS 3,498E403 OVOTHICMI/CAD 2-4 .845E0G0 OELEIZ= 6.577t-04 TTERP 2
KMUCE (GMI 4.700E-0% YNOE= 6.846E-05 XNOE(PPHI= 6,554F 401 SOOTEIGHI= L.174E-05

a

THETACDES.CADS -10,00  PEATHDS S.OLTESIL VECHSI= 1.HO1F002 TCKES 1.195E003 ULCALZGMI= 7.703E401 PHEE= 3.5082-c1
MPOTHCATHZCAD= 2eQG3E4S3 WOICOL/CAD2-5.397E 499 MF IOTHIGH/CAD= 5.8622-03 SMODTHIGH/CAD® 1.123E-02 SMVOTHIGH/CA)= 1.213E-63
OVNDTHICEL/CAIS 2,205£400 OORDTHIGAL/GAD= 2,.879E-01 SVUNTHIGM/CA)= 2,351E-03 SOBDYHIGH/CAD= 2.471£-03 SLOOTHIGM/CAI= 9.GL4E-03
GLCGELIS L.BOICHCA SURF UCGHD= Z2.G716-C2 SFBOTHIGH/EAI= G4.B27E-03 HKNETICALI=-2.675E402 RRICHS-ATH/GN-KI= 2.d61€600
UCHAR(CH/SECI=® 3.629E403 OVOTHIGM3/CAD =~4.412E400 DFLEI2Z= 4.72GE-D4 TTERPC=

HANOE(GMDS 34323004 YNOE2 1.336F-04 XNOE(PPMD= 1.2776 +02 SOOTEIGH)= 2.255E-05

THETS(DES.CAds -9.00 PEATHDS S.3L7EGd VICH3I= 2.6H0FF02 TEKIs 1.222E003 UICAL/GHI= 7.752E+01 PHEE= 3,910c-01
OPUTHEGTHFGADS 2o5L3E6G2 WPCCAL/CADS-5.226€ 450 OFJOTHIGH/CA)= 4.6620-93 SHODTHIGH/CANs 1.293E-02 SHVOTHIGH/CAI= 1.293€-03
CICCTHICOL/CAIS 2.47RE4CG DOROTHICAL/GAD= 3.442E-01 SVONTHIGH/CAI= 3.794E-04 SOBOTHIGM/GA)= 2.866€-03 SLODTH(GM/CAD= 5.503€-03
SUMILICALDS 2.266F4G1 SUMFACGHIS 2,795£-G2 SFANTH(GH/CAI= 3,245€-03 RKNET(CALD=-2.927E+92 RRICHI-ATH/GM-K)= 2.864E000
USMARICH/SECDS 3.755E493 DVDTHICHI/CA) =-2.977E405 OLLEIZ= 1.088E-33 ITERPC= 2

XMNCEIGMIS SQA2IE-G4 NOES 2.33H0-24 XANOE(PPMD= 2.251E 402 SOOTECGMD= 3.754E-05

THETACDEG.CAYS BLOC PEATHDS S.623EHUL VIGCHSI= 266020402 TIKI= 1.2600403 ULCAL/GH)= 7.786EF01 PHEE= &.252e-01
COITHIATHZCADS 2.G7CEVSS WOCCAL/CADS-4.968E FLO DFIDTHIGH/CAD= 5.AS2E-03 SMODTHIGN/CADS 1.482E-02 SMVOTHIGH/GAI= 1.482E-03
OCEDIMICAL/EAI= 2.H1BE6DD OARDTHICAL/EAI> $.921E-01 SVADTHIGH/CAD= 3.901E-0% SNOOTHIGH/CAD= 3.185E-23 SLONTHIGH/CAI= 5.492E-03
SYMLL(CALYS 2.567EOCL SUMFAIGHDS 3,153C+02 SFADTHIGH/CAD= 3,575b-03 WKMETICALD=-2.977E002 RRICM3-ATM/GN-K)= 2.057E 000
UCHAOIEM/SECDS T.R76ECG3 OVOTHICMI/CAD =-3,54CEH00 OFLELZ= 1-.6uKE-93 TTERPCH 2

MMNGEIGHDS 8.695F-04 YNDE= 3.6u9E-O4 XNOZ(PPM)> 3.299402 SOOTEIGH= 5.518E-05

VHETLIGEG CADE 7.00 PEAT) = S.930E91 VICM3D= 145598002 TEKD= 1,2976653 UICAL/GHD® 7.799L401 PHEES 4,595: -C1
OPOTHCATH/CAdS 2.478EFIG WPICAL/CAD 2-4 6352450 OFINTHIGN/CAD= 5.8622-G3 SMODTHIGH/CAI® 1,682E-92 SMVOTHIGH/CAD= 1,652E-03
OICITHICAL/CAS= 2.7H0F +90 OOPOTHICAL/CAD= &.43ZF-01 SVADTHIGM/CAD= 4.006E-04 SOMDTHIGH/C: 3.305E-03 SLODTHIGH/CAD= 5.442E-03
SUMMACALY= 2.947E693 SUMFUIGMD= §.523L-02 SFADTHIGH/EAD= 3.706£-03 WKNETICAL)=-3.0246902 RRICHS-ATH/GN-K)= 2.8706 400
UCHARICH/SECI= 369945403 OVOTHICH3/CAD=-3,101E490 OELEIZ= 2.4920-03 ITERPCS 2

XHNIE(GHY= 1.179E-G3 YNOF= W.7IAE-0% XNOECPPHD= 4.495E002 SCOTEIGHS= 7.510E~65

THETA(TEG.CAD= 6.00 = PEATHD= G.237F 401 VICH3D= 165460002 TCKI= 1. 3S6E403 ULCAL/GHI= 7.792E+C1 PHEE= 4.937
GPIENCATI/COD= 2.458F 409 WPECAL/CA)=-6 2346660 OFJDTHIGH/CA)= 5.8822-03 SMODTHIGH/CAD= 1,090E-G2 SHVDTHIGH/CA)= 1489CE-C3
OICOTHICAL/COI= 2.960E400 DAROTHICAL/CAD= 4.969E-01 SVUDTHIGH/GAD= 4 .1U8E-0% SODDTHIGH/CAD= 3.50GE-03 SLONTHIGH/CAD= 5.472E-03
SUMOLICALDS S.227E4CL SUMFACGHD= 3.9150-62 SFIOTHIGH/CAD= 3,915€-G3 HKNET(CALI=-3.066E402 RRICMS-ATN/GN-K)= 2.873E690
UCMEKICM/SECDS &S07E903 OVDTHICM3/C4) =-2,6616000 OFLEIZ= 3.270E-93 ITERPI
XANCE (GHD? 1.532F-03 YNOLS G.123F-04 KUOC(PPMDs 5.827F402 SOOTEIGMD® 9.702E-05

THe TACIEG-C8 75.00 PUATMD= G.535E031 VACMSD= 1.5266002 TEKd= 1.3746403 ULCAL/GMD= 7.751E001 PHEES 5.27901
OMDIHIATH/CADS 2,61KE460 MPICEL/CAIZ-3.764E4U0 OFIOTH(GH/CAI= 5.882E-I3 SMOOTHIGH/CAD® 2.165E-02 SHVOTHIGH/C! 2.105€-c3
THOCOLACAIS S.0S9F1TC DGRDTNICAL/CAD> 5,5486~B1 SVADTHIGH/CAD= 4. 207£-06 SOUDTHIGH/CAI= 3,653£-03 SLODTHIGH/CAD= 5.462E-93
SLUCSLIS F,G9GLACL RUMFRIGMD= 4. 322E-02 SFOOTHIGH/CAD= 4.072E-03 WKNET(CALDS-3.105E402 RRICHI-ATM/GH-KI= 248766090
UCHERUCH/SECIE Ge2iKEOO3 OVATHICHS/CA)=-2,219€400 NELEIZ= 3.674i-03 ITERPCS 2

XANCEIGHI= 1.6866F-03 YNOC= 7.444F-04 XNOECPPMI= 7.076002 SOOTI(GHD= 1.206F-64



SYMBOL DESCRIPTION

ITERPC No. of iterations per step

XMNG@E Over-all mass of NO (gm) per power stroke

YNG@E Mass averaged NC mass fraction

XNOE NC concentration (ppm)

S@OTE Over-all mass of soot (gm) emission per power stroke

(4) Engine Performance and Exhaust Pollutants

The computed engine performance parameters and pollutants in the exhaust
gas are printed at the end of engine cycle (i.e. 4 = THEV@) or at the end of
integration (i.e. 6@= THEMAX). They are:

SYMBOL DESCRIPTION

ILM.E.P. Indicated mean effective pressure (psi)

L.H.P. Indicated engine horse-power (hp)

I.S.F.C. Indicated specific fuel consumption (lb /hp-hr)

S2OT Exhaust soot concentration (mg/m? at standard conditions)
NITRIC Exhaust nitric oxide concentration (ppm)

@XIDE

20


ScPARAM
INOL
INO2
Tsoor
Te8uG
NTER
NPRT
NAXIL
NSHEL
NSZG
Icrope
TRATEP
Tspcyp

TOTAIL
NITST
OTHMIN
OTHMED
OTHMAX
THEMAX
EPsIt
EpsIp

SEND -

ls
te
1s
a,
SQ,
2s
209
Ye
3,
Os
O
Qs .

Oy

SO,
O.1e+dt,
O.25E+Ct,
O.SEtO1,
0,35+02,
G.2E-02,
0.2E-02,

26


(x4

4CIVER DE=13.67 ,PEWH7 OF? PLATE My CPR17 Dy SLHO 2032 TH=S ONO eSFME15COay

CPCHL=B.
CINS=47.
APT

wereee Aye

ECERR AY
TEPAC DY

2eSCR=1.UTHEV
oXMFJ=00100, TF

224.4, TRETJ= 26

€GIVEN 7.

ERULS, TOR TTF ASF Toe
CKFV=92,9ES HEF LE -2FU, 790, FING

DEGFIOC, MF. SHTPL (CASE

130.060, THEVO=140. eXMAT=2.4009,EGR=0e0,THET I=
S5% eg WH7 Oy XXH12 60, CEN=4G67,5URFTH27 ag XMUFE1L7F H-2,
TREE = 177.6, CPOFL=0.51 CPFV=E9.232,CVFVA0.22,

Ley DNCZ7=C.GZ254,RNCZ7=F

200 ) tavee

20,

sCVFGA=A00.¢


se

TMETACDEG,COd= 4.92 PEATMD= HAZTEGGL VACMSD= 1.4360L0u2 TEK)= 1,613£003 ULCAL/GMD= 7.687E001 PHEES 5. 621€-01
OPNTH(ATHZCADS 2eSAPE CGD WPACAL/CAD=~3.223E 000 GFIOTH(GH/CAD= S.ARZE-03 SHODIHIGH/CAD® 2.324E-02 SHVOTHIGH/CAD= 2.324E-C3
CACOTHICRL/CADS 3.175E400 DORDTHICAL/CA)= H.0BOF-01 SVHUTHIGH/CAD= 4.300E-04 SMEDTHIGH/CA)= 3.743E-03 SLODTHIGH/CA)2 5.452E-03
SUMIL ICAL = 3.964EO61 SUMFZ(GMI= &.7395-02 SFBOTHIGH/CA= &.173E-03 WKNETICALI=-3.136E002 RRICMS-ATH/GH-K}= 2.879F 090
UCHCKISHSSECD= Ge I25E003 OVOTHICM3/C4) =-3.777£+00 OFLE12= 4.578F-03 ITERPCS 2

MMMOS(GH)= 2.271E-03 YNOE= 9.C4BE~GH XNOE(PPHDS B.591E+U2 SOOTEIGHD= 1.444E-06

THETEEDEG CAD GO PEATHD= 7.197E%u1 VICH3)= 1.GAGEHO2 TIKD> 1.4S1E403 UICAL/GHD= 7.587E001 PHEE= 5.963E-01
OPOTHCATHACON= 2.706F4II WPICAL/CAD=-2.015E490 OFIOTHIGH/CAD= 5,662E-93 SMNDIN(GH/CADE '2.54KE*02 SMVOTHIGH/CAD= 2.54KE-C3
DISIIMICAL/EDIS F.30HEOIG MOSOTHICAL/CAI= 6.652E-D1 SVBDIHIGH/CA)> 4. 388£-D4 SOBDTHIGH/CAD= 3.863E-D3 SLODIHIGH/CAD= S.444£-03
SUMCLICALD= 4.3626 001 SUMFAIGMD= 5.1695-02 SFODTHIGH/CAD= 4. 3u1C-03 NKNET(CALD=-3.162€¢02 RRICH3-ATH/GH-K)= 2.882E000
UCHARUCAZSECD= Geb29ECOT DVOTHIGHI/EAD=-1,333E005 OFLELZ= 5.170%-03 ITERPC= 2 .

HMNGOGH= 2,702E+93 YNOC= B.C74E-C$ XNOE(PPMI= 1.019£493 SOOTECGMD= 1,702E-04

THETAUNEG.CADS 2.60 — PEATHD = 7.3778631 VICMID= LehGGECO2 TIKI= 1.490E403 ULCAL/GMD= 7,535E001 PHEES 5. 963E-01
OPUTHEATH/CAD= 1,6996E499 WPCCAL/EAD==2 9416400 NFINTHEGHZCAD= 0. SHODTH(GH/CA) = 2.767E-02 SMVOTHIGNH/CA= 2.767E-63
DICOTHICALZEAD= 36433E400 OOROTHICAL/E A= 7.234k-OL SVEOTHIGH/CAI= 4.4 72E-0% SOUDTHIGH/CAD= 4.O73E~03 SLOOTH(GH/CADS 5.435E-03
SUMCLICALD= G.777LOCL SUEFICGHD= S.621c-02 SFANTHIGM/CAD= &,521t-03 HKNETICALD=-3.181E002 RRICHS-ATM/GN-Kb= 2.086£+00
UTHAIEICH/SECIE GGILECCS OVOTHICH3/CA) =-8.892E-G1 OFLELZ= 6,031f-93 TTLRPC= 2

MMMOE (GMD = 3,209€-93 YNOES 1.230E-03 XMOECPPHD= 1.166603 SOOTE(GM)= 1.968E-04

THETACIEG.CADS «4400 © POKATMD = 766222051 VICHII= 1.452642 TIKD= 165246403 UMCAL/GMD= 7.513E¢0L PHEE= 5, 963-31
DOOTHUATA/CAD= Le72CECG0 WPCCAL/CAD=-1.295E459 MFIJOTHIGH/CAD= 0. SMHONTHIGH/CAD= 2.923E-G2 SMVOTHIGH/C4)= 2.923E-63
DICCTHCCAL CANS 26555E€4I0 IAROTHICAL/CAD= 7.H38E-GL SVGUTHIGH/CADS 4.551E-04 SOQUTHISH/CAD= 4.1G3E-03 SLOOTHIGH/CAI= 5.427E-03
SUMCLICALIS S.2116661 SUMFIIGH)= 6.I77E-02 SFUDTHIGH/CAD= &.558L-03 WKNCTICALD=-3,293E402 RRICH3-ATH/GN-KI= 2.39CE +00
UCHORICH/SECH= 4.53TEEGS OVOTHICHI/CAN=-G.447E-01 OLLELZ= 8.554E-05 ITERPC= 1

X4NGTIGMDS J.47LE-93 VHOES 1.375F-03 XNOELPPMD= 1,3C1£403 SOOTE(GMI= 2.2218-04

THETACNES COE 1.66] PUATMES 7.84GEOSL VICHSI= 1.462F4G2 TCKI= 1,5670003 UCCAL/GHS= 7,972E001 PHEE= 5.963C-01
OPSTHCATUZCOD= G.E29FeG1 WOLCAL/CAY= &.225E~C1 OFJOTHIGH/CAD= Cy SMMDIM(SH/CA)= 3.070E-G2 SMVOTHIGH/C4)> 2.370E-C3
Casa tetcaL say 7266420 ODOKOTHICAL/CAN= d.413E-01 SVBOTH(GH/CAI= 0. SOMOTH(GH/CAI= 2.790€-03 SLOOFHIGY/CAD= a.
SUNOLICALD= Ge11GEOG1 SUMFAIGN)= 6.685£-G2 SFODIH(GH/CAD= 2.79Gh-d3 WKNETICAL)=-3,193£402 RRICHI-ATH/GH-KI= 2.894E600
UCHEDICH/SECI= G.530E0CS OVOTHICMI/CAD= G4K7E-O1 OFLEIZ= 6.30GL-03 ITERPC= 1

YMNCE(G4)= Get4IE-93 YNGE= 1.637F-03 XNOECPPHD= 1,546F403 SOUTE(GMD= 2.622E-04

G.Chre 3.60 PEATMD = 7.915€401 VICM3)= 1.48GEFd2 TKKI= 1.5970403. UICAL/GMI= 6.540E401 PHEE= 5, 9032-51

CHOPACATMs/CAD= 2oPOLE@GL WRECAL/CAP= 2.1295 000 NFINIMIGH/CAI= 0, SMOOTH(GH/CAI= 3.168E-02 SMVOTHIGH/CA)= 3,168E-C3
OTCUTHIGAL/CMNS 3.803F OIG OOKOTHIGAL/CAI= 8,607E-01 SVANTHIGH/CAD= 0. SOBDIH(SH/CA)= 2.377E-03° SLOOTHIGN/CAI= co
SUMGLICALDS 7.366L091 SUMFAIGH= 7.174 2-c2 SFNOTH(GH/CAD= 2.3776-93 HKNET(CAL)=~3,159E002 RRICH3-ATH/GH-KI= 2.8956 000
UCHANICM/SEC DS GeSTEF EIS OVOTHICHS/CAD= 1.333E400 ATLELZ= S.HuhE-03 TTERPCE 1

YANGEIGAD= WH25E-93 YNOL= £.A25E-03 XNOF(PPMD= 1.725E4C3 SOOTEIGMD= 2.673F-04

THF TO(DEG 6 cade GeG2  PEATHDS 769017491 VICHSD= 2.496LeU2 TKD= 1.6098 903 UCCAL/GM)= 622186061 PHEE= 5.903E-91
OPIVHUOTM/CAIE@2,6239F-G1 WPICOLZCOY= 2,979E900 NFINIMIGM/CAD= 0. SMQDIMIGH/CA)= $,199£-02 S#VOTHIGN/CAD= 3,199F-03
CICTIMICAL/CAL= 3.A23E09G AGROTHICAL/CAl= 4.719€-01 SVBLTNIGH/CAI= Oy SOBOIHIGH/CAI= 2.234E-03 SLDOTHIGH/CAI= 2.

SUMCLICALI= 7.515E461 SUMF3IGMH)= 7,398£-02 SFADIH(GM/CA)= 2.254£-03 WKNET(CALY=-3.129F+02 RRICM3~ATM/GH-K)= 2.895E 400
UCHBE(SH/SECD= GeSSLEOCS OVOTHICHI/CAD= 1.777E400 DELELZ= 4,877E-03 ITERPC= 1
XMNOEIGHD= &.B1Z2E~U3 YNOE= 1.897E-03 XNOE(PPH)= 1.790E403 SOUTEIGM)= 2.964E-%



below.

format.
under:

IV. DESCRIPTION OF PROGRAM OUTPUT

The regular printed output of the computer program is described
The description is presented in the order of its appearance.

(1) Printout of Program Input

The program input group (III) is printed out under namelist output
Descriptions of these variables has been given in Section 3

NAMELIST GIVEN
I.D. Card
NAMELIST CPARAM
NAMELIST EPARAM
NAMELIST GIVEN

(2) Initial Conditions

The printout of namelist INITL gives the initial conditions of the engine

and the charge in it. The elements of NAMELIST INITL are described below.

SYMBOL

THETAI
THEMAX

UA
HA
XMA@

DESCRIPTION

Initial C.A. for integration (or compression) _

Maximum or final C.A.; the integration is terminated when
6 = THEMAX

Initial pressure (atm) .
Initial charge temperature (Px)

Initial charge volume (om)

Specific internal energy of charge (cal/gm)
Specific enthalpy of charge (cal/gm)

Total mass of initial charge

7


FDFCK TPROP

SUBROUTINE TEROPCTX PX yVyCVe COU sH
Cxevee esse THIS RCUTINE CALCULATES --
Cerreeeeeee EQUILIBRIUM THERMODYNAMICS OF Ccmousticr PRODUCTS
Ceeeeseeeeee U(P,T,PHE) AND R(P,T,PHE) AND THEIS OERIVATIVES
” COMMCN/FUELP/ Hs XXoYFLyZ2FLyYSTCAsZSTOAZAPI,FLEN»FAROSs AMWFE,SURFT

1 + XMUF yCEN,OXST,HCOM

CCMMCN/RANGS/ TRANG1, TRANG2,PRANGI »PRANGZsPHERG1,PHEFG2

COMMCN /PORVIZ PUPF,PUPT PROF ,PUPF »PRFET .PRPP

COMMCN /HPCOE/S AA1(6),B881(E) ,f01011) »RRICK), SA2CE) ,BE206) fO2C11)

SHW) GAR oR ePHEX)

1 yRR2(8)
Pi= PX

Ti= TX

PHE= PHEX

IF (P¥.GT.PRANG2) Pi= PRANG2

IF (TX.GT.TRANG2) T1= TRANG2

IF (PRKEX.GT.PHERG2) PHE= PHERG2
TIS1000./TL

Te= T1*T1
Tis T2*ti
Ta= T3*T1
TS= Tutti

PHE? = PHE*PHE

PHE3= PHE*PHE2

QALCOGP= ALOG(F1)

IF (FPRE.GT.1.0) GO TO 2

1 CONTINUE
C FUEL LEAN

A= AALCLVFAALCZVFTLFAAL(AV¥FT2 FAALIG) FTI FAAL(SVFETH +AAL(69FTS
B= 381(1)#8B102)*T14881(03)*T2 #BB1CGI*TI +BR1(5I*TY #28116) *T5
ots COL(1+OCL(2) FPHE+ ON1(3) FPHES

OP= (£01(4)+001(15)*PHF* COL(E) #PHES)D FTI

G30=¢(001(714C04(8) #PHES DAS) PRES) + {CO1(10)+001011) *PHE) FTT)
DF= O307ALOGP

XXX=91#D2403

IF (XXX .LTe-2006) XXX=-200.

O= EXPCXXX)

UX= A- BFPHEt D

XXX= CRELC3P+ RRLCGP*TI# RR105) #ALOGPIFPHE+ RRAC6I+ RRIC7I*TI+
1RR1(8)*ALCGP

IF (XXXeLT.-200.)9 XxXX=-200.

C= EXPOXXX)

SX= RRLCLD+ SRL(2)*PHE + C6

HXSUX#RXFTL

AP= AAL(2)42.0FAAL(3)*T1+ 3.07 AAL(H)FT2+ 4 .OFAAL(SIFTI+ S.OFAAL(E)

1 *TG
BPs BEL(2)+2.0*BRL(SP*TIt SCF ESL (GIFT Z+ GG*ERI(S)*TS +5. 0¥AB1(6)
1 eT4

O2P= -(001(4)+ DDL (5) *PHE+ NC1(6) FPHEI) F100C./T2
O3P= -(001(10)+ DO1¢11) #PHE) #1000.74T2* ALOGF
oP = C*(92P+03P)
CVX= AP-OPFPHE+0P
CPX=CVX4#RX
CONVE?T UNIT TG PER MASS CF MIXTURE
XBASE= 1.0+FAROS*PHE
Fi= 1.0/X@ASE

Uz UXFFL
H= HX*FL
R=RX*F 4

AMH= 1.98779
92...


FOFCK TEOHZZ
SUPROUTINE RADHZZ
© POTS GCF RANTATIVE H.T. FROM INOIVIQUAL HCT ZCNE TC THE WALL
© SASEN CN LOCAL ZONE TEMPERATUFE
COMMON/INTEGRZ INONL,INC2,ISOCT,IWETSTNPRT,ICTATL yNTER NAXIL,
1 NSHEL»NSZG,NITST,IGUGyII,IIT
CCMMONSCONVES XJsATI OTR aXKy SIGs PeRUsCIOTH sO PHOT, CTROTsPT
COMMCNJPPATETY ONODTI(1NC) »INGNT2(1L00),ONSOT1(100),OmMSDT2(1 00)»
t ORMOTH(100) ,OLGOTH(100),XMODF (109
COMMCR/TPROSTZ P,TEM (100) ,USM (100) .FRM (100) ,VEM(100),CPEM(19C),

t AMWBM (109) 5RAMC100) sPREC1CO)sAFROLOC) yXMFAM( 100),
2 . XMPAM (100)

COMMCMZENGEC/S DE RFWyPLyVC SC oS oUST Ch y PF My TK yONOZZyNNOZZyDECKL,
1 VFJO

THES TW¥*L

ALFAH=0.9

ALFAS=0.25

XLT=SS+S¢

1.570795227*0E*0E
ASLY F.LGLS9F#OFFXLT
ALPH= (GH*ALFAHFASLVIFALFAS) /(AHFASLVTD
co 1 I= 1,1!
SuCrI=0.9
IF (VQM(T) .GT. 1.G&-10) RHOFI= XMFBM(I)/VOM(T)
RCP= S.7FA*RKCFL
V7E= YVRMCT) /NNOZZ
XMS ES TEVZE/SPT
PL=YX¥**(0.3333323)
AZE= 261415EFPL*PL
FPSA= ALPH*(1.0-EXP(-PCPFFL))
XYX= CTAM CT) FF GHTHG) F SIG
GPE= NNOZZ*EPSA*AZEFXXX
SEMOTHCT)= CREFCTOTH
1 CONTIMUF
RETUPR
FRO

79


#NECK BOATA
RLOCK DATA
FOLYNGCMIAL COEFFICIENTS FOR THERMODYNAMICS CATA OF CURDH(1.7N) COMBUSTI
PRODUCTS UCPyTsPHE) IN CAL/EM-OF AIR ANC R(PyTyPHE) IN CAL/GM=OF -AL2~K

PIS TN AT, T IS IN Ky PHE IS EQUIVALENCE RATIO
CrRRTY 470=1.71

COMMON /HFCOE/S AAL(6),B810€),COL(11),kR1(8), AAZ(E), E8206), 002011)

oPR208)

Qa0

1:
© LEAN PRODUCTS
DATA AAL
1 / ~S.S5G49233EE+01, 7e83I7475RE-02, 1.6734745E-04, -1.1376329E-07,4
2 GO737BLGEM115 ~SeGTBIB4GEW15 /
DATA BBL
L/S 7e3BG65276E+02, FW.7BL747EE-O2, -1.0774E0ZE-9%, 6.1112624E-08,
2 T1.G435962E-11, 9.9489198E-16 /
DATA OD1

1 / 1.37CH2C2E+O1, -1.,1017790E +00, -S.7155SC9E-01, -2.67R04ESE+01,
2 7.2G48177E+00, 143090986E+00, S.O189712E-01, -S5.2147343E-01,5
z 5 .3444392E-02, -2.6396E27F +00, 1,.60€8921E+00

DATA PRL

1 4 6.8789266E-92, 4.2293620F-02, -S.83612E8E-01, 6.6435220E +00,
2 1.1172543E-01, 2.3055596E+00, ~2.5753610F #01, -4.0113958E-O17

C RICH FROOLCTS

DATA AA2

1 / -R.SSSOSISE+O2, 2.5318868E-G1, -1.61766C4E-05, -1.0563951E-07,4
2 7.2508183E-11, ~1.2100101E-14 7

DATA eR2

Lf ste 79EUKTE+O2, 6.598E551F-01, -1.16E2944E-03, 7.1690345E-07,
2 TLB777U4TER10, 1.79C5017E-14 /

DATA No?

1 4 Ga7714237E400, Gs3547107E 400, -7.25€01Z4E-02, -9.6825140E 400,
2 TLS7H7997E401, 2.1 2GINFE-01, -E.U44S7ECEHO0, 7.2392234E400,
3 5.2083600F-03, 2.2071211E+01, -2.5775327E+01 7%

DNATA PR2

1 4 4.324E505FE-02, 1.250CE87E-02, 2.7144593E-01, -3.8400357E-01,
2 -3.,3R507845-03, ~4.0003995E4+00, 2.9256778E-O1, -1-79734EFE-02/
END

64


INPUT
Je Engine description
@ Intake air, fuel and CCR
properti
Je Contret parameters
je _Empiricat parameters
INITIALIZATION
© Physical constants &
parameters

© Residual gas

Properties of initial cylinder
content

© Set zeros

VOLUME
¢ Cylinder volume (ELM;

J

-EUEL SPRAY

Dropsize distribution (SIZED)
Alr/apray entrainment (SPRAY)
Evaporator & ignition (EVAP)

Wet bulb temperature (WETBT)

Yes

© Alr/product mixing
(amt

EQUILIBRIUM COMPUTATION
1 themodynamte properties (TPR?
@_Species distribution (SPECY2

DIFEUSION FLAME
@ Droplet burning rate (FLAM2)
@_Hetergeneous NO rate (HICTNGR)

Co

Figure 1

HEAT TANSr te
© Convection (CONHT)
© Radiation (RADIIZD

HEMOGENCOUS NO RATE
2 NC kinetics (HOMNCZ)

SCOT FORMATION.
© Soot kinetics (S@@TRZ)

rary

oar

oe

CYCLE THERMCDYNAMICS
Integrating energy equation
Mass conservation equations
NC and soot conservation equations
Determining Zonal
Temperature and a/F
NO and soot
Detine Overna
Temperature
Pressure
Mass of NO and soct

QUIPUT
© Overall properties
© Zonal detsiis

© Engine performance
© Dxhaust pollutants



ze

TRETAIVEG.CAIS =75.00 PEATHD= S.101E400 VEGHSD= 161508003 TIKI= S.O9MECUZ UICAL/GMD= 4.548E400 PHEEZ 1.458E-02
OPUTHIATY/CAIS Pel31F 902 AOUCM/GADE-1.493E600 DFIOTHIGH/CA)2 Oe SHODTHIGH/CAD® O. SNVDTHIGN/CAD= 0,
DICHUM (CLEC ADS-1.755E902 OROTACCALJC ADs 3,082E-07 SVPOTHIGH/CAI= QO. SOMDIH(GH/CAD® 0, SLODTHIGN/CAD=
SUMCLICALDS-1.774E4C1 SUMFIIGHI= O. SFBOTHIGM/CAD= Oy WKNETI(CALD=-5.293E¢01 RRICHS“ATH/GH-KI= 2.840E400
UCHAVITM/SECIS 2e43LE4GS OVOTHICH3/CAI =-2.100E02 DELEL2> 7.619E-64 ITERPC=

XYNCEIGMD® G.4G3E“UA YNOF= 2.6150-4A XNOS(PPMD= 2,516€-42 SOOTLIGHD= 1.582E-22

THTTAODES GAPS -65,09 — PEATHD= GO79E 400 VECHID= D.419F002 TEKI= SeG92E0U2 UMCAL/GMD® L.285E4O1 PHEES 1. 468E-02
OPOTHEATA/CAIZ L6C4-E-GL WRICAL/CAD=-2.902E400 OFJOTHIGN/CAIS Of SHOOTHIGH/CAl = 0. . SMVOTHIGN/G4)= 0,
QACOTHICAL/CAD= &.7E7E-02 MOR OTN (CAL/CAD = Ls7T14E-05 SVBOTHIGH/CAI= 0. SOBOTHIGH/CAI= Ge SLOOTHIGH/CAD=
SUMGLIROLDe-$.733F 401 SURFALGMD =

SFAOTHIGM/CA)= 0, WKNETICALD=-7,082E¢01 RRACMS-ATM/GH-KI= 2.640E000
UCHAEICH/SECDS 26427E003 OVNTHEC Sica e-2.ounevod OELELZ= 1.298b-33 ITERPCE 2

XMNOF(GHI= 1.999E-G3 YNOEZ 4. 461F=44 XNOE(PPMI= 4.290E-36 SOOTEIGMD= 2.8566-21

THETACQEG.COD= -55.00 PAAIHD= 5.619600 VICHSI= 726316402 TEKI= 5.968£402 UICAL/GH)= 2.072E+01 PHEE= 1. 486-62
OMVTHCLTM/CAI= LePIZE-GL WPICAL/CAD=-2 448E 400 MFINTHIGH/CAD= 0. SHOOTHIGH/CAD= 0. SMVNTHIGH/CA) = 0.
OISUTHICAL/CAD® 14729E 01 ATNTHICAL/CAD= &.566E-05 SVBOTHIGH/CAI= 0. SORDTHIGH/CAD= 0, SLODTHIGH/CAD= Ge
SUMOLICOLD=-1.54N6461 SUMFOIGHD= 0. SFADTHIGH/CAI= 0. WKNETI(CALD=-9,385E¢01 RRICHS~ATH/GH-K)= 2.840E000

UCHSCICH/SECI= 26400E403 OVOTHICMI/CAD =1.914E401 OFLEL2Z= 1,.565£-93 ITERPCS 2
MAMIEUGHES 2.299639 YNOFS 8B, 930E-60 XNOELPPMD= 8,568F-3% SODTEIGHD= 5,668F-20

THETAUGEG.CADS -65,00 . PUATHD= A6160C000 VECHSI= 5.6170002 TEK)= 6.5516002 ULCAL/GMD= 3.161E¢01 PHEES 1, 4682-92

DERTHCATH/CAIS 2e899E-01 WO(CAL/CAD=-3,158E400 OFIDTHIGH/CAD= 0. SHODTH(GH/CA) = 0, SMVOTHIGH/CA)= 0.
DIG THICAL/CAD® 3.112E-21 MRATHICAL/CADS 9.5u2E-0% SVNOIHIGH/CAI= Oo SOBDTHIGH/CAD= 0, SLOOTHIGH/CAD= Ge
SUMPLICSLD=-1.514E 421 SUMFA(GMI= 0. SFOOTHIGH/CAD= 0. WKNET(CAL)=-1,235E¢02 RRICHS-ATH/GH-KI= 2.600E400

UCHEE(CHSSECD= 2463526453 OVOTHIGH3/CAD =-1.70CE401 OLLELZ= 1,914E-93 ITERPC= 2
MMICE(GAD= 5.222E-35 YNOES 2.323£-35 MNOCIPPMD= 2,2340-29 SOOTLIGHD= 1.205€-18

THETAGDES COPS -35,09 PATH= L.2GTF OL VICHSD= 4.G55E4G2 TIKI= 76262E402 UICAL/GHD= 4.565001 PHEES 1, 4582-32
CRDTHIATMsCANS 5.608F=01 WPICAL/CQ)=-6.0410000 OFINTHIGM/CAI= 0. SHODIHIGH/CAD= Oo SMVDTHIGH/CAD= Oe
DGCOTAICAL/CAS= §.955E-01 JOFOTHICAL/CAI= 1.826-04 SYHOTHIGM/CAD= 0. SOBDTHIGH/C4)= 0. SLODIH(GH/CAD= 06
SUMILICALDS-A.613ECIC SUMFACUMD= 2, SFNOTHIGH/CAI= 0. WKNET(CALD==1,616£402 RRICHS-ATH/GH-KI= 2.800F 400

UCHLEICH/SECHS 2429AE903 OVOTHICMS/CA) =91.G11E001 OFLEIZ= 2.407E-03 TTERPC
YHNIFIGMI= L.572E-35 VNOE= 6. 3A2F=$1 XNOE(POMD= G.137E-25 SONTLIGMD= 2.651E-17

2

THETOCOES CAD= -25,.90 PUATHD= 2.012F6G2 VICMID= 2,616F 002 TEK)= B.O98EFI2 UICAL/GM)= G.2G6EOL PHERE 1, 408-02

DOTHEATH/ChIS B6E33E~GL WOLCAL/CO)=-6.949E 40 OFINTHIGH/CAD= 0, SHODTH(GH/CAI= Of SMVOTIIGH/CAD= O.
OICUTHICALACAD= 7.939E-01 OIWNTHICAL/TAD= 3.337E-04 SVAOTHIGH/Cad= 0. SOBDTHIGH/CAI= 0, SLODIN(GMH/CAIS Oe
SUMCLICALE==1,4656600 SUMFUIGHI= 9. SFODTHIGM/CAT= 6. WKNET(CALI=-2,090€002 RRICHS-ATH/GH-KI= 2,86CE +00

ULHAPZICH/SECD= 262198403 OVOTHICN3/68) =-1,95AE001 OLLEL G99E-03 TTERPCS 1
XMNCEIGM) = 2.455E-26 YNOE= 9.9H7E-27 XNOE(PPM)= 9,585£-21 SQOTEIGM)= 5.182E-16

TRETA(DEGCAD= -19.00 — PUATHD= 2.6765¢01 VICHSD= 2.251F 402 TIK)= B.G11E402 UICAL/GH)= 7.206401 PHEE= 4. 889E-c2
COITHEATA/UEDS 162078460 WICAL/CAI=-5.325E400 OFIDTHIGH/CAI= 5.8622-03 SHNOTHIGM/CAD= 9. SMVOTHIGH/CL)= Cy
DTSOTAICLESCAY= 14215E 00d DOFDINICAL/EAD= 5.125E-04 SVADTHIGM/CAI= Oe SMAOTHISH/CAI= O. SLOCTHIGH/CAD= de
SUACL(CALD= &.691E470 SUMFRIGM= 0. SFROTH(GH/CAD= 0. WKNET(CALD =-2.405E002 RRISHS-ATH/GH-KI= 2.846E006

UCHLOICH/SECI= 24156E493 OVOTHICHI/CAD=-8,211€400 OELEIZ= 3.674E-03 ITERPCS 1
KMNOTIGIDS S6H21E-24 NOES 1. 39t-24 XNOE(PPMD= 1.356£-18 SOOTEIGMD= 2.622E-15


GO TC 16
C CPITICAL CNIGN SKIN ASSUMFTICN
45 XL= CPLY(TL~-TLOD
1h CPLXL= CPRASSXL
Re= CPLYL*(TS-TL+ QCP¥YC2Z84)
ALCG2= ALOGI1,3+82)
YFL= 160-(1.0+¥9285)701.0482)
BLOGG= ALCG(1,0+YFL/(1.9-YFL))
BAA= ALOG2/4LOG4
C FLOME TEMPERATURE
JFe (TS4TL FYC23S/SYFL+CCP*YO28U V7TL.O4YG2BG/YFL)
TETAS TE/T8
XLHS= PUAX FEZ
OR B= -0.05*(SMAX-8L)
R3= OMAX-C2B
BO 4 ITS9=1,50
RA= Oh4 ORF
IF (RP.LE.RL) GO TO 12
PNEMP= Da/(RR-2L)
YMQOTS SLFRORMO*AL CG?
¥XX= RORMOEADA
RPFQLS XXXAUXXX1.0)
PFs RLFSFOL
XXXK= XMOOT/EF
TF (X¥X%.GT.2900) XXXX=200.
ERF= EXOTXXxX¥)
XXXX= XMNOT/IB
IF (X¥XY¥.GT.209.) XXXxX= 200,
ESI= EXP ($X¥XX)
OPTS (PMAX-=FFI/30.
fT= RF- Nel
SuM=0.0
NO 2 ITR=1,30
RT= <+DOT
Eee EXP C-XMDOTZ°L)
KTTGO= (EQF-ER)*RIFRIZ (CC ERF-EQ) FTE TB *(ER@ERB))
IF (XITGH.LT.O.M) XTTGO=9.0
St cuv+ XITGO
2 CONTINUF
XFHS= 3.0*SUMFDOT
CHK= (XLHS<XRHSI/XLHS
IF (CITEF.LEQ.1) GO TO 4
TF (ARSICKK) LT.C.01) GO TO F
CHKX= CHKLFCHK
TF (GRKX.LT.C.9) GO TO 5
& CONTINUE
CKK1= CHK
GO TO 1
5 CONTINUE
QR= Re ORR
NFB= £.5*D28
CONTINUE
ITE(6,161) PsTA,TLs»TLO,YC28,RLGFMAY
WETTE (6,100) ITER,XLHS,XRHS,CHK,ARF,TFE,RB PT
sTae
© IN TSeANSTENT M-S OIFFUSTCN FLAME CAN NOT GE ESTABLISHED
42 CONTINUE
XL= YL2O8* ((TCRITF<TLCIS (TCOITE-298.0)**0.38
HXKS 160+ CEMARFLTB-TLOIV/XL
MPMLITIA £2.56F V7OH2]FRHCDFOL*ALOG (XXX)

71


14 CONTINUS
XMDECTI= ¥YMOF (CT) + XNGICI9L) #XXXMCF *XMCFL
41 CONTINUE
OMLIT(T,Ld= CMLOTIT
XKK(T,L)= (3.6366197723/F0EN/RL) FOMLOTI
OLIOT= CLENT+ IMLOTLI*XNAT(IVL)
10 CONTINUE
C OVERALL DIFFUSIVE ZONAL BURNING PATE
DLENTE(TI= ALENTFOTITH
TF (COLBOTHETIFOTHETA .GT.XMFLI(I)) DLBCTHITI= XMFLICI)/DTHETA
20 SCNTINUE
100 FCR4ATCIX,I5,1P12E19. 7)
LOL FORMATOS/1X,% ee0e FOPFOR IN FLAME ITERATION «oe * eI P7E10.3
RETURN
END

73


#NECK HFTNCR
SUBROUTINE HETNOR(P,YNC8 ,AMW,RHCDsNITST,SCNCR,CNOOTID
© HETERGENOUS NC FCRSMATTION FATE
RIMENSION Yi509 52050) ,WOOTL CSO) »XKZ(°50),XK 3050) pSAVEL(S19,
1 SAVE2(51),WCOTAV(50)
COMMON/FLAMESY AMWMP,RMAX,YN2L,TF,RF»SXMOCT SRL RS ,TOF (SO) »YFCF (50)
1 YO2OF(50),YN20F (5G) ,YNCOF (50) »WOOT(50)4R80F (50)
NITST1I= NITST- 1
DO 7 T=1,NITST
wOCT1(T)=0.0
XK2(T)= 0.0
XK3(T)= 0.0
IF (YC20F(I).LE.0.0) GO Tel 7
C SATE CCNSTANTS
RT=1.9877 TOF (1)
¥KFL= L-26E14*EXP(-75400./RT)
S.12ELZFEXP(-400./8T)
1.33E10* TOF (1) FEXP(-7080./PT)
XK 3.2E9F TNE CT) FEXP(-39190./RT)
XKEQS= 5.0FEXP(-58300./RT)
XKi= CNOP*S.0S5E-4*(XKF1*XKECZ)
WOOTL(I)= XK1* (CAMWMP/TOF (IV) *¥1.5) FYN20F (1) *SG8TCYO2Z0F(T))
XK2Z0T)= 1.9955* (XKBL*xXKI2/ (XKFL*XKF2))/ (YN20F (1) *YO2DF (19)
XKI(TI= 1.14297 (XKE1/XKF 2) /YC2DF (1)
7 CONTINUE
© INTEGFATING NO CCNSERVATION ECUATION
OY1L= 2.0E-4
IF (YNO8.GT22.06-3) DY1= G.1*YNOB
Yi(1)= 2.0€-4- D1
WOCT(1)= 0.9

CHKUPP= 1.0619
CHKLOW=-1.0€10

2 CONTINUE
Yitto= Yitiy+ OVt

& CONTINUE

© IMITIAL CONDITIONS
IF (ITER.GE.NITST) GO TO 4
Y2C1)= YMDOTFY1¢1)
YNONF(1)=Y101)
eo 1 Ts1,NITstt
~ RI= O.S* (POF (1+41)4+ROF (II)
CPT= RDF (T+13-R20F (1)
RI2= RI*RT
YACT#1d= YLCTI+ C¥2¢T)/R12)F0RT
TF CY1 CT#4) LT. -1.06-29) YiCI+i)= -1.0E-20
TF CVYACT#1).6T.100) YAICI+19= 1.0
Yeoor(L#1)= yacr+d)
C PSODUCTION SATES

CORVL= 1.60- XK2CT #1) *(YNCOF(CIT4#1) #YNOCF (T4199
FORe 1.0% XK ICT tL) FYNONF (T4419
CORR= CNPRA/CORR]
WOOT(T+#1)= WOCT1(I+1) #CORP
WOCTAYCTIY= CoS*CNDOT (TY #WOCT(T+1))
Y2CT#1)= Y2CT) + COXMDOTFYZCTI PT 2)-CRIZFWOGTAV(II/RHOC) I*DRT

1 GONTINUE
ITER= ITER+ 1
SAVELCITEPY= YLCL)

74


TCRITF ©

CPFL
CPFYV
CVF
OHFV
HEFL
HCON
SEND

O.S67E4C3,
O.S51E+00,
0.332E+004
0.32E+00,
0.92865E+02,
70.26439E+03y

Oe,

29


DELPHE= PHETS(IOPHE)- PHETBC(IOPH1)
DPi= (P-PTBCICP1))/0ELP

OTi= (T-TTACIOTidy/OELT

OPHE1= (PHE-PHETS (IOPH1) )/CELPHE
OP2= (PTBCTOP) =P) /NELP

NT2= (TTBCIOT)-1T)/DELT

DPHE2= (PHETA(IOPHE)- PHE) /NELPHE
CALL TRIPIT(XH2TB, XH2)

CALL TRIPIT(XH20TS8,XH20)

CALL TRIPIT(XCOTS,XCO)

CALL TRIPIT(xC92TBA,XCO2)

CALL TRIPIT(XN2TBs XN2)

CALL TRIPIT(XNOTB,XNO)

GALL TRIPIT(XGHTB,XOH)

CALL T@IPIT(XHTB, XH)

CALL TRIPTT(XCTBy XO)

CALL TRIPIT{UXNTB XN)

CALL TRIOITCXC2TB,X02)

CALL TRIPITUXCHGTByXCH4)

AMMC=26 *XH2#28 6 FXN2418 o FXHZ0417 0 FXOHF1 oe FXHETL a XNCHLGS FXNt3 24 * X02
1416. 7XO4G46FxC024+28. FXCO+16. FXCHG
SETURN

EMO

84


SUBROUTINE EVAP

The model requires that the premixed air-fuel packages are tracked
throughout their lifetime and this operation is handled in the subroutine
EVAP. The vaporization rate of each droplet size groups and the air-fuel
vapor ratio of each package are defined in Lagrangian coordinates. The
onset of ignition in each air-fuel package is also monitored.

SUBROUTINE FINJIN

The rate of fuel injection is calculated by calling this subroutine.

SUBROUTINE FLAMEZ

This subroutine computes quasi-steady flame structure by constructing
the diffusion field around a confined burning droplet with finite space. Flame
temperatures, flame zone thickness, fuel, Oo, and temperature distributions
are determined. Subroutine HETNOR is called in FLAMEZ to calculate the
production rate for nitric oxide during droplet burning. The nitric oxide
produced during droplet burning forms the heterogeneous contribution to the
total emission.

UBROUTINE HETNOR

The nitric oxide formation rate within the flame structure defined by
FLAME is determined in this subroutine by integrating the finite-rate NO
conservation equation. This two-point boundary value problem is solved
by a trial and error method. An initial guess of the NO concentration at the
droplet surface is made. The NO production equations which are based on
the Zeldovich mechanism are integrated and the calculated NO concentration
at the end points is checked with the given value at the end point.

SUBROUTINE HOMNOZ (HOMNOR)

The homogeneous contribution to the total amount of NO produced
(which was not included in HETNOR) at the zonal NO formation rate is
determined by this subroutine. The Zeldovich chemical kinetic mechanism
is also used as the basis for the calculation of the homogeneous NO
formation. Subroutine HOMNCR is used for that of the "air" zone.


II. PROGRAM DESCRIPTIONS

The main body of the program is the routine MAIN which simulates
the cycle thermodynamics of a reciprocating engine and also manipulates
the major input output mechanisms of the program. The cycle thermodynamic
routine is used as the skeleton and the model components, coded as
separated subroutines, are linked together by CALL statements, thus forming
the complete program. The flow chart is shown in Figure l.

A brief description of the subroutines included in the program are
given below.

SUBROUTINE AMIX

This subroutine estimates the rate of mixing between the “air” zone
and the burning zones.

BLOCK DATA

This subroutine presents the coefficients for both the internal energy
polynomial, U(P, T, 9), and the gas constant polynomial, R (P,T,¢), for
lean (9 s 1.0) and rich (> 1.0) equilibrium combustion products. The
coefficients are obtained for a specific fuel by curve-fitting the calculated
results from the NASA equilibrium program (Ref. 2).

SUBROUTINE CONHT

The gross rate of convective heat transfer from the entire charge
within the cylinder to the cylinder walls is computed in this subroutine.
The calculation is based upon the mean charge temperature and the character-
istic charge Reynolds No.

FUNCTION ERF

This subroutine estimates the error function of a given argument x, i.e.,

2 x 2
erf (x) =

-C
d
nf J


SYMBCL DESCRIPTION

AMWA Molecular weight of charge content
FARA Fuel-air mass ratio of charge content
PHEA Equivalence ratio of charge content
FAR@S Stoichiometric fuel-air ratio

XMEGR Mass of EGR (gm)

XMAEGR Mass of pure air contained in XMEGR (gm)
XMFEGR Mass of fuel contained in XMEGR, (gm)
XMRD Trapped or residual gas (gm)

XMARD Mass of pure air contained in XMRD (gm)
XMFRD Mass of fuel contained in. XMRD (gm)
AMWF Molecular wt. of fuel

@XST Stoichiometric oxygen-fuel mass ratio
X@2A Oxygen mole fraction in charge

FDEN Liquid fuel density (gm/om?)

HHV High heating value of fuel (cal/gm)
HCOM Heat of combustion (cal/gm)

HFFL Heat of formation of fuel (cal/gm)

HFJ Enthalpy of injected fuel (cal/gm)

VC Engine clearance volume (cm?)

VDIS Engine displacement (cm)

VFI Initial velocity of fuel injection (cm/sec)

(3) Printout of Current Information at Each Specified Print Station

The program will printout averaged over-all charge properties and rate
information at the specified print station controlled by NAMELIST CPARAM and
every ignition step (from the onset to the end of zonal ignition). The elements

of the printout are described below.

18


SUBROUTINE TPROPZ (TPRCP)

This is a subroutine to evaluate the thermodynamic properties of
the equilibrium combustion products. The internal energy and gas constant
(R) are calculated as functions of P, T, and» by polynomial. Expressions
and coefficients of these polynomials are specified in DATA BLOCK. The
partial derivatives (of U and R) with respect to P, T, and %, are also
computed in subroutine TPROPZ. Subroutine TPROP is used to evaluate the

thermodynamic properties of the "air" zone.

SUBROUTINE TRIPIT

Subroutine TRIPIT is a three-dimensional linear interpolation routine
called for in subroutines SPECY and SPECYZ.

SUBROUTINE VOLM

This subroutine gives the instantaneous volume of engine cylinder.

SUBROUTINE WETBT

This subroutine estimates wet-bulb temperature of a vaporizing fuel
drop at high ambient pressures and temperatures. The Cox-Antonie equation
is used to evaluate the equilibrium partial vapor pressure.


Ve Gt .5*PFTL/PL
CVECVXFFL
Cr PXTFS
GhvM= CP/CY
C PARTIAL DEPIVATIVES OF U AND R W.ReTe Poly 2NO FRE
O1PF 00112)+3.0*001(3) #PHKF2
Derr (901°) +3.0*001 (6) FPHE 2) FTI
OzPS =C (ONTO +3.0*N0119) HOHE 2) + NOL CLL) *TIDFSLOGP
CPF= C*(DiPF+C2PF+02PF)
PFLPF= -FL FF LFF AROS
PUPS FLFC (930791)
PUPT = CV
FUPF= FL (-8+CPF) + PFAPFFUX
PEPD= FLFCE(RPL(S) FPHEFRP1 (99) /P1
PROTS -FLFCF(PRL (4) FPHE+ROL(7)9F1000./T2
FOXPF= RO4(2)+ C¥CRRL(3)+ RALC4)¥FTI+ FR1I(5) FALCGP)
PEPF= FLIFPRXPFt PFAPFERX
RETURN
2 CONTINUE
© FUEL OICH
A= MAVSCLIFAAZS(2IFTAFANSCIVFTZ FAA2(GIFTI #OA22(5)FTH +AA2 (6) FTS
P= APACLI+ATZ(2VFTA+BESCIIFT2 FAB2(GIFTZ +ERZ(Sd*¥TH 4PB21E) FTS
Mls PN261)94+902(2)FPHF + COZZI FPHEZ
Des COC2(4) +002 (5) *PHE+ COZ(E) FOHES) FTI
NICHC(NH2(77+4+OC0208FFPHE+ OO209) #PHEZ) + (CO2C10) +4NO2]¢11) *PHE) FTI)
Oss OF0*4L0GP
XKX=914024N72
TF (XXX .LT.-200.) XXK=-260.
N= EXPCYXY)
UX= A- PEPHE+ OD
XXX= (PO2SCZP+ QOZTG)FTI+ CRIS) FALCGP)FOKE+ RI2(KhD + FR2C7) FTI +
LRE2(BY*ALOGP
IF (XXX eLTe- 2004) XXX=-20%%
FXF (XXX) .
5 RR2(1)4 KR2C2)¥PHE +
HYSUX4+RX*¥TY
AP= ARQC2V #2 OFAN2[](ZVFTL+ SB.0FAA2(4)FT2+ 4oOFDAS(S)PFT3+ 5.0F A026)

1 FT4
Pps BP2(2)+2.0FEB2ACTIFTL+ F.CFRR2CHDET2+ 4OFERZISV*TZ +5.0*882(6)
1 *T4

~{(O0N2(4)+ NDSCS) FPHE+ 902 (6) *PHES) 71000472
~(0N2(10)+ O002(11)*PHE) F1000./T2*% ALCGP
0 C¥IO2PeD IP)
CVX= OP-BP*CHF+eNP
CPX=EWX +X
CONVERT UNIT 10 PFS MASS OF MIXTURE
XRASE= 1. 0+F ASOS*PHE
Fi= 1.0/XRDSE
U= UY*F1

Gave CP/OY
C FBPTIAL NFRIVATIVES CF U AND R Wa®eT. PyTy ANC PRE
C1PF = 002(2)+3.0*0N2(2) FPHEZ
Ceer (90205) +3.0*902(6) =PHE?) FTT
NSPE =0(9652(8)+3.9*9ON2(9) *PHE?) 4+ CGSC14)* TT) #SLOGF

93


GO Ta 6

@ CONTINUE
PHEL= PHETS(4)
IDPHE= 2

& CONTINUE
Iopi= TOP-1
IOti= IOT-1

TOPH INPHE-1

OFL PTACTIOF)= PTACIGP1)

DELT= TTB(IOT)- TTBCTOTI)

QELPHE= PHETBC(IDPHE) = PHETBCIDPH1)

OP1= (P1-PT8(19P1))/ DELP

NT1= (V1i-TTSCIOTAII/DELT

DPHEL= (PHE1-PHETBCTDPH1)) /DELPHE

OP2= (PTBCIDP)-P1)/DELP

OT2= (TTBCIDT) -T1)/DELT

DPHE2= (PHETEC(IDPHE)— PHEL)/DELPHE

COLL TRIPIT(XK2TB,XH2(T))

CALL TRIPIT(XH20TB,XH20(T))

CALL TRIPIT(XCOTB,XCO(TI))

CALL TRIPITIXCO2THyXCO2{TI)

CALL TRIPIT(XN2TS,XN2(T)9

CALL TRIPITCXNOTA,XNOCID)

CALL TRIPITCXCHTByXOH(IDD

CALL TRIPTTCXHTS,XHCT)D

CALL TRIPIT(XOTBsXOCT))

CALL TRIPIT(XNTSeXN(I9)

CALL TRIPIT(xC2TB,XO2(T))

GALL TRIPIT(XCHGTR yXCHG(T))

AMWBMCTI= 26¥XH2(T)+ 28. ¥*XN2(IT)+ 18.*XH2C(ID+ 17.*XOH(I)+ 1.*XH(T)
1 + FOoFXNOCTP + 14.FXNCT)# 32. 7XO2(TI+ 16.7 KOT 446F
2 XCC2CT) + 28.*XCC(I)& 166*XCH4CT)

19 CONTINUF
RETURN
END


FOS 1.0-FTA22
f= (FCO*FR) FETA
ETA= ET&- G/NGDE
IFCETASLT.O.C) ETA=0.00C1
IF CETALGT.1L.66) ETA=0.999°
S CONTINUES
WOITE (6,103) JIsTsGeIQsETA
sToe
Cc AVF, VELOCITY OISTRIBUTICNS--- FREE JET REGION
& CONTINUE
ETOI(IyJ+id= ETA
AB FCH L.0-(OSFCETALT(T UHL dV FETAL CIT 9 Jd0) F105
F(l,J)= RHOFOLFFC
ETAT(T.J) #9
ETAT(TI,J+1)*8
PL¥(R2*R2-R1*R1)
VEAP OUF/RHCF(T,J) /OAREA
VE(Ls dbs VRAPALVF IO
AFT(TyJ)= (1.0/VS(T,J))-1.0
CET(L,J= BHCR/OHOF ITS)
6c To 9
7 CONFTRUS
G WALL JET REGION
TF (JsEC.NSHELD GA TO 17
GM= FLOAT( I) FOUFY (2.0 FFT FC RX4EW) FERVCL FREGFCLD
co 10 ITEC=1,29
FTOL= ETA
ETA2=(7,9/2.5) # (ETA**#2.5)
ETAZ= (7.9/4.0) FCETAF FG)
ETAM= (1.0/5.5) #(ETA*F5.5)
G= CTAL-ETA2tETAZ-ETAG- CC
CKK= ARS(G/GC)
TF (CHK ALT.1.068-5 9 GO TO 11
ETA3 ETAF*1,5
FO= 1.0-ETA32
CGDE= FCR*S
ETA= ETA- G/CGNE
TFCETAWLT.O.0) ETA=C.0001
TFCET&L.GT.1.69) ETA eae
10 CONTINUE
WRITE (AG193) Jel yGeQ0,ETS
sT9e
14. CONTTNUF
© B/E, VELCCITY NISTRIPUTTCNS--- WALL JET REGICN
ETAT(I,J+1)= ETA
17 FO= L.O-(OSFCETAL CTs Uti) +E TAL (Ly Jdd) F105
FHOF (I,J)= RHOFCLYFC
Rg ETAT(I,J)*3
R2= ETATCT, +1) *8
MARCOS 2.0FPTF(Ox¥+BW) # (22-84)
VGAR2= NUF/PHOF (I,J) DAREA
VS(T,d)= VBAR2/VF IC
AFT(T,J)= RHCE/RHOF IT, UJ)
9 CONTTNUF
YS(LsJ)= 0.57 (81462)
CYS(I,J)= 82-81
IF (J.L4.1) €C To 2
G THTAGPATY FOR PENETRATION TIME
VOAR= C.SFCVS (I-12, 5) #VSCT, Sd) FVE UC
PTIM(T,J)= PTIM(I-t,4)+ (OX/VEAR)

a

oao

89


TCRITF 0.567E4035

CPFL = 0651E+06e
CPFYV = 0.332E+009
CVF = Ue32E+005
OHFV = 0.92365E+02,

HFFL = -0.2E439E +035
HCOM = 0604
SEND

24


__|EPA-460/3-76-008-b
| March 1976
|
| STUDY
| ON OXIDES OF NITROGEN
| | AND CARBON FORMATION
IN DIESEL ENGINES -
COMPUTER PROGRAM
: USER’S MANUAL
LT
fe \
—( rey,

U.S. ENVIRONMENTAL PROTECTION AGENCY
Office of Air and Waste Management
Office of Mobile Source Air Pollution Control
Emission Control Technology Division
Ann Arbor, Michigan 48105


2 CONTINUE
XS(Ih. = X
BBB(I)= B
21-CONTINUE
‘EFFECT OF SWIRL CN ENTRINNMENT
USBAR= 0s25°CBOWL*CMEGA
ER=. 100+ SEFJ*(USBAR/VFUC)
00 21 J=1,NSHEL
O00 21 I=1,NAXIL
AFT(TyJ)= ERPAFT(I9J)
21 CONTINUE
© DETAIL OUTPUT
IF (IDTAIL.EC.0) GO TO 22
WRITE (6,116) (I,I=1,NSZG)

WRITE (65117) (SZBARCI) ,I=1,NSZG)
WRITE (65115) € GDSZ(I),I=1,NSZ6)
WRITE (6598)

WRITE (6,107). (XSI 3 y1T=1,20)
WRITE (6,106) € BBBC(I) ,I=1,20)
BO 20 J=1,NSHEL
WRITE (6,98)
WRITE (65108) € YSCI,J),T=1,20)
WRITE (6,114) ( OYS(I,J),I=1,20)
WRITE (6,1099 (RHOF(T,J),1=1,20)
WRITE (65119) (AFT(I,J) 5721220)
WRITE (6,110) (VS (I oJ) 9151528)
WRITE (6,111) (PTIM(L,J) yI=1520)
20 CONTINUE
22 CONTINUE
98 FORMAT(/)
103 FORMAT(/1X,* JET SPRAY LCCUS ITERATICN GOES NCT CONVERGE *y2I55
1 1P5€10.3)
106 FORMAT(1X,* @ *,20F6.2)
107 FORMAT(1X,* X *,20F6.3)
108 FORMAT(1X,*% Y *,20F6.4)
109 FORMAT(1X,*ROF*,20F6 24)
116 FORMAT(1X,* V *,20F6.4)
111 FORMATC(1X,*TIM*,20FE 04)
112 FORMAT(1X,*BEFORE TRANSITION*,2X41P9E12.3)
113 FORMAT(1X,*AFTER TRANSITIOCN*,2X%51P9E12.3)
1144 FORMAT(1X,*DY *,20F6.4) .
115 FORMAT(1X,*GROUP MASS FRACTICN*®,5€12.3)
116 FCRMAT(/1X,*CROPLET SIZE GROUP *95(7X,I59)
117 FORMAT(1X,*MEAN GROUP SIZE *,5E12.3)
118 FORMAT (/27X,1HB,9X »GHRHOFCL 4 6X5 SHVCL 4 SX y SHXSM)
'119 FORMATC1Xs9AFT* + 20F6.~ 2)
RETURN: .
“END

90


G INTEGRATION sTeP SIZE CONTROL
TF (THETAI.GE.THELJS) GTHETA= OQTHMIN
IF (CTHETAT.SLTeTHETJ) LAND. CCTHETAL+OTHETA).ST.THEIU)) OTHETA=
1 THEIJ-THETAT
IF (SUMFS.GT.XMFI95 ) OTHETA= OTHMED
IF ((THETAL#+OTHETA).GTeTHEVO) DTHETA= THEVO-THETAT
OTHET2= G.S5*OTHETA
DT= OTHETA*DTOTH
OT2= G.5*0T
C CHARACTERISTIC VELOCITIES
UP =OSOTH*OTHOT
XXX= USIFUST+UR FUP
YT= SQRTCXXX)
XXXX= (XM AO*XXX +SUMF JFVF JO FVE JO) ZC XMAD# SUMF J)
UCH4R= SORT UXXXX)
CG CHARATERISTIC RENOLOS NO.
XMU = ‘4.25E-5*(T **1.5)0/(T #110.)
REPBAR= RHOFUPBARFIE/XMU
REI= RHO*UT#DE/XMU
REJ= PHOFUCHAR*DE/XMU
C*** TOTAL CONVECTIVE Hols
CALL CONHT(REPBAR ,T,0QCOTH,CA,CB)

SLOOT
aosoT
XMSOTU= 0.0
IF (THETAI.LT. THETU) GO TO 86
IF (I@UG.NE.9) IPRT=0
IF (CTHETAT-OTHETA) .GE.THEIY) GO TO 3
S ESTIMATE SPRAY TSAJECTORIES
OELP= PINJ= 9
CALL SPRAY (RHODA ,IMEGA, DELP 4SEF Jy COCHG)
C INITIAL FUEL PROPERTIES
CALL WETBTIP,TO,AMWA ,AMWF,4,8,CeTFLyDHFVP)
UFL= UFFL+ CPFL*(TFL-298.)
TFG= TFL
HFL= UFL+ P/FOEN*FATOXS
HFV= HFL+ DHFYP
UFG= HFV= RFGTFG

SMAPB= 6.0
svaot o.9
DFGOTH= 0.9
DFVITH= 9.0
OFJOTH= 6.9

0.0

a.0

TFLB= TFL
UFL UFL
TFGO= TFG

IF CICRITALEN.0) GO TO 17
HFOS=. UFCOIT
TFOS= TCRITF
60 10.3
17 CONTINUF
HEOS= HFL
TFOS= TFL
3 CONTINUE
S SATE OF FUEL INECCTION
IF (THETAT.LT.THEET)
{CALL FINJN(THETAL,OFJOTH,XMFJ)
IIOLo= [1
© ESTIMATE THE VAPORIZATION OF FUEL SPRAY
IF (KOUNT.GT. NSS ) GO TO 25

50


When it is unnecessary to modify the parameters given
in Tables 4, 5 and 6, a $CPARAMS card starting at
Column 2 must be included. It should also be noted
that when NSHEL and DTHMIN are specified, the
constraint:

(DINJ/DTHMIN) x NSHEL = 100
must be satisfied.

(d) Namelist EPARAM

The empirical parameters, which are preset, required by

the program are provided with this set of input cards.

These parameters which are preset prior to each individual
case need only to be input if it is desired to alter the preset
values. Therefore, when the alterations are required, they
must be made in every case.

If modification of the parameters listed in Table 7 is not
required, a SEPARAMS card starting at Column 2 must be
included. .

(e) Namelist GIVEN

The variables in this set of input are the same as that of (a).
This set of inputs serve as a modification or as an addition
to set (a). When neither a modification nor an addition
complement is desired a $GIVENS card must be provided.

The order of card input is illustrated in Figure 2.

10


Ill. DESCRIPTION OF THE PROGRAM INPUT
\

The data input to the diesel engine program consists of three distinct
groups:

(1) Tape Input (Species Tables): The mole fraction distributions,
X; (P,T,»), of the equilibrium combustion products is read in
as TAPE 8 from either a magnetic tape or magnetic disc. This
input is required only once for all stacked cases having the
same kind of fuel. Cases using different fuels cannot be
stacked.

(2) Block Data Input (Thermodynamic Properties): The internal
energy, U (P,T,%), and gas constant, R(P, T, #), of the
equilibrium combustion products are calculated from fitted
polynomials (detail, see Ref. 1). The coefficients of these
polynomials form an integral part of the program via the
BLOCK DATA subroutine. The BLOCK DATA has to be revised
whenever a different fuel is used.

(3) Card Input (Engine Geometry, Operating Conditions and Air-
Fuel Physical Properties): Five sets of input cards are required
in this group. All of them except the case identification card
are input via the standard namelist format.

The five sets of cards are described below according to their
input order. The definitions, format, and preset values of
the variables in each set of input cards are described in
functional subgroups. The order of the variables within each
set of input data is arbitrary.

(a) Namelist GIVEN

This set of input cards details the common constants of engine
geometry, engine operation, and fuel physical properties. The
variables in the GIVEN namelist are presented in Tables 1,

2 and 3.

(b) Case Identification Card

A single card contains the case identification information in
Columns 1 through 80. When no identification is desired a
blank card must be provided..

(c) Namelist CPARAM

This set of card inputs provides the logic flow control parameters
for the program. All the parameters in this set of data are

preset to nominal values and only require to be input if it is
desired to alter the preset values. These preset values apply to
each individual case, therefore, all charges must be made -
whenever required.


FOFCK ATX
SUBPOUTINE AMTX(XMAD,SUMEM NCHA SATAEM AFP CMR OTH sy DMVCTHSCMIX,
1 FyMTcr)
© THIS FOUTINE ESTIMATES THE RATF CF MIXING BETWEEN ATR ANG SUPNING ZONES
MIMENSICN OMGSCTH(100) »BTLAM(190) ,OFR(100),OMVCTHI190)
COMMON/SINTEGR/S INOLsINO2,TSOCT,IWETBT sNPRT,ICTAIL yNTERyNAXILy

1 NSHELNSZ2G,NTTST,I°UG,II,ITIT

CCOMMCEZENGECS DE, REWyRL VO 4SC oS yUST yCRy RPMs TWyONOZZ yNNOZZ,0RCHL
1 vFJc.

CES O.5¥*0E

AP= 3.,1U150%RE FRE
V= VC+ APFS
XLBAR= VZAP
YLCHAS= (FE*XL3ARI/2.9/(REtXL 9A)
COMTX= OMIX#UCHARSXLCHARSEPM
IF (XMAM.LT.O.0) XMAA=0.0
XXX=XMAA*SUMOMS (SUMBM+XMAAD
DCMADTHECOMIX*XXX
SUMBTA= 9.9
90 1 LT=1,IT
SUMBTS= SUMETA+ BIARM(T)
41 CONTINUE
NO 2 l=i,1r
DMOITHCL)= C8TABM (CT) /SUMETA) FOMADTH
OVVOTECTY= CVIOLTFOMOOTE(T)
? CONTINUE
BETURI
ENO

63


ro)

oO

oO a ao

oOo aA a a

oO

Qo

NSHEL= 4&
NAXIL= 2
NITST= 5
TeuG= ¢@
IWErBT=
ICRITA=
NSTROK=
INTEGRATION
DTHMINS
OTHMED=
OTHMAX=
THEMO&X= THEVO
RELATIVE ERROR ROUNDS
EPSIT= 2.0E-3
EPSIP= 2.6E-3
EPSINO= 1.9
EPSIST= 1.0
PRE-SET EMPIRICAL SONSTANTS
ATR/PRODUCT MIXING
CMIX=7.8F-3
CYIJI= set
HETEROGENEOUS BRUNING RATE
CBURN= 2.9
XTLYOZ= 6.02
VAPORIZATION
cEVAP= §.276
HUTIM= 2.0E-4
IGNITION DELAY
ATGNO= 1.460E-4

2
a

P SIZE (CIN DEGREE CA.)

VNR UFOG

EMPIRICAL PARAMETERS ACCOUNT FOR SWIRL EFFECT

SEFJ= 4.3
NOZZLE OISCHARGE
COCHG= 1.6
IGNITION DELAY
OTPAHY= 2.0E-4
HEAT CONVECTION
Ch= 9.26
c5= 0.75
DROPLET
RHOD= 3.0E-4
CPBAR= 2.35
NO KINETICS
GNOR= 5.37
SOOT KINETICS
SCEN=2.0
OSOOT=1.0E-4
AsoOTa= 3.52
ASOTP= 3.52

COEFFICIENTS FOR VAPOR PRESSURE (COX-ANTONIE
A= 4.091
R= 1324.77
C= 53.74

OHEFG= 3.9
WRITE (6,167)

INPUT CASE IDENTIFICATION
PEAD (5,99) (109 (7T),1=1,8)
IF (ECE(S)) 999,34

34 CONTIN
WETTE (6,99) (IDD(1),1=1,8)

ON

EQe)

SPRAY GROWTH

MCOTFICATION OF FONTROL PARAMETERS THROUGH NAMELIST INPUT

46


FPECK SPRAY
SUIROLTINE SPRAY {RHOR,CMEGA,CELPs SEF JyCOCHG)
2 TRIS FCUTTNE STMULATES THE STEADY STATE FUEL JET SPRAY ANC DEFINES THE
C TROJECTGFIES CS FUEL ELEMENTS AND A/F VIA SIMILARTTY ASSUMFTTON
DIMENSION BAA(20)
CCMMONZINTEGR/S INO1,INC2,ISCCT,IWETET sNPRT,ICTAILsNTER NAXIL,

1 NSHEL,NSZG,NITST,I@UG,ITI,III

COMMCHR/SCA/S THEMAX,THEIJ,THEEIG, DINJ» THENEsTEEVCyTHEVCOTING,
1 CTHUIN,OTHMEN »ITHMAX »CTHETA,OT

CCHUMON/SENGEOS NE, RFWyRLyVCoSCySyUSTsCReRPMy TK eONCZZyNNOZZ,OROWL
1 VFJO

COMMONS FUELF/W eo XX yYFLsZ2EL YSTOA, ZSTOR, APT »FCEN »FAROS y AMWF SUPFT
1 aXMUF yCENOXST»HOOM

COMMONS SETOZS OFF yXS0 20) 4 ¥S(2994) pVS( 2094) yRECE (2094) pAFT (2094) 5
PTIM( 2054) ,ETAI (20,5) sOYS(2054) 5 GOSZ(3),SZRAR(3),
0025544 3) 9 XSS(25 44) pYSS (2554) gLYSS(25 54) sTAFVI25 94) y
AFV (2554) sXMFVS (2594) yXMSOT( 2554) ¢TO(2594) POV 25 944
INIG(2544) a
COMMON/SCONVF/S XJpATOSNTRy XK y SIG, 9 ,RUSCTOTH,CTEOT,DIROT PI
RAYELTST/GIVENS FOEN,SURFT,XMUF DELP sRHOF yONCZZ »OBONL VF JOeCIF
1 sCWF yCUF XI ,BTyNSHEL»NAXILGNS7ZG sCUFoOX,XWALL BC, YMAX
XIFL= 1.3446
0.9195
16559272539
0.4209
6.33115
0.2709
7= 0.2275
XT2= 2.9774
YI3= 2.3122
XTG= 1.9854
XIS= 1.74872
XTe= 1.57383
XI7= 1.437
CIF=0.135* (1.9+RHOP/FOEN)
CJF=0.175
eMrsg.206
REONL= 9.5 *0R0WL
RNCZ «SFDNCZZ
ANOZ PT*RNCZZ*ONCZZ
XJF= FOEN*VFJUC*ANOZZ*COCHG
XJM= XUFFV FIO
QuF= XJF/NSHEL
C GEFINING THE ATMFNSION CF INITTAL REGION
PIF= PI/XTE3
XT= (SORT (PIF) 1,0) #RNCZZ/CTF
ATs CIFFXT#RNCZZ
c XWALL= (QNOWL#CUFFXI-8107(1.0+CUF)
XWALL= PROWL
BO= @T+¢XWALL=-XT) #CUF
XMOX= XWALL+ (FT *OB0WL/NNCGZZ)-20
Dx= XMAY/INAXTL=1)
C CALCULATE DROPSTZE OISTRIBUTION
COLL SIZED COELPyRHOR, FCEN, SURFT, XMUF ,0N072Z,S274%,GDSZ,NSZG)
TF (TOTAIL.NE.0) WATTELEGTVEND
MO 3 J= 1, NSHEL
PTIM( 1, JD= O66
3 CONTINUE
(1) = XI
1 T=1,NAXIL

foawe

87


Te

THETEIVEG.CADS =245.00 PEATHDS 1.96666 30 VICHSD= 2.319FOO3 TIKD= So46REC02 UICALZGMD
DPOTHUATHsCADS 4eFSZESTI WPLEALZCQ) 2=2.1285-01 NF UOTHIGH/CANS 0.
GISHTHIGALZCAD2-4,357E-01 OIPNTHICAL/CAS-2.90SE-05 SVBNTHIGH/CAI= 0.
SUMCLICALD=-2,279E000 SUMFAIGHD= 9 SFONTHIGH/ZA)= De
UCHADICA/SECD 2 2.167E403 OVOTHICH3/CAD=-9,273E400 OCLES2=-8.230£-05 ITERPC= 1
WANDE GHD S.O38E-75 YNOFS 2,209F-75 XNOGECPPM)= 2.2010-b9 SOOTE(GHD= 8,160€-31

243196001
SHOOTH(GH/CA) © 0
SOBDTHIGH/G:

THETACIEG.CADs 135.60 PEATHD= 1.152660 VICM3DS 2.2136403 TAKDS S.642E802 ULCAL/GHD =-2,037E001
OPITHUATH/CAD= Ly,LGE “G2 WPKCAL/EAD=-3.U47E-01 MF IOTHIGH/CAD= 0. SHOOTHIGH/CA)=.0
DICFTHINBL/CAI=-3,9GKE-01 OOROTHICAL/CAI2-2.769£-05 SVHOTHIGH/CAI= 9. SORDTHIGH/CO)= 04
SUMCLACALD==6 6220696 SURFALGII= 0. SFNOTHIGHZEAD= 0.
USHEO(SM/SECD= Ze1ASECIS OVOTHICHS/CADS-1,.183E601 MELEAZ= 1.041€-04 ITERPC= 1
YANOEIGNIS 6.546C°71 YHOC= 2.657E-71 XNOECPPMD® 2,555E-65 SOOTE(GHI= 1.502E-29

TQCNEG.CADS 9125.60 — PEATMD= 242856690 VECHSP= 2.UASFOGS TEKI= 3.825602 UMGALSGMD=~-1.739E 401
QPATHCATMZCAD= L,423E202 WPLCAL/TAd EG 1326-61 NFINTHIGH/CAD= 0, SHOOTH(GH/CAD= 0
OGCeTH(CAL/CAD EH 3.257E-71  MOROTH(CAL/CAN=-2,598E-05 SVEDTHIGH/CAI= 0. ‘SNBOTHIGH/C!
SUITE (0402 2-9.5150002 SUMFAIGH)= 0. SFODTHIGH/CAD= Oe
UGHLEMCU/SER DS 26231C 00S OVOTHICHIZCAD s-1.425E401 MELEIZ= 2.171E-0% TTERPCS 1
WMNCE ISHS G.555F~67 YHOES 21,849F-67 XHOLCPPMD= 1.77KE-61 SOOTEIGMI= 2. 7198-28

UMETACIEG.CADS 9115,00  PEATMD= 164595050 VICNSD= 1.9296003 TEKI= 4,02Z2E402 ULGAL/GHI =-1 4116401
OPATHESTM/CEDS UASTESIZ WPECAL/CAD==5 WOIE-OL MF UNTNCGH/CAD= 0. SHODTHIGH/CAD= 0.
CICSIMICAL/CRES-2.759F=92 IAKNTNICAL/E MD =~2.2756£-05 SVANTHIGH/CAI= 0. SOOOTH(GH/CA)= OD.
MUL(SALIS<$.2565402 SUMFSIGM= Oy SFANINIGH/ZAT= 0.
UCNAPIGA/SFC IS 2.226423 OVOTMIGM3/CAD=-1,647E+01 DELEL2= 3.206¢-C4 ITERPC= 1
AG4)= 2.595F-63 YNOE= 1,054F-63 XNOECPPMDS 1,013E-57 SOOTEIGMI= 3.7562-27

THETSNDEG.CAI= ~105.06
SPDT NETH E
DCN TH ICEL

PUATMD= 1.691690 VICHS)= 1.755103 TIKI W240E 402 UICAL/GHD=~1.04GECOL
Ye 2,U97E<C2 WP(CAL/CAD=-7 .GOCE~G1 NFIOTHIGH/CAD= O. SMDDTHIGH/CA)= OQ.
APsH2.1G5E-92 DOROTHICAL/CAI=-2.07H5-05 SVEDTHIGH/CAI= J. SULOTHIGH/CAD= 0.
LeHaSeeli SUMFICCHIZ a. SFUOTHIGH/CAS= Oe
De 2.333053 OVOTMIGNS/CRd=-1.83BE 901 DILFIZ= &.SH3F-O% TTERPC= 1
GNDE 1.429E=59 YHOrs S.BG1F-h0 XNOE(PPND= 5,579F-54 SOOTFIGHD= 5,061E-26

THETBIDEG «CA

“95005 PERTH = 2.008F909 VICH3D= 165036403 TEKI= WobBSESO2 UICAL/GHD
UPOIACATN/COI= S.435E=-92 WPLCAL/CAY=~9.076E-91 OFIOTHIGM/CAI= 0.
DACPTACCOL/CAYZ-1,550E-01 MOPOTHICALZCA)=-1.640E-05 SVBDTHIGH/CAI= 0.
CUCCSLES-2 555451 SUMFALGHD= 0, SFANTHIGM/CAD= 0.
UCHEPISHASECDS 2.379E0)3 OVOTHIC43/CR)=-2,985E FUE DEL:
YONGE IGMS A .h92E-56 YNO

6.2396 000
SMNOTHIGM/CAD=
SOROTH(G H/C!

25 S.ut7i-de 1VERPCS 1
= 3.579E-56 YNOF(PPMD= 3,394E-50 SOOTE(GMI= 6.95CE-25

THETECUEG.TAD= 85.50 PEATNDS 2,454E 000 VICHSD= 16559005 TCKI= 4.768002 ULNALZGMD=-2,320E000
OPNTHCATsCAD= UOB57E=02 WPLCAL/TAY=-1,161F 060 UFININCCH/CAI= 0, SMDOTHESH/CAD= 0.
OWATAICAL ACA O2KE-02 ONGOTHICAL/CAD =-9.MBLE-06 SVBOTHIGH/CAd= O46 SODOTHESH/CAD= O04
SUMGLICELD S274 22051 SUMFSCGMD= 9, SFODTHIGH/CAD= 06
GOMIDEAMAGEOAS BeOLTEATT OYOTH(AMI/CAD=-2,076E00t OFLELZ= 6.501E-04 TTERPC= t
XANCE (GHD = G.UIGE-G2 HOES 2.637E-52 KNOELPPHD= 2.535c-46 SOOTECGHD= 1.007E-23

PHEE® 1, 468C-02
SMVOTHIGH/CAD = Co
SLOOTHIGN/CAS= 0.

=
WXNETICAL) =-1,069E690 aRtens-aTHsGHeK 2-800F 000

PHEES 2, WG8E-02
SMVDTHIGH/CAD= 0,
SLOOTHIGHZCADS Ue

WKNET (CAL) 2=3.879E000 RRICMS-AIM/GN-K)= 2.346E 096

PHEES 1, 4n8E-o2
SMVOTHIGH/CA)= &.
SLOGTHIGH/CAd= oO.

Oe
WKNETICAL) 2-7.727E+00 RRICHI-ATH/GH-K)= 2.540E696

PHEES 1, 4b82-02
SMVOTHIGH/GAI= Ce
SLOOTHIGH/CAI= Oe

WKNET (GALI == 21,2B3E602 RRUCHS~ATM/GM-KI= 2. a4CE OOS

PHEES 1,4682-22
SMVOTHIG¥/CA)= 0.
SLOOTH(GH/CAD= Ge

WKNET(GALI=~1.947E 401 RRICH3~ATM/GN-KI= 2,04GE400

PHEES 144585 -02
SMVOTH(GH/CAd= G6
SLODTHIGH/CAD= Oe

WKHETICALD =-2,602E901 RR(ICHS-ATH/GH-K I= 2.54GE 060

PHEES 1, 4680-62
SHVDTHIGH/CS1= 0.
SLOOTH(GH/CADE G,

HKNET(CALD=-3.895E001 RR(ICHS-ATM/GH-K}= 2,d40E000


IFRT=9

CALL EVAPIP,T, TA, TFL yO PA, OMEGA TIM, RHOAyXXXXyXMAFA,KOUNT CEVOP,
1 AIGND,X4SOTU)

XMAA= XMO-XMAFO

DFGOTH= .CXXXX-XMFGI/DTHETA
IF (XMAA.LT.O.9) XMAA=0.0
IF (XMAA.LE.0.9) KOUNT=NSS
IF (21et » GO TO 25
IGEND= IGEND+ 1
IF(IIOLO.c&Q.IF) GO TO 25

c INITIAL CONDITIONS FOR NEW ZONE

IGEND= 9

ITOLO1= TIOLO+ 1

99 21 I=II1o0L01,11

AMWBM(T )= amya

XMAPBI= XMBM(T)

DAPITH(I) = XMAPBI/OTHETA
xO2(I) = X924

XNONEQ(I )= XNOANE

YNO(T)= YNOA

XMNOCT J= XMASFDLFYNOCT)

TaM(I)= TA
UEM(I) = ua
H38M(T)= HA

CP84(T)I= CPA

PHE(T) = PHEA

ETAG= Led/(1.0+ AFRO)

ALPA= 4 ETaAa

XMFBM (LT XMBM CT) FET AG
XHPaM(T XM3i4(T) FAL PA
XMFLICID= OLOOTHCT)*OTHETS
BTASM(IT )= BTAG*(XMBM(IT )/X4A)
GEM(I d= GA*(XMEMUT 27X44) **0.5666667
BTAA= BTAA- 3TABN (I)

Gh= GA-G8N(1)

YSCE)= XMSOOTCIIV/XMBM(T)

suas SMOPB+ XMBM(T)

SLOOTH= SLOOTH+ OLINTH{T)
SY80TH= SVBDTHt DV30TH(I)

IF (IBUG.EQ.0) SO TO 21

WRITE (6,111) [Ty AF ROL)» XMFG,PHE (I), TSM(L) ,USM(T) ,RBM(T) ,CPSN(T)
1 ,OVBOTH(T),XMFLICI) .
WEITE (6,112)

69 20 L=1,NS2G

WRITE (6,169) Ly DICLT gLd,XNOT(T 4L),O0MFLICT 4b)

20 CONTINUE
24 CONTINUE

SUMBM= SUMBM+ STMQP3

XMAS XMO- SMAPB

IF (XMA.LT.G.9) XMO=6.0

XMAAS XMAB- SMAPB

IF (XMAS.LT.O.0) XMAG=0.9

IF (X#4A.LE.9.0) KOUNT= NSS
XMFA= XMAFET AA

XMPAN= XMAFALOA

XMNOA= XMA#YNOG

XMSOTA= XMOFYSA

OPB30TH= SMAPS/OTHETA

$1


TTNITL
THETAT
THCMAX
P

T

v

UA

HA
XMAO
AMNA
FARA
PHEA
FAROS
xMEGO
XMAEGR
XMFEGR
XMRO
XMARD
XMF 2D
AMNE
OXxsT
xO24
FOEN
HHV
HCOM
HEFL
HFS
vo
vors
VF JO
SEND

0

ow uw ow 8

-O.15E +03,

0.3€+02,

O.1C+ot,
0.333554271A8909E+03,
0,23621320794931E +04,
-0.24477070064609E +02,
-0,11905387510382E+01,
0.24536R63750245E +01,
9.2835G04302659RF +02,
3.162556003768316-92,
0.1487747340597S5E-O14
0.69873063933353E-01,
O.Cy

O6Ce

O.c,
0663086375924526E-01,
0.6%552920023545E-01,
0,.2523455000981€-02,
0696E+G2y,
0.33202666666667E +014
0.25853426216717E+0Cy
G.85193285988741E +00,
0.1696S555555556E +054
0.10269305555556E +05,
70264398 +03,
70.2313562613658E+03,.
3.18599813601726E403,
0.233597017627615 +04,
0.20%3994172532RE +055

30


*DECK FLAME7
SUBROUTINE FLAMEZ(TLsTLC,CPL,CNOP)
C COMPUTATION OF GUASI-STEANY THIN FLAME STSUCTURE (WITHIN A FINITE SPACE)
COMMON/INTEGSS INOL,ZINC2,ISOCT,IWETBT »NPRTsICTAILsNTERsNAXIL,
NSHEL,NSZG,NITST,IGUG,II,III
{MON/OCNYF/ XJyATO,CTRAXKySIG,R,RU,CTCTH,CTHOT yOTROT yPI
COMMCN/SPECYIZ = XH20100) oe XH2C (106) »XCCC1CO)4XCO20109) yXN2(100)5

1 XNOC100) oXO4(100) XH (100) ,XC (100) XN(100),xX02¢100)
2 »XCH4 (100) ,XNCNEQ(100)
COMMCN/PRATET/ DNOOT1 (100) ,ONCOT2(100) ,OMSOT1(100) ,OMSDT2(100),
1 ORMOTH(100),DL80TH(100),XMOF (100)
COMMCN/TPRCPIY PyTOM(100) ,UBM(100) pHEMCICG) »VEM(100) sCPBM(100),
1 AMWBM (100),R936(100),PHE (100) ,AFR(100),XMFBM(100),
2 XMPAM (100)

4 COMMON/SEAL THEMAX, THE TJy THEIGsDINJy THERE, THEVO,THEVCSDTINGy
CTHMIN,DTEMED, OTHMAX sOTHE TA,CT
* conmonso20rcy HUTIM,OTPHY,XTLYC2,CPBARyRHOD,TCRITF yXL2985
4 a XMBM (100) ,CVEDTHI10G) »PLONTH(1609 ,XMFLI(1000,4
2 OLC100,3) ,0MFLI (10052) sCMLOT(100,3) sX¥NCI (10093) 5
3 XKK (100,3)
COMMON/SEUELP/W yg XX 4 VF Xs ZFLyYSTCA,ZSTOA,API yFEEN yFAPOS yAMWF, SURFT
1 aXMUF yCEN,OXST»HCOM
COMMGN/SFLAMES/ AMWMPyRMAX,YN2L TF a RF y XMOCT RL» RB, TDF (50) »YFDF (50) y
1 YO2DF(50),YN2NF (50), YNCOF (56) ,WOCT(S50) »ROF (50)
CCP= HCOM/CFEAR
NITST1= NITST- 1
NITST NITST/S
NITSTZ= NITST?2# 2
00 20 I=1,II
XMOF(T)= 6.0
CLENTHITI= 9.9
ONCNTZ(T)= 9.9
NO 19 L=1,NS2G6
XKKIT,LI= 0,9
SOMLOT(I,LI= 9.0
19 CONTINUE
CLEIT= 0.9
OMLOTT= 0.9
IF (XMFLI(I)LE.0.0) GO TO 20
YO28= XO2U1)*32./AMWBM(T)
YN23= XN2C1)*28./AMWRM(TD
YNC3= XNCNEQ(I) #30./AMWBM(T)
Tes TEM(T)
AMWMP= AMWAM(T) #P
YO28J= YC28/0XST
XYXMOF= 12.56637062*AMWMP/82.054
© PRORLET SPACING
SUMNOT= 0.0
00 9 L=1,NS7E
IF (OT(T.L).GT20.0) SUMNOT= SUMNDI+ XNDICT4L)
CONTINUE
VMAX= VBMCT)/SUMNOT
QmMax= (C.23872*VMAX) ¥# C0, 333332
CO 10 L=1,NS7G
TF (OIUI,LY.LE}M.0) GO TO 10
RL= O.5*O0T(TsL)
IF (Y¥C28.LTAXTLYO2) GO TO 12
. IF (CCTCOTTE-TLY.LE1.0E-€) GC TO 15
G SUBCRITICAL ASSUMFTION
XLS XL298* (CTCRITF=TL IZ CTCRITF-298.9) © *0.38

Ro)

70


SUBROUTINE RADH2Z (RADHZ2)

This subroutine evaluates the zonal rate of radiative heat transfer.

Subroutine RADHZ is used for that of the "air" zone.

SUBROUTINE SIZED

This subroutine determines the drop-size distribution of the fuel jet
by an upper-limit distribution function based upon the Weber number and the
Reynolds number of the jet.

SUBROUTINE SOOTRZ(SOOTR)

This subroutine calculates the zonal soot formation rate both prior
to and after ignition. It also determines the soot oxidation rate based upon the
zonal pressure, temperature and equivalence ratio. Subroutine SCOTR is used for
that of the "air" zone.

SUBROUTINE SPECYZ (SPECY) :

Zonal equilibrium species distribution is obtained via table look up
for specified independent variables P, T, and ¢%. Subroutine TRIPIT is called
to perform the three-dimensional linear interpolation. An error message will
be printed out whenever any independent variable lies outside the range of
the table. Subroutine SPECY is used for that of the "air" zone.

+ SUBROUTINE SPRAY
This subroutine calls subroutine SIZED to determine the drop-size
distribution of the fuel spray. The fuel jet is simulated by a steady jet

(free and.wall jets) theory. Therefore, the trajectories and over-all air-fuel
ratios of fuel elements are defined in Eulerian spatial coordinates.

SUBROUTINE TAPEIN

This subroutine is called at the very beginning of MAIN to the input
species distribution tables from TAPE 8.


V. SAMPLE CASE

The program input and output of a sample case is presented in the
following pages to facilitate program checkout.

21


X= XT#(I-1) "0X
TF(X.GTeXWALL) GO TO 8
C FREE JET REGTON
BI+CUF*(X-XT)
XJF/(RHOR* BFR FXIF3)
c -EBFVF JO
XXX¥= BBFPO=4 0 *CC
VCL= 2.5*(-S38+SORT(XXX))
IF (VCL.GT. VFJO) VCL= VFUC
RHOFCL= XJF/{(VCL*® B¥B*XIF3)
AFCL= RHOR/RHOFOL
60 TO &
® CONTINUE ‘
IF (€X-0X9.GT.XWALL) GO TO 12
C WALL TRANSIENT ZCNE
BW= 68
xO= X=-0X
RHOCLP= XI2*xXJF/(BWFAWFVCL*XIF2*X13)
XXX= XT2% (RHORFXIFE+RHOPCL*XIF7) J CXITF2* (RHOR FXIG+RHOCLP*XI7))
VOLP= VCL*SGRT XXX)
BOP= XIF2*BWFVCL/(XT2*VCLF)
XJMP= VCLPFVCLPFBWFBOPF(RHCRFXIGtRHOCLPFXISY
IF (T@UG.SEG.0) GO TO 14
WRITE (6,113)
WPITE (6.112) 8yRHOFCL,VEL,XUM
WRITE (6,113) 80P,RHOCLP,VCLP.XJMP
14 CONTINUE
vcl= VCLP
= BCP
RHOFCL= RHOCLF
f WALL SEGICN
12 CONTINUE
RX= X=xX0
B= 29 + CWFFRX
AR= (RX+AW) FBFRHOR*X IL
P= XIS/XIS*XUF
CC=-XJMP
XXX= BB*BR- 4.0*AAFCC
VCL= -BRtSORT (XXX) /(2.07AA)
RHOFCL= XJF/IVCL¥(RX4BW) FB*X13)
AFCL= RHOR/RHCFCL
Cc WPITE (6,104) [,85XP29XN1UsXN2ZU sXMMoUST,VCL,RHCFCL yAFCL
& CONTINUE
ETAICT,19= 0.0
ETATCI,NSHEL+19= 1.0
ETA=0.1
00 2 J=1,NSHEL
TF (X.GT.. XWALLY GO TO 7
CG FREE JET REGION
IF (JeEQsNSHFLY GO TO 16
GN=(FLOAT (YS) VFOUF/ (2.0 7PTFEHOFCLFVCL*R*O)
Oo 5 ITEP=1,20
ETAL= G.5*(ETAFETA)
ETA?= 0.85714*(ETA**3.5)
ETAZ= 0.67 CF TASS)
ETAG= G.15384*(ETA**6.5)
G= ETAL-ETA24ETASZ-ETAG -f0
CHK= ABS(EsCR)
IF (CRK.LY.1.06-5 2 GO TO E
ETA32= FTA**1.5

88


2

oO

fs!

IF (Ispcy?.€9.9) GO To 12
WRITE (6,110)
TF (XMA.GTeGed) WRITE (6,109) IID, XH2AyXH20HyXCOA,XCO2MyXN24,XNOAy
1 XOHA,XHQ,XOA,XNA,XO24, XCHUA
no 11 Is1,11
WRITE (6,199) T,XH2(T) XH2011),XCOCT) »XCO2Z(T) »XN2(¢1),XNOC(I) »XOHEI),
4 XHCT)  XO(I) sXNCT) X02 (1) »XCH4(T)
11 CONTINUE
12 CONTINUES
TF (IFATES.ER.0) GO TO 54
WRITE (6,114)
OO 14 I= 1,1T
WMQITE(6,199) LT, OTSOTHCI) sOMDOTH(I) sOMVOTH(I) ,OVBOTH(T) »CLBOTH(T)
1 ,ONODTH( I), ONOITL(T), ONOOT2 (1), D0MSOTH(T) ,OMSOTL(1),0MSOT2(1)
14 CONTINUE
S51 CONTINUS
IF CINROPP.EN.9) GO TO 77
WRITE (6,115) (LyL=1sNSZG)
00 15 I=1,11
WEITE (6,109) Ty, (OIL LI, XNOT(TyL) sOMFLI(ITyL) »IMLOT(I5t) pL=tyNSZ5)
15 CONTINUE
77 CONTINUE
IF (IGEND.GE.3) KOUNT= NSS
TF (KOUNT.EQ. NSS) KOUNT= KOUNT#1
TPRT=IPRT +4
TF (THETAI.LT. THEMAX) GO TO 2

: INDICATED ENGINE PERFORMANCE

VorS &,102366-24V0IS
. 37 .0SS75*SUMNK/VNTS1/2
MEP*VDTS1*RPM/ (395090. *NSTROK/2.)

60 o#XMF J FRPM/ (XTHP FNS TPOKS 2.) /453.59

xTSFC

~ SQOT CONCENTPATION (MG/M3STP)

PSTP= 1.56

TSTA= 390.

VSTP= (PFV/T)*CTSTP/PSTP)

XISOOT= SOOTE*1.UE3/ (VSTP*1.,6E-5)

WRITE (65116) XIMEP,XIHPs XISFC,XISOOTsXNOC

/PROGENE TO THE NEXT CASE

50 TO 1
FOOMAT STATEMENTS
97 FORMAT (1X,*ITESATION NOES NOT CONVERGE TO THE CORRECT TEMPFs
1 215,1P7F10.37)
98 FORMAT (H1)
99 FORMAT (8419)
199 FORMAT(SE10.3)
101 FOR“ATI//1X%y FTHETAIDEG CAI=*#,F 8.2, 3X, *P(ATMI=*,4PE10.352X5 FVICM3)
LEF,210 03, 2X, FT (KIS, F106 5, 2%» FU(CAL/GM) =*,E10.3,2X 5 *PHEE=*,E19.3/
22XyFOPOTH COTM/C AI =F, E1063, 2X%y FWO(CAL/CA) =*#,E106352Xy FOF UNTH(GM/CAD
S=%,E1063, 2X5 FSH OOTH(GM/CAV=* 015 6352%5 F#SMVOTHIGM/CAD =F, E19. 3/
42X_*DOCDTH(CAL/ CAD HF yE1062,2X5 FGQROTH(CAL/CA) =*,£106342X5 *SVBOTHIG
SM/CAI =F, £10. 3,2Xe *SOBDTH(GM/CA)=*, E1063 52Xy*SLOOTHIGM/CA) =F,
6 £10. 27/2X 5 *SUMILICAL) =*%, £10.55 2% 9* SUMFS (GM) =*,E10.3,2Xy*SFBOTHIGH/
TORY HF, ELV. Fy 2X4 FWKNET (CAL) = 4F 10 652% 5 FRRICMI-ATM/GM-K) =*, 610.37
32% yFUCHS® (CM/SESI=*%,E19.3,2% ,*OVOTHICH3/CAI=*, E1023,
D2Xs*OELSL2=* Eid. 3,2X, FT TESOCH*,IS/
F245 XMNOE (GMI=* E106 35 2X5 FYNOES* 9015 0 342K, FXNIE(PPMI=F E1563, 2X5
1 *SOOTE(GMI=*,E19.3)

102? FORMAT U/L Xs FSYSTEM* yk y*MASS (GM) Fo 4X9 FT (CK) #y 8X, FUTICAL/SGM) * 4 2X,
L *PHE* 19X_e* BETA*S7Xy*XMDF (GM) F y4Xy ¥XNONEQ* 9 5X_*XMSOOT (GM) *42Xy

2 FXMFLICGMDF 2X FA TRIGM) #,1X,* FUEL (GM) *)

61


FOECK SPECYZ

SUBROUTINE SFECYZ

© TRIS FCUTINE CALCULATES THE FCUTLIQ°IUM SFECIES DISTRIBUTICN OF A GIVEN
C STATE I.€. (PsT,FHE) FCR A SPFCIFIC FUEL FROM TRIFLE TASLE INTERPCLATION

1

1
2

1

1

nN

Ae

t
?

COMMCMZINTEGRS INQL,INC2+ISOCT,IWETBT»NFEFT,ICTAILeNTERsNAXIL,
NSHEL4NSZG,NITST,IFUG,II,III

COMMCN/SPECYIZ = XH20100) »XF20(100) sXCCC1060)exXCO20100) sXN201 000,
XNC (100) »XCH(100) ,XH(109) ,XC (100) 4XN(1009,X02(100)

+ XCHG(100) XKCNECCIOC)

COMMON/PRATETZ ONOOTL (106) ,ONNOT2] (100) ,OMSOT1(100) ,OMSOT2{(1C0),
ORMOTH(100),0LR0TH(160),xMDF (100)

COMMCN/STSROPIS ©, TBM (100) -USM(106) »HEM (100) »VEM(100) »CPBM(100),
AMWBM(100),R3M(100),PKE (100) ,AFR(100) »XMFOM(100),
XMPAM{100)

“COMMCN/SPTAB/ NETRNTTByNPHTE,PTAC7) pVTE(13) sPHETBU 11)
XH2TBC7 913511) pXH2CTE (7513511) eXCOTS( 713,11) 6
XCO2TBC7,12411) sXN2TBI7 413411) sXNCTEC7 413511) 4
XOHTR (7413511) oXHTE(7413511) pXCTS(7 413,11) 5
XNTSC7 5134410 eXC2TBRC7y12y 42) sXCR4TIU7y 13511)

CC¥MCN/YXCON/ TOP,IOT,IOPHE,

TLP1,IOT1, TOSH1, OFLF,CELT,CELFRE,CP1,0T1,0°HE 1,
0P2,0T2,0PHEZ

Fis P

O06 1 IQP=2yNFT2

TF(P1.GE.FTSCIOP=1) sAND. PL.LLE.PTBC(ICP)) GC TO 2

CCNTIAUE

WPOITE(6,100) FL

FOOMOT(* ERROR eeee F CUT OF RANGE *, 1PE10.3)

IF(PA.LLT.PTO2(1)) GO TO 7

Pi= PTB(NFTR)

InP= KPT

GO TO 2

CONTINUE

Piz PTR(L)

Taps 2

CONTINUE

de 10 [=1,11

Ti= TeM¢(L)

GPREL= PHEIT)

QO 3 IOT=e@,NTTa

IF (T4.5F.TYA(IOT=-1) «AND. TI.LE~TTB(ICT)) GCG TO 4

SONTINUS

WOTTEC6.104) Tt

FORMAT(® FRROR seeese T OUT OF RANGE * 41FF1IC.3)

IF (T1.LT. TIP(1)) GO TO eB

Tis TIBCNTTA)

TOT= NTT?

Go TOG

CONTINUE

Tis TTeR(1)

Trt= 2

CONTINUE

D0 S TOPHE= 2, NPHTR

TF (PRE 4.GE.PHETRCIDPHE~-1) .AND. PHELLE*PHETECIOPHE)) GO TO E

CONTINUE

WEITE(6,1092) PHEL

FORMAT(* FRECR seeee PHE CUY-OF FANGE *,1PE1C.3)

If (PRE1.LT.PRET@ (12) GO To ¢

PHEL= PHETRCNPHTAD

TOPHE= MPHTS

85


TECHNICAL REPORT DATA

read Instructions on the reverse before completing)

aa

(Please
7. REPORT NO. 2

EPA-460/3-76-008-b

[3. RECIPIENT'S ACCESSIONNO.

4, TITLE ANO SUBTITLE
Study on Oxides of Nitrogen And Carbon Formation

5. REPORT DATE

‘Issued May 1976

In Diesel Engines - Computer Program User's Manual

l6. PERFORMING ORGANIZATION CODE

7. AUTHOR(S)

C. J. Kau and T. J. Tyson

ja. PERFORMING ORGANIZATION REPORT NO.

19. PERFORMING ORGANIZATION NAME ANO ADORESS
Ultrasystems, Inc.

2400 Michelson Drive

Irvine, California 92715

10. PROGRAM ELEMENT NO,

TT CONTRA [ANT NO.

68-01-0436

12, SPONSORING AGENCY NAME AND AOORESS
U.S. Environmental Protection Agency

"Finat Repore” sue’ "732.50Ne"*7q

Office of Air and Waste Management
Office of Mobile Source Air Pollution Control
Ann Arbor, Michigan 48105

114, SPONSORING AGENCY CODE

15, SUPPLEMENTARY NOTES
Co-sponsor:

New York, 10020
[ie AastTRAGT

Coordinating Research Council, 30 Rockfeller Plaza, New York,

16. AGSTAA

tion diesel engines has been developed and tested. The

which control fuel/air mixing.

The model has been tested against results obtained with
engine.
meters on NO, emissions were obtained.
not satisfactory.

A mathematical model describing heat release and pollutant formation in direct injec-

rical constants which can be tuned to fit the requirements of a particular engine.
Sensitivity studies indicate that the model is most responsive to those constants
Numerical experiments strongly suggest that diffusion
flames modelled by spherical droplet flames are unsuitable for this type of system.

Reasonable predictions of the influence of engine design and operation para-
However, predictions of smoke emissions were

model includes several empi-

a single cylinder diesel

7. KEY WORDS AND DOCUMENT ANALYSIS:

IN ENDED TERMS |c. COSATI Field/Group

Diesel Engines
Combustion
Emission

Nitric Oxide (NO)
Nitrogen Oxides

IBUTION STATEMENT

Computer Program
Fuel Sprays
Internal Combustion
Engine

Fuel Consumption

OESCRIPTORS: >. IOENTIFIERS/OPE!
Diffusion Flames Mathematical Modeling Mobile Sources
Air Pollution Soot

Exhaust Gas Recirculation
Divided Chamber Engine

19. SECURITY CLASS (This Report)

21-07 (Recipro-
cating Engines)
21-02 (Combus-
tion)

‘BT, NO. OF PAGES
106

Release Unlimited

20. SECURITY GLASS (This page)

22. PRICE

EPA Form 2220-1 (9-73)

100:


FQECK TAPTIN
SUBROUTINE TSPEIN
C SPESTES TABLE INPUT FROM A TAPE
COMMCN/SPTAB/ NPTBNTTA,NFHTO,PTEC7) »ITSC12) ,PRETBI 11),

* KMH2TAC7 yA %e14) 4 XH2CTAC7,1 34114) pXCOTAC7e13,11)5
1 XCO2TBC7 a1 7011) os XN2TSIC7,13,11) oXNOTO (7413551175
2 XOHTBC7, 12,11) ,XHTRC7o12511),XCTRC7,13011)4

z XNTBE7,13,21),xXO2TB(7412411) »XCHUTI(7,13, 11)
NPT?2= 7

NTTS9= 13

MPHTE= 11

C SET ZER0

NO 84 T=1,12043

PTR(T)= 0.9
8G CONTINUE

DC 31 IPHTS=1, NPHTA

NG $1 TPTR=1,NPTB

READ (8,121) PTBCIPTB) »PHETS(TPHTE)

2EAI 18,122) (TTB(I),1=1,NTTED

TE (OHETS(TOHTS) LE60.05) GC TO 82

PEA9 (8,122) (XH2TACTETByIyTPHTS) yT=1,NTTS)
PEAD (8,122) (X4ZOTACTPTS,I, 1PHTB) ,I=tyNTTAD
BEAD (8,122) € XCOTRCTPTS, 1, TPHTE) yI=1,NTTED
DO (8,122) (XCO2T3CTPTEy Ty IPHTO) yI=1yNTT2)
a? PONT IMU

FFA (84122) ( XN2TS(IPTA,T,IPHTB) yI=t,NTTV

- PEAD (8,122) ( XNOTA(TPTA,I,TSHTB) ,I=1,NTTED

TF(OHETPCTPHTA) LE.0.05) GC TO 83

269 (8,127) ( XYOHTBCIPTS,I,ICHTE) »I=1,NTTSD

REAN (8,122) ( YHTSCIPTA,I,IPHTB) ,I=1,NTTE)
8% CONTINUE

FLAN (8,122) € XOTM(TPTE,I,IPHTB) gI=19NTTP

2FAN (8,122) ( XNTBCTPTS,Ly1PHTB8),I=4,NTTE)

PES (8,122) € XO2ZTRCIPTE,I,IPHTE) ,I=1,NTTE

TF(OHEFECIPHTA) LE.0.05) GO TO 34

PEND (8,122) (XCHTS¢TPTA,I,IFHTE) »T=1yNTT2)

81 CONTTINUF
171 FORMAT (1¥,1P3E10.3)
122 FORMAT(HX,1P13E10.3
PETUPN
ENO

gl


FCECK TEIPIT
SUPROUTINE TRIPIT{(G,QT)

G TRIPLE INTERPOLATION CF A TABLE FUNCTION Q«P,T,PRKE)
DIMENSION G(7,13,11)
COMMCN/XXCON/ IOP,IOT,IOPHES

1 LOP1,INT1,INPH1,0ELP »DELT »DELFHEsDP1»y0T1,0PHE1,
2 0P2,0T2,0PHE2

QJLiKi= Q€ IOP, TIDT,IOPHE) *OP1+ QC ICP1,1I0T,TOPHE) *0P2

Quik = Q¢ IOP,IOT,IOPH1)#0°1+ QC IOFL,1IOT,IOPH1) ¥DP2

CUKt = QC TOP,INT1,ICPHE) *OP1+ GC L[OP1,IOT1»IOPHE)*OP2

QuK = Q¢( TOP,INT1,IDPH1)*#0P1+ Q( ICP1,I0T1,1IDPH1) *0P2

G2= GILKLFDT1+ QIK1FDT2
C1i= GJLKFOTL+ QUK*OT2

or = C2*DPHE1* Q1L*O0PHE2
RETURN

END

*DECK VCLM
SUBRCUTINE VOLM(THETAT, OSOTH,CVDTH.V)
C#¥Fee THIS ROUTINE CALCULATES THE TOTAL VOLUME AT CURRENT CA
COMMONSENGEOZS DESRFWeRLAVG_SCyS,UST,CReRPM,TWyONOZZyNNOZZ,0BOWL,
4 VFJO
CCMMON/LONVEZS XJeATD OTR eXKySIG,R»RU,CTOTH,CTEOT ,OTROT»PI
XxX2= THETAT*#OTR
SINT=SIN(XX2)
COST= COS (XX2)
RLR= RFW/RL
XXX= 1. 0-€RUR*SINT) ¥¥*2
XXX= SQRTCXXX)
= RFWH(1.0-COST #(41.0-XXXI/RLR)
OSDTH= @FW *SINT *£1.0+CCST#RLR/XXX) = FNTR
OVITH= 0.7853981E35* DE *DE*OSOTH
V= VC+ O.7A5S39B1LES5*DE FOES
RETURN __
FNO

98


7
€ cece

2
1

15

39
102
101
102
163
104
105
106
108

CONTINUE
LECTING THE SURNING SYSTEMS (I.E. DEFINING BURNING RATED
IDIS(1,wW=3

IIs I+

KOQUNT= KOUNT #1

XMSCOT(TI)= XMSOT(I,J) #NNOZZ

XMBM(TI)S NNCZZ*AFTIFOMFIJT

AVBOTECIID= NNOZZ*XMEVS (Ty J) FOMFIUT/OTHETA

XMFLICIID= 0.0

99 10 K= 1,NSZ26

OI(LI,K)= Olly JeK)

SMFLICII,K= 0.9

XNDIC(IIyK)= 0.9

IF (DICIIgK).LE.O.0) GO TO 10

DMFLI(II,K¥= NNOZZ* (DI (114K) /SZBAP (K)) ¥*3*GOSZ(K) FOMFIUT
XNOECIIgK)= €-O7DMFLICIT, KY/CFDENSPI*COI(IT, Kd#e3))
XMFLICIID= XMFLICIIO+ OMFLICII,K)

CONTINUE

CLOOTHCIIV= XMFLICITI/OTHETS

IF (IOTAIL.NE.O) WRITE(6,101)
1 Tod REV (Td) sXMFVS (Ty J) »XMOMCIIDy TIMI) DIG) XMSOT(T J)
CONTINUE

CONTINUE

XMFG= NNOZZ*XMFGL

XMAFA= NNCZ7*XMAFAL

XMSOTUS NNOZZ*XMSOT1

TIM= TIM#+CT

IF (18UG.EG.9 } RETUPN

NO 15 J=1,NSHEL

WETTE (6,99)

WRITE (65102) ( XSSCTII-T+1,J) ,I=1,1I1)

NPITE (65103) C YSSCIII-I+t1,J),l=1,110)

WRITE (65109) COYSSCIII-L4+ty,J) sI=1,1I1

WRITE(6,194) (XMFVS(III-I+1,u),I=t,111)

WRITE (6,105) ( AFV(III-I#1,J),T=1,11I)

WRITE (65108) CTAFV(III-I+t1,J),T=1,TIL)

CONTINUE

WRITE (6599)

FORMAT(/)

FORMATCAX,* T=*,T3,* J=*,I3," ELEMENT TMPINGES A NEGHBOR JET*
FORMAT (1X_¥%----- SYSTEM %,213,% IS BUFNNEC *,1P6E12.3)
FORMAT(1X,* X *,1P10E12.23),

FORMAT(IX,® Y ¥%,1P10612.3

FORMAT(AX,*X¥FV *,1P10E12.3

FORMATCIX,*4/FV #,4P10E12.2

FORMAT(1X,*TAVF *,1P10612.3)

FORMAT(1X,* DY *,1°10£12.3)

RETURN

ENC

68


VI, PROGRAM STATISTICS

The computer program is coded in standard FORTRAN IV, and the
version is currently compiled and optimized in FIN (FORTRAN Extended)
mode and ran ona CDC 6600 computer. The computer statistics (based on the

sample case presented in preceding section ) are shown in the Table below.

Compilation Time 71 sec.
Run Time 128 sec.
Core Used 120 K (octal)
Tape Required 1 (7 track)
Compilation I/O Time 98 sec.
Run I/O Time 6 sec.

39


FCECK SOOTR
SUBROUTINE SCOTP(PyTBReAFVeYOZFyV oT, ILIGIsXMSCT »VSMy DMSDTs
4 mMsoTia,oMsot2)
© THIS RCUTINF ESTIMATE THE SCCT FCRMATICN ANC DESTSUCTION FATE

4 AMWA gSFOA,FATA,PHEA
COMMON SFUELP/Wy XX yg YF Ly ZFL ay YSTCAs ZSTOA SAFI yFCENSFAROS yAMNF, SURFT
4 »XMUF CEN, CXSTHOOM .

SOMMCH/D2CPO/Z HUTIM,OTFERY,XTLYC2 sCPRAR,RHOOSTCRITF yXL298,
XMBM (109) ,0VANTH(109) ,CLOOTH (109) »XMFLIC1G0)4
O1€109,3) ,OKFLI (100, 2) pOMLOT (100039 4XNOT(100,3),
XKK (19042)
COMMCN /SCOTPY SOEN,GSCOT, ASCTP,ASOTAsXMSCCT (100
IFCIOIGT.GT.2) GO TO 2
& VAPORIZING ZONE
1 CONTINUE
PFuS F/(1.0+4FV)
PREEV= PHEA*(1,0+(1, 2+AFROAI/AFY)
GMSDT= ASCTPF(YBM*VAM/VIF(T/TAY *(PHEEV**3) FPFUFEXP(-20C00/T®)
RETUPN
fC RUCNING ZONE
2 CONTINUE
CG SCOT FGQMATTON RATE
PREB= 1.0/(AFVFFARGS)
DMSNT1L= ASOTA*(VAMFVSM/VIFCT/T 3) F(PHFS*42) FP FEXP(-20090./T8)
G SOOT CXIQIZATION PATE

ane

DOMMON SAIRPLS XMATyXME Us XMEGRSXMPD yXMAOsXMARC »XMAEGRyXMFROyXMFEGR

DMSOT2= &.51E4* (XMSOTFYOZB*O J /(SCEN*OSCOT*SERT(T8)) “EXE (-1956500/

4 78)
© NET FECOUCTION RATE
OMSOT4- OMSOT2

81


VII. REFERENCES

Kau, C.J., Tyson, T.J. and Heap, M.P. "Study of the Oxides
of Nitrogen and Carbon Formation in Diesel Engines". Final
Report Phase II and III, Ultrasystems, Inc., June, 1975.

Gordon, S. and McBride, B.J. "Computer Program for Calculation of
Complex Chemical Equilibrium Compositions,. Rocket Performance,
Incident and Reflected Shocks, and Chapman-Jouguet Detonation",
NASA SP-273, NASA Lewis Research Center, 1971.

41


31

7

REF

XMFLB= XMFLO- XMFOO

XMFG = XMFV- XMFPB

TF (XMAWLT.C41.0E-5*XMA0)) XMA= 0.9
IF (XMPL.LT.9.0) XMFL= 0.0

IF (XMFG.LT. (1.08 -S5*XMFJD) XMFG= 0.0
Ga= 0.0

GFG= 6
BTAA= 9
8TOFG=
IF (IIeLE.9) GO TO 29
OO 30 I=1, II
BTASM(T)=0.9

GEM(I)= 0.9

CONTINUE

CONTINUE

SUMR=0.€

SUMG=6
SUMU=C.%

SUMRT=0.8

IF(XMG.LE.9.9) S50 TO SL

Tas TAO# COTADTI+OTAODTH) #OTHET2

XMNOA= XMNOAOt (ANON TO#ANOOTH) #OTHET2

IF (XHNOA.LT.G.0) KMNOA= 6.9

XMSOTM= XMSTA0+ (AMSOTOFAMSOTH) FOTHET2

IF (TA.LT.TRONGL WOR. TA.GT.TRANG2) WRITE(5,10938) THETAI,TA,PHEA,P,
4 CTADTH,DTADTO, TAO

IF (PHEA.LLT.PHEPGL .OR. PHEA.GT.PHERG2) WRITE(5,168) THETAI,TA,P,
1 PHEA,XMFO,XMPAA

YNOQ= XMNOA/XMA

XNOANE= (AMWA/ 30.) *YNOA

YSQ= O69

IF (XMO8,6T.9.90) YSA= XMSOTS/XMAA

QTAQ=XMAF RASTA

G (TASTH) FCSTAR **(.6666567))

SUMS= SUMS+ GA

SUM2= SUM2+ XMAFRO

SUMU= SUMU+ XM4A*YU4a

SUMRT= SUMRT+ OTAA

SCNTINUE

IF (XMFG.LE.G.0) GO TO 33

TFG= TFGO+ (TFSOTOFTFGOTH) FD THET2?)

UFG= UFLBP+ (TFG-TEPF)*CVFYV

IF (TFG.LT.U.9) WRITE (6,100) TFG, TFGOTH, AAFGs3BFG,0P0TH, OF VOTH
1 aXMFG

RT AFG=XMF G*REGFTFG

GFG= (TFG-TW)*(3TOFG** (, 6666667) )

SUMS= SUNG+ GFS

SUMR= SUMR+ XMFG#RFG

SUMU= SUMU+ XMFSPUFG

SUMRT= SUMRT+ STAFG

CONTINUE
SUMU= SUMU+ XNFL*UFL+ XMFLB*UFLB

IF (II.LT.1) GO TO 73

oe 72 I=1,1I

IF (I.GT.ITIOLD) GO TO 67

TAM(I) = TSMO(T) + (OTGOTOCI) +0T30TH(I)) FOTHET2
50 TC
SONTINUE
INE THE TEMPEQTURRE CALCULATION FOR NENLY IGNITEO SYSTEMS
QAMAVX= 9.5% (93"(1)+03M0(T))

57


SEPARAM

cA = 0.26€00,
ce = O.75E+00,
NOR = 0,537E+01,

cvior = G.1E+00,
CMIx = 0.78E-02»
SOEN = O.2E+014

ascot = 0.1€-03,
ASOTA = 0.352E+0ty
ASOTP = C.352E +01,
HUTIM = 0.2€-33,
OTPHY = 9.2&-92,
SEFU = O.435+Ct,

COCHG = O.1E+%1,

CEvap = G.276€400,
AIGND = 0.146€-03,
RHID = (0, 3E=03,
CPBAR = 0.35E400y
XTLYOZ = 0.2€-%1,
SEND -

27


IF (IIT.EQ.1)
IXX=IPRT/NPRT
IXX=1XX#NPRT
IF € IPRT.NE.IXX ) GO TO 77
19 CONTINUES
PRINT OUTPUT
CVER ALL ENERGY BALANCE
Eis SuMU-£O
€2= HF U*SUMFJ- WKNET= SUMOL
DeLEL12= (E1- £2)/(XMFS*HCOM)
C AVERAGED POLLUTANT FORMATION
SUMMW= XMAFAMWA+ XMFG*AMWF
SUMOAE= XMPAA
SUMPFE= XMFO+ XMFL# XMFLB+ XMFG
IF (II.LE.0) GO TO 23
06 16 I= 1,11
SUMMW= SUMMW+ XM3M(T) *#AMWBMEI)
SUMPAE= SUMPAE+ XMPAM(T)
SUMPFE= SUMPFE+ XMFBM(T)

GO To 19

aoa

16 CONTINUE
23 CONTINUE
PHEE= SUMPFE/SUMPAE/FAROS
AMWE= SUMMW/SUMM
XMNOE/ SUMM
XNOE= (AMWE/30.)*YNOE*1.056
C GROSS OUTPUT
WPS. = 0.57 (WP + WPO ’
OPOTHB= 0.5*(OPDTH + DPOTHO)
OFJOTR= 2.5*(O0F JDTH+ OF JOTO
NACOTS= C.5*(BACDTH+ DGCDTO)
DGROTR= £.5*(DQ20TH+ NQROTO)
SMOOTR= 0.5*(SMODTH+ SMDDTO
SMVOTB= 9.5*(SMVOTH+ SMVOTO)
SVBOTR= 0.5*(SVRDTH+ SVBDTO
SOBDTB= 5.5*(SDANTH+t SOBOTO)
SLOOTB= 0.5*(SLODTH+ SLODTO
SFSO9T8= SvB80Ta+ SOBDTS

WEITE(6,101) THET&I,PyVyTsU,PHEEy
NP OTH38, WPB,OFJOTB,SMOOTS,SMVOTBy
OQCDTB, OQROTS»SVEDTB,SDBOTB,ySLOOTB,
SUMNQL »SUMFB,SFBDTS,WKNET,RR,y
UCHAR,NVDTHsNELEL2,ITERPC,
XMNOEsYNOE,XNOE,SOOTE
INEWP= INEWP+ 1
IF (€ CINEWP/7) *7.EQ.INEWP)
IF (IOTAIL.EQ.0) GO TO 77
_€& DETAIL OUTPUT.
WRITE (6,162).
XXXXX=  XMFA+ XMPAA
WPITE (65105) XMA »TAyUNsPHEAy3TAA, XMNOAY
1 XMF A
IF¢THETAT.LT. THEIJ) GO TO 77
WRITE (6,104) XMFG,TFG »UFG,PHEFG y»STAFG
WRITE (6,193) XMFL TFL ,UFL
WRITE (6,103) XMFLByTFLB,UFLB
IF (II.LT. 1) 509 19 77
90 78 T=1,IT
WPITE (6,166)
4 XNONE Q(T) sXMSOOT (I) sXMFLI (IT), XMPAM(T) »XMFSM(T)
78 CONTINUE

vr une

oAND. CIOTAIL.EQ.0))

60

WRITE (6598)

XNOANEy XMSOTASXXXXXoXMPARY

Ty XMBM(I) yTSMCT) »UIMC(T) »PHE (CT) sBTASM(T) XMOF (ID,


FMECK FINUN
SUB2CUTINE FINUNOTHETAT, OF UOTH, XMF J)
© TRIS RCUTINE GIVES CONSTANT FUEL INJECTION RATE

c ASSUMING CONSTANT INJECTICN 2ATE
COMMCN/SCAS TKEMAX THE TS eTHEIGsDINJoTRENS g TREEVOsTHEVCSOTINGG
1 QTHMIN, OTHMEDyOTHMAXsCTHETA,OT

OF JDTH=90.0
TF CCTRETAT LY .THETU) CRe(THETAT GESTTRETU+OINS)0) RETURN

DFUATHS YMEU/OTNG
SETURN
ENO

69


NOMENCLATURE |

Specific heat

Gas binary diffusivity

Pressure

Gas constant

Temperature

Specific internal energy

Error function

Equivalence ratio based on fuel-air mass ratio
Crank angle

Density


a= Tt
RHO= XMAO/Y
CALL TPROP(TA,> VXA,CVA,CPAsUA,H By AMWA,GAMA,RA,PHEAD
CALL SFECY(P,TA,PHEA,&MWA,XH24,XHZOA yXCON,XCOZAyXNZAyXNOAXOHA,
4 YHA, ¥Y OAs XNO,XO20,XCHGA)
YO24= XO24*32.74MWA
YNOA= XNOA*(33¢/74MWA)
XMNOA= XMA*YNOA
XMOANE= XNOA
XMNOAG= XMNOO
YSA= 6.6
XMSOTA= 0.6
XMSTAOQ= XMSOTA
SUMM= XMAD
EO= XMagFUA
DTHETA= OTHMAX
C PRINT INITIAL CONDITIONS
WRITE (6,INITUD
WSITE (6,98)
SET ZERO
GFG=0.6

in START THE INTESSATION
2 CONTINUE
XMAd= XMO

49


a

fone)

Qo

oO

ro)

READ (5,0PARAM)
WRITE (6,0 PARAM)

MODIFICATION OF EMPIRICAL PAPAMETERS THROUGH NAMELIST INPUT
READ (S,EPAR AMD
WRITE (6,5 PARAM)

MODIFICATION OF ENSINE PARAMETERS AND FUELING PROPERTIES THROUGH INPUT
READ (5,GTVEN) .
WRITE (6,GIVEN
OTROT= FPM* 3.14159/30.

OTHO «0 *RPM
ATOTH=1.9/9THOT

MEAN SWIRL AND MEAN PISTON VELOCITIES
UST= 6.25*0E*OMEGA
UPBAR= RFWFREPM/15,

FUEL COMPOSITION CiW)H(XX)
THeOFETICAL (0),(%) COMPOSITION FOR STOICHIOMETRIC BURNING OC(YSTOA) N(Z5TO4:
YSTOA= 2. 0* A459. 5*XX
ZSTOA= 3.764*YSTOA
FUEL SPECIFICATION °*
MOLECULAR WEIGHTS OF FRESH AIR,FUEL AND THEORETICAL COMBUSTION P®ODUCTS
AMNWAT= 22 7
MHFS WEL2 60 4+XX
XL298= OHFV
REGSRU/OM WF
CALCULATE STOICHIOMETRIC F/A RATIO REFERENCE TO PURE AIR
FAROSSAMWFES (CW 456 25FXX) R03 2.0430754*28.6))
OXST= 0.252/FAR0S
DENSITY GF LIQUID FUEL (GM/CM3)
IF (FOEN.LE.9.0)
LFOEN= 141 57 (OOTH131 25)71.6939%32.6179/05.96243* 1.9E-3
FUEL INJECTION VELICITY :
ININJ= OINJ+ 6.5
DINJ= IDINJ
THEEIU= THETU+I ING
DTINJ= DINJ*OTOTH
XJFS XMFI/OTINI/NNOZZ
ANOZZ= 3.25*P1*ON07Z *0N0Z2
VFJO= XJF/(F DOZEN *ANOZZ*COCHG)

NC. OF 3URNING 20! TO BE TREATES
NXXX=(DINJ+0.00 14D THMIN) /OTHMIN
NSS= NSHELFNXXX
XMFJOS= 09,95*XMFJ
XMFJO9= 0,99¥*XMFJ
IF (HCOM.NE. 3.0) GO TO 39

APT FORMULA FOR HIGH HEATING VALUE OF FUEL
HHV= (19659-44640 .*#(G>T-10.0)/1.3F AMWE

HEAT OF COM3USTION OR LaHeV.
HCOM= HHV 70.5*XX*DHVH20
aetna CONVERT ALL THE HEATS TO CAL/GM
HRV SHH V/A MWE
HCOM=HCOM /AMWF
Za CONTINUE

HEAT OF FORMATION (SAL/GM) OF FUEL AT 298 K
IF (HFFL.NE.d.9) GO TO 4G
HEELS HOOMFAMWNE + WEHFCO2+ 0.5*XX *HFH29G
HFFL= HFFL/AMWF

“4O CONTINUE

I.e. GF LIQUIN FUEL AT 8.°.%

UFFL= HFEL= ATOXJ/FOEN
UFL3P= CPFL*(T39F-298.9) 4UFFL

47


6C TO it
2 CONTINUE
C BURNING RATE
OMLOTI= 12.5663 7062FRHCO*XMNCT
TF (INO2.EG.0) GO TO 11
C MASS FRACTION CF N2 ON SURFACE
WXXG= WH 0625*XX
XXX1= GGLO*(HFEWXXG*XCO2(TIZXC2CT))
XXX2= L2690F(G.S5FXKX#WXYGFXHZOCII/SXC2CT))
XXX3= 2B. O*CWXXG*XN2ACTISXC2CIID
XXX=  (OXXXL*XXX2IVAXXXG
YMN2L= (1.0-YFLIZ(120+XXX)
XXXX= XMOOT/AL
IF (XXXX.GT.200.) XXxXxX=Z00.
ERL= FXPC-XXXX)
ERLF= EPL- ERF
E°L8= ERL- ERB
ESFR= ERF=- ERB
TFMTLE= (TF=TL)/ERLF
YFLE= YFL/ERLE
TEMTBE= (TF -T8)SESFA
YO2R YO28/ERFR
YN2L8E= CYN2L=YN28)/ERLB
C COMPUTATICN OF DISTRIBUTIONS AROUND A BURNING ORCFLET
DORT= (RFSRLI/NITST2
RI= RL-DRT
“99 3 IR=1,NITST
IF (IP.EQ.NITST39 ORI= (RA-2FISINITSTI-NITST2)
RI= RT4+D9T
ROFCIRV= RI
XXXX= XMDOTZRI
TF CXXXXeGTeZ600) XXXX= 200.
EQ= EXPI=XXXX)
IF (RI.GT.RF) GO TA &
TEMPESATURE DIST2TRUTION --- FUEL SIDE
TOF (IQ )= TE+ TEMTLE* CERF-ERD
YFOF(IRI= YFLE*(ER<EPF)
VO2ZOF(I°)= 060
GO TO 7
C TEMPERATURE DISTRISUTION --- CXIOANT SIDE
6 CONTINUE
TNECIRI= TS+ TEMTSE* (ER-EPAR)
YO2DFCIP)= YC28E*(ERF-ER)
YFOFC(IRI= 0.0
CN? OTSTRIGFUTION
7 CONTINUE
YN2DF (TRID= YN28+ YN2LSE*(ER-ERA)
8 CONTINUE
C CALCULATE HETERGENOUUCTION PATE APOUNT A FURNING CROPLET:
IF (OMLOTT.LT.1.0&-15) GO TO 14
CALL HETNOR(P,YNO8sAMWEM(T),RHOOsNITST,CNCRyENOOTID
NeCOT2(TI= CNCOTZ(1)+ CNCOTI®XNOI(I,L)
C MASS INVCLVEQ IN HETERGGENEOUS FCRMATION
XMDFTI= 0.0
RO 14 ITR=NITST3,NTITST
RT= O.SF(RDF(T2+ ROFCIR=1))
8l2s FL¥et
NRI= RDOFC( TRY + QOFCI9-1)
TOEAVE O.S*UTOE(IS)+ TOF CIR-1))
XMOSTS XMOFI+ RE2FORT/INFAV

°

72


Table 4.

Pollutant Computation Control

FORTRAN FORMAT | PRESET
SYMBOL DESCRIPTION S
IN@1 Homogeneous NO rate calculation is to be I l
performed if ING@1# 0
IN@2 Heterogeneous NO rate calculation is to be I 1
performed if INO2 # 0
ISSZ@T | Soot kinetic calculation is to be performed if I L
= IS@OT #0
Table S. Numerical Stability Control
NTER Max. No. of iterations allowed for stable solution] I 50
NAXIL No. of axial fuel jet stations (max. = 20) I 20
NSHEL No. of conical fuel jet shells (max. =4) I 4
NSZG No. of droplet size groups(max. =3) I 3
NITST No. of spatial points in droplet field intergration IT 50
(max. = 50)
DTHMIN Intergration step size (C.A,) used during fuel F 1.0
injection and combustion
DTHMED | Integration step size (C.A.) used after combustion| F 2.5
DTHMAX Integration step size (C.A.) used before injection | F $.0
THEMAX Maximum crank angle to be integrated to F THEVC
EPSIT Relative error bound on averaged temperature F 2.0 E-3
calculation :
EPSIP Relative error bound on pressure calculation F 2.0 E-3

13



© GROSS CYLINNES PROPERTIES
& AVERAGES DENSITY AND INTERNAL ENERGY
RHO= SUNM/V
U=SUMU/SUMCKM
c AvEeza GAS CONSTANTS IN GM3-4TM/G4-K
RR =(SUMR/SUMNM) #RORU
5 AVERAGED PRESSURE
B= SUMRT/V#RORU
AVEPAGED TEMPERATURE
T= PEV/(SUMMFRRD
QVER=OLL POLLUTANTS
XMNOE= XMNOA
SOOTE= XMSOTO+ XMSOTU
IF (II.LE.0) GO TO 6G
00 39 I=1,II
SOOTE= SOOTE+ XMSOOT (I)
XMNOE= XMNOE+# XMNO{I)
59 CONTINUE
65 CONTIN
IF tITERP
> accuracy CHE
CHKT= TTX
IF (ABS(CHKT) «ST. EPSIT*TO) GO TO 57
CHKP= P=PX
IF (QES(CHKP) .GT. EPSIP#PO) GU TO 57
CHKXNO= X¥MNOE= X4NOEX
IF (XMNOE.GT+S.3E-4 «AND. AIS(CHKNO) 6GT.EPSINJ*XMNOE) GOTO 57
CHKSOT= SOOTE- 300
IF (SOOTE.GT.1.9E-5 .AND. ASS(CHKSOT).GT.EPSIST*SOOTE) GO TO 57
se TO Se
£7 CONTINUF :
& CALCULATION O9£S NOT PASS THE ACCURACY CHECK
IF (ITEP°C.LT.NTER) GO TO 64
NOITE (6,129) THETAL,P,TA,TEG sy DTANTH, TFGUTH,OPITH
WEITE (6,119) ITERPC, TAO +TA,XMNOAO +XMNOS
WOTTE (6,119) ITERPC,TFGO +TFG
IF (II.LE.9) GO TO 56
DG 45 I=1,IT
WEITE (6,119) I, TBMO(I) + TBM(I) yXMNOO (1) +¥MNO(T)
1,OT30TO(I) ,DTADTH(L) »ANODTI(I) ONODT2 UT)
£© CONTTNUE
&6 CONTINUE
STOP
64 CONTINUE
PX P
Tx= T
XMNOEX= YMNOE
SOOTEX= SOOTE
ITE2PC= ITEPPCE 1
GO TO @4
54 CONTINUES
SUMFB= XMFP8+ XMFO8
IF (SUMFE.GT.XM£J99) KOUNT= NSS#4
SUMMK=SUMNK+ &3S(HP) *DTHETA
SUMIL= SUMOL+ COQ@CITH+ DOCITOF OQRDTH OGFOTO) FOTHET2
WKNETS= WKNET+# (AP +90) FOTHET2
C DETAILED OUTPUT STEP SIZE CONTOOL
IF (COPDTH*OPOTHO) «LT. 0.0) GO TO 10
OTINTH=(T=TO)/N THETA
IF (COTIOTHFOTITHO) «LT. Ged) GO TO 19

59

a

a

£9.90) GO TO 64



FOECK CCNYT
SUBRCUTINE CCNHTCRET,TI,OCCITH,SA,COS)

C THIS POUTINE FSTIMATES OVER=OLL CONVECTIVE HEAT TRANSFER TO AN UNIFORM

© TEMPERTUREO CYLINDER WALL BASED CN MASS AVERAGED GAS TEMPERATURE
COMMON/GONVEZS XJsATDOTA A XK,SIG sR sPUSCTOTHsDTHOT,OTROT
COMMCN/ENGEOS DE, RFWe RL y VO gSCeSeUST CR yg RE MyTWyONGZZyNNOZZeDBGWL
4 WF JO

s+Sc

Oh FXK*(RETFFCRI/CE

2.7 PS 39816 35* OE DE

AP=AH

ASLVI= 2.1415¢*0F*XLI

ST= ARF AP+ASLVI

OQCG= KWFAI*(TI-TW )

SOCITH=ANGFOTLTH

RETURN

END

FOECK FPF
FUNSTION FSF (7)

ER2QP FUNCTICN EVALUATION
X=?

=1.0

FUX.LT.0.04=-1.6

x

(X#*2)

TH1.0/01.4+.32759114X)

EMIX25EXPOXT)

ESFSA* (1. EMTX2*T#(.25482959Z+4T * (= .28UGSE7IEFTR(L.G2LUL37G14T* (<1,

1453452027+T#1.0514605429)))9)

RETURN

aD)

a

S

Ke

65


(SGIVEN
LO SEFAFAM

SCPARAM and Case

Ist Case ILD. Card

SGIVEN

Ist Case

Preset Engine and Fuel Paramete:s

Figure 2. Input Data

15


TABLE OF CONTENTS

ACKNOWLEDGEMENTS

NOMENCLATURE

I.
II.
Ill.

(1)
(2)
VI.
VI.
VIII.

INTRODUCTION

PROGRAM DESCRIPTIONS

DESCRIPTION OF THE PROGRAM INPUT
DESCRIPTION OF THE PROGRAM OUTPUT
SAMPLE CASE

INPUT

OUTPUT

PROGRAM STATISTICS

REFERENCES

APPENDIX -- FORTRAN LISTING OF DIESEL ENGINE
PROGRAM

ili

Page No.
i
Vv
1
3
9

17
al
22
23
39
41
43


Ze

TMETAUNES.CADS 97,09 — PEATHD= S.7UGECT WECHTD= 2.095F402 TEKD= 1.620E0C3 UMMALJGHD® W.SC3E-91 PHFES 5.958E- a1
RPVMEATHPCLI2—2.253F 090 HPCCALSCAD® L.Od3ECCL OF IOTHIGH/CADS 0, SMODIM(GH/CA)® 2,.693E-G2 SHVOTHIGH/CA)= 2.093E-03.
DISNTHICALZCAIS 342526069 DOROTHICAL/CAD= 6.767E-01 SVBOTHIGM/CAI= 0. SOBOTHIGH/CAD= 7.084E~04 SLODTHIGH/CAI= Ge
SUNTLICALDS 1,32 G2 SUMFS(GHD= 9,163E-02 SFBDTHIGK/EAD= 7.08GL-04 WKNETICALI=-2.146Et02 RRICHI-ATH/GH-KD= 2.59KE 006
UCHLOICM/SECI= GS4AEOG3 OVOTHIGHS/CAD = 7.349006 DELEIZ= 2.6170-03 ITERPCE

WMMOE(GHI= 5.054E-33 YNOES 2.3G3E-03 KNOZ(PPH)= 2.175E4C3 SOOTECGMI= 3.275604

THETACIEG.C
GROTH ETYsCa

19.50 PCATMD> S.270E¢01 VICHNI= 2.251£002 TAK) 2,6056003 UICAL/GH)2-1.972E401 PHEES 5,903z-G1
2e2GSEVOT WOCCMLZCAD= L.CKIECOL OFJOTHIGH/CAD= 0. SHDOTH(GH/C4)= Z2:540E-02 SMVDTH(SH/CAI= 2.560E-63

GUENTHICAL/CAD= 3.098F090 ONPOTHICAL/CS)= A.2A7E-G1 SVROTHIGM/CAI= O. SONDIHIGH/CAD= 5,497E-04 SLODTHIGN/CAD= v,
SUICLICALDS 1.601f692 SUMFUCGHI= 9.280602 SFONTHIGH/CAD= 5.497E-04 WKNETICALI=-1.936£002 RR(CMI-ATM/GH-K)= 2.692600
UCHAP(CM/SECI= 4.55Z2E603 OVOTHIRMS/CAD= A.2U1E600 OELELZ= 2.329%-03 ITERPC= 1

XMNCEUGHD= 5.922803 YNOE= 2.317E-93 XNOECPPMI= 2.1886 953 SOOTEIGMD= 3.268£-04

THETAUVEG.CADE 21.00 — PEATMD= GARUTEFOL VECMSD= 2.424002 TEKI= 1,597E403 ULCAL/GHDS-2.185E+01 PHEES 5. 90.

at

OPITHCATM/TAD=-2.106E409 WRI AL/CADS 1.656001 OFJOTHIGH/CAI= 6. SHODTH(GH/CA)= 2.3882-C2 SMVDTHIGH/CA)= 2.356E-53
COOGTHICGL CANS 2.944F 400 ICKOTHICAL/CA)= 5.822E-G1 SVBCTHIGM/CAI= O. SOBOTHISM/CAI= &291E-04 SLOOTHICH/CAd= Ge
SUISLICELI= 1.472E652 SUMFRIGMD= 9, 3720-02 SFUDTHIGH/CAD= 4.291£-0% WKNETICALD=-1.7266002 RRICHI-ATH/GH-K)= 2.8918 006
UGHEPICM/EECIS HES7ECIS OVOTHICMS/CAD= 9,047E400 OfLes2= 2.036f-03 TTERPCE 1

XMNOTGHI= 5.966E-33 YNOL= 2.3258-93 XNOECPPMD= 2.1976042 SOOTECGHD= 3.200€-04

THETACNEG.CADS 23669 " PEATHD= &.HU9E091 VICM3D= 2.612602 TEKI= 1,569E003 UMCAL/GMD=~-3.271F 401 PHEE= 5.963t-51
OPITMEATSZCR)E-2,97CE OGD WPRECAL/CAD= 2.0587 401 OFJOTHIGH/CAD= 0. SMDOTHEGM/CA) = 2.239E-02 SHVOTHIGH/C3)= 2.2350-33
OnCoratcaLscey= GhECII ANROTHICAL/OAD= 5, 3916-01 SVMOTHIGHZCAD= 0. SOMOTHIGH/CAI= 3.429E=C4 SLODIHIGN/CEI=
L.SUTERO? SUNFVEGMI= G.445E-02 SFOOTIIGHSTADS 3.429-94 WKNETICALD=<16524E402 RR(CMS-ATH/G4-K I= 2.890E00G
TISH/SECDS G.562E4G3 OVDTMICHI/CAD= 4.AO7EHOG OLLEL2Z= 2,095(-93 TTERPCE 1

AMIGE GMD = 5,955E-53 YNOF= 2.329E-33  ANDECPPMD= 2.252703 SOOTEIGHD= 3.253F-04

THETACMEG.CAY= «25,00 — PEATMD= &GGZEOU1 VECMID= 2.N1GEF02 TIK)= 1,549E003 UICAL/GHD=-4.368E601 PHZE= 5. 903E-O1
MEATH(ATY/CAI==2.834E9I5 WPICAL/CA)= L,045E 901 DFINTHIGH/CAD= 0. SMDDIM(GH/CAD= 2.097E-02 SMVOTHIGN/CED= 2,097E-63
OTCITHICAL/TOD= 2eH5CF OIC DOLDTHICAL/CAD= 4.9BNE-OL SVAOTHIGH/CAI= 2. : SOMOTHIGH/CA= 2,879E-0% SLODIHIGH/CAD= Ce
SUNCLECSLDS 1.636462 SUMFS(GHI= 9,554E-02 SFROTHIGH/CAD= 2,879. -04 WKNETICALD=~1,305€002 RRICMS-ATM/GN-K)= 2.837E005
USHORICM/SECD= G5ERF O53 OVOTHICMS/CAD= 1,0580001 DELTA2= 1,5560-03 THERPC= 1
ANCE GMD = 5.961E-03 = 2e33LE-O3 XNOELPPMD= 2.205403 SOOTLIGNI= 3.246E-Ch

THETACDEG. CADE — 3G4Ce — PUATMD= 3,252C4U1 VICM3D= 3.391002 TIKI= 1.49KE003 ULCAL/GND=-6.94BE001 PHEE= 5.9632-01
OPATHEATN/CAD==2.521E450 WRICAL/EAD= 9.961E 400 DFIDINIGM/CAD= SHODIHIGH/CAD= 1.837E-G2 SNVDTHIGM/CAI= 1.837E-03
OGCCIAICALSCAIs 2o41KE 6G DOPOTHICAL/CAD= 4.1A7E-01 SVBOTHIGH/CAI= 0. SOSOTHIGH/CAD= 2.466€-0% SLOOTH(GH/CAD= Cy
QLICALDS 1,751E462 SUMFICGHD= 9,633E-02 SFANTHIGH/CAD= 2.460£-04 WKNETICALI=-7.983E¢01 RRICH3-ATH/GH-K}= 2.8S4E 006
UCHAEICH/SECDE G57HFOCS OVNTHICHS/CAD= A.262E4dL OFLCI2= 5.082¢-04 ITERPCE 1

MANGE (GHD S6G5HE-0% YNDT= 2.S3GE-D3 ANDECPAMD= 2,2U7L093 SOOTEILMD= 5.2298-04

sree INOTCATED ENGINE PERFOPMACE AUD ENGINE EMISSION -----

TMeCePoCPSIIS 7.4836 001
T.HePoCHPD= — 2.025E631
LiS.FeCe(LB/HP=HQI= 4 692KE-01
SOOTIMG/M3STPD= — -1.454E002
NITFIC OXIDE (PPMD= —2.207E +03



#DECK SOOTR7 :
SUBROUTINE SOOTRZCTAV,V)

G THIS RCUTINE ESTIMATE THE SOCT FCEMATION ANC DESTRUCTION RATE
COMMCNZINTEGRY INO1,INC2,TSOCT,IWETBT yNPRTsIQTAIL sNTERgNAXILy

t NSHEL »NSZGNITST,IPUGSII,IIT
COMMCN/SPECYIZ = XH20100) +XH2C(100) 4XCC (100) ,xCC2(100) »XN2(100),
1 XNO(160).XOH(100),XH(100),x0 (100) »XN(1009,xX02(100)
2 » XCH4 (10-0) »XNCNEG (100)
COMMCN/PRATET/ ONODT11109) ,DNCOT2(100) ,OMSOT1(100),OMSOT2Z(100),
t DRMDTHT190) sOLEDTH(100),XMDF (100)
COMMCN/TPROPT/ P, TBM (100) ,URM (100) »HOAMC100),VEM(1009,CPBM(100),
1 AMWEM (100),83M (100) sPHE C100) sAFR(100)4XMFEMC LOO)»
2 XMPAM (100)
COMMON /AIRP/ XMAIyXMF Jy XMEGRyXMRDgXMAO,XMARD yXMAEGRyXMFRO sg XMFEGRy
1 AMWA,AFRA,FARA,PHEA
COMMCN/FUELP/Wy XX a VFL y ZFL 9 VSTCAs ZSTOSD yg API» FLEN  FARCS pAMWF y SUSFT
t »XMUF »CEN,OXSTyHCOM

COMMON/ORCPC/ HUTIM,DTPHY,XTLYO2,CPBAR,RECDsTCRITF yXL29By
XMBM(100),C0VEDTH(100),CLOOTHI{100) »XMFLIC100),
DI(10052) sOMFLI (10053) pOMLOT (10053) »XNDT (10043),
XKK (10053)
COMMON /SOOTPY SDEN,DSODT, ASCTP,ASOTAsXMSCCT (100)
oo 1 I=1,1T
OMSOTL(TI= 0.0
OMSOT2(TI= 0.9
IF (ISOOT.EQ.9) GO TO 1
T8= TRM(T)
yO28 = XO2(12*32./4MWBM(T)
AFOVAL= XMPAM(TY/(XMFBM(TI4XMFLICID)
PHES= 1.0/(AFOVAL*FAROS)
CG SOOT FORMATION RATE
BMSOTLOTI= ASCTA*(VBMCI) FVBMCTD/V) (TAV/T8)* (PHEB**3) FPF
1 ExP(-20000./T8)
© SCOT OXTOCIZATION PATE
OMSOT2Z(T)= 6.51E4*(XMSCOTI(ID*#YC28*P IV (SOEN*OSOOT*SCRT(TADY*
1 EXP (-19650./T8)
1 CONTINUE
RETUSN
EMO

Wwe

82


eeeee AVF=24., THEIJ= 29 DEG.BTIC, MED. SHIRL (CASE O60 ) deaee

25


ACKNOWLEDGEMENTS

The CAPE-20-71 Steering Committee of the Coordinating Research Coun-
cil and the EPA originally defined the project and guided the work to it's final

outcome,

The ultimate value of this project is largely due to the Committee,

without whose help and criticism it would have been difficult to complete.
The authors sincerely acknowledge their gratitude to the following individuals
who have served on the Steering Committee during the lifetime of the project:

F. J. Hills, Committee Chairman
J. L. Bascunana
N. G. Beck

T. C. Belian

J. E. Bennethum
G. L. Borman
W. L. Brown

J. C. Hoelzer
G. D. Kittredge
G. A. Lavoie
D. F. Merrion
P.C. Meurer
J. M. Perez

S. M. Shahed
A. D. Tuteja

A. V. Wilson
A. E. Zengel

Mobil Research & Development Corp.
Environmental Protection Agency
International Harvester

Coordinating Research Council
General Motors Research Laboratory
University of Wisconsin

Caterpillar Tractor

International Harvester

Environmental Protection Agency

Ford Motor Company

General Motors-Detroit Diesel Allison
International Harvester

Caterpillar Tractor

Cummins Engine

General Motors-Detroit Diesel Allison
Cummins Engine

Coordinating Research Council


CMSCVX*EY
CPAM(I)= CPX¥FL
GAM= CPRM(T)/CV

C PARTIAL CEPIVATIVES OF U AND R WeReTs Pals AND PHE

10

O1Pe = PO2(2)+3.0G*C0213) FPHES?

D2PF = (002(5)+%.9*ND2 (6) FPHE2) FTL

CIPE =((002(8)+3.0*002(9)*PHE2) + CO2C11)* TI) FAL CGP
NEES F*(DIPFED2PF FOP)

PFIOFS -FAI*F1*F AROS

PUPPBMCT)= F1FN*(D30/P1)

ouprev(t) = Cv

PUPFEM(I)= FLIF(-P+DOF) + PFAPF TUX

PFPCRM(TI= FLFC#(PR2(5) *PHELFRR2(8)) SFL

PRETRM(TI= “FIFO (RAZ (4) FPHEL+QP 207) FL OCOSTE
PRXOF= PO2(2)+ CFERAZ(3)+ WQ2(G)FTI+ FR205) FALCGP)
POPFBM(I)= FA*PRXPFH PFLFF ROX

CONTINUE

RETURN

ENO

97


25 CONTINUE
DFVOTH= OFGOTH* SVBOTH
82 CONTINUES

ITERPC= 9

Po= P

To= T

TAO= TA

TFGO= TFG
TFLO= TFL
WPO= WP

HAQ= HA

HFGO= HFG
HFOSO= HFOS
OPOTH OPOTH
OTOTH OTIOTH
XmD00 XMD
XMVO XMV

XMFVO = XMFV
XMFOO= XMFA
XMAPSO= XMAPQ

XMPAA XMPAA
XMFLO XMFLOD
XMF PB XMFP3
XMFOB XMFOB

XMNOAO= XMNOA
XMSTAQ= XMSOTA
SUMFJQ= SUMFJ
IF (II.LE.9) GO TO 63
90 62 I=1,1T
TaMO(L)= TEM(T)
XMNOOCI)= XMNOCT)
XMSOTO(I) = XMSOOTIID
XMF3MO(T) = XMF3MCT)
XMPAMO(I) = XMPAM(T)
DO 32 L=1.NSZG
OLo(I,Ly= OCT.)
IF (OL(IT,L).LE.9.6) XNDI(T,LI= 0.6
82 CONTINUES
IF (1.GT.ITOLO) GO TO 62
DAPOTH(T)= 0.3
DVBOTH(I)= 0.9
OLOITH(T)= 0.6
62 CONTINUE
63 CONTINUE
C IF NO NEW SYSTEM IS CREATED ,SPARE THE FIRST ITERATION
IF (¢ THETALT.GT.THEVC «AND. THETAI.LT«THEIJ) 20%,
1 {IGEND.GT.3)) GO TO 49
€41 SONTINUE
SUMi= 9.0
ded
0.9
0.9
SUM5= 0.6
OTLOTH= 0.9
SOBOTH= 9.9
SMOOT 1]
SMVOT

$2


ee

THETACDEG.CADs 18.06 — PEATMI= 2,798E001 VICH3I= 2.1716602 TEKIs R.GRKEFOZ ULCAL/GHD® 7.298E01 PHEES 8. 3116-92

CPUTNCATH/CADE 162256600 WPICAL/EAD==5.300E 900 UFJOTHIGH/CA)= &.8A2E-03 SHUDTHIGH/CAD= 0. SHVDTHIGH/CAD= 0.
DCLUTHIGAL/CAD= 2.2E3FOG0 OORDIMICGL/CO)= S.354E-06 SVUOTHIGH/CAD= 36 SDBOTHIGH/CAI= Ca SLODTH(GW/CAd= Gy
SUMLLICALD2 5.356009 SUMFA(GHI= Ge SFBOTHIGH/CAI= 0, WKNETICAL)©92.458£402 RRICMS-ATH/GH-K)= 2.840600
UCHA>ICHATECHS 2.35KE0GS OVOTHICMI/CAD S@7 8020400 OELCIZ= S.6n7b-03 LTERPCS 1

XHNOL (GH:

Se921E-26 NOES 2-404E-24 XNOECPPHD= 2.312E-18 SOUTEIGNDS 3.2640-15

THETAINEGS

DE 227,00 — PEATMD = 2.926692 VINMID= 26G95F902 TIKI= B.755E002 WICAL/GHDS 7.3BGEOO1 PHERS 16173£-01

OPISHCATH/CED= 16255E400 WPKCAL/CAD=-5.256E 400 OFIDTH(GH/CAD= $.84ZE-03 SNOOTHIGH/CADE 0. SMVOTHICH/CAI= O4
OCCOTHMICOL/CAD= Le T19E9GO OMANTHICALZCAD= S.SA1E-04 SVADTN(UM/CAI= Os SOBOTHIGN/CAD= 0. SLOUTHIGH/CAD= Go
SUvOL(DALD= 7.26GF+u0 SUMFAIGHD= 0. SFONTHIGH/CAD= 0. WKNET(CALD

2e5LOE*O2 RRICHS“ATH/GH-K)= 2.840E 090
UCHAUICH/SECD= 2,562E093 OVOTHIGMS/CA) =-7.369b400 OLLCIZ= 3.679£-03 ITERPC= 1.

KANDECGMD= 1.676E-23 YHOE= G.3PHE-24 KNOE(PPMI= &.202E-18 SOOTE(GAI= 420766915

THITAUTEG.CADS 16.00 PEATMD= $.0525001 VICHSD= 2.0236+02 TEKI= B,825F002 UCAL/GMD= 7.467EC01 PHEED 1.515£-61
CPOTHCATHZCAD= 1424184060 WOLCAL/CAD=-5.1892 400 OFJOTHIGH/GAI= 5.882E-03 SMODTHIGH/CAD= 0. SMVOTHIGH/CA) = Oy
OICOTHICAL/CAI= 14357E430 OOROTHICAL/CAD= 5.808E-04 SVAUTHIGH/CAI= 0. SODDTHIGH/CAD= 0. SLOOTH(GH/CAD= 0.
SUVOL(CELI= B.622E660 SUMFRIGMD= 0. SFROTH(GH/SAI= 0, WENET(CALD=-2,562E¢02 RRICMS-ATH/GH-KI= 2,2396400
UCHAP(SH/SECI= 24742E463 OVOTHICM3/CA) =-6.972E400 O&LE1Z= 3.765F-03 ITERPC= 1

MANGE (GMD 1.6991E+22) YNOFZ O.c77F-24 XNOE(PPMD= 7.7AZE-18 SOUTEIGMI= 5.046E~-15

THETS(OEG.08 715.00 PeATH

312A OL VECMSD= L.956F002 FCKI= B.B9LEFO2 ULCAL/GHD= 7.544401 PHEES 1.857E-01

MOMTHCATM/CADS 2,391E600 WRICAL/CAIS-5.G9RE 460 DFJOTHIGH/CAI> 5.AA2E-03 SHODTHIGH/CA)= 0. SHVOTHIGH/CAl= 0.
DICOTMICAL/CAD= 1,466£ 600 OOROTHICAL/CA)= 6.G52E-04 SVHOTHIGH/CAI= 0. SOBOTHIGH/CAI= 0. SLOOTHIGH/CAD= 0.
SUMQLICALDS 3.903E001 SUMFACGMI= 0. SFANTHIGH/CAD= Cy WKNETIGALD=-2.613E002 RRICHI“ATM/GN-KI= 2.837E400

UCULEISH/SECH= 2.910€4C3 OVOTHIC43/CAD =-6.553E006 OiLEAZ= 3.721k-03 LTERPCS 1
MAHCE(G4D= 3.669E°23  YHOE= 1,4A7E-23 XNOE(PPHI> 1.436£-17 SOOTEIGMI= 6.222E-15

THETAQIEG GANS +1460)  PUATHDS 3.323801 VACHSS LeROPENGZ TCKI= B.953E402 UIGAL/GM)= Toot3EeGL PHEES 2. 20GE-d1
QHITHCATM/CED= 143106400 WRC ALSCADH-6,941£ 0G OFIOTHICM/CAD= S.ABZE-03 SHDDIM(GH/CAD= 2. SMVOTHIGM/C4)= 0.
OQCOTHICALJCAD= 1,455E06G MOROTHICAL/C4I= 6.254F-04 SVADTHIGH/CAI= 0. SOBOTHIGH/CA)= 0. SLOOTH(GM/CAD= Cy
SUNCLICELDS 1.468C451 SUMFAIGH= 0. SFANTHIGH/CAD= 0, WKHED(CAL}=-2,663E402 RRICM3-ATH/GN-KD= 2.534E090
UDMES(CM/SECDE 3,G5EE003 OVOINICMS/CR) =-6.930F 000 OLLEL2Z= 3,738f-03 TTERPC 1

FANOE(GHD= G.643E-23 YHOES 2.6A9E-23 XNOECPPHI= 2.604£-17 SOOTCIGH= 7.6116615

THITSEXEG. CADE -13,60  PEATMDS 3,70GC UL VECMSDS LeHS3C0O2 TEKI= 9.d62E002 ULCAL/GND= 7.598E401 PHEE= 2. 5422-01
OPATHCATH/COP= L,754E400 WPKCAL/CAI=-5.070E 460 OFSOTHIGH/CAD= 5.682€-03 SMOGTH(GH/CA) = 3.016E-03 SMVOTHIGH/CAI= 3.518E-34
DINDIN(CAL/CAI= 1.503F 690 OUKMTHICAL/CAD= &.ABSE-O2 SVROTHIGH/CAD= 5.759E-03 SOADTHIGH/CAI= &.162E-04 SLOOIH(GH/CAD= 5.CGEE-J3
SUNOQLICOLD= 2.3ILEOCL SIMFACGMDS T.A77E-3  SFBOTHIGH/C ADS 7.177E-03 HKNETICALD2-2. 7146002 RRISM3~ATM/GN-KI= 2.dU2E400
UCHERICHASECI= 3.219€663 OVOTHICNI/CA sTGwC+OO OCLEL2= 1,698C-03 ITERPC= 2

KMNCE(GMDS 2.055E-05 YNOE= 8.3936-C6 XNOE(PPHD= B.GLar+GO SOOTEIGHY= 1.204E-C6

THETBISEGCAD= 912,00 PLATMD= G.227E001 VIGMSD= Le779EHO2 TIKI= 1.00KE403 UCCAL/GMD= 7.6Z3E401 PHEE= 2.8942-O1
OPOTHEAT“/CEIE 2.25ZE+00 WRICAL/DAD=-5 511646) OFIOTHIGH/CAD= 5.382E-03 SMODTHIGH/CAD= 5.913E-03 SHVOTHIGH/CAI= 5.913L-o6
DCLDIMILAL/CAY2 1.905E%I0 OORDTHICAL/CAD= 1.4595E-01 SVBOTHIGM/CAD= 5.820£-03 SDBDIHIGH/CA)= B.430E~0% SLDITHIGM/CAD= 5.945E-63
SUMGLICALD® 16498E0G1 GUEFOIGHD= 26394E-02 SFANTHIGH/CAD= 6,663E-03 HKNET(COL)=-2. 7676002 RRICMS-ATH/GM-K}= 2.85CE005
USHOZICH/SECD= F.361E9TS OVOTHICNI/CA) 2-S.270E000 OELELZ= 7.9272-04 ITERPC=

KMNGE(GHDS 7.734E-05 YNOE= 3.1202-05 XNOECPPMD= 2.9972+01 SOOTE(GHI= 5.116E-06



20

oO

©

oO

oO

oO

oO

20

UNI

NOMELIST/CPARAM/ INI1,INO2,1SGOT,

IBUG

1
“2 os NTHMIN, OTHMEDsOTHM4 Xy THEMAX
3 sEPSTT,EPSIP,EPSINO,EPSIST

NAMELTST/EP A447 C4,C8,CNOR,

Sad

MAMELTST /INITL/

i
2
CONSTANTS

R=32.654

RU
THERMAL CONDUCTIVITY GAL/CM-SEC-K

1.987

XK= 3.4155E-4

S-3 CONSTANT IN CAL /SEC-CM2-K4

STG= 1.55 38F-12

AT OF FORMATION (CALVMOLE) OF CO2,

aHUTIM,OTPHY, SEFU
2 »RHOD yCPEAR,XTLYO?
THE TAI, THEMAX 4? ,
FAROS,XMEGRyXMAE

HFC )2=
HFH2
OHVH20=
HFH20G=
CONVERSION
PI= 3.1
XJ= 4.1
a7Tn=
OT2= 3.
RORU= R
ATOXJ=
INPUT EQUIL
CALL TA
SPECIES TAB
PRANGL=
PRANG2=
TRANGL=
TRANG 2=
PHE®Gi=
PHE2G2=

1.5

a4 .ISER
$.3 63
LU.5ES
HFH20+
FACTORS
41592654
BS74E7
32556
14159/130.
SRS

&TO/XJ

I3kIUM COMPOSITION TALES
PEIN
LE RANGES

PTa (41)

PTRCNPT 2)

TT3C4)

TTSCNTT 3)

PHET2ACL)

PHETB(NPHT 3)

QHVH20

CVINTsoMIXs SVEN, DSOOTSASOTAyASOTP
»COCHG,»CEVAP,CBUSN,AIGND

T9VyUAyHAyXMAD,AMWA
GR XAFEGP oe XMROe XMARD» XMFZODVAMWE,
OXSTyXO28s5OEN, HHV, HCOM,HFFL HFS,» VCeVOIS, VF JO

UNIVERSAL GAS CONSTANTS IN CMZ<ATM/G-MOLE-K

S4L GAS CONSTANTS IN CAL/G-MOLE-K

H20 AT 298 K

FROM T4PE3

PPE-SET ENGINE PARAMETERS AND FUELING PROPERTIES

HFFL= 9
XMRO= 0
FOE D)
HOOM= 0
READ (5
WRITE (6

1 CONTINU
PeE-SET CON

INTEGERS

INOt= 4
INO
Isaor

6

ae

+c

of
+GIVEN)
»GIVEN)
iE

TROL PARAMETERS

NPRT= 2

IOTaIe
IZATE
T2°9PP

SPlYP=
NTE
NSZ

rome

ww

45

sNTESyNPRTsMNAXILNSHEL
yNSZGyIDROPP,IRATEP, ISPCYPyIWETBT yLOTAILyNITST

oFORQ,PHEA,


*DECK TPROP?Z
SUBROUTINE TPEOPZ
Ceeteeeeese THIS RCUTINE COLEULATES --------
C#Feeereeee FCUTLIARP TUM THERMOSYNAMICS OF CCMRUSTICN PROOUCTS
Creer eeeee UIP,ZT,SHE) AND PR(P,T,PRE) AND THEIR CERIVATIVES
GOMMCN/INTEGR/S INC1,TNC2,TSOCT,IWETBT »NPEFT,ICTAIL »NTERyNAXIL,

1 NSHEL»NSZG,NITST,ISUC,IT,III
COMMCM/TPROPIZ PoTOM (100) ,U3M( 100) sHEME100) pVOM(1COd,CPEMI1000,
1 AMWBM (100) ,R3"(100) sPHE C100) ,AFR(109),X¥MFBM(100),4
2 XMPAM (100)
CAMMON/FUFLP/WyXXsYFL4ZELYSTOAsZSTOA SAFI yFEEN pFAROS s AMWF y SURFT
1 aXMUF sCENSOXSTsHCOM

COMMCH/PANGSS TRANG1sTRANG2Z,FRANG!I yPRANGZ yPEERGL,PHERSG2
COMMON/FOPVIZ/S PUPPBM(100) sFUPTBM(100),PS5PFEM (100) ,PUPFEM(100),

1 POPTIM (100) ,PRFPBM(100

COMMCN /HECOE/ AA1(6), E22 (E) »DD1 (11) ,82108), 60206) ,89216) ,0N2011)
1 +RR2C8)

Pi= Pp

TF (P1L.ST.PRONG2) PL= SRANG2

Oc 10 I=t,It

TiL= TAMC(T)

CHEL= PHE(T)

TE (T1.GT.TRONG2) Ti= T2ANG2

TS (PRE1.GT.PHESG2) PHEL= PHERG2
TISIMLHOW/ Tt

T2s Ti*T4
T?= T2*Tt
Tus T3*T4
TS= T4FTL
PHE2 = PHELSCHEY
PHE3= PHEL*PRE2

ALOGP= SLOGTF1)
IF (OKE1.CT.1.9) GO TO 2
41 CONTINUE
© FUEL LEAN
GQ= AALCADFAAL(2VFTLEOALC IDET? EAALCGINATS FAAL (SIFTS #AA1 (6) TS
BALCAFABLC2|.FTAFABLCSIFT2?] 4 9BLCGVFTS +881(5) 474 4281 (6) FTS
ot COLCL#COL (2) FPHE1+ NO1(3) PRES
N2= CNNLCG) FONTS) FOHEL+ NOL CE) FPREZDFTI
NSO=C(NOL(7IFON1L (CBVFPHEL+ OOLIGD*FRES) + (0011004001011) FPHEL) #TT)
N3= D3CFALOGP
¥XX=01402+03
IF, (XXX.LT.-200.9 XXX=-209.
D2 €XP (XXX)
UX= Ae @FPHELt 9
XXX CREL C34 BLOG) FTI+ PRICS)FALOGFIFPHEL# 2R1(6)+ RRAC7PFTI+
LREL (8) FALOGP
IF XXXL Te-2006) XXX=-2006
Or EXP (XXX)
RX= ORL(L)# RELC22F°OHF 1+ C
UX4PXFTL
ALLC2VF2.9FAAL(SVFTLF J.C AALCYIFT2+ G.CHOAL(SDFTI+ S.OFAAL(E)
1 *T4
FELC]V +2 .0F8R1(SVFTL4 B.OFSRLCHIFT2+ GOFEBL(SDFTS 45.07 82116)
Th
Des -€N01(4)+ OCLOSI*PHEL+ COLCEI*PRES) *i0CC./T2
O3P= -(791(10) 8 COLC11) FEHEL)IF1000./T2* ALOGE
ne f¥(G2P+C Py
CV¥s AS-28%—PHE14DP
COXSC UK +OY

w

95


SGIVEN
OE
REN
Re

cr

sc

TH
RPM
sco
ONOZZ
NNOZZ
DeOWL
THEVG
THEVO
xXMAT
EGR
XMRO
OMEGA
NSTROK
THEIJS
ONS
XMF J
PINS
TFS

w

xx
FOEN
CEN
API
xMUF
SUPFT
TePF

“

o

a

ry

"

0.1397E +02,
0.762€901)
Co 30GAEHO2s
Q.17E402y
0.2032E+005
0.SE003%
O.15E404,
Q.1€+01,
0.254E-01,
6
0.82E+01y
-0.15€+03,
G-14E +03,
G.2HE+01s
ety

6.0,
0.6€403,

&
=0.26402,
G.17E402y
0.164005
0.231403,
0.35E+03,
O.7E+0Ly
0.126402
0.0,
0.u67E+02,
0.3458402,
0.12E-C1y
0.278402,
0.37766 +03,

28


CONVF2T UNIT TC PER MASS OF MIXTURE
XBASE= 1.0+FA80S#PHE1
Fiz 1.0/XBASE
UBM(T)= UX*F4
HAM(T)= HX*FL
RAM(T)= 2X9F4
AMWBM(T)= 1.987/RBM(T)
Vx= G41.6°RBM(T)*TIsP4
CVSCVXFF 4
GPEM(I)= CPX*F1
GAM= CPBM(T)/CV
C PARTIAL CERIVATIVES OF U AND R WeRsT. PyTy AND PRE
DIPF = 001(2)+3.0*D01(3) #PKE2
D2PF = (0N1(5)+3.07001 (6) *PHE2) #TI
D3PF =C€(N01(8)+3.04N01(9) *PHE2) + COLC14)* TI) FALOGP
OPF= C*(D1PF+C2PFHDSPF)
PFIPF= -F1*F1*F AROS
PUPPAMITI= F1°D*(D30/P1)
PupTa¥(1) = CV
PUPFAM(T)= FLIF(-B#DPF)+ PFLPFFUX
FRPPAM(T)= FLFC#(2R1(5) FPHELERR1 (BDI SPL
FRPTEM(T)= -FL®C* (R91 (4) #PHEL*RA1(72994000./T2
PRXOF= PRA(2)+ CF(REL(3)+ RALCHIFTIF RR1(5) FALOGP)
PRPFRM(I)= F18PRXPF+ PFLPFERX
co To 10
2 MCNTINUE
© FUEL RICH
A= ANZC1)FAA2C2IFTAFANZ(IIFTZ HAAZ(HIFTE FAAZ(SVETH #AAZ(G) FS
B= AB2S(1)F+BE2(2)FTLEBAZ(3IFT2 +BE2(4) FTI +BEZ(SVFTH 498216) *TS
D 02(194+002(2) *PHEL+ CO2(3) *PHES
(0D2(6) #00215) *PHEL® DN2(6) #PHES) #TI
D30=((0N2{7) +002(8)*PHE1+ DOZ(9) PHES)+ (002 (10) 400211) *PHEI)#TI)
D%= O30*ALOGP
XXX=D14#024D3
TF (XXX LT. +2006) XXX=-200.
Ds EXP(XXX)
Ux= A= 8*PHEL4 D
XAX= (ROZCID+ RAVCGI*TI+ RRZ(5)*ALOGPI*PHEL+ RRZ(EV+ REZ(TIETIFE
1RR2(8)* ALOGP
IF (XXX6LT 22006) XXX==200.
C= EXP IXXX)
QX= RR2(1)+ &RZ(2)*PHEL+ C

H YeRKFTL

APS AAZ(2)+2.0*AAZ(3IFT1+ Z.OFAAZ(HIFTZ# GOFAAZ(S)*T3+ S.OFAA2(6)
1 FT

BP= A2Z(2742.0FBB2(3VFT1+ F.G*AB2(GI*T2+ G.O*GB2(5)*TS 4+5.0%082(63
1 *Th ;

D2P= -(DN2(4)+ OD2(5)*PHE1+ OD2(6)*PHES)F10C0./T2
DIP= -€002(10)+ O02(11)*PHELIF1000.7T2* ALCEGP
np = OF (N2P+03P)
CVX= OP-RP#PKE1+0P
CPX=CVX#ERX
CONVERT UNIT TO PER MASS OF MIXTUPE
YBASE= 1.0+FAROS*PHEL
Fi= 1,0/X84SE
UeM(T)= UXFF1
HEMOT HXFF 4
FPM(I)= OX*FL
AMWAM(TY= 1.987/REM(T)
VX= GLEFRBM(T) FTAP 1

96


Table 1. Engine Geometry

ceipsien DESCRIPTION FORMAT | PRESET
DE Engine bore (cm) F * NA**
RFW Crank radius (cm) F NA
RL Rod length (cm) F NA
CR Compression ratio F NA
sc Engine clearance (cm) F NA
TW Mean engine wall temperature (Cx) F NA
SCR Supercharge (atm) F NA
DNZZZ Diameter of fuel nozzle (cm) F NA
NNG@ZZ No. of fuel nozzles F NA
THEVC Clesing angle of intake valve (C.A.) F NA
THEVG Opening angle of exhaust valve (C. A.) F NA
DBOWL Diameter of bowl on piston head (cm) F NA
NSTROK No. of strokes per power cycle [eee 4
Table 2. Engine Operation
XMAI Mass of fresh air/powerstroke (gm) . F NA
XMRD Mass of residue gas/powerstroke (gm) . F Option*
OMEGA Charge swirl angular velocity (rad. /sec) F NA
THE! Fuel injection time (C.A.) F ; NA
DIN Duration of fuel injection (C.A.) F NA
XMFI Mass of fuel injection/powerstroke (gm) P NA
PINJ Fuel line pressure (atm) F NA
TF] Fuel line temperature (°k) F NA
EGR % of E.G.R, (base on total charge) F NA
RPM Engine speed (rpm) F NA

+ Floating point
* This value will automatically be calculated in program if not specified.
** Not applicable

wee Integer

ll


11 CONTINUE
WEITE (6,159) IJ
6a To ?
12 CONTINUE
InMis 18-4
opT= PTIM(IN,J) = PTIM(TOMi,4)
APTA=C(ETIM(IO,J)-TIMI) /NPT
OPT2=(TIMT<STIM(L OMS yd) OPT
XSS(ToJ) = XSCTOM1)) «-FDPTL1 + XSCIO) FOPT2
YSS(T,J)= -YSCIOML,JI#09T1 # YS(I0,0)*0PT2
AYSS(T,JD= AYSCIOML,J)FDPT1 + OYS(TC,J) *OPT2
REOFT= FPHOFCIOM1, J) FOPTLF2HCF (CIG,I) FOPT2
AFTI= AFT(ION1,J) *DPT1it AFT(ID,J) *OPT2
VSSIJ= VSCTOML, JPFNPTL + VS(TO,u)FOPT2
WSSIv= VSSTUFVEIO
XMFVS(I4J)= 6.0
IF (TIMI.LT. HUTIM) GC TG 17
S VAPORIZATION FUNCTION
00 3 L=1,NS7G
IF (9(IyJd,l) L020) GO TO 18
X¥XL= 1.04 CROAR*(TA-TLIZXL
XKKO= 2.0*PT*¥RHODFALOCG(XXX4)
CSS= XSS(TsJ)
IF (2SS.GT.PRONL) &SS= 2ECHL
USIT= RSS*CMEGA
QEN= SHCAFUSTTFEN (Ty Jol )7XMUM
CREV= 1.0+ CEVAP*SGST (RED) #PRI
XKKK= CPEVFXKKO
OC2NT= -4.0*XKKK/OT/SFOEN
Otis OT
IF CCTIMT-9T).LTSHUTIMD PTL= TIMI- HUTIM
OP= MCTysoLd*2 + DOZOT OTL
TF (M2.LT.9.9) D2=0.0
O(T,J,L)= Ssert(p2) ~
18 AXXX= GOS7ELIE(LI-( DCT, SL) 252302 (LD) FF3)
XMFVS(I,J)= XMEVS(Tyu)# OXXX
3 SONTINUE
SAVEAFS AFTIZKMFVS{(I yu)
AFV(I,J)= SAVEAF
TAFV(TyJd= (CTLEAFVIT,J)#TADZ (1. GFAP VL, Jd)
G SCOT FOPMATION
AFVJ= AFV(I4J)
TAFVJ= TAFV(T,J)
VIJ= 2.0°T*YSS (I,J) #OYSS (I,J) FVSSTU*ET
IOTGJ= 6
vows VIU/VC
CALL SCOTPIP,TAFVIsSAFVUs YC2,VCyT ,LOICIyXMSCT (Ty J) VIS OMSDT,
1 OVSOT1,O0MSOT2)
YMSOTI{T J) = XMSOT(Iys)+ CMSOT#OT
© SELF-IGNITION CRITERIA (AFTES SHIPINSKIyETsAL«)
PXX= 0.57 (8+PO(T, Jd)
TXX= G.5FT #TO(Tyd))
IG= OTGNO = *EXP(UEGS./TXXI/ (CENGO* (PXX**O.396))
IF (OIG+LE.TIMI) 50 TO 7
17 CONTINUE
C UNPURDEN FUEL VAPAR
XMEGL= XMPG44¢ XMPVSCT SJ)" COMFIUT
XMAFAL= YMOF O14 AFTIFOMFIJT
XMSOT1= XMSCTL+ XMSOT(I,J)
«7 79 2

67


FOFCK SPATY
SUBPOUTINE SPECY( PX, TX sPHEXs AMM gXH2 4 XH2G o XCC¥XCO2 9 XN29XNO2 XCH
ty XHyXOy¥NyxXS2yXCH4)
C THIS SCUTINE CALCULATES TEE EGUILIARIUM SFECIES CIST°I@UTION CF & GIVEN
© STATE T.€-6 (PyT,PHE) FOR A SPECIFIC FUEL FRCM TRIPLE TABLE INTERPCLATION
COMMCR/SSTAR/S NPTBNTTG,NPHTByOTE(7),TTBC13) sPHETS(11) 4
XH2TB07413411) 9 ¥HZOTBU7 413411) sXCOTRI7 13911) gy
XOO2ZTIC7 13,11) eXN2TRI7,1 7,11) pXNOTBU7e13,11)¢
XOHTOC7 13,11) oXHTBC7e13511) XCTEC7 13,11) 5
XNTPC7 423,110 ,XO2TEC7g 12944) »XCH4TB(74134511)
COMMCN/YXCON/ IDP,ICT,IOPHE,
TDP1L,IOTL,ICPH1, SFL PsDELT,DELPRE,9P1,0T1,0PHE1,
OP2,0T2,C0PKE2

AM

Nie

OHES BHOX
DO 1 IOP=2yNeT3
TF (CP.LGE.PTAC(IOP=-13 JANN. P.LE.PTSCICFIY GE TC 2
4 CONTINUES
WRITE(6,100) F
100 FORMATC*® ESRCRese. P OUT CF RANGE *, 1PE10.2)
IF (P.LT.PTI(1)) GO TC 7
f= PTECNPTA)
IOP= NPT?
Go TO 2
7? CONTINUE
Fs PT2(4)
Ipe= 2
2 CONTINUF
oO 3 [GT=2.NTTS
IF (T.GE.TTIR(IOT<$1) .AND. TLE.TTACICT)) GC TC 4
2 CONTINUE
WPITS(F, 194) T
LOL FOOMAT(® ERROR .eoee T SUT OF RANGE * 4 1PE1C.3)
IF (T.LT. TTE(1)) GO TO #
T= TTBCNTTA)Y
Tots NTTP
Go mT! &
CONTINUE
T= TTECL)
IOT= 2
G CONTINUES
90 5 IOPH 2, NPHTS
TF (PRE.GEsPHETA{IDPHE-1) eAKC. PHE.LE.PHETECIDOHE)) GO TO F
5 CONTINUE
WRITE(6,102) SHE
192 FORMATOC* FRROR .weees PHE CUT OF RANGE *,1PE1C.3)
IF (PFE.LT.PRET3¢1)) GO TC 9
PHE= PHETRCNPHTS)
IMPHE= NPHTE
GS 10 6
9 CONTINUE
PHE= PHET3(4)
TQOHE= 2
& CONTINUE
THL= TOP-1
TPTi= TOT-4
TPH. TOPHE@1
CFLP= PTA(ICE)= PTRC(ICP1)
OFLT= TTACINT)= TIRCTOTL)

oo

83


c

Cc

ro}

ao

6

BURNING SYS

IF (II 20) 59 73 79

S

MIXING PATE OF OT%-3URNING SYSTEM

13

CALL AMIX (XM44, SUM34,UCHAR BTABM AFR yOMONTHs OMVOTH y»CMIX,CVIOT
00 18 I=1,II

SMDOTH= SMDOTH+ DMODTH(I)

SMVDTH= SMVOTHt OMVOTH(T)

CONTINUE

ETOA= 1.0/(1.9+4F R48)

ALPA= 1.0- ETAA

XAAOTHS = ALPA*SMOOTH

XFAOTH= - ETQA *SMOOTH

FOUILIGRIUM CALCULATION

SPE

55
or
FST

Est

2

CALL TPPOPZ

CIES DISTRIBUTION
CALL SPECYZ
00 55 I=1,I1T
VeM(T)= v*BTaaM (1)

SONTINUE
FFUSICN SURNING RATE OF THE SYSTEM
IMATE HETEROGENE JUS NO PRODUCTION RATE

POLL FLAMEZCTEOS, TFL »CPFL,CNOR)
IM@TE HCMOGENEOUS NG PRODUCTION RATE
CALL HOMNOZ{(CNOR)

ESTIMATE SOOT FORMATION RATE

CALL SOOTFZIT,V)

ATs OF RADTATIVE H.T. FOR EACH BURNED SYSTEM
GALL RADHZZ

oc 71 [=1,11

SOBOTH= SOBOTH+ OL3OTHIT)

ETASM= L.O7(1.0+AF RTI)

ALPSM= 1.C-ETABM

XFIOTH (TD) = (OLSOTH(T)+ OVBITHIT)+ ETASFOMDDTH(I)- ETASM*OMVOT
LH(T))
XAITQTH(I) = (OLPA*DMODTH(T)- ALPBM*OMVOTH(I))

XFAQTH= XFANTH+ ETABM*OMVOTH(T)

XOAQTH= XAANTH+ ALPOM*FOMVDTH(T)

ONOITH(T)= (ONOITL(IY+ONONT2¢1)) *OTOTH + YNOA*OMODTH(T)
1 - YNOCI)*OMVOTH(T)

NET PRODUCTION RATE

ONSOT= OMSOTI(TI- OMSOT2(T)

OMSITH(T) = OMSOT#OTOTH+ YSA*OMDOTH(I)- YS(I) #OMVOTH(T)
OGROTH=0Q20TH+O WMOTHC(T)

DCMOTH (I) = GBM(I) #0QCOTH

Q2T= CUFLBSUFL) FOLOOTRI(IT)+ (HF OS-UFLB)*OLODTH(I)

OEM(I)= OSMOTH( T+ OCMOTH(T) + Q2T

OPHOTH(T) = (XMPAM(T) #XFIOTH(CI)-XMFBM(T)FXAIDTH(III/S(FAROS#XMPAMCET)
1 **2)

XX4=  ( H4-H3MCI)) FOMOOTH(T)+ (HEOS-HBM(T)) *OLBOTH(T)

1 + CHFG -HOM(I)) #OVENTHUT)

X¥O= XMAM (LT) *#(R 3M (1) #TOM(I)/7P-PUPPBM (I) -TSM(T) FPRPPBMII))

X¥B= XMOM (TY *(PUPFBEM(I) +TBM(T) #PRPFBM(I))

XXG= XM3M(T) #( OUP TIM (I) +TSMC1) FPRPTBM(I) +R 3M (7)

BUMCT) = CXK4L-XX3* NPHOTHOT) -98M(L)) /XXC

AAM(T) = XXA/XX3

SUM1= SUM1+ PBYCI)FTOM(T) * (OMDOTH(T) 40L 80TH(I) tOVBOTH(T) -OMVOTH(T)
1)
SUM2= SUM24+ XM3MCT)* T3407) PROF SM(T) *OPHOTH(T)

SUM 3+ 88M (TL) FXMEM(T) #(RSMCID +TBM(L) FPRPTBM(I))
SUMG+ XM3M(TIFTOMCT)#9RPOQM(TD

53



TF (Il «LE.9) GO TO 38
oC $7 I=1,1I
XX4= XX4+ CHIM CT) -HA) FOMVOTH(T)
SNODTH= SNOOTH+ YNO(T)*#OMVOTH(T)
SMSOTH= SMSOTH+ YSCI)FOMVOTH(I)
27 SONTINUE
38 CONTINUE
PHAOTH= (XMPAS*#XFADTH= XMFA*XAADTH)/ (FAROS*#XMPAA*XMPAD)
AAAS XX1/XX2
BBA= (-XX3*PHADTHH+XX4=Q4) /XX2
SUML= SUM1+ RA*T&*(SMVOTH=SMODTH)
SUM SUM2+ XMA*FTAFPRPFFPHAOTH
SUM SUM3+ 334 *XMA* (RA+TAFPRPT)
SUMG= SUM44+ XMAFT AFP RPP
SUMS= SUNS+ AANFXMAF(RAFTO*PRPT)
BNOOTA= 0.6
IF (INO1.EQ9.9) GO TO 26
CALL HOMNOR (P,TA,XNOAyXN24 4XO24, XNOANE »CNOR,XNODTAI
YNOOTA= XNOOTA*(30./AMWA)
ONODTA= XMAFYNODTA
CONTINUE
ANOITH= ANOOTA*OTOTH+ SNOOTH= YNO&*SMDOTH
YC24= XO2Z8*32./4MHE
OMSOTA=0.0
IF (XMAM.LE.d.9 602. ISOOT.EQ.3) GO TO 32
VAA= XMIAFVASXMA
IDIGI= 3
CALL SOOTR(P,TA ,AFRAVYG2A,Vy Ts TIIGTs XMSOTA,VAN,SOMSOTA,AMSOTA
1 sAMSOT2)
32 CONTINUE
AMSITH= OMSOTA*FOTOTH+ SMSNTH- YSA*SMOOTH
5 PRESSURE GERTVATIVE
WR= P*OVOTHFATOXS
OPOTHS (-AP+SU414+SUM2+SUM3)/ (VFA TOXS=-SUMG=SUMS)
TF COPDTH.GT.2G4) ARITE(E,109) ITERPC WP yVyATIXI,
1 SUM1,SUM2, SUM3,S5UM45 SUMS
TEMPERATURE DERIVATIVES
OTAOTH= ASA*DSOTH+ 5BA
TFGOTH= AAFG#DPOTHt BEFS
IF (IL «LE.2) GO TO 6
DO 7 T=1,IT
OTSOTH(I) = AAMCI) *OPOTH+ 8B4(T)
7 CONTINUE
CONTINUES
IF (ITERPC.NE.9) GO TO 50
~ PREDICTORS. 66 PREDICTORS
49 CONTINUE
OTADTO= OTADNTH
ANOOTG= ANODTH
AMSOTO= AMSOTH
TFGOTO= TFGOTH
IF (II.LE.0) GO TO 54
00 52 I=1,IT
QBMO(T)= O8M(T)

0
oa

oO

OMONTC(T) = OMIOTH(T)
OMVOTO (I) = DMVOTH(T)
OL3NTO(I)Y= OLBITHCT)
BTRITO(I)= OT8ITHCT)
ANOOTO (I) = ONOOTHC(T)
IMSITC (I) = OMSOTH(T)

$5


FOECK WETAT
SUIRQUTINE WETST(PsTSsAMWAGAMWE yAg RoC yTWyXL)
5 ESTIMATE STAPLE WET-BULB TEMPERATURE (N-FEPTENE)D
COMMON/PROSCS HUTIM,OTPHY,XTLYC2sCPEAR,RHOD,TICRITF yXL298y
XMBM (100) ,OVENTH(100) ,CLOOTHC1CG) ,XMFLI(1C0),
910100, 3),OMFLI(1005 3) ,CMLOT(100,5) »XNCI (100 3)»
XKK (10543)
OTW= O.MS*(TCRITF-298,)
TWw= 208.- OTW
O80 1 ITF8=1,59
Th= TH+ ATW
IF (TH.GE.TCRITF) GO TO 3
XLS XL29URFCOTLOTTF<THI/(TCRITFE-298.)) ¥¥0.38
CPRAQ*(TS-TNIZXL
Xx h- 2/7 (THK-C)
PFL= 10.**XXX
XXX= PFELFAMNE
XXXSOXXN+ (P=PEL) FAMWA)D
YFL/(1.C-YFL)
C (XPHS=XLESI/XRHS
T© (ASS (CHK).LT.O.01) GO TO 2
IF (ITEREQ.1) CHKI= CHK
IF (CRKK*CHK1 LT. C20) SC TO 2
Ge Tl 4
3 CONTIBUE
Tye TWwe OTW
OTWs §.570TH
1 CONTINUE
WITTE CE, 102) TTER Py TB oT Wy XQESyXLESSCTW
100 FORMATOSX,*eeeee ERFOR IN WETAT ITERATION. oo *¥ oT 3e1P7E10.3)
ST9e
? CONTINUE
RFTUGK
END

ane

99


Table 6.

Output Control

rSAtne DESCRIPTION FORMAT | PRESET
NPRT Print every NPRT intergration steps I 2
IDTAIL Print detailed zonal information if IDTAIL #0 I 0
IDR@PP Print detailed droplet information if IDROP # 0 I
& IDTAIL # 0
ISECYP Print detailed species information of each zone if I 0
ISPECY #-0 & IDTAIL #0
ERATEP Print detailed rate information if IRATEP #0 & I 0
IDTAIL #4 0
Table 7
eyenee DESCRIPTION FORMAT PRESET
CA Preconstant for rate for heat convection rate F 0.26 -
CB Reynolds No.exponent heat convection rate F 0.75
CNOR Multiplier of N@ rate constant F 5.37
CVIDI Vitiation/dilution ratio F 0.1
CMIX Air/product dilution coefficient F 7.8E-3
SEF] Coefficient for swirl effect on jet entrainment F ; 4.3
CEVAP Crossflow effect on fuel vaporization rate F 0.276
RHOD PD in droplet calculation ¢m/cm-sec) F 3.0 E-4
CPBAR c, in droplet calculation(cal/gm) F 0.35.
XTLY@2 Droplet extinction limit (oxygen mass fraction) F 0.02
HUTIM Droplet heat-up time (sec) F 2.0 E-4
SDEN Density of soot particle (am/em >) F 2.0
DSOOT Assumed mean soot particle size (cm) F 1.0 E-4 +
ASOTA Pre-constant for soot formation rate in burning F 3.52
zone -
ASOTP Pre-constant for soot formation rate in F 3.52 1
vaporizing zone

14