Program EDDYBL

More documents

Recommendations

Info

422 APPENDIX D . PROGRAM <strong>EDDYBL</strong>ICPU = 2 : VAX FortranSpecial Comments : None .Compile and Link :for instllink instlfor eddybllink eddyblfor seteblfor edgefor graficfor initilfor ioeblfor main0for misclink setebl,edge,grafic,initil,ioebl,main0,miscICPU = 3 : SUN Fortran . . . SUN/OS or MS-DOS/SP-1Special Comments : Using the -03 option yields the maximum degree ofoptimization .Compile and Link :f77 instl .f -03 -o instlf77 eddybl.f -03 -o eddyblf77 setebl .f edge .f grafic .f initil .f ioebl .f main0.f misc.f -03 -o seteblICPU = 4 : Cray Fortran . . . UNICOSSpecial Comments : None .Compile and Link :cf77 -o instl instl .fcf77 -o eddybl eddybl .fcf77 -o setebl setebl .f edge .f grafic .f initil .f ioebl.f main0 .f misc .fICPU = 5: Silicon Graphics IrisSpecial Comments : None .Compile and Link :f77 -o instl instl .ff77 -o eddybl eddybl .ff77 -o setebl setebl .f edge.f grafic .f initil .f ioebl.f main0.f misc .f
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 .
Page 3 and 4: D.2. GETTING STARTED QUICKLY 393are
Page 5 and 6: D.2 . GETTING STARTED QUICKLY 395sh
Page 7 and 8: D.3. INSTALLING SETEBL 397D .3Insta
Page 9 and 10: D .5 . RUNNING A GENERAL CASE 399D
Page 11 and 12: D .5 .RUNNING A GENERAL CASEthe sta
Page 13 and 14: D .5 . RUNNING A GENERAL CASE 403ho
Page 15 and 16: D.5 . RUNNING A GENERAL CASE 405Tur
Page 17 and 18: D .5 . RUNNING A GENERAL CASE 407pl
Page 19 and 20: D .5 . RUNNING A GENERAL CASE 409di
Page 21 and 22: D .5 . RUNNING A GENERAL CASE 41 1l
Page 23 and 24: D.7. EDDYBL OUTPUT PARAMETERS 413Co
Page 25 and 26: D.8 . PROGRAM PLOTEB : PLOTTING UTI
Page 27 and 28: D.8 . PROGRAM PLOTEB : PLOTTING UTI
Page 29 and 30: D .9 . ADAPTING TO OTHER COMPILERS/
Page 31: D.10. COMPILE AND LINK COMMANDS 421
Page 35 and 36: D .11 . ADDITIONAL TECHNICAL INFORM
Page 41 and 42: D.11 . ADDITIONAL TECHNICAL INFORMA

Program EDDYBL

You also want an ePaper? Increase the reach of your titles

Delete template?

Save as template?