TALK=T;RUN( 1, 1)
 
    GROUP 1. Run title and other preliminaries
TEXT(S-A 1D DEVELOPED PIPE FLOW
TITLE
 
  DISPLAY
  The case considered is the 1d solution of fully-developed
  turbulent flow in a circular pipe at a Re=1.E5. The solution
  is performed by use of the single-slab solver with a specified
  mass flow rate. Calculations are made with the high-Re
  form of the Spalart-Allmaras model with wall functions.

  This model predicts the following results compared to experimental
  data:
 
          f        vdp
  SA      0.018    3.96
  Data    0.018    3.75
 
  where f =8.*(wstar/wbulk)**2 is the friction factor, and
  vdp=(wcl-wbulk)/u* is the velocity defect parameter.
	
  In-Form is used to print the f, vdp and Nu to the PHIDA, RESULT
  and inforout files.
 
  ENDDIS
 
 
REAL(DELT,US,DIAM,WIN,REY,TKEIN,EPSIN,MIXL,FRIC,DPDZ,MASIN,DTF)
REAL(YPLUS_W)
DIAM=0.1;WIN=1.0; REY=1.E5;FRIC=1./(1.82*LOG10(REY)-1.64)**2
DPDZ=FRIC*RHO1*WIN*WIN/(2.*DIAM);US=WIN*(FRIC/8.)**0.5
TKEIN=0.25*WIN*WIN*FRIC;MIXL=0.09*0.5*DIAM
EPSIN=TKEIN**1.5/MIXL*0.1643
 
    GROUP 3. X-direction grid specification
CARTES=F;XULAST=0.1
    GROUP 4. Y-direction grid specification
YPLUS_W=40.0
 
ENUL=WIN*DIAM/REY;DELT=2.*YPLUS_W*ENUL/US;
NREGY=2; REGEXT(Y,0.5*DIAM)
IREGY=1;GRDPWR(Y,29,0.5*DIAM-DELT,0.8);IREGY=2;GRDPWR(Y,1,DELT,1.0)
 
    GROUP 5. Z-direction grid specification
ZWLAST=0.1*DIAM
DTF=15.*ZWLAST/WIN
 
    GROUP 7. Variables stored, solved & named
SOLVE(W1,ENTI);STORE(ENUT,LEN1);SOLUTN(W1,P,P,P,P,P,N)
STORE(STRS,YPLS,SKIN)
TURMOD(SPALART-ALLMARAS)
 
KELIN=1
 
 
 
REY;FRIC;US;DPDZ
 
    GROUP 11. Initialization of variable or porosity fields
FIINIT(W1)=WIN
FIINIT(ENTI)=1.4025E-4
 
 
    GROUP 13. Boundary conditions and special sources
WALLCO=GRND2
WALLCO
WALL(WALLN,NORTH,1,1,NY,NY,1,NZ,1,1)
 
 
MASIN=RHO1*WIN*DIAM*DIAM/80.
 
IENUTA = 21
DISWAL
 
  ** activate single-slab solver
FDSOLV(FLOW,MASIN)
 
    GROUP 15. Termination of sweeps
SELREF=T;RESFAC=1.E-5
LSWEEP=30;LITHYD=10
 
    GROUP 16. Termination of iterations
    GROUP 17. Under-relaxation devices
VARMIN(W1)=1.E-10;WALPRN=T
RELAX(W1,FALSDT,DTF)
RELAX(ENTI,FALSDT,0.5*DTF)
SPEDAT(SET,GXMONI,CLASSIC,L,T)
 
    GROUP 22. Spot-value print-out
ISWPRF = 1 ;ISWPRL = 100000
ITABL=3;IYMON=2;NPLT=2;NZPRIN=1;NYPRIN=1
IYPRF=1;TSTSWP=-1
 
   ** compute flow output parameters
 
(stored of FRIC is 8.*STRS/(:WIN:*:WIN:))
(stored of UTAU is STRS^0.5)
(stored of VDP is (W1[&1&]-WIN)/UTAU[&NY&])
(STORED of USTR is UTAU[,NY,])
(STORED of UPL is W1/USTR)
(STORED of YPL is ((YVLAST-YG)*USTR)/(1E-6))
 
(make ffac is 0.0)
(make ustar is 0.0)
(make vdus  is 0.0)
(store1 of ffac at walln is fric)
(store1 of ustar at walln is utau)
(store1 of vdus at walln is vdp)
   ** print to inforout file
(print of f  is ffac)
(print of u* is ustar)
(print of (ucl-ub)/u* is vdus)
 
DISTIL=T
EX(W1  )=   1.050E+00
EX(YPL )=   1.054E+03
EX(UPL )=   2.235E+01
EX(USTR)=   4.699E-02
EX(VDP )=   3.959E+00
EX(UTAU)=   1.576E-03
EX(FRIC)=   5.889E-04
EX(DWDY)=   1.304E+01
EX(LTLS)=   3.882E-04
EX(WDIS)=   1.786E-02
EX(VOR1)=   1.304E+01
EX(SKIN)=   1.646E-04
EX(YPLS)=   1.322E+00
EX(STRS)=   7.361E-05
EX(ENUT)=   1.493E-04
EX(ENTI)=   1.494E-04
 
 LIBREF = 0
STOP