LOAD(N110)
TEXT(KO_2D FLOW THROUGH A ORIFICE PLATE :T302
TITLE
  DISPLAY
  This case is identical to library case N110, except that the
  calculation is performed with the k-omega turbulence model
  rather than the k-e model. The case considered is 2d turbulent
  axisymmetric incompressible flow through an orifice plate
  located in a pipe.
  ENDDIS
TURMOD(KOMODL)
WALLCO=GRND2
CONPOR(ORIF,0,VOLUME,1,NX,-#2,NY,-#5,-#5)
REAL(GOMEG)
DO JJ=1,NY
+ GYP=YFRAC(JJ)*YVLAST;GY=.5*(GYP+GYM);GYDR=GY/PR
+ GYDR2=GYDR*GYDR;GYDR4=GYDR2*GYDR2
+ GLM=0.14-0.08*GYDR2-0.06*GYDR4;GLM=GLM*PR
+ GYDR3=GYDR2*GYDR;GKI=1.+2.*GYDR/3.+10.*GYDR3/3.;GKI=GKI*US2
+ GEPI=0.1643*GKI**1.5/GLM;GOMEG=GEPI/(0.09*GKI)
+ COVAL(IN:JJ:,OMEG,ONLYMS,GOMEG)
+ GYM=GYP
ENDDO
PATCH(T1,NWALL,1,NX,NY,NY,1,NZ,1,1)
COVAL(T1,W1,LOGLAW,0.0)
COVAL(T1,KE,LOGLAW,LOGLAW);COVAL(T1,OMEG,LOGLAW,LOGLAW)
RELAX(KE,FALSDT,TSTEP); RELAX(OMEG,FALSDT,TSTEP);KELIN=0
FIINIT(OMEG)=FIINIT(EP)/(0.09*FIINIT(KE));STORE(EP)