Program EDDYBL

D .11 . ADDITIONAL TECHNICAL INFORMATION 423D.11 Additional Technical InformationThe program uses the conventional Levy-Lees variables [see Hayes andProbstein (1959)] and much of the program notation follows that of Harrisand Blanchard (1982) . The numerical procedure is the Blottner (1974) variablegrid method augmented with an algorithm devised by Wilcox (1981b)to permit large streamwise steps . Section 7.3 of the main text provides anin-depth discussion of the algorithm . This section first presents the governingequations for mean-flow properties and all turbulence-model equationsimplemented in the program. Then, the transformed, nondimensional formof the equations is presented for the k-w and k-c models .D .11 .1Mean-Flow EquationsThe equations governing conservation of mass, momentum and mean energyfor all models are the same . For compressible two-dimensional (j = 0) andaxisymmetric (j = 1) boundary layers, the program uses body-orientedcoordinates (s, n), where s is arc length and n is distance normal to thesurface . The equations are as follows .TS (Pu) + i an= 0 (D .3)_8u _8u _ _dP 1 _a _auPlc as +pv an ds + ri an ~ (P an + PT (D .4)Pu Ts +pv an =u ds+'U (aaiiza 'j 'U 'UTn) + Pc + r an I ( PrL + PIT) On(D .5)The perfect gas law is used as the equation of state and the fluid is assumedcalorically perfect so thatP= and h = CPT (D .6)In Equations (D .3) through (D .6) : u and v are streamwise and normalmass-averaged velocity components ; P, P and h are fluid density, pressureand enthalpy ; p and PT are molecular and eddy viscosity ; r is specificReynolds shear stress ; e is turbulence dissipation rate ; PrL and PIT arelaminar and turbulent Prandtl numbers ; T is mass-averaged temperature ;R is the perfect gas constant ; and Cp is specific heat at constant pressure .