TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 3.6, Date 21/06/05
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
************************************************************
IRUNN = 1 ;LIBREF = 2
************************************************************
Group 1. Run Title
TEXT(1D Evaporation of a Water Jet in a Duct )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,10,1,1)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(136) =VAPO ; NAME(137) =MASS
NAME(138) =HEAT ; NAME(139) =SPH1
NAME(140) =KOND ; NAME(141) =ENUL
NAME(142) =DEN1 ; NAME(143) =MOMX
NAME(149) =TEM1
* Solved variables list
SOLVE(P1 ,U1 ,VAPO,TEM1)
* Stored variables list
STORE(MOMX,DEN1,ENUL,KOND,SPH1,HEAT,MASS)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
************************************************************
Echo InForm settings for Group 7
inform7begin
(stored var tdec is tem1-273.)
inform7end
************************************************************
Group 8. Terms & Devices
TERMS (TEM1,Y,Y,Y,Y,N,Y)
NEWRH1 = T
NEWENL = T
DENPCO = T
************************************************************
Group 9. Properties
PRESS0 = 1.000000E+05
RHO1 = GRND5
RHO1A = 0.000000E+00 ;RHO1B = 3.484300E-03
RHO1C = 7.142900E-01
ENUL = GRND2
ENULA =-4.946790E-06 ;ENULB = 4.583940E-08
ENULC = 8.097400E-11
CP1 = 1.150000E+03
ENUT = 2.000000E-03
DRH1DP = GRND5
DVO1DT = 3.330000E-03
PRNDTL(TEM1) = -2.630000E-02
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(TEM1) = 1.051000E+03
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 (GENPAT ,CELL ,0,0,0,0,0,0,1,1)
COVAL (GENPAT ,U1 , FIXFLU , GRND )
COVAL (GENPAT ,VAPO, FIXFLU , GRND )
COVAL (GENPAT ,TEM1, FIXFLU , GRND )
PATCH (GENMAS ,CELL ,0,0,0,0,0,0,1,1)
COVAL (GENMAS ,P1 , FIXFLU , GRND )
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 150
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
ENDIT (TEM1) = 1.000000E-04
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 7.000000E-01)
RELAX(U1 ,FALSDT, 1.000000E-02)
RELAX(VAPO,LINRLX, 1.000000E-01)
RELAX(MASS,LINRLX, 1.000000E-01)
RELAX(HEAT,LINRLX, 1.000000E-01)
RELAX(DEN1,LINRLX, 1.000000E-01)
RELAX(MOMX,LINRLX, 1.000000E-01)
RELAX(TEM1,FALSDT, 1.000000E+00)
************************************************************
Group 18. Limits
VARMAX(P1 ) = 1.000000E+10 ;VARMIN(P1 ) =-8.000000E+04
VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06
VARMAX(DEN1) = 1.000000E+10 ;VARMIN(DEN1) = 1.000000E-05
VARMAX(TEM1) = 1.000000E+10 ;VARMIN(TEM1) = 2.730000E+02
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
ASAP = T
PARSOL = F
L($G001
GENTR = T
*------------------------------------------------------
* GENTRA GROUP 1: Particle physics
*------------------------------------------------------
* Particle type - 60
GPTYPE = 60
* Gravity components in GENTRA Cartesian system
GGRAX = 0.000000E+00 ;GGRAY = 0.000000E+00
GGRAZ = 0.000000E+00
* Buoyancy forces
GBUOYA = F ;GSURPR = F
* Stochastic model of turbulence
GSTOCH = F
* Data for vaporising particles
* Data for melting/solidifying particles
GDRAG = GRND1 ;GKONC = 2.630000E-02
GNUSS = GRND1
GKONV = GRND1 ;GCPCON = 1.150000E+03
GCPVAP = GRND1
GMWCON = 2.890000E+01 ;GMWVAP = 1.800000E+01
GLIQST = 2.731500E+02
GCPLIQ = 4.131800E+03 ;GHLIQD = GRND1
GVAPST = GRND1
GLATVP = GRND1 ;GSTPRE = GRND1
GDTRCT = 5.000000E-06
*-------------------------------------------------------
* GENTRA GROUP 2: Boundary conditions for particles
*-------------------------------------------------------
* Inlet-data file-name
GINFIL ='Q1'
*_XP__UP__DI_LDEN_MDOT_TEM_(NUM)
1.0 10. 10.e-5 1000. 2.5 298.
* Wall treatment, and rest coefficient if appropriate
GWALLC = 4
* Porosity threshold
GPOROS = 0.000000E+00
*-------------------------------------------------------
* GENTRA GROUP 3: Numerical controls
*-------------------------------------------------------
* 1st GENTRA sweep; frequency of calls
GSWEP1 = 1 ;GSWEPF = 1
* Maximum Lagrangian time-step; time step size multplier
GDTMAX = 1.000000E-03 ;GRTFRC = 2.000000E-01
* Min # of t-steps per cell; max # of t-steps; timeout
GLAGTS = 5 ;GSTEMX = 10000
GTIMMX = 1.000000E+02
*-------------------------------------------------------
* GENTRA GROUP 4: Output controls
*-------------------------------------------------------
* Restart-file, history-file and frequency for output
GRSFIL ='NONE'
GHFILE ='GHIS'
GOUTFR = 10
* The identifier of the individual history and
trajectory files
GH1STC ='H'
GT1STC ='T'
* The first, last particles and the interval for
writing history and trajectory files
NGWSTR = 0 ;NGWEND = 0 ;NGWINT = 1
GSWOUT = 1
* Particle number for residence-time calculation
GRESTI = 0
L(GENSET
************************************************************
Group 20. Preliminary Printout
ECHO = T
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 8 ;IYMON = 1 ;IZMON = 1
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NXPRIN = 1
IXPRF = 1 ;IXPRL = 8
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,0.000000E+00,-1.000000E+00,0.000000E+00)
GVIEW(UP,0.000000E+00,0.000000E+00,1.000000E+00)
> DOM, SIZE, 3.000000E+01, 5.000000E+00, 5.000000E+00
> DOM, MONIT, 2.250000E+01, 2.500000E+00, 2.500000E+00
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> GRID, RSET_X_1, 10, 1.000000E+00
> GRID, RSET_Y_1, 1, 1.000000E+00
> GRID, RSET_Z_1, 1, 1.000000E+00
> OBJ, NAME, INLET
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 5.000000E+00, 5.000000E+00
> OBJ, GEOMETRY, cube
> OBJ, ROTATION24, 9
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, VOLUFLOW, 7.541000E+01
> OBJ, TEMPERATURE, 1.051000E+03
> OBJ, GENTRA_EXIT, 1.000000E+00
> OBJ, NAME, B4
> OBJ, POSITION, 3.000000E+01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 5.000000E+00, 5.000000E+00
> OBJ, GEOMETRY, cube
> OBJ, ROTATION24, 9
> OBJ, GRID, NO
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, SAME
> OBJ, COEFFICIENT, 1.000000E+03
STOP
Energy balance for gas phase
----------------------------
25kg/s of hot gas enters at 798degC into a duct with a water
spray of 2.5kg/s at 25degC. The specific heat of the gas is
1150 J/kgK. The expected vapour mass fraction at the exit is
0.0909 and the expected exit temperature is 538degC.
Energy inflow rate of gas = 25*1150*1071 = 30.79 MW
Energy to heat droplets = 2.5*4183*(370-298) = 0.75 MW
Energy to vaporise droplets = 2.5*2.45E6 = 6.13 MW
Energy convected by evaporating
droplets into gas phase = 2.5*1880*370 = 1.74MW
Energy outflow rate of gas = 30.79 + 1.74 - 0.75 - 6.13
= 25.65 MW
Gas Exit temperature = 25.65E6/(1150*27.5) = 811 K
Hence, the expected gas exit temperature = 538degC.
distil=t
EX(P1 )=1.207E-01;EX(U1 )=2.564E+00;EX(VAPO)=9.082E-02
EX(MASS)=2.500E-01;EX(HEAT)=4.542E+05;EX(SPH1)=1.150E+03
EX(KOND)=2.630E-02;EX(ENUL)=8.572E-05;EX(DEN1)=4.289E-01
EX(MOMX)=2.500E+00;EX(TEM1)=8.123E+02;EX(TDEC)=5.393E+02