************************************************************ Q1 created by VDI menu, Version 3.5, Date 20/10/03 CPVNAM=VDI;SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings DISPLAY FLAME HOLDERS IN AN AFTERBURNER The case is aimed to demonstrate the 'cut-off' technique of PARSOL method to represent the inclined walls of combustor and flameholders on Cartesian grid. The demonstration case considers the turbulent two-stream diffusion flame in a 2D model of afterburner. Two flameholders placed in the flow are of the open-mouth-type V-gutters. The task is to calculate the temperature and composition of combustion gases along with all related field distributions. ENDDIS ************************************************************ IRUNN = 1 ;LIBREF = 290 ************************************************************ Group 1. Run Title TEXT(Flame holders in an afterburner ) ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,60,60,1,7.620000E-05) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T * Non-default variable names NAME(140) =KOND ; NAME(141) =SPH1 NAME(142) =DEN1 ; NAME(143) =FUEL NAME(144) =TMP1 ; NAME(145) =PROD NAME(146) =OXID ; NAME(147) =MIXF NAME(148) =EPKE ; NAME(149) =LEN1 NAME(150) =ENUT * Solved variables list SOLVE(P1 ,U1 ,V1 ,H1 ,FUEL,MIXF) * Stored variables list STORE(ENUT,LEN1,EPKE,OXID,PROD,TMP1,DEN1,SPH1) STORE(KOND) * Additional solver options SOLUTN(P1 ,Y,Y,Y,N,N,N) SOLUTN(H1 ,Y,Y,Y,N,N,Y) TURMOD(KEMODL) ************************************************************ Group 8. Terms & Devices TERMS (H1 ,N,Y,Y,Y,Y,Y) NEWRH1 = T ************************************************************ Group 9. Properties PRESS0 = 1.000000E+05 ;TEMP0 = 2.730000E+02 RHO1 = GRND6 RHO1A = 1.600000E+01 ;RHO1B = 2.900000E+01 RHO1C = 3.200000E+01 ENUL = 1.544000E-05 CP1 = GRND10 CP1A = 1.500000E+03 ;CP1B = 1.500000E+03 CP1C = 1.500000E+03 ;CP1D = 0.000000E+00 TMP1 = GRND8 TMP2A = 5.405406E-02 ;TMP2B = 4.900000E+07 TMP2C = 0.000000E+00 DVO1DT = 3.410000E-03 PRNDTL(H1 ) = 7.000000E-01 PRT (EP ) = 1.314000E+00 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(H1 ) = 4.000000E+04 ;FIINIT(FUEL) = 0.000000E+00 FIINIT(PROD) = 0.000000E+00 ;FIINIT(OXID) = 0.000000E+00 FIINIT(MIXF) = 1.000000E+00 ;FIINIT(EPKE) = 1.000000E+00 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 PATCH (CHSO ,PHASEM,0,0,0,0,0,0,1,1) COVAL (CHSO ,FUEL, GRND9 , GRND9 ) EGWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 300 RESFAC = 1.000000E-03 ************************************************************ Group 16. Terminate Iterations ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX, 7.000000E-01) RELAX(U1 ,FALSDT, 1.000000E-02) RELAX(V1 ,FALSDT, 1.000000E-02) RELAX(KE ,LINRLX, 3.000000E-01) RELAX(EP ,LINRLX, 3.000000E-01) RELAX(H1 ,FALSDT, 1.000000E-02) RELAX(DEN1,LINRLX, 1.000000E-02) RELAX(FUEL,FALSDT, 1.000000E-02) RELAX(MIXF,FALSDT, 1.000000E-02) KELIN = 3 ************************************************************ Group 18. Limits VARMAX(P1 ) = 1.000000E+10 ;VARMIN(P1 ) =-5.000000E+04 VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06 VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06 VARMAX(TMP1) = 3.000000E+03 ;VARMIN(TMP1) = 0.000000E+00 VARMAX(PROD) = 1.000000E+00 ;VARMIN(PROD) = 0.000000E+00 VARMAX(OXID) = 1.000000E+00 ;VARMIN(OXID) = 0.000000E+00 VARMAX(MIXF) = 1.000000E+00 ;VARMIN(MIXF) = 0.000000E+00 ************************************************************ Group 19. EARTH Calls To GROUND Station USEGRD = T ;USEGRX = T GENK = T ASAP = T CHSOA = 5.405406E-02 ;CHSOB = 1.000000E+00 CSG1 ='SW' ************************************************************ Group 20. Preliminary Printout ECHO = T ************************************************************ Group 21. Print-out of Variables OUTPUT(DEN1,Y,N,Y,N,Y,Y) OUTPUT(TMP1,Y,N,Y,N,Y,Y) OUTPUT(PROD,Y,N,Y,N,Y,Y) OUTPUT(OXID,Y,N,Y,N,Y,Y) ************************************************************ Group 22. Monitor Print-Out IXMON = 22 ;IYMON = 41 ;IZMON = 1 NPRMON = 100 NPRMNT = 1 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 NXPRIN = 1 NYPRIN = 1 IPLTF = 1 ;IPLTL = 100 NPLT = 1 ISWPRF = 1 ;ISWPRL = 100000 No PATCHes used for this Group ************************************************************ Group 24. Dumps For Restarts NOWIPE = T GVIEW(P,0.000000E+00,0.000000E+00,1.000000E+00) GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00) > DOM, SIZE, 1.200000E+00, 5.000000E-01, 2.000000E-01 > DOM, MONIT, 4.300000E-01, 3.375000E-01, 1.000000E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, SNAPSIZE, 5.000000E-03 > GRID, RSET_Z_1, 1, 1.000000E+00 > OBJ, NAME, OUTLET > OBJ, POSITION, 1.200000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 5.000000E-01, 2.000004E-01 > OBJ, CLIPART, cube12t > OBJ, ROTATION24, 1 > OBJ, GRID, 2 > OBJ, TYPE, OUTLET > OBJ, PRESSURE, 0.000000E+00 > OBJ, TEMPERATURE, SAME > OBJ, COEFFICIENT, 1.000000E+03 > OBJ, TURBULENCE, SAME , SAME > OBJ, NAME, IN3 > OBJ, POSITION, 0.000000E+00, 4.300000E-01, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 7.000000E-02, 2.000004E-01 > OBJ, CLIPART, cube3t > OBJ, ROTATION24, 1 > OBJ, GRID, 2 > OBJ, TYPE, INLET > OBJ, PRESSURE, 0.000000E+00 > OBJ, VELOCITY, 1.000000E+01, 0.000000E+00, 0.000000E+00 > OBJ, TEMPERATURE, 2.700000E+01 > OBJ, INLET_FUEL, 1.000000E+00 > OBJ, INLET_MIXF, 1.000000E+00 > OBJ, TURB-INTENS, 5.000000E+00 > OBJ, NAME, B11 > OBJ, POSITION, 2.000000E-01, 3.300000E-01, 0.000000E+00 > OBJ, SIZE, 1.500000E-01, 1.500000E-02, 2.000004E-01 > OBJ, ROT-ANGLE, -2.200000E+01, 0.000000E+00, 0.000000E+00 > OBJ, ROT-MODE, 1 > OBJ, CLIPART, cube > OBJ, ROTATION24, 1 > OBJ, GRID, 2 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 103 > OBJ, INI_TEMP, 2.700000E+01 > OBJ, NAME, B12 > OBJ, POSITION, 2.000000E-01, 3.300000E-01, 0.000000E+00 > OBJ, SIZE, 1.600000E-01, 1.500000E-02, 2.000004E-01 > OBJ, ROT-ANGLE, 2.200000E+01, 0.000000E+00, 0.000000E+00 > OBJ, ROT-MODE, 1 > OBJ, CLIPART, cube > OBJ, ROTATION24, 1 > OBJ, GRID, 2 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 103 > OBJ, INI_TEMP, 2.700000E+01 > OBJ, NAME, B7 > OBJ, POSITION, 4.400001E-01, 1.800000E-01, 0.000000E+00 > OBJ, SIZE, 1.600000E-01, 1.500000E-02, 2.000004E-01 > OBJ, ROT-ANGLE, 2.200000E+01, 0.000000E+00, 0.000000E+00 > OBJ, ROT-MODE, 1 > OBJ, CLIPART, cube > OBJ, ROTATION24, 1 > OBJ, GRID, 2 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 103 > OBJ, INI_TEMP, 2.700000E+01 > OBJ, NAME, B8 > OBJ, POSITION, 4.400001E-01, 1.800000E-01, 0.000000E+00 > OBJ, SIZE, 1.520000E-01, 1.500000E-02, 2.000004E-01 > OBJ, ROT-ANGLE, -2.200000E+01, 0.000000E+00, 0.000000E+00 > OBJ, ROT-MODE, 1 > OBJ, CLIPART, cube > OBJ, ROTATION24, 1 > OBJ, GRID, 2 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 103 > OBJ, INI_TEMP, 2.700000E+01 > OBJ, NAME, B13 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 5.719998E-01, 2.800000E-01, 2.000000E-01 > OBJ, CLIPART, wedge1 > OBJ, ROTATION24, 9 > OBJ, GRID, 2 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 102 > OBJ, INI_TEMP, 2.700000E+01 > OBJ, NAME, B9 > OBJ, POSITION, 6.280003E-01, 3.080000E-01, 0.000000E+00 > OBJ, SIZE, 5.719998E-01, 1.920000E-01, 2.000000E-01 > OBJ, CLIPART, wedge1 > OBJ, ROTATION24, 23 > OBJ, GRID, 2 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 103 > OBJ, INI_TEMP, 2.700000E+01 > OBJ, NAME, IN1 > OBJ, POSITION, 0.000000E+00, 2.800000E-01, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 7.000000E-02, 2.000004E-01 > OBJ, CLIPART, cube3t > OBJ, ROTATION24, 1 > OBJ, GRID, 2 > OBJ, TYPE, INLET > OBJ, PRESSURE, 0.000000E+00 > OBJ, VELOCITY, 1.000000E+01, 0.000000E+00, 0.000000E+00 > OBJ, TEMPERATURE, 2.700000E+01 > OBJ, INLET_FUEL, 1.000000E+00 > OBJ, INLET_MIXF, 1.000000E+00 > OBJ, TURB-INTENS, 5.000000E+00 > OBJ, NAME, IN2 > OBJ, POSITION, 0.000000E+00, 3.500000E-01, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 8.000000E-02, 2.000004E-01 > OBJ, CLIPART, cube3t > OBJ, ROTATION24, 1 > OBJ, GRID, 2 > OBJ, TYPE, INLET > OBJ, PRESSURE, 0.000000E+00 > OBJ, VELOCITY, 2.000000E+01, 0.000000E+00, 0.000000E+00 > OBJ, TEMPERATURE, 2.700000E+01 > OBJ, TURB-INTENS, 5.000000E+00 > OBJ, NAME, B14 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 1.200001E+00, 0.000000E+00, 2.000000E-01 > OBJ, CLIPART, cube13 > OBJ, ROTATION24, 1 > OBJ, TYPE, PLATE > OBJ, SURF_TEMP, 0.000000E+00, 2.700000E+01 > OBJ, NAME, B15 > OBJ, POSITION, 0.000000E+00, 5.000000E-01, 0.000000E+00 > OBJ, SIZE, 1.200001E+00, 0.000000E+00, 2.000000E-01 > OBJ, CLIPART, cube13 > OBJ, ROTATION24, 1 > OBJ, TYPE, PLATE > OBJ, SURF_TEMP, 0.000000E+00, 2.700000E+01