TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2018, Date 08/01/18 CPVNAM=VDI; SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings DISPLAY Library case 640: Tutorial: A duct flow (Immersol) ENDDIS ************************************************************ IRUNN = 1 ;LIBREF = 197 ************************************************************ Group 1. Run Title TEXT(A duct flow (Immersol) ) ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,35,1,47) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(138)=QRZ ;NAME(139)=QRX NAME(140)=PRPS ;NAME(141)=WDIS NAME(142)=WGAP ;NAME(143)=T3 NAME(144)=SCAT ;NAME(145)=EMIS NAME(146)=TEM1 ;NAME(147)=EPKE NAME(148)=DEN1 ;NAME(149)=EL1 NAME(150) =ENUT * Solved variables list SOLVE(P1,U1,W1,T3,TEM1) * Stored variables list STORE(ENUT,EL1,DEN1,EPKE,EMIS,SCAT,WGAP,WDIS) STORE(PRPS,QRX,QRZ) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) SOLUTN(TEM1,Y,Y,Y,N,N,Y) TURMOD(KECHEN) ************************************************************ Group 8. Terms & Devices TERMS (T3 ,N,N,Y,N,Y,N) NEWRH1 = T NEWENL = T ************************************************************ Group 9. Properties PRESS0 =1.01325E+05 ;TEMP0 =273. * Domain material index is 2 signifying: * Air using Ideal Gas Law, STP SETPRPS(1, 2) DISWAL DRH1DP = GRND5 DVO1DT =3.33E-03 PRNDTL(T3)=-0.0263 ;PRNDTL(TEM1)=-0.0263 PRT(KE)=0.75 ;PRT(EP)=1.15 EMISS =1. ;SCATT =0. ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(P1)=0. ;FIINIT(PRPS)=-1. FIINIT(WDIS)=0.05 ;FIINIT(T3)=20. FIINIT(SCAT)=0. ;FIINIT(EMIS)=1. FIINIT(TEM1)=20. No PATCHes used for this Group INIADD = F ************************************************************ Group 12. Convection and diffusion adjustments No PATCHes used for this Group ************************************************************ Group 13. Boundary & Special Sources No PATCHes used for this Group EGWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 500 RESFAC =1.0E-04 ************************************************************ Group 16. Terminate Iterations LITER(P1)=200 ;LITER(KE)=10 LITER(EP)=10 ;LITER(T3)=50 LITER(TEM1)=50 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,0.5 ) RELAX(U1 ,LINRLX,0.5 ) RELAX(W1 ,LINRLX,0.5 ) RELAX(KE ,LINRLX,0.5 ) RELAX(EP ,LINRLX,0.5 ) RELAX(LTLS,LINRLX,1. ) RELAX(T3 ,LINRLX,0.25 ) RELAX(TEM1,LINRLX,0.25 ) KELIN = 3 ************************************************************ Group 18. Limits VARMAX(T3)=3000. ;VARMIN(T3)=-204.75 VARMAX(TEM1)=3000. ;VARMIN(TEM1)=-204.75 ************************************************************ Group 19. EARTH Calls To GROUND Station GENK = T PARSOL = F CONWIZ = T IENUTA = 2 ISG62 = 1 SPEDAT(SET,FACETDAT,NUMOBJ,I,7) SPEDAT(SET,GXMONI,PLOTALL,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=0.1221 ;EX(U1)=1.799 EX(W1)=0.05239 ;EX(KE)=0.04272 EX(EP)=1.821 ;EX(LTLS)=9.06E-03 EX(QRZ)=9.113E+05 ;EX(QRX)=5.84E+04 EX(PRPS)=5.253 ;EX(WDIS)=0.06917 EX(WGAP)=0.3156 ;EX(T3)=1746. EX(SCAT)=0. ;EX(EMIS)=0.996 EX(TEM1)=1660. ;EX(EPKE)=13.09 EX(DEN1)=336.299988 ;EX(EL1)=1.521E-03 EX(ENUT)=1.541E-04 ************************************************************ Group 21. Print-out of Variables OUTPUT(WDIS,Y,N,N,N,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 18 ;IYMON = 1 ;IZMON = 25 NPRMON = 100000 NPRMNT = 1 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 ISWPRF = 1 ;ISWPRL = 100000 No PATCHes used for this Group ************************************************************ Group 24. Dumps For Restarts GVIEW(P,0.,-1.,0.) GVIEW(UP,0.,0.,1.) GVIEW(VDIS,0.72) GVIEW(CENTRE,0.5,0.5,0.25) > DOM, SIZE, 1.000000E+00, 1.000000E+00, 5.000000E-01 > DOM, MONIT, 5.000000E-01, 5.000000E-01, 2.677807E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > GRID, AUTO, T T T > GRID, RSET_X_1, -12, 1.100000E+00,G > GRID, RSET_X_2, -11, 1.100000E+00,G > GRID, RSET_X_3, -12, 1.100000E+00,G > GRID, RSET_Y_1, 1, 1.000000E+00 > GRID, RSET_Z_1, -13, 1.100000E+00,G > GRID, RSET_Z_2, -3, 1.100000E+00,G > GRID, RSET_Z_3, -15, 1.100000E+00,G > GRID, RSET_Z_4, -3, 1.100000E+00,G > GRID, RSET_Z_5, -13, 1.100000E+00,G > DOM, INI_AMB, YES > DOM, INI_BUOY, YES > OBJ, NAME, INLET > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, TO_END, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cube3t > OBJ, TYPE, INLET > OBJ, PRESSURE, 0. > OBJ, VELOCITY, 1. ,0. ,0. > OBJ, TEMPERATURE, 727. > OBJ, KE_IN, 1.000000E-03 > OBJ, EP_IN, 9.000000E-03 > OBJ, NAME, OUTLET > OBJ, POSITION, AT_END, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, TO_END, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cube12t > OBJ, TYPE, OUTLET > OBJ, PRESSURE, 0. > OBJ, TEMPERATURE, 727. > OBJ, COEFFICIENT, 1000. > OBJ, TURBULENCE, SAME , SAME > OBJ, NAME, WALL1 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, TO_END, TO_END, 0.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cube13 > OBJ, TYPE, PLATE > OBJ, SURF_TEMP, 0. ,727. > OBJ, EMISSIVITY, 0.9 > OBJ, NAME, WALLH > OBJ, POSITION, 0.000000E+00, 0.000000E+00, AT_END > OBJ, SIZE, TO_END, TO_END, 0.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cube13 > OBJ, TYPE, PLATE > OBJ, SURF_TEMP, 0. ,2227. > OBJ, EMISSIVITY, 0.9 > OBJ, NAME, PLT2 > OBJ, POSITION, 3.000000E-01, 0.000000E+00, 1.300000E-01 > OBJ, SIZE, 4.000000E-01, 1.000000E+00, 2.000000E-02 > OBJ, GEOMETRY, cube14 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 103,COPPER at 27 deg C > OBJ, EMISSIVITY, 0.9 > OBJ, NAME, B6 > OBJ, POSITION, 3.000000E-01, 0.000000E+00, 3.500000E-01 > OBJ, SIZE, 4.000000E-01, 1.000000E+00, 2.000000E-02 > OBJ, GEOMETRY, cube14 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 111,STEEL at 27 deg c (C = 1%) > OBJ, EMISSIVITY, 0.9 STOP