Program EDDYBL
Program EDDYBL
Program EDDYBL
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
426 APPENDIX D . PROGRAM <strong>EDDYBL</strong>The compressibility function F(Mt )is given by_ 0, Mt < MtoF(Mt) - { Mt - Mi , Mt > Mto(D .26)where Mt - 2k/a 2 , a is the speed of sound, and Mt . is given byMt, = 1/4 (D .27)D.11 .3 k-c Model EquationsFor the k-E model the equations for k and E are :_8k _8k _ _8u1 _8 aklpu pv as + an-pl pe an - + -r-3- an [Ij (p + PT/O'k) (D .28)an_aF _aE E _a ii E 2 1 a _aEpe as +pv an = f1CEikPTan-f2CE2P k +PE+r~ an [~ (4 -f PT/U,) an](D .29)whereand the eddy viscosity isE = E + Eo (D .30)pT = Ci f ;,Pk2/E (D .31)<strong>Program</strong> <strong>EDDYBL</strong> includes six low-Reynolds-number versions of thek-E model . The models differ in the form of the damping functions ft,fl, f2, eo, E, in the values of the closure coefficients, and in the surfaceboundary condition imposed on E. The damping functions depend uponone or more of the following three dimensionless parameters .+ k2 _ k i/2 n _ u .rnReT = Ry, y - (D .32)Ev,v vThe damping functions, closure coefficients and surface boundary conditionon E for the six models are as follows .