Program EDDYBL

426 APPENDIX D . PROGRAM EDDYBLThe 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)Program EDDYBL 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 .