TALK=T;RUN( 1, 1)
#cls
display
1-D PARSOL HEAT CONDUCTION
Ý------------------Ý\\\\\\Ý--------------------------Ý
air AL air
z= 0. 0.55 0.75 1.4
BC.T=0 [k] T=1.2[K]
Analytical Solution
-----------------------------------------------------------
iz = 1 2 3 4 5 6 7
T = 0.05 0.15 0.25 0.35 0.45 0.55 0.55
-----------------------------------------------------------
iz = 8 9 10 11 12 13 14
T = 0.55 0.65 0.75 0.85 0.95 1.05 1.15
-----------------------------------------------------------
Note: This case was used to test PARSOL efficiency, the explicit
method was demonstrated so slow convergence. Implicit
method shows significant improvement.
Also, this case can verify the links settings between
1) cut-cell to non-cut solid cell
2) cut-cell to non-cut fluid cell
enddis
#pause
************************************************************
Q1 created by VDI menu, Version 3.4, Date 05/11/01
CPVNAM=VDI;SPPNAM=Core
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(Parallel plates with uniform temperature)
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,1,1,14)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(148) =PRPS ; NAME(149) =TEM1
* Solved variables list
SOLVE(TEM1)
* Stored variables list
STORE(PRPS)
* Additional solver options
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
SETPRPS(1, 0)
ENUT = 0.000000E+00
DRH1DP = 1.000000E-05
DVO1DT = 3.410000E-03
PRNDTL(TEM1) = -2.580000E-02
EMISS = 0.000000E+00 ;SCATT = 0.000000E+00
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(PRPS) = -1.000000E+00 ;FIINIT(TEM1) = 1.000000E-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
No PATCHes used for this Group
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 50
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(TEM1,FALSDT, 1.000000E+07)
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
ASAP = T
PARSOL = T
************************************************************
Group 20. Preliminary Printout
distil=t; dsttol=1e-03
EX(PRPS)= 8.071E+00
EX(TEM1)= 5.928E-01
ECHO = T
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 1 ;IZMON = 4
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
IzPRF = 1 ;IzPRL = 14
ISWPRF = 1 ;ISWPRL = 100000
IPROF = 3
PATCH (MYPLOT ,PROFIL,1,1,1,1,1,14,1,1)
PLOT(MYPLOT ,TEM1, 0.000000E+00, 0.000000E+00)
************************************************************
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.000000E+00, 1.000000E+00, 1.400000E+00
> DOM, MONIT, 5.000000E-01, 0.500000E-01, 3.500000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> GRID, RSET_X_1, 1, 1.000000E+00
> GRID, RSET_Y_1, 1, 1.000000E+00
> GRID, RSET_Z_1, 14, 1.000000E+00
> OBJ, NAME, B3
> OBJ, POSITION, 1.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 1.000000E+00, 1.400000E+00
> OBJ, CLIPART, cube11
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, PLATE
> OBJ, ADIABATIC, 0.000000E+00, 0.000000E+00
> OBJ, NAME, B5
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 1.000000E+00, 1.400000E+00
> OBJ, CLIPART, cube11
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, PLATE
> OBJ, ADIABATIC, 0.000000E+00, 0.000000E+00
> OBJ, NAME, B7
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 1.000000E+00, 0.000000E+00
> OBJ, CLIPART, cube13
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, SURF_TEMP, 0.000000E+00, 0.000000E+00
> OBJ, NAME, B8
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 1.400000E+00
> OBJ, SIZE, 1.000000E+00, 1.000000E+00, 0.000000E+00
> OBJ, CLIPART, cube13
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, SURF_TEMP, 0.000000E+00, 1.200000E+00
> OBJ, NAME, B6
> OBJ, POSITION, 0.000000E+00, 0.000000E+00,5.500000E-01
> OBJ, SIZE, 1.000000E+00, 1.000000E+00,2.000000E-01
> OBJ, CLIPART, cube14
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 100
> OBJ, ADIABATIC, 0.000000E+00, 0.000000E+00
M